如何将以网格间隔收集的点数据转换为 r 中的地理参考数据集？

Question

我有这个数据集：https://www.dropbox.com/s/k06n9l05t25r6x2/newdata.csv?dl=0

（示例）

"","row","col","flagrv"
"1",2361,530,2
"2",2378,531,2
"3",2360,531,2
"4",2355,531,2
"5",2363,532,2
"6",2359,532,2
"7",2368,533,2
"8",2367,533,2
"10",2359,533,2

如果我使用此代码进行绘图：

gs.pal <- colorRampPalette(c("blue", "green","yellow","orange","red"),bias=1,space="rgb")
ggplot(data=ndata,aes(x=col,y=row,color=flagrv)) + 
  geom_point(size = 0.01)+
  scale_colour_gradientn(name = "Scale",colours = gs.pal(5))+
  xlab('Longitude')+
  ylab('Latitude')+
  theme_bw()+
  theme(line = element_blank())+
  theme(legend.position = c(.93,.20),panel.grid.major = element_line(colour = "#854440"))+
  ggsave("test.png",width=10, height=8,dpi=300)

我们得到这个数字：

现在，问题是我没有 Lat-Long 值。我想覆盖州边界但不能使用 Maps 包。有人建议我使用 gdal，但我不知道如何。您能否告诉我如何将其映射到 Lat-Long 域，以便我可以轻松地对其进行操作。

编辑：

我从别人那里了解到我可以使用这个：

gdal_translate -a_srs EPSG:4269 FILE.asc FILE.tif

#

答案 1 的错误

Error: unexpected ']' in "spdf = SpatialPointsDataFrame(coords, all_data[, c("flagrv"]"

然后我把代码改成：

spdf = SpatialPointsDataFrame(coords, all_data[, c("flagrv")]) 

但现在我有这个错误：

Error in validObject(.Object) : invalid class “SpatialPointsDataFrame” object: invalid object for slot "data" in class "SpatialPointsDataFrame": got class "integer", should be or extend class "data.frame"

Answer 1

在至少不知道数据集的投影和基准面的情况下（但希望有更多信息，例如分辨率和范围），没有简单的方法可以做到这一点。如果这是一个派生地图，请尝试查找用于生成它的内容。 有了这些信息，您就可以使用 raster 包中的投影函数来定义数据集的投影。

编辑（根据提供的其他信息，有一个可行的解决方案）： 这是一个可行的解决方案，因为数据集的左下角坐标为 24.55，-130，行/列之间的间距为 0.01 度，投影为 nad83。请注意，提供的元数据信息是错误的，因为最小纬度值不是 20 度，但可以从最南端（基韦斯特）估计为 24.55。

#load dataset 
all_data=(read.csv('new_data.csv',header=T, stringsAsFactors=F))
res=0.01 #spacing of row and col coords pre-specified
#origin_col_row=c(0, 0) 
origin_lat_lon=c(24.55, -130) 
all_data$row=(all_data$row)*res+origin_lat_lon[1] 
all_data$col=(all_data$col)*res+origin_lat_lon[2]

#now that we have real lat/lon, we can just create a spatial dataframe
library(rgdal)
library(sp)
coords = cbind(all_data$col, all_data$row)
spdf = SpatialPointsDataFrame(coords, data=all_data) #sp = SpatialPoints(coords)
proj4string(spdf) <- CRS("+init=epsg:4269") 

r 似乎在尝试绘制那么多点时窒息，所以为了检查答案是否有意义，我将数据集保存为 shapefile 并将其绘制在 arcgis 上：

writeOGR(spdf,"D:/tmp_shapefile4.shp", "flagrv", driver="ESRI Shapefile")

我设法使用 ggplot2 和下面的代码来绘制它，请耐心等待，因为绘制它需要一段时间：

df=as.data.frame(spdf)
library(ggplot2)
ggplot(data=df,aes(x=col,y=row,color=flagrv))+ 
   geom_point(size = 0.01)+
  xlab('Longitude')+
  ylab('Latitude')