function I = Stitching_jq(I1,I2,H)
[height1,width1,c1]=size(I1);
[height2,width2,c2]=size(I2);
if c1~=c2
error('Image Channel Error!');
end
N=floor(H(2)*height1+H(5)*width1+H(8)-15);
for j=1:width1-N
for i=1:height1
I(i,j,:)=I1(i,j,:);
end;
end;
for j=width1-N+1:width1
for i=1:height1
for k=1:c1
x=H(1)*i+H(4)*j+H(7);
y=H(2)*i+H(5)*j+H(8);
x1=floor(x);
y1=floor(y);
x2=x1+1;
y2=y1+1;
if x<1|x>=height2|y<1|y>=width2
continue;
else
P1=I2(x1,y1,k);
P2=I2(x1,y2,k);
P3=I2(x2,y1,k);
P4=I2(x2,y2,k);
B=(1-x+x1)*(1-y+y1)*P1+(x-x1)*(1-y+y1)*P3+(1-x+x1)*(y-y1)*P2+(x-x1)*(y-y1)*P4;
d=width1-j;
I(i,j,k)=d/N*I1(i,j,:)+(1-d/N)*B;
end
end;
end;
end;
for j=width1+1:width1+width2
for i=1:height1+height2
for k=1:c2
x=H(1)*i+H(4)*j+H(7);
y=H(2)*i+H(5)*j+H(8);
x1=floor(x);
y1=floor(y);
x2=x1+1;
y2=y1+1;
if x<1|x>=height2|y<1|y>=width2
continue;
else
P1=I2(x1,y1,k);
P2=I2(x1,y2,k);
P3=I2(x2,y1,k);
P4=I2(x2,y2,k);
I(i,j,k) =(1-x+x1)*(1-y+y1)*P1+(x-x1)*(1-y+y1)*P3+(1-x+x1)*(y-y1)*P2+(x-x1)*(y-y1)*P4;
end;
end
end;
end
figure; imshow(uint8(I));
评论0