C语言使用LibTIFF库计算DEM图像填挖方体积

本代码使用C语言和LibTIFF库解析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:\桌面

编译代码，输入以下命令：

gcc dem.c -o dem -ltiff

其中，-ltiff 表示链接 libtiff 库。 4. 运行代码，输入以下命令：

./dem D:\kjfx\dem.tif

其中，D:\kjfx\dem.tif 是 dem 图像的路径。 5. 程序会输出总填方体积和总挖方体积。

**注意：**在编译和运行代码时，需要确保 libtiff 库已经正确安装并配置好了环境变量。

环境配置：

如果你的系统还没有安装 libtiff 库，你需要先安装它。具体安装方法取决于你的操作系统。一般来说，你可以使用包管理器进行安装，例如在 Ubuntu 上使用 sudo apt install libtiff-dev。

安装完成后，你需要将 libtiff 库的路径添加到系统的环境变量中，这样编译器才能找到它。具体操作方法也取决于你的操作系统。

其他说明：

代码中 ROWS、COLS 和 PIXEL_SIZE 变量需要根据你的 DEM 图像的实际尺寸进行调整。
代码只计算了图像中所有点的填挖方体积，并没有考虑地形坡度等因素的影响。
代码中的计算公式假设图像的坐标系是平面直角坐标系，并且单位是米。
如果你的 DEM 图像使用的是其他坐标系或单位，你需要根据实际情况修改代码中的计算公式。

希望本代码能够帮助你计算 DEM 图像的填挖方体积。如果你有任何疑问，请随时留言。