| 查看: 2569 | 回复: 1 | ||
xiashaoyan新虫 (正式写手)
|
[求助]
normxcorr2函数求互相关量最大值位置的问题 已有1人参与
|
|
看了一些求基准图像和配准图像的互相关值的最大值对应的位置，我现在看到有两种方法， 第一种是使用MATLAB里面自带的normxcorr2函数，我看了mathwork网站上有一个求互相关量最大值位置的程序：(参考网址为：http://cn.mathworks.com/help/images/ref/normxcorr2.html) onion = rgb2gray(onion); peppers = rgb2gray(peppers); Perform cross-correlation and display result as surface. c = normxcorr2(onion,peppers); [ypeak, xpeak] = find(c==max(c(  );yoffSet = ypeak-size(onion,1); xoffSet = xpeak-size(onion,2); 可是有一个疑问是，为什么在求出互相关量最大值位置后还要减去基准图像的尺寸呢？(也即为什么yoffSet = ypeak-size(onion,1)？？xoffSet = xpeak-size(onion,2);？？) 第二种方法是采用两副图像进行互相关运算，计算的程序，我找了网上的程序为： (参考网址为：http://www.mathworks.com/matlabc ... /dftregistration.m) buf1ft是基准图像，buf2ft是配准待图像 [m,n]=size(buf1ft); CC = ifft2(buf1ft.*conj(buf2ft));%先计算buf1ft与buf2ft共轭的傅立叶变换，得到互相关函数的频谱，然后求结果ifft2得到互相关函数 %% 计算互相关函数的最大值及最大值的位置 [max1,loc1] = max(CC); [max2,loc2] = max(max1); rloc=loc1(loc2); cloc=loc2; CCmax=CC(rloc,cloc); md2 = fix(m/2); nd2 = fix(n/2); if rloc > md2 row_shift = rloc - m - 1; else row_shift = rloc - 1; end if cloc > nd2 col_shift = cloc - n - 1; else col_shift = cloc - 1; end 这里我同样有一个，就是为什么求出互相关量最大值的位置后还要减去一个1或者减去1+基准图像的尺寸呢？(也即：row_shift = rloc - m - 1？？或者row_shift = rloc - 1)，还有为什么rloc > md2时，就得 row_shift = rloc - m - 1，其中的具体含义是什么呢？ 两种方法计算最大值对应的位置，我思索了半天想，最后计算的结果表示什么意思，是表示配准时两幅图像的偏差量吗？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
国社科又开始会评了，不知道这次命运如何
已经有13人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有13人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
2楼2015-03-10 10:36:26
30