C++ TIFF 图像处理：计算填挖方体积

本代码使用 C++ 和 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;
}

运行步骤：

打开命令行窗口（Windows）或终端（Mac/Linux）。
进入代码所在的目录，例如：
```
cd D:\kjfx\Project9\Project9
```
编译代码，命令为：
```
gcc dem.cpp -o dem -ltiff
```
解释一下这个命令：
- "gcc" 是编译器的命令，用于将源代码编译成可执行文件。
- "dem.cpp" 是源代码文件的名称。
- "-o dem" 表示将编译生成的可执行文件命名为 "dem"。
- "-ltiff" 表示链接 TIFF 库，这里假设已经安装了 TIFF 库。
运行程序，命令为：
```
./dem D:\kjfx\dem.tif
```
解释一下这个命令：
- "./dem" 表示运行当前目录下的 "dem" 可执行文件。
- "D:\kjfx\dem.tif" 是 DEM 图像文件的路径和名称。
程序会输出总填方体积和总挖方体积，单位为立方米。

代码说明：

该代码使用 TIFF 库来打开和读取 TIFF 文件。
ROWS、COLS 和 PIXEL_SIZE 是 DEM 图像的尺寸和像素大小。
代码逐行读取 DEM 图像数据，并根据每个像素的高度值计算填方体积或挖方体积。
最终输出总填方体积和总挖方体积。

注意：

确保已经安装了 TIFF 库。
代码中的 dem.cpp 文件和 dem.tif 文件应替换成你的实际文件路径。
代码中假设 DEM 图像的像素大小为 5 米，如果你的 DEM 图像的像素大小不同，请修改代码中的 PIXEL_SIZE 常量。

更多信息：

TIFF 库文档：https://www.libtiff.org/
DEM 图像格式：https://en.wikipedia.org/wiki/Digital_elevation_model