Файл Shapefile и векторный файл точно не перекрываются при построении графика

Question

Я хочу нарисовать шейп-файл над растровым файлом в R, но я не могу сделать их идеально наложенными: растровое изображение поворачивается на несколько градусов против часовой стрелки. Это проблема с проекцией?

Пожалуйста, рассмотрим следующую 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) 

Answer 1

Видимо, в пакете rgdal произошла ошибка. Моя проблема была решена после того, как я обновил ее до версии 0.8-12.

Answer 2

Это, как представляется, проблема с проекцией. Вам нужно узнать точную проекцию обоих наборов данных и преобразовать один из наборов данных в проекцию другого набора данных.