Web Mercator Projeksiyonunda Daha İyi Mesafe Ölçümleri

Katmanlarımı bir sql-uzamsal etkin-SDE-coğrafi veritabanında (Geometri türü, Web Mercator-3857) depolayan bir ESRI yığını ile çalışıyorum.

Ben bir web haritalama uygulaması inşa ediyorum, bu yüzden varsayılan olarak, fayans 3857 web mercator da bulunmaktadır.

Saklanan procs aracılığıyla, bir kullanıcının konumundan (web mercator da Koordinatlar) çeşitli katmanlara mesafeleri sorgulamak için STDistance kullanın.

Sorun, web mercatorunun bozulmasından dolayı, mesafe kireçlerimin gittikçe azalması, ekvatordan daha uzak hale getirilmesidir.

Katmanlarımı sql-uzamsal-coğrafya (geometri yerine) tipinde saklamayı düşündüm, ancak:

• Mesafe sorgularımın daha uzun süreceğini tahmin ediyorum (küresel yüzeyde mesafe kireçlenmeleri)
• çok fazla veriyi yeniden içe aktarmam gerekecek
• arcgis hizmetleri anında projeksiyon yapmaları gerektiği kadar hızlı olmayacak

Google haritalarına gidersem ve bir mesafe hesaplaması yaparsam, geri gönderilen mesafe, Kuzey / güney bölgelerinde bile çok daha doğru olur, bu nedenle Google'ın web mercator projeksiyonunun neden olduğu bozulmayı düzeltmesi gerektiğini varsayıyorum.

O zaman sorum: "Doğru" mesafeyi elde etmek için web mercator projeksiyonunda yapılan mesafe hesaplarına uygulanabilecek basit bir faktör değeri var mı?

Kısa mesafeler için, hesaplanan mesafeyi cos(lat)çarpan projeksiyonunun ölçeği enlemin sekantıyla orantılı olduğundan (sekant 1/cos) hesaplayabilirsiniz . Ayrıca bkz. Http://en.wikipedia.org/wiki/Mercator_projection#Mathematics_of_the_projection

@ Jeremiah-england sayesinde Zeyilname : yukarıdaki düzeltme gerçek Mercator projeksiyonları söz konusu olduğunda doğru olsa da, Web Mercator (EPSG: 3857) bir Mercator değildir. EPSG buna "Sözde Mercator" diyor. Sorun, eliptik bir model olan WGS84'ü kullanması ve küresel mercator hesaplamaları (Google'ın daha hızlı oldukları için kullandığı) kullanarak yansıtmasıdır. Mesafelerinizi 1/cos(phi)web Mercator ile ölçeklendirirseniz , ekvatorda yaklaşık% 0,6 indirim yapacaksınız. Daha fazla bilgi için Noel Zinn'in Web Mercator'daki sunumuna bakın.

Yukarıda belirtilen sunuma göre, web mercator koordinatlarından daha doğru mesafeleri hesaplamak için aşağıdaki yöntem kullanılabilir. Verilen dx- yatay koordinat farkı (WE yönü) ve dy- dikey koordinat farkı (SN yönü):

e = 0.081819191
adjustedX = dx * cos(lat) / sqrt(1 - e^2 * sin(lat)^2)
adjustedY = dy * cos(lat) * (1 - e^2) / pow(1 - e^2 * sin(lat)^2, 3/2)
adjustedDistance = hypot(adjustedX, adjustedY)

Bu ayar ile oran arasındaki oran ekvatordan kutuplara cos(lat)kadar SN yönünde daha büyüktür . WE yönü oranı ekvatorda başlar ve kutuplarda büyür .0.99331.003411.0034

Bu düzeltmenin hala sadece Dünya yüzeyinin düzlemsel geometrisinin kabul edilebildiği kısa mesafeler için oldukça iyi çalıştığını unutmayın.

GEOGRAPHYKüresel veriler üzerinde doğru sonuçlar arıyorsanız , verilerinizi tekrar formatta saklama seçeneğinizi düşünürüm .

