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气候数据获取

8.4 根据经纬度,解析行政地址

library(baidumap),申请百度地图API 参考1
参考2