กำลังคำนวณขอบเขตภาพ / รอยเท้าของภาพถ่ายดาวเทียมโดยใช้เครื่องมือโอเพนซอร์สหรือไม่?

ฉันจำเป็นต้องสร้างรูปหลายเหลี่ยมของภาพแรสเตอร์เลเยอร์เดียวหลายรูปไม่ใช่กล่องขอบเขต / ขอบ แต่พื้นที่ที่ไม่มีค่า Nodata ดังแสดงที่นี่: การสร้างรูปร่างไฟล์แสดงรอยเท้าของ Rasters? .

ในการตอบคำถามข้างต้นปลั๊กอิน Image Boundaryได้ถูกกล่าวถึง แต่ฉันไม่พบมันใน QGIS 1.8.0 Lisboa ที่ทำงานบน Ubuntu

เครื่องมือยังคงมีอยู่หรือไม่?

ถ้าไม่มีมีวิธีการทำเช่นนี้โดยใช้เครื่องมือ opensource ไม่ว่าจะเป็น gdal, R, QGIS, GRASS หรือคล้ายกันไม่ใช่กับ ArcMap?

คุณสามารถเพิ่มที่เก็บเก่าที่มีปลั๊กอินในที่เก็บอย่างเป็นทางการ

http://pyqgis.org/repo/contributed

ปลั๊กอินขอบเขตภาพมีตัวเลือกสำหรับการคำนวณพิกเซลที่ถูกต้อง แต่ภาพต้องเป็นฉากเต็มรูปแบบวิธี CBERS หรือ Landsat ที่กระบวนการนี้คำนวณมุมทั้งสี่แรก ภาพจำเป็นต้องมีค่าตัวเลขด้วยค่าศูนย์สำหรับพื้นที่ที่ไม่มีภาพ (สแกนโดยเซ็นเซอร์รับสัญญาณดาวเทียม)

ผู้แต่ง Image Boundary

ในที่สุดฉันก็พบวิธีที่จะทำสิ่งนี้:

ขั้นตอนที่ 1: gdalwarp -dstnodata 0 -dstalpha -of GTiff foo1 foo2

สิ่งนี้ทำสองสิ่งสำคัญ: กำหนดปลายทางไม่มีข้อมูล (นอกขอบเขต) เป็น 0 และสร้างแถบอัลฟา

ขั้นตอนที่ 2: gdal_polygonize.py foo2 -b 2 -f "ESRI Shapefile" foo3

ขั้นตอนที่สองใช้แถบอัลฟา (แบนด์ 2) สร้างในขั้นตอนที่ 1 และสร้างไฟล์รูปร่างจากวงดนตรีนั้น

สิ่งนี้สามารถสคริปต์ได้อย่างง่ายดายในสคริปต์ทุบตีหากคุณมีภาพจำนวนมากเพื่อสร้างโครงร่างที่แน่นอนสำหรับ  

ฉันใช้gdal_translateตามที่แนะนำโดยโครงการ GDAL

gdal_translate -b mask -of vrt -a_nodata 0 test.tif test.vrt
# Note the  -a_nodata 0 doesn't seem to work when the mask is input, so do another pass
gdal_translate -b 1 -of vrt -a_nodata 0 test.vrt test2.vrt
gdal_polygonize.py -q  -8 test2.vrt -b 1 -f "ESRI Shapefile" testdata.shp

Image Boundary pluginไม่ทำงานสำหรับฉันดังนั้นฉันจึงใช้วิธีการเดียวกันกับ GDAL อย่างไรก็ตามมันใช้งานได้สำหรับฉันหลังจากเปลี่ยนขั้นตอนแรกเป็น:

ขั้นตอนที่ 1 : gdalwarp -srcnodata 0 -dstalpha -of GTiff foo1 foo2

ฉันกำลังทำงานกับวง Landsat8 (ที่ไม่มีข้อมูล = 0) และเมื่อใช้-dstnodataฟังก์ชั่นที่ฉันได้รับ:

band1กับไม่มีข้อมูล = 'ไม่มีข้อมูล'
band2 (อัลฟาวงดนตรี) = '255' สำหรับทั้งฉาก / ขอบเขต

ในขณะที่-srcnodataฟังก์ชั่นฉันได้รับ:

band1กับไม่มีข้อมูล = 'ไม่มีข้อมูล'
band2 (alpha วงดนตรี) ไม่มีข้อมูล = 'ไม่มีข้อมูลและข้อมูลที่ถูกต้อง = '255' ซึ่งก็จะช่วยให้สารสกัดจากรูปหลายเหลี่ยมสำหรับข้อมูลที่ถูกต้องเท่านั้น

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

ขอบเขตของภาพอยู่ล้าสมัยแล้ว คุณสามารถใช้ปลั๊กอิน Footprint http://pt.slideshare.net/LuizMotta3/qgis-ibama-imagefootprint