学习网址：https://blog.csdn.net/SUN5_The_answer/article/details/135131492
library(raster)
library(terra)
library(ncdf4)
occur <- read.csv("sample information1.csv", header=T, row.names = "Code")
objects(occur)
bio1 <- raster("wc2.1_10m_bio_1.tif")
crs.wgs <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"
crs(bio1) <- crs.wgs
occ <- SpatialPointsDataFrame(
occur[,c('Lon','Lat')],
proj4string=CRS(crs.wgs),
data=occur
)
lo_bio1 <- raster::extract(bio1,occ)
for (i in 1:19) {
var = raster(paste0("wc2.1_10m_bio_",i,".tif"))
var_val = raster::extract(var, occ)
occur = cbind(occur, var_val)
colnames(occur)[i+2] <- paste0("bio",i)
}
write.csv(occur, 'bio_val.csv', quote = FALSE)
Rt提取气候因子
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 工具准备 ArcGIS Desktop下载安装 软件的下载和安装自行解决，可参考文末提供的参考资料。 Worldc...
- 写在前面 昨日，我推出了两个 TBtools 小插件，主要用于快速基于输入的蛋白序列和启动子序列，预测一个初级的基...
