Mengidentifikasi Sistem Koordinat Shapefile saat Tidak Diketahui?

Saya memiliki Shapefile tetapi sistem koordinatnya tidak diketahui, dan tidak ada file * .prj.

Bagaimana saya bisa mengidentifikasinya sekarang?

Selalu ada metode "brute force":

1. Ambil layer dengan sistem koordinat yang dikenal yang seharusnya overlay dengan layer Anda yang tidak dikenal.

2. Sekarang buat beberapa tebakan yang terpelajar tentang proyeksi apa yang bisa dilakukan oleh lapisan yang tidak dikenal. (UTM, Plat Carree, dll). Proyeksikan layer sistem koordinat yang diketahui ke dalam setiap proyeksi hingga Anda menemukan sebanyak mungkin yang cocok dengan layer yang tidak diketahui.

Saya mempromosikan komentar oleh Mapperz dan Brandon Copeland dengan menambahkan jawaban yang menggunakan teknik mereka.

1. Jika Anda memiliki ArcGIS Desktop, tambahkan beberapa data referensi yang mencakup area yang sama ke ArcMap. Data referensi ini harus memiliki definisi sistem (proyeksi) koordinat yang valid.
2. Tambahkan data dalam sistem koordinat yang tidak diketahui.
3. Teliti sistem koordinat apa yang mungkin digunakan di bidang ini. Satu tempat untuk memeriksa adalah Registry EPSG .
4. Atur sistem koordinat data frame ke kemungkinan. Jika data-dengan-coordsys-dikenal dengan data-tanpa-coordsys-dikenal, Anda mungkin telah mengidentifikasi sistem koordinat.

Teknik ini memanfaatkan fakta bahwa ArcMap dapat mengonversi data ke sistem koordinat baru dalam memori. Data-tanpa-diketahui-coordys tidak dapat dikonversi sehingga hanya ditampilkan. Dengan memproyeksikan data yang diketahui dalam memori, kita dapat melihat dengan cepat koordinat apa yang dimilikinya dalam berbagai sistem koordinat yang memungkinkan.

Karena pertanyaan ini tidak pernah menjadi tua, saya membangun situs yang menggunakan metode Brute Force . Jika Anda menyeret shp + shx yang zip ke peta itu akan memetakannya di setiap sistem koordinat yang tersedia di PostGIS. Dengan asumsi Anda tahu seperti apa "benar" itu, Anda dapat memperbesar ke area itu dan klik pada poligon untuk mendapatkan file .prj dari epsg.io.

Ada dua tautan hebat dari Esri yang menjelaskan tentang hal ini:

• Identifikasi sistem koordinat terproyeksi yang tidak diketahui menggunakan ArcMap
• Identifikasi referensi spasial, proyeksi, atau sistem koordinat data

Di pekerjaan terakhir saya, saya memiliki lapisan data Geologi (poligon) bernama "FSU_Geol.shp". Bos saya memberikannya kepada saya dan meminta saya untuk mencari tahu beberapa hal. Pertama-tama, dia menyerahkan formulir ini oleh klien, dan tidak ada file .prj, jadi dia ingin saya mengetahuinya. Dia juga ingin saya mencari tahu apa kategori Geologi itu. Saya bisa terus dan terus, tetapi mari kita memotong ke pengejaran .... Saya googled itu, dan saya berakhir di SINI .

Akronim dalam konvensi penamaan adalah singkatan dari "permukaan Geologi Bekas Uni Soviet", dan pencarian Google membawa saya "langsung" ke sumbernya (USGS). Semua yang saya perlu ketahui tentang shapefile ini adalah di link teratas yang saya tekan. Saya tidak mengatakan bahwa Google dapat menemukan apa saja & segalanya , tetapi saya baru lulus dari universitas, dan hanya mengambil "kesempatan dalam kegelapan", dan melihat umpan balik yang saya dapatkan!

Dalam contoh lain, seseorang di kantor telah mengunduh sejumlah shapefile menggunakan fungsi batch. Saya tidak ingat nama paket perangkat lunak saat ini, tetapi file .prj hilang. Saya hanya pergi ke metadata di ArcCatalog, dan saya benar-benar menemukan jawaban saya di sana. Sekali lagi, jangan mengutip saya, tapi saya pikir ada URL ke sumber di sana.

Sejak hari itu saya "mulai" dengan pencarian google jika tidak ada petunjuk di metadata!

Tanyakan pada produser.

Jika Anda tahu CRS umum untuk wilayah geografis Anda, Anda dapat mencobanya. Tapi bertanya itu lebih baik.

Saya sangat suka alat web ini: http://projfinder.com/ . Periksa file Anda dan cari beberapa koordinat. Zoom ke sekitar tempat yang sama di dunia dan biarkan menebak.

