qt界面实现C+++GDal生成植被覆盖度影像代码

以下是一个基于 Qt 和 GDAL 的简单示例程序，可以生成植被覆盖度影像：

#include <QCoreApplication>
#include <QDebug>
#include <gdal_priv.h>

int main(int argc, char *argv[])
{
    QCoreApplication a(argc, argv);

    // 打开输入数据集
    GDALDataset *inDs = (GDALDataset *) GDALOpen("input.tif", GA_ReadOnly);
    if (inDs == NULL) {
        qDebug() << "无法打开输入数据集！";
        return 1;
    }

    // 获取输入数据集的信息
    int width = inDs->GetRasterXSize();
    int height = inDs->GetRasterYSize();
    int bandCount = inDs->GetRasterCount();
    qDebug() << "输入数据集信息：";
    qDebug() << "宽度：" << width << "像素";
    qDebug() << "高度：" << height << "像素";
    qDebug() << "波段数：" << bandCount;

    // 创建输出数据集
    GDALDriver *driver = GetGDALDriverManager()->GetDriverByName("GTiff");
    GDALDataset *outDs = driver->Create("output.tif", width, height, 1, GDT_Float32, NULL);
    if (outDs == NULL) {
        qDebug() << "无法创建输出数据集！";
        return 1;
    }

    // 设置输出数据集的投影和地理信息
    outDs->SetProjection(inDs->GetProjectionRef());
    outDs->SetGeoTransform(inDs->GetGeoTransform());

    // 计算植被覆盖度
    GDALRasterBand *inBand = inDs->GetRasterBand(1);
    GDALRasterBand *outBand = outDs->GetRasterBand(1);
    float *inData = new float[width * height];
    float *outData = new float[width * height];
    inBand->RasterIO(GF_Read, 0, 0, width, height, inData, width, height, GDT_Float32, 0, 0);
    for (int i = 0; i < width * height; i++) {
        if (inData[i] >= 0.5) {
            outData[i] = 1.0;
        } else {
            outData[i] = 0.0;
        }
    }
    outBand->RasterIO(GF_Write, 0, 0, width, height, outData, width, height, GDT_Float32, 0, 0);

    // 释放内存
    delete[] inData;
    delete[] outData;

    // 关闭数据集
    GDALClose(inDs);
    GDALClose(outDs);

    return 0;
}

这个程序假设输入数据集为 1 个波段的 TIFF 影像，每个像素的值表示植被覆盖度的比例，输出数据集也是 1 个波段的 TIFF 影像，每个像素的值为 0 或 1，分别表示无植被和有植被。具体实现方式是将输入数据集中大于等于 0.5 的像素值设为 1，将小于 0.5 的像素值设为 0。