赤道上的方位角是否等于不在赤道上的方位角？

我想完全理解方位角的概念，并遇到一些不一致（或者可能是我的错误）。赤道上的方位角是否等于不在赤道上的方位角？

我向你展示了一些不匹配的例子，希望有人能够解释我是如何工作的。

我在EPSG中显示坐标：900913，在PostGIS中使用我自己的JavaScript函数。

MY FUNCTION

/* Difference between the two longitudes */ 
var dLon = lon2 - lon1; 
/* Y value */ 
var y = Math.sin(dLon) * Math.cos(lat2); 
/* X value */ 
var x = Math.cos(lat1) * Math.sin(lat2) - Math.sin(lat1) * Math.cos(lat2) * Math.cos(dLon); 
/* Calculates the azimuth between the two points and converts it to degrees */ 
var angle = Math.atan2(y, x)/Math.PI * 180;

实施例

/* Same Y, not on the equator */ 
Point A: (-81328.998084106, 7474929.8690234) 
Point B: (4125765.0381464, 7474929.8690234) 
Result in PostGIS: 90 degrees 
Result in my JS function: 74.232 degrees 

/* Same Y, on the equator */ 
Point A: (-81328.998084106, 0) 
Point B: (4125765.0381464, 0) 
Result in PostGIS: 90 degrees 
Result in my JS function: 90 degrees

我明白，在赤道上，方位角是90（或270），用于一条水平线。认为如果你画一条水平线稍微偏向赤道北（或南），那么方位角不再是90度。但是... PostGIS告诉我，当我们有相同的Y时，总是90度。

此外，这个calculator还显示，当Y！= 0时，水平线的方位角不是90度（不在赤道上）。

它是如何正确的？

感谢

来源

2014-08-27 joaorodr84

[如何计算PostGIS中两点之间的方位角？]（http://stackoverflow.com/questions/25526684/how-to-calculate-the-azimuth-between-two-points-in- postgis） – 2014-08-28 17:19:57

这个问题似乎是题外话，因为它不是关于编程（尽管它使用编程）。请参阅帮助中心的[我可以询问哪些主题]（http://stackoverflow.com/help/on-topic）。也许[Geography Stack Exchange]（https://gis.stackexchange.com/）会是一个更好的地方。 – jww 2014-08-28 17:53:50

@jww，谢谢。我知道GIS.StackExchange更好。我尝试了两种方法，而且这里比较好。所以，非常感谢回答的人！ :)我将来会更加小心。 – joaorodr84 2014-08-28 18:06:48

在你的榜样，你已经使用EPSG：900913，这是平面的，预计在米。这意味着公式使用将是ATAN2，这永远是90时的纬度相同，与式为：

azimuth = atan2(y1-y2, x1-x2)

而第二部分将始终为0，得到的方位90.所以，使用平面坐标，是的，对于相同纬度的坐标对，方位角总是相同的，这就是为什么在使用EPS：900913时，Postgis总是给出相同的答案。

如果您切换到地理数据类型，并因此使用测地坐标，则不再是这种情况。

例如：

select degrees( 
    st_azimuth(
    st_makepoint(0, 10)::geography, 
    st_makepoint(90, 10)::geography));

给出80.1318065在PostGIS中，您链接计算器页面上给出了80.139。

随着x /经度越接近，对于给定的纬度，数值越接近90。例如，

select degrees( 
    st_azimuth(
    st_makepoint(0, 10)::geography, 
    st_makepoint(1, 10)::geography));

现在给出POSTGIS 89.9131737和和89.333的在线计算器（稍微差异）。

所有这些都是由于公式现在说明了曲率的事实，所以两个向量的投影之间的角度不会超过90度，除了在赤道上。

看看Wikipedia azimuth文章中有关球体版本的等式。这应该很容易在JavaScript中进行编码，并且应该与Postgis的地理类型给出类似的答案。

来源

2014-08-27 15:56:12

Hi @JohnBarça。我以错误的方式使用PostGIS。谢谢。正如你所看到的，我已经在使用[来自维基百科的功能]（http://upload.wikimedia.org/math/2/3/8/238805bcce98ded92525289f49b3d6f9.png），稍作修改。或者你的意思是[第二个]（http://upload.wikimedia.org/math/e/5/e/e5ecad955f1b1f7e2f84c3dd36eb4296.png）？ – joaorodr84 2014-08-27 16:40:33

我的意思是第二个，它代表了扁圆球体，从http://upload.wikimedia.org/math/7/6/e/76e90f1c786b852709ce30d9504cb19a.png开始。如果你使用它，你应该更接近你在Postgis中使用地理的值。 – 2014-08-27 16:43:55

非常感谢。 :) – joaorodr84 2014-08-27 16:56:26

赤道上的方位角是否等于不在赤道上的方位角？

回答

相关问题