从栅格堆栈(R)中提取具有相关值的经纬度

时间:2018-10-19 09:54:20

标签: r

我经常想探讨现代气候变量的模式,因为它们与我的古数据中出现的问题有关。在这种情况下,我对Tmin的分布感兴趣,因为它与Tmax的子集有关。我使用的气候数据是来自Worldclim的19个Bioclim变量的栅格堆栈。我不熟悉使用Raster Stacks,并且想将特定数据提取到一个数据框中,以便我可以更直观地使用它。在其他场合,我也希望为其他目的做类似的事情,这要求将最终数据集导出为.csv,因此我有多个原因想要导出数据。

我可以使用以下代码读取数据并提取我感兴趣的气候值。这给了我两列数据帧,分别为Tmax和Tmin。我还需要另外两列,即与每对数据关联的纬度和经度。到目前为止,我仅发现了一些示例,人们希望按照给定的坐标进行提取,并且他们用于提取数据的通话中的经纬度可以绑定到数据帧,但是我正在提取给定范围内的所有数据。

到目前为止,我的代码:

wc_bio <-getData("worldclim",var="bio",res=2.5)
minmax <- wc_bio[[c(5,6)]]
names(minmax) <- c("tmax","tmin")

#设置要从(北半球)提取的地图的范围

`extN <- extent(-180,180,0,90)`

#将tmax和tmin的值提取到数据帧中

df <- extract(minmax,extN, df=TRUE)

此外,这些是栅格堆栈wc_bio的特征

class       : RasterStack 
dimensions  : 3600, 8640, 31104000, 19  (nrow, ncol, ncell, nlayers)
resolution  : 0.04166667, 0.04166667  (x, y)
extent      : -180, 180, -60, 90  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 
names       : bio1, bio2, bio3, bio4, bio5, bio6, bio7, bio8, bio9, bio10, bio11 ... 

和minmax

class       : RasterStack 
dimensions  : 3600, 8640, 31104000, 2  (nrow, ncol, ncell, nlayers)
resolution  : 0.04166667, 0.04166667  (x, y)
extent      : -180, 180, -60, 90  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 
names       : tmax, tmin 

有人知道如何提取或以其他方式生成提取的气候值的纬度和经度矢量吗?

0 个答案:

没有答案