24小时热门版块排行榜

返回列表

472381565

新虫 (小有名气)

应助: 3 (幼儿园)
金币: 607.6
散金: 14
帖子: 133
在线: 19.5小时
虫号: 2082313
注册: 2012-10-24
性别: MM
专业: 遥感机理与方法

[求助] 给这段MATLAB程序中没加注解的地方加上注解。

I1 = double(imread('X6.jpg'));
[h1 w1 d1] = size(I1);
I2 = double(imread('X7.jpg'));
[h2 w2 d2] = size(I2);
%显示图像并在显示的图像上取坐标

figure; subplot(1,2,1); image(I1/255); axis image; hold on;
title('first input image');
[X1 Y1] = ginput2(2); %在第一幅图像上确定2个点
subplot(1,2,2); image(I2/255); axis image; hold on;
title('second input image');
[X2 Y2] = ginput2(2); %在第二幅图像上取2个点

Z  = [ X2'  Y2' ; Y2' -X2' ; 1 1 0 0  ; 0 0 1 1 ]';
xp = [ X1 ; Y1 ];
t  = Z \ xp;
a  = t(1);
b  = t(2);
tx = t(3);
ty = t(4);

T = [a b tx ; -b a ty ; 0 0 1];

cp = T*[ 1 1 w2 w2 ; 1 h2 1 h2 ; 1 1 1 1 ];
Xpr = min( [ cp(1,

0 ] ) : max( [cp(1,

w1] );
Ypr = min( [ cp(2,

0 ] ) : max( [cp(2,

h1] );
[Xp,Yp] = ndgrid(Xpr,Ypr);
[wp hp] = size(Xp);

X = T \ [ Xp(

Yp(

ones(wp*hp,1) ]';

clear Ip;
xI = reshape( X(1,

,wp,hp)';
yI = reshape( X(2,

,wp,hp)';
Ip(:,:,1) = interp2(I2(:,:,1), xI, yI, '*bilinear');
Ip(:,:,2) = interp2(I2(:,:,2), xI, yI, '*bilinear');
Ip(:,:,3) = interp2(I2(:,:,3), xI, yI, '*bilinear');

offset = -round( [ min( [ cp(1,

0 ] ) min( [ cp(2,

0 ] ) ] );
Ip(1+offset(2):h1+offset(2),1+offset(1):w1+offset(1),

= double(I1(1:h1,1:w1,

);

figure; image(Ip/255); axis image;
title('mosaic image');

回复此楼