CN106023097B - 一种基于迭代法的流场图像预处理方法 - Google Patents

一种基于迭代法的流场图像预处理方法 Download PDF

Info

Publication number
CN106023097B
CN106023097B CN201610307839.9A CN201610307839A CN106023097B CN 106023097 B CN106023097 B CN 106023097B CN 201610307839 A CN201610307839 A CN 201610307839A CN 106023097 B CN106023097 B CN 106023097B
Authority
CN
China
Prior art keywords
image
information
flow field
scale structure
particle
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
CN201610307839.9A
Other languages
English (en)
Other versions
CN106023097A (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.)
Huazhong University of Science and Technology
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201610307839.9A priority Critical patent/CN106023097B/zh
Publication of CN106023097A publication Critical patent/CN106023097A/zh
Application granted granted Critical
Publication of CN106023097B publication Critical patent/CN106023097B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/77—Retouching; Inpainting; Scratch removal
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20172—Image enhancement details
    • G06T2207/20192—Edge enhancement; Edge preservation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种基于迭代法的流场图像预处理方法。该方法可以根据尺度信息，利用高斯滤波特点将微小粒子信息完全滤除，然后对处理后的图像利用可迭代的边缘恢复算法，将模糊化的大尺度结构信息不断恢复，恢复后的大尺度结构信息就是流场图像的噪声信息。最后，利用图像背景相减法即可得到实际粒子图像信息。按照本发明实现的流场图像预处理方法可以根据粒子尺寸信息设置距离权重来实现微小粒子去除，并且采用基于迭代方法的边缘恢复获得背景噪声图像具有较好的收敛特性，采用减去噪声背景的方法可以很好地消除由于光照、粒子不均等引起的噪声，因此最终获得的粒子图像具有很好的边缘特性且对比度较高。

Description

