如何在与栅格相同的投影中获得点坐标

时间:2017-05-20 15:32:27

标签: r spatial raster r-raster

我正在尝试在栅格上绘制一些点(具有OSGB36投影),但点x和y位置与绘制的栅格的程度不同。如何使点与光栅处于相同的坐标系/范围内,以便它们出现在光栅上。

栅格图层crs:

extent      : 420000, 480000, 440000, 5e+05  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +datum=OSGB36 +units=m +no_defs +ellps=airy +towgs84=446.448,-125.157,542.060,0.1502,0.2470,0.8421,-20.4894 

绘制栅格:

enter image description here

然而,点(也在OSGB6中)的x / y位置是:

54, -1.3

关于我可能做错了什么/如何解决它的任何想法?

谢谢!

1 个答案:

答案 0 :(得分:1)

您的观点似乎是在wgs84地理坐标中,而不是在预测的Osgb6中。您可能错过了在Arcgis中指定的crs。你需要修改你的观点的crs。我假设您的栅格名为r

pt <- data.frame(x=54,y= -1.3)
coordinates(pt) <- ~x+y
projection(pt) <- "+init:epsg=4326"
pt_osgb <- spTransform(pt, CRS(projection(r)))

顺便说一句,您知道在使用函数locator(sp=TRUE)时可以获得R中某点的坐标。

抱歉,我使用的是智能手机,我无法测试这个脚本,但我认为这可能就是原因......