| 查看: 791 | 回复: 3 | ||
472381565新虫 (小有名气)
|
[求助]
给这段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'); |
回复此楼» 猜你喜欢
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
472381565
新虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 607.6
- 散金: 14
- 帖子: 133
- 在线: 19.5小时
- 虫号: 2082313
- 注册: 2012-10-24
- 性别: MM
- 专业: 遥感机理与方法
2楼2013-06-01 16:42:49
3楼2013-06-02 07:30:51
 |
4楼2013-06-08 00:19:28
投票:
30