CN108917632B - 一种高效率高精度数字图像相关位移后处理方法 - Google Patents
一种高效率高精度数字图像相关位移后处理方法 Download PDFInfo
- Publication number
- CN108917632B CN108917632B CN201810461488.6A CN201810461488A CN108917632B CN 108917632 B CN108917632 B CN 108917632B CN 201810461488 A CN201810461488 A CN 201810461488A CN 108917632 B CN108917632 B CN 108917632B
- Authority
- CN
- China
- Prior art keywords
- displacement
- template
- digital image
- dic
- post
- 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.)
- Expired - Fee Related
Links
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000012805 post-processing Methods 0.000 title claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 30
- 238000012937 correction Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 abstract description 14
- 238000000691 measurement method Methods 0.000 abstract description 2
- 230000003287 optical effect Effects 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000010297 mechanical methods and process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000009666 routine test Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/24—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种高效率高精度数字图像相关位移后处理方法,包括计算各计算模板覆盖区域内的高斯点坐标、计算各计算模板覆盖区域内高斯点位移、高斯点位移修正步骤。本发明通过对传统DIC获得的位移及位移梯度进行后处理,能够在付出极小计算代价的条件下将位移精度提高一阶,是一种兼具高精度和高效率的测量方法;本发明将插值后高斯点位移的平均值作为最终位移测量结果可以有效抑制随机误差,具有较强的鲁棒性;本发明的测量方法仅需对传统DIC获得的位移及其梯度进行简单数学运算,不需要改动DIC核心算法,编程复杂度低,易于实现。
Description
技术领域
本发明涉及一种数字图像相关位移后处理方法,尤其涉及一种高效率高精度数字图像相关位移后处理方法,属于光测力学技术领域。
背景技术
数字图像相关法(Digital image correlation,DIC)是一种基于现代图像处理技术、通过匹配变形前后图像中具有强辨识度散斑的灰度值进行计量的光测力学方法。因其具有对测试环境要求宽松、能够获得全场位移和应变等优点,目前DIC已作为理想的变形测量手段应用于多个领域。
基于国内外学者在DIC亚像素算法及其影响因素方面的研究成果,目前DIC算法在均匀变形测量上已能够满足常规试验的精度要求。但在非均质材料的非均匀变形测量上,测量偏差较大,无法满足非均质材料高精度的全场测量要求。造成测量偏差增大的主要原因是基于常应变假设的一阶计算模板形函数与真实位移场的欠匹配。为提高匹配性,现有的解决方案是减小计算模板尺寸或提高计算模板形函数的阶次,但随着计算模板尺寸的减小,位移随机误差急剧增大甚至无法收敛,而提高形函数阶次又将显著增加DIC的计算复杂度和对图像噪声的敏感性,无法满足高效率和高鲁棒性的测量需求。
因此,寻求一种能同时满足精度、效率和鲁棒性的方法是目前应用DIC进行非均匀变形测量的关键技术问题。
发明内容
本发明要解决的技术问题是提供一种高效率高精度数字图像相关位移后处理方法。
为解决上述技术问题,本发明采用的技术方案是:
一种高效率高精度数字图像相关位移后处理方法,在待测区域划分网格点,设定网格步长,网格点坐标为(xi,yi),i=1,2,...,n,n为网格点数量,以每一个网格点为中心取其邻域作为计算模板,模板尺寸为M×N,网格点位移(ui,vi)和位移梯度由传统DIC方法获得;包括以下具体步骤:
进一步,还包括高斯点位移修正步骤:
采用上述技术方案所取得的技术效果在于:
1)本发明通过对DIC算法获得的位移及位移梯度进行后处理,能够在付出极小计算代价的条件下将位移精度提高一阶,是一种兼具高精度和高效率的测量方法;
2)本发明将插值后高斯点位移的平均值作为最终位移测量结果可以有效抑制随机误差,具有较强的鲁棒性;
3)本发明的测量方法仅需对DIC算法获得的位移及其梯度进行简单数学运算,不需要改动DIC核心算法,编程复杂度低,易于实现。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1是本发明的流程图;
图2是本发明的实施例1中高斯点位置示意图;
图3是本发明的实施例1中测量结果与DIC算法测量结果的对比图。
具体实施方式
实施例1:
参考图1-图3,一种高效率高精度数字图像相关位移后处理方法,在待测区域划分网格点,设定网格步长,网格点坐标为(xi,yi),i=1,2,...,n,n为网格点数量,以每一个网格点为中心取其邻域作为计算模板,模板尺寸为M×N,网格点位移(ui,vi)和位移梯度由传统DIC方法获得;包括以下具体步骤:
进一步,还包括高斯点位移修正步骤:
在本实施例中,采用Zhou Peng等,见Zhou P,Goodson K E.Subpixeldisplacement and deformation gradient measurement using digital image/specklecorrelation[J].Optical Engineering,2001,40(8):1613-1621.提出的模拟散斑算法生成变形前后散斑图。其中,图像尺寸为200×400像素,散斑数量为3200个,散斑大小为5像素。变形场为u(x,y)=sin(2πx/100),v(x,y)=0。
以变形前散斑图作为参考图像,在待测区域划分网格点,本实施例中网格步长取3×15像素,网格点数量n=536,网格点坐标为(xi,yi),i=1,2,...,n。以每一个网格点为中心取其邻域作为计算模板,计算模板尺寸41×41像素,计算模板中像素点在变形图像中的像素坐标与参考图像中的像素坐标之间的映射关系采用下式表示的一阶形式:
进一步地,采用本发明的方法对网格点,即计算模板中心点位移及其梯度做进一步后处理,获得具有二阶精度的高斯点位移结果。
本发明的核心内容是对DIC算法获得的位移及位移梯度数据进行后处理,通过将计算模板中心点位移外推至高斯点位移来提高位移精度,因此,本发明要求保护的范围并不局限于一阶形函数,对于计算模板上的二阶或更高阶形函数,可以通过增加高斯点个数,采用本发明的步骤获得更高精度的位移结果。
参考图3,从图中可以看出,本发明的数字图像相关位移后处理方法能够显著提高位移测量精度。此外,由于本发明的测量方法仅需对DIC算法获得的位移及其梯度进行简单数学运算,计算复杂度低,具有较高计算效率。
Claims (2)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810461488.6A CN108917632B (zh) | 2018-05-15 | 2018-05-15 | 一种高效率高精度数字图像相关位移后处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810461488.6A CN108917632B (zh) | 2018-05-15 | 2018-05-15 | 一种高效率高精度数字图像相关位移后处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108917632A CN108917632A (zh) | 2018-11-30 |
CN108917632B true CN108917632B (zh) | 2020-06-02 |
Family
ID=64403853
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810461488.6A Expired - Fee Related CN108917632B (zh) | 2018-05-15 | 2018-05-15 | 一种高效率高精度数字图像相关位移后处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108917632B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111986152B (zh) * | 2020-07-17 | 2024-03-15 | 南京航空航天大学 | 一种数字图像相关变形测量数据处理时长预测方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4836067B2 (ja) * | 2005-05-23 | 2011-12-14 | 日立造船株式会社 | 構造物における変形量計測方法 |
US8107694B2 (en) * | 2006-07-20 | 2012-01-31 | Ultrasound Medical Devices, Inc. | Method of tracking speckle displacement between two images |
CN102221341B (zh) * | 2011-03-16 | 2013-01-23 | 中国人民解放军国防科学技术大学 | 基于随机并行梯度下降优化技术的快速数字图像相关测量方法 |
CN102967261B (zh) * | 2012-11-10 | 2015-04-15 | 华中科技大学 | 一种基于数字散斑相关法的激光位移测量方法 |
CN103729849A (zh) * | 2013-12-31 | 2014-04-16 | 南京航空航天大学 | 一种数字图像变形初值的计算方法 |
CN103761768A (zh) * | 2014-01-22 | 2014-04-30 | 杭州匡伦科技有限公司 | 一种三维重建的立体匹配方法 |
CN104197846B (zh) * | 2014-09-18 | 2017-03-15 | 中国科学院光电技术研究所 | 一种基于特征图样反射的反射镜位移测量方法 |
CN105139346B (zh) * | 2015-07-09 | 2019-05-14 | Tcl集团股份有限公司 | 一种数字图像处理方法和数字图像处理装置 |
CN106202801B (zh) * | 2016-07-22 | 2019-07-09 | 中国人民解放军空军预警学院 | 一种天基光学空间碎片监视图像中目标运动拖尾仿真方法 |
CN107945216B (zh) * | 2017-11-10 | 2019-10-11 | 西安电子科技大学 | 基于最小二乘估计的多图像联合配准方法 |
-
2018
- 2018-05-15 CN CN201810461488.6A patent/CN108917632B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN108917632A (zh) | 2018-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103292701B (zh) | 基于机器视觉的精密器件在线尺寸测量方法 | |
CN111288967B (zh) | 一种基于机器视觉的远距离高精度位移检测方法 | |
CN109360246B (zh) | 基于同步子区搜索的立体视觉三维位移测量方法 | |
CN107798698B (zh) | 基于灰度修正与自适应阈值的结构光条纹中心提取方法 | |
CN107328502B (zh) | 一种锚杆托盘载荷可视化数字成像方法 | |
CN110455222B (zh) | 一种高精度旋转角测量方法、装置及设备 | |
CN110906875B (zh) | 一种孔径测量的视觉处理方法 | |
CN109373912B (zh) | 一种基于双目视觉的非接触式六自由度位移测量方法 | |
CN103994732A (zh) | 一种基于条纹投影的三维测量方法 | |
CN103632338A (zh) | 一种基于匹配曲线特征的图像配准评估法 | |
CN104657999B (zh) | 一种基于核函数的数字图像相关方法 | |
CN111256607B (zh) | 一种基于三通道标志点的变形测量方法 | |
CN108917632B (zh) | 一种高效率高精度数字图像相关位移后处理方法 | |
CN113012143B (zh) | 一种基于二维数字图像相关法的试件质量检测方法 | |
CN107240077B (zh) | 一种基于椭圆构像偏差迭代矫正的视觉测量方法 | |
CN114820820B (zh) | 一种电脑后壳的3d检测方法 | |
CN113838120B (zh) | 一种基于数字图像相关的大转角运动位移测量方法 | |
CN107328365B (zh) | 一种柔顺定位平台位移的鲁棒测量方法 | |
CN113532308B (zh) | 数字图像相关中带初值的岭回归应变测量方法 | |
CN115512343A (zh) | 一种圆形指针式仪表校正和读数识别方法 | |
CN112927299B (zh) | 标定方法、装置及电子设备 | |
CN114705216B (zh) | 一种用于三维视觉测量系统的二次标定方法 | |
CN116500851B (zh) | 3d成像系统及其标定方法、装置、设备和存储介质 | |
CN110060305B (zh) | 一种高精度简化式线阵相机标定方法 | |
CN114299079B (zh) | 一种面向稠密点云数据的发动机叶片截面线数据获取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200602 |