24小时热门版块排行榜    

查看: 2421  |  回复: 0

xiashaoyan

新虫 (正式写手)

[求助] 配准图像时,求互相关函数的最大值位置的疑问

看了一些求基准图像和配准图像的互相关值的最大值对应的位置,我现在看到有两种方法,
第一种是使用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的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 xiashaoyan 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
信息提示
请填处理意见