r이있는 지정된 그리드가있는 세계지도를 계획하십시오.

Question

경도 / 위도 그리드가있는 세계지도를 플롯하기 위해 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 속성을 사용하는 방법을 모르겠습니다!

Answer 1

maps mapproject 패키지를 보았습니까?원하는 해상도로 해안선을 매우 빨리 그릴 수 있으며 투영 및 그리드가 함께 제공됩니다.Object 중 하나 (expand.grid(long,lat))에 투사를 적용 할 수 있습니다.Antartica의 예 :

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)

.