我正在尝试计算多个物种栅格文件(Maxent 输出)的栅格宽度,并使用 for 循环执行此操作。计算宽度的每个循环我也打算将结果保存为 csv 文件,但我遇到了问题。
这是我到目前为止的代码:
data <- read.csv("orchid_occurrences.csv", sep = ";", stringsAsFactors = FALSE)
spp_to_model <- unique(data$species)
orchid <- paste0("./orchid_data/",spp_to_model,"_avg.asc")
for(spp_to_model in orchid){
spp_raster <- raster(spp_to_model)
spp_breadth <- raster.breadth(spp_raster)
write.table(spp_breadth, spp_to_model,"_breadth.csv", row.names = FALSE)
}
我稍微调整了一下以适合我在此处共享的发生数据和栅格文件的示例:https ://drive.google.com/file/d/1oA78dnZBSd9B_9EJuaWxg2isK-9lfQ5Q/view?usp=sharing
提前致谢!