C++ TIFF 文件处理示例：计算填挖方体积

该代码示例使用 C++ 和 libtiff 库读取 TIFF 文件，计算 DEM 数据的填挖方体积。代码示例包含详细的注释，并提供运行说明。

#include <stdio.h>
#include <stdlib.h>
#include <tiffio.h>

#define ROWS 327
#define COLS 486
#define PIXEL_SIZE 5

int main(int argc, char* argv[]) {
    TIFF *tif;
    float *data;
    uint32 width, height;
    uint16 bits_per_sample, samples_per_pixel;
    tsize_t strip_size;
    unsigned long long total_fill_volume = 0, total_cut_volume = 0;
    int i, j;

    if (argc < 2) {
        fprintf(stderr, "Usage: %s <input file>
", argv[0]);
        exit(1);
    }

    tif = TIFFOpen(argv[1], "r");
    if (!tif) {
        fprintf(stderr, "Error opening TIFF file
");
        exit(1);
    }

    TIFFGetField(tif, TIFFTAG_IMAGEWIDTH, &width);
    TIFFGetField(tif, TIFFTAG_IMAGELENGTH, &height);
    TIFFGetField(tif, TIFFTAG_BITSPERSAMPLE, &bits_per_sample);
    TIFFGetField(tif, TIFFTAG_SAMPLESPERPIXEL, &samples_per_pixel);

    if (width != COLS || height != ROWS || bits_per_sample != 32 || samples_per_pixel != 1) {
        fprintf(stderr, "Invalid TIFF format
");
        exit(1);
    }

    strip_size = TIFFStripSize(tif);
    data = (float*) _TIFFmalloc(strip_size);
    if (!data) {
        fprintf(stderr, "Error allocating memory
");
        exit(1);
    }

    for (i = 0; i < height; i++) {
        TIFFReadEncodedStrip(tif, TIFFComputeStrip(tif, i, 0), data, strip_size);
        for (j = 0; j < width; j++) {
            float elevation = data[j];
            if (elevation > 0) {
                total_fill_volume += (PIXEL_SIZE * PIXEL_SIZE * elevation);
            } else {
                total_cut_volume += (PIXEL_SIZE * PIXEL_SIZE * (-elevation));
            }
        }
    }

    printf("Total fill volume: %llu cubic meters
", total_fill_volume);
    printf("Total cut volume: %llu cubic meters
", total_cut_volume);

    _TIFFfree(data);
    TIFFClose(tif);

    return 0;
}

运行说明

确保已安装 libtiff 库。
将代码保存为 dem.cpp 文件，并将 dem.tif 文件放置在与 dem.cpp 相同的目录下。
在命令行中进入该目录，并使用以下命令进行编译和运行：

gcc dem.cpp -ltiff -o dem
./dem dem.tif

其中，-ltiff 表示链接 libtiff 库，-o dem 表示输出可执行文件名为 dem，dem.tif 是要处理的 dem 图像文件名。

代码说明

代码使用 libtiff 库读取 TIFF 文件，并获取图像的宽度、高度、位深和像素样本数。
代码检查 TIFF 文件的格式是否符合预期，并分配内存存储数据。
代码逐行读取 TIFF 文件中的数据，并计算每个像素的填挖方体积。
最后，代码打印总的填挖方体积。

注意

代码中的 ROWS、COLS 和 PIXEL_SIZE 定义了 DEM 图像的尺寸和像素大小，需要根据实际情况进行修改。
代码假设 TIFF 文件存储的是单通道浮点数数据。
代码只计算了水平方向的填挖方体积，如果需要计算垂直方向的体积，需要进行相应的修改。

总结

该代码示例演示了如何使用 C++ 和 libtiff 库读取 TIFF 文件，并计算 DEM 数据的填挖方体积。代码示例包含详细的注释，并提供运行说明，方便用户学习和使用。