CN117314791B - Infrared image cold reflection noise correction method based on Butterworth function fitting - Google Patents

Infrared image cold reflection noise correction method based on Butterworth function fitting Download PDF

Info

Publication number
CN117314791B
CN117314791B CN202311588699.3A CN202311588699A CN117314791B CN 117314791 B CN117314791 B CN 117314791B CN 202311588699 A CN202311588699 A CN 202311588699A CN 117314791 B CN117314791 B CN 117314791B
Authority
CN
China
Prior art keywords
image
infrared image
cold reflection
information
fitting
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
CN202311588699.3A
Other languages
Chinese (zh)
Other versions
CN117314791A (en
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.)
Changchun University of Science and Technology
Original Assignee
Changchun 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 Changchun University of Science and Technology filed Critical Changchun University of Science and Technology
Priority to CN202311588699.3A priority Critical patent/CN117314791B/en
Publication of CN117314791A publication Critical patent/CN117314791A/en
Application granted granted Critical
Publication of CN117314791B publication Critical patent/CN117314791B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10048Infrared image

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Image Processing (AREA)
  • Transforming Light Signals Into Electric Signals (AREA)

Abstract

本发明提出了基于巴特沃斯函数拟合的红外图像冷反射噪声矫正方法,首先设计平滑窗口对红外图像数据进行平滑处理,去掉多余的场景信息同时保留冷反射信号,并将其逆变换可以最大限度保留冷反射周围的信息,同时将其以外的信息进行抑制消除,使其图像场景信息达到最小;然后设计二维巴特沃斯函数,初始化需要拟合的参数,利用最小二乘法进行拟合,最后提取拟合函数的参数作为红外图像中冷反射的特征,利用二维巴特沃斯函数对映射后形成的曲面进行拟合,使用二维巴特沃斯函数的几何特征表示红外图像中冷反射的特征,通过与红外图像差分从而去除图像中的噪声,能够清楚地将图像与背景区分开,提高红外图像质量。

The present invention proposes a cold reflection noise correction method for infrared images based on Butterworth function fitting. First, a smoothing window is designed to smooth the infrared image data, remove redundant scene information while retaining the cold reflection signal, and inversely transform it to maximize the The information around the cold reflection is retained to the maximum extent, while the information other than it is suppressed and eliminated to minimize the image scene information; then a two-dimensional Butterworth function is designed, the parameters that need to be fitted are initialized, and the least squares method is used for fitting. Finally, the parameters of the fitting function are extracted as the characteristics of cold reflection in the infrared image, the two-dimensional Butterworth function is used to fit the surface formed after mapping, and the geometric characteristics of the two-dimensional Butterworth function are used to represent the characteristics of the cold reflection in the infrared image. Features, through difference with the infrared image to remove noise in the image, can clearly distinguish the image from the background and improve the quality of the infrared image.

Description

基于巴特沃斯函数拟合的红外图像冷反射噪声矫正方法Infrared image cold reflection noise correction method based on Butterworth function fitting

技术领域Technical field

本发明属于红外图像去噪技术领域,具体地,涉及基于巴特沃斯函数拟合的红外图像冷反射噪声矫正方法。The invention belongs to the technical field of infrared image denoising, and specifically relates to a method for correcting cold reflection noise of infrared images based on Butterworth function fitting.

背景技术Background technique

红外图像已经被广泛的应用于光学测量、目标识别、国防军工等一系列关系国计民生的方面,然而,红外探测器在成像过程中不可避免的会受到噪声的影响,冷反射现象作为红外图像中的一种缺陷,不仅影响红外图像的可视觉性,更会影响红外图像的目标探测、识别跟踪等后续处理应用。Infrared images have been widely used in a series of aspects related to the national economy and people's livelihood, such as optical measurement, target recognition, national defense and military industry. However, infrared detectors will inevitably be affected by noise during the imaging process. The cold reflection phenomenon is one of the main problems in infrared images. A defect that not only affects the visibility of infrared images, but also affects subsequent processing applications such as target detection, identification and tracking of infrared images.

目前采用非均匀性校正技术分为基于标定的方法以及基于场景的校正方法,这些方法对噪声有一定的抑制作用,但是对于由冷反射引起的噪声通过处理后会产生严重的伪影现象,对提高红外图像的质量的作用并不明显。The current non-uniformity correction technology is divided into calibration-based methods and scene-based correction methods. These methods have a certain inhibitory effect on noise, but the noise caused by cold reflection will produce serious artifacts after processing, which will cause serious artifacts. The effect of improving the quality of infrared images is not obvious.

发明内容Contents of the invention

针对上述问题,本发明提出了基于巴特沃斯函数拟合的红外图像冷反射噪声矫正方法,具体为基于巴特沃斯曲面拟合进行红外图像去噪,能够抑制冷反射所带来的影响,提高在低温环境下红外图像的信噪比,从而去除图像中的噪声。In response to the above problems, the present invention proposes a method for correcting cold reflection noise in infrared images based on Butterworth function fitting. Specifically, infrared image denoising is based on Butterworth surface fitting, which can suppress the impact of cold reflection and improve The signal-to-noise ratio of infrared images in low temperature environments, thereby removing noise in the images.

