8 环境
8.1 经纬度格式转化
支持向量,将度分秒格式转为°格式
a<-"12°35′65″"
DD(a)
## [1] 12.60138889
8.2 WC气候因子提取
参考代码,wt为包含long,lat(°格式的经纬度的data.fram) 将bio1.bil格式的WC网站下载文件,加载library(raster)包,利用raster函数读取
TP<-function(wt){
mat<-raster("bio1.bil")
map<-raster("bio12.bil")
wt$long<-as.numeric(wt$long)#将经度°转为数字格式
wt$lat<-as.numeric(wt$lat)
wt$long<-signif(wt$long,10)#保留小数位
wt$lat<-signif(wt$lat,10)
cells<-cellFromXY(mat,wt[,c(3,4)])#提取单元格位置
map1<-map[cells] #降雨信息
mat1<-mat[cells]/10 #温度信息,除以10
wt$mat<-mat1
wt$map<-map1
return(wt)
}
8.3 spocc包
文档链接 支持GBIF数据获取,Leaflet 绘图,WC气候数据获取