R Shapefile和向量文件don'；绘制时，不能精确重叠

我想在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)

库（光栅）
图书馆（rgdal）
#下载自http://biogeo.ucdavis.edu/data/gadm2/shp/ITA_adm.zip
shape\u file=“ITA\u adm1.shp”
#下载自http://sedac.ciesin.columbia.edu/data/set/gpw-v3-population-density/data-download
#背景地理：国家，意大利；数据属性：网格
pop_density_file=“w001001.adf”
意大利地图这似乎是一个投影问题。您需要找出两个数据集的精确投影，并将其中一个数据集转换为另一个数据集的投影。
显然，
rgdal
包中有一个bug。我的问题在更新到版本0.8-12后得到解决
 我怎样才能找到他们的投影？如果我在QGIS中打开它们，它们会很好地重叠…也许在QGIS属性中有关于投影的说法？否则，您需要从数据供应商处获得该信息。从QGIS中的图层属性来看，它们似乎具有相同的CRS EPSG:4326-WGS 84。我还尝试导出它们，指定了CRS，但仍然没有重叠。它们对我来说完全重叠-请注意，如果绘制
italy\u map\u dens>0
以获得双色光栅，则更容易看到。我不知道为什么在这两种情况下都使用xlim和ylim进行绘图-使用xlim、ylim和add=TRUE没有意义。确保您拥有最新的rgdal，并告诉我们您使用了哪些光栅文件-数据集中有一些。人口密度数据需要注册和登录。你可以不用登录就下载。登录详细信息下面有一个按钮。