24小时热门版块排行榜    

查看: 2493  |  回复: 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,其中的具体含义是什么呢?

两种方法计算最大值对应的位置,我思索了半天想,最后计算的结果表示什么意思,是表示配准时两幅图像的偏差量吗?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cctt126

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
图像做类似卷积运算前,都要考虑边界问题。
图像的四边要‘扩充’配准图像的半径。
所以你得到的最大值位置是扩充后图像的坐标位置,需要减去扩充的尺寸。
在我用c写了图像卷积算法以后,才知道这个道理的。
数学,英语,编程,没一样让我省心。。。
2楼2015-03-10 10:36:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiashaoyan 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 271材料工程求调剂 +6 .6lL 2026-03-18 6/300 2026-03-19 15:41 by haoshis
[考研] 材料专硕274一志愿陕西师范大学求调剂 +8 薛云鹏 2026-03-13 8/400 2026-03-19 15:36 by haoshis
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 一志愿南昌大学,327分,材料与化工085600 +3 Ncdx123456 2026-03-19 3/150 2026-03-19 13:18 by houyaoxu
[考研] 0703化学 305求调剂 +4 FY_yy 2026-03-14 4/200 2026-03-19 05:54 by anny19840123
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
信息提示
请填处理意见