แยก Raster จาก Raster โดยใช้ Polygon shapefile ใน R

ฉันใหม่สำหรับ R และใช้แพ็คเกจแรสเตอร์ ฉันมีปัญหาในการแยกรูปหลายเหลี่ยมจากไฟล์แรสเตอร์ที่มีอยู่ ถ้าฉันใช้

extract(raster, poly_shape)

ฟังก์ชั่นบนแรสเตอร์มันมักจะสร้างรายการที่มีข้อมูล สิ่งที่ฉันต้องการคือการแยกไฟล์แรสเตอร์อื่นที่ฉันสามารถโหลดด้วย ArcGIS อีกครั้ง หลังจากอ่านอีกเล็กน้อยฉันคิดว่าฟังก์ชั่นการครอบตัดคือสิ่งที่ฉันต้องการจริงๆ แต่เมื่อลองใช้ฟังก์ชั่นนี้

crop(raster, poly_shape)

ฉันได้รับข้อผิดพลาดนี้:

Error in .local(x, y, ...) : extents do not overlap
In addition: Warning message:
In intersect(extent(x), extent(y)) : Objects do not overlap

ไฟล์ raster และ poly_shape เหมือนกันทั้งสองฟังก์ชั่น คุณช่วยบอกฉันทีว่าอะไรที่นี่ผิด มันถูกต้องไหมที่ฟังก์ชั่นการครอบตัดสร้างอีกแบบแรสเตอร์และไม่ใช่รายการ?

แก้ไข : ฟังก์ชั่นขอบเขต () ไม่ทำงานสำหรับฉัน ฉันยังคงได้รับข้อผิดพลาดเดียวกัน แต่ฉันแน่ใจว่าชุดข้อมูล 2 ชุดทับซ้อนกัน! กับ

extract(raster, poly_shape)

ฉันได้รับข้อมูลที่ถูกต้องจากมัน เช่นเดียวกับรายการและไม่ใช่แรสเตอร์อย่างที่ฉันต้องการ ฉันเพิ่งโหลดชุดข้อมูลใน ArcGIS มาก่อนและพวกมันก็พอดีมากดังนั้นฉันจึงไม่ได้ตรวจสอบการฉาย ตอนนี้ฉันพยายามแล้ว

projection(raster) # "+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs"
projection(poly_shape) # "+proj=utm +zone=32 +ellps=GRS80 +units=m +no_defs"

และคุณจะเห็นได้ว่าการคาดการณ์ไม่เหมาะสม ฟังก์ชั่นแตกไฟล์ดูเหมือนจะสามารถแปลงไฟล์โดยอัตโนมัติในวิธีที่ถูกต้อง ฉันรู้ว่าเพราะฉันทำต่อไปนี้:

• ฉันตัดส่วนที่แน่นอนของรูปหลายเหลี่ยมที่ฉันแยกใน R ออกมาใน ArcGIS ด้วย
• ฉันคำนวณผลรวมของค่าทั้งหมดของรูปหลายเหลี่ยม R ที่แยกแล้ว (รายการ)
• ฉันคำนวณผลรวมของเซลล์แรสเตอร์ทั้งหมดที่ฉันตัดออกใน ArcGIS

2 มีผลลัพธ์เดียวกันแน่นอนดังนั้นฉันเดาว่าข้อสรุปควรเป็นว่าฟังก์ชันการแยกทำงานได้ถูกต้อง ตอนนี้ฉันมี 2 ตัวเลือกฉันเดา:

1. ฉันต้องการวิธีที่จะนำ Raster ออกจากรายการที่แยกอีกครั้งหรือ
2. ชุดข้อมูล 2 ชุด (raster + poly_shape) จำเป็นต้องใช้การปฏิเสธแบบเดียวกันและฟังก์ชันการครอบตัดควรทำงาน

คุณอยากแนะนำให้ทำอะไรที่นี่?

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

โปรดจำไว้ว่า R ไม่สามารถฉายภาพได้ทันทีดังนั้นโปรดตรวจสอบการคาดการณ์ของคุณ คุณสามารถตรวจสอบว่าขอบเขตของคุณทับซ้อนกันโดยใช้ฟังก์ชัน "scope ()"

นี่คือตัวอย่างของการปลูกพืชโดยใช้ขอบเขตรูปหลายเหลี่ยมจากนั้นปิดบังแรสเตอร์ที่เกิดขึ้นโดยใช้รูปหลายเหลี่ยมที่ "rasterized"

# Add required packages
require(raster)
require(rgeos)
require(sp)

# Create some data using meuse 
data(meuse)
  coordinates(meuse) <- ~x+y
    proj4string(meuse) <- CRS("+init=epsg:28992")
data(meuse.grid)
  coordinates(meuse.grid) = ~x+y
    proj4string(meuse.grid) <- CRS("+init=epsg:28992")
      gridded(meuse.grid) = TRUE    
        r <- raster(meuse.grid) 
          r[] <- runif(ncell(r))

# Create a polygon
f <- gBuffer(meuse[10,], byid=FALSE, id=NULL, width=250, 
                         joinStyle="ROUND", quadsegs=10)   

# Plot full raster and polygon                       
plot(r)
  plot(f,add=T)

# Crop using extent, rasterize polygon and finally, create poly-raster
#          **** This is the code that you are after ****  
cr <- crop(r, extent(f), snap="out")                    
  fr <- rasterize(f, cr)   
    lr <- mask(x=cr, mask=fr)

# Plot results
plot(lr)
  plot(f,add=T)

สิ่งที่ฉันค้นหาจริงๆคือmask()ฟังก์ชั่น

mask(raster, poly_shape)

ทำงานได้โดยไม่มีข้อผิดพลาดและส่งคืนสิ่งที่ฉันค้นหา

ขอบเขตทำงานได้ดี ... ฉันคิดว่าขอบเขต Xmin, Xmax, Ymin และ Ymax ของคุณนั้นแตกต่างจาก X และ Y ของแรสเตอร์ของคุณ - เช่นตั้งอยู่ตรงกันข้าม