一种基于迭代法的流场图像预处理方法
技术领域
本发明属于图像处理技术领域，具体地，涉及一种基于迭代法的流场图像预处理方法。
背景技术
流场图像在进行实际粒子图像测速算法之前，由于图像本身的一些噪声，都要进行图像的预处理操作。预处理的基本目的是通过图像滤波和增强技术以减小噪声、激光反射、粒子聚团等所造成的速度误测，其根本目的就是获得更强对比度的粒子图像。图像预处理可以提高粒子图像的质量，减小后期算法的消耗，较好的预处理算法可以很好地提高后期流场处理的精度。
流场中粒子图像具有以下特点，粒子尺寸相对均匀，分布散乱无具体分布规律，实际背景噪声更不可避免，为提高图像处理的质量，对预处理算法速度要求较高。
在现有技术中，一般通过高斯滤波的方式进行图像处理，在二维高斯滤波中，目标点上的像素值是通过如下方式获得的：对所在位置周围一定范围内的像素值分别赋予不同的高斯权重值，并在加权平均后得到当前的结果，但此种方法仅考虑到了位置权重信息，却没有考虑到像素间的相似程度信息，因此单纯使用高斯滤波模板通常会使整张图片变得模糊，应用在流场预处理算法的应用中具有一定局限性。由此，亟待涉及一种针对流场图像的预处理方法，来提高图像处理质量。
发明内容
针对现有技术的不足，本发明提出了基于迭代的、距离权重和像素相似信息权重、背景相减法的流场图像预处理方法，旨在去除流场图像中的噪声，最终达到增强粒子数据对比度的目的。
为实现以上发明目的，本发明采用以下技术方案：
本发明提出一种基于迭代法的流场图像预处理方法，其特征在于，该方法包括如下步骤：
(1)采集流场原始图像函数I；
(2)利用高斯滤波对所述原始图像函数I进行处理进而去除小尺度结构信息获得输出图像函数O，其中所述小尺度结构信息为低于4像素的图像信息；
(3)利用双边滤波函数迭代获取大尺度结构信息恢复函数Gn+1，其中n为迭代次数，所述迭代次数满足设定的阈值精度则停止迭代，其中所述大尺度结构信息为高于15像素的图像信息；
(4)获得最终图像信息为J＝I-Gn+1。
进一步地，所述步骤(2)的滤波过程为：
其中，用来实现归一化，N(α)为α点邻域的像素集合。
进一步地，所述步骤(3)的迭代过程为：
其中，用来对结果归一化，N(α)为α点邻域的像素集合，其中σs和σr分别为基于空间距离和像素间相似程度的权重参数。
进一步地，在所述步骤(2)中，其中可通过调节σs，消除不同尺寸的小尺度结构信息。
本发明具有以下有益效果：
1.利用高斯滤波算法对小尺度轮廓信息进行消除，因为流场中示踪粒子的尺寸大小较为均匀，因此具有很好的滤除效果；
2.基于迭代的边缘恢复，对小尺度轮廓信息完全无影响，且可以通过原始图像信息对模糊化的大尺度轮廓信息恢复；
3.基于迭代的流场图像预处理方法具有很好的收敛性，且收敛于输入图像信息。
附图说明
图1是按照本发明实现的方法流程示意图；
图2是本发明实现的方法具体实施过程效果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
实际上在现有的图像处理领域中，双边滤波在实际处理时不仅考虑了图像像素在空间距离上的信息，同时考虑到了像素间的相似程度信息，因而可以保持原始图像的大体分块信息，进而保持边缘信息。
基于背景相减原理的预处理算法在图像的预处理中具有较为广泛的应用，对于去除背景噪声具有很好的实际应用效果。
本发明主要是在以上几种算法的研究基础上，提出了一种基于迭代的流场图像预处理方法。可很好在保持粒子轮廓信息的基础上，达到较好效果地滤除流场图像中非粒子图像的大尺度结构信息。
在本实施例中，图像数据通过行列方向的二维数组来表示，每个数组元素的值代表该点的像素灰度值，以图像左上角为图像坐标的原点，垂直向下方向为Y正方向。
总体而言，本发明所描述的流场图像预处理方法主要分为如下步骤：
(1)去除小尺度结构信息；
(2)大尺度结构信息迭代恢复；
(3)图像背景相减。
本发明所提出的流场图像预处理方法中，对于高斯滤波的尺度一般是根据标准差来决定的，实际的流场图像处理中实际粒子大小在8个像素大小左右，所以在本算法涉及的步骤中，进行图像处理时，将高于15像素的图像信息称为大尺度结构信息，低于4像素的为小尺度结构信息。
下面对本发明提出的流场图像预处理方法作详细说明，具体步骤如下：
(1)去除小尺度结构信息
第一步先通过高斯滤波滤除小尺度结构信息，尺度算子通过加权平均的形式表述，在公式中采用I为输入图像，O为输出图像，α和β分别表示图像坐标系中的像素坐标，σs表示标准差，决定高斯滤波器的宽度，滤波器函数表达式为：
其中，用来对结果归一化，N(α)为α点邻域的像素集合。
使用这个滤波器可以将小的边缘信息全部都滤除，而大轮廓边缘信息被模糊化，注意，此步骤中通过对σs的调节，即通过不同尺寸粒子的消除设置不同的标准差来完成调节，可以很好地适用于不同尺寸的小尺寸结构信息粒子的消除，将尺度小于σs的粒子信息全部滤除，使该预处理算法在粒子尺寸不同时中也有很好的适用性。
(2)大尺度结构信息恢复
基于迭代方法的大轮廓边缘恢复是该方法的关键之一。在该过程中，图像G被不断迭代更新，用Gn+1表示第n+1次的迭代输出图像，G1就是公式(1)中的O，即公式(1)中高斯滤波处理后的输出图像数据。第n次的迭代结果Gn+1的是通过双边滤波函数获得的，该函数的输入由最初原始图像I和上一次的迭代输出结果Gn，该双边滤波函数公式为：
其中，用来对结果归一化，N(α)为α点邻域的像素集合，其中σs和σr分别为基于空间距离和像素间相似程度的权重参数。
利用上述式子进行迭代处理，直到达到所需精度终止算法迭代，例如流场测量互相关算法中，我们算法迭代至0.1pixel精度终止运行，从而获得最终图像函数Gn+1以供下一步处理。
同时，在公式(2)中，由于小尺度结构信息已经被公式(1)完全去除，因此||Gn(α)-Gn(β)||约等于0，这也导致基于像素间相似程度的权重对小的轮廓信息完全不起作用。也就是说，之后的所有的迭代操作所具有的边缘恢复对小的轮廓信息都不起作用，小的轮廓信息将无法被恢复。
在大尺度结构信息迭代恢复的过程中，模糊化的大尺度结构信息被逐渐恢复，最终输出的结果为恢复具有大尺度结构信息的图像，小尺度结构信息完全被滤除。因此，获得的图像数据为去除粒子信息的噪声数据。
(3)背景相减，获得粒子图像
在PIV图像中，由于示踪粒子尺寸相对较小且尺寸相似，因此利用公式(1)的高斯滤波将示踪粒子轮廓从图像中滤除。此时，获得的图像就是去除了所有小的轮廓信息(粒子图像)和部分大轮廓信息的数据图像，将模糊后的一些噪声信息重新进行边缘恢复，最终公式(2)中所获得的图像数据就是滤除所有示踪粒子的噪声信息。下面，就需要利用原始图像数据减去边缘恢复后图像噪声数据，即可获得示踪粒子轮廓信息，具体如下公式：
J＝I-Gn+1 (3)
J即为最终获得的流场图像中的示踪粒子图像信息。
如图2所示，是按照本发明实现的流场迭代图像处理方法中的图像处理的具体执行之后的图像效果图，可以看出在经过了上述的高斯滤波和迭代之后获得的图像能够获得较好的效果，并且能够实现。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (2)

