既存の投影グリッド (137x108 グリッド ポイント) にマップをプロットする良い方法を見つけようとしています。ライブラリ フィールドと as.image-image.plot は非常に簡単で分かりやすいことがわかりましたが、表示されているボックスの解像度には完全に満足していません。変更しようとしましnrowたncolが、臨界値 (100+) を超えると、マップに白い隙間が表示されます。以下に例を示しnetcdf4ます。ファイルを読み取るにはライブラリが必要です。
アドバイスや代替アプローチは大歓迎です。
#Elevation data can be downloaded in the following link,
#http://www.ecad.eu/download/ensembles/data/Grid_0.44deg_rot/elev_0.44deg_rot_v9.0.nc.gz
#INSTALLING-LOADING NECESSARY LIBRARIES
install.packages(fields)
install.packages(ncdf)
library(fields)
library(ncdf)
#READING FILES
orog.nc <- open.ncdf("elev_0.44deg_rot_v9.0.nc")
orog <- get.var.ncdf( orog.nc, "elevation")
orog[orog <= -99]=NA
x <- get.var.ncdf( orog.nc, "Actual_longitude")
y <- get.var.ncdf( orog.nc, "Actual_latitude")
xy <- data.frame(lon=c(x),lat=c(y))
set.panel(1,3)
#PLOTING LOW RES
look <- as.image(c(orog), x=xy, na.rm=T)
image.plot(look, col=terrain.colors(100))
world( add=TRUE)
#PLOTING "HIGH RES"
look <- as.image(c(orog), x=xy, na.rm=T, nrow=100, ncol=100)
image.plot(look, col=terrain.colors(100))
world( add=TRUE)
#PLOTING TOO "HIGH RES"
look <- as.image(c(orog), x=xy, na.rm=T, nrow=150, ncol=150)
image.plot(look, col=terrain.colors(100))
world( add=TRUE)
アップデート
tim riffe の建設的なコメントを使用すると、実際にはデータの解像度がより細かくなり (0.22deg - 272x214 など)、白いグリッド ポイントの効果に問題を起こすことなく、コードでより多くの nrow と ncol を設定できます。nrow-ncol しきい値はどのように定義されるのでしょうか?
#Install Libraries and read data as the previous post,
#but instead download 0.22deg file for fine resolution.
#http://www.ecad.eu/download/ensembles/data/Grid_0.22deg_rot/elev_0.22deg_rot_v9.0.nc.gz
set.panel(1,2)
#PLOTING "LOW RES"
look <- as.image(c(orog), x=xy, na.rm=T)
image.plot(look, col=terrain.colors(100))
world( add=TRUE)
#PLOTING "HIGH RES"
look <- as.image(c(orog), x=xy, na.rm=T, nrow=175, ncol=175)
image.plot(look, col=terrain.colors(100))
world( add=TRUE)