本发明通过以下技术方案实现:The present invention is realized through the following technical solutions:

一种基于巴特沃斯函数拟合的红外图像冷反射噪声矫正方法:A method for correcting cold reflection noise in infrared images based on Butterworth function fitting:

步骤一、通过传感器读取红外图像Step 1. Read the infrared image through the sensor ;

步骤二、设计平滑窗口对步骤一读取的红外图像进行平滑预处理;去掉多余的场景信息、并将其逆变换以最大限度保留冷反射周围的信息,同时将其以外的信息进行抑制消除,使其图像场景信息达到最小;Step 2. Design a smoothing window for the infrared image read in step 1. Perform smoothing preprocessing; remove excess scene information and inversely transform it to retain the information around the cold reflection to the maximum extent, while suppressing and eliminating information other than it to minimize the image scene information;

步骤三、构建二维巴特沃斯函数,初始化需要拟合的参数;Step 3: Construct a two-dimensional Butterworth function and initialize the parameters that need to be fitted;

步骤四、利用最小二乘法对数据进行拟合,最后提取拟合函数的参数作为红外图像中冷反射的特征;Step 4: Use the least squares method to fit the data, and finally extract the parameters of the fitting function as the characteristics of cold reflection in the infrared image;

步骤五、利用二维巴特沃斯函数对映射后形成的曲面进行拟合,使用二维巴特沃斯函数的几何特征表示红外图像中冷反射的特征,通过与红外图像差分从而去除图像中的噪声,完成红外图像冷反射噪声矫正。Step 5: Use the two-dimensional Butterworth function to fit the surface formed after mapping, use the geometric characteristics of the two-dimensional Butterworth function to represent the characteristics of cold reflection in the infrared image, and remove the noise in the image by difference with the infrared image. , complete the cold reflection noise correction of infrared images.

进一步地,在步骤二中,设计平滑窗口对红外图像的局部区域进行处理,以减少突出的噪声点或异常值;Further, in step two, a smoothing window is designed to process local areas of the infrared image to reduce prominent noise points or outliers;

所述预处理具体采用下式:The pretreatment specifically adopts the following formula:

预处理算法为:The preprocessing algorithm is: ;

其中,表示的是图像信息,/>表示的是均值图像的最大灰度值或者为的最大值,/>表示的信息像素/>与冷反射信号中心/>之间的距离即,/>表示的是均值图像中的标准差也就是抑制图像中场景信号的程度;in, Represents image information,/> It represents the maximum gray value of the mean image or The maximum value,/> Represented information pixels/> with cold reflection signal center/> The distance between ,/> It represents the standard deviation in the mean image, which is the degree of suppression of the scene signal in the image;

对红外图像进行预处理,达到平滑图像的目的,平滑后的图像定义为/>to infrared images Perform preprocessing to achieve the purpose of smoothing the image. The smoothed image is defined as/> .

进一步地,在步骤二中,通过所述逆变换来重建被平滑过程中丢失的细节,确保冷反射信号的重要特征不会因处理而丢失;Further, in step two, the inverse transformation is used to reconstruct the details lost during the smoothing process to ensure that important features of the cold reflection signal will not be lost due to processing;

并通过信息抑制来消除图像中与冷反射无关的信息,确保平滑后的图像只有关键的信息被保留。And use information suppression to eliminate information irrelevant to cold reflection in the image to ensure a smoothed image Only critical information is retained.

进一步地,在步骤三中,Further, in step three,

定义二维巴特沃斯函数Define the two-dimensional Butterworth function ;

式中表示幅度,/>表示属性,/>表示属性的均值,/>表示属性的标准差。in the formula Indicates the amplitude,/> Indicates attributes,/> Represents the mean value of the attribute,/> Represents the standard deviation of the attribute.

进一步地,在步骤四中,Further, in step four,

均初始化为1;利用最小二乘法矩阵运算,调整/>知道拟合精度达到阈值,获取拟合后的巴特沃斯拟合参数,所述巴特沃斯拟合参数视为红外图像中冷反射特征参数。Will are initialized to 1; use least squares matrix operation to adjust/> When the fitting accuracy reaches the threshold, the fitted Butterworth fitting parameters are obtained, and the Butterworth fitting parameters are regarded as cold reflection characteristic parameters in the infrared image.

进一步地,在步骤五中具体包括:Further, step five specifically includes:

利用步骤四中通过最小二乘法得到的巴特沃斯函数作为数学模型,拟合步骤二映射后形成的曲面,以捕捉和表示冷反射信号的几何特征;Use the Butterworth function obtained by the least squares method in step 4 as a mathematical model to fit the surface formed after mapping in step 2 to capture and represent the geometric characteristics of the cold reflection signal;

并使用二维巴特沃斯函数的几何特征表示红外图像中冷反射的特征:通过将拟合后的巴特沃斯函数与原始红外图像进行差分,去除图像中的噪声。And use the geometric characteristics of the two-dimensional Butterworth function to represent the characteristics of cold reflection in the infrared image: by difference between the fitted Butterworth function and the original infrared image, the noise in the image is removed.

