方格网法计算土方量matlab代码

以下是方格网法计算土方量的 MATLAB 代码：

% 输入参数：x,y为横纵坐标，z为高程
% x,y,z必须为向量，且长度相等
% dx,dy为方格大小
% 输出参数：V为土方量
function V = earthwork(x,y,z,dx,dy)
    % 计算x,y方向上的网格数
    nx = ceil((max(x)-min(x))/dx);
    ny = ceil((max(y)-min(y))/dy);
    % 初始化土方量为0
    V = 0;
    % 循环计算每个网格的土方量
    for i = 1:nx
        for j = 1:ny
            % 计算网格的四个顶点
            x1 = min(x) + (i-1)*dx;
            x2 = min(x) + i*dx;
            y1 = min(y) + (j-1)*dy;
            y2 = min(y) + j*dy;
            % 判断顶点是否在数据范围内
            if x1<min(x) || x2>max(x) || y1<min(y) || y2>max(y)
                continue
            end
            % 判断顶点是否在数据点内
            idx = x>=x1 & x<=x2 & y>=y1 & y<=y2;
            if ~any(idx)
                continue
            end
            % 计算网格内的土方量
            z1 = min(z(idx));
            z2 = max(z(idx));
            V = V + (z1+z2)/2*dx*dy;
        end
    end
end

使用方法：

1. 将上述代码保存为一个.m文件，如earthwork.m。

2. 在MATLAB命令窗口中输入以下代码：

x = [1 2 3 4 5 6 7 8 9 10];
y = [1 1 1 1 1 2 2 2 2 2];
z = [0 0 0 0 0 1 1 1 1 1];
dx = 1;
dy = 1;
V = earthwork(x,y,z,dx,dy)

其中，输入参数x、y、z为示例数据点的坐标和高程，dx、dy为方格大小，输出参数V为土方量。实际使用时，需根据实际数据进行修改。