Menghasilkan kontur dari DEM menggunakan GDAL?

16

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=-9999

Saya 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?

— Nathan W
sumber

1

Untuk apa Anda menggunakan kontur? Apa statistik ketelitian posisi dan nilai dari DEM?

— whuber

DEM dihasilkan dari LiDAR jadi saya akan mengatakan mereka cukup bagus. Yah saya ingin satu set untuk dilihat publik untuk 1m dan 5m dan 250mm dengan di rumah untuk proposal desain.

— Nathan W

Untuk proposal desain dan pekerjaan in-house Anda lebih baik menggunakan DEM itu sendiri: dibandingkan dengan seperangkat polyline kontur, itu cocok untuk hampir semua perhitungan. Memang, untuk tontonan umum, Anda sering dapat melakukannya dengan baik dengan menumpangkan fitur-fitur utama (stream, ridges, puncak, dll.) Pada DEM yang berbukit. Apakah Anda benar-benar perlu membuat kontur sama sekali?

— whuber

Saya pribadi setuju dengan @whuber, tetapi saya mendapatkan pushback untuk menggunakan DEM karena perbedaan budaya antara pengguna akhir yang berbeda. Insinyur (setidaknya di sektor pasar saya) enggan menggunakan apa pun kecuali kontur. Selain itu, kontur dapat dihamparkan pada citra udara lebih mudah daripada DEM.

— DPierce

14

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 .

— Julien
sumber

1

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).

— pengguna23715
sumber

0

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

— christian.gobel
sumber