获取步骤五的结果,提取作为含有冷反射噪声的红外图像/>的特征,实现对红外图像冷反射噪声的矫正。Get the result of step five and extract As an infrared image containing cold reflection noise/> characteristics to achieve the correction of cold reflection noise in infrared images.

一种基于巴特沃斯函数拟合的红外图像冷反射噪声矫正系统:A cold reflection noise correction system for infrared images based on Butterworth function fitting:

所述矫正系统包括:图像处理模块、拟合模块和矫正模块;The correction system includes: an image processing module, a fitting module and a correction module;

所述图像处理模块读取红外图像;并设计平滑窗口对读取的红外图像/>进行平滑预处理;去掉多余的场景信息、并将其逆变换以最大限度保留冷反射周围的信息,同时将其以外的信息进行抑制消除,使其图像场景信息达到最小;The image processing module reads the infrared image ;And design a smooth window for the read infrared image/> Perform smoothing preprocessing; remove excess scene information and inversely transform it to retain the information around the cold reflection to the maximum extent, while suppressing and eliminating information other than it to minimize the image scene information;

所述拟合模块构建二维巴特沃斯函数,初始化需要拟合的参数;利用最小二乘法对数据进行拟合,最后提取拟合函数的参数作为红外图像中冷反射的特征;The fitting module constructs a two-dimensional Butterworth function, initializes the parameters that need to be fitted; uses the least squares method to fit the data, and finally extracts the parameters of the fitting function as the characteristics of cold reflection in the infrared image;

所述矫正模块利用二维巴特沃斯函数对映射后形成的曲面进行拟合,使用二维巴特沃斯函数的几何特征表示红外图像中冷反射的特征,通过与红外图像差分从而去除图像中的噪声,完成红外图像冷反射噪声矫正。The correction module uses the two-dimensional Butterworth function to fit the curved surface formed after mapping, uses the geometric characteristics of the two-dimensional Butterworth function to represent the characteristics of cold reflection in the infrared image, and removes the cold reflection in the image by difference with the infrared image. Noise, complete the cold reflection noise correction of infrared images.

一种电子设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述方法的步骤。An electronic device includes a memory and a processor. The memory stores a computer program. When the processor executes the computer program, the steps of the above method are implemented.

一种计算机可读存储介质,用于存储计算机指令,所述计算机指令被处理器执行时实现上述方法的步骤。A computer-readable storage medium is used to store computer instructions that implement the steps of the above method when executed by a processor.

本发明有益效果Beneficial effects of the invention

本发明结合现有的红外噪声图像去噪方法,从噪声分布的角度对红外图像进行分析。利用二维巴特沃斯函数对映射后形成的曲面你和,使用二维巴特沃斯函数的几何特征表示红外图像数据特征。This invention combines existing infrared noise image denoising methods to analyze infrared images from the perspective of noise distribution. The curved surface formed by mapping the two-dimensional Butterworth function is used to represent the infrared image data characteristics using the geometric characteristics of the two-dimensional Butterworth function.

本发明方法对红外图像进行二维巴特沃斯曲面拟合并进行去噪处理。相比于其他红外图像分析方法,本发明能够明显提高红外图像的信噪比,从而去除图像中的噪声,提高红外图像的质量。The method of the present invention performs two-dimensional Butterworth surface fitting on infrared images and performs denoising processing. Compared with other infrared image analysis methods, the present invention can significantly improve the signal-to-noise ratio of the infrared image, thereby removing noise in the image and improving the quality of the infrared image.

附图说明Description of the drawings

图1为本发明的方法流程图。Figure 1 is a flow chart of the method of the present invention.

具体实施方式Detailed ways

下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some of the embodiments of the present invention, rather than all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts fall within the scope of protection of the present invention.

一种基于巴特沃斯函数拟合的红外图像冷反射噪声矫正方法:A method for correcting cold reflection noise in infrared images based on Butterworth function fitting:

步骤一、通过传感器或者其他设备读取红外图像Step 1. Read the infrared image through a sensor or other device ;

步骤二、设计平滑窗口对步骤一读取的红外图像进行平滑预处理;去掉多余的场景信息、并将其逆变换以最大限度保留冷反射周围的信息,同时将其以外的信息进行抑制消除,使其图像场景信息达到最小;使得原始图像中的冷反射信号更为明显。Step 2. Design a smoothing window for the infrared image read in step 1. Perform smoothing preprocessing; remove redundant scene information and inversely transform it to retain the information around the cold reflection to the maximum extent, while suppressing and eliminating information other than it to minimize the image scene information; making the cold reflection in the original image The reflected signal is more obvious.

在步骤二中,设计平滑窗口对红外图像的局部区域进行处理,以减少突出的噪声点或异常值;In step two, a smoothing window is designed to process local areas of the infrared image to reduce prominent noise points or outliers;

所述预处理具体采用下式:The pretreatment specifically adopts the following formula:

预处理算法为:The preprocessing algorithm is: ;

其中,表示的是图像信息,/>表示的是均值图像的最大灰度值或者为的最大值,/>表示的信息像素/>与冷反射信号中心/>之间的距离即,/>表示的是均值图像中的标准差也就是抑制图像中场景信号的程度;in, Represents image information,/> It represents the maximum gray value of the mean image or The maximum value,/> Represented information pixels/> with cold reflection signal center/> The distance between ,/> It represents the standard deviation in the mean image, which is the degree of suppression of the scene signal in the image;

