Saya memiliki kelas fitur (kabupaten Carolina Selatan, jadi area geografis yang cukup besar) di NAD83 SC State Plane. Itu perlu ditransformasikan ke proyeksi kedua (NAD83 UTM 17), kemudian diubah kembali ke aslinya. Saya akan menggunakan alat Proyek Esri untuk mencapai ini.
Dapatkah transformasi ganda ini menyebabkan pergeseran lokasi koordinat poligon, dan seberapa jauh - sentimeter, meter, kilometer?
Jawaban:
Saya tidak tahu mesin proyeksi mana yang digunakan ArcGis, tetapi pertanyaan yang sangat menarik juga untuk proj.4. Jadi saya mencobanya untuk menguji mesin proyeksi proj.4 dalam lingkungan GNU-R. Saya menggunakan sudut NAD 83 - UTM 17 dan EPSG 26917 dan memproyeksikan ulang 10000 dan 1000000 kali secara berulang dan menghitung perbedaannya dengan nilai awal.
Inilah hasilnya:
Tampaknya kesalahan "proyeksi ulang" berada dalam kisaran sentimeter untuk 10.000 putaran.
Dan berkembang menjadi kesalahan dalam rentang meter jika Anda menjalankan loop 1000000 kali.
Ini skripnya.
Tes lebih lanjut dalam lingkungan statistik harus mudah. Penjelasan skrip dan kode untuk lingkungan linux tersedia di github.com/bigopensky .
sumber
WGS84 Bounds: -84.0000, 24.0000, -78.0000, 83.0000
adalah wilayah yang tepat untuk diminati. Apakah saya melakukan kesalahan?Esri memiliki mesin proyeksi sendiri.
Sebagian besar metode proyeksi dan transformasi geografis / datum berperilaku baik ketika digunakan dalam bidang yang sesuai minat. Jika Anda terlalu jauh di luar zona UTM, Mercator melintang tidak selalu 'terbalik' (konversikan ke lintang-bujur). Proyeksi yang digunakan untuk seluruh dunia mungkin memiliki beberapa masalah di atau sekitar kutub atau +/- 180 meridian atau 'anti-meridian' (meridian yang berlawanan dengan pusat sistem referensi koordinat yang diproyeksikan).
Saya berlari 4 poin yang jatuh di luar Carolina Selatan melalui mesin proyeksi Esri. Untuk stress test 1k atau 10k atau 1M poin, saya harus memberi kode sesuatu karena tes serupa saya yang ada hanya melakukan 'round-trip' - diproyeksikan ke geografis untuk diproyeksikan. 32133 adalah NAD 1983 State Plane South Carolina (meter). 26917 adalah NAD 1983 UTM zona 17 Utara.
Jadi Anda dapat melihat kami memiliki dua poin yang kembali pada 10e-09.
Penanganan di ArcGIS dipersulit oleh fakta bahwa ada referensi spasial. Referensi spasial mencakup sistem koordinat ditambah beberapa nilai penyimpanan dan analisis. Secara default, sistem koordinat yang menggunakan meter disimpan dengan presisi sepersepuluh milimeter, 0,0001.
Pengungkapan: Saya bekerja untuk Esri.
sumber
Saya pikir ini adalah kasus di mana Anda perlu menguji alur kerja yang diusulkan terhadap beberapa fitur titik uji, yang mudah untuk menambahkan bidang koordinat XY.
Bandingkan nilai-nilai XY poin awal Anda dengan yang telah Anda proyeksikan / ubah (betapapun berkali-kali), dan Anda akan menghitung perbedaannya.
sumber