Bukan alat (saya tidak tahu ada yang akan membiarkan Anda melakukan itu), tetapi periksa balasan @ mkennedy untuk Bagaimana saya mengubah titik ini ke WKID 4326? . Dia menjelaskan bagaimana dia sampai pada referensi spasial yang benar. SpatialReference.org dan kesabaran akan menjadi teman Anda.

Selain itu, Esri menyediakan a panduan tentang cara menebak sistem koordinat (meskipun saya lebih suka metode mkennedy jika Anda tahu lebih banyak tentang data).

Jawaban terbaik yang saya temukan untuk pertanyaan ini adalah non-teknis: cari tahu dari mana data Anda berasal. Lembaga dan organisasi cenderung konsisten dengan penggunaan proyeksi mereka. Tahu itu datang dari DOT negara Anda? Lihatlah sisa data mereka dan lihat apa yang dikatakannya kepada Anda. Tidak tahu dari mana asalnya? Tebakan yang berpendidikan kemungkinan besar akan mengirim Anda ke jalan yang benar.

Setidaknya itu membuat mengatasi masalah dengan brute force sedikit lebih bisa dilakukan!

Serius sudah ketinggalan zaman pada saat ini, tetapi Werner Flacke dan Birgit Klaus memposting Cari Proyeksi pada ArcScripts pada 2007. Sayangnya, saya tidak menemukan kode sumber di sana. Ini berbasis VBA sehingga hanya dapat digunakan di ArcGIS Desktop v9.2 dan mungkin 9.3. Itu tidak termasuk dua shapefile dengan bidang yang diminati oleh EPSG Geodetic Parameter Dataset yang dapat digunakan untuk mempersempit kemungkinan.

Kalkulator Geografis Biru Marmer dan Transformer Geografis memiliki alat bantu pemulihan sistem yang terkoordinasi.

Coba program ogrinfo yang disediakan sebagai bagian dari GDAL.

Lihat Bagaimana cara mengakses metadata Shapefile menggunakan OGR?

Jadi ogrinfo mungkin tidak memberi Anda informasi proyeksi tanpa file .prj, tetapi masih merupakan alat yang berguna untuk membantu Anda menyelidiki daftar proyeksi potensial.

Sebagai contoh:

Geometry: Polygon
Feature Count: 269 
Extent: (320000.000000, 505000.000000) - (323000.000000, 511000.000000) 

Tanggapan ini menunjukkan kepada saya bahwa shapefile menggunakan sistem referensi berdasarkan meter dan bukan pada derajat.

Dengan asumsi Anda tahu dari mana data berasal, Anda sekarang memiliki daftar proyeksi yang lebih pendek.

Praktik terbaik lainnya mungkin untuk:

• Tanyakan pada orang-orang yang memberi Anda data di tempat pertama, bahkan jika mereka tidak tahu persis mereka harus memiliki daftar proyeksi yang mungkin bisa dilihat pada contoh pertama.
• Lihat apakah ada metadata terkait.

Mungkin lebih baik untuk memulai dari beberapa posisi yang diketahui, daripada mencoba dan kasar dari awal.

Untuk menambahkan tanggapan lain dalam posting ini, saya akan menambahkan yang berikut:

Cara: Mengidentifikasi sistem koordinat yang tidak diketahui yang diproyeksikan menggunakan ArcMap

Berikan cek khusus pada NAD1927

Jika koordinat dalam derajat desimal, seperti antara bujur -180 dan +180, dan lintang -90 dan +90, identifikasi Sistem Koordinat Geografis (datum) yang digunakan untuk data. Sebelum versi 9.2, ArcMap memberikan GCS_Assumed_Geographic_1 ke data, secara default. Ini menempatkan data pada datum NAD_1927

Tambahan

Jika data di Amerika Serikat dan menunjukkan sejauh mana koordinat di sebelah kiri desimal adalah 6, 7 atau 8 digit, data mungkin diproyeksikan ke State Plane atau sistem koordinat UTM.

FAQ: Dasar-Dasar Proyeksi: Apa yang perlu diketahui oleh profesional GIS

Halaman ini menyediakan beberapa pemeriksaan untuk menentukan proyeksi yang penting dalam membantu mencoba mengidentifikasi sistem koordinat yang tidak diketahui.

Konsep-konsep berikut ini mendasar untuk memahami penggunaan proyeksi peta di ArcGIS. Harap perhatikan bahwa topik proyeksi sangat luas, dan artikel ini tidak dapat melakukan lebih dari menyentuh beberapa topik penting.

1. Sistem koordinat, juga dikenal sebagai proyeksi peta, adalah sebutan yang sewenang-wenang untuk data spasial. Tujuannya adalah untuk menyediakan dasar bersama untuk komunikasi tentang tempat atau wilayah tertentu di permukaan bumi. Masalah paling kritis dalam berurusan dengan proyeksi peta adalah mengetahui apa proyeksi itu dan memiliki informasi sistem koordinat yang benar terkait dengan suatu dataset.