对红外图像进行预处理,达到平滑图像的目的,平滑后的图像定义为/>to infrared images Perform preprocessing to achieve the purpose of smoothing the image. The smoothed image is defined as/> .

在步骤二中,通过所述逆变换来重建被平滑过程中丢失的细节,确保冷反射信号的重要特征不会因处理而丢失;In step two, the inverse transformation is used to reconstruct the details lost during the smoothing process to ensure that important features of the cold reflection signal will not be lost due to processing;

并通过信息抑制(通过特定的滤波器或者处理技术来实现)来消除图像中与冷反射无关的信息,确保平滑后的图像只有关键的信息被保留,以便更好地突出冷反射信号。And through information suppression (implemented through specific filters or processing techniques) to eliminate information unrelated to cold reflection in the image to ensure a smooth image Only critical information is retained to better highlight cold reflection signals.

步骤三、构建二维巴特沃斯函数,初始化需要拟合的参数;Step 3: Construct a two-dimensional Butterworth function and initialize the parameters that need to be fitted;

在步骤三中,In step three,

定义二维巴特沃斯函数Define the two-dimensional Butterworth function ;

式中表示幅度,/>表示属性,/>表示属性的均值,/>表示属性的标准差。in the formula Indicates the amplitude,/> Indicates attributes,/> Represents the mean value of the attribute,/> Represents the standard deviation of the attribute.

步骤四、利用最小二乘法对数据进行拟合,最后提取拟合函数的参数作为红外图像中冷反射的特征,实现对红外图像噪声特征的提取;Step 4: Use the least squares method to fit the data, and finally extract the parameters of the fitting function as the cold reflection features in the infrared image to extract the noise features of the infrared image;

在步骤四中,在拟合过程中,通过最小化拟合函数与实际数据之间的误差调整巴特沃斯函数的参数,以使该函数更好地拟合红外图像中的冷反射信号,最小二乘法会尝试不断地调整巴特沃斯函数的参数,直到达到使误差最小化的状态;完成拟合后,最终得到的巴特沃斯函数的参数将被视为红外图像中冷反射特征的表示。In step four, during the fitting process, the parameters of the Butterworth function are adjusted by minimizing the error between the fitting function and the actual data, so that the function can better fit the cold reflection signal in the infrared image, with a minimum The square method will try to continuously adjust the parameters of the Butterworth function until it reaches a state that minimizes the error; after the fitting is completed, the parameters of the final Butterworth function will be regarded as a representation of the cold reflection features in the infrared image.

均初始化为1;利用最小二乘法矩阵运算,调整/>知道拟合精度达到阈值,获取拟合后的巴特沃斯拟合参数,所述巴特沃斯拟合参数视为红外图像中冷反射特征参数。Will are initialized to 1; use least squares matrix operation to adjust/> When the fitting accuracy reaches the threshold, the fitted Butterworth fitting parameters are obtained, and the Butterworth fitting parameters are regarded as cold reflection characteristic parameters in the infrared image.

即在步骤四中,通过最小二乘法,巴特沃斯函数的参数已经被调整以最佳拟合原始红外图像中的冷反射信号。这个巴特沃斯函数是一个数学模型,其中的参数是通过拟合过程得到的。That is, in step four, through the least squares method, the parameters of the Butterworth function have been adjusted to best fit the cold reflection signal in the original infrared image. This Butterworth function is a mathematical model whose parameters are obtained through a fitting process.

步骤五、利用二维巴特沃斯函数对映射后形成的曲面进行拟合,使用二维巴特沃斯函数的几何特征表示红外图像中冷反射的特征,通过与红外图像差分从而去除图像中的噪声,能够清楚地将图像与背景区分开,提高红外图像质量,完成红外图像冷反射噪声矫正。Step 5: Use the two-dimensional Butterworth function to fit the surface formed after mapping, use the geometric characteristics of the two-dimensional Butterworth function to represent the characteristics of cold reflection in the infrared image, and remove the noise in the image by difference with the infrared image. , can clearly distinguish the image from the background, improve the quality of the infrared image, and complete the cold reflection noise correction of the infrared image.

在步骤五中具体包括:Step five specifically includes:

利用步骤四中通过最小二乘法得到的巴特沃斯函数作为数学模型,拟合步骤二映射后形成的曲面,以捕捉和表示冷反射信号的几何特征;Use the Butterworth function obtained by the least squares method in step 4 as a mathematical model to fit the surface formed after mapping in step 2 to capture and represent the geometric characteristics of the cold reflection signal;

并使用二维巴特沃斯函数的几何特征表示红外图像中冷反射的特征:通过将拟合后的巴特沃斯函数与原始红外图像进行差分,去除图像中的噪声。And use the geometric characteristics of the two-dimensional Butterworth function to represent the characteristics of cold reflection in the infrared image: by difference between the fitted Butterworth function and the original infrared image, the noise in the image is removed.

