sql - 线串的平均方位

Question

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

例如，给定这个线串:

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



起点(线串中的第一个点)位于 map 底部(从墨西哥开始)。

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

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

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

我不确定该方法是否正确，或者是否有更好的方法(内置函数或算法)

Answer 1

我得到的角度是向北，但我的北角是 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