c实现读取并计算dem数据的填方挖方

由于DEM数据是以栅格形式存储的，因此可以通过读取DEM数据的方式来计算填方挖方。

首先，需要读取DEM数据的头文件信息，包括栅格的行列数、坐标系、分辨率等信息。然后，读取DEM数据中每个栅格的高程值，并根据高程值的变化来计算填方挖方。

具体实现步骤如下：

读取DEM数据头文件信息，包括栅格的行列数、坐标系、分辨率等信息。
根据DEM数据的行列数，创建一个二维数组存储DEM数据中每个栅格的高程值。
读取DEM数据中每个栅格的高程值，并将其存储到二维数组中。
遍历二维数组，计算填方挖方。
对于每个栅格，计算其四周栅格的高程值差（即该栅格高程值与周围栅格高程值的差值），并根据高程值差的正负来判断是填方还是挖方。
将填方和挖方的体积分别累加到对应的变量中。
输出填方和挖方的体积。

下面是一个简单的C语言代码实现：

#include <stdio.h>

#define ROWS 10
#define COLS 10

int main()
{
    // 读取DEM数据头文件信息
    int rows = ROWS; // DEM数据的行数
    int cols = COLS; // DEM数据的列数
    double cellSize = 1.0; // DEM数据的分辨率
    double xllCorner = 0.0; // DEM数据左下角的X坐标
    double yllCorner = 0.0; // DEM数据左下角的Y坐标

    // 创建二维数组存储DEM数据
    double dem[ROWS][COLS];

    // 读取DEM数据
    for (int i = 0; i < rows; i++) {
        for (int j = 0; j < cols; j++) {
            scanf("%lf", &dem[i][j]);
        }
    }

    // 计算填方挖方
    double fillVolume = 0.0; // 填方体积
    double cutVolume = 0.0; // 挖方体积

    for (int i = 1; i < rows - 1; i++) {
        for (int j = 1; j < cols - 1; j++) {
            double diff1 = dem[i][j] - dem[i-1][j]; // 上方栅格高程差
            double diff2 = dem[i][j] - dem[i][j-1]; // 左方栅格高程差
            double diff3 = dem[i][j] - dem[i+1][j]; // 下方栅格高程差
            double diff4 = dem[i][j] - dem[i][j+1]; // 右方栅格高程差

            if (diff1 >= 0 && diff2 >= 0 && diff3 >= 0 && diff4 >= 0) {
                fillVolume += cellSize * cellSize * dem[i][j];
            } else if (diff1 <= 0 && diff2 <= 0 && diff3 <= 0 && diff4 <= 0) {
                cutVolume += cellSize * cellSize * dem[i][j];
            }
        }
    }

    // 输出填方挖方体积
    printf("Fill volume: %lf\n", fillVolume);
    printf("Cut volume: %lf\n", cutVolume);

    return 0;
}

注意，上述代码仅适用于二维平面上的DEM数据，对于三维DEM数据需要进行适当修改。此外，代码中未考虑DEM数据的边界情况，实际应用中需要对边界情况进行特殊处理