这是因为拟合后的函数主要捕捉了冷反射信号的特征,而与冷反射无关的信号被去除,从而提高了图像的质量。通过上述差分操作,去除了噪声,使得图像更为清晰。拟合的巴特沃斯函数的几何特征有助于将冷反射信号与图像背景区分开,使得目标信号更为突出。综合上述步骤,整体的目标是提高红外图像的质量,同时通过拟合巴特沃斯函数和差分操作,完成对冷反射噪声的矫正。This is because the fitted function mainly captures the characteristics of the cold reflection signal, while signals unrelated to the cold reflection are removed, thereby improving the quality of the image. Through the above differential operation, noise is removed, making the image clearer. The geometric characteristics of the fitted Butterworth function help distinguish the cold reflection signal from the image background, making the target signal more prominent. Based on the above steps, the overall goal is to improve the quality of infrared images, and at the same time complete the correction of cold reflection noise through fitting Butterworth functions and differential operations.

完成这些步骤后,得到的红外图像应当更清晰,冷反射信号更为凸显,同时噪声得到有效的去除。After completing these steps, the resulting infrared image should be clearer, the cold reflection signal should be more prominent, and the noise should be effectively removed.

获取步骤五的结果,提取作为含有冷反射噪声的红外图像/>的特征,实现对红外图像冷反射噪声的矫正。Get the result of step five and extract As an infrared image containing cold reflection noise/> characteristics to achieve the correction of cold reflection noise in infrared images.

一种基于巴特沃斯函数拟合的红外图像冷反射噪声矫正系统:A cold reflection noise correction system for infrared images based on Butterworth function fitting:

所述矫正系统包括:图像处理模块、拟合模块和矫正模块;The correction system includes: an image processing module, a fitting module and a correction module;

所述图像处理模块通过传感器或者其他设备读取红外图像;并设计平滑窗口对读取的红外图像/>进行平滑预处理;去掉多余的场景信息、并将其逆变换以最大限度保留冷反射周围的信息,同时将其以外的信息进行抑制消除,使其图像场景信息达到最小;使得原始图像中的冷反射信号更为明显。The image processing module reads infrared images through sensors or other devices ;And design a smooth window for the read infrared image/> Perform smoothing preprocessing; remove redundant scene information and inversely transform it to retain the information around the cold reflection to the maximum extent, while suppressing and eliminating information other than it to minimize the image scene information; making the cold reflection in the original image The reflected signal is more obvious.

所述拟合模块构建二维巴特沃斯函数,初始化需要拟合的参数;利用最小二乘法对数据进行拟合,最后提取拟合函数的参数作为红外图像中冷反射的特征,实现对红外图像噪声特征的提取;The fitting module constructs a two-dimensional Butterworth function, initializes the parameters that need to be fitted; uses the least squares method to fit the data, and finally extracts the parameters of the fitting function as the characteristics of cold reflection in the infrared image to realize the analysis of the infrared image. Extraction of noise features;

所述矫正模块利用二维巴特沃斯函数对映射后形成的曲面进行拟合,使用二维巴特沃斯函数的几何特征表示红外图像中冷反射的特征,通过与红外图像差分从而去除图像中的噪声,能够清楚地将图像与背景区分开,提高红外图像质量,完成红外图像冷反射噪声矫正。The correction module uses the two-dimensional Butterworth function to fit the curved surface formed after mapping, uses the geometric characteristics of the two-dimensional Butterworth function to represent the characteristics of cold reflection in the infrared image, and removes the cold reflection in the image by difference with the infrared image. Noise, it can clearly distinguish the image from the background, improve the quality of the infrared image, and complete the cold reflection noise correction of the infrared image.

一种电子设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述方法的步骤。An electronic device includes a memory and a processor. The memory stores a computer program. When the processor executes the computer program, the steps of the above method are implemented.

一种计算机可读存储介质,用于存储计算机指令,所述计算机指令被处理器执行时实现上述方法的步骤。A computer-readable storage medium is used to store computer instructions that implement the steps of the above method when executed by a processor.

本申请实施例中的存储器可以是易失性存储器或非易失性存储器,或可包括易失性和非易失性存储器两者。其中,非易失性存储器可以是只读存储器read only memory,ROM、可编程只读存储器programmable ROM,PROM、可擦除可编程只读存储器erasablePROM,EPROM、电可擦除可编程只读存储器electrically EPROM,EEPROM或闪存。易失性存储器可以是随机存取存储器random access memory,RAM,其用作外部高速缓存。通过示例性但不是限制性说明,许多形式的RAM 可用,例如静态随机存取存储器static RAM,SRAM、动态随机存取存储器dynamic RAM,DRAM、同步动态随机存取存储器synchronous DRAM,SDRAM、双倍数据速率同步动态随机存取存储器double data rate SDRAM,DDR SDRAM、增强型同步动态随机存取存储器enhanced SDRAM,ESDRAM、同步连接动态随机存取存储器synchlink DRAM,SLDRAM和直接内存总线随机存取存储器direct rambus RAM,DR RAM。应注意,本发明描述的方法的存储器旨在包括但不限于这些和任意其它适合类型的存储器。The memory in the embodiment of the present application may be a volatile memory or a non-volatile memory, or may include both volatile and non-volatile memories. Among them, the non-volatile memory can be read only memory, ROM, programmable ROM, PROM, erasable programmable read-only memory erasablePROM, EPROM, electrically erasable programmable read-only memory. EPROM, EEPROM or flash memory. Volatile memory can be random access memory, RAM, which is used as an external cache. By way of illustration, but not limitation, many forms of RAM are available, such as static RAM, SRAM, dynamic RAM, DRAM, synchronous DRAM, SDRAM, double data rate synchronous dynamic random access memory double data rate SDRAM, DDR SDRAM, enhanced synchronous dynamic random access memory enhanced SDRAM, ESDRAM, synchlink dynamic random access memory synchlink DRAM, SLDRAM and direct rambus RAM ,DR RAM. It should be noted that the memory of the method described herein is intended to include, but is not limited to, these and any other suitable types of memory.

