C语言实现DEM数据填方挖方计算

由于DEM数据是以栅格形式存储的，因此可以通过遍历每个栅格点，计算该点的高程值与基准高程值的差值，进而得到该点的填方或挖方量。

具体实现步骤如下：

1. 读取DEM数据，并获取基准高程值。

2. 遍历每个栅格点，计算该点的高程值与基准高程值的差值。

3. 如果差值大于0，则表示该点需要填方，填方量为差值乘以栅格面积。

4. 如果差值小于0，则表示该点需要挖方，挖方量为差值乘以栅格面积。

5. 统计所有填方量和挖方量，得到总的填方量和挖方量。

6. 输出填方量和挖方量。

代码实现：

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

#define ROWS 10
#define COLS 10
#define CELL_SIZE 10.0 //栅格面积

int main()
{
    float dem[ROWS][COLS]; //DEM数据
    float base_height = 0.0; //基准高程值
    float fill_volume = 0.0; //填方量
    float cut_volume = 0.0; //挖方量

    //读取DEM数据
    //此处省略读取过程，直接赋值
    for(int i=0; i<ROWS; i++)
    {
        for(int j=0; j<COLS; j++)
        {
            dem[i][j] = i+j;
        }
    }

    //获取基准高程值
    base_height = dem[0][0];

    //遍历每个栅格点，计算填方量和挖方量
    for(int i=0; i<ROWS; i++)
    {
        for(int j=0; j<COLS; j++)
        {
            float diff = dem[i][j] - base_height;
            if(diff > 0)
            {
                fill_volume += diff * CELL_SIZE;
            }
            else if(diff < 0)
            {
                cut_volume += (-diff) * CELL_SIZE;
            }
        }
    }

    //输出填方量和挖方量
    printf('Fill volume: %.2f\n', fill_volume);
    printf('Cut volume: %.2f\n', cut_volume);

    return 0;
}