3

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)
4

2 に答える 2

2

どうやらrgdalパッケージにバグがあったようです。バージョン 0.8-12 に更新したところ、問題は解決しました。

于 2013-11-21T20:42:50.807 に答える
1

投影に問題があるようです。両方のデータセットの正確な射影を見つけて、データセットの 1 つを他のデータセットの射影に変換する必要があります。

于 2013-11-21T07:44:26.623 に答える