LANDSAT 8 GRASS'ta pansharpen nasıl?

Haritalama için QGIS'i bir süre kullandıktan sonra, pansharpening QGIS'de doğrudan uygulanamayacağı için GRASS'a girmeye 'zorlandım'. Yapıyı az çok anladıktan sonra, i.fusion.brovey ile bazı denemeler yaptım. 15m çözünürlükte raster üretebildim, ancak QGIS'e yüklerken orijinal 30m rasterlerde olduğu gibi güzel renk kontrastını elde edemiyorum / geri alamıyorum.

Daha önceki bir gönderide birisi bunu doğruladı ve landsat için başka bir işlev kullanmayı önerdi: 'i.his.rgb' (kullanmadığım i.landsat.toar kullandıktan sonra). İ.his.rgb'yi kullanmadan önce min / max'ı 1-255'e dönüştürmeye başlar ... Orijinal rasterlerimin 0-65535 aralığı vardır. Benim için işe yaramadı ...

Sadece bir pansharp yapmak için uğraştıktan birkaç gün sonra çıktım, sıkıştım ... Neden sadece i.pansharp uygulamak başkalarının önerdiği gibi Landsat için çalışmıyor? Birisi Landsat 8 görüntülerini pansharpening yapmak için daha basit bir iş akışı önerebilir mi? Biraz tavsiye bekliyorum, şimdiden teşekkürler ...

Bir göz atın i.landsat.rgb - gerçekleştirir otomatik dengeleme muhtemelen pansharpening çalıştırmadan önce, LANDSAT görüntüler için renk. Ayrıca, kanalların dijital sayılarını i.landsat.toar ile atmosfer üstü radyasyona veya yansıtma oranına dönüştürmeyi de düşünebilirsiniz . Ayrıca bkz http://grasswiki.osgeo.org/wiki/LANDSAT

BTW: orijinal LANDSAT verileri için 0-65535 aralığında olması tuhaftır, çünkü 8 bit kanal (0-255) olarak iletilirler. Bu yüzden önce giriş verilerinizi kontrol edin.

Güncelleme: 10/2013

Mevcut GRASS GIS 7 ile basitleştirilmiş bir prosedür için GRASS GIS 7: RGB kompozitleri ve pan keskinleştirme içindeki Landsat 8 verilerini işleme sayfasına bakın

Teşekkürler Michal ve Markus. Sonunda belirtilen fonksiyon brov ile pansharpening yapabildim. GRASS'ta gerçekleştirildikten sonra, RGB rasterlerini Uint16 veri türü olarak gdal GTiff'e aktardım. Markus, birkaç Landsat 8 görüntüsünü (11 bandın her birinin tiff'i) indirip kullanıyorum, hem doğal / yanlış renkte güzel görsel sonuçlar elde ediyorum. Bildiğim kadarıyla, orijinal landsat tiffleri Uint16 olarak 0 ila 65535 arasında değişen değerlerle geliyor (ya da bir şey eksik miyim?) ... GRASS LANDSAT bağlantısı için teşekkürler ...

Aşağıdaki gibi bir makaledeki formülleri takip ederseniz, pan keskinleştirici FIHS veya Brovery Dönüşümü oldukça basittir: http://www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

GRASS GIS ve pansharpening Landsat 8 kullandığınızı varsayarsak.

Bölgenizi pan 8 tif kullanarak ayarlayın - bunu yaparsanız ayarlanacak 15m çözünürlüğe sahip olmak önemlidir.

Ağırlıksız Brovery için bunu sırasıyla Mavi, yeşil ve kırmızı için Harita Hesaplayıcıya girin:

band8 / (band2 + band3 + band4) / 3 * band2

band8 / (band2 + band3 + band4) / 3 * band3

band8 / (band2 + band3 + band4) / 3 * band4

Yukarıdaki çıktıların her biri için r.colors komutunu çalıştırın ve histogram eşitleme ve renk tablosunu gri olarak ayarlayın.

R.composite komutunu çalıştırın, uygun bantları seçin (kırmızı = 4, yeşil = 3, mavi = 2)

Resmi kontrol et.

Açıkçası bu işlem, atmosferik düzeltme, renk gerilmesi ve bant ağırlıkları için yukarıdaki makalede tartışıldığı gibi daha fazla düşünülerek yapılabilir.

Pansharpening yapmak için Monteverdi'yi de deneyebilirsiniz . Bu araç, giriş olarak normal GeoTIFF'leri tüketebilir (ancak, ayrı çok bantlı bantların bir çok bantlı rasterle birleştirilmesi gerekir) ve makul şekilde iyi performans gösterir.

Bu sorunla mücadele ettikten sonra, SCP'yi kullanarak QGIS'de (2.18.21 benim sürümüm) keskinleştirmenizi sağlayan basit bir çözüm buldum. Bu yöntemi içeren SCP güncellemesi burada açıklanmaktadır https://fromgistors.blogspot.com/2015/07/major-update-semi-automatic-44.html .

Bu yapmak kolaydır ve bir cazibe gibi koştu! Son çıktıyı kontrol ettim ve pikseller pankromatik 8 bandıyla eşleşerek 15 m çözünürlükte ölçülüyor. Merak eden biri varsa - kılavuzda açıklandığı gibi Brovey Transform yöntemini kullanır (bkz. Https://media.readthedocs.org/pdf/semiautomaticclassificationmanual/latest/semiautomaticclassificationmanual.pdf ).

Janos sitesine hoş geldiniz. İdeal olarak sizinki gibi cevapların neden iyi yöntemler olduklarına dair bir açıklaması olmalıdır. ..

Landsat 8 resimlerinin taranması için daha kolay yöntem http://www.geosage.com/highview/download.html . Benim tavsiyem, dene. Sadece bir tıklama ....

Pansharpenig için Landsat 8 resimleri için GUI spektral transformatörü kullandım http://www.geosage.com/highview/download.html Çok iyi. BR Janos