วิธีการปรับความคมชัดของภาพ Landsat สำหรับการจัดประเภทภาพใน GRASS


11

ฉันต้องการเรียนรู้วิธีการจัดหมวดหมู่ฉาก Landsat ที่ไม่มีการควบคุมโดยไม่ใช้i.cluster> i.maxlikใน GRASS โดยใช้ภาพความละเอียด 15 ม. ที่คมชัดขึ้น (ตัวอย่างที่ให้ไว้ในWikiใช้ความละเอียด 30m ทั่วไป)

ฉันพยายามใช้i.pansharpenในการสร้างภาพแพนรุนแรงขึ้นครั้งแรก แต่i.pansharpenสามารถส่งออกเพียง 3 ช่องทางที่สามารถรวมกันใช้หรือd.rgb r.compositeเท่าที่ฉันทราบกระบวนการจัดประเภทภาพต้องใช้คลื่นความถี่หลายหน่วยที่สมบูรณ์ 1 ถึง 7 ฉันจะสร้างแถบแยก (วงดนตรี 1 ถึงวงดนตรีที่ 7) ได้อย่างไรโดยใช้ภาพความละเอียด 15 เมตร 8 วงก่อนที่จะป้อนลงใน กระบวนการจำแนกภาพ?

ฉันพบกระดาษที่แสดงให้เห็นว่าพวกเขาทำสิ่งนี้ได้อย่างไร โดยพื้นฐานแล้วพวกเขาใช้การวิเคราะห์องค์ประกอบหลักเพื่อผสานวงดนตรี 30m-multispectral กับวง 15m-panchromatic คำพูดที่ถูกต้องจะเป็น:

"วิธีการนี้จะปรับภาพ multispectral ขนาด 30 เมตรเป็น 15 เมตรก่อนจากนั้นคำนวณส่วนประกอบหลักทั้งหกของรูปถ่ายแบบหลายส่วน ความละเอียด) ถูกลดขนาดให้ตรงกับฮิสโตแกรมขององค์ประกอบหลักตัวแรกของภาพขนาด 30 เมตรและองค์ประกอบแรกจะถูกแทนที่ด้วยแถบความยาวเปลี่ยนรูปแบบ rescaled นี่เป็นธรรมเนื่องจากองค์ประกอบหลักตัวแรกแทนความสว่างโดยรวมในลักษณะที่คล้ายกับ แถบสเปกตรัมกว้างของภาพ panchromatic หลังจากการเปลี่ยนชิ้นส่วนประกอบทั้งหกจะถูกเปลี่ยนกลับไปเป็นพื้นที่ข้อมูลดั้งเดิมเพื่อปรับปรุงความละเอียดเชิงพื้นที่ "

ก่อนอื่นกระดาษไม่พบอัลกอริทึม / สูตรใด ๆ ฉันไม่รู้ว่าจะเปลี่ยนคำพูดข้างต้นให้เป็นสูตรทางคณิตศาสตร์ที่สอดคล้องกันได้อย่างไร ฉันรู้ว่าฉันสามารถใช้i.pansharpenกับอัลกอริทึม PCA แทน Brovey หรือ IHS ปกติ - แต่ยัง - เอาท์พุทจะมีเพียง 3 ช่องทางของสีแดงสีน้ำเงินและเขียว - ซึ่งฉันเศร้าไม่มีเงื่อนงำเกี่ยวกับวิธีการใช้สำหรับการจำแนกภาพ ..

ดังนั้นก่อนที่ฉันจะพยายามแตกหัวพยายามเขียนอัลกอริธึม PCA ใหม่ด้วยตนเองใครช่วยฉันชี้ไปที่วิธีที่ง่ายและดีกว่าในการจำแนกภาพบนภาพ Landsat ที่คมชัดขึ้นหรือไม่? ฉันหมายถึง - ควรมีวิธีที่ง่ายกว่านี้ใช่ไหม ฉันรู้สึกว่าฉันขาดอะไรง่าย ๆ ไป

หากเส้นทางเดียวที่เหลือคือการเขียนสคริปต์ของตัวเองพวกคุณช่วยชี้ให้ฉันดูบางสิ่งที่นั่นคล้ายกับตัวอย่างของสิ่งที่ฉันพยายามจะทำหรือไม่?

ความช่วยเหลือใด ๆ ที่ชื่นชมมาก!


ก่อนอื่นมันไม่ชัดเจนสำหรับฉันหลังจากอ่านย่อหน้าที่ 1 ในหน้า 5 หากขั้นตอนที่ 2 ของวิธีการนำเสนอ (PCA สำหรับทุกแบนด์ยกเว้นเทอร์มอลความร้อน) ใช้resampledเป็นวง 15m หรือ 30m ดั้งเดิม จากขั้นตอนที่ 3 (การจับคู่ฮิสโตแกรมของ Pan-band ที่ใช้เป็นข้อมูลอ้างอิงสำหรับพีซีเครื่องแรกที่มีความละเอียดเชิงพื้นที่ ... ?) ฉันเดาว่าขั้นตอนที่สอง (PCA) ถูกนำไปใช้กับวงดนตรีดั้งเดิม (30m) ในขั้นตอนที่ 4 Pan-band ที่ได้รับการปรับปรุงจะแทนที่พีซีเครื่องที่ 1 ดังนั้นอาจเป็นขั้นตอนที่ 2 ที่ใช้กับแบนด์ 15m! - และสุดท้าย PCA หลังจะได้ชุดรูปภาพที่ได้รับการปรับปรุง มันเป็นอย่างนั้นเหรอ?
Nikos Alexandris

พวกเขาใช้ PCA สองครั้งครั้งเดียวสำหรับชุดเดิมของวงดนตรีหกวง (30 ม.) และอีกครั้งสำหรับวงดนตรีที่มีความยาวตั้งแต่ 15 เมตรถึงหนึ่งครั้ง?
Nikos Alexandris

คำตอบ:


4

i.cluster ต้องการแรสเตอร์อย่างน้อยสองตัว ดังนั้นสามเอาต์พุตจาก i.pansharpen จะเพียงพอ


คุณหมายถึงฉันจะต้องใส่เอาท์พุททั้ง 3 จากi.pansharpenในi.group? อืมไม่เคยคิดเลยว่า (ฉันคิดว่าi.groupต้องใช้แรสเตอร์เป็นวงแยกต่างหาก) .. คุณภาพของการจัดหมวดหมู่ที่ใช้แรสเตอร์เพียง 3 ตัวควรแตกต่างจากการใช้ช่วงทั้งหมดใช่ไหม? ไม่เป็นไรฉันจะลองก่อน ขอบคุณ @Vladimir Naumov!

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