GRASS'ta görüntü sınıflandırması için Landsat görüntüleri nasıl keskinleştirilir?


11

Pan keskinleştirilmiş 15m çözünürlüklü görüntüler kullanarak GRASS'ta i.cluster> kullanarak Landsat sahnesinin denetimsiz bir sınıflandırmasını nasıl yapacağınızı öğrenmek istiyorum i.maxlik( Wiki'sinde verilen örnek normal 30m çözünürlükleri kullanıyor).

i.pansharpenÖnce pan-netleştirilmiş görüntüler oluşturmak için kullanmayı denedim , ancak i.pansharpenyalnızca d.rgbveya kullanılarak birleştirilebilen 3 kanal çıktısı verebilirim r.composite. Bildiğim kadarıyla, görüntü sınıflandırma işlemi 1 ila 7 arasında eksiksiz çok bantlı bantlar gerektirir. 15m çözünürlüklü bant 8 görüntüsünü kullanarak pan ile keskinleştirilen ayrı bantları (bant 1 ila bant 7) nasıl üretebilirim. görüntü sınıflandırma süreci?

Bunu nasıl yaptığını gösteren bir makale buldum ; temel olarak bir şekilde 30m-multispektral bantları 15m-pankromatik bantla birleştirmek için bir tür Temel Bileşen Analizi kullandılar. Kesin teklif:

"Yöntem ilk olarak 30 metrelik multispektral görüntüyü 15 metreye yeniden örnekliyor. Daha sonra multispektral görüntünün altı temel bileşenini de hesaplıyor (termal bant, kurucu çözünürlüğünden dolayı silindik). çözünürlük), 30 metrelik görüntünün ilk ana bileşeninin histogramıyla eşleşecek şekilde yeniden ölçeklendirilir ve ilk bileşen, yeniden ölçeklendirilmiş pankromatik bant ile değiştirilir.Bu, ilk ana bileşen genel parlaklığı, Pankromatik görüntünün geniş spektral bandı. Değiştirildikten sonra altı bileşen orijinal veri alanına geri dönüştürülerek uzamsal çözünürlükte bir iyileşme sağlanıyor. "

İlk olarak, kağıt hiçbir algoritma / formül göstermedi. Yukarıdaki alıntıyı karşılık gelen bir matematiksel formüle nasıl çevireceğim hakkında hiçbir fikrim yok. Her i.pansharpenzamanki Brovey veya IHS yerine PCA algoritması ile kullanabileceğimi fark ettim - ama yine de - çıktı sadece 3 kanal kırmızı, mavi ve yeşil olacak - ne yazık ki görüntü sınıflandırması için nasıl kullanılacağına dair hiçbir fikrim yok ..

Yani, elle yeni bir PCA algoritması yazmaya çalışırken kafamı kırmaya çalışmadan önce, herhangi biri pan keskinleştirilmiş Landsat görüntülerde görüntü sınıflandırması yapmanın daha kolay ve daha iyi bir yolunu göstermeme yardımcı olabilir mi? Yani - daha kolay bir yol olmalı, değil mi? Basit bir şeyi kaçırmış gibi hissediyorum.

Geriye kalan tek yol kendi senaryomu yazmaksa, lütfen beni orada yapmaya çalıştığım şeyin bir örneğine benzeyen bir şeye yönlendirebilir misiniz?

Herhangi bir yardım çok takdir!


İlk olarak, sunulan yöntemin 2. adımı (termal bant hariç tüm bantlarda PCA), yeniden örneklenmiş 15m bantlara veya orijinal 30m bantlara sahipse, sayfa 5'teki 1. paragrafı okuduktan sonra benim için net değil . 3. adımdan (referans olarak Pan-bandının histogram eşleşmesi 1. PC wihch uzamsal çözünürlüktedir ...?), Sanırım 2. adım (PCA) orijinal bantlara (30m) uygulandı. 4. adımda, gelişmiş Pan-bandı 1. PC'nin yerini alır - bu nedenle, muhtemelen 2. adım 15m-bantlarına uygulanır! - ve son olarak, geriye doğru bir PCA gelişmiş bir görüntü seti elde eder. Öyle mi?
Nikos Alexandris

İki kez PCA, biri orijinal altı bant (30m) seti için ve bir kez de yeniden örneklenmiş-15m bant seti için mi uygulanır?
Nikos Alexandris

Yanıtlar:


4

i.cluster için en az iki raster gerekir. Yani, i.pansharpen'den üç çıktı yeterli olacaktır.


Yani 3 çıktıyı da i.pansharpeniçine koymak zorunda i.groupmıyım? Hmm, bunu hiç düşünmemiştim ( i.groupayrı bantlarda raster gerektirdiğini düşündüm ) .. Sadece 3 raster kullanan sınıflamanın kalitesi, tüm aralığı kullanmaktan farklı olmalı, değil mi? Boş ver, önce deneyeceğim. TeşekkürlerVladimir Naumov!

Bunun hile yapması gerektiğini düşünüyorum! Basit bir şeyi kaçırdığımı biliyordum :)
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.