用于检测 R 中空间自相关的坐标投影问题答案

【问题标题】：Issue with coordinate projection for detecting spatial autocorrelation in R用于检测 R 中空间自相关的坐标投影问题
【发布时间】：2017-03-19 15:33:10
【问题描述】：

我们有一个数据集，其中包含纬度和经度坐标以及属性信息，每个都在自己的单独列中，以数字形式存储。这些坐标已根据地理坐标系 WGS 1984 进行了地理编码。我们知道我们的数据中存在显着的空间自相关，我们希望使用“sp”包在气泡图中可视化。我们正在对其他在线示例进行建模，例如：https://beckmw.wordpress.com/2013/01/07/breaking-the-rules-with-spatial-correlation/。但是，当我们尝试在“sp”中使用坐标命令时，我们不断收到错误消息：

代码示例：坐标

if (nchar(projargs) == 0) projargs

我们无法加载我们的直接代码，因为它很敏感并且托管在无法访问互联网的虚拟环境中。有没有人知道为什么会发生这种情况？我们查看了 proj4 包，但不知道如何指定投影系统（或者这甚至是我们得到的错误？）。如果有人知道 R 中的任何其他包或可视化空间自相关的方法，那也将不胜感激。

【问题讨论】：

尝试使用不敏感的数据重新创建错误消息，并在可重现的示例中提供这些信息。为什么不直接使用coords = SpatialPoints(data.frame(lead$X, lead$Y), CRS("+init=epsg:4326"))

标签： r geospatial spatial projection sp

【解决方案1】：

您的代码有点“奇怪”：似乎您正在尝试构建一个仅包含坐标的数据集。 AFAIU，你可能需要这行的东西：

data <- data.frame(lead$X, lead$Y, lead$Z)

，其中 Lead$Z 对应于您要检查的通用“变量”，然后

coordinates(data) <- c('X','Y')`
proj4string(data) <- "+init=epsg:4326"

，它应该给你一个适当的“SpatialPointsDataframe”和lat-lon WGS84地理坐标（第一行也可以被删除，你将把所有变量保留在spatialpointsdataframe的数据中）。

HTH

【讨论】：