在上述实施例中,可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用软件实现时,可以全部或部分地以计算机程序产品的形式实现。所述计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行所述计算机指令时,全部或部分地产生按照本申请实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中,或者从一个计算机可读存储介质向另一个计算机可读存储介质传输,例如,所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线例如同轴电缆、光纤、数字用户线digital subscriber line,DSL或无线例如红外、无线、微波等方式向另一个网站站点、计算机、服务器或数据中心进行传输。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质例如,软盘、硬盘、磁带、光介质例如,高密度数字视频光盘digital video disc,DVD、或者半导体介质例如,固态硬盘solid state disc,SSD等。In the above embodiments, it may be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented using software, it may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer instructions are loaded and executed on the computer, the processes or functions described in the embodiments of the present application are generated in whole or in part. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions may be stored in or transmitted from one computer-readable storage medium to another, e.g., the computer instructions may be transferred from a website, computer, server, or data center Transmission to another website, computer, server or data center through wired means such as coaxial cable, optical fiber, digital subscriber line, DSL or wireless means such as infrared, wireless, microwave, etc. The computer-readable storage medium may be any available medium that can be accessed by a computer or a data storage device such as a server, data center, etc. that contains one or more available media integrated. The available media may be magnetic media such as floppy disks, hard disks, magnetic tapes, optical media such as high-density digital video discs (DVD), or semiconductor media such as solid state disks (SSD).

在实现过程中,上述方法的各步骤可以通过处理器中的硬件的集成逻辑电路或者软 件形式的指令完成。结合本申请实施例所公开的方法的步骤可以直接体现为硬件处理器执行完成,或者用处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器,闪存、只读存储器,可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器,处理器读取存储器中的信息,结合其硬件完成上述方法的步骤。为避免重复,这里不再详细描述。During the implementation process, each step of the above method can be completed through the integrated logic circuit of the hardware in the processor or instructions in the form of software. The steps of the methods disclosed in conjunction with the embodiments of the present application can be directly implemented by a hardware processor for execution, or can be executed by a combination of hardware and software modules in the processor. The software module can be located in random access memory, flash memory, read-only memory, programmable read-only memory or electrically erasable programmable memory, registers and other mature storage media in this field. The storage medium is located in the memory, and the processor reads the information in the memory and completes the steps of the above method in combination with its hardware. To avoid repetition, it will not be described in detail here.

应注意,本申请实施例中的处理器可以是一种集成电路芯片,具有信号处理能力。在实现过程中,上述方法实施例的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器可以是通用处理器、数字信号处理器DSP、专用集成电路ASIC、现场可编程门阵列FPGA或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本申请实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成,或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器,闪存、只读存储器,可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器,处理器读取存储器中的信息,结合其硬件完成上述方法的步骤。It should be noted that the processor in the embodiment of the present application may be an integrated circuit chip with signal processing capabilities. During the implementation process, each step of the above method embodiment can be completed through an integrated logic circuit of hardware in the processor or instructions in the form of software. The above-mentioned processor may be a general-purpose processor, a digital signal processor DSP, an application-specific integrated circuit ASIC, a field programmable gate array FPGA or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components. Each method, step and logical block diagram disclosed in the embodiment of this application can be implemented or executed. A general-purpose processor may be a microprocessor or the processor may be any conventional processor, etc. The steps of the method disclosed in conjunction with the embodiments of the present application can be directly implemented by a hardware decoding processor, or executed by a combination of hardware and software modules in the decoding processor. The software module can be located in random access memory, flash memory, read-only memory, programmable read-only memory or electrically erasable programmable memory, registers and other mature storage media in this field. The storage medium is located in the memory, and the processor reads the information in the memory and completes the steps of the above method in combination with its hardware.

以上对本发明所提出的基于巴特沃斯函数拟合的红外图像冷反射噪声矫正方法,进行了详细介绍,对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。The cold reflection noise correction method for infrared images based on Butterworth function fitting proposed by the present invention has been introduced in detail above, and the principles and implementations of the present invention have been explained. The description of the above embodiments is only used to help understand the present invention. The method and its core idea; at the same time, for those of ordinary skill in the field, there will be changes in the specific implementation and application scope according to the idea of the present invention. In summary, the contents of this specification should not be understood as Limitations on the invention.

Claims (7)

