R中的削波光栅
我正在绘制美国东北部的地图。 地图背景需要是海拔图或年平均温度图。我有两个来自Worldclim.org的光栅,它们给了我这些变量,但我需要将它们裁剪到我感兴趣的状态。关于如何做到这一点的任何建议。 这就是我到目前为止所做的:R中的削波光栅,r,maps,raster,clip,R,Maps,Raster,Clip,我正在绘制美国东北部的地图。 地图背景需要是海拔图或年平均温度图。我有两个来自Worldclim.org的光栅,它们给了我这些变量,但我需要将它们裁剪到我感兴趣的状态。关于如何做到这一点的任何建议。 这就是我到目前为止所做的: #load libraries library (sp) library (rgdal) library (raster) library (maps) library (mapproj) #load data state<- data (stateMap
#load libraries
library (sp)
library (rgdal)
library (raster)
library (maps)
library (mapproj)
#load data
state<- data (stateMapEnv)
elevation<-raster("alt.bil")
meantemp<-raster ("bio_1.asc")
#build the raw map
nestates<- c("maine", "vermont", "massachusetts", "new hampshire" ,"connecticut",
"rhode island","new york","pennsylvania", "new jersey",
"maryland", "delaware", "virginia", "west virginia")
map(database="state", regions = nestates, interior=T, lwd=2)
map.axes()
#add site localities
sites<-read.csv("sites.csv", header=T)
lat<-sites$Latitude
lon<-sites$Longitude
map(database="state", regions = nestates, interior=T, lwd=2)
points (x=lon, y=lat, pch=17, cex=1.5, col="black")
map.axes()
library(maps) #Add axes to main map
map.scale(x=-73,y=38, relwidth=0.15, metric=T, ratio=F)
#create an inset map
# Next, we create a new graphics space in the lower-right hand corner. The numbers are proportional distances within the graphics window (xmin,xmax,ymin,ymax) on a scale of 0 to 1.
# "plt" is the key parameter to adjust
par(plt = c(0.1, 0.53, 0.57, 0.90), new = TRUE)
# I think this is the key command from http://www.stat.auckland.ac.nz/~paul/RGraphics/examples-map.R
plot.window(xlim=c(-127, -66),ylim=c(23,53))
# fill the box with white
polygon(c(0,360,360,0),c(0,0,90,90),col="white")
# draw the map
map(database="state", interior=T, add=TRUE, fill=FALSE)
map(database="state", regions=nestates, interior=TRUE, add=TRUE, fill=TRUE, col="grey")
#加载库
图书馆(sp)
图书馆(rgdal)
图书馆(光栅)
图书馆(地图)
图书馆(mapproj)
#加载数据
state光栅包中的函数crop允许您使用范围对象或可以计算范围的对象来切割(子集)另一个对象。软件包示例:
r <- raster(nrow=45, ncol=90)
r[] <- 1:ncell(r)
e <- extent(-160, 10, 30, 60)
rc <- crop(r, e)
r查看上的或MapGAM
软件包。一个可以帮助你自己建造。另一个已经有了这样的功能,你可以进行逆向工程。难道没有办法用“地图”来填充光栅吗?或者甚至在光栅中使用覆盖功能?@mnel:您对解决此问题是否有任何意见?交叉发布: