0001 function [new_coast,slen]=join_cst(coast,tol);
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022
0023
0024
0025
0026
0027
0028
0029
0030
0031
0032
0033
0034 coast=fixcoast(coast);
0035
0036 plot(coast(:,1),coast(:,2));set(gca,'aspectratio',[nan 1]);
0037 figure(gcf);
0038
0039 iseg=find(isnan(coast(:,1)));
0040 nseg=length(iseg)-1;
0041 irem=ones(nseg,1);
0042 z=coast(:,1)+sqrt(-1)*coast(:,2);
0043 zstart=z(iseg(1:nseg)+1);
0044 zstop=z(iseg([2:(nseg+1)])-1);
0045 zends=[zstart(:) zstop(:)];
0046
0047 znew=nan+sqrt(-1)*nan;
0048 k=0;
0049
0050
0051 while(~isempty(find(irem))),
0052
0053 ii=find(irem);
0054 id0=ii(1);
0055 id=id0;
0056 z0=z(iseg(id)+1);
0057 zind=[(iseg(id)+1):(iseg(id+1)-1)];
0058 line(real(z(zind)),imag(z(zind)),'erasemode','none','color','white')
0059 zc=z(iseg(id+1)-1);
0060
0061 irem(id)=0;
0062
0063 found_next=1;
0064 tried_reverse=0;
0065 while (found_next),
0066 ii=find(irem);
0067 nrem=length(ii);
0068 dgap=abs((ones(size(zends(ii,:))))*zc-zends(ii,:));
0069 id=find(dgap<tol);
0070 if(~isempty(id)),
0071 idt=find(dgap==min(dgap(:)));
0072 idt=idt(1);
0073 if(idt>nrem),
0074 id=ii(idt-nrem);
0075 zi=[(iseg(id+1)-1):-1:(iseg(id)+1)];
0076 else
0077 id=ii(idt);
0078 zi=[(iseg(id)+1):1:(iseg(id+1)-1)];
0079 end
0080 ni=length(zi);
0081 line(real(z(zi)),imag(z(zi)),'erasemode','none','color','white')
0082 irem(id)=0;
0083 zc=z(zi(ni));
0084 zind=[zind zi];
0085 else
0086 if(~tried_reverse),
0087 zc=z0;
0088 nzind=length(zind);
0089 zind=zind([nzind:-1:1]);
0090 tried_reverse=1;
0091 else
0092 found_next=0;
0093 end
0094 end
0095 end
0096 k=k+1
0097 znew=[znew; nan+sqrt(-1)*nan; z(zind)];
0098 end
0099 znew=[znew; nan+sqrt(-1)*nan];
0100 iseg=find(isnan(znew));
0101 nseg=length(iseg)-1;
0102 slen=(iseg([2:(nseg+1)])-iseg(1:nseg)-2);
0103 [y,isort]=sort(slen);
0104 isort=flipud(isort);
0105 zind=[];
0106
0107 for i=1:nseg;
0108 zind=[zind (iseg(isort(i)):iseg(isort(i)+1))];
0109 end
0110 znew=znew(zind);
0111 new_coast=[real(znew) imag(znew)];
0112 new_coast=fixcoast(new_coast);
0113 iseg=find(isnan(new_coast(:,1)));
0114 nseg=length(iseg)-1;
0115 slen=(iseg([2:(nseg+1)])-iseg(1:nseg)-2);