CN104637041B - 一种基于参照特征的宽幅织物图像获取与拼接方法 - Google Patents

一种基于参照特征的宽幅织物图像获取与拼接方法 Download PDF

Info

Publication number
CN104637041B
CN104637041B CN201410854510.5A CN201410854510A CN104637041B CN 104637041 B CN104637041 B CN 104637041B CN 201410854510 A CN201410854510 A CN 201410854510A CN 104637041 B CN104637041 B CN 104637041B
Authority
CN
China
Prior art keywords
image
reference picture
fabric
harris
pixel
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.)
Active
Application number
CN201410854510.5A
Other languages
English (en)
Other versions
CN104637041A (zh
Inventor
田秋红
陈本永
孙政荣
姚佳宝
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Guangyu New Material Technology Co Ltd
Original Assignee
Zhejiang Sci Tech University ZSTU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang Sci Tech University ZSTU filed Critical Zhejiang Sci Tech University ZSTU
Priority to CN201410854510.5A priority Critical patent/CN104637041B/zh
Publication of CN104637041A publication Critical patent/CN104637041A/zh
Application granted granted Critical
Publication of CN104637041B publication Critical patent/CN104637041B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Analysis (AREA)

Abstract

本发明公开了一种基于参照特征的宽幅织物图像获取与拼接方法。将织物水平放织物载台上,其中间水平放上参照织物,再移走参照织物,过程中分别获得左、右图参照图像和原始左、右图图像;将左、右图参照图像采用多尺度Harris角点检测方法得到各自Harris角点集;采用归一化的相关测度方法进行初始匹配,得到初始匹配点集;采用RANSAC方法进行初始匹配特征点求精和单应矩阵求解,得到最终单应矩阵;对原始左、右图图像进行投影变换配准,然后拼接融合,完成织物图像的获取与拼接。本发明实现循环周期强或无明显特征的织物图像的无缝拼接,拼接精度高,可精确计算出待拼接原始织物图像的相对几何位置关系,具有较强的抗干扰能力。

Description

一种基于参照特征的宽幅织物图像获取与拼接方法
技术领域
本发明涉及一种图像处理方法,尤其是涉及一种基于参照特征的宽幅织物图像获取与拼接方法,循环周期性较强或无明显特征。
背景技术
在研究纺织品织物特性过程中,如织物悬垂性、织物褶皱性、织物有无疵点等,受硬件条件限制使得一次拍摄只能实现局部图像的采集,特别是对于宽门幅织物。因此,通常需要采用双目相机分别采集两幅具有一定重叠区域的局部图像,再进行无缝拼接,形成一幅完整的织物图像,进而实现织物特性研究,故图像拼接技术是织物特征研究的首要前提。
目前图像拼接方法主要有模板匹配、基于网格的匹配和基于特征的匹配。模板匹配先在参考图像的重叠部分设置模板,用此模板在目标图像中搜索具有相同(或相似)值的对应区域,从而确定重叠范围,实现图像拼接。该方法在重叠部分足够大时,匹配精度较高,但计算量也随着模板的增大而迅速增大,很难满足在线实时拼接的需要,且易受光照、变形等因素的影响;基于网格匹配是在搜索过程中进行粗略和精确两步匹配。该方法虽然在运算速度上较模板匹配有所改善,但是在粗略匹配过程中,移动的步长难以确定,有可能将所取网格在目标图像匹配过程中被划分开,造成匹配中无法取出与原网格完全匹配的最佳网格,难以达到精确匹配;基于特征的匹配首先提取两帧图像重叠局部的特征,然后用相似性准则匹配特征。其中Harris角点提取方法由于不受摄像机姿态及光照的影响,是目前效果最好的方法。但该方法不适用于图像内容单一特征较少的情况。
发明内容
针对织物特性测试等技术领域需要高精度的织物图像拼接,本发明的目的在于基于参照特征的宽幅织物图像获取与拼接方法,利用参照织物提取Harris角点和计算单应矩阵完成原始待拼接图像的高精度无缝拼接,是一种非接触式的织物图像自动拼接方法,解决了现有拼接方法针对循环周期性较强或没有明显特征的织物图像的拼接失真问题。
本发明解决其技术问题所采用的技术方案步骤如下:
1)先将被拼接织物水平放在织物载台上,获取占有被拼接织物左侧三分之二的图像,将该图像进行归一化处理得到被拼接织物的原始左图图像I1
2)在被拼接织物的中间水平放上参照织物,获取置有参照织物的左图参照图像I1',左图参照图像I1'与原始左图图像I1大小相同;水平向左移动织物载台,获取置有参照织物的右图参照图像I2',并使得左图参照图像I1'和右图参照图像I2'共同有至少参照织物三分之一的重叠区域;
3)再保持被拼接织物位置不变,移走参照织物,获取被拼接织物的原始右图图像I2,原始右图图像I2与右图参照图像I2'大小相同;
4)对于左图参照图像I1'和右图参照图像I2',均采用多尺度Harris角点检测方法分别得到左图参照图像I1'的左图Harris角点集S1和右图参照图像I2'的右图Harris角点集S2
5)采用归一化的相关测度方法对两个角点集S1和S2进行Harris角点的初始匹配,得到初始匹配点集S;
6)采用RANSAC方法对初始匹配点集S进行初始匹配特征点求精和单应矩阵的求解,得到初始匹配特征点求精结果集合Si’,并得到最终单应矩阵H';
7)根据得到的最终单应矩阵H',对原始左图图像I1和原始右图图像I2进行投影变换配准,再拼接并利用加权平均方法消除拼接缝,进行图像的融合,完成织物图像的获取与拼接。
所述的步骤4)中左图参照图像I1'和右图参照图像I2'的两个参照图像均采用以下多尺度Harris角点检测方法进行处理:
4.1)图像的尺度空间按高斯窗的尺度δm采用以下公式1进行计算,尺度δm在0.5~5的范围内以0.5的步长变化并依次计算:
式中,Ix'、Iy'分别为参照图像在x方向和y方向的偏导,G(x,y,δm)为高斯核,是尺度变换的唯一变换核,m为尺度计算的序数,Lxm)、Lym)、Lxym)分别表示x方向、y方向和xy方向的高斯尺度空间;
4.2)由上述步骤4.1)的公式1计算得到不同尺度δm上的矩阵M:
以及采用以下公式3计算得到在不同尺度δm上的响应函数C(x,y):
C(x,y)=Det[M(x,y)]-βTrace2[M(x,y)] (3)
式中,Det[M(x,y)]为矩阵M的行列式,Trace[M(x,y)]为矩阵M的迹,β是常数因子;
4.3)如果响应函数C(x,y)大于响应函数阈值T1,则响应函数C(x,y)所对应的像素点(x,y)为角点;为避免在不同尺度下手动调整阈值,响应函数阈值T1可由图像中可能存在的角点总数n0代替;
4.4)对于尺度δm以0.5的步长在0.5~5的范围内变化的每个尺度,采用非极大抑制法抑制多响应的影响,非极大抑制窗的尺寸为3×3,将响应函数C(x,y)的值按照从大到小降序排列,将各个尺度下的左图参照图像I1'和右图参照图像I2'中可能存在的角点作为候选角点;
4.5)对于尺度δm以0.5的步长在0.5~5的范围内变化的第二个尺度δm=1开始到最大尺度δm=5,依次遍历各个尺度,查找是否有当前尺度下的候选角点在前一尺度中出现,如果没有,保留该候选角点;否则,剔除候选角点;由此得到左图参照图像I1'的左图Harris角点集S1和右图参照图像I2'的右图Harris角点集S2
所述的步骤5)的左图参照图像I1'和右图参照图像I2'的两个参照图像具体采用以下方法进行初始匹配:
5.1)在左图参照图像I1'中,以左图Harris角点集S1的任一Harris角点X1(x1,y1)为中心,使得能建立大小为3×3的搜索窗口;
5.2)对于左图参照图像I1'和右图参照图像I2',均采用以下公式4依次遍历各个Harris角点分别计算两个参照图像灰度值的标准差stdd(Xk):
式中:为参照图像中Harris角点Xk(x,y)处的灰度平均值,Ik'(x+i,y+i)表示参照图像中像素点Xk(x+i,y+i)处的灰度值,k为左图Harris角点集S1和右图Harris角点集S2中的Harris角点序数,Xk为左图Harris角点集S1和右图Harris角点集S2中的Harris角点,i、j分别表示参照图像的x方向、y方向的像素点序数,N为参照图像的像素点总数;
5.3)分别遍历两个参照图像中的各个像素点,采用以下公式5计算两幅参照图像之间的协方差:
其中,X1和X2分别为左图Harris角点集S1和右图Harris角点集S2中的Harris角点;
5.4)采用以下公式6的归一化自相关函数cor(X1,X2)进行归一化自相关计算:
5.5)当归一化自相关函数cor(X1,X2)大于自相关阈值T2时,则该归一化自相关函数cor(X1,X2)对应的Harris角点为候选匹配点,由所有候选匹配点组成初始匹配特征点集S。
所述的步骤6)具体过程如下:
6.1)将初始匹配特征点集S里面的匹配点对作为样本,根据置信概率α和所有样本的数量计算得到迭代次数t,采用以下步骤6.2)~6.4)进行循环迭代计算;
6.2)从初始特征匹配点集S中随机选取l对样本点对,利用这2l个匹配点采用以下公式7计算得到两幅参照图像之间的临时单应矩阵Hv
式中,(Xv,Yv)是左图参照图像I1'中的Harris角点,(xv,yv)是右图参照图像I2'中的Harris角点;v为当前循环次数的序数,第一次计算的序数v为1;
6.3)由临时单应矩阵Hv采用最小二乘法,采用以下公式8计算初始特征匹配点集S剩余各个匹配点对的极线距离dv
若极线距离dv小于极线距离阈值T3,则将该极线距离dv所对应的匹配点对加入第v匹配点集Sv,该极线距离dv加入第v距离集合Dv
6.4)将当前循环次数的序数v加1,若加1后的循环次数的序数<t,t为迭代次数,则跳转到步骤6.2),否则进行下一步骤6.5);
6.5)选取所有匹配点集中匹配点对最多的匹配点集作为初始匹配特征点求精结果集合Sv’,所对应的临时单应矩阵Hv作为最终单应矩阵H'。
所述的步骤7)具体过程如下:建立纯黑色的背景图像I12,背景图像I12各个像素点的RGB像素值均为(255,255,255),计算最终单应矩阵H'和原始右图图像I2的乘积,将原始右图图像I2投影变换到背景图像I12中,背景图像I12中重叠区域外对应像素值按左图和右图的像素值不变,重叠区域处的像素值采用公式9的加权平均方法进行图像融合,由此获得完整的拼接织物图像:
式中,I1(x,y)为原始左图图像I1重叠区域的像素值;I2(x,y)为原始右图图像I2重叠区域的像素值;d1为像素点(x,y)距原始右图图像左边缘的垂直距离;d2为像素点(x,y)距原始左图图像右边缘的垂直距离;I12(x,y)表示融合后的像素点(x,y)的像素值。
所述的步骤6.1)中迭代次数t采用以下公式10计算得到:
t=wα (10)
式中,w为初始匹配特征点集S中特征点数量,α为置信概率。
所述方法用于对宽幅织物图像的获取与拼接。
本发明具有的有益效果是:
(1)本发明采用Harris角点检测和RANSAC方法,实现循环周期强或无明显特征的织物图像的无缝拼接,拼接精度高,进而保证宽幅织物图像性能的准确测试。
(2)依据参照特征得到的Harris角点和单应矩阵,可精确计算出待拼接原始织物图像的相对几何位置关系,具有较强的抗干扰能力,成本低,易于实用化。
附图说明
图1为本发明方法的流程图。
图2为本发明方法实施装置结构框图。
图3为本发明实施例的原始左图图像和原始右图图像。
图4为本发明实施例的左图参照图像和右图参照图像。
图5为本发明实施例的Harris角点检测效果图。
图6为本发明实施例的最终拼接效果图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
本发明方法的实施装置如图2所示,宽幅织物放在织物载台上,CCD相机位于宽幅织物的正上方并朝向宽幅织物,CCD相机经图像采集卡连接计算机,织物载台连接有用于控制左右移动的步进电机,步进电机经步进电机控制单元连接PC机。
本发明图像处理的坐标系是以图像的左上角点为原点,横向为x方向,纵向为y方向,斜对角为xy方向。
如图1所示,本发明方法的实施例如下:
1)被拼接织物采用44英寸的宽幅织物,先将被拼接织物水平放在织物载台上,利用CCD相机(型号:OK_AC1300)获取占有被拼接织物左侧三分之二的图像,经图像采集卡(型号:OK_RGB20B)输入计算机(型号:Lenovo s520),将该图像进行归一化处理得到宽高为800p×530p的被拼接织物的原始左图图像I1,如图3(a)。
2)在被拼接织物的中间水平放上参照织物,利用CCD相机获取置有参照织物的左图参照图像I1',如图4(a),左图参照图像I1'与原始左图图像I1大小相同,宽高均为800p×530p。
由步进电机控制单元驱动织物载台,水平向左移动织物载台,利用CCD相机获取置有参照织物的右图参照图像I2',如图4(b),并使得左图参照图像I1'和右图参照图像I2'共同有具有参照织物三分之一的重叠区域。
3)再保持被拼接织物位置不变,移走参照织物,利用CCD相机获取被拼接织物的原始右图图像I2,如图3(b),原始右图图像I2与右图参照图像I2'大小相同。
4)对于左图参照图像I1'和右图参照图像I2',均采用多尺度Harris角点检测方法分别得到左图参照图像I1'的左图Harris角点集S1和右图参照图像I2'的右图Harris角点集S2
4.1)图像的尺度空间按高斯窗的尺度δm采用公式1进行计算,尺度δm在0.5~5的范围内以0.5的步长变化并依次计算:
4.2)由公式1计算得到不同尺度δm上的矩阵M,并采用公式3计算得到在不同尺度δm上的响应函数C(x,y),其中常数因子β一般取值在0.05~0.5之间,常数因子实施例取值为0.4。
4.3)如果响应函数C(x,y)大于响应函数阈值T1,则响应函数C(x,y)所对应的像素点(x,y)为角点;实施例的响应函数阈值T1由图像中可能存在的角点总数n0代替;
4.4)对于尺度δmδi以0.5的步长在0.5~5的范围内变化的每个尺度,尺度共计十个,分别为0.5、1、1.5、2、2.5、3、3.5、4、4.5和5,采用非极大抑制法抑制多响应的影响,非极大抑制窗的尺寸为3×3,将得到的响应函数C(x,y)的值按照从大到小降序排列,将各个尺度下的左图参照图像I1'和右图参照图像I2'中可能存在的角点作为候选角点;
4.5)对于尺度δm以0.5的步长在0.5~5的范围内变化的第二个尺度δm=1开始到最大尺度δm=5,依次遍历各个尺度,查找是否有当前尺度下的候选角点在前一尺度中出现,如果没有,保留该候选角点;否则,剔除候选角点;这一步保证了在尽量小的尺度下提取角点,提高了角点的定位精度;
由此得到左图参照图像I1'的左图Harris角点集S1和右图参照图像I2'的右图Harris角点集S2,如图5(a)和5(b)所示。
5)采用归一化的相关测度方法对两个角点集S1和S2进行Harris角点的初始匹配,得到初始匹配点集S。
5.1)在左图参照图像I1'中,以左图Harris角点集S1的任一Harris角点X1(x1,y1)为中心,使得能建立大小为3×3的搜索窗口,并建立以该角点X1(x1,y1)为中心3×3的搜索窗口;
5.2)对于左图参照图像I1'和右图参照图像I2',均采用公式4依次遍历各个Harris角点分别计算两个参照图像灰度值的标准差:
5.3)分别遍历两个参照图像中的各个像素点,采用公式5计算两幅参照图像之间的协方差cov(X1,X2):
5.4)采用公式6的归一化自相关函数cor(X1,X2)进行归一化自相关计算:
5.5)当归一化自相关函数cor(X1,X2)大于自相关阈值T2时,T2实施例取0.8,则该归一化自相关函数cor(X1,X2)对应的Harris角点为候选匹配点,由所有候选匹配点组成初始匹配特征点集S。
6)采用RANSAC方法对初始匹配点集S进行初始匹配特征点求精和单应矩阵的求解,得到初始匹配特征点求精结果集合Si’,并得到最终单应矩阵H'。
6.1)将初始匹配特征点集S里面的匹配点对作为样本,根据置信概率α和所有样本的数量计算得到迭代次数t,置信概率α取α=0.96,采用以下步骤6.2)~6.4)进行循环迭代计算;
上述迭代次数t采用公式10计算得到,实施例中公式10的w取值w=61。
6.2)从初始特征匹配点集S中随机选取l对样本点对,实施例中l为4,故利用这八个匹配点采用公式7计算得到两幅参照图像之间的临时单应矩阵Hv
6.3)由临时单应矩阵Hv采用最小二乘法,采用公式8计算初始特征匹配点集S剩余各个匹配点对的极线距离dv
极线距离阈值T3取T3=2,以为例,依据公式8得其极线距离dv=1.0e+0.006*0.1425dv小于极线距离阈值T3,则将该极线距离dv所对应的匹配点对加入第v匹配点集Sv,该极线距离dv加入第v距离集合Dv
6.4)将当前循环次数的序数v加1,若加1后的循环次数的序数<t,t为迭代次数,实施例中迭代次数t=61×0.96≈59;则跳转到步骤6.2),否则进行下一步骤6.5);
6.5)选取所有匹配点集中匹配点对最多的匹配点集作为初始匹配特征点求精结果集合Sv’,所对应的临时单应矩阵Hv作为最终单应矩阵H',实施例中H'的匹配点对数为37。
7)根据得到的最终单应矩阵H',对原始左图图像I1和原始右图图像I2进行投影变换配准,再拼接并利用加权平均方法消除拼接缝,进行图像的融合,完成织物图像的获取与拼接,实施例中H'的取值为
由于原始左图图像I1和原始右图图像I2的位置关系和置有参照织物的左图参照图像I1'和右图参照图像I2'的位置关系一致,故利用左图参照图像I1'和右图参照图像I2',计算出的单应矩阵H',同时也反映了原始图像间的投影映射关系。
具体建立大小为1600p×530p和RGB像素值为(255,255,255)的背景图像I12,计算最终单应矩阵H'和原始右图图像I2的乘积,将原始右图图像图像I2投影变换到背景图像I12中,背景图像I12中重叠区域外对应像素值按左图和右图的像素值不变,重叠区域处的像素值采用公式9的加权平均方法进行图像融合,保留更多的图像细节信息:
采用上述方法对左图参照图像I1'和右图参照图像I2'的拼接融合结果如图6(a)所示,原始左图图像I1和原始右图图像I2的拼接融合结果如图6(b)所示。
由此,本实施例最终采用了加权平均融合方法生成无缝拼接图像,实现了对于宽幅织物图像的获取与拼接。
上述具体实施方式用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。

