4

如何计算代表非直线线串几何的方位角？

例如，给定这个线串：

LINESTRING(-100.06372406847015 25.4583895045113,-95.14184906847015 30.197590679284158,-100.67895844347015 33.335334871825495,-95.40552094347015 36.576044042090466,-100.59106781847015 39.14304403636714,-95.40552094347015 41.61974326920709,-100.15161469347015 43.56067943577098,-95.66919281847015 45.87106856382069,-110.08325531847015 48.44030667059785,-85.38598969347015 48.73100364391479)

您可以在这里可视化：http: //arthur-e.github.io/Wicket/sandbox-gmaps3.html

在此处输入图像描述

起点（线串中的第一个点）位于地图底部（从墨西哥开始）。

我希望这条线的平均/整体方位指向北方（~0度）

大多数对线方位角的计算只是采用线串的起点和终点并以这种方式计算。这对于只有 2 个点的直线是正确的。

我相信对线串执行此操作的方法是计算每个 POINT 对的方位以及它们的平均值？

我不确定该方法是否正确或是否有更好的方法（内置函数或算法）

4

2 回答 2

2

我得到的角度是向北，但我的北方是 90 度，我从你的角度得到的角度是 84 度左右，所以向北有点东。注意：这可能不正确，但只是一个尝试。

我首先使用 split 获取单独的点，然后使用 xml 获取每个东和北点，然后重新排列数据，以便我可以在同一行上看到每个点与下一个点。然后计算增量，然后计算每个点之间的方位角，最后是平均值。希望能帮助到你。

declare @data nvarchar(max)='-100.06372406847015 25.4583895045113,-95.14184906847015 30.197590679284158,-100.67895844347015 33.335334871825495,-95.40552094347015 36.576044042090466,-100.59106781847015 39.14304403636714,-95.40552094347015 41.61974326920709,-100.15161469347015 43.56067943577098,-95.66919281847015 45.87106856382069,-110.08325531847015 48.44030667059785,-85.38598969347015 48.73100364391479'
declare @points table (seq int,x float,y float)

    ;with points as (select ROW_NUMBER() over (order by (select 1) desc) seq,* from string_split(@data,','))
    ,xmldata as (select points.seq,CONVERT(XML,'<Points><Point>'+ REPLACE(points.value,' ', '</Point><Point>') + '</Point></Points>') AS xmldataPoints
    from points),pointXY as (
    SELECT seq,
        xmldata.xmldataPoints.value('/Points[1]/Point[1]','float') AS [x],
        xmldata.xmldataPoints.value('/Points[1]/Point[2]','float') AS [y]
    FROM xmldata
    ),nextPoint as (
    select *,LEAD(x,1,null) over (order by seq) x2,LEAD(y,1,null)  over (order by seq) y2 from pointXY
    ),delta as (
    select *,(x2-x) dNorth,(y2-y) dEast from nextPoint
    ),bearing as (
    select * ,
    DEGREES(  IIF(dEast=0,
        IIF(dNorth<0,PI(),0),
        IIF(dEast<0,(-aTan(dNorth / dEast) + PI() / 2.0+ PI()),(-aTan(dNorth / dEast) + PI() / 2.0)))) bearing
    from delta
    )
    select AVG(bearing) from bearing

结果

84.5262691250142
于 2019-07-29T16:52:06.817 回答
0

当您像您一样排列线段时，每个线段代表一个向量，总线是一个合成向量。您可以找到向量的角度，找到起点和终点之间的角度；

制图方位角（十进制度）可以在平面上已知 2 个点的坐标（制图坐标）时计算：

在此处输入图像描述

来源：https ://en.wikipedia.org/wiki/Azimuth

于 2019-08-02T22:15:40.740 回答