| 查看: 758 | 回复: 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'); |
回复此楼» 猜你喜欢
关于Filecode分析方法
已经有6人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有9人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
FileCode能看出啥？
已经有21人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有35人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
3楼2013-06-02 07:30:51
472381565
新虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 607.6
- 散金: 14
- 帖子: 133
- 在线: 19.5小时
- 虫号: 2082313
- 注册: 2012-10-24
- 性别: MM
- 专业: 遥感机理与方法
2楼2013-06-01 16:42:49
