方法一:R语言 (可直接提取并导出数据,快慢依电脑性能和网速而定...)
soil_world在R中的使用指南(可以直接调用SG的数据)
代码如下 ///感谢zihao师兄鼎力相助
install.packages("geodata")
library(geodata)
library(raster)
library(sp)
library(rgdal)
# Get soil grid layers (tif files)
gph <- soil_world(var="phh2o", depth=5, path=tempdir())
#Image visualization
plot(gph, ylim = c(-60,90), xlim = c(-180,180),
col = colorRampPalette(c("#3288BD","#66C2A5","#FEE08B","#D53E4F"))(250))
# create RasterLayer
ph <- raster(gph)
plot(ph)
# Get soil grid points
data = read.csv("point.csv")
coordinates(data) = c("longitude", "latitude")
phd <- extract(x = ph, y = data)
datapoint <-data.frame (data,phd)
write.csv(datapoint,"E:/Seoyeol/Paper/4ENG2/Data/ph.csv")#生成csv文件
ps:一些栅格数据的R语言处理方法,下面这个博客有介绍