gis - PROJ.4 库和 OSGB36

Question

希望你一切都好

我正在尝试使用 proj.4 库将纬度/经度坐标转换为 OSGB36 x 和 y。

有没有其他人成功地做到了这一点？我需要填充 srcPrj4String 和 destPrj4String 变量，例如

string srcPrj4String = "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs";
string destPrj4String = "+proj=utm +zone=11 +ellps=GRS80 +datum=NAD83 +units=m";

但我无法弄清楚 destPrj4String 应该是什么与 OSGB36 - 我知道数据应该是 +datum=OSGB36，但我尝试的一切都不起作用

有任何想法吗？

提前谢谢了

莱迪

Answer 1

谷歌搜索出现 this来自曼彻斯特大学地球科学学者约翰史蒂文森博士 - 如果有人这样做，他应该做对。这是一个报价。

问题是进入 OSGB36 需要投影和
datum conversion .之前 October 2007 , 项目只携带
出投影，从而导致较大的偏移。你可以检查
如果您通过运行“proj -v”或查看您的
.epsg 文件:

cat /usr/share/proj/epsg | grep -A 1 "British National Grid" 

# OSGB 1936 / British National Grid 
<27700> +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 
+y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs  <> 

新版本有 +datum=OSGB36。

如果您有旧版本，则可以通过将以下行替换为以下内容来更正它:

+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.999601 +x_0=400000 +y_0=-100000 
+ellps=airy 
+towgs84=446.448,-125.157,542.060,0.1502,0.2470,0.8421,-20.4894 +units=m 
+no_defs <> 

一个复杂的问题是 OSGB36 是 slightly distorted关于 GPS
预测(例如 WGS84 和 ETRS89)。这个偏移很小，并且
仅对更高精度的测量很重要。许多搜索关于
OSGB36 偏移量调出与此相关的页面。如果你想
也弥补了这一点，
您可以下载 nadgrid 文件和 use it .对于我的数据，这移动了
点大约 1 m。