2. Ketika proyeksi peta pertama dibuat, diasumsikan, secara keliru, bahwa bumi itu datar. Kemudian anggapan itu direvisi, dan bumi dianggap sebagai bidang yang sempurna. Pada abad ke-18, orang mulai menyadari bahwa bumi tidak bulat sempurna. Ini adalah awal dari konsep spheroid kartografi.

3. Untuk mewakili lokasi di permukaan bumi dengan lebih akurat, pembuat peta mempelajari bentuk bumi (geodesi) dan menciptakan konsep spheroid. Kemudian sistem koordinat geografis (GCS) dirancang, yang meliputi datum, satuan ukuran, dan meridian utama. Datum menghubungkan spheroid ke bagian tertentu dari permukaan bumi. Datum terbaru dirancang agar sesuai dengan seluruh permukaan bumi dengan baik.

4. Datum yang paling umum digunakan di Amerika Utara adalah:

• NAD 1927 (Datum Amerika Utara 1927) menggunakan spheroid Clarke 1866
• NAD 1983 (Datum Amerika Utara 1983) menggunakan spheroid GRS 1980
• WGS 1984 (World Geodetic Survey 1984) menggunakan spheroid WGS 1984

Spheroids yang lebih baru dikembangkan dari pengukuran satelit dan lebih akurat daripada yang dikembangkan oleh Clarke pada tahun 1866. Istilah 'sistem koordinat geografis' dan 'datum' digunakan secara bergantian, tetapi seperti disebutkan di atas, GCS mencakup unit pengukuran datum, spheroid, satuan dan meridian utama.

5. Koordinat untuk perubahan data tergantung pada datum dan spheroid yang menjadi dasar koordinat tersebut, bahkan jika mereka menggunakan proyeksi dan parameter peta yang sama.

Misalnya, koordinat geografis di bawah ini adalah untuk satu titik yang terletak di dalam kota Bellingham, Washington, menggunakan 3 datum berbeda:

Kode: DATUM X -ordinate Y-Coordinate NAD_1927 -122.466903686523 48.7440490722656 NAD_1983 -122.46818353793 48.7438798543649 WGS_1984 -122.46818353793 48.7438798534299

6. Prinsip manajemen data yang baik adalah mendapatkan parameter proyeksi dari sumber data yang menyediakan data. Jangan membuat perkiraan berpendidikan tentang proyeksi data, karena database GIS yang tidak akurat akan menjadi hasilnya. Parameter yang diperlukan adalah sebagai berikut:

• Proyeksi
• Satuan ukuran
• ZONE (untuk UTM)
• Zona FIPS (untuk Pesawat Negara)
• Datum

Parameter lain mungkin diperlukan, tergantung pada proyeksi. Misalnya, proyeksi Albers dan Lambert memerlukan parameter berikut:

• Pararel standar 1, dalam derajat, menit dan detik (DMS)
• Pararel standar 2 (DMS)
• Meridian pusat (DMS)
• Lintang asal proyeksi (DMS)
• Arah timur dan satuan ukuran
• Arah utara dan satuan pengukuran
• Pergeseran-X dan satuan ukuran
• Pergeseran-Y dan satuan pengukuran

Saya memiliki masalah yang sama ketika saya sedang mengerjakan proyek saya (DEM). Salah satu DEM adalah koordinat yang tidak diketahui dan tidak dapat diproyeksikan, jadi saya melakukan 'georeferensi' menggunakan alat Georeferencing di ArcMap . Setelah Georeferencing, Anda akan mendapatkan file prj.file Anda.

Lihatlah situs web ini:

http://www.egger-gis.at/shapefile-projectionfinder/

Saya adalah pengembang alat gratis ini. Mungkin alat ini dapat membantu Anda menemukan dan menetapkan proyeksi shapefile Anda dengan benar. Ini didasarkan pada gagasan http://projfinder.com/ oleh Aaron Racicot.

Anda juga dapat mencoba solusi khusus ini untuk:

Austria (Jerman): https://www.data.gv.at/anwendungen/checkaustrianprojection/

Australia: https://maegger.github.io/map_australia.html

UTM - Zona: https://maegger.github.io/map_utm.html

Anda dapat membedakan antara koordinat geografis dan kartografi, lat / lon pertama dalam derajat, North / East kedua dalam meter. Jika file tersebut berbentuk Anda dapat membuang koordinat geometri dengan shapelib-1.3.0 dan melihat nilainya kemudian. ..

Anda bisa melakukannya dengan Fiona.

    import fiona as f
    a = fiona.open("C:\QGIS_ShapeFile\qgis\shafile_XXX.shp")
    print(a.crs)

Anda akan mendapatkan jawabannya sebagai

    {'init': 'epsg:4326'}