[go: up one dir, main page]

CN110298806B - Infrared image enhancement method and system - Google Patents

Infrared image enhancement method and system Download PDF

Info

Publication number
CN110298806B
CN110298806B CN201910601556.9A CN201910601556A CN110298806B CN 110298806 B CN110298806 B CN 110298806B CN 201910601556 A CN201910601556 A CN 201910601556A CN 110298806 B CN110298806 B CN 110298806B
Authority
CN
China
Prior art keywords
band
sub
infrared image
enhanced
representing
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
CN201910601556.9A
Other languages
Chinese (zh)
Other versions
CN110298806A (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.)
Shanghai Rongjun Technology Co ltd
Original Assignee
Shanghai Rongjun Technology Co ltd
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 Shanghai Rongjun Technology Co ltd filed Critical Shanghai Rongjun Technology Co ltd
Priority to CN201910601556.9A priority Critical patent/CN110298806B/en
Publication of CN110298806A publication Critical patent/CN110298806A/en
Application granted granted Critical
Publication of CN110298806B publication Critical patent/CN110298806B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/10Image enhancement or restoration using non-spatial domain filtering
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration using local operators
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/136Segmentation; Edge detection involving thresholding
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10048Infrared image
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20056Discrete and fast Fourier transform, [DFT, FFT]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20064Wavelet transform [DWT]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging

Landscapes

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

Abstract

The invention provides an infrared image enhancement method and system, wherein the method comprises the following steps: acquiring a stationary wavelet transform result and a saliency map corresponding to the infrared image; wherein the stationary wavelet transform result comprises: a low-pass sub-band and a detail sub-band; carrying out contrast increasing processing on the low-pass sub-band and carrying out linear enhancement processing on the detail sub-band to obtain an enhanced low-pass sub-band and an enhanced detail sub-band; performing low-pass filtering and threshold segmentation processing on the saliency map to obtain a processed saliency map; performing dot multiplication on the enhanced detail sub-bands and the processed saliency map to obtain fused detail sub-bands; and performing stationary wavelet inverse transformation based on the enhanced low-pass sub-band and the fused detail sub-band to obtain an enhanced infrared image. The invention can effectively improve the infrared image enhancement processing efficiency and the real-time property.

Description

