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

由于缺乏具体的DEM数据和计算方法，以下是一般的读取DEM数据和计算填方挖方的C语言代码示例：

读取DEM数据：

FILE *fp;
fp = fopen("dem.txt", "r");
if (fp == NULL) {
    printf("Failed to open file.\n");
    return 1;
}
int rows, cols;
fscanf(fp, "%d %d", &rows, &cols);
float **dem;
dem = (float **)malloc(rows * sizeof(float *));
for (int i = 0; i < rows; i++) {
    dem[i] = (float *)malloc(cols * sizeof(float));
    for (int j = 0; j < cols; j++) {
        fscanf(fp, "%f", &dem[i][j]);
    }
}
fclose(fp);

计算填方挖方：

float fill_volume = 0.0, cut_volume = 0.0;
float unit_width = 1.0, unit_length = 1.0, unit_height = 1.0;
for (int i = 0; i < rows - 1; i++) {
    for (int j = 0; j < cols - 1; j++) {
        float a = dem[i][j], b = dem[i][j+1], c = dem[i+1][j+1], d = dem[i+1][j];
        float avg_height = (a + b + c + d) / 4.0;
        float diff = avg_height - unit_height;
        if (diff > 0) {
            fill_volume += diff * unit_width * unit_length;
        } else if (diff < 0) {
            cut_volume += (-diff) * unit_width * unit_length;
        }
    }
}
printf("Fill volume: %.2f cubic meters.\n", fill_volume);
printf("Cut volume: %.2f cubic meters.\n", cut_volume);

这里的DEM数据是以文本文件形式存储的，第一行包含行数和列数，后面是每个格点的高程值。计算填方挖方时，将DEM数据按照每个网格计算平均高程，然后与给定的单位高度比较，如果高于单位高度则计算填方体积，否则计算挖方体积。最终输出填方体积和挖方体积。注意，这里的单位宽度、单位长度和单位高度是可以根据具体情况设定的