Claims (7)

1.一种基于参照特征的宽幅织物图像获取与拼接方法,其特征在于该方法的步骤如下:
1)先将被拼接织物水平放在织物载台上,获取占有被拼接织物左侧三分之二的图像,将该图像进行归一化处理得到被拼接织物的原始左图图像I1
2)在被拼接织物的中间水平放上参照织物,获取置有参照织物的左图参照图像I1',左图参照图像I1'与原始左图图像I1大小相同;水平向左移动织物载台,获取置有参照织物的右图参照图像I2',并使得左图参照图像I1'和右图参照图像I2'共同有至少参照织物三分之一的重叠区域;
3)再保持被拼接织物位置不变,移走参照织物,获取被拼接织物的原始右图图像I2,原始右图图像I2与右图参照图像I2'大小相同;
4)对于左图参照图像I1'和右图参照图像I2',均采用多尺度Harris角点检测方法分别得到左图参照图像I1'的左图Harris角点集S1和右图参照图像I2'的右图Harris角点集S2
5)采用归一化的相关测度方法对两个角点集S1和S2进行Harris角点的初始匹配,得到初始匹配点集S;
6)采用RANSAC方法对初始匹配点集S进行初始匹配特征点求精和单应矩阵的求解,得到初始匹配特征点求精结果集合Si’,并得到最终单应矩阵H';
7)根据得到的最终单应矩阵H',对原始左图图像I1和原始右图图像I2进行投影变换配准,再拼接并利用加权平均方法消除拼接缝,进行图像的融合,完成织物图像的获取与拼接。
2.根据权利要求1所述的一种基于参照特征的宽幅织物图像获取与拼接方法,其特征在于:所述的步骤4)中左图参照图像I1'和右图参照图像I2'的两个参照图像均采用以下多尺度Harris角点检测方法进行处理:
4.1)图像的尺度空间按高斯窗的尺度δm采用以下公式1进行计算,尺度δm在0.5~5的范围内以0.5的步长变化并依次计算:
L x ( &delta; m ) = G ( x , y , &delta; m ) &CircleTimes; ( I x &prime; ) 2 L y ( &delta; m ) = G ( x , y , &delta; m ) &CircleTimes; ( I y &prime; ) 2 L x y ( &delta; m ) = G ( x , y , &delta; m ) &CircleTimes; ( I x &prime; I y &prime; ) 2 - - - ( 1 )
式中,Ix'、Iy'分别为参照图像在x方向和y方向的偏导,G(x,y,δm)为高斯核,是尺度变换的唯一变换核,m为尺度计算的序数,Lxm)、Lym)、Lxym)分别表示x方向、y方向和xy方向的高斯尺度空间;
4.2)由上述步骤4.1)的公式1计算得到不同尺度δm上的矩阵M:
M = L x ( &delta; m ) L x y ( &delta; m ) L x y ( &delta; m ) L y ( &delta; m ) - - - ( 2 )
以及采用以下公式3计算得到在不同尺度δm上的响应函数C(x,y):
C(x,y)=Det[M(x,y)]-βTrace2[M(x,y)] (3)
式中,Det[M(x,y)]为矩阵M的行列式,Trace[M(x,y)]为矩阵M的迹,β是常数因子;
4.3)如果响应函数C(x,y)大于响应函数阈值T1,则响应函数C(x,y)所对应的像素点(x,y)为角点;’
4.4)对于尺度δm以0.5的步长在0.5~5的范围内变化的每个尺度,采用非极大抑制法抑制多响应的影响,非极大抑制窗的尺寸为3×3,将响应函数C(x,y)的值按照从大到小降序排列,将各个尺度下的左图参照图像I1'和右图参照图像I2'中可能存在的角点作为候选角点;
4.5)对于尺度δm以0.5的步长在0.5~5的范围内变化的第二个尺度δm=1开始到最大尺度δm=5,依次遍历各个尺度,查找是否有当前尺度下的候选角点在前一尺度中出现,如果没有,保留该候选角点;否则,剔除候选角点;
由此得到左图参照图像I1'的左图Harris角点集S1和右图参照图像I2'的右图Harris角点集S2
3.根据权利要求1所述的一种基于参照特征的宽幅织物图像获取与拼接方法,其特征在于:所述的步骤5)的左图参照图像I1'和右图参照图像I2'的两个参照图像具体采用以下方法进行初始匹配:
5.1)在左图参照图像I1'中,以左图Harris角点集S1的任一Harris角点X1(x1,y1)为中心,使得能建立大小为3×3的搜索窗口;
5.2)对于左图参照图像I1'和右图参照图像I2',均采用以下公式4依次遍历各个Harris角点分别计算两个参照图像灰度值的标准差:
s t d d ( X k ) = &Sigma; i = - N i = N &Sigma; j = - N j = N &lsqb; I k &prime; ( x + i , y + j ) - I &OverBar; k &prime; ( X k ) &rsqb; 2 ( 2 N + 1 ) &times; ( 2 N + 1 ) - - - ( 4 )
式中:为参照图像中作为Harris角点的像素点Xk(x,y)处的灰度平均值,Ik'(x+i,y+j)表示参照图像中像素点Xk(x+i,y+j)处的灰度值,k为左图Harris角点集S1和右图Harris角点集S2中的Harris角点序数,Xk为左图Harris角点集S1和右图Harris角点集S2中的Harris角点,i、j分别表示参照图像的x方向、y方向的像素点序数,N为参照图像的像素点总数;
5.3)分别遍历两个参照图像中的各个像素点,采用以下公式5计算两幅参照图像之间的协方差:
cov ( X 1 , X 2 ) = &Sigma; i = - N i = N &Sigma; j = - N j = N &lsqb; I 1 &prime; ( x 1 + i , y 1 + j ) - I &OverBar; 1 &prime; ( X 1 ) &rsqb; &CenterDot; &lsqb; I 2 &prime; ( x 2 + i , y 2 + j ) - I &OverBar; 2 &prime; ( X 2 ) &rsqb; ( 2 N + 1 ) &CenterDot; ( 2 N + 1 ) &CenterDot; ( 2 N + 1 ) &CenterDot; ( 2 N + 1 ) - - - ( 5 )
其中,X1和X2分别为左图Harris角点集S1和右图Harris角点集S2中的Harris角点;
5.4)采用以下公式6的归一化自相关函数cor(X1,X2)进行归一化自相关计算:
c o r ( X 1 , X 2 ) = cov ( X 1 , X 2 ) s t d d ( X 1 ) &CenterDot; s t d d ( X 2 ) - - - ( 6 )
5.5)当归一化自相关函数cor(X1,X2)大于自相关阈值T2时,则该归一化自相关函数cor(X1,X2)对应的Harris角点为候选匹配点,由所有候选匹配点组成初始匹配特征点集S。
4.根据权利要求1所述的一种基于参照特征的宽幅织物图像获取与拼接方法,其特征在于:所述的步骤6)具体过程如下:
6.1)将初始匹配特征点集S里面的匹配点对作为样本,根据置信概率α和所有样本的数量计算得到迭代次数t,采用以下步骤6.2)~6.4)进行循环迭代计算;
6.2)从初始特征匹配点集S中随机选取l对匹配点对,利用l对匹配点对采用以下公式7计算得到两幅参照图像之间的临时单应矩阵Hv
X v Y v 1 = H v x v y v 1 - - - ( 7 )
式中,(Xv,Yv)是左图参照图像I1'中的Harris角点,(xv,yv)是右图参照图像I2'中的Harris角点;v为当前循环次数的序数,第一次计算的序数v为1;
6.3)由临时单应矩阵Hv采用最小二乘法,采用以下公式8计算初始特征匹配点集S剩余各个匹配点对的极线距离dv
d v = X v Y v 1 H v x v y v 1 - - - ( 8 )
若极线距离dv小于极线距离阈值T3,则将该极线距离dv所对应的匹配点对加入第v匹配点集Sv,该极线距离dv加入第v距离集合Dv
6.4)将当前循环次数的序数v加1,若加1后的循环次数的序数<t,t为迭代次数,则跳转到步骤6.2),否则进行下一步骤6.5);
6.5)选取所有匹配点集中匹配点对最多的匹配点集作为初始匹配特征点求精结果集合Sv’,所对应的临时单应矩阵Hv作为最终单应矩阵H'。
5.根据权利要求1所述的一种基于参照特征的宽幅织物图像获取与拼接方法,其特征在于:所述的步骤7)具体过程如下:建立纯黑色的背景图像I12,背景图像I12各个像素点的RGB像素值均为(255,255,255),计算最终单应矩阵H'和原始右图图像I2的乘积,将原始右图图像I2投影变换到背景图像I12中,背景图像I12中重叠区域外对应像素值按左图和右图的像素值不变,重叠区域处的像素值采用公式9的加权平均方法进行图像融合,由此获得完整的拼接织物图像:
I 12 ( x , y ) = I 1 ( x , y ) d 1 + I 2 ( x , y ) d 2 d 1 + d 2 - - - ( 9 )
式中,I1(x,y)为原始左图图像I1重叠区域的像素值;I2(x,y)为原始右图图像I2重叠区域的像素值;d1为像素点(x,y)距原始右图图像左边缘的垂直距离;d2为像素点(x,y)距原始左图图像右边缘的垂直距离;I12(x,y)表示融合后的像素点(x,y)的像素值。
6.根据权利要求4所述的一种基于参照特征的宽幅织物图像获取与拼接方法,其特征在于:所述的步骤6.1)中迭代次数t采用以下公式10计算得到:
t=wα (10)
式中,w为初始匹配特征点集S中特征点数量,α为置信概率。
7.根据权利要求1所述的一种基于参照特征的宽幅织物图像获取与拼接方法,其特征在于:所述方法用于对宽幅织物图像的获取与拼接。
CN201410854510.5A 2014-12-31 2014-12-31 一种基于参照特征的宽幅织物图像获取与拼接方法 Active CN104637041B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410854510.5A CN104637041B (zh) 2014-12-31 2014-12-31 一种基于参照特征的宽幅织物图像获取与拼接方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410854510.5A CN104637041B (zh) 2014-12-31 2014-12-31 一种基于参照特征的宽幅织物图像获取与拼接方法