Infrared image enhancement method and system
Technical Field
The invention relates to the technical field of image processing, in particular to an infrared image enhancement method and system.
Background
The image enhancement technology is mainly used for improving the recognition capability of an image, selectively highlighting interesting features in the image and subtracting redundant features. Due to the characteristics of the infrared imaging technology, the infrared image generally has the problems of low contrast, blurred edge and unobvious details, so that an image enhancement technology is required to improve the quality of the infrared image.
At present, the infrared image enhancement methods suitable for the market are mainly two types, namely an image enhancement method based on a space domain and an image enhancement method based on a transform domain. The spatial domain image enhancement method is to directly process the pixels in the image, and is basically based on the gray mapping transformation of the image, and the type of the mapping transformation used depends on the purpose of enhancement. The spatial domain image enhancement method mainly comprises gray level transformation, Histogram Equalization (HE), smoothing, sharpening and the like. The transform domain enhancement method firstly converts an image in an image space into other spaces in a certain form, then performs image enhancement processing by using the special property of the space, and finally converts the image into the original image space, thereby obtaining an enhanced image. Common variations are wavelet transforms, fourier transforms, discrete cosine transforms, and the like.
However, the space complexity and the time complexity of the conventional image enhancement method are high, and the real-time requirement of the system is difficult to meet.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide an infrared image enhancement method and system for improving the infrared image enhancement processing efficiency and instantaneity.
In a first aspect, the present invention provides an infrared image enhancement method, including:
s1: acquiring a stationary wavelet transform result and a saliency map corresponding to the infrared image; wherein the stationary wavelet transform result comprises: a low-pass sub-band and a detail sub-band;
s2: carrying out contrast increasing processing on the low-pass sub-band and carrying out linear enhancement processing on the detail sub-band to obtain an enhanced low-pass sub-band and an enhanced detail sub-band;
s3: performing low-pass filtering and threshold segmentation processing on the saliency map to obtain a processed saliency map;
s4: performing dot multiplication on the enhanced detail sub-bands and the processed saliency map to obtain fused detail sub-bands;
s5: and performing stationary wavelet inverse transformation based on the enhanced low-pass sub-band and the fused detail sub-band to obtain an enhanced infrared image.
Optionally, the performing, in S2, a contrast increasing process on the low-pass sub-band includes:
performing CLAHE contrast enhancement processing on the low-pass sub-band, wherein the processing algorithm is as follows:
L′=(1-α)L+α·CLAHE(L)
wherein: l denotes a low-pass subband, L' denotes an enhanced low-pass subband, α denotes a control strength, clahe (L) denotes a limit contrast adaptive histogram equalization, and · denotes a multiplication operation.
Optionally, the processing algorithm for performing linear enhancement processing on the detail sub-bands in S2 is as follows:
Figure BDA0002118647060000021
Figure BDA0002118647060000022
Figure BDA0002118647060000023
v=2.5s/t·M
wherein: s represents the coefficient magnitude in the transform domain, M represents the maximum coefficient magnitude, and t represents the first enhancement parameter; hhRepresenting horizontal direction detail sub-bands, sign (·) representing sign function operation, b representing a second enhancement parameter, v representing an intermediate parameter, tanh (·) representing hyperbolic tangent function operation, c representing a third enhancement parameter, exp (·) representing exponential function operation,
Figure BDA0002118647060000024
representing the enhanced horizontal subband, HvShowing the sub-bands of detail in the vertical direction,
Figure BDA0002118647060000025
representing the enhanced vertical detail subband, HdThe detail sub-bands in the diagonal direction are shown,
Figure BDA0002118647060000026
representing the enhanced diagonal detail sub-bands.
Optionally, in S4, the significance map is multiplied by each detail subband element by element to obtain a calculation formula of the fused detail subband, where the calculation formula is as follows:
Figure BDA0002118647060000027
Figure BDA0002118647060000028
Figure BDA0002118647060000029
wherein, H'hRepresents the horizontal detail sub-band, H ', of the final enhancement'vRepresents the vertical-direction detail sub-band, H ', of the final enhancement'dThe final enhanced diagonal detail subband is represented, S represents the saliency map, and o represents the dot product operation of the matrix.
In a second aspect, the present invention further provides an infrared image enhancement system, comprising: the device comprises a memory and a processor, wherein the memory stores executable instructions of the processor; wherein the processor is configured to perform the infrared image enhancement method of any one of the first aspects via execution of the executable instructions.
Compared with the prior art, the embodiment of the invention has the following beneficial effects:
according to the infrared image enhancement method and system, a stationary wavelet transform result and a saliency map corresponding to an infrared image are obtained; wherein the stationary wavelet transform result comprises: a low-pass sub-band and a detail sub-band; carrying out contrast increasing processing on the low-pass sub-band and carrying out linear enhancement processing on the detail sub-band to obtain an enhanced low-pass sub-band and an enhanced detail sub-band; performing low-pass filtering and threshold segmentation processing on the saliency map to obtain a processed saliency map; performing dot multiplication on the enhanced detail sub-bands and the processed saliency map to obtain fused detail sub-bands; and performing stationary wavelet inverse transformation based on the enhanced low-pass sub-band and the fused detail sub-band to obtain an enhanced infrared image. The embodiment of the invention utilizes the stationary wavelet transform of the saliency map, can carry out multi-scale and multi-direction decomposition on the image, well represents the detail information of the image, and is beneficial to the detail enhancement of the infrared image; in addition, the saliency map can be very robust to infrared images containing noise; the efficiency is higher, the real-time is good.
Drawings
Other features, objects and advantages of the invention will become more apparent upon reading of the detailed description of non-limiting embodiments with reference to the following drawings:
FIG. 1 is a schematic diagram of an infrared image enhancement method according to an embodiment of the present invention;
FIG. 2(a) is a first infrared image;
FIG. 2(b) is a first infrared image after AMSR enhancement processing;
FIG. 2(c) is a first infrared image after BPDHE enhancement processing;
FIG. 2(d) is a first infrared image after CRM enhancement processing;
fig. 2(e) is a first infrared image after the EFF enhancement processing;
FIG. 2(f) is a first infrared image after homomorphic filter enhancement processing;
FIG. 2(g) is a first infrared image after Natural Factor enhancement processing;
FIG. 2(h) is a first infrared image after enhancement processing according to an embodiment of the present invention;
FIG. 3(a) is a second infrared image;
FIG. 3(b) is a second infrared image after AMSR enhancement processing;
FIG. 3(c) is a second infrared image after BPDHE enhancement;
FIG. 3(d) is a second infrared image after CRM enhancement processing;
fig. 3(e) is a second infrared image after the EFF enhancement process;
FIG. 3(f) is a second infrared image after homomorphic filtering enhancement processing;
FIG. 3(g) is a second infrared image after Natural Factor enhancement processing;
fig. 3(h) is a second infrared image after enhancement processing according to an embodiment of the present invention.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the invention, but are not intended to limit the invention in any way. It should be noted that it would be obvious to those skilled in the art that various changes and modifications can be made without departing from the spirit of the invention. All falling within the scope of the present invention.
Fig. 1 is a schematic diagram illustrating an infrared image enhancement method according to an embodiment of the present invention, and as shown in fig. 1, first, a stationary wavelet transform and a saliency map calculation are respectively performed on an input image; secondly, respectively carrying out contrast enhancement and nonlinear enhancement on the low-pass sub-band and the detail sub-band after the stationary wavelet transform; meanwhile, low-pass filtering and threshold segmentation are carried out on the acquired saliency map; then, carrying out dot multiplication on the result after nonlinear enhancement and the result after threshold segmentation to obtain a new detail sub-band; and finally, performing inverse stationary wavelet transform on the low-pass sub-band after contrast enhancement and the new detail sub-band to obtain an enhanced infrared image.
With reference to fig. 1, the infrared image enhancement method provided in this embodiment may include:
step 1: acquiring a stationary wavelet transform result and a saliency map corresponding to the infrared image; wherein the stationary wavelet transform result comprises: a low-pass sub-band and a detail sub-band.
In this embodiment, the infrared image I is subjected to stationary wavelet transform to obtain a low-pass sub-band L and a detail sub-band Hh、HvAnd Hd(ii) a And performing significance calculation on the infrared image I to obtain a significance map S of the infrared image I.
Specifically, first, a two-dimensional fourier transform is performed on the infrared image I, and a natural logarithm is taken. The specific calculation formula is as follows:
Figure BDA0002118647060000041
P(f)=φ(F(I))
L(f)=log(A(f))
wherein F represents a two-dimensional Fourier transform,
Figure BDA0002118647060000042
and phi represents amplitude operation and phase operation, respectively, p (f) and l (f) represent phase spectrum and log spectrum, respectively, a (f) represents the phase of I after fourier transform, and f (I) represents the fourier transform of I.
Further, a 3 × 3 mean filter h is applied to the amplitude part of the spectrumn(f) Filtering is carried out, and the calculation formula is as follows:
V(f)=L(f)*hn(f)
wherein: v (f) represents the result of the amplitude mean filtering of the spectrum, hn(f) Mean filter is represented, convolution operation is represented.
Further, the spectral residual r (f) is calculated as follows:
R(f)=L(f)-V(f)
further similarly, the significance map is obtained by taking an index of the log spectrum residual error and performing inverse Fourier transform, and the calculation formula is as follows:
Figure BDA0002118647060000051
wherein:
Figure BDA0002118647060000052
showing the significance map obtained by the inverse transformation, F-1Representing an inverse fourier transform.
However, experiments have found that obtaining a saliency map of an image is not directly useful for enhancement. Because the calculated saliency map is either too bright or too dark, the dynamic range is small. Therefore, the segmented binary image can be obtained as required through corresponding threshold processing. Since the binary significance map is abrupt at the boundary. Therefore, the invention makes it smooth enough by low-pass filtering it many times, and finally enhances the smooth result.
According to the embodiment of the invention, the logarithm of the amplitude part of the image subjected to the discrete Fourier transform is taken, and the initial saliency map is calculated.
Figure BDA0002118647060000053
Then, the initial saliency map is compared
Figure BDA0002118647060000054
Threshold segmentation is performed to make it a binary image, which is then smoothed using an averaging filter to the final directly usable saliency map S.
Figure BDA0002118647060000055
Wherein G is3×3Mean filter, G, representing a filter kernel size of 3 x 3meanWhich is indicative of the mean value filter,
Figure BDA0002118647060000056
graph representing significance
Figure BDA0002118647060000057
The binary segmentation result represents a convolution operation.
Step 2: and carrying out contrast increasing processing on the low-pass sub-band and carrying out linear enhancement processing on the detail sub-band to obtain an enhanced low-pass sub-band and an enhanced detail sub-band.
In this embodiment, CLAHE contrast enhancement processing is performed on the low-pass sub-band, and the processing algorithm is as follows:
L′=(1-α)L+α·CLAHE(L)
and step 3: and carrying out low-pass filtering and threshold segmentation processing on the saliency map to obtain a processed saliency map.
The steps 2 and 3 are not separated, and can be carried out simultaneously.
And 4, step 4: and performing dot multiplication on the enhanced detail sub-band and the processed saliency map to obtain a fused detail sub-band.
In this embodiment, the processing algorithm for performing linear enhancement processing on the detail subband is as follows:
Figure BDA0002118647060000058
Figure BDA0002118647060000059
Figure BDA0002118647060000061
v=2.5s/t·M
and 5: and performing stationary wavelet inverse transformation based on the enhanced low-pass sub-band and the fused detail sub-band to obtain an enhanced infrared image.
In this embodiment, the enhanced detail is sub-band H'h、H'v、H'dPerforming inverse stationary wavelet transform with the enhanced low-pass sub-band L' to obtain an enhanced infrared image
Figure BDA0002118647060000062
In the embodiment, a stationary wavelet transform result and a saliency map corresponding to an infrared image are obtained; wherein the stationary wavelet transform result comprises: a low-pass sub-band and a detail sub-band; carrying out contrast increasing processing on the low-pass sub-band and carrying out linear enhancement processing on the detail sub-band to obtain an enhanced low-pass sub-band and an enhanced detail sub-band; performing low-pass filtering and threshold segmentation processing on the saliency map to obtain a processed saliency map; performing dot multiplication on the enhanced detail sub-bands and the processed saliency map to obtain fused detail sub-bands; and performing stationary wavelet inverse transformation based on the enhanced low-pass sub-band and the fused detail sub-band to obtain an enhanced infrared image. The embodiment utilizes the stationary wavelet transform of the saliency map, can carry out multi-scale and multi-direction decomposition on the image, well represents the detail information of the image, and is beneficial to the detail enhancement of the infrared image; in addition, the saliency map can be very robust to infrared images containing noise; the efficiency is higher, the real-time is good.
Simulation experiments prove that the method has clear details, higher contrast, better visual effect and better objective evaluation index for the infrared image enhancement problem containing noise, and the method has good real-time performance and is an effective and feasible infrared image enhancement method.
Specifically, as can be seen from the enhancement result graphs shown in fig. 2(b) to 2(h) and fig. 3(b) to 3(h), the enhanced image of the present invention not only has improved overall contrast, but also has very clear edges and surface textures of the target object.
In addition, in order to better illustrate the superiority and advancement of the present invention, the objective quality of the enhancement results obtained by the technique of the present invention and those obtained by other methods were evaluated by using the commonly used 6 typical objective evaluation indexes for image enhancement. The 6 evaluation indexes are respectively as follows: GI-F, SSEQ, ES, ACG, EME, and AVC, and in addition to SSEQ, higher values of these evaluation indexes indicate better enhanced image quality. The average value of the objective evaluation index of the two sets of experimental images is shown in table 1.
TABLE 1
Figure BDA0002118647060000063
Figure BDA0002118647060000071
As can be seen from Table 1, the 6 objective evaluation indexes obtained by the enhancement result of the invention are superior to those of other methods, so that the invention can effectively improve the definition and detail information of the image.
In conclusion, the enhanced image obtained by the stationary wavelet transform infrared image enhancement method based on the saliency map has the advantages of good visual effect, rich detail information, high contrast and high efficiency.
The invention is not limited to the examples, and any equivalent changes to the technical solution of the invention by a person skilled in the art after reading the description of the invention are covered by the claims of the invention.
The embodiment of the present invention further provides an infrared image enhancement system, including: the device comprises a memory and a processor, wherein the memory stores executable instructions of the processor; wherein the processor is configured to perform the above-described infrared image enhancement method via execution of the executable instructions.
Optionally, a memory for storing a program; a Memory, which may include a volatile Memory (RAM), such as a Random Access Memory (SRAM), a Double Data Rate Synchronous Dynamic Random Access Memory (DDR SDRAM), and the like; the memory may also comprise a non-volatile memory, such as a flash memory. The memory 62 is used to store computer programs (e.g., applications, functional modules, etc. that implement the above-described methods), computer instructions, etc., which may be stored in one or more memories in a partitioned manner. And the computer programs, computer instructions, data, etc. described above may be invoked by a processor.
The computer programs, computer instructions, etc. described above may be stored in one or more memories in a partitioned manner. And the computer programs, computer instructions, data, etc. described above may be invoked by a processor.
A processor for executing the computer program stored in the memory to implement the steps of the method according to the above embodiments. Reference may be made in particular to the description relating to the preceding method embodiment.
The processor and the memory may be separate structures or may be an integrated structure integrated together. When the processor and the memory are separate structures, the memory, the processor may be coupled by a bus.
In addition, embodiments of the present application further provide a computer-readable storage medium, in which computer-executable instructions are stored, and when at least one processor of the user equipment executes the computer-executable instructions, the user equipment performs the above-mentioned various possible methods.
Computer-readable media includes both computer storage media and communication media including any medium that facilitates transfer of a computer program from one place to another. A storage media may be any available media that can be accessed by a general purpose or special purpose computer. An exemplary storage medium is coupled to the processor such the processor can read information from, and write information to, the storage medium. Of course, the storage medium may also be integral to the processor. The processor and the storage medium may reside in an ASIC. Additionally, the ASIC may reside in user equipment. Of course, the processor and the storage medium may reside as discrete components in a communication device.
The present application further provides a program product comprising a computer program stored in a readable storage medium, from which the computer program can be read by at least one processor of a server, the execution of the computer program by the at least one processor causing the server to carry out the method of any of the embodiments of the invention described above.
Those of ordinary skill in the art will understand that: all or a portion of the steps of implementing the above-described method embodiments may be performed by hardware associated with program instructions. The program may be stored in a computer-readable storage medium. When executed, the program performs steps comprising the method embodiments described above; and the aforementioned storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.

