如何使用 R 一次性将多个位置的 GeoTiff 数据提取/转换为 csv？

Question

在 WorldClim 网站中，有 2021-2100 年期间 SSP [126、245、370 和 585] 的未来气候数据，变量的不同分辨率，例如来自不同模型的 tn、tx、pr 和 bc 。同样，也有历史数据。

我计划通过一次针对一个模型中的一个变量来管理存储空间。当我从模型 MIROC6 以 30 秒分辨率下载数据（例如，变量“pr”）后，是否可以使用 R 将多个特定位置的数据从 GeoTiff 文件提取到 csv 中？

所需的经纬度位置为 -[28.03,88.71]、[27.37,88.76]、[27.71,88.56]、[27.33,88.61] 和 [27.90,88.53]。

来自 WorldClim 网站的参考链接：

1. 2021-2100 年 30 秒分辨率的未来数据：https://www.worldclim.org/data/cmip6/cmip6_clim30s.html
2. 1970-2000 年 30 秒历史气候数据：https://www.worldclim.org/data/worldclim21.html
3. 1960-2018 年历史月度天气数据，分辨率为 2.5 分钟：https://www.worldclim.org/data/monthlywth.html#

谢谢你， 问候， 个人电脑 研究学者， CRDIST，印度理工学院克勒格普尔。 联系方式：[电子邮件受保护]

Answer 1

正如对您的问题的评论中提到的，这不是适合您的问题类型的论坛，但我将在这里给您一些提示：

• 据我所知，您无法从 WorldClim 下载小区域的数据。您需要下载完整的数据集，正如您已经观察到的那样，这是大量数据。
• 您正在查找特定位置的数据，因此查找 30 弧秒数据是没有意义的。相反，您应该获取纬度 26 - 30 度和经度 86 - 92 度的边界框的数据。这超出了您的要求，但请继续阅读。
• 您可以从气候数据存储下载所有 CMIP5 或 CMIP6 模型、实验和变量的小区域数据。这需要一些努力，因为您必须单独执行每个组合，但您可以同时暂存多个请求，然后将它们下载为较小的文件。请注意，您还可以获得这样的每日数据。许多生物气候变量可以根据每日 tmin、tmax 和 pr 计算出来，或者您也可以从 CDS 下载它们，作为气候学并用于多个模型和实验 RCP4.5 和 RCP8.5（例如基于 CMIP5 数据）。
• 这些模型的原始分辨率为 0.5 到 2 度，因此您应该指定更大的边界框。获得数据后，您必须对几个像素进行插值以找到每个位置感兴趣的像素。请注意，不同的模型具有不同的网格，因此您必须对每个模型重复插值。
• 既然你要开始攻读博士学位，就要做好充分的准备。不要使用 CSV 文件，然后使用 Excel 或类似的东西，而是投资学习 R。