碎片拼接

 

1.问题描述：

碎片拼接

2.部分程序：

%cd piecejoint

a=tic;
%%
global allP
global margl margr
global ld
global pn pw ph

allP=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.仿真结论：

C74