Lagi

Kepadatan data LiDAR, pengembalian per meter persegi


Bagaimana cara menghitung raster yang menunjukkan jumlah titik balik per meter persegi menggunakan SAGA GIS atau perangkat lunak sumber terbuka lainnya?

Tujuan saya adalah untuk mendapatkan topeng yang akan menunjukkan tempat (sel) di mana tidak ada hit dan akan membantu saya dalam menafsirkan data yang diturunkan dari LiDAR.


Saya mengembangkan GIS gratis dan open-source yang disebut Whitebox Geospatial Analysis Tools (dapat diunduh di sini) yang memiliki fungsionalitas analisis ekstensif untuk memproses data LiDAR. Whitebox berisi alat khusus untuk menghitung kepadatan titik file LiDAR LAS yang disebut Titik Kepadatan LiDAR.

Alat ini sangat spesifik untuk LiDAR, mengambil satu atau lebih file LAS sebagai input. Itu juga dapat menghitung kepadatan poin pengembalian pertama saja, poin pengembalian terakhir saja atau semua pengembalian. Selain itu, jika data LiDAR Anda berisi informasi klasifikasi, Anda dapat mengecualikan titik berdasarkan jenis klasifikasinya (misalnya bangunan atau titik vegetasi).

Ada alat kedua, yang disebut Hex-binning, yang juga akan melakukan analisis tipe kepadatan pada file LAS (ShapeFiles tipe titik juga), tetapi akan menampilkan kisi heksagonal vektor (ShapeFile). Ini bisa menjadi cara yang sangat efektif untuk memvisualisasikan cakupan titik kepadatan tinggi (lihat blog tentang hex-binning di sini).


Anda mungkin ingin menggunakan Fusion. Ini adalah perangkat lunak gratis untuk pemrosesan dan visualisasi point cloud.

1- Lihat ke manual untuk program perintah yang disebutKatalog. Ini mengembalikan statistik deskriptif dari titik cloud.

Yang Anda inginkan adalah sakelar Katalogkepadatan: luas, min, maks. Deskripsi manual mengatakan:

Membuat gambar untuk semua file data yang menunjukkan kepadatan kembali untuk area yang diwakili oleh setiap piksel. area adalah area piksel, min adalah kerapatan titik minimum yang dapat diterima per satuan luas, dan maks adalah batas atas untuk rentang kerapatan yang dapat diterima. Sel dengan densitas titik yang berada dalam kisaran min-max berwarna hijau, sel dengan densitas titik di bawah minimum berwarna merah, dan sel dengan densitas di atas maksimum berwarna biru.

Anda perlu mengatur min dan max yang sesuai sehingga akan menyoroti informasi yang Anda inginkan. Misalnya: jika Anda ingin melihat piksel dengan pengembalian kurang dari 1, setel min sama dengan 1. Piksel dengan pengembalian kurang dari satu per meter persegi akan muncul merah.

Baris perintah akan mirip dengan yang berikut, dengan mempertimbangkan piksel 5 meter (25 m²).

Katalog /density:25,1,5 pointcloud.las

Jika Anda memerlukan bantuan lebih lanjut untuk menjalankan baris perintah ini, di sini saya memberikan contoh terperinci tentangnya.

Ini adalahkatalog / kepadatanoutput, diambil dari satu sampel data yang disediakan oleh Fusion:


Gambar RGB dari pemandangan di atas dapat dilihat di sini.

2- Ada pilihan lain, lebih fleksibel dari Katalog.

Gunakan perintahGridMetrikuntuk mengekstrak jumlah pengembalian per piksel. Atur ukuran piksel dengan cara yang lebih cocok untuk analisis. Outputnya akan menjadi file nilai yang dipisahkan koma (csv) di mana setiap baris akan sesuai dengan piksel.

Contoh (dengan ukuran piksel sama dengan 5):

gridmetrics /noground /nointensity 2 5 4800K_gridmetrics lda_4800K_data.las

/tanpa tanahadalah sakelar yang digunakan bersama dengan kumpulan data yang dinormalisasi,/nointensitasadalah sakelar untuk tidak menghasilkan statistik intensitas LiDAR (akan mempercepat pemrosesan), dan2adalah jeda ketinggian untuk menghitung statistik penutup.

Kemudian, gunakan perintahCSV2Griduntuk mengubah file .csv menjadi file format raster ASCII (.asc), dan Anda akan dapat membukanya dalam perangkat lunak GIS yang berbeda.

Contoh (dimana variabel bunga 'total jumlah pengembalian' ada di kolom 5 dari sebelumnyagridmetricskeluaran):

csv2grid 4800K_gridmetrics_all_returns_elevation_stats.csv 5 4800K_return_per_pixel.asc


Untuk memeriksa jumlah titik per sel grid dari titik cloud LIDAR di GRASS, Anda dapat menggunakanr.in.xyzmodul. Modul ini membuat kisi dari awan titik menggunakan parameter "metode" yang dengannya Anda memilih cara menggabungkan titik saat mengimpor. Jika Anda memilihmetode = nmaka raster yang dihasilkan akan berisi jumlah poin untuk setiap sel. Dengan mengaturmetode = rata-rataAnda mendapatkan nilai rata-rata poin di setiap sel grid.

Detail lebih lanjut di halaman pemrosesan GRASS Wiki LIDAR.