CN103824301A - 一种数字散斑相关方法中测度函数改进方法 - Google Patents
一种数字散斑相关方法中测度函数改进方法 Download PDFInfo
- Publication number
- CN103824301A CN103824301A CN201410088944.9A CN201410088944A CN103824301A CN 103824301 A CN103824301 A CN 103824301A CN 201410088944 A CN201410088944 A CN 201410088944A CN 103824301 A CN103824301 A CN 103824301A
- Authority
- CN
- China
- Prior art keywords
- content
- octagon
- subdomain
- target
- respect
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000002596 correlated effect Effects 0.000 claims abstract description 17
- 238000013507 mapping Methods 0.000 claims abstract description 8
- 230000000875 corresponding effect Effects 0.000 claims description 33
- 230000006835 compression Effects 0.000 claims description 14
- 238000007906 compression Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 5
- 230000003287 optical effect Effects 0.000 abstract description 4
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000011524 similarity measure Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明涉及光学测量、图像匹配领域,特别是涉及一种数字散斑相关方法中测度函数改进方法。现有的相关测度公式中,图像模板子域f和目标子域g的形状都为正方形且尺度和方向相同,当目标图像中的内容相对于原始图像有变形时,以上述方法在目标图像中截取的任何一个子域g都不可能和模板子域f含有完全相同的内容,这将导致相关测度函数曲面的主峰不够突出并降低图像匹配的精度。本发明将模板子域形状定为正八边形,将目标子域的形状定为与模板子域具有映射关系的八边形,这将使二者尽量包含相同的内容,并使相关测度函数曲面的主峰更突出,主峰中心更尖锐,明显提高图像配准精度。
Description
技术领域
本发明涉及光学测量、图像匹配领域,特别是涉及一种数字散斑相关方法中测度函数改进方法。
背景技术
图像匹配(图像识别)大体上可分为两种情况:基于对象的识别和非对象识别。所谓对象是一些像素集合,该集合具有特定的图像特征(灰度变异、轮廓、纹理等)。基于对象的识别是相对容易的,其匹配测度函数一般都具有一个陡峭的单峰。但有时需要做比对的图像中没有明显的对象或由于图像过于模糊(或模板过小)使模板中无法定义足以识别的对象,例如岩石力学实验中试件的表面图像由一些随机的灰度斑点组成(一般称为散斑图)。对两幅含有相同内容的散斑图进行匹配时,由于单独的散斑不含有足以分辨的信息量,无法视作对象加以搜索,只有通过增大比对窗口,使窗口内的散斑增多,直到其具有足以分辨的信息量才能实现匹配。这种匹配是一种非对象识别。在光学测量领域,图像匹配函数一般采用某种相关测度计算,所以将这种图像匹配方法称作数字散斑相关方法-DSCM(Digital SpeckleCorrelation Method)。
基于区域灰度的图像匹配可作如下数学描述:
设I为原始图像,f为其中截取的图像模板。T为目标图像,g为其子域。在将模板和目标子域比较时,对应点处的灰度值分别为f(i,j)和g(i+u,j+v),(u,v)为g与f的像间相对位移矢量,相应可将g表示为g(u,v)。根据某相关算法,可建立图像间的相关测度。计算g与f的相关测度值,记作C(u,v)。基于图像灰度相关的测度有多种算法,针对具体情况,采用合适的算法,可以使相关测度函数C(u,v)获得较好的图形表现。若某相关测度函数曲面有一个明显的单峰,其峰值所对应的子域与模板f恰好包含相同内容。
基于区域灰度的各种配准方法,其差异主要表现在相似度量准则的不同,基于不同的度量准则,发展出了很多相关公式,金观昌在《计算机辅助光学测量(第2版)》(北京:清华大学出版社,2007)第146页将现有的大多数相关公式列举如下:
C1=-∑∑|f-g|
C2=-∑∑|f-〈f〉-(g-〈g〉)|
C3=-〈|f-g|〉
C4=-〈(f-g)2〉
C5=〈(f·g)〉
以上各式中,f的g分别表示在原始图像中截取的图像模板和目标图像中截取的目标子域,〈f〉和〈g〉分别表示f和g的均值。
以上各式中,f和g两个子域的形状都为正方形且尺度和方向相同。当目标图像中的内容相对于原始图像有变形时,以上述方法在目标图像中截取的任何一个子域g都不可能和模板子域f含有完全相同的内容,这将导致相关测度函数曲面的主峰不够突出并降低图像匹配的精度。
发明内容
本发明提供了一种数字散斑相关方法中测度函数改进方法,通过将模板子域形状定为正八边形,将目标子域的形状定为与模板子域具有映射关系的八边形,使二者尽量包含相同的内容,改善相关测度函数曲面的主峰范围内的曲面形态,提高配准精度。
为实现上述目的,本发明所采用的技术方案是:一种数字散斑相关方法中测度函数改进方法,包括如下步骤:
(1)设I为原始图像,T为目标图像,f为I中截取的模板子域,g为T中截取的目标子域;
(2)将模板子域f的形状取为正八边形,目标子域g的形状取为和f具有映射关系的八边形;
(3)根据f与g的映射关系,用插值的方法求取g与f的节点网格对应的网格节点处的图像灰度值,将该组插值获得的图像灰度值与f的灰度数据比较,得到相关测度。
进一步地,所述的步骤(2)中,获取目标子域g的形状的方法为:初步计算或预测T中内容相对于I中内容的变形规律,根据这个变形规律对对应于f的正八边形进行变形处理,变形处理后获得的八边形作为T中目标子域g的形状。
所述的根据这个变形规律对对应于f的正八边形进行变形处理,进一步包括:
(1)若T中内容相对于I中内容无任何变形,目标子域g的形状为正八边形,其大小和方向与f完全相同;
(2)若T中内容相对于I中内容仅有横向拉伸或压缩,目标子域g的形状为所述的对应于f的正八边形经同样横向拉伸或压缩后获得的八边形;
(3)若T中内容相对于I中内容仅有纵向拉伸或压缩,目标子域g的形状为所述的对应于f的正八边形经同样纵向拉伸或压缩后获得的八边形;
(4)若T中内容相对于I中内容仅有刚性旋转,目标子域g的形状为所述的对应于f的正八边形经同样刚性旋转后获得的正八边形;
(5)若T中内容相对于I中内容仅有横向剪切变形,目标子域g的形状为所述的对应于f的正八边形经同样横向剪切变形后获得的八边形;
(6)若T中内容相对于I中内容仅有纵向剪切变形,目标子域g的形状为所述的对应于f的正八边形经同样纵向剪切变形后获得的八边形;
(7)若T中内容相对于I中内容包含以上两种以上的变形,目标子域g的形状为所述的对应于f的正八边形经同样混合变形后获得的八边形。
进步地,求取g与f的节点网格对应的网格节点处的图像灰度值的插值算法包括:
(1)双线性插值法;
(2)双立方插值法;
(3)双三次样条插值法。
本发明达到的有益效果:使相关测度函数曲面的主峰更突出,主峰中心更尖锐,明显提高图像配准精度。
附图说明
图1原始图像和目标图像的关系;
图2正八边形的目标子域;
图3经横向拉伸或压缩的目标子域;
图4经纵向拉伸或压缩的目标子域;
图5经刚体旋转的目标子域;
图6经横向剪切变形的目标子域;
图7经纵向剪切变形的目标子域。
具体实施方式
下面结合附图对本发明做进一步详细说明。
本发明的一种数字散斑相关方法中测度函数改进方法,包括如下步骤:
(1)设I为原始图像,f为其中截取的图像模板子域,T为目标图像,g为其目标子域,在将模板子域f和目标子域g比较时,对应点处的灰度值分别为f(i,j)和g(i+u,j+v),(u,v)为g与f的像间相对位移矢量,相应可将g表示为g(u,v),根据某相关算法,可建立图像间的相关测度C,计算g与f的相关测度,记作C(u,v),当(u,v)为变量时,函数C(u,v)称为相关测度函数,简称测度函数,现有的测度函数算法获取模板子域f和目标子域g的方法都为截取正方形子域,如图1所示。将这些测度函数算法做以下改进:将模板子域f的形状取为正八边形,目标子域g的形状都取为和f具有映射关系的八边形。
(2)根据f与g的映射关系,用插值的方法求取g的与f的节点网格对应的网格节点处的图像灰度值,将该组插值获得的图像灰度值与f的数据比较计算相关测度。
进一步地,获取目标子域g的形状的方法为:初步计算或预测T中内容相对于I中内容的变形规律,根据这个变形规律对对应于f的正八边形进行变形处理,变形处理后获得的八边形作为T中目标子域g的形状,进行的变形处理进一步包括:
(1)若T中内容相对于I中内容无任何变形,目标子域g的形状为对应于f的正八边形,其大小方向与f完全相同,如图2所示;
(2)若T中内容相对于I中内容仅有横向拉伸或压缩,目标子域g的形状为对应于f的正八边形经同样横向拉伸或压缩后获得的八边形,如图3所示;
(3)若T中内容相对于I中内容仅有纵向拉伸或压缩,目标子域g的形状为对应于f的正八边形经同样纵向拉伸或压缩后获得的八边形,如图4所示;
(4)若T中内容相对于I中内容仅有刚性旋转,目标子域g的形状为对应于f的正八边形经同样刚性旋转后获得的正八边形,如图5所示;
(5)若T中内容相对于I中内容仅有横向剪切变形,目标子域g的形状为对应于f的正八边形经同样横向剪切变形后获得的八边形,如图6所示;
(6)若T中内容相对于I中内容仅有纵向剪切变形,目标子域g的形状为对应于f的正八边形经同样纵向剪切变形后获得的八边形,如图7所示;
(7)若T中内容相对于I中内容包含以上两种以上的变形,目标子域g的形状为对应于f的正八边形经同样混合变形后获得的八边形;
进一步地,求取g的与f的节点网格对应的网格节点处的图像灰度值的插值算法包括:
(1)双线性插值法;
(2)双立方插值法;
(3)双三次样条插值法。
Claims (4)
1.一种数字散斑相关方法中测度函数改进方法,其步骤为:
(1)设I为原始图像,T为目标图像,f为I中截取的模板子域,g为T中截取的目标子域;
(2)将模板子域f的形状取为正八边形,目标子域g的形状取为和f具有映射关系的八边形;
(3)根据f与g的映射关系,用插值的方法求取g与f的节点网格对应的网格节点处的图像灰度值,将该组插值获得的图像灰度值与f的灰度数据比较,得到相关测度。
2.根据权利要求1所述的数字散斑相关方法中测度函数改进方法,其特征在于,所述的步骤(2)中,获取目标子域g的形状的方法为:初步计算或预测T中内容相对于I中内容的变形规律,根据这个变形规律对对应于f的正八边形进行变形处理,变形处理后获得的八边形作为T中目标子域g的形状。
3.根据权利要求2所述的数字散斑相关方法中测度函数改进方法,其特征在于,所述的根据这个变形规律对对应于f的正八边形进行变形处理,进一步包括:
(1)若T中内容相对于I中内容无任何变形,目标子域g的形状为正八边形,其大小和方向与f完全相同;
(2)若T中内容相对于I中内容仅有横向拉伸或压缩,目标子域g的形状为所述的对应于f的正八边形经同样横向拉伸或压缩后获得的八边形;
(3)若T中内容相对于I中内容仅有纵向拉伸或压缩,目标子域g的形状为所述的对应于f的正八边形经同样纵向拉伸或压缩后获得的八边形;
(4)若T中内容相对于I中内容仅有刚性旋转,目标子域g的形状为所述的对应于f的正八边形经同样刚性旋转后获得的正八边形;
(5)若T中内容相对于I中内容仅有横向剪切变形,目标子域g的形状为所述的对应于f的正八边形经同样横向剪切变形后获得的八边形;
(6)若T中内容相对于I中内容仅有纵向剪切变形,目标子域g的形状为所述的对应于f的正八边形经同样纵向剪切变形后获得的八边形;
(7)若T中内容相对于I中内容包含以上两种以上的变形,目标子域g的形状为所述的对应于f的正八边形经同样混合变形后获得的八边形。
4.根据权利要求1所述的数字散斑相关方法中测度函数改进方法,其特征在于,求取g与f的节点网格对应的网格节点处的图像灰度值的插值算法包括:
(1)双线性插值法;
(2)双立方插值法;
(3)双三次样条插值法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410088944.9A CN103824301A (zh) | 2014-03-12 | 2014-03-12 | 一种数字散斑相关方法中测度函数改进方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410088944.9A CN103824301A (zh) | 2014-03-12 | 2014-03-12 | 一种数字散斑相关方法中测度函数改进方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103824301A true CN103824301A (zh) | 2014-05-28 |
Family
ID=50759342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410088944.9A Pending CN103824301A (zh) | 2014-03-12 | 2014-03-12 | 一种数字散斑相关方法中测度函数改进方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103824301A (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102722890A (zh) * | 2012-06-07 | 2012-10-10 | 内蒙古科技大学 | 基于光流场模型的非刚性心脏图像分级配准方法 |
CN102968788A (zh) * | 2012-10-25 | 2013-03-13 | 武汉大学 | 一种基于规则格网面元的波段配准方法 |
CN103534726A (zh) * | 2011-05-17 | 2014-01-22 | 苹果公司 | 用于全景摄影的位置传感器辅助的图像配准 |
CN103544711A (zh) * | 2013-11-08 | 2014-01-29 | 国家测绘地理信息局卫星测绘应用中心 | 遥感影像的自动配准方法 |
CN103593675A (zh) * | 2013-11-25 | 2014-02-19 | 中国人民解放军国防科学技术大学 | 一种基于对数极坐标变换的静脉匹配方法 |
-
2014
- 2014-03-12 CN CN201410088944.9A patent/CN103824301A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103534726A (zh) * | 2011-05-17 | 2014-01-22 | 苹果公司 | 用于全景摄影的位置传感器辅助的图像配准 |
CN102722890A (zh) * | 2012-06-07 | 2012-10-10 | 内蒙古科技大学 | 基于光流场模型的非刚性心脏图像分级配准方法 |
CN102968788A (zh) * | 2012-10-25 | 2013-03-13 | 武汉大学 | 一种基于规则格网面元的波段配准方法 |
CN103544711A (zh) * | 2013-11-08 | 2014-01-29 | 国家测绘地理信息局卫星测绘应用中心 | 遥感影像的自动配准方法 |
CN103593675A (zh) * | 2013-11-25 | 2014-02-19 | 中国人民解放军国防科学技术大学 | 一种基于对数极坐标变换的静脉匹配方法 |
Non-Patent Citations (2)
Title |
---|
张军辉等: "《数字散斑相关法应用于非接触应变的测量方法》", 《理化检验-物理分册》 * |
王伟等: "《利用数字图像相关性实现岩土自动跟踪》", 《应用科学学报》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110443836B (zh) | 一种基于平面特征的点云数据自动配准方法及装置 | |
CN103727930B (zh) | 一种基于边缘匹配的激光测距仪与相机相对位姿标定方法 | |
EP3091479B1 (en) | Method and apparatus for fingerprint identification | |
CN102402685B (zh) | 基于Gabor特征的三马尔可夫场SAR图像分割方法 | |
CN101556695A (zh) | 一种图像匹配方法 | |
CN101950364A (zh) | 基于邻域相似度和阈值分割的遥感图像变化检测方法 | |
CN102360503B (zh) | 基于空间贴近度和像素相似性的sar图像变化检测方法 | |
CN102903109B (zh) | 一种光学影像和sar影像一体化分割配准方法 | |
CN104361590A (zh) | 一种控制点自适应分布的高分辨率遥感影像配准方法 | |
CN103886606B (zh) | 基于联合广义伽玛分布参数的sar图像分割方法 | |
CN105160686A (zh) | 一种基于改进sift算子的低空多视角遥感影像匹配方法 | |
CN103871039A (zh) | 一种sar图像变化检测差异图生成方法 | |
CN111325184B (zh) | 一种遥感影像智能解译与变化信息检测方法 | |
CN110704559A (zh) | 一种多尺度矢量面数据匹配方法 | |
CN107392953B (zh) | 基于等高线的深度图像识别方法 | |
CN105488798B (zh) | 基于点集对比的sar图像相似性度量方法 | |
CN117274627A (zh) | 一种基于图像转换的多时相积雪遥感图像匹配方法及系统 | |
CN101894368B (zh) | 基于测地距离的半监督sar图像水域分割方法 | |
CN103325104B (zh) | 基于迭代稀疏表达的人脸图像超分辨率重建方法 | |
CN117474765B (zh) | 基于参考影像纹理转移的dem超分辨率重建系统 | |
CN103310456B (zh) | 基于Gaussian-Hermite矩的多时相/多模态遥感图像配准方法 | |
CN104200460A (zh) | 基于图像特征和互信息的图像配准方法 | |
CN111080649A (zh) | 一种基于黎曼流形空间的图像分割处理方法及系统 | |
CN103854285A (zh) | 基于随机投影和改进谱聚类的sar图像地物分割方法 | |
CN103824301A (zh) | 一种数字散斑相关方法中测度函数改进方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140528 |
|
WD01 | Invention patent application deemed withdrawn after publication |