DEM数据转成ascii码之后，求栅格体积（土方量）的C代码或C++代码

以下是计算DEM数据转换为ASCII码之后的栅格体积（土方量）的C++代码示例：

#include <iostream>
#include <fstream>
#include <sstream>
#include <string>

using namespace std;

int main() {
  // 读取文件
  ifstream inFile("DEM_data.asc");
  if (!inFile) {
    cerr << "Error: cannot open file!" << endl;
    return 1;
  }

  // 读取DEM数据的头文件信息
  string line;
  getline(inFile, line); // 忽略第一行
  getline(inFile, line);
  int numCols, numRows;
  double xllCorner, yllCorner, cellSize, noDataValue;
  stringstream ss(line);
  ss >> numCols >> numRows >> xllCorner >> yllCorner >> cellSize >> noDataValue;

  // 计算栅格体积
  double volume = 0.0;
  double elevation;
  for (int i = 0; i < numRows; ++i) {
    for (int j = 0; j < numCols; ++j) {
      inFile >> elevation;
      if (elevation != noDataValue) { // 如果不是缺失值
        volume += elevation * cellSize * cellSize; // 计算栅格体积
      }
    }
  }
  inFile.close();

  // 输出结果
  cout << "The volume of the DEM data in ASCII format is " << volume << endl;

  return 0;
}

其中，"DEM_data.asc"是DEM数据转换为ASCII码后的文件名。代码中，首先读取了DEM数据的头文件信息，包括栅格列数、行数、左下角坐标、栅格单元大小和缺失值等信息。然后，对于每个栅格单元，读取其高程值，如果高程值不是缺失值，则计算栅格体积，最后输出结果。