Yanıtlar:
Her zaman "kaba kuvvet" yöntemi vardır:
Bilinmeyen katmanınıza binmesi gereken bilinen bir koordinat sistemine sahip bir katman alın.
Şimdi, bilinmeyen katmanın hangi projeksiyon olabileceği konusunda eğitimli tahminler yapın. (UTM, Levha Taşıyıcı, vb.) Bilinen koordinat sistem katmanınızı, her bir projeksiyonda, bilinmeyen katmana mümkün olduğunca uygun bir katman bulana kadar yansıtın.
Mapperz ve Brandon Copeland'ın tekniklerini kullanan bir cevap ekleyerek yorumlarını tanıtıyorum.
Bu teknik, ArcMap'in verileri bellekteki yeni bir koordinat sistemine dönüştürebilmesinden yararlanır. Bilinmeyen bir coordsys verisi dönüştürülemiyor, bu yüzden sadece görüntüleniyor. Bilinen verileri hafızaya yansıtarak, çeşitli koordinat sistemlerinde hangi koordinatlara sahip olduğunu hızla görebiliriz.
Bu soru hiç eskimez olduğundan, Brute Force yöntemini yapan bir site yaptım . Eğer sıkıştırılmış bir shp + shx'i haritaya sürüklerseniz, onu PostGIS'te bulunan her koordinat sisteminde eşler. "Doğru" neyin neye benzediğini bildiğinizi varsayarsak, o bölgeyi yakınlaştırıp poligonu tıklayarak .prj dosyasını epsg.io'dan edinin.
Bu konuda detaylara giren Esri'den iki harika bağlantı var:
Son işimde, "FSU_Geol.shp" adında bir jeoloji verisi (poligon) tabakası vardı. Patronum bana verdi ve birkaç şey bulmamı istedi. Öncelikle, bu şekil dosyasını müşteri tarafından teslim aldı ve .prj dosyası yoktu, bu yüzden bunu çözmemi istedi. Ayrıca Jeoloji kategorilerinin ne olduğunu bulmamı istedi. Devam edip gidebilirdim ama kovalamacayı kesmemize izin verdi .... Bunu googledim ve BURADA bittim .
Adlandırma sözleşmesindeki kısaltma, "Eski Sovyetler Birliği'nin yüzey jeolojisi" nde duruyordu ve Google araştırması beni "doğrudan" kaynağa yönlendirdi (USGS). Bu şekil dosyası hakkında bilmem gereken her şey, vurduğum en üst bağlantıdaydı. Google’ın her şeyi bulabildiğini söylemiyorum , ama üniversiteden yeni çıktım ve sadece "karanlıkta bir atış" çektim ve aldığım geribildirimlere baktım!
Başka bir durumda, ofisteki bir kişi toplu iş işlevini kullanarak birkaç şekil dosyası indirmişti. Şu anda yazılım paketinin adını hatırlamıyorum, ancak .prj dosyaları eksikti. ArcCatalog'daki metadataya girdim ve cevabımı orada buldum. Yine, bana teklif etme, ama oradaki kaynağın bir URL'si olduğunu düşünüyorum.
O günden beri meta verilerde ipucu yoksa google aramaya "başladım"!
Bu web aracını gerçekten çok seviyorum: http://projfinder.com/ . Dosyanızı kontrol edin ve bazı koordinatlara bakın. Dünya üzerinde yaklaşık olarak aynı yere zum yapın ve tahmin edin.
Değil bir araç (bunu yapmana izin verirdim varlığı birinin bilmiyorum) ama check out @ mkennedy yanıtını için 4326 WKID için bu noktayı dönüştürmek istiyorsunuz nasıl? . Mekansal referansa nasıl ulaştığını açıklıyor. SpatialReference.org ve sabır arkadaşlarınız olacak.
Ek olarak, Esri bir koordinat sisteminin nasıl tahmin edileceğine dair bir rehber sunar (veri hakkında biraz daha fazla şey biliyorsanız mkenedy yöntemini tercih etsem de).
Bu soruya bulduğum en iyi cevap teknik değil: verilerinizin nereden geldiğini bulun. Ajanslar ve kuruluşlar, projeksiyonları kullanmalarıyla tutarlı olma eğilimindedir. Senin devletin DOT'tan geldiğini biliyor musun? Verilerinin geri kalanına bak ve sana ne söylediğini gör. Nereden geldiğini bilmiyor musun? Eğitimli bir tahmin, sizi doğru yoldan aşağı gönderir.
En azından kaba kuvveti olan problemin üstesinden gelmek biraz daha mümkün!
Bu noktada ciddi olarak modası geçmiş, ancak Werner Flacke ve Birgit Klaus , 2007'de ArcScripts'te Projeksiyon Bul'u yayınladı . Ne yazık ki kaynak kodun orada olduğunu sanmıyorum. VBA tabanlıdır, bu nedenle yalnızca ArcGIS Desktop v9.2 ve muhtemelen 9.3'te kullanılabilir. EPSG Geodetic Parameter Dataset'in olasılıkları daraltmak için kullanılabilecek ilgi alanları olan iki şekil dosyası içerir.
Blue Marble Geographic Calculator ve Geographic Transformer, sistem kurtarma araçlarını koordine ediyor.
GDAL'in bir parçası olarak verilen ogrinfo programını deneyin .
Bkz Öğr kullanarak Şekil dosyası meta nasıl erişilir?
Bu nedenle, ogrinfo size .prj dosyası olmadan projeksiyon bilgisini vermeyebilir, ancak potansiyel projeksiyonların listesini araştırmanıza yardımcı olmak için hala yararlı bir araçtır.
Örneğin:
Geometry: Polygon
Feature Count: 269
Extent: (320000.000000, 505000.000000) - (323000.000000, 511000.000000)
Bu cevap bana, şekil dosyasının derecelere göre değil, sayaçlara dayanan bir referans sistemi kullandığını gösteriyor.
Verilerin kabaca nerede olduğunu bildiğinizi varsayarsak, şimdi olası projeksiyonların daha kısa bir listesine sahipsiniz.
Diğer en iyi uygulamalar şunlar olabilir:
Muhtemelen sıfırdan denemek ve kaba kuvvet kullanmak yerine, bazı bilgilerden başlamak daha iyidir.
Bu yayındaki diğer yanıtları eklemek için aşağıdakileri eklerdim:
Nasıl Yapılır: ArcMap kullanarak bilinmeyen, öngörülen bir koordinat sistemini belirleme
NAD1927'de özel bir kontrol sağlayın
Koordinatlar -180 ve +180 boylamları ve -90 ve +90 enlemleri gibi ondalık derecelerde ise, veriler için kullanılan Coğrafi Koordinat Sistemini (datum) tanımlayın. 9.2 sürümünden önce, ArcMap varsayılan olarak verilere GCS_Assumed_Geographic_1 atar. Bu, verileri NAD_1927 referans noktasına yerleştirir.
Ek olarak
Veriler ABD’de ise ve ondalık basamağın solundaki koordinatların 6, 7 veya 8 basamaktan oluştuğunu gösteriyorsa, veriler muhtemelen Devlet Düzlemine veya UTM koordinat sistemlerine yansıtılır.
SSS: Projeksiyon Temelleri: GIS uzmanının bilmesi gerekenler
Bu sayfa, bilinmeyen bir koordinat sistemini tanımlamaya çalışırken önemli olan projeksiyonları belirleme konusunda birkaç kontrol sunar.
ArcGIS'te harita projeksiyonlarının kullanımını anlamak için aşağıdaki kavramlar temeldir. Lütfen projeksiyonun konusunun oldukça geniş olmasına rağmen, bu makalenin birkaç önemli konuya değinmekten fazlasını yapamayacağını unutmayın.
Harita projeksiyonları olarak da bilinen koordinat sistemleri, mekansal veriler için keyfi tanımlamalardır. Amaçları, dünya yüzeyindeki belirli bir yer veya alan hakkında iletişim için ortak bir temel sağlamaktır. Harita projeksiyonları ile ilgili en kritik sorun, projeksiyonun ne olduğunu bilmek ve veri setiyle ilgili doğru koordinat sistemi bilgisine sahip olmaktır.
İlk harita projeksiyonları tasarlandığında, yanlış bir şekilde, dünyanın düz olduğu varsayılmıştır. Daha sonra varsayım revize edildi ve dünyanın mükemmel bir küre olduğu varsayıldı. 18. yüzyılda, insanlar dünyanın tam olarak yuvarlak olmadığını fark etmeye başladılar. Bu kartografik sferoid kavramının başlangıcıydı.
Harita yapıcılar, dünya yüzeyindeki yerleri daha doğru bir şekilde temsil etmek için, dünyanın şeklini (jeodezi) inceledi ve sferoid kavramını yarattı. Daha sonra, bir veri, ölçü birimleri ve bir ana meridyen içeren coğrafi koordinat sistemleri (GCS) geliştirilmiştir. Veri, bir sferoyu dünya yüzeyinin belirli bir bölümüne bağlar. Son veriler, tüm dünyanın yüzeyine iyi uyacak şekilde tasarlanmıştır.
Kuzey Amerika'da en yaygın kullanılan veriler:
• Clarke 1866 sferoidini kullanan NAD 1927 (Kuzey Amerika Datum 1927)
• GRS 1980 sferoidini kullanan NAD 1983 (Kuzey Amerika Datum 1983)
• WGS 1984 sferoidini kullanan WGS 1984 (Dünya Jeodezi Anketi 1984).Daha yeni sferoidler, uydu ölçümlerinden geliştirilmiştir ve 1866'da Clarke tarafından geliştirilenlere göre daha doğrudur. 'Coğrafi koordinat sistemi' ve 'datum' terimleri, birbirlerinin yerine kullanılır, ancak yukarıda belirtildiği gibi, bir GCS, bir veri, sferoid, ölçü birimlerini içerir. ve birinci sınıf bir meridyen.
- Verilerin koordinatları, aynı harita projeksiyonunu ve parametrelerini kullanıyor olsalar bile, bu koordinatların dayandığı referans ve küreye bağlı olarak değişir.
Örneğin, aşağıdaki coğrafi koordinatlar, Washington eyaletinin Bellingham şehrinde bulunan 3 farklı veriyi kullanan tek bir nokta içindir:
Kod: DATUM X-Koordinat Y-Koordinat NAD_1927 -122.466903686523 48.7440490722656 NAD_1983 -122.46818353793 48.7438798543649 WGS_1984 -122.46818353793 48.7438798534299
- İyi bir veri yönetimi ilkesi, projeksiyon parametrelerini veri sağlayan veri kaynağından elde etmektir. Verilerin projeksiyonu hakkında bilinçli bir tahminde bulunmayın, çünkü yanlış bir GIS veri tabanı ortaya çıkar. Gerekli parametreler şunlardır:
• Projeksiyon
• Ölçü birimleri
• ZONE (UTM için)
• FIPS bölgesi (Devlet Düzlemi için)
• DatumProjeksiyona bağlı olarak başka parametreler de gerekebilir. Örneğin, Albers ve Lambert projeksiyonları aşağıdaki parametreleri gerektirir:
• 1. standart paralel, derece, dakika ve saniye cinsinden (DMS)
• 2. standart paralel (DMS)
• Merkezi meridyen (DMS)
• Çıkıntıların başlangıç noktası (DMS)
• Yanlış doğu ve ölçü birimleri
• Yanlış kuzey ve ölçü birimleri
• X-shift ve ölçü birimleri
• Y-shift ve ölçü birimleri
Bu web sitesine bakın:
http://www.egger-gis.at/shapefile-projectionfinder/
Bu ücretsiz aracın geliştiricisiyim. Belki bu araç, şekil dosyanızın doğru izdüşümünü bulmanıza ve tanımlamanıza yardımcı olabilir. Aaron Racicot tarafından http://projfinder.com/ fikrine dayanmaktadır.
Aşağıdakiler için bu özel çözümleri deneyebilirsiniz:
Avusturya (Almanca): https://www.data.gv.at/anwendungen/checkaustrianprojection/
Avustralya: https://maegger.github.io/map_australia.html
UTM - Bölgeler: https://maegger.github.io/map_utm.html
Fiona ile yapabilirsin.
import fiona as f
a = fiona.open("C:\QGIS_ShapeFile\qgis\shafile_XXX.shp")
print(a.crs)
Olarak cevap alacaksınız
{'init': 'epsg:4326'}