Triangulación de Delaunay restringida de un mapa geográfico

Cree una triangulación de Delaunay restringida utilizando un mapa del perímetro de Estados Unidos.

Cargue un mapa del perímetro de los Estados Unidos contiguos.

load usapolygon

Defina un límite de borde entre dos puntos sucesivos que conforman el límite poligonal y cree la triangulación de Delaunay. Esta triangulación abarca un dominio acotado por la envolvente convexa del conjunto de puntos. El conjunto de datos contiene puntos de datos duplicados; es decir, dos o más puntos de datos tienen la misma ubicación. Los puntos duplicados se rechazan y delaunayTriangulation vuelve a establecer el formato a los límites en consecuencia.

nump = numel(uslon);
C = [(1:(nump-1))' (2:nump)'; nump 1];
dt = delaunayTriangulation(uslon,uslat,C);

Warning: Duplicate data points have been detected and removed.
 The Triangulation indices and constraints are defined with respect to the unique set of points in delaunayTriangulation.

Warning: Intersecting edge constraints have been split, this may have added new points into the triangulation.

Filtre los triángulos que se encuentran dentro del dominio del polígono y represéntelos.

io = isInterior(dt);
patch(Faces=dt(io,:),Vertices=dt.Points,FaceColor="r")
axis equal
axis([-130 -60 20 55])
title("Constrained Delaunay Triangulation of usapolygon")