Claims (4)

1. An infrared image enhancement method, comprising:
s1: acquiring a stationary wavelet transform result and a saliency map corresponding to the infrared image; wherein the stationary wavelet transform result comprises: a low-pass sub-band and a detail sub-band;
performing stationary wavelet transform on the infrared image I to obtain a low-pass sub-band L and a detail sub-band Hh、HvAnd Hd(ii) a Carrying out significance calculation on the infrared image I to obtain a significance map S of the infrared image I;
specifically, first, two-dimensional fourier transform is performed on the infrared image I, and a natural logarithm is taken, and a specific calculation formula is as follows:
Figure FDA0003532170280000011
P(f)=φ(F(I))
L(f)=log(A(f))
wherein F represents a two-dimensional Fourier transform,
Figure FDA0003532170280000012
phi represents amplitude operation and phase operation respectively, P (f) and L (f) represent phase spectrum and logarithm spectrum respectively, A (f) represents amplitude after Fourier transformation to I, F (I) represents Fourier transformation to I;
using a 3 × 3 mean filter h for the log spectrum L (f)n(f) Filtering is carried out, and the calculation formula is as follows:
V(f)=L(f)*hn(f)
wherein: v (f) represents the result of the amplitude mean filtering of the spectrum, hn(f) Representing an averaging filter, representing a convolution operation;
calculating a spectrum residual R (f) according to the following formula:
R(f)=L(f)-V(f)
logarithm is taken on the amplitude part of the image after discrete Fourier transform, and an initial saliency map is calculated:
Figure FDA0003532170280000013
then, the initial saliency map is compared
Figure FDA0003532170280000014
Performing threshold segmentation to make the image into a binary image, and then smoothing the binary image to a finally directly usable saliency map S by using an average filter:
Figure FDA0003532170280000015
wherein G is3×3Mean filter, G, representing a filter kernel size of 3 x 3meanWhich is indicative of the mean value filter,
Figure FDA0003532170280000016
graph representing significance
Figure FDA0003532170280000017
Binary segmentation results, representing convolution operations;
s2: carrying out contrast increasing processing on the low-pass sub-band and carrying out linear enhancement processing on the detail sub-band to obtain an enhanced low-pass sub-band and an enhanced detail sub-band;
s3: performing low-pass filtering and threshold segmentation processing on the saliency map to obtain a processed saliency map;
s4: performing dot multiplication on the enhanced detail sub-bands and the processed saliency map to obtain fused detail sub-bands;
s5: and performing stationary wavelet inverse transformation based on the enhanced low-pass sub-band and the fused detail sub-band to obtain an enhanced infrared image.
2. The infrared image enhancement method according to claim 1, wherein the performing contrast increasing processing on the low-pass sub-band in S2 includes:
performing CLAHE contrast enhancement processing on the low-pass sub-band, wherein the processing algorithm is as follows:
L′=(1-α)L+α·CLAHE(L)
wherein: l denotes a low-pass subband, L' denotes an enhanced low-pass subband, α denotes a control strength, clahe (L) denotes a limit contrast adaptive histogram equalization, and · denotes a multiplication operation.
3. The infrared image enhancement method according to claim 1, wherein the processing algorithm for performing linear enhancement processing on the detail sub-bands in S2 is as follows:
Figure FDA0003532170280000021
Figure FDA0003532170280000022
Figure FDA0003532170280000023
v=2.5s/t·M
wherein: s represents the coefficient magnitude in the transform domain, M represents the maximum coefficient magnitude, and t represents the first enhancement parameter; hhRepresenting horizontal direction detail sub-bands, sign (·) representing sign function operation, b representing a second enhancement parameter, v representing an intermediate parameter, tanh (·) representing hyperbolic tangent function operation, c representing a third enhancement parameter, exp (·) representing exponential function operation,
Figure FDA0003532170280000024
representing the enhanced horizontal subband, HvShowing the sub-bands of detail in the vertical direction,
Figure FDA0003532170280000025
representing the enhanced vertical detail subband, HdThe detail sub-bands in the diagonal direction are shown,
Figure FDA0003532170280000026
representing the enhanced diagonal detail sub-bands.
4. An infrared image enhancement system, comprising: the device comprises a memory and a processor, wherein the memory stores executable instructions of the processor; wherein the processor is configured to perform the infrared image enhancement method of any of claims 1-3 via execution of the executable instructions.
CN201910601556.9A 2019-07-04 2019-07-04 Infrared image enhancement method and system Active CN110298806B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910601556.9A CN110298806B (en) 2019-07-04 2019-07-04 Infrared image enhancement method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910601556.9A CN110298806B (en) 2019-07-04 2019-07-04 Infrared image enhancement method and system

