Saya memiliki cloud point yang koordinatnya berkaitan dengan sistem koordinat lokal. Saya juga memiliki titik kontrol tanah dengan nilai GPS. Bisakah saya mengonversi koordinat lokal ini ke sistem koordinat global menggunakan PROJ.4 atau perpustakaan lain?
Setiap kode dalam Python untuk masalah yang disebutkan di atas akan sangat membantu.
python
coordinate-system
gps
proj
ground-control
pengguna18953
sumber
sumber
Jawaban:
Anda tampaknya ingin melakukan transformasi affine antara sistem koordinat lokal Anda dan sistem koordinat georeferensi.
Affine mentransformasikan secara mendasar semua sistem koordinat dan dapat diwakili oleh persamaan matriks di bawah ini.
Namun, Anda memiliki masalah dua langkah.
Proj.4 unggul di # 2: mentransfer antara sistem koordinat georeferensi dengan matriks transformasi yang dikenal. Setahu saya tidak bisa digunakan untuk menemukan matriks transformasi dari data titik. Namun, Anda dapat melakukan semuanya dengan mudah dengan menggunakan aljabar linier ringan (inversi matriks kuadrat-terkecil) di Numpy. Saya telah menggunakan versi kelas ini untuk mengurangi data dari beberapa studi lapangan:
Ini dapat digunakan sebagai berikut:
projected_coordinates
sekarang di WGS84, UTM, atau sistem koordinat apa pun yang direkam oleh GPS Anda. Fitur utama dari metode ini adalah dapat digunakan dengan sejumlah titik pengikat (3 atau lebih) dan mendapatkan akurasi semakin banyak titik pengikat digunakan. Anda pada dasarnya menemukan yang paling cocok melalui semua poin dasi Anda.sumber
Fungsi transformasi rumput melakukan persis apa yang Anda butuhkan, meskipun tidak berbasis python atau proj seperti yang diminta:
http://grass.osgeo.org/grass65/manuals/g.transform.html
sumber
Selalu lebih mudah untuk mengidentifikasi sistem koordinat lokal, seperti yang kami lakukan di sini:
Proyeksi stereografik WGS84 ellipsoid di pesawat [python]
GDAL sekarang dapat mengubah data vektor menggunakan poin GCP.
sumber
Saya terjebak dalam masalah yang sama beberapa minggu yang lalu, saya menemukan skrip python yang dapat membantu. Solusi asli dari sini
sumber