วิธีการ pansharpen LANDSAT 8 ใน GRASS?


9

เมื่อใช้ QGIS มาระยะหนึ่งในการทำแผนที่ฉันถูก 'บังคับ' ให้ขุดลงใน GRASS เนื่องจาก pansharpening ไม่สามารถใช้งานได้โดยตรงใน QGIS หลังจากทำความเข้าใจกับโครงสร้างมากขึ้นหรือน้อยลงฉันได้ทำการทดลองกับ i.fusion.brovey ฉันสามารถสร้างแรสเตอร์ที่มีความละเอียด 15 ม. แต่เมื่อโหลดลงใน QGIS ฉันไม่สามารถบรรลุ / ฟื้นความคมชัดของสีที่ดีดังที่ปรากฏในแรสเตอร์ 30 ม. ดั้งเดิม

ในโพสต์ก่อนหน้านี้มีคนยืนยันว่าและเสนอให้ใช้ฟังก์ชันอื่นสำหรับ landsat: 'i.his.rgb' (หลังจากใช้ i.landsat.toar ซึ่งฉันไม่ได้ใช้) เขาเริ่มแปลง min / max เป็น 1-255 ก่อนใช้ i.his.rgb ... rasters ดั้งเดิมของฉันมีช่วง 0-65535 มันไม่ได้ผลสำหรับฉัน ...

ฉันลาออกหลังจากผ่านไปหลายวันที่ต้องดิ้นรนเพื่อทำแพนเทอร์ฉันรู้สึกว่าติดขัด ... ทำไมการสมัคร i.pansharp ไม่ควรใช้กับ Landsat ตามที่คนอื่นแนะนำ? บางคนสามารถแนะนำเวิร์กโฟลว์ที่ตรงไปตรงมามากขึ้นเพื่อดำเนินการ pansharpening ของภาพ Landsat 8 หรือไม่? รอคำแนะนำบางอย่างขอบคุณล่วงหน้า ...

คำตอบ:


7

ลองดูที่ i.landsat.rgb - ทำการปรับสมดุลสีอัตโนมัติสำหรับภาพ LANDSAT ซึ่งอาจเป็นไปได้ก่อนที่จะเรียกใช้การแพนเพนชัน นอกจากนี้คุณยังอาจพิจารณาการแปลงตัวเลขดิจิตอลของแต่ละช่องเพื่อความกระจ่างใสด้านบนของบรรยากาศหรือการสะท้อนกับi.landsat.toar ดูเพิ่มเติมที่http://grasswiki.osgeo.org/wiki/LANDSAT

BTW: การมีช่วง 0-65535 สำหรับข้อมูล LANDSAT ดั้งเดิมนั้นแปลกเนื่องจากมีการส่งเป็น 8 บิต (0-255) ดังนั้นตรวจสอบข้อมูลอินพุตของคุณก่อน

อัปเดต: 10/2013

สำหรับขั้นตอนที่ง่ายขึ้นด้วย GRASS GIS 7 ปัจจุบันให้ดูหน้าป้อนคำอธิบายรูปภาพที่นี่การประมวลผลข้อมูล Landsat 8 ใน GRASS GIS 7: คอมโพสิต RGB และการลับคมกระทะ


