วิธีการคำนวณมุมที่สองบรรทัดตัดกันใน PostGIS

ฉันต้องการคำนวณมุมระหว่างสองบรรทัดที่พวกเขาตัดกันใน PostGIS

จุดเริ่มต้นสำหรับการคำนวณมุมใน PostGIS นั้นน่าจะเป็นST_Azimuth - แต่นั่นก็เป็นจุดที่ป้อนเข้า ความคิดแรกของฉันคือการหาจุดสิ้นสุดของเส้นตัดกันและทำการคำนวณแบบอะซิมุ ธ นั่นไม่ดีพอเพราะคุณสมบัติส่วนใหญ่ของเส้นไม่ตรงและฉันสนใจในมุมที่จุดตัด ดังนั้นสิ่งที่ฉันเกิดขึ้นคือการดำเนินการซ้อนกันที่ทำตามขั้นตอนต่อไปนี้:

1. ระบุจุดแยกทั้งหมดระหว่างตารางคุณลักษณะสองบรรทัด
2. สร้างบัฟเฟอร์ขนาดเล็กมากรอบ ๆ จุดตัด
3. ระบุจุดที่คุณสมบัติของเส้นตัดกันด้านนอกของบัฟเฟอร์ (ใช้จุดแรกถ้ามีมากกว่าหนึ่ง - ฉันสนใจเฉพาะจริงๆว่ามุมนั้นใกล้กับ 0, 90 หรือ 180 องศา)
4. คำนวณ ST_Azimuth สำหรับสองจุดนั้น

SQL แบบเต็มมีความยาวในการโพสต์ที่นี่ แต่ฉัน gisted ที่นี่ถ้าคุณสนใจ (โดยวิธีการมีวิธีที่ดีกว่าที่จะดำเนินการมากกว่าเขตข้อมูลทั้งหมดจะลงงบกับ?)

ผลลัพธ์ไม่ถูกต้องดังนั้นฉันจึงทำสิ่งที่ผิดอย่างชัดเจน:

แก้ไขฉันเพิ่มการคำนวณใน EPSG ใหม่: 3785 และผลลัพธ์จะแตกต่างกันเล็กน้อย แต่ก็ไม่ถูกต้อง:

คำถามของฉันอยู่ที่ข้อบกพร่องในกระบวนการนี้ ฉันเข้าใจผิดว่า ST_Azimuth ทำอะไรได้บ้าง มีปัญหา CRS หรือไม่ มีอะไรอีกบ้างไหม? หรืออาจจะมีวิธีที่ง่ายกว่านี้มากในการทำเช่นนี้?

ฉันมีความศักดิ์สิทธิ์ มันค่อนข้างธรรมดา ฉันกำลังปล่อยข้อมูลสำคัญหนึ่งชิ้นสำหรับ PostGIS เพื่อคำนวณมุมที่ถูกต้อง

สิ่งที่ฉันคำนวณคือมุมระหว่างจุดสองจุดที่ตัดกันด้านนอกของบัฟเฟอร์เล็ก ๆ ในการคำนวณมุมของจุดตัดฉันต้องคำนวณมุมทั้งสองระหว่างจุดทั้งสองที่ด้านนอกของบัฟเฟอร์และจุดตัดของคุณลักษณะสองบรรทัดแล้วลบออก

ฉันอัปเดตSQLเต็มรูปแบบแต่นี่คือบิตเด่น

SELECT
    ...
    abs
    (
        round
        (
            degrees
            (
            ST_Azimuth
            (
                points.point2,
                points.intersection
            )
            -
            ST_Azimuth
            (
                points.point1,
                points.intersection
            )
        )::decimal % 180.0
        ,2
    )
)
AS angle
...
FROM
points 

ฉันเพิ่งต้องคำนวณสิ่งเดียวกัน แต่ตัดสินใจด้วยวิธีที่ง่ายกว่าและน่าจะเร็วกว่า

เพื่อหาคะแนนพิเศษสำหรับการคำนวณแบบราบฉันเพียงตรวจสอบความยาวด้านหลังสี่แยก (หรือหลังจากนั้นในกรณีที่เกิดขึ้นได้ยากในจุดเริ่มต้นของบรรทัด) โดยใช้ST_Line_Locate_PointและST_Line_Interpolate_Point :

abs(degrees( 
  ST_Azimuth (
    intersection, 
    ST_Line_Interpolate_Point(
      line1, 
      abs(ST_Line_Locate_Point(line1, intersection) - 0.0001)
    )
  )
  -
  ST_Azimuth (
    intersection, 
    ST_Line_Interpolate_Point(
      line2, 
      abs(ST_Line_Locate_Point(line2, intersection) - 0.0001)
    )
  )
))

permyriad นั้นเป็นกฎเกณฑ์และเพื่อผลลัพธ์ที่สอดคล้องกันมากขึ้นมันจะดีกว่าถ้าใช้ offset แบบสัมบูรณ์ สำหรับตัวอย่างเช่นตรวจสอบ 20m ล่วงหน้าคุณจะเปลี่ยน 0.0001 เป็น20/ST_Length(line1)และ20/ST_Length(line2)ตามลำดับ