用c语言编写一个代码，要求：（1）能够读取dem数据；（2）对读取的dem数据进行填方挖方算法

好的，下面是一个简单的用C语言实现读取DEM数据并进行填方挖方算法的示例代码：

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

#define WIDTH 100
#define HEIGHT 100
#define CELL_SIZE 10

typedef struct {
  float elevation;
} Cell;

typedef struct {
  int width;
  int height;
  Cell cells[WIDTH][HEIGHT];
} Map;

// 读取DEM数据
void read_dem(Map* map, const char* file_name) {
  FILE* file = fopen(file_name, "rb");
  if (file == NULL) {
    fprintf(stderr, "Error: Cannot open file '%s'\n", file_name);
    exit(1);
  }

  // 读取DEM的宽度和高度
  int width, height;
  fread(&width, sizeof(int), 1, file);
  fread(&height, sizeof(int), 1, file);
  if (width != map->width || height != map->height) {
    fprintf(stderr, "Error: Invalid DEM dimensions\n");
    exit(1);
  }

  // 读取DEM中每个格子的高度值
  for (int y = 0; y < height; y++) {
    for (int x = 0; x < width; x++) {
      float elevation;
      fread(&elevation, sizeof(float), 1, file);
      map->cells[x][y].elevation = elevation;
    }
  }

  fclose(file);
}

// 计算填方挖方的结果
float calculate_cut_and_fill(Map* map) {
  float cut = 0.0f;
  float fill = 0.0f;

  for (int y = 0; y < map->height - 1; y++) {
    for (int x = 0; x < map->width - 1; x++) {
      Cell* c00 = &map->cells[x][y];
      Cell* c10 = &map->cells[x + 1][y];
      Cell* c01 = &map->cells[x][y + 1];
      Cell* c11 = &map->cells[x + 1][y + 1];

      float a = c00->elevation;
      float b = c10->elevation;
      float c = c01->elevation;
      float d = c11->elevation;

      float avg = (a + b + c + d) / 4.0f;

      if (avg > a) {
        cut += (avg - a) * CELL_SIZE * CELL_SIZE;
      } else {
        fill += (a - avg) * CELL_SIZE * CELL_SIZE;
      }
    }
  }

  return fill - cut;
}

int main() {
  Map map = { WIDTH, HEIGHT };

  // 读取DEM数据
  read_dem(&map, "dem.bin");

  // 计算填方挖方
  float result = calculate_cut_and_fill(&map);

  printf("Result: %.2f cubic meters\n", result);

  return 0;
}

这段代码定义了一个`Map`结构体来表示地图，其中包含了宽度、高度和每个格子的高度值。`read_dem`函数用来从二进制文件中读取DEM数据，`calculate_cut_and_fill`函数用来计算填方挖方的结果。在`main`函数中，我们首先读取DEM数据，然后调用`calculate_cut_and_fill`函数计算结果，并输出到控制台中。