C语言代码：使用TIFF库读取DEM文件并计算填方和挖方体积

该程序使用C语言和TIFF库读取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;
}

如何使用该代码

确保您的系统已安装TIFF库。
将代码保存为一个名为dem_volume.c的文件。
使用以下命令编译代码：
```
gcc -o dem_volume dem_volume.c -ltiff
```
运行程序，并将您的DEM文件作为参数传递给它：
```
./dem_volume <input_file>
```
例如，如果您的DEM文件名为dem.tif，则命令应为：
```
./dem_volume dem.tif
```
程序将输出填方和挖方总体积。

注意

该程序假设您的DEM文件是32位单精度浮点型数据，并且每个像素只有一个样本（单通道）。
代码中的ROWS、COLS和PIXEL_SIZE常量需要根据您的DEM文件进行调整。
如果您的DEM文件不是TIFF格式，则需要将其转换为TIFF格式才能使用该程序。
该程序仅计算DEM文件中的填方和挖方总体积，并不考虑其他因素，如地形坡度等。

参考资料