利用C++对dem数据进行填方挖方算法

填方挖方算法是指根据地表高程数据（如DEM数据）计算地形表面被填方或挖方的体积。以下是一个利用 C++ 对 DEM 数据进行填方挖方算法的基本步骤：

1. 读取 DEM 数据文件并解析高程数据。
2. 定义填方挖方的高度阈值，计算出 DEM 数据中每个像素点的填方或挖方高度（即 DEM 数据中该点高度减去阈值）。
3. 对于每个像素点，根据其填方或挖方高度计算其对总体积的贡献。
4. 对所有像素点的贡献值进行累加，得到填方或挖方的总体积。

以下是一个 C++ 示例代码，用于计算 DEM 数据的填方挖方体积：

#include <iostream>
#include <fstream>
#include <vector>

using namespace std;

const double THRESHOLD = 10.0; // 高度阈值

int main()
{
    // 读取 DEM 数据文件
    ifstream fin("dem.txt");
    int rows, cols;
    fin >> rows >> cols;
    vector<vector<double>> data(rows, vector<double>(cols));
    for (int i = 0; i < rows; ++i) {
        for (int j = 0; j < cols; ++j) {
            fin >> data[i][j];
        }
    }
    fin.close();

    // 计算填方挖方体积
    double fill_volume = 0.0, cut_volume = 0.0;
    for (int i = 0; i < rows - 1; ++i) {
        for (int j = 0; j < cols - 1; ++j) {
            double height1 = data[i][j] - THRESHOLD;
            double height2 = data[i][j+1] - THRESHOLD;
            double height3 = data[i+1][j+1] - THRESHOLD;
            double height4 = data[i+1][j] - THRESHOLD;
            double avg_height = (height1 + height2 + height3 + height4) / 4.0;
            double volume = (avg_height >= 0.0 ? 1.0 : -1.0) * avg_height;
            if (avg_height >= 0.0) {
                fill_volume += volume;
            } else {
                cut_volume += volume;
            }
        }
    }

    // 输出结果
    cout << "Fill volume: " << fill_volume << endl;
    cout << "Cut volume: " << cut_volume << endl;

    return 0;
}

在该示例代码中，我们首先从 `dem.txt` 文件中读取 DEM 数据，并将其存储在二维向量 `data` 中。然后，我们遍历 DEM 数据中的每个像素点，计算出其填方或挖方高度，并根据其高度贡献计算出填方或挖方的总体积。最后，我们输出填方和挖方的总体积。