1. An infrared image cold reflection noise correction method based on Butterworth function fitting is characterized by comprising the following steps of:
step one, reading an infrared image I through a sensor;
step two, designing a smoothing window to carry out smoothing pretreatment on the infrared image I read in the step one; removing redundant scene information, inversely transforming the scene information to keep information around cold reflection to the maximum extent, and simultaneously suppressing and eliminating other information to minimize the image scene information;
in the second step, a smooth window is designed to process the local area of the infrared image so as to reduce the prominent noise point or abnormal value;
the pretreatment specifically adopts the following formula:
the preprocessing algorithm is as follows:
wherein f (x, y) represents image information, max { f } represents the maximum gray value of the mean image or the maximum value of f (x, y), and the distance between the information pixel (x, y) represented by R (x, y) and the cold reflection signal center (a, b), i.e., R (x, y) = | (x-a, y-b) | 2 Delta represents the standard deviation in the mean image, i.e., the extent to which the scene signal in the image is suppressed;
preprocessing an infrared image I to achieve the purpose of smoothing the image, wherein the smoothed image is defined as I';
in the second step, reconstructing the details lost in the smoothed process through inverse transformation, so as to ensure that important characteristics of the cold reflection signals are not lost due to processing;
the information irrelevant to cold reflection in the image is eliminated through information inhibition, so that only key information of the smoothed image I' is ensured to be reserved;
step three, constructing a two-dimensional Butterworth function, and initializing parameters to be fitted;
fitting the data by using a least square method, and finally extracting parameters of a fitting function as characteristics of cold reflection in the infrared image;
fitting the curved surface formed after mapping by using a two-dimensional Butterworth function, using the geometric features of the two-dimensional Butterworth function to represent the features of cold reflection in the infrared image, and removing noise in the image by differentiating the geometric features with the infrared image to finish the correction of the cold reflection noise of the infrared image.
2. The method according to claim 1, wherein: in the third step of the process, the process is carried out,
definition of two-dimensional Butterworth function
Wherein A represents amplitude, mu 12 Representing the mean, sigma of the attributes xy Representing the standard deviation of the properties.
3. The method according to claim 2, characterized in that: in the fourth step of the process, the process is carried out,
will A, mu 12xy All initialized to 1; adjusting A, mu by least square matrix operation 12xy And until the fitting precision reaches a threshold value, acquiring the fitted Butterworth fitting parameters, wherein the Butterworth fitting parameters are regarded as cold reflection characteristic parameters in the infrared image.
4. A method according to claim 3, characterized in that:
the fifth step specifically comprises:
fitting the curved surface formed after the mapping in the step two by using the Butterworth function obtained by the least square method in the step four as a mathematical model so as to capture and represent the geometric characteristics of the cold reflection signals;
and using the geometric features of the two-dimensional butterworth function to represent the features of cold reflection in the infrared image: the noise in the image is removed by differentiating the fitted Butterworth function with the original infrared image;
obtaining the result of the fifth step, extracting (A, mu) 12xy ) As a feature of the infrared image I containing the cold reflection noise, correction of the cold reflection noise of the infrared image is realized.
5. An infrared image cold reflection noise correction system based on Butterworth function fitting is characterized in that:
the corrective system includes: the device comprises an image processing module, a fitting module and a correction module;
the image processing module reads an infrared image I; a smoothing window is designed to carry out smoothing pretreatment on the read infrared image I; removing redundant scene information, inversely transforming the scene information to keep information around cold reflection to the maximum extent, and simultaneously suppressing and eliminating other information to minimize the image scene information;
designing a smooth window to process a local area of the infrared image so as to reduce prominent noise points or abnormal values;
the pretreatment specifically adopts the following formula:
the preprocessing algorithm is as follows:
where f (x, y) represents image information, M represents a maximum gray value of the mean image or a maximum value of f (x, y),the distance between the information pixel (x, y) represented by R (x, y) and the cold reflection signal center (a, b), i.e., R (x, y) = | (x-a, y-b) || 2 Delta represents the standard deviation in the mean image, i.e., the extent to which the scene signal in the image is suppressed;
preprocessing an infrared image I to achieve the purpose of smoothing the image, wherein the smoothed image is defined as I';
reconstructing lost details in the smoothed process through inverse transformation, so as to ensure that important characteristics of the cold reflection signal are not lost due to processing;
the information irrelevant to cold reflection in the image is eliminated through information inhibition, so that only key information of the smoothed image I' is ensured to be reserved;
the fitting module constructs a two-dimensional Butterworth function, and initializes parameters to be fitted; fitting the data by using a least square method, and finally extracting parameters of a fitting function as characteristics of cold reflection in the infrared image;
the correction module fits the curved surface formed after mapping by using a two-dimensional Butterworth function, uses the geometric features of the two-dimensional Butterworth function to represent the features of cold reflection in the infrared image, and removes noise in the image by differentiating the features with the infrared image, so that the correction of the cold reflection noise of the infrared image is completed.
6. An electronic device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor implements the steps of the method of any of claims 1 to 4 when the computer program is executed.
7. A computer readable storage medium storing computer instructions which, when executed by a processor, implement the steps of the method of any one of claims 1 to 4.
CN202311588699.3A 2023-11-27 2023-11-27 Infrared image cold reflection noise correction method based on Butterworth function fitting Active CN117314791B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311588699.3A CN117314791B (en) 2023-11-27 2023-11-27 Infrared image cold reflection noise correction method based on Butterworth function fitting

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311588699.3A CN117314791B (en) 2023-11-27 2023-11-27 Infrared image cold reflection noise correction method based on Butterworth function fitting

