Kami baru saja menerima sejumlah besar DEM di tempat kerja dan saya ingin menghasilkan kontur dari mereka. DEM memiliki resolusi 1m dan ukuran 1kmx1km.
Output dari gdalinfo:
Driver: AAIGrid/Arc/Info ASCII Grid
Files: 380000_6888000_1k_1m_DEM_ESRI.asc
Size is 1000, 1000
Coordinate System is `'
Origin = (380000.000000000000000,6888000.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Corner Coordinates:
Upper Left ( 380000.000, 6888000.000)
Lower Left ( 380000.000, 6887000.000)
Upper Right ( 381000.000, 6888000.000)
Lower Right ( 381000.000, 6887000.000)
Center ( 380500.000, 6887500.000)
Band 1 Block=1000x1 Type=Float32, ColorInterp=Undefined
NoData Value=-9999Saya tahu saya bisa menggunakan gdal_contour untuk menghasilkan kontur ( posting blog saya pada topik ) tapi saya bertanya-tanya apa beberapa praktik terbaik untuk menghasilkan kontur.
Apakah ada aturan yang harus Anda ikuti untuk mendapatkan hasil maksimal dari kontur tetapi tidak membuat barang jadi atau kehilangan terlalu banyak informasi?
Katakanlah saya ingin menghasilkan tiga set kontur:
- 250mm
- 1 m
- 5 m
Apakah ada yang harus saya lakukan pada DEM sebelum setiap set?
Apakah memposting garis adalah cara yang baik untuk pergi atau memuluskan raster adalah pilihan yang lebih baik?
Jawaban:
Aturan kartografi untuk mewakili lega karena kontur disajikan dalam buku terkenal Imhof pada representasi bantuan, bab C . Beberapa aturan ini diberikan pada halaman wikipedia ini . Rekomendasi utama saat menyederhanakan kontur untuk mempertahankan karakteristik utama medan.
Menghaluskan kontur secara independen tidak mencegahnya tumpang tindih: Disarankan untuk menghaluskan DTM terlebih dahulu. Perataan Gaussian tradisional dengan parameter yang sesuai tergantung pada resolusi target memungkinkan untuk menghapus detail kecil. Kekurangannya adalah smoothing DTM mengisi lembah dan depresi, dan meratakan pinggiran dan puncak. Menggunakan algoritma filter douglass-peucker seperti dalam makalah ini dapat menjadi solusi. Ada juga sejumlah metode berdasarkan penggunaan algoritma pemulusan kontur yang dibatasi oleh jaringan drainase atau kerangka. Akhirnya, untuk mencegah kontur tumpang tindih di bagian yang ceroboh, dimungkinkan untuk menghapusnya secara lokal atau juga merusaknya .
sumber
Saya ingin komentar kedua @ whuber. Analisis Kuantitatif selalu lebih baik dari DEM secara langsung dan Analisis Visual sering lebih baik bila dilakukan dari Hillshade daripada kontur.
Untuk menjawab pertanyaan secara langsung: Di ArcGIS saya akan menggunakan Focal Statistics atau Aggregate [Spatial Analyst Toolbox] untuk memperlancar garis kontur yang dihasilkan. Karena kontur adalah fitur analisis visual, jumlah smoothing akan bervariasi sesuai kebutuhan Anda. Jadi, Anda harus bereksperimen dan melihat mana yang paling cocok untuk proyek Anda.
Menghaluskan garis setelah membuatnya tidak bekerja tetapi sedikit kikuk dibandingkan dengan memodifikasi raster terlebih dahulu. Satu pembersihan pasca kontur yang sering saya lakukan adalah memilih garis dengan panjang tertentu (mis. Panjang <10 ') dan menghapusnya. Ini rids data "bising" sedikit dari kontur tertutup (yaitu loop kecil) yang tidak mungkin untuk mencerminkan karakter sebenarnya dari permukaan yang dimodelkan oleh kontur.
Hal-hal lain yang mungkin Anda lihat di [pengguna ArcGIS] adalah:
Membuat Mosaik Raster dan / atau Statistik Fokus Mosaik Raster dan menggunakan model untuk menghasilkan kontur untuk seluruh dataset.
Membuat kontur 3D untuk penggunaan AutoCAD.
Waspadai area data DEM yang tidak biasa (mis. Bentangan luas tanah dataran rendah yang membutuhkan interval kontur yang berbeda untuk secara akurat mewakili permukaan, area vegetasi lebat yang memberikan nilai Bare Earth palsu, area relief vertikal - tebing, dll).
sumber
Ada cara mudah menggunakan gdal_contour. Setelah mengatur semua opsi di jendela dialog Anda kemudian dapat mengedit baris perintah dan sebagai gantinya "-i interval" Anda dapat menggunakan level tetap "-fl level". Seperti gambar di bawah. Anda dapat memeriksa opsi lain di sini http://www.gdal.org/gdal_contour.html
sumber