ฉันมาเจอบทความที่คุณเชื่อมโยงจาก neteler.org แต่คำแนะนำไม่ได้ผลเลย (เพราะรูปแบบจุดลอยตัวของ rasters ที่สร้างด้วยi.landsat.toarดังนั้นฉันจึงสิ้นสุดคำแนะนำจากที่นี่: gis.stackexchange.com / a / 39239/9857
Juan

3
โปรดทราบว่าคำแนะนำนั้นได้รับมาหลายปีแล้วและไม่สะท้อนรหัส GRASS GIS ล่าสุด BTW: โมดูล pansharpening ที่ยอดเยี่ยมจริงๆคือ addon ใหม่: i.fusion.hpf - การผสมผสานความละเอียดสูงแบบ Panchromatic และความละเอียดต่ำของข้อมูล Multi-Spectral โดยใช้เทคนิคการเพิ่มตัวกรอง High-Pass ให้ดูgrass.osgeo.org/grass70/manuals/ addons / i.fusion.hpf.html
markusN

4

ขอบคุณ Michal และ Markus ในที่สุดฉันก็สามารถทำ pansharpening กับฟังก์ชั่นที่ระบุ brov เมื่อดำเนินการใน GRASS ฉันส่งออก RGB rasters ไปยัง gdal GTiff ด้วยประเภทข้อมูล Uint16 Markus ฉันได้ดาวน์โหลดและใช้ภาพ Landsat 8 หลายภาพ (tiff ของแต่ละวง 11) เพื่อให้ได้ภาพที่ดีทั้งในสีธรรมชาติ / สีเท็จ สำหรับความรู้ของฉัน Tiff ดั้งเดิมของ Landsat มาเป็น Uint16 ซึ่งมีค่าอยู่ระหว่าง 0 ถึง 65535 (หรือฉันขาดอะไรไป?) ... ขอบคุณสำหรับลิงก์ GRASS LANDSAT ...


3
ใช่รูปภาพ Landsat 8 ถูกจัดเก็บเป็น Uint16 ชนิด
webrian

3

การเหลา FIHS หรือการแปลง Brovery ค่อนข้างง่ายหากคุณทำตามสูตรจากบทความเช่น: http://www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

สมมติว่าคุณกำลังใช้ GRASS GIS และ pansharpening Landsat 8

ตั้งค่าภูมิภาคของคุณโดยใช้ pan 8 tif - สิ่งสำคัญคือต้องมีความละเอียด 15m ซึ่งจะถูกตั้งค่าหากคุณ

สำหรับ Brovery ที่ไม่มีตุ้มน้ำหนักให้ป้อนสิ่งนี้ใน Map Calculator สำหรับสีน้ำเงินเขียวและแดงตามลำดับ:

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

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

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

เรียกใช้ r.colors สำหรับแต่ละผลลัพธ์ข้างต้นและตั้งค่าการปรับสมดุลฮิสโตแกรมและตารางสีเป็นสีเทา

เรียกใช้ r.composite เลือกแถบที่เหมาะสม (แดง = 4, เขียว = 3, น้ำเงิน = 2)

ตรวจสอบภาพ

เห็นได้ชัดว่ากระบวนการนี้สามารถทำได้โดยพิจารณาเพิ่มเติมเกี่ยวกับการแก้ไขบรรยากาศการยืดสีและตุ้มน้ำหนักตามที่กล่าวไว้ในบทความด้านบน


2

คุณยังสามารถลองMonteverdiเพื่อดำเนินการ pansharpening เครื่องมือนี้สามารถใช้ GeoTIFF ปกติเป็นอินพุต (อย่างไรก็ตามต้องแยกรวมวง multispectral เป็น multiband raster เดียว) และทำงานได้ดีพอสมควร


2

หลังจากดิ้นรนกับปัญหานี้ฉันพบวิธีแก้ปัญหาง่ายๆที่ช่วยให้คุณสามารถปรับความคมชัดใน QGIS (2.18.21 เป็นรุ่นของฉัน) โดยใช้ SCP การปรับปรุงกับ SCP ผสมผสานวิธีการนี้จะอธิบายไว้ที่นี่https://fromgistors.blogspot.com/2015/07/major-update-semi-automatic-44.html

มันง่ายที่จะทำและมันก็วิ่งได้อย่างมีเสน่ห์! ฉันตรวจสอบผลลัพธ์สุดท้ายและพิกเซลวัดที่ความละเอียด 15 ม. - จับคู่กับวง 8 panchromatic หากใครอยากรู้อยากเห็น - ใช้วิธีการแปลง Brovey ตามที่อธิบายไว้ในคู่มือ (ดูhttps://media.readthedocs.org/pdf/semiautomaticclassificationmanual/latest/semiautomaticclassificationmanual.pdf )

ขั้นตอนที่ 1 สำหรับเครื่องมือ pansharpen ของ SCP ขั้นตอนที่ 3-4 สำหรับ SCP pansharpen


1

ยินดีต้อนรับสู่เว็บไซต์ Janos คำตอบที่ดีเช่นคุณควรมีคำอธิบายว่าทำไมพวกเขาถึงมีวิธีการที่ดี ..

สำหรับ pansharpening ของ Landsat 8 ภาพวิธีการง่ายขึ้นhttp://www.geosage.com/highview/download.html แนะนำของฉันลองมัน เพียงคลิกเดียว ....Spectraltransformer


-2

สำหรับ pansharpenig ฉันใช้ GUI สเปกตรัมหม้อแปลงสำหรับ Landsat 8 ภาพป้อนคำอธิบายรูปภาพที่นี่ http://www.geosage.com/highview/download.html ดีมาก สาขา BR Janos


2
ยินดีต้อนรับสู่เว็บไซต์ Janos คำตอบที่ดีเช่นคุณควรมีคำอธิบายว่าทำไมพวกเขาถึงมีวิธีการที่ดี ลิงก์คำตอบเท่านั้นมีแนวโน้มที่จะหมดอายุหากลิงก์นั้นหายไป
แอรอน
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.