Saya mencari metode untuk memproses gambar penginderaan jauh dan mengekstrak area mahkota pohon individu dari gambar.
Saya memiliki citra areal panjang gelombang visual, dan data Lidar dari daerah tersebut. Lokasi yang dimaksud adalah daerah gurun, jadi tutupan pohonnya tidak sepadat hutan. Resolusi citra udara adalah 0,5 kaki kali 0,5 kaki. Resolusi Lidar adalah sekitar 1 x 1 kaki. Baik data visual dan lidar berasal dari dataset Pima County, Arizona. Contoh tipe citra udara yang saya miliki ada di akhir tulisan ini.
Pertanyaan ini deteksi Single Tree di ArcMap? tampaknya menjadi masalah yang sama, tetapi tampaknya tidak ada jawaban yang bagus di sana.
Saya dapat memperoleh klasifikasi yang wajar dari tipe-tipe vegetasi (dan informasi tentang persentase tutupan keseluruhan) di area tersebut dengan menggunakan klasifikasi Iso Cluster di Arcmap, tetapi ini memberikan sedikit informasi tentang masing-masing pohon. Yang paling dekat saya dengan apa yang saya inginkan adalah hasil melewati keluaran klasifikasi isocluster melalui fitur Raster ke Polygon di Arcmap. Masalahnya adalah bahwa metode ini bergabung dekat oleh pohon menjadi satu poligon.
Sunting: Saya mungkin harus memasukkan beberapa detail tentang apa yang saya miliki. Kumpulan data mentah yang saya miliki adalah:
- Data las lengkap, dan taster raster dihasilkan darinya.
- Citra visual (seperti contoh gambar yang ditampilkan, tetapi mencakup area yang jauh lebih luas)
- Pengukuran langsung secara manual terhadap sebagian pohon di daerah tersebut
Dari ini saya telah menghasilkan:
- Klasifikasi tanah / vegetasi.
- DEM / DSM raster.
sumber
Jawaban:
Ada banyak literatur tentang deteksi mahkota individu dalam data spektral dan lidar. Metode bijak, mungkin mulai dengan:
Falkowski, MJ, AMS Smith, PE Gessler, AT Hudak, LA Vierling dan JS Evans. (2008). Pengaruh tutupan kanopi hutan konifer pada akurasi dua algoritma pengukuran pohon individu menggunakan data Lidar. Canadian Journal of Remote Sensing 34 (2): 338-350.
Smith AMS, EK Strand, CM Steele, DB Hann, SR Garrity, MJ Falkowski, JS Evans (2008) Produksi peta struktur-spasial vegetasi dengan analisis per objek perambahan juniper dalam foto udara multi-temporal. Canadian Journal Remote Sensing 34 (2): 268-285
Jika Anda tertarik pada metode Wavelet (Smith et al., 2008), saya memilikinya dengan kode Python tetapi, ini sangat lambat. Jika Anda memiliki pengalaman Matlab, ini adalah tempat penerapannya dalam mode produksi. Kami memiliki dua makalah di mana kami mengidentifikasi ~ 6 juta hektar perambahan juniper di Oregon bagian timur menggunakan metode wavelet dengan citra NAIP RGB-NIR sehingga, terbukti dengan baik.
Baruch-Mordo, S., JS Evans, J. Severson, JD Naugle, J. Kiesecker, J. Maestas, dan MJ Falkowski (2013) Menghemat belibis bijak dari pohon: Solusi proaktif untuk mengurangi ancaman utama terhadap kandidat spesies Konservasi Biologis 167: 233-241
Poznanovic, AJ, MJ Falkowski, AL Maclean, dan JS Evans (2014) Penilaian Akurasi Algoritma Deteksi Pohon di Juniper Woodlands. Rekayasa Fotogrametri & Penginderaan Jauh 80 (5): 627–637
Ada beberapa pendekatan yang menarik, dalam dekomposisi objek umum, dari literatur ruang matematika keadaan terapan menggunakan proses Gaussian multiresolusi untuk mendekomposisi karakteristik objek lintas skala. Saya menggunakan tipe-tipe model ini untuk menggambarkan proses multi-skala dalam model ekologis tetapi bisa diadaptasi untuk mendekomposisi karakteristik objek gambar. Menyenangkan, tetapi sedikit esoteris.
Gramacy, RB, dan HKH Lee (2008) Bayesian mengamati model proses Gaussian dengan aplikasi untuk pemodelan komputer. Jurnal Asosiasi Statistik Amerika, 103 (483): 1119–1130
Kim, HM, BK Mallick, dan CC Holmes (2005) Menganalisis data spasial non-stasioner menggunakan proses Gaussian piecewise. Jurnal Asosiasi Statistik Amerika, 100 (470): 653-668
sumber
Untuk membuat DHM, kurangi DEM dari DEM, ini bisa dilakukan di Esri Raster Calculator atau GDAL_CALC . Ini akan menempatkan semua ketinggian Anda di 'bidang permainan level'.
Sintaks (Pengganti jalur lengkap untuk DEM, DSM & DHM):
DHM sebagian besar adalah 0 (atau cukup dekat), yang Anda nilai dari nodata Anda. Dengan Raster Calculator atau GDAL_CALC Anda dapat mengekstraksi nilai lebih dari nilai arbitrer berdasarkan jumlah kebisingan yang Anda amati dalam DHM. Tujuan dari ini adalah untuk mengurangi kebisingan dan menyoroti hanya mahkota vegetasi - dalam contoh di mana dua 'pohon' berdekatan ini harus dibagi menjadi dua gumpalan yang berbeda.
Sintaks (Pengganti jalur penuh untuk Binary & DHM dan nilai yang diamati untuk Nilai):
Sekarang dengan baik GDAL_CALC atau Esri IsNull membuat raster biner, yang dapat dipoligonisasi dengan GDAL_Polygonize atau Esri Raster to Polygon .
Untuk memperbaiki poligon menghapus poligon terlalu kecil dan kemudian membandingkannya dengan band RGB mencari tanda tangan, di Esri alat Statistik Zonal akan membantu. Kemudian Anda dapat membuang poligon yang jelas tidak memiliki statistik yang benar (berdasarkan eksperimen dan data Anda, saya tidak bisa memberi Anda nilai).
Ini akan membuat Anda mencapai akurasi sekitar 80% dalam merencanakan setiap mahkota.
sumber
eCognition adalah perangkat lunak terbaik untuk itu, saya melakukan itu dengan menggunakan perangkat lunak lain tetapi eCognition lebih baik. Berikut ini referensi literatur tentang subjek:
misalnya http://www.mdpi.com/1424-8220/14/12/22643
Selain itu:
mis. http://www.nrcresearchpress.com/doi/abs/10.1139/x05-030#.VJmMb14gAA
sumber
Saya memiliki masalah yang sama beberapa tahun yang lalu. Saya punya solusi yang tidak memerlukan data LAS yang difilter atau data tambahan lainnya. Jika Anda memiliki akses ke data LiDAR, dan dapat menghasilkan DEMs / DSMs / DHMs (DEM selanjutnya, saya tidak memperdebatkan semantik nomenklatur model permukaan) dari hasil yang berbeda, skrip berikut mungkin berguna.
Script arcpy menelan 3 DEM dan mengeluarkan poligon hutan dan titik pohon. 3 DEM harus memiliki resolusi spasial yang sama (yaitu 1 meter) dan luasan, dan mewakili pengembalian pertama, pengembalian terakhir, dan tanah kosong. Saya memiliki parameter yang sangat spesifik untuk ekstraksi sayuran, tetapi parameternya dapat diubah sesuai dengan kebutuhan lainnya. Saya yakin prosesnya dapat ditingkatkan, karena ini adalah upaya serius pertama saya pada skrip python.
sumber
Saya memposting ini sebagai jawaban karena batas panjang dalam komentar, tidak ada harapan untuk kredit :). Kuas yang sangat luas, asalkan Anda punya DEM.
Jumlah maksimum grup dalam proses = jumlah pohon di dalam setiap poligon. Kriteria tambahan, misalnya jarak antara 'pohon' di dalam poligon mungkin membantu ... Perataan DEM menggunakan kernel juga merupakan opsi.
sumber