Saya membantu memastikan bahwa perangkat lunak sumber terbuka dapat menangani datum baru Australia dengan tepat, lihat situs web ICSM untuk perincian tentang proyek GDA2020.
Sekarang, QGIS sudah mendapatkan definisi GDA2020 termasuk, melalui GDAL, saya mengerti.
Contoh sistem referensi koordinat GDA2020 adalah ini:
+proj=utm +zone=55 +south +ellps=GRS80 +units=m +no_defs
Dan jika Anda melihat GDA94 CRS, didefinisikan seperti ini:
+proj=utm +zone=55 +south +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
Seperti yang Anda lihat, ini sangat mirip.
Sekarang, kedua CRS didefinisikan sama persis, tetapi, ada pergeseran koordinat dalam GDA94 ke GDA2020 sekitar 1,5 m ke arah timur laut. (Ada file grid shift dalam format NTv2 yang akan segera siap yang akan memungkinkan transformasi yang tepat, tapi bukan itu pertanyaannya.)
Tapi, jika Anda mengonversi antara GDA94 dan GDA2020 sekarang , menggunakan QGIS, tidak ada perubahan dalam koordinat. Ini pada dasarnya hanya label secara berbeda.
Haruskah transformasi 7 parameter sederhana diimplementasikan dalam Proj.4 atau alat sumber terbuka lainnya yang merupakan transformasi default (meskipun, tidak sempurna) antara GDA94 dan GDA2020?
Atau apakah hanya alat yang akan selalu melakukan perubahan?
Bagaimana ini harus ditangani?
(Dan saya hanya ingin mencatat lagi bahwa mengubah menggunakan kisi adalah ideal, dan itu ditangani dalam beberapa cara termasuk plugin QGIS ini .)
Jawaban:
Jika Anda mencari di database EPSG untuk
GDA94
CoordinateTransformation, Anda mendapatkan:EPSG:1150
GDA94 ke WGS84 (1) yang memiliki nilai semua-nolEPSG:8048
GDA94 ke GDA2020 (1) dengan 7 nilai yang diberikan oleh @ user30184Jadi adalah aman untuk mengambil mereka untuk GDA2020 ke WGS84 (mengurus tanda dan unit!) Sampai pergeseran grid baru diterbitkan. Itu akan mendapatkan nomor kode transformasi baru.
Saat ini terdapat kode Transformasi
EPSG:8049
ITRF2014 ke GDA2020 (1) yang menyatakan bahwa keduanya sama untuk saat ini, dengan nilai kenaikan tahunan. Jadi Anda bisa menggunakan kerangka waktu ITRF juga.sumber
Kamu bertanya:
FAQ di http://www.icsm.gov.au/gda2020/faq.html menginformasikan:
Setelah ICSM mempublikasikan 7 parameter, parameter transformasi kesamaan Anda dapat mulai menggunakannya sebagai
Saya sepertinya sudah diterbitkan di http://www.icsm.gov.au/gda2020/InterimReleaseNoteV1.0.pdf .
Anda dapat mencoba dengan parameter + towgs84 ini tetapi saya ingat bahwa Proj.4 mungkin menginginkan beberapa parameter dengan tanda terbalik.
Membuat tiket Proj.4 ketika parameter tersedia secara resmi dapat mempercepat proses dengan Proj.4 tetapi ketika database EPSG diperbarui dan Proj.4 mulai menggunakan database baru itu, perubahan dapat terjadi secara otomatis. Tergantung sedikit pada bagaimana GDA2020 akan diimplementasikan dalam database EPSG dan apakah algoritma baru diperlukan atau jika itu hanya pertanyaan tentang menambahkan parameter towgs84.
sumber
Masalahnya di sini adalah asumsi dalam pertanyaan bahwa proj.4 benar melaporkan dua CRS (sistem referensi koordinat) sebagai sama. Mereka tidak. String proj.4 yang dikutip bukan definisi CRS. Mereka dihasilkan dari definisi CRS, dan string proj.4 bukan gambaran lengkap. Definisi registri EPSG memberi kami informasi tambahan yang kami butuhkan untuk memahami apa yang sebenarnya terjadi.
Ini berasal dari pandangan dunia bahwa WGS-84 adalah datum global 'dan' proj historis.4 telah menggunakannya sebagai perantara ketika mengkonversi antara datum. Masalahnya, WGS-84 akan didefinisikan ulang setiap beberapa tahun ( kami berada di G1762 sekarang, disejajarkan dengan ITRF-08 ) karena diselaraskan kembali dengan perubahan dalam kerangka referensi ITRF, dari mana GDA juga diturunkan.
Ini telah menyebabkan pintasan dan asumsi ini dimasukkan ke dalam perilaku proj, meskipun dalam versi terbaru ini mulai berubah.
Melacak implikasi perubahan pada kerangka referensi, dan ketika mereka berubah bukanlah masalah besar sementara GPS konsumen memiliki ketelitian> 5 juta tetapi waktu berubah dan akurasi sub-meter membutuhkan alat yang memperhitungkannya dengan benar.
Jadi, untuk menjawab pertanyaan tersebut, kita perlu melacak kerangka data dan referensi apa yang menjadi dasar GDA94 dan GDA2020 CRS dan kemudian melihat transformasi apa yang tersedia.
EPSG:7844
GDA2020 2D Geografis CRS (Lat / Lon), dariEPSG:7843
GDA2020 3D Geographic CRS (L / L / H), dariEPSG:7842
GDA2020 3D Geocentric (ECEF X / Y / Z), yang semuanya digunakanEPSG: 1168 mendefinisikan kerangka referensi jangkar:
Jika Anda melakukan hal yang sama untuk GDA94 Anda akan melihat bahwa kerangka referensi adalah ITRF92, selaras pada 01/01/1994.
Jika mentransformasikan antara data ITRF02 / 14 dan GDA94 / GDA2020, datum diselaraskan, dan transformasi di antara mereka
null
hanya pada tanggal penyelarasan zaman. Pada dasarnya itulah yang dikatakan oleh string-string itu. Untuk kenyamanan, kami biasanya tidak suka harus terus-menerus mengubah koordinat yang kami simpan, jadi lebih mudah untuk melangkah mengubah pergeseran di antara mereka setiap beberapa tahun dan menerima tingkat kesalahan.Untuk sebagian besar aplikasi yang membutuhkan akurasi> 1m, itu cukup baik.
Tetapi kenyataan tidak melangkah berubah setiap beberapa tahun dan jika kita menginginkan transformasi yang lebih akurat, kita perlu mempertimbangkan penyimpangan jarak yang tergantung waktu dari datum sebelum / setelah penyelarasannya. Ini adalah transformasi 4D daripada 3D.
Transformasi antara GDA2020 dan WGS-84 atau ITRF2014 dijelaskan dalam:
Jika mentransformasikan antara GDA94 dan GDA2020, semuanya lebih sederhana karena kita hanya perlu mengetahui perbedaan antara kerangka referensi. Semacam. Ada lebih dari satu, dan yang benar untuk digunakan tergantung pada kapan, bagaimana dan di mana data dirujuk ke GDA94. Ini merupakan upaya untuk menghilangkan kesalahan karena metode yang kurang disempurnakan yang digunakan di tahun 90-an.
Ini adalah:
Untuk memahami dalam keadaan apa ini harus digunakan, baca manual teknis GDA2020
sumber
Membangun jawaban sebelumnya definisi proj4 terlihat seperti ini:
+proj=longlat +ellps=GRS80 +towgs84=-0.06155,0.01087,0.04019,-0.0394924,-0.0327221,-0.03289790,0.009994 +no_defs
Anda kemudian dapat menggunakan ini di salah satu zona grid proyeksi standar dengan hanya menambahkan parameter towgs84. misalnya
+proj=utm +zone=55 +south +ellps=GRS80 +towgs84=-0.06155,0.01087,0.04019,-0.0394924,-0.0327221,-0.03289790,0.009994 +units=m +no_defs
Untuk mendapatkan angka yang tepat dari bagian 3.1 dari spesifikasi, pertama-tama Anda membalikkan tanda parameter rotasi (seperti yang dibahas dalam bagian 2.2.1), tetapi kemudian membalikkan semuanya karena parameter dalam spesifikasi adalah transformasi dari WGS / GDA94 dan kami ingin transformasi ke WGS untuk definisi proj4. Jadi pada dasarnya semuanya kecuali rotasi dalam spesifikasi memiliki tanda terbalik.
Satu-satunya hal lain yang harus diperhatikan adalah untuk proj4 skalanya adalah parameter terakhir.
Purist akan menyarankan menggunakan pendekatan grid shift NtV2 tetapi file-file ini sangat besar dan saya telah menemukan bahwa di atas memberikan akurasi lebih dari 5cm menggunakan data sampel untuk Victoria. Saya juga menginginkan solusi yang akan bekerja dengan proj4js.
sumber