Publications (2)

Publication Number Publication Date
CN110298806A CN110298806A (en) 2019-10-01
CN110298806B true CN110298806B (en) 2022-04-12

Family

ID=68030473

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910601556.9A Active CN110298806B (en) 2019-07-04 2019-07-04 Infrared image enhancement method and system

Country Status (1)

Country Link
CN (1) CN110298806B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111709898B (en) * 2020-06-20 2023-05-23 昆明物理研究所 Infrared image enhancement method and system based on optimized CLAHE
CN112308114A (en) * 2020-09-24 2021-02-02 赣州好朋友科技有限公司 Method and device for sorting scheelite and readable storage medium
CN116580290B (en) * 2023-07-11 2023-10-20 成都庆龙航空科技有限公司 Unmanned aerial vehicle identification method, unmanned aerial vehicle identification device and storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105184752A (en) * 2015-09-23 2015-12-23 成都融创智谷科技有限公司 Image processing method based on wavelet transform
CN108389158A (en) * 2018-02-12 2018-08-10 河北大学 A kind of infrared and visible light image interfusion method
CN109191417A (en) * 2018-09-11 2019-01-11 中国科学院长春光学精密机械与物理研究所 It is detected based on conspicuousness and improves twin-channel method for self-adaption amalgamation and device
US10249032B2 (en) * 2010-04-23 2019-04-02 Flir Systems Ab Infrared resolution and contrast enhancement with fusion

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2431919A1 (en) * 2010-09-16 2012-03-21 Thomson Licensing Method and device of determining a saliency map for an image
CN104050638A (en) * 2014-06-12 2014-09-17 杭州电子科技大学 Saliency method infrared small target enhancing method combined with scale optimization
CN107862666A (en) * 2017-11-22 2018-03-30 新疆大学 Mixing Enhancement Methods about Satellite Images based on NSST domains

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10249032B2 (en) * 2010-04-23 2019-04-02 Flir Systems Ab Infrared resolution and contrast enhancement with fusion
CN105184752A (en) * 2015-09-23 2015-12-23 成都融创智谷科技有限公司 Image processing method based on wavelet transform
CN108389158A (en) * 2018-02-12 2018-08-10 河北大学 A kind of infrared and visible light image interfusion method
CN109191417A (en) * 2018-09-11 2019-01-11 中国科学院长春光学精密机械与物理研究所 It is detected based on conspicuousness and improves twin-channel method for self-adaption amalgamation and device

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Image Contrast Enhancement by Contourlet Transform;Ehsan Nezhadarya 等;《48th International Symposium ELMAR-2006》;20060609;81-84 *
Infrared image enhancement through saliency feature analysis based on multi-scale decomposition;Jufeng Zhao 等;《Infrared Physics & Technology》;20131121;第62卷;86-93 *
Infrared super-resolution imaging using multi-scale saliency and deep wavelet residuals;Gunnam Suryanarayana 等;《Infrared Physics & Technology》;20181231;第97卷;177-185 *
一种基于平稳小波域的红外图像增强方法;龚昌来;《激光与红外》;20130620;第43卷(第6期);703-707 *
基于视觉显著性与对比度增强的红外图像融合;张承泓 等;《红外技术》;20170525;第39卷(第5期);421-426 *
小波变换和自适应变换相结合的图像增强方法;于天河;《哈尔滨理工大学学报》;20181231;第23卷(第6期);100-104 *

