0

経度/緯度グリッドで世界地図をプロットするために、maptools パッケージの wrld_simpl を使用しようとしています。

たとえば、経度と緯度を含む netCDF ファイルがあり、グリッド上に配置できるすべてのポイントを含むマトリックスを作成します。

lat <- get.var.ncdf(nc,"lat")
long <- get.var.ncdf(nc,"lon")
pts <- SpatialPoints(expand.grid(long,lat), CRS(proj4string(wrld_simpl)))

ここで、座標がグリッド上にある国をプロットしたいと思います。

wrld_simpl 属性の使い方がわかりません!

4

1 に答える 1

0

maps mapprojectパッケージを見ましたか？必要な解像度で海岸線を非常に迅速に描画でき、投影とグリッドが付属しています。オブジェクトの 1 つにプロジェクションを適用できます ( expand.grid(long,lat))。南極の例:

library("mapproj")
m <- map( "world", "Antarctica", plot=FALSE)
map("world", "Antarctica", proj="conic", param=-90, fill=TRUE, col="gray")
map.grid(m, nx=4, ny=6, col="black")

front.lim <- data.frame(x=seq(-180, 180, length.out=1000),
                        y=rep(-62, 1000))
front.lim <- mapproject(front.lim$x, front.lim$y, projection="conic", parameters=-90)
lines(front.lim$x, front.lim$y, col="red", lwd=2)

既に ncdf 行列をプロットしている場合は、単純に を追加できmap(add=TRUE)ます。x と y の範囲を読み取り、この領域に境界線を描画します。

plot(c(-10, 45), c(30, 60), type="n", xlab="", ylab="") # image(ncdf)
map(add=TRUE, fill=TRUE, col="gray")
abline(v=seq(-10, 45, 7.5), lty=2)
abline(h=seq(-30, 60, 7.5), lty=2)
于 2013-11-18T20:11:06.787 に答える