Mengidentifikasi Sistem Koordinat Shapefile saat Tidak Diketahui?


Jawaban:


49

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.


1
ini juga metode saya. Saya menggunakan dua perangkat lunak yang berbeda untuk "membuktikannya". Anda juga dapat membaca dengan teliti pertanyaan ini untuk beberapa arahan dasar pada proyeksi. gis.stackexchange.com/questions/2769/…
Brad Nesom

6
Anda dapat menyederhanakan langkah 2 sedikit dengan hanya mengubah proyeksi dataframe (ArcMap), dan membiarkan proyek yang dikenal memproyeksi ulang dengan cepat.
Brandon Copeland

13
Perhatikan bahwa Anda tidak boleh memproyeksikan data. Anda harus mendefinisikan proyeksi. Memproyeksikan data mengubah koordinat. Menentukan proyeksi hanya memberikan proyeksi ke data itu dan memberi tahu perangkat lunak itu tentang data yang sudah ada. Ini adalah dua alat berbeda di ArcMap, "Tentukan Proyeksi" adalah yang ingin Anda gunakan.
Sean

Saya lebih suka membaca nilai koordinat, saya biasanya bisa tahu apakah itu derajat desimal, meter atau kaki (lebih sulit dibedakan) dan kemudian pindah ke sistem terbaik. Mulai dengan WGS 84 jika muncul desimal, UTM jika meter, Nyatakan bidang jika kaki. Saya pernah punya satu yang tidak ada yang bisa memperbaikinya (untuk anomali gravitasi). Itu pesawat negara bagian AS dengan unit linier ...... kilometer. Saya bisa tahu dari fakta bahwa itu dekat saya memiliki sistem yang benar tetapi negara itu kecil, ahh kilometer adalah unit linier.
Jika Anda tidak tahu- hanya GIS

33

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.


4
Catatan samping: Dalam QGIS memproyeksikan kembali dalam memori disebut "on-the-fly".
underdark

3
Ini disebut 'on-the-fly' di ArcGIS juga! Saya terus menemukan banyak orang yang tidak mengerti frasa itu, jadi saya mencoba menggunakan 'di-memori' sebagai gantinya.
mkennedy

Apakah tautan ini masih valid? Tampaknya sedikit lebih berbelit-belit daripada yang saya ingat, namun ada epsg-registry.org
GISKid

1
@GISKid Anda benar. Saya pikir situs web epsg-registry.org setelah tanggal ketika saya menulis jawaban ini. Saya akan memperbarui.
mkennedy

1
Sumber hebat lain: epsg.io
GISKid

29

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.


15
ide yang luar biasa, tapi saya pikir situs web itu tidak benar-benar berfungsi
Fábio Dias


22

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!


18

Tanyakan pada produser.

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


Satu-satunya masalah adalah saya tidak bisa bertanya padanya ... Ini juga pikiran pertama yang saya miliki.
com

Kemudian, Anda harus pergi brute force seperti kata brenth.
underdark

18

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.


Ini sebenarnya alat terbaik yang saya temukan. Tebak dalam ArcMap sebenarnya tidak berguna dalam jangka panjang atau jika Anda tidak tahu apa yang seharusnya. Tidak bisa membantu dengan sistem koordinat homebrew dan kekurangan satu sama sekali, tetapi siapa yang bisa!
Martin

Sayang sekali tidak lagi berfungsi - ada yang tahu siapa yang harus dihubungi?
sventechie


13

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!


11

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.


Terima kasih atas tipnya Melita. Kode ada di sana (* .frm). Lisensi adalah domain publik. PrjFinder adalah pendamping buku "Bekerja dengan Proyeksi dan Transformasi Datum di ArcGIS" (2006).
matt wilkie

4

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.


1
Ini tidak menjawab pertanyaan yang diajukan.
Dan C

1
Jika tidak ada file prj, ayam ogrinfo tidak akan mengembalikan hasil untuk informasi sistem koordinat.
Devdatta Tengshe

4

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.

  1. 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

  1. 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


Apa yang harus saya lakukan jika saya memiliki layer grid yang hanya menunjukkan sistem koordinat yang tidak diketahui dan saya tidak dapat menetapkan apa pun untuknya? Terima kasih.
AndrewLebron

2

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.


0

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


0

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. ..


0

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'}

Apakah Anda yakin itu akan berfungsi tanpa file .prj? bagaimana Fiona tahu?
JR

Tidak tanpa file .prj, itu tidak akan berfungsi
Sandip Nepal
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.