将 NAD83 UTM 转换为纬度和经度时出现问题

Question

编辑 因此，我从 Jenner 在 VBforums 发布的解决方案显然是用 VB 编写的，我使用在线转换器将其移植到 C#。翻译过程中丢失了一些东西，这就是为什么它偏离了 10 英里。

我猜我只是误解了 Proj.Net 讨论板上的算法的用途，这就是为什么他们没有按照我的意愿行事。

我将在两天内结束这个问题，届时我可以标记我的答案，除非有人提供了很棒的东西。

---

我在将 UTM 转换为纬度和经度时遇到问题。例如，我有以下 NAD83 UTM 坐标：

Easting: 686029.702258

Northing: 3581213.621173

zone: 15

我在 vbforums 不久前给了我一个位于我预期的位置以南约 10 英里处。我在此处找到了一个简单的数学解决方案，给了我意想不到的结果。

double[] inverseMercator (double x, double y) {
     double lon = (x / 20037508.34) * 180;
     double lat = (y / 20037508.34) * 180;

     lat = 180/Math.PI * (2 * Math.Atan(Math.Exp(lat * Math.PI / 180)) - Math.PI / 2);
     return new double[] {lon, lat};
}
double[] toPoint = inverseMercator (686029.702258, 3581213.621173);

我得到以下结果：

lat: 30.602349476368449

long: 6.1627096689832594

我在同一线程中使用 D_Guidi 提供的 Proj.Net 解决方案得到类似的结果。

使用在线转换器，我能够得到更接近我所期望的结果：

lat：32.35238307052292

long：-91.0230710652583

任何人都可以阐明我做错了什么吗？

编辑 - 更喜欢 .NET 解决方案或易于转换的东西

相关线程

Answer 1

您是否尝试过“proj”中包含的cs2cs二进制文件（投影，http://trac.osgeo.org/proj/ ）？我相信这可以准确地将 UTM 转换为东向/北向和纬度/经度。问题是简单的数学公式并不能准确、足够地解释地球的形状。

Answer 2

因此，我从 Jenner 在 VBforums 上发布的解决方案显然是用 VB 编写的，我使用在线转换器将其移植到 C#。翻译过程中丢失了一些东西，这就是为什么它偏离了 10 英里。

我猜我只是误解了 Proj.Net 讨论板上的算法的用途，这就是为什么他们没有按照我的意愿行事。