R에서 다각형 shapefile을 사용하여 래스터에서 래스터 추출

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

파일 래스터와 poly_shape는 두 기능에서 동일합니다. 여기서 무엇이 잘못 될 수 있는지 말해 줄 수 있습니까? 자르기 기능이 목록이 아닌 다른 래스터를 만드는 것이 옳습니까?

편집 : extent () 함수가 작동하지 않습니다. 여전히 같은 오류가 발생합니다. 그러나 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"

투영이 맞지 않는 것을 볼 수 있습니다. 추출 기능은 파일을 올바른 방식으로 자동 변환 할 수있는 것 같습니다. 나는 다음을했기 때문에 알고 있습니다.

• ArcGIS에서도 R로 추출한 다각형의 정확한 부분을 잘라냅니다
• 추출 된 R 다각형 (목록)의 모든 값의 합계를 계산했습니다.
• ArcGIS에서 잘라낸 모든 래스터 셀의 합계를 계산했습니다.

2는 동일한 결과를 가지므로 추출 기능이 올바르게 작동했다는 결론을 내릴 것입니다. 이제 두 가지 옵션이 있습니다.

1. 추출 된 목록에서 래스터를 다시 가져 오는 방법이 필요합니다.
2. 2 개의 데이터 세트 (래스터 + poly_shape)는 동일한 프 젝션을 사용해야하며 자르기 기능이 작동해야합니다

여기서 무엇을 제안 하시겠습니까?

추출 기능이 정확하게 동작합니다. 자르기 함수가 다각형의 범위를 사용하도록 강제 한 다음 객체를 마스크하여 다각형 영역을 나타내는 정확한 래스터를 반환 할 수 있습니다. 오류가 계속 발생하면 실제로 데이터가 겹치지 않음을 의미합니다.

R은 "즉석"투영을 수행하지 않으므로 투영을 확인하십시오. "extent ()"함수를 사용하여 익스텐트가 겹치는 지 확인할 수 있습니다.

다음은 다각형 범위를 사용하여 자르고 "래스터 화 된"다각형을 사용하여 결과 래스터를 마스킹하는 예입니다.

# 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와 다릅니다. 즉, 반대로 설정됩니다.