Also Published As

Publication number Publication date
CN110298806A (en) 2019-10-01

Similar Documents

Publication Publication Date Title
Chen et al. Gaussian-adaptive bilateral filter
Jain et al. A survey of edge-preserving image denoising methods
Ahn et al. Block-matching convolutional neural network for image denoising
CN114820352B (en) Hyperspectral image denoising method, device and storage medium
Bartyzel Adaptive kuwahara filter
CN105046664B (en) A kind of image de-noising method based on adaptive EPLL algorithms
CN110298806B (en) Infrared image enhancement method and system
CN112508810A (en) Non-local mean blind image denoising method, system and device
US9569684B2 (en) Image enhancement using self-examples and external examples
Dai et al. Entropy-based bilateral filtering with a new range kernel
US9443286B2 (en) Gray image processing method and apparatus based on wavelet transformation
Kamble et al. Performance evaluation of wavelet, ridgelet, curvelet and contourlet transforms based techniques for digital image denoising
CN105740876A (en) Image preprocessing method and device
CN102542542A (en) Image denoising method based on non-local sparse model
CN104574293A (en) Multiscale Retinex image sharpening algorithm based on bounded operation
CN114155161B (en) Image denoising method, device, electronic equipment and storage medium
CN104657951A (en) Multiplicative noise removal method for image
CN106815818A (en) A kind of image de-noising method
CN111192204A (en) Image enhancement method, system and computer readable storage medium
Shahdoosti et al. Combined ripplet and total variation image denoising methods using twin support vector machines
CN104616259A (en) Non-local mean image de-noising method with noise intensity self-adaptation function
Zhao et al. Constant time texture filtering
Walha et al. Handling noise in textual image resolution enhancement using online and offline learned dictionaries
CN105303538A (en) Gauss noise variance estimation method based on NSCT and PCA
CN102339460B (en) Adaptive satellite image restoration method

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