3

很抱歉为此打扰您，但可能是 5 个小时我在这个问题上很生气，但我无法解决它。

我有一个包含大约 37,000 个条目的数据集。它们中的每一个都有自己的 Lon 和 Lat 坐标值。检查整体值，它们的范围分别如下：纬度（-54.4871,70.66344）和经度（-177.375,178.4419）。这是绝对合理的。

我使用 ArcGIS 创建了一个包含 37 千个点的 shapefile：一切正常。

然后我需要使用 R 处理这些数据，我用于代码的命令是（maptools 包）：

cells <- readShapeSpatial('RES',IDvar="id_obj", 
                          proj4string=CRS("+proj=longlat +datum=WGS84"))

但 R 给出了一个错误：

有效性方法中的错误（as（object，superClass））：给予不合格数据的地理 CRS：2.76663393422e+145

（我不知道这个数字来自哪里，它不是我数据集的一部分......）

阅读此博客上的其他帖子似乎原因应该是 lon 或 lat 的数据无效，但正如我上面提到的，我的数据集并非如此。

我尝试创建不同的shapefile，第一个没有投影，使用了几个投影（WGS84 Mercator，web mercator ...），但错误总是一样的......

谢谢你的帮助。

4

1 回答 1

6

底线是您的 shapefile 似乎已损坏。

点 shapefile 有两个主要部分，一个coords包含点坐标的部分，以及一个包含“属性”数据的数据部分（关于点的信息，如您的案例中的地区和国家/地区）。您的 shapefile 在数据部分也有 Lon 和 Lat，但它们不匹配：

library(rgdal)
setwd("<directory with shapefile...>")
map <- readOGR(dsn=".", layer="test")

range(map@data$Lat)
# [1] -54.48708  70.66344

range(map@coords[,2])
# [1]  -5.448708e+01  2.766634e+145

重投影涉及转换coords部分中的信息，这就是它失败的原因。

这是一种解决方法，但破解 SpatialPointsDataFrame 不是一个好主意：

map@coords <- as.matrix(map@data[c("Lon","Lat")])
map@bbox   <- rbind(range(map@coords[,1]),range(map@coords[,2]))
wgs.84 <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"
proj4string(map) <- CRS(wgs.84)

library(ggplot2)
gg <- data.frame(map@coords)
ggplot(gg) + 
  geom_point(aes(x=Lon,y=Lat), size=1, alpha=0.5, colour="blue") + 
  coord_fixed()


mercator <- "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs"
map.mercator <- spTransform(map,CRS=CRS(mercator))
gg <- data.frame(map.mercator@coords)
ggplot(gg) + 
  geom_point(aes(x=Lon,y=Lat), size=1, alpha=0.5, colour="green") + 
  coord_fixed()

我建议您重新创建 shapefile 并重试。

于 2014-01-29T21:35:57.373 回答