inkmlbb

fname = 'Wet_Skagit_BB.kml';
x = geog(:,1);
y = geog(:,2);

0001 function [pts] = inkmlbb(fname,x,y)
0002 %fname = 'Wet_Skagit_BB.kml';
0003 %x = geog(:,1);
0004 %y = geog(:,2);
0005 
0006 
0007 triangle = '/usr/local/bin/triangle';
0008 
0009 [latb,lonb,dumz] = read_kml(fname);
0010 [ktri,kx,ky] = tri_from_bndry(lonb,latb,triangle);
0011 
0012 box = [min(kx),max(kx),min(ky),max(ky)];
0013 pts1 = inbox(box,x,y);
0014 
0015 mark = zeros(length(x),1);
0016 
0017 dims = size(ktri);
0018 
0019 for ii=1:length(pts1);
0020 i = pts1(ii);
0021 for j=1:dims(1);
0022   xtri = kx(ktri(j,1:3));
0023   ytri = ky(ktri(j,1:3));
0024   if(isintriangle(xtri,ytri,x(i),y(i)) ) ; mark(i) = 1; end;
0025 end;
0026 end;
0027 
0028 pts = find(mark==1);
0029 
0030 %figure
0031 %patch('Vertices',[kx(:),ky(:)],'Faces',ktri,...
0032 %       'Cdata',kx,'edgecolor','k','facecolor','k');
0033 %hold on;
0034 %plot(x(pts1),y(pts1),'r+');
0035 %plot(x(pts),y(pts),'g+');
0036 %axis equal
0037 %colorbar
0038