gdal c++ 提取ndvi

GDAL是一个功能强大的地理数据抽象库，用于读取、写入和处理各种栅格数据格式。在使用GDAL提取NDVI（归一化植被指数）时，我们可以按照以下步骤进行操作。

首先，我们需要引入GDAL库，并声明需要使用的数据集和波段。

#include "gdal_priv.h"
#include "cpl_conv.h"

int main()
{
    //声明需要使用的数据集和波段
    GDALDatasetH hDataset;
    GDALRasterBandH hRedBand, hNirBand;
    //...
    return 0;
}

接下来，我们需要打开影像文件，并获取红光和近红外波段。

hDataset = GDALOpen("input_image.tif", GA_ReadOnly);
hRedBand = GDALGetRasterBand(hDataset, 1);
hNirBand = GDALGetRasterBand(hDataset, 2);

然后，我们需要读取红光和近红外波段的数据，并计算NDVI。

int width = GDALGetRasterBandXSize(hRedBand);
int height = GDALGetRasterBandYSize(hRedBand);
float* redData = (float*) CPLMalloc(sizeof(float)*width*height);
float* nirData = (float*) CPLMalloc(sizeof(float)*width*height);
GDALRasterIO(hRedBand, GF_Read, 0, 0, width, height, redData, width, height, GDT_Float32, 0, 0);
GDALRasterIO(hNirBand, GF_Read, 0, 0, width, height, nirData, width, height, GDT_Float32, 0, 0);

float* ndviData = (float*) CPLMalloc(sizeof(float)*width*height);
for (int i = 0; i < width*height; i++) {
    if (redData[i] + nirData[i] == 0) {
        ndviData[i] = 0;
    } else{
        ndviData[i] = (nirData[i] - redData[i]) / (nirData[i] + redData[i]);
    }
}

最后，我们可以将计算得到的NDVI数据保存到新的影像文件中。

GDALDriverH hDriver = GDALGetDriverByName("GTiff");
GDALDatasetH hOutputDataset = GDALCreateCopy(hDriver, "output_ndvi.tif", hDataset, 0, NULL, NULL, NULL);
GDALRasterBandH hOutputBand = GDALGetRasterBand(hOutputDataset, 1);
GDALRasterIO(hOutputBand, GF_Write, 0, 0, width, height, ndviData, width, height, GDT_Float32, 0, 0);

上述代码片段展示了如何使用GDAL C API提取NDVI。通过逐像素计算红光和近红外波段之间的比值，我们可以得到NDVI并将其保存为新的栅格影像文件output_ndvi.tif。