我正在开展一个大型项目,该项目涉及在城市之间创建节点和链接网络,然后将这些数据网格化到预设网格上。我已经将R栅格库用于此过程。
目前,栅格化数据与我预期的数据并不完全匹配(即应该覆盖城市的网格单元与城市实际所在的位置不匹配)。
cities http://i58.tinypic.com/2rwowuo.jpg
暗红色网格单元应该位于黑点之上,但它们似乎都向西移动了一个网格单元。
我怀疑在投影方面存在一些我遗漏的问题,并且' crs'光栅库中的选项,但我无法弄清楚我做错了什么,以及如何解决它。
更新:我不认为我可以挑出一大块解决问题的代码,但我的栅格流程如下:
rasterize(SpLin, fine.raster, field=*a number*, crs="proj=longlat", background=0)
图像来自于此(略微修剪):
pC <- ggplot(aes(x=long, y=lat, fill=value), data=dat_grid) + # dat_grid is the rasterized data
geom_tile() + scale_fill_gradient(low="white",high="red") +
geom_polygon(data=maine_state,aes(x=long,y=lat,group=group),colour="grey",fill="white",alpha=0) +
... (trimmed)
geom_point(data=points.df,aes(long,lat),size=2) # the black dots marking the cities