voronoi泰森多边形设置边界且求面积

function vorarea
%UNtitleD Summary of this function goes here
%多边形面积
%   Detailed explanation goes here
clear;clc;
m=5;

x=[4228325869];
y=[8423433769];
%  [vc]=voronoin（[x;y]‘）   % 不设置边界时的情况，v是所有的voronoi边的顶点坐标，c{i}就是第i个点周围的voronoi边的顶点下标
pbound=Polyhedron（[0 0;0 100;100 100;100 0]）;
VV=mpt_voronoi（[x; y]‘bound‘pbound）;
for iPoly = 1:m
    % make points in order
    k = convhull（VV.Set（iPoly）.V（:1）VV.Set（iPoly）.V（:2））;
    % make clockwise
    [pXpY] = poly2cw（VV.Set（iPoly）.V（k1）VV.Set（iPoly）.V（k2））;