Publications (2)

Publication Number Publication Date
CN117314791A CN117314791A (en) 2023-12-29
CN117314791B true CN117314791B (en) 2024-02-20

Family

ID=89281400

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311588699.3A Active CN117314791B (en) 2023-11-27 2023-11-27 Infrared image cold reflection noise correction method based on Butterworth function fitting

Country Status (1)

Country Link
CN (1) CN117314791B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103592367A (en) * 2013-10-23 2014-02-19 江苏大学 Portable poultry egg quality detection device and method
CN109035362A (en) * 2018-06-11 2018-12-18 西安电子科技大学 Cold emission removing method based on cold emission strength model
CN109191387A (en) * 2018-07-20 2019-01-11 河南师范大学 A kind of Infrared Image Denoising method based on Butterworth filter
CN109829861A (en) * 2018-12-29 2019-05-31 西安电子科技大学 A kind of cold emission suppressing method based on wavelet decomposition
CN111738933A (en) * 2020-05-15 2020-10-02 南京邮电大学 Method and device for noise suppression in infrared digital holographic reconstruction
CN113326855A (en) * 2021-06-22 2021-08-31 长光卫星技术有限公司 Night lamplight remote sensing image feature extraction method based on two-dimensional Gaussian surface fitting

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015133130A1 (en) * 2014-03-06 2015-09-11 日本電気株式会社 Video capturing device, signal separation device, and video capturing method
US10740889B2 (en) * 2017-12-29 2020-08-11 Huizhou China Star Optoelectronics Technology Co., Ltd. Method and system for detection of in-panel mura based on hough transform and gaussian fitting

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103592367A (en) * 2013-10-23 2014-02-19 江苏大学 Portable poultry egg quality detection device and method
CN109035362A (en) * 2018-06-11 2018-12-18 西安电子科技大学 Cold emission removing method based on cold emission strength model
CN109191387A (en) * 2018-07-20 2019-01-11 河南师范大学 A kind of Infrared Image Denoising method based on Butterworth filter
CN109829861A (en) * 2018-12-29 2019-05-31 西安电子科技大学 A kind of cold emission suppressing method based on wavelet decomposition
CN111738933A (en) * 2020-05-15 2020-10-02 南京邮电大学 Method and device for noise suppression in infrared digital holographic reconstruction
CN113326855A (en) * 2021-06-22 2021-08-31 长光卫星技术有限公司 Night lamplight remote sensing image feature extraction method based on two-dimensional Gaussian surface fitting

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于傅里叶变换的森林红外图像增强算法研究;崔帅;刘波;丁德红;;红外技术(第01期);全文 *

Also Published As

Publication number Publication date
CN117314791A (en) 2023-12-29

Similar Documents

Publication Publication Date Title
CN107301661B (en) High-resolution remote sensing image registration method based on edge point features
CN108830818B (en) Rapid multi-focus image fusion method
CN109242888B (en) Infrared and visible light image fusion method combining image significance and non-subsampled contourlet transformation
US10650260B2 (en) Perspective distortion characteristic based facial image authentication method and storage and processing device thereof
JP6446374B2 (en) Improvements in image processing or improvements related to image processing
WO2023193401A1 (en) Point cloud detection model training method and apparatus, electronic device, and storage medium
CN112598708B (en) A hyperspectral target tracking method based on four-feature fusion and weight coefficient
CN109509163B (en) A method and system for multi-focus image fusion based on FGF
CN112950685B (en) Infrared and visible light image registration method, system and storage medium
WO2023019555A1 (en) Cell fluorescence image thresholding method and system, terminal, and storage medium
CN106934398B (en) Image de-noising method based on super-pixel cluster and rarefaction representation
CN113421205A (en) Small target detection method combined with infrared polarization imaging
CN111695575A (en) Weld image feature point extraction method based on improved mean filtering method
CN109064402B (en) Single Image Super-Resolution Reconstruction Method Based on Enhanced Non-Local Total Variational Model Prior
CN106651923A (en) Method and system for video image target detection and segmentation
Zhang et al. Depth enhancement with improved exemplar-based inpainting and joint trilateral guided filtering
CN117314791B (en) Infrared image cold reflection noise correction method based on Butterworth function fitting
CN115601569A (en) A method and system for optimal matching of heterogeneous images based on improved PIIFD
CN112927169B (en) A Noise Removal Method for Remote Sensing Image Based on Wavelet Transform and Improved Weighted Kernel Norm Minimization
CN109949337A (en) Moving target detection method and device based on Gaussian mixture background model
CN106778822B (en) Image straight line detection method based on funnel transformation
CN106651781B (en) A kind of image noise suppression method of Laser active illuminated imaging
CN113470074B (en) Self-adaptive space-time regularization target tracking method based on block discrimination
CN116596801A (en) Image non-local mean denoising method and device
CN116596958B (en) Target tracking method and device based on online sample augmentation

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