天天看点

碎片拼接

1.问题描述:

碎片拼接

2.部分程序:

%cd piecejoint
a=tic;
 %%
 global allP
 global margl margr
 global ld
 global pn pw phallP=brief_pretreat;
 [margl,margr]=extractmargin; lmark=brief_markline_accurate;
 % save('a','lmark')
 % load('a')
 %%%%%%%%%%%%%%%%%%%
 pind=zeros(11,19);
 tick=[0 4.5 10 20 26 30 36 40 47.5 50 57.01 60];
 for k=1:11
     pind(k,:)=find(tick(k)<lmark&lmark<tick(k+1));
 end
 rerow=[4 6 8 10 2 3 5 7 9 11 1];
 pind=pind(rerow,:);
 %%
 %%%%%%%%%%%%%%%%%%%%%%%%%%55
 P=zeros(11,19);
 for k=1:11
     P(k,:)=arrangemethod3_1(pind(k,:));
 end
 %writecomP(P,1)
 %showP(P)
 %%%%%%%%%%%%%%%%%%%%%%%%%%
 %%
 upind=[];
 dnind=P(2,:);
 ppind=pind([1,7],:);
 ppind=ppind(:);
 P(1,:)=arrangemethod4(upind,ppind,dnind);
 %showP(P(1:2,:))upind=P(6,:);
 dnind=P(8,:);
 ppind=pind([1,7],:);
 ppind=ppind(:);P(7,:)=arrangemethod4(upind,ppind,dnind);
 %showP(P(6:8,:))
 %show2line(ppind(1:19)-1,ppind(20:38)-1)
 %%%%%%%%%%%%%%%%%%%5
 showP(P)
 toc(a)      

3.仿真结论:

d3