"추출"기능은 어떻게 다른 예측을 처리합니까?

동일한 크기의 사각형 격자 셀을 사용하여 래스터에서 가중 평균 추출을 수행하려면 extract() 함수를 사용해야합니다. 내 폴리곤 그리드는 UTM21n이고 래스터는 GCS WGS84 데이텀 D입니다. 래스터를 extract()에 사용하기 전에 재 투영해야합니까? 아니면 함수가 제대로 처리할까요?"추출"기능은 어떻게 다른 예측을 처리합니까?

출처

2013-10-11 Herman Toothrot

SpatialPolygonshere에 대해 extract의 소스 코드를 찾을 수 있습니다.

setMethod('extract', signature(x='Raster', y='SpatialPolygons'), 
function(x, y, fun=NULL, na.rm=FALSE, weights=FALSE, cellnumbers=FALSE, small=FALSE, df=FALSE, layer, nl, factors=FALSE, sp=FALSE, ...){ 

    px <- projection(x, asText=FALSE) 
    comp <- .compareCRS(px, projection(y), unknown=TRUE) 
    if (!comp) { 
     .requireRgdal() 
     warning('Transforming SpatialPolygons to the CRS of the Raster') 
     y <- spTransform(y, px) 
    } 
...

extract는이에 설명되지 않는다는 사실에도 불구하고, 사실 투사 자체 (래스터의 돌출부에 SpatialPolygon의 투사를 변경)을 수행에서하는 것을 제안 : 코드는 다음 코드로 시작 도움말 페이지.

출처

2013-10-11 07:32:12 plannapus

+1 annoyingly 나는이 기능이 문서화되어 있지 않기 때문에 정확한 상황을 처리하기 위해 추출기 주위에 래퍼 함수를 작성했습니다. 알았어, 고마워! –

흥미 롭습니다. 프로젝션이 폴리곤 레이어에 설정되어 있지 않으면 어떻게되는지 궁금합니다. –

설명서에는 자동 재 투영에 대한 언급이 없습니다. 그래서, 나는 함수가 이것을하지 않는다고 가정하는 것이 바람직하다고 생각합니다. 따라서 extract에 전화하기 전에 자신을 재 투척해야합니다.

출처

2013-10-11 06:49:19

"추출"기능은 어떻게 다른 예측을 처리합니까?

답변

관련 문제