matlab实现计亩，多边形面积计算，相交面积计算

function [multi_outer_points, multi_inner_points] = multi_polygonFit(multi_points, w) num_lines = length(multi_points); multi_outer_points = cell(1, num_lines); multi_inner_points = cell(1, num_lines); for j = 1:num_lines points = multi_points{j}; n = size(points, 1); outer_points = zeros(n, 2); inner_points = zeros(n, 2); for i = 1:n curr_point = points(i, :); if i == 1 prev_point = points(n, :); else prev_point = points(i-1, :); end if i == n next_point = points(1, :); else next_point = points(i+1, :); end % 对于第一个和最后一个点，只考虑一个方向向量 if i == 1 dir_vec = next_point - curr_point; elseif i == n dir_vec = curr_point - prev_point; else % 对于其他点，计算前一个方向向量和后一个方向向量 dir_vec_prev = curr_point - prev_point; dir_vec_next = next_point - curr_point; % 使用平均法向量 dir_vec = (dir_vec_prev + dir_vec_next) / 2; end % 计算法向量 normal_vec = [-dir_vec(2), dir_vec(1)]; % 单位化法向量 normal_vec = normal_vec / norm(normal_vec); % 计算外部和内部点 outer_points(i, :) = curr_point + w/2 * normal_vec; inner_points(i, :) = curr_point - w/2 * normal_vec; end multi_outer_points{j} = outer_points; multi_inner_points{j} = inner_points; end end