Saya memiliki sejumlah besar data dalam bentuk CSV yang berisi koordinat UTM sebagai X
dan Y
dan nilai ketinggian sebagai Z
informasi. Saya perlu mengubah data ini menjadi DEM sebagai GeoTIFF untuk analisis lebih lanjut. Dalam hal ini, jumlah yang sangat besar berarti 16 m. garis, dengan satu titik X
, Y
dan Z
per baris. Poin didistribusikan secara merata, oleh karena itu interpolasi tidak diperlukan; setiap titik hanya perlu dikonversi menjadi sel raster.
Data asli datang tanpa pemisah, dengan lebar kolom tetap. Saya sudah menemukan cara mengubah sintaks file untuk menggunakan pemisah alih-alih lebar tetap dan menghilangkan semua karakter spasi, menggunakan sed editor teks stream . Dari sini, biasanya alur kerja saya akan mengimpor data ke ArcGIS dengan membuat kelas fitur dari X
, Y
dan Z
data dan sebagai langkah kedua, mengubah titik shapefile menjadi GeoTIFF, menggunakan alat Point to Raster . Namun, file yang saya miliki saat ini terlalu besar untuk proses ini.
Alih-alih alur kerja yang dijelaskan di atas, saya mencari alternatif yang efisien dan menemukan GDAL. Namun, dalam gdal_translate
, format terdekat yang didukung yang dapat saya temukan dalam daftar tipe file yang didukung, adalah kisi ASCII tetapi tidak ada XYZ yang dipisah koma. Kesulitan lain adalah, bahwa saya memiliki koordinat UTM , sementara sebagian besar contoh tampaknya menggunakan koordinat derajat desimal. Namun, saya harus tetap berada dalam sistem UTM (atau setidaknya, GeoTIFF keluaran saya harus dalam sistem koordinat UTM).
Jadi saya mencari cara untuk mengubah CSV XYZ menjadi GeoTIFF, menggunakan GDAL , tetapi sejauh ini tidak dapat menemukan contoh yang berhubungan dengan masalah yang tepat ini. Saya akan sangat senang untuk beberapa petunjuk atau bahkan contoh kode.
Jawaban:
Anda dapat melakukan ini menggunakan GDAL, secara langsung mendukung format XYZ . Tidak masalah jika koordinat Anda adalah UTM, gdal_translate akan ditampilkan dalam sistem koordinat yang sama.
Jadi untuk mengkonversi ke GeoTIFF semudah:
Lihatlah dokumen GeoTIFF untuk opsi keluaran (seperti kompresi) dan dokumen gdal_translate untuk info penggunaan lebih lanjut. Secara khusus, Anda harus menentukan apa sistem koordinat dengan
-a_srs
parameter.Dipisahkan koma / spasi dan lebar kolom tetap, dengan dan tanpa baris tajuk didukung.
Satu-satunya gotcha yang saya ketahui adalah:
File harus diurutkan dengan benar (oleh Y, lalu X).
sumber
-a_srs EPSG:12345
Data Terbuka DGM200 Jerman: https://gdz.bkg.bund.de/index.php/default/open-data/digitales-gelandemodell-gitterweite-200-m-dgm200.html
Unduh file xyz: https://daten.gdz.bkg.bund.de/produkte/dgm/dgm200/aktuell/dgm200.utm32s.xyzascii.zip
Konversikan xyz ke file GeoTIFF:
sumber