1.一种基于迭代法的流场图像预处理方法，其特征在于，该方法包括如下步骤：
(1)采集流场原始图像函数I；
(2)利用高斯滤波对原始图像函数I进行处理进而去除小尺度结构信息获得输出图像数据O，其中所述小尺度结构信息为低于4像素的图像信息；
(3)利用双边滤波函数迭代获取大尺度结构信息恢复函数Gn+1，其中n为迭代次数，G1是步骤(2)中的输出图像数据O，所述迭代次数满足设定的阈值精度则停止迭代，其中所述大尺度结构信息为高于15像素的图像信息；
(4)获得最终图像信息为J＝I-Gn+1；
所述步骤(2)的滤波过程为：
其中，用来实现归一化，α和β为图像坐标系中的像素坐标，N(α)为α点邻域的像素集合；
所述步骤(3)的迭代过程为：
其中，用来对结果归一化，N(α)为α点邻域的像素集合，其中σs和σr分别为基于空间距离和像素间相似程度的权重参数。
2.如权利要求1所述的基于迭代法的流场图像预处理方法，其特征在于，在所述步骤(2)中，其中通过调节σs，消除不同尺寸的小尺度结构信息。
CN201610307839.9A 2016-05-11 2016-05-11 一种基于迭代法的流场图像预处理方法 Active CN106023097B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610307839.9A CN106023097B (zh) 2016-05-11 2016-05-11 一种基于迭代法的流场图像预处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610307839.9A CN106023097B (zh) 2016-05-11 2016-05-11 一种基于迭代法的流场图像预处理方法

Publications (2)

Publication Number Publication Date
CN106023097A CN106023097A (zh) 2016-10-12
CN106023097B true CN106023097B (zh) 2019-04-12

Family

ID=57100306

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610307839.9A Active CN106023097B (zh) 2016-05-11 2016-05-11 一种基于迭代法的流场图像预处理方法

Country Status (1)

Country Link
CN (1) CN106023097B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106908218B (zh) * 2017-03-23 2018-11-30 北京尚水信息技术股份有限公司 河道表面流场的处理方法
CN107977935A (zh) * 2017-11-15 2018-05-01 合肥中科奔巴科技有限公司 一种掌静脉图像纹理滤波增强方法
CN109919853B (zh) * 2019-01-07 2023-04-25 南京航空航天大学 一种基于粒子特征检索的piv图像预处理方法
CN111398625B (zh) * 2020-03-19 2022-04-12 西安理工大学 一种物理模型试验中的测速方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7146059B1 (en) * 2003-03-05 2006-12-05 Massachusetts Institute Of Technology Method of performing fast bilateral filtering and using the same for the display of high-dynamic-range images
CN102789636A (zh) * 2012-08-01 2012-11-21 中山大学 基于多尺度联合双边栅格平滑模型的图像细节增强方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7146059B1 (en) * 2003-03-05 2006-12-05 Massachusetts Institute Of Technology Method of performing fast bilateral filtering and using the same for the display of high-dynamic-range images
CN102789636A (zh) * 2012-08-01 2012-11-21 中山大学 基于多尺度联合双边栅格平滑模型的图像细节增强方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Horn–Schunck optical flow applied to deformation measurement of a birdlike airfoil;Gong Xiaoliang等;《Chinese Journal of Aeronautics》;20151015;第28卷(第5期);第1305-1315页
基于变分原理的光流估计方法研究;谭君;《中国优秀硕士学位论文全文数据库信息科技辑》;20160315(第03期);摘要，第23-32页

Also Published As

Publication number Publication date
CN106023097A (zh) 2016-10-12

Similar Documents

Publication Publication Date Title
Guo et al. LIME: Low-light image enhancement via illumination map estimation
Wang et al. Dehazing for images with large sky region
CN105608676B (zh) 一种视频图像的增强方法及装置
CN108230264B (zh) 一种基于ResNet神经网络的单幅图像去雾方法
CN106023097B (zh) 一种基于迭代法的流场图像预处理方法
CN108564549B (zh) 一种基于多尺度稠密连接网络的图像去雾方法
Matakos et al. Accelerated edge-preserving image restoration without boundary artifacts
Yin et al. Highly accurate image reconstruction for multimodal noise suppression using semisupervised learning on big data
CN105844601A (zh) 一种基于双边滤波和多尺度Retinex算法的矿井图像增强方法
Zhao et al. Single image fog removal based on local extrema
Yan et al. Method to Enhance Degraded Image in Dust Environment.
CN102819827A (zh) 一种基于灰度分割的自适应矩匹配条带噪声去除方法
CN109949256B (zh) 一种基于傅里叶变换的天文图像融合方法
Yu et al. Image and video dehazing using view-based cluster segmentation
CN107169932A (zh) 一种适用于中子成像系统图像的基于高斯‑泊松混合噪声模型的图像复原方法
CN105719254B (zh) 一种图像降噪方法及系统
CN106709876B (zh) 一种基于暗像元原理的光学遥感图像去雾方法
Zhang et al. Dehazing with improved heterogeneous atmosphere light estimation and a nonlinear color attenuation prior model
CN102339460B (zh) 卫星图像自适应复原方法
CN103337055A (zh) 一种基于梯度拟合的文本图像去模糊方法
CN111079689A (zh) 一种指纹图像增强方法
CN116342443A (zh) 一种近红外与可见光图像融合方法及系统
Wang et al. Adaptive Bright and Dark Channel Combined with Defogging Algorithm Based on Depth of Field
Jyothi A robust and efficient pre processing techniques for stereo images
Liu et al. The validity analysis of the non-local mean filter and a derived novel denoising method

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