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

Question

我正在尝试在栅格（具有 OSGB36 投影）上绘制一些点，但点 x 和 y 位置与绘制的栅格的范围不同。我怎样才能让这些点与光栅处于同一坐标 system/extent，以便它们出现在光栅上。

栅格图层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

关于我可能做了什么 wrong/how 来修复它的任何想法？

谢谢！

Answer 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 中某个点的坐标

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

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

How to get point coordinates in same projection as raster

r

raster

spatial

r-raster