Bagaimana cara melakukan klasifikasi deret waktu LANDSAT tanpa pengawasan menggunakan Rumput?

8

Saya menggunakan GRASS untuk melakukan klasifikasi tanpa pengawasan pada serangkaian waktu gambar LANDSAT (res temporal bulanan, lebih dari 30 tahun), untuk membandingkan perubahan tutupan lahan di area tertentu.

Saya khawatir jika saya menjalankan i.clusteruntuk setiap peta dalam seri waktu tanda tangan yang digunakan untuk kelas mungkin sedikit berbeda untuk setiap peta, yang memengaruhi validitas perbandingan. Apakah ini masalahnya? Jika demikian, apakah ada cara untuk memastikan bahwa tanda tangan kelas yang digunakan pada setiap peta dalam deret waktu identik di seluruh?

Ide awal saya adalah menjalankan i.clusterpada peta tertentu di mana semua kelas terwakili dengan baik, kemudian gunakan sigfile yang dihasilkan dari ini sebagai input sigfile i.maxlikuntuk semua peta.

grass remote-sensing classification landsat si_2012
sumber
1
Salah satu cara untuk mendapatkan tanda tangan yang umum untuk koleksi grid terkait adalah memindahkan semuanya ke lokasi yang berbeda dan membuat mosaik mereka (dengan beberapa pita nilai nol di antara mereka semua): membuat tanda tangan untuk itu. Ini bekerja di sekitar masalah potensial yang salah satu gambar mungkin tidak menampilkan berbagai kelas penuh.
whuber
1
Terima kasih @whuber. Ini terdengar seperti ide yang bagus, namun, area yang saya lihat terdiri dari 11 ubin LANDSAT yang dirubah menjadi satu. Membubarkan ini dan merapikannya bersama-sama akan menghasilkan file besar dan memakan waktu lama. Sebaliknya, dapatkah saya membuat grup gambar dari semua titik dalam rangkaian waktu dan menghasilkan tanda tangan untuk ini?
si_2012
1
Kumpulan data Anda tampaknya sangat besar dan mungkin itu akan menjadi tugas yang sulit untuk memperbaiki atmosfer ditambah radiometrik menormalkan gambar Anda. Yang terakhir, jika berhasil, mungkin mendorong peningkatan yang signifikan terhadap perbandingan silang radiometrik.
Nikos Alexandris
Saran bagus @ Nikos - Saya akan menulis skrip menggunakan paket R 'landsat' untuk melakukan identifikasi PIF untuk normalisasi relatif ubin. Saya akan memperbarui posting ini ketika saya mendapat hasil.
si_2012
1
Saya tidak menyarankan bahwa pencocokan histogram ( i.histo.match in grass 7 ) sesuai dengan kebutuhan Anda, tetapi Anda mungkin ingin melihatnya juga.
Nikos Alexandris

Jawaban: