2013-11-21 3 views
3

R 파일의 래스터 파일 위에 모양 파일을 플롯하려고하지만 완전히 겹치게 할 수 없습니다. 래스터가 반 시계 방향으로 몇도 회전 한 것처럼 보입니다. 투영에 문제가 있습니까? 형상 파일과 벡터 파일이 플롯 할 때 정확하게 겹치지 않음

enter image description here

는 다음 MWE 투사에 문제가 나타나지 않습니다

library(raster) 
library(rgdal) 

# Download from http://biogeo.ucdavis.edu/data/gadm2/shp/ITA_adm.zip 
shape_file = "ITA_adm1.shp" 
# Download from http://sedac.ciesin.columbia.edu/data/set/gpw-v3-population-density/data-download 
# Setting Geography: Country, Italy; Data Attributes: Grid 
pop_density_file ="w001001.adf" 

italy_map <- readOGR(dsn = shape_file, layer = "ITA_adm1") 
italy_map_dens <- raster(pop_density_file) 

colPal <- colorRampPalette(c("white", "red"))(500) 

par(mar=c(0,0,0,0)) 
plot(italy_map_dens, xlim = c(6.70, 18.32), ylim = c(35.2, 47.6), axes=FALSE, box=FALSE, legend=FALSE, col=colPal) 
plot(italy_map, xlim = c(6.70, 18.32), ylim = c(35.2, 47.6), border="grey", add=TRUE) 
+0

인구 밀도 데이터는 등록과 로그인이 필요합니다. 응. – Spacedman

+0

로그인하지 않고 다운로드 할 수 있습니다. 로그인 세부 정보 아래에 버튼이 있습니다. – CptNemo

답변

2

분명히 rgdal 패키지에 버그가있었습니다. 내 문제는 버전 0.8-12로 업데이트 한 후 해결되었습니다.

1

을 고려하시기 바랍니다. 두 데이터 세트의 정확한 투영법을 찾고 데이터 세트 중 하나를 다른 데이터 세트의 투영법으로 변환해야합니다.

+0

투영법을 어떻게 알 수 있습니까? QGIS에서 열면 잘 겹칩니다 ... – CptNemo

+0

아마도 QGis 속성에서 투영과 관련하여 어떤 것이 있습니까? 그렇지 않으면 데이터 공급 업체로부터 해당 정보를 가져와야합니다. –

+0

QGIS의 레이어 속성에서 동일한 CRS EPSG : 4326 - WGS 84를 갖고있는 것처럼 보입니다. 또한 CRS를 지정하여 내보내려고했지만 여전히 겹치지 않았습니다. – CptNemo