5
เพิ่มความเร็วของการครอบตัดหน้ากากและแยกแรสเตอร์โดยรูปหลายเหลี่ยมใน R
ฉันกำลังแยกพื้นที่และร้อยละของการใช้ที่ดินประเภทต่าง ๆ จากแรสเตอร์ตามขอบเขตหลายเหลี่ยมหลายพัน ฉันพบว่าฟังก์ชั่นการดึงข้อมูลทำงานได้เร็วขึ้นมากถ้าฉันวนซ้ำรูปหลายเหลี่ยมและครอบตัดแต่ละอันจากนั้นปิดบังแรสเตอร์ให้มีขนาดเท่ากับรูปหลายเหลี่ยมที่เฉพาะเจาะจง อย่างไรก็ตามมันค่อนข้างช้าและฉันสงสัยว่าใครมีคำแนะนำในการปรับปรุงประสิทธิภาพและความเร็วของรหัสของฉันหรือไม่ สิ่งเดียวที่ฉันได้พบที่เกี่ยวข้องกับการนี้คือการตอบสนองนี้โดยโรเจอร์ Bivand ที่แนะนำให้ใช้GDAL.open()และGDAL.close()เช่นเดียวกับและgetRasterTable() getRasterData()ฉันดูสิ่งเหล่านี้ แต่เคยมีปัญหากับ gdal ในอดีตและไม่รู้จักดีพอที่จะรู้วิธีนำไปใช้ ตัวอย่างที่ทำซ้ำได้: library(maptools) ## For wrld_simpl library(raster) ## Example SpatialPolygonsDataFrame data(wrld_simpl) #polygon of world countries bound <- wrld_simpl[1:25,] #name it this to subset to 25 countries and because my loop is set up with that variable ## Example RasterLayer …