24小时热门版块排行榜    

查看: 2595  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 3/150 2026-11-09 21:05 by FChKX8meu1ds
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
请填处理意见