Bir tabloda iki uzamsal alana sahip olmanızı engelleyen hiçbir şey yoktur - biri Mercator'da bir GEOMETRYtür ve diğeri WGS84'te bir GEOGRAPHYtür olarak (en azından SQL Server'da, ArcSDE hakkında emin değilim).

Her iki alanı da orijinal verilerinizle dolduracak basit bir coğrafi işlem komut dosyası oluşturabilmeniz gerekir. Düzenleme ve sık güncellemeler varsa, bu bir seçenek olmayabilir.

İki alana sahip olduğunuzda, hızlı görüntüleme için Mercator'unuzu kullanmaya devam edebilir ve mesafeler için kullanıcı tarafından girilen noktaları Enlem / Boylam'a dönüştürebilirsiniz.

Bunun iki büyük avantajı vardır:

• ayrıca kullanıcı sorgularına dayalı olarak daha doğru alanlar ve uzunluklar elde edebilirsiniz
• ölçümlerle ilgilenen her sorgu için özel kod eklemeyi unutma konusunda endişelenmenize gerek yoktur

Sorgu hızları daha karmaşık olacaktır, ancak bir kullanıcı tarafından fark edilmeyebilir. Bir çözüme karar vermeden önce bir özellik sınıfıyla test etmek isteyebilirsiniz. Ayrıca kullanıcı tarafından girilen noktaları Mercator'dan dönüştürmeniz gerekir (bu önemsiz olmalı ve tarayıcıda yapılabilir).

Coğrafya türüyle bile, yine de bir hata payı var:

Coğrafya yöntemleri için hata toleransı 1.0e-7 * kapsamları kadar büyük olabilir

Gönderen MSDN

ESRI'dan alternatif bir öneri, geometrileri bir ArcGIS Sunucusuna göndermeyi ve sonucun geri verilmesini içeren bir "Geometri servisi" kullanmaktır. Bir web hizmeti kullanarak herhangi bir darboğaz sorunu beklemiyorsanız, bu çok etkili bir yaklaşım olabilir. Aynı yaklaşımı bir Silverlight uygulamasında kullandım.

İşte ESRI'dan orijinal bir blog girişi: http://blogs.esri.com/Dev/blogs/arcgisserver/archive/2010/03/05/Measuring-distances-and-areas-when-your-map-uses-the -Mercator-projection.aspx

Blogda basitleştirilmiş bir JavaScript örneği var ve burada da tam bir örnek uygulama var: http://serverapps.esri.com/javascript_examples/compare_measurements.htm

Google Earth'ün yaptığı gibi, geometriyi yerel WGS84 bölge projeksiyonuna veya Güney Amerika'daki SIRGAS gibi gelişmiş bir WGS84 projeksiyonuna dönüştürebilirsiniz.

Hemen hemen tüm "bölgelendirilmiş" projeksiyonun koordinatları için http://www.spatialreference.org adresine bakabilir ve bunları dönüştürmek için bir tablo vb. Yapabilirsiniz.

Bir masa bölgeleri hayal ediyoruz

|   minx     |    miny    |    maxx    |   maxy     |  srid  |
+------------+------------+------------+------------+--------+
|234567.34314|234567.34334|234567.34334|234567.34334|  1234  |

Küresel Mercator'da (Web Mercator) saklanan tüm minx, miny, maxx, maxy değerleri. bu yüzden postgis'i örnek olarak kullanacağım.

SELECT ST_Distance(
         ST_Transform( -- transform/reproject
            ST_SetSRID(geom_line, 3857) -- geom_line with forced srid assignation to web mercator
            , ( -- here we get the first SRID from the spatial position of geom_line
                 SELECT  srid
                 FROM    zones
                 WHERE   ST_Contains(
                           ST_MakeBox2d(ST_Point(minx,miny),ST_Point(maxx,maxy))
                           , geom_line
                         )
                 LIMIT 1
            ))

Umarım bu faydalı olur, iyi günler.

OpenLayer'ların elipsoid üzerindeki iki EPSG: 4326 (WGS84) noktası arasındaki mesafeyi hesaplamak için bir yardımcı yöntem vardır. OpenLayers açık kaynak olduğundan, nasıl uygulandığını burada görebilirsiniz: https://github.com/openlayers/openlayers/blob/release-2.12/lib/OpenLayers/Util.js#L750

OpenLayers ve Measure kontrolünü kullananlar için, bu hesaplamayı kullanmak için litteral seçeneklerde jeodezik etkinleştirin.