สร้างแวดวงใน PostGIS หรือไม่

ฉันใช้ PostGIS 1.5.2 กับรูปทรงเรขาคณิตใน SRID: 900913 ฉันต้องการสร้างแวดวงโดยใช้รายการคะแนนเป็นจุดศูนย์กลางด้วยรัศมี 600 กิโลเมตร ฉันใช้คำค้นหานี้:

INSERT INTO circles (geom) (
   SELECT ST_Buffer(point, 600000, 'quad_segs=8') 
   FROM points
);

แต่วงกลมที่สร้างขึ้นมีรัศมีไม่ 600 กิโลเมตร (รัศมีอยู่ใกล้กับความยาวนี้ แต่ไม่แน่นอน)

มีวิธีอื่นในการสร้างแวดวงใน PostGIS หรือไม่

หมายเหตุ: ข้อมูลที่แสดงนั้นมาจากสเปน การฉายภาพที่ถูกต้องคือ 4326 แต่ลูกค้าใช้ Google rasters ดังนั้นฉันจึงจัดเก็บข้อมูลไว้ใน 900913 เพื่อหลีกเลี่ยงการถอดแบบและเพิ่มประสิทธิภาพ

ลองสิ่งนี้:

 SELECT ST_Transform(geometry( 
            ST_Buffer(geography( 
                ST_Transform( point, 4326 )), 
                600000)), 
            900913) FROM points`

สิ่งนี้พลิกเป็นภูมิศาสตร์แล้วใช้การเลือก inbuilt SRID เพื่อ (ย้อนกลับ) ในรูปเรขาคณิตที่มีบัฟเฟอร์ระนาบที่ดีจะถูกเรียกใช้แล้วพลิกกลับ ปัญหาเกี่ยวกับวิธีการของคุณใน Mercator คือ Mercator ไม่รักษาระยะห่าง การใช้การฉายภาพในพื้นที่ที่เหมาะสมยิ่งขึ้นจะให้ผลลัพธ์ที่ดีกว่าและนั่นคือสิ่งที่เกิดขึ้นในวิธีการข้างต้น

คุณสามารถเพิ่มค่า quad_seg เพื่อให้ได้วงกลมที่แม่นยำยิ่งขึ้น เพื่อให้ได้วงกลมจริงคุณต้องทำงานกับรูปทรงเรขาคณิตโค้ง แต่ฉันไม่รู้ว่าซอฟต์แวร์ใดที่สามารถแสดงได้

ความแม่นยำของ srid 90013 นั้นก็แย่มากเช่นกันเนื่องจากเป็นภาพที่ครอบคลุมทั่วทั้งโลก

คุณจะได้ผลลัพธ์ที่แม่นยำยิ่งขึ้นด้วยการฉายภาพในพื้นที่

ขึ้นอยู่กับที่แวดวงของคุณกำลังถูกสร้าง พวกเขาอยู่ใกล้เส้นศูนย์สูตรหรืออยู่ใกล้กับเสาหรือไม่

ลองดูแผนที่นี้ คุณคิดว่า Antartica หรือกรีนแลนด์นั้นใหญ่จริงหรือ นั่นคือการฉายภาพที่คุณใช้ใช่ไหม

ฉันอยากจะแนะนำให้คุณอ่านเอกสาร USGS ฉบับนี้อย่างรวดเร็วเกี่ยวกับการฉายภาพโดยเฉพาะตารางด้านล่างซึ่งจะช่วยให้คุณทราบว่าการคาดการณ์แบบใดที่เหมาะกับอะไร

และหลังจากนั้นฉันคิดว่าในที่สุดฉันก็ควรตอบคำถามของคุณ :)

สิ่งที่ Nicklas พูดนั้นเป็นคำแนะนำที่ดี มีการประมาณการเฉพาะสำหรับพื้นที่ท้องถิ่นของคุณที่จะทำงานได้ดีขึ้นหรือไม่

มิฉะนั้นคุณอาจต้องการดูประเภท PostGIS Geographyใหม่ อย่างไรก็ตามคำตอบที่เหมาะสมที่สุดขึ้นอยู่กับว่าข้อมูลของคุณอยู่ที่ไหน

อัปเดต: ตั้งแต่ตอนนี้เรารู้ว่าข้อมูลของคุณอยู่ในสเปนคุณดูที่เก็บไว้ในเครื่องฉายภาพในพื้นที่เช่นUTM Zone 31Nทำการดำเนินการของคุณโดยใช้ข้อมูลนั้นแล้วฉายกลับไปที่ Google Web Mercator หรือไม่

คุณสามารถใช้ newish SQL / MM ตอนที่ 3ชนิดเรขาคณิต CIRCULARSTRING และ / หรือ CURVEPOLYGON

อย่างไรก็ตามขอเตือนว่าการสนับสนุนสำหรับประเภทนี้มี จำกัด ทั้งด้วยฟังก์ชั่นภายในและโปรแกรมภายนอก คุณสามารถใช้ST_CurveToLineเพื่อช่วยให้เห็นภาพได้

นอกจากนี้มันเป็นบิตนอกวาง ( ณ PostGIS 2.0 SVN) ที่SELECT ST_Area('CURVEPOLYGON(CIRCULARSTRING(2 1, 1 2, 0 1, 1 0, 2 1))'::geometry)เป็นเพียงการประมาณคร่าวๆของπ (เปรียบเทียบ 3.14033115695475 ถึง 3.14159265358979 จากpi())