ฉันกำลังใช้ QGIS ฉันต้องการคลิปเลเยอร์แรสเตอร์แรสเตอร์โดยใช้เลเยอร์ขอบเขตผู้ดูแลระบบที่เป็นข้อมูลเวกเตอร์ อย่างไรก็ตามเครื่องมือการประมวลผลทางภูมิศาสตร์ดูเหมือนจะใช้งานได้กับข้อมูลเวกเตอร์เท่านั้น
ฉันจะคลิปชั้นฝนนี้ได้อย่างไร
ฉันกำลังใช้ QGIS ฉันต้องการคลิปเลเยอร์แรสเตอร์แรสเตอร์โดยใช้เลเยอร์ขอบเขตผู้ดูแลระบบที่เป็นข้อมูลเวกเตอร์ อย่างไรก็ตามเครื่องมือการประมวลผลทางภูมิศาสตร์ดูเหมือนจะใช้งานได้กับข้อมูลเวกเตอร์เท่านั้น
ฉันจะคลิปชั้นฝนนี้ได้อย่างไร
คำตอบ:
ติดตั้งปลั๊กอิน GDAL จากนั้นใช้เครื่องมือ Clipper
หากคุณมีความสนใจใช้งูใหญ่เป็นเอกสารที่ดีที่มีอยู่ใน GeospatialPython.com, ที่นี่
และแหล่งที่มา clipraster.py เป็นที่นี่
กระบวนการ:
การตัดแรสเตอร์เป็นชุดของการคลิกปุ่มอย่างง่าย ๆ ในแพ็คเกจซอฟท์แวร์เชิงพื้นที่ระดับสูง ในแง่ของการประมวลผลภาพเชิงพื้นที่นั้นจริงๆแล้วมีขนาดใหญ่มากและมีอาเรย์หลายมิติ การสำรวจระยะไกลที่ง่ายที่สุดคือการดำเนินการทางคณิตศาสตร์ในอาร์เรย์เหล่านี้เพื่อดึงข้อมูลจากข้อมูล เบื้องหลังนี่คือสิ่งที่ซอฟต์แวร์กำลังทำอยู่ (ให้หรือทำตามขั้นตอน):
- แปลง vectorfilefile เป็นเมทริกซ์ซึ่งสามารถใช้เป็นหน้ากาก
- โหลดภาพเชิงพื้นที่ลงในเมทริกซ์
- โยนเซลล์รูปภาพใด ๆ ออกนอกขอบเขตรูปร่าง
- ตั้งค่าทั้งหมดนอกขอบเขต shapefile เป็นค่า NODATA (null)
- ตัวเลือก: ดำเนินการยืดฮิสโตแกรมบนภาพเพื่อสร้างภาพที่ดีขึ้น
- บันทึกภาพผลลัพธ์เป็นแรสเตอร์ใหม่
clipraster.py
ลิงค์เสีย พิจารณาอัปเดต Tks
ฉันขอแนะนำให้ใช้gdalwarp
เพราะคุณสามารถเพิ่มประสิทธิภาพของกระบวนการ
ตัวอย่างเช่นหากคุณต้องการตัดแรสเตอร์inraster.tif
ด้วย Shapefile extent.shp
:
gdalwarp -cutline extent.shp -crop_to_cutline -of GTiff -dstnodata 255 inraster.tif inraster_cropped.tif -co COMPRESS=LZW -co TILED=YES --config GDAL_CACHEMAX 2048 -multi
พารามิเตอร์สองตัวสุดท้ายช่วยให้คุณสามารถเร่งกระบวนการโดยใช้:
ฉันไม่แน่ใจว่าจะทำอย่างไรใน QGIS อย่างไรก็ตามคุณสามารถใช้ซอฟต์แวร์ GIS โอเพ่นซอร์สอื่นเช่น GRASS หรือ GDALWarp เพื่อตัดภาพแรสเตอร์ของคุณ ดูคำตอบ keymirror และ geographika ด้านล่าง
สำหรับโซลูชันที่ใช้ ArcGIS Desktop คุณอาจต้องการดูการตัดแรสเตอร์โดยใช้รูปหลายเหลี่ยมหรือรูปหลายเหลี่ยม