WorldClim 应用

worldClim官网

https://www.worldclim.org/data/worldclim21.html

此网站包含了世界各地的气候数据：

我的理解是，它把整个地球表面割成了许多相等大小的小方块，分辨率有4种：

30秒这个，代表1平方千米大小面积的方块，其它依次增大。分的小块越小，总的数据也就越大，我下载的是2.5分钟大小的数据，大约1.2G。

这些数据如何下载应用呢？我的目的是对应到我的取样样本地点的气候数据，我有取样点的经纬度，那如何得到相应气候数据呢？

我打算用已有的R包来实现，主要用到的R包是：

library(rgdal)

library(raster)

安装这两个R包，也是废了功夫啊，我是安装在我的mac笔记本上的，集群上我更安不上。

需要安装gdal这个库，Mac用brew install gdal安装。这个应该是需要sudo权限，所以集群上我安不上。（brew这儿稍后写一下）

工具安装好了，下面就是如何用了（下面都是R命令）。

w = getData('worldclim', var='bio', res=2.5) 

此命令下载分辨率为2.5分钟的数据，可以指定下载位置，我就下载到当前文件夹了，getData是raster这个R包的命令，具体raster用法见https://cran.r-project.org （CRAN官网找这个包）

ll<-read.table("**/**_cluster_location_k4_v2.txt", header=T, sep="\t")

此命令读取经纬度信息,"**/**_cluster_location_k4_v2.txt"文件格式如下，只给前三列就可以。

extract(w, ll[c(2,3)])

此命令用于提取气候数据。

此数据可以结合样本SNP数据，进行RDA分析，找到环境相关的突变位点。RDA分析我后面额外写。

其它可以参考的相关文章：

https://www.gis-blog.com/r-raster-data-acquisition/

https://www.jianshu.com/p/c34fb73f5a57

https://www.modb.pro/db/541343

欢迎留言讨论！