C语言读取TIFF格式DEM数据并计算填挖方体积

本文提供了一个简单的C语言程序，可以读取TIFF格式的DEM数据，获取每个像素的高程值，并计算填挖方体积。

程序代码：

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

#define DEM_FILE 'dem.tif'
#define PIXEL_SIZE 5.0

int main()
{
    TIFF* tif = TIFFOpen(DEM_FILE, 'r');
    if (!tif) {
        printf('Failed to open DEM file %s
', DEM_FILE);
        return 1;
    }

    uint32_t width, height;
    TIFFGetField(tif, TIFFTAG_IMAGEWIDTH, &width);
    TIFFGetField(tif, TIFFTAG_IMAGELENGTH, &height);

    float* data = malloc(sizeof(float) * width * height);
    TIFFReadFloat(tif, data);

    float fill_volume = 0.0;
    float cut_volume = 0.0;

    for (int y = 0; y < height; y++) {
        for (int x = 0; x < width; x++) {
            float elevation = data[y * width + x];
            if (elevation < 0.0) {
                cut_volume += (0.0 - elevation) * PIXEL_SIZE * PIXEL_SIZE;
            }
            else {
                fill_volume += elevation * PIXEL_SIZE * PIXEL_SIZE;
            }
        }
    }

    printf('Fill volume: %f
', fill_volume);
    printf('Cut volume: %f
', cut_volume);

    free(data);
    TIFFClose(tif);

    return 0;
}

解释：

包含头文件： 程序首先包含必要的头文件，stdio.h 用于标准输入输出，stdlib.h 用于内存分配，tiffio.h 用于操作TIFF文件（需要安装LibTIFF库）。
定义常量： DEM_FILE 定义DEM文件路径，PIXEL_SIZE 定义像素大小（单位：米）。
打开DEM文件： 使用 TIFFOpen 函数打开DEM文件，并检查是否成功打开。
获取图像信息： 使用 TIFFGetField 函数获取DEM图像的宽度和高度。
读取DEM数据： 使用 TIFFReadFloat 函数将DEM数据读取到一个浮点型数组 data 中。
计算填挖方体积： 遍历每个像素，根据高程值计算填挖方体积。如果高程值小于0，则计算挖方体积；如果高程值大于0，则计算填方体积。
输出结果： 打印计算得到的填方和挖方体积。
释放内存和关闭文件： 使用 free 函数释放 data 数组占用的内存，使用 TIFFClose 函数关闭DEM文件。

注意事项：

该程序假设DEM文件中的高程值单位为米，因此计算出的体积也是以立方米为单位。如果DEM文件中的高程值单位不是米，需要将计算出的体积乘以相应的系数进行转换。
程序假设DEM数据存储为浮点型。如果数据类型不同，需要修改代码中相应的数据类型和函数调用。