Publications (2)

Publication Number Publication Date
CN104637041A CN104637041A (zh) 2015-05-20
CN104637041B true CN104637041B (zh) 2017-06-23

Family

ID=53215747

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410854510.5A Active CN104637041B (zh) 2014-12-31 2014-12-31 一种基于参照特征的宽幅织物图像获取与拼接方法

Country Status (1)

Country Link
CN (1) CN104637041B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106296580A (zh) * 2015-06-03 2017-01-04 无锡天脉聚源传媒科技有限公司 一种图像拼接的方法及装置
CN105447864B (zh) * 2015-11-20 2018-07-27 小米科技有限责任公司 图像的处理方法、装置及终端
CN105806485B (zh) * 2016-03-21 2017-09-26 浙江大学 一种用于纺织品自动色彩采集及校准的方法
CN106651903B (zh) * 2016-12-30 2019-08-09 明见(厦门)技术有限公司 一种运动物体检测方法
CN109099835A (zh) * 2018-07-02 2018-12-28 广西电网有限责任公司北海供电局 基于手机影像的电力线夹测量方法及其系统
CN109345513B (zh) * 2018-09-13 2021-06-01 红云红河烟草(集团)有限责任公司 一种具有烟包姿态计算的卷烟包装缺陷检测方法
CN109598675B (zh) * 2018-11-13 2023-03-10 北京交通大学 多重复纹理图像的拼接方法
CN111755084B (zh) * 2020-06-22 2021-03-16 常州市新创智能科技有限公司 一种织物铺层检验方法
CN113609943B (zh) * 2021-07-27 2024-05-17 东风汽车有限公司东风日产乘用车公司 手指指静脉识别方法、电子设备及存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101866482A (zh) * 2010-06-21 2010-10-20 清华大学 基于摄像设备自标定技术的全景图拼接方法和装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9182228B2 (en) * 2006-02-13 2015-11-10 Sony Corporation Multi-lens array system and method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101866482A (zh) * 2010-06-21 2010-10-20 清华大学 基于摄像设备自标定技术的全景图拼接方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于ASIFT的无缝图像拼接方法;杨小辉等;《计算机工程》;20130228;第39卷(第2期);241-244,249 *

