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

Question

我试图在栅格（具有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 

绘制的栅格：

这些点的x / y位置（也在OSGB6中）为：

54, -1.3

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

谢谢！

Answer 1

您的观点似乎是在wgs84地理坐标中，而不是在预计的Osgb6中。 您可能在Arcgis中错误指定了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中点的坐标。

抱歉，我在使用智能手机，无法测试此脚本，但是我认为这可能是原因...