PROJ.4 库和 OSGB36

一切顺利

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

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

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

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

有任何想法吗？

提前谢谢了

leddy

谷歌搜索出现this曼彻斯特大学地球科学学者约翰·史蒂文森 (John Stevenson) 博士的成果——如果有人能做到的话，他应该能做对。这是一个引用。

问题是转向 OSGB36 需要投影和数据转换。之前2007年10月, proj 只携带 出投影，从而导致较大的偏移。你可以检查 如果您通过运行“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 是稍微扭曲关于 GPS 投影（例如 WGS84 和 ETRS89）。这个偏移量很小，并且 仅对于更高精度的测量很重要。很多搜索关于 OSGB36 偏移量会显示与此相关的页面。如果你想 也弥补这一点 您可以下载 nadgrid 文件并use it。对于我的数据，这移动了 点相差约 1 m。