与R光栅化过程略有不匹配

时间:2014-11-24 20:29:12

标签: r projection raster

我正在开展一个大型项目,该项目涉及在城市之间创建节点和链接网络,然后将这些数据网格化到预设网格上。我已经将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

0 个答案:

没有答案