Also Published As

Publication number Publication date
CN104637041A (zh) 2015-05-20

Similar Documents

Publication Publication Date Title
CN104637041B (zh) 一种基于参照特征的宽幅织物图像获取与拼接方法
CN106485275B (zh) 一种实现盖板玻璃与液晶屏定位贴合的方法
CN108717531B (zh) 基于Faster R-CNN的人体姿态估计方法
US11051000B2 (en) Method for calibrating cameras with non-overlapping views
CN109961399B (zh) 基于图像距离变换的最佳缝合线搜寻方法
Tuytelaars et al. Noncombinatorial detection of regular repetitions under perspective skew
CN105096307B (zh) 在成对的立体图像中检测物体的方法
CN103679702B (zh) 一种基于图像边缘矢量的匹配方法
CN102834845A (zh) 用于多摄像装置校准的方法与装置
CN107424160A (zh) 通过视觉系统查找图像中线的系统和方法
Surh et al. Noise robust depth from focus using a ring difference filter
CN104008542B (zh) 一种针对特定平面图形的快速角点匹配方法
CN108369650A (zh) 将校准图案的图像内的候选点标识为该校准图案的可能特征点的方法
CN103136525B (zh) 一种利用广义Hough变换的异型扩展目标高精度定位方法
CN107993258A (zh) 一种图像配准方法及装置
CN105989603A (zh) 机器视觉图像传感器校准
CN103426186A (zh) 一种改进的surf快速匹配方法
CN109598234A (zh) 关键点检测方法和装置
CN104599258A (zh) 一种基于各向异性特征描述符的图像拼接方法
CN110232387A (zh) 一种基于kaze-hog算法的异源图像匹配方法
CN108052909B (zh) 一种基于心血管oct影像的薄纤维帽斑块自动检测方法和装置
CN102930251A (zh) 两维收藏品数据收录与甄别的装置和方法
CN110288511A (zh) 基于双相机图像的最小化误差拼接方法、装置、电子设备
CN109815823A (zh) 数据处理方法及相关产品
CN106296587A (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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20200310

Address after: 310012 Room 309, floor 3, building 5, 198 Jiaogong Road, Xixi street, Xihu District, Hangzhou City, Zhejiang Province

Patentee after: Hangzhou Guangyu New Material Technology Co., Ltd

Address before: 310018, No. 928, No. 2, Poplar Street, Jianggan Economic Development Zone, Hangzhou, Zhejiang

Patentee before: ZHEJIANG SCI-TECH University

TR01 Transfer of patent right