PostGIS'de daire oluşturma?

PostGIS 1.5.2 kullanıyorum, geometrileri SRID: 900913. 600 kilometrelik yarıçapı olan merkez olarak bir nokta listesi kullanarak bir daire oluşturmam gerekiyor. Bu sorguyu kullanıyorum:

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

Ancak oluşturulan çevrelerin 600 kilometre yarıçapı yoktur (yarıçap bu uzunluğa yakındır, ancak tam olarak değildir).

PostGIS'te çevreler yaratmanın başka yöntemleri var mı?

NOT: Temsil edilen bilgiler İspanya'dandır. Doğru projeksiyon 4326'dır, ancak müşteri Google raster kullanmaktadır, bu nedenle yeniden yapılanmaları önlemek ve performansı artırmak için 900913’de veri depolarım.

Bunu dene:

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

Bu, coğrafyaya kayar, daha sonra, iyi bir düzlemsel tamponun çalıştığı geometriye geri dönmek için (ironik olarak) yerleşik SRID seçimini kullanır ve sonra geriye kayar. Mercator'daki yaklaşımınızın sorunu, Mercator'ın mesafeyi korumamasıdır. Daha uygun bir yerel projeksiyon kullanmak daha iyi sonuçlar verir ve yukarıdaki yöntemde olan budur.

Daha doğru bir daire elde etmek için quad_seg değerini artırabilirsiniz. Gerçek bir daire elde etmek için eğri bir geometri tipi ile çalışmak zorundasınız, ancak hangi yazılımı gösterebileceğini bilmiyorum.

90013 sridinin doğruluğu da çok kötü, çünkü tüm dünyayı kapsayan bir projeksiyon.

Yerel bir projeksiyon ile daha doğru bir sonuç alacaksınız.

Çevrelerinizin nerede oluşturulduğuna bağlıdır. Ekvatorun yanında mı yoksa kutuplara mı yakın?

Bu haritaya bir göz atın . Antartica veya Grönland’ın gerçekten bu kadar büyük olduğunu mu düşünüyorsun? Kullandığınız projeksiyon bu, değil mi?

Bu USGS belgesini projeksiyonlarda , özellikle de aşağıdaki tabloda size hangi projeksiyonların neyin iyi olduğu konusunda hızlı bir fikir veren hızlı bir şekilde okumanızı tavsiye ederim .

Ve hepsinden sonra, nihayet sorunuza cevap vermeliyim :)

Nicklas'ın söylediği şey iyi bir tavsiyeydi. Yerel bölgeniz için daha iyi çalışacak belirli bir projeksiyon var mı?

Aksi takdirde, yeni PostGIS Coğrafya türüne bakmak isteyebilirsiniz . Bununla birlikte, en uygun cevap, verilerinizin bulunduğu yere bağlıdır.

GÜNCELLEME: Artık verilerinizin İspanya'da olduğunu biliyoruz, UTM Zone 31N gibi yerel bir projeksiyonda saklamayı , işlemlerinizi kullanarak ve ardından Google Web Mercator'a yansıtmayı düşündünüz mü?

Newish SQL / MM Part 3 geometri tipini CIRCULARSTRING ve / veya CURVEPOLYGON kullanabilirsiniz.

Bununla birlikte, bu tür desteğin hem iç işlevlerde hem de dış programlarda sınırlı olduğu konusunda uyarılmalıdır. Görselleştirmeye yardımcı olmak için ST_CurveToLine kullanabilirsiniz .

Ayrıca, bir biraz kapalı (PostGIS 2.0 svn arasında gibi) SELECT ST_Area('CURVEPOLYGON(CIRCULARSTRING(2 1, 1 2, 0 1, 1 0, 2 1))'::geometry)yalnızca kaba bir yaklaşımdır tt . (3.14033115695475 ila 3.14159265358979'dan karşılaştırın pi()).