在坐标系 WGS 74 和 WGS 84 之间转换

Question

我正在对 GSHHS 数据库进行栅格化，主要是将海岸线多边形和河流线转换为栅格。

http://www.soest.hawaii.edu/pwessel/gshhg/

河流和海岸数据库是两个不同的文件。

我注意到河流和海岸之间数百米的位置错位，而这些点显然应该对齐。我在README 中注意到的一件事是，海岸线数据库使用 WGS84 坐标，而河流数据库是使用 WGS 72 从其他来源生成的。区别应该在于本初子午线的偏移和地球主轴尺寸的差异。

我在互联网上搜索了关于这两个集合之间的转换，但没有找到。

我需要的答案：

• 如何在它们之间进行转换？

或者

• 如何解决 GSHHS 数据库中的错位问题？

Answer 1

以下是将 WGS 72 坐标转换为 WGS 84 坐标的公式和参数

公式

• Δφ" = (4.5 cos φ) / (a sin 1") + (Δf sin 2φ) / (sin 1")（单位 = 弧秒）
• Δλ" = 0.554（单位 = 弧秒）
• Δh = 4.5 sin φ + a Δf sin2 φ - Δa + Δr（单位 = 米）

---

• Δφ = Δφ" x 1"（单位 = 弧度）
• Δλ = Δλ" x 1"（单位 = 弧度）

参数

• Δf = 0.3121057 x 10-7
• a = 6378135 米
• Δa = 2.0 m
• Δr = 1.4 m
• 1" = pi / (3600.0 * 180.0) 弧度

要获得 WGS 84 坐标，请将使用 WGS 72 坐标计算的 Δφ、Δλ、Δh 变化添加到 WGS 72 坐标（分别为 φ、λ、h）。纬度为正北，经度为正东（0° 到 180°）。

结果

• φ WGS84 = φ + Δφ
• λ WGS84 = λ + Δλ
• h WGS84 = h + Δh

Documentation[p105-106]

Answer 2

您可以使用 proj4:https://trac.osgeo.org/proj/wiki/man_proj。支持许多椭球标识符。您可以使用 -le 选项开关获取列表。要从 wgs74 转换，您可以使用“towgs84”选项开关。

Answer 3

一个权威来源是国家地理空间情报局 (NGA)。他们在http://earth-info.nga.mil/GandG/geotrans/ 发布了开源的Geotrans 程序。这提供了 GUI，批处理。我希望内部有一个您可以作为开发人员调用的 API。