poli.m

来自「麻省理工学院的人工智能工具箱,很珍贵,希望对大家有用!」· M 代码 · 共 27 行

%Cargamos lospologonos de los paises
STR=load('worldlo.mat','POpatch');

%Dentro hay una variable que me da el nombre.
%   Lo usaremos para saber cual es Spain
[names{1:length(STR.POpatch)}]=deal(STR.POpatch.tag);
i=strmatch('Spain',names);
%Ahora que lo sabemos extraemos sus poligonos
XB=STR.POpatch(i).long;
YB=STR.POpatch(i).lat;

%Definimos cual es el rectangulo que queremos
XC=[-11 5 5 -11 -11]';
YC=[34 34 45 45 34]';

%Y con todo esto definimos obtenemos la mascara.
[XA,YA]=polybool('-',XC,YC,XB,YB,'cutvector');

[X,Y]=meshgrid(linspace(-11,5,100),linspace(34,45,100));
Z=peaks(100);
pcolor(X,Y,Z);
shading('interp')
hold on
%y la pintamos para ver si es lo que queremos
h=patch(XA,YA,'w');
set(h,'EdgeColor','none')
plot(XB,YB,'k')