30米空间分辨率中国年度耕地栅格数据集(1986-2021)


在这里插入图片描述

数据下载地址

点击这里下载数据

数据指标说明

一、数据介绍
精确、详细且及时的耕地范围信息对于粮食安全保障和环境可持续性至关重要。然而,由于农业景观的复杂性和足够训练样本的缺乏,在大范围下进行高时空分辨率的耕地动态监测仍然具有挑战性,尤其是对于农业用地发生剧烈变化的地区。在这里,我们开发了一种新颖而高效的年度耕地制图框架,该框架集成了时间序列Landsat影像、自动化训练样本生成以及机器学习和变化检测技术。我们将所提出的方案实施到谷歌地球引擎的云计算平台上,并首次以30米空间分辨率生成了中国年度耕地数据集(CACD)。结果表明,我们的方法能够跟踪不同农业区域的动态耕地变化。CACD年度地图和变化地图的像素尺度F1分数分别为0.79±0.02和0.81。更进一步,跨产品间的精度评估、与统计数据的相关性以及空间细节方面显示出CACD相比于其他数据集的精确性和稳健性。根据我们的估计,从1986年到2021年,中国的耕地总面积扩大了30, 300平方公里(1.79%),在2000年之前经历了增长,但在2000年至2015年期间总体下降,之后略有恢复。耕地扩张集中在西北部,而东部沿海地区经历了大量的耕地流失。此外,我们观察到在研究期间至少有419, 342平方公里(17.57%)的耕地至少发生了一次弃耕。一致而高分辨率的CACD数据集可以被用于各种研究应用中,以促进可持续农业利用和粮食生产进展。
数据容量:21.25GB
数据时间:1986-2021
数据格式:tif
二、数据作者引用方式
引用方式
Tu, Y., Wu, S., Chen, B., Weng, Q., Gong, P., Bai, Y., Yang, J., Yu, L., and Xu, B.: A 30 m annual cropland dataset of China from 1986 to 2021, Earth System Science Data Discussion. [preprint], https://doi.org/10.5194/essd-2023-190, in review, 2023.

项目备注

本公众号资源部分自主研发,部分来源于网络。
来源于网络的资源版权归出版社或原作者所有,仅供学习,请勿用于商业。

数据下载地址

点击这里下载数据

  • 2
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
使用 GDALDataset::RasterIO 函数来降低栅格数据集分辨率,可以通过设置读取和写入的像素数来实现。假设原始数据集分辨率为 $w_0 \times h_0$,要将其降低到 $w_1 \times h_1$,可以先计算出每个像素的宽度和高度比例: ```c++ double scaleX = (double)w_0 / w_1; double scaleY = (double)h_0 / h_1; ``` 然后,可以创建一个大小为 $w_1 \times h_1$ 的缓冲区,将原始数据集中的数据按比例写入到缓冲区中: ```c++ // 打开原始数据集 GDALDataset *pSrcDS = (GDALDataset *) GDALOpen("path/to/src/dataset", GA_ReadOnly); // 创建输出数据集 GDALDriver *pDriver = GetGDALDriverManager()->GetDriverByName("GTiff"); GDALDataset *pDstDS = pDriver->Create("path/to/dst/dataset", w_1, h_1, pSrcDS->GetRasterCount(), GDT_Float32, NULL); // 计算每个像素的宽度和高度比例 double scaleX = (double)w_0 / w_1; double scaleY = (double)h_0 / h_1; // 设置读取和写入的像素数 int nXSize = (int) ceil(scaleX); // 保证至少读取一个像素 int nYSize = (int) ceil(scaleY); // 创建缓冲区 float *pBuf = (float *) CPLMalloc(nXSize * nYSize * sizeof(float)); // 逐行读取原始数据集的数据,并按比例写入到缓冲区中 for (int y = 0; y < h_1; y++) { for (int x = 0; x < w_1; x++) { // 计算在原始数据集中的位置 int srcX = (int) (x * scaleX); int srcY = (int) (y * scaleY); // 读取数据 CPLErr err = pSrcDS->RasterIO(GF_Read, srcX, srcY, nXSize, nYSize, pBuf, nXSize, nYSize, GDT_Float32, pSrcDS->GetRasterCount(), NULL, 0, 0, 0); if (err != CE_None) { // 处理错误 // ... } // 计算平均值 float sum = 0; for (int i = 0; i < nXSize * nYSize; i++) { sum += pBuf[i]; } float avg = sum / (nXSize * nYSize); // 写入数据 err = pDstDS->RasterIO(GF_Write, x, y, 1, 1, &avg, 1, 1, GDT_Float32, pSrcDS->GetRasterCount(), NULL, 0, 0, 0); if (err != CE_None) { // 处理错误 // ... } } } // 释放缓冲区 CPLFree(pBuf); // 关闭数据集 GDALClose(pSrcDS); GDALClose(pDstDS); ``` 在上述代码中,我们首先打开原始数据集,然后创建大小为 $w_1 \times h_1$ 的输出数据集。接下来,我们计算出每个像素的宽度和高度比例,并根据比例设置读取和写入的像素数。然后,我们创建一个大小为 $nXSize \times nYSize$ 的缓冲区,逐行读取原始数据集的数据,并按比例写入到缓冲区中。对于每个输出像素,我们计算出其对应的原始像素的平均值,并将其写入到输出数据集中。最后,我们释放缓冲区,并关闭数据集。 需要注意的是,上述代码仅适用于单波段数据集,如果处理多波段数据集,则需要在读取和写入数据时指定要处理的波段。同时,如果原始数据集分辨率不能整除目标分辨率,则需要对每个输出像素进行加权平均,权重为其覆盖的原始像素的面积。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Android安卓科研室

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值