CN112258430A - Universal correction method for remote sensing image radiation nonuniformity - Google Patents
Universal correction method for remote sensing image radiation nonuniformity Download PDFInfo
- Publication number
- CN112258430A CN112258430A CN202011187801.5A CN202011187801A CN112258430A CN 112258430 A CN112258430 A CN 112258430A CN 202011187801 A CN202011187801 A CN 202011187801A CN 112258430 A CN112258430 A CN 112258430A
- Authority
- CN
- China
- Prior art keywords
- image
- nsct
- low
- band
- radiation
- 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.)
- Granted
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000012937 correction Methods 0.000 title claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 12
- 238000001914 filtration Methods 0.000 claims abstract description 10
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 238000005286 illumination Methods 0.000 abstract description 3
- 238000003384 imaging method Methods 0.000 abstract description 3
- 230000009466 transformation Effects 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration by the use of local operators
-
- G06T5/80—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/194—Segmentation; Edge detection involving foreground-background segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/25—Determination of region of interest [ROI] or a volume of interest [VOI]
-
- 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/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- 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/20024—Filtering details
Abstract
A universal correction method for remote sensing image radiation unevenness relates to the technical field of image processing, and solves the problem of radiation unevenness caused by factors such as illumination conditions, atmospheric conditions, sensor equipment, imaging conditions and the like of the existing remote sensing image; step two, carrying out NSCT on the obtained low-frequency subband imageLowpassPerforming Gaussian filtering to obtain an image background image, and removing the image background image; performing linear restoration enhancement processing on the image radiance value based on the local interesting area radiance information; step four, replacing the low-pass sub-band image, and performing NSCT reconstruction to obtain a corrected wave band shadowAn image; the invention can well inhibit the phenomenon of nonuniform radiation in the image by combining non-subsampled contourlet transformation, Gaussian filtering and a linear enhancement method based on local characteristics, and simultaneously, the change of radiation correction on the real radiation information of the ground object is greatly reduced.
Description
Technical Field
The invention relates to the technical field of image processing, in particular to a universal correction method for a remote sensing image with non-uniform radiation.
Background
At present, most of correction methods for remote sensing image radiation unevenness perform corresponding compensation correction based on a distribution mode of radiance, and have certain applicability. However, in practical situations, the radiation unevenness of the image has diversity, and there is a phenomenon that the radiation brightness changes irregularly and irregularly, and it is difficult to fit with an appropriate mathematical model.
In summary, the conventional radiation correction method for the radiation nonuniformity of the remote sensing image generally has certain limitations, and is difficult to adapt to the radiation nonuniformity in reality.
Disclosure of Invention
The invention provides a universal correction method for remote sensing image radiation nonuniformity, aiming at solving the problem of nonuniform radiation caused by factors such as illumination conditions, atmospheric conditions, sensor equipment, imaging conditions and the like of the existing remote sensing image.
A universal correction method for remote sensing image radiation nonuniformity is realized by the following steps:
firstly, carrying out multi-scale decomposition on an original image by adopting an NSCT (non-subsampled Contourlet transform) method;
dividing an original image according to wave bands, and performing multi-scale decomposition on each wave band of the original image by adopting an NSCT (non-subsampled Contourlet transform) method; obtaining a low-frequency sub-band image NSCT of each wave bandLowpassAnd a plurality of directional band-pass sub-band images NSCTBandpassSte；
Step two, performing NSCT on the low-frequency subband image obtained in the step oneLowpassPerforming Gaussian filtering to obtain an image background image, and removing the image background image; the method specifically comprises the following steps:
NSCT of the low-frequency subband image by adopting Gaussian filteringLowpassProcessing to obtain non-uniform background image InfoBackgroundAnd the non-uniform background image information InfoBackgroundEliminating to obtain ground object radiation low-frequency image InfoFeatureLow；InfoFeatureLow＝NSCTLowpass-InfoBackground；
Performing linear restoration enhancement processing on the image radiance value based on the local interesting area radiance information; the method specifically comprises the following steps:
the ground object radiation low-frequency image Info obtained in the step twoFeatureLowSelecting a local region of interest, and fitting a low-frequency sub-band image NSCT in the local region of interest by a least square methodLowpassAnd ground object radiation low-frequency image InfoFeatureLowAnd the linear coefficient obtained by fitting is used for radiating the low-frequency image Info of the ground objectFeatureLowCorrecting to obtain corrected low-frequency image InfoFeatureLowCor；
Replacing the low-pass sub-band image, and performing NSCT reconstruction to obtain a corrected band image;
adopting the corrected low-frequency image Info obtained in the third stepFeatureLowCorReplacement of low frequency subband images NSCTLowpassCombining the multi-direction band-pass sub-band image NSCT of the step oneBandpassSteNSCT reconstruction is carried out to obtain each wave band of the corrected image, and each wave band is combined to realize general correction of the original image.
The invention has the beneficial effects that: the general correction method can well inhibit the phenomenon of radiation nonuniformity in the image and simultaneously reduce the change of radiation correction on the real radiation information of the ground object to a greater extent by combining non-subsampled Contourlet (NSCT), Gaussian filtering and a local feature-based linear enhancement method. In addition, the method is not required to distinguish the type of the radiation unevenness of the image, and is applicable to images with various radiation unevenness conditions.
The universal correction method can effectively inhibit the phenomenon of uneven image radiation, simultaneously inhibits the change of radiation correction to the real radiation information of the ground object, has good universality to various types of uneven radiation, and provides a universal solution for the problem of uneven image radiation.
Drawings
FIG. 1 is a flow chart of a general correction method for remote sensing image radiation non-uniformity according to the present invention;
FIG. 2 is a graph of the luminance effect of an original remote sensing image;
fig. 3 is a luminance effect diagram after radiation homogenization correction.
Detailed Description
The invention is further illustrated with reference to the accompanying drawings and specific examples.
The main process of the invention comprises four parts of NSCT multi-scale decomposition, background information estimation and elimination by means of Gaussian filtering, ROI characteristic-based linear recovery enhancement and NSCT image reconstruction, and the detailed process is shown in figure 1.
Firstly, carrying out multi-scale decomposition on an original image wave band 1 by adopting NSCT;
NSCT is a two-dimensional image representation method combining non-downsampling pyramid decomposition and a non-downsampling direction filter bank, has the characteristics of multi-resolution, locality, multi-directionality, anisotropy, translation invariance and the like, and can effectively decompose an original two-dimensional image into a series of multi-scale band-pass sub-band images and a low-pass sub-band image. Here, NSCT is used to perform multi-scale decomposition on the original image band 1 to obtain 1 low-frequency subband image and several directional band-pass subband images, which are respectively denoted as NSCTLowpassAnd NSCTBandpassSte。
Secondly, based on the low-pass sub-band image obtained in the previous step, obtaining an image background image by means of Gaussian filtering, and removing the low-frequency sub-band image NSCT obtained in the previous stepLowpassThe method comprises the background information of the image and part of the low-frequency radiation information of the ground objects, wherein the image is processed by adopting Gaussian low-pass filtering, and then the non-uniform background image Info in the original image can be extractedBackgroundRemoving the uneven background image from the low-frequency sub-band image, and leaving the ground object radiation low-frequency image InfoFeatureLowThe following formula:
InfoFeatureLow＝NSCTLowpass-InfoBackground
and thirdly, based on a local Region of Interest (ROI) radiance image, after removing the background image in the previous step of linear restoration and enhancement of the radiance value of the image, the radiance information of the image is lost, so that linear restoration and enhancement processing of radiance needs to be performed on the low-frequency image of the ground feature obtained after processing. Because the radiation of the whole image is not uniform, an ROI with relatively uniform radiation and real color is selected in the image, and NSCT in the ROI range is fitted through a least square methodLowpassAnd InfoFeatureLowAnd the linear coefficient obtained by fitting is used for the whole InfoFeatureLowCorrecting to obtain InfoFeatureLowCor。
Fourthly, replacing the low-pass sub-band image, carrying out NSCT reconstruction to obtain a corrected band image, and converting the Info into a new entryFeatureLowCorReplacement NSCTLowpassCombined NSCTBandpassSteNSCT reconstruction is performed to obtain an image band 1 in which the radiation non-uniformity phenomenon is suppressed by correction.
And (3) performing similar processing on other residual wave bands, repeating the steps 1-4 on the wave band 2 and the wave band 3 of the original image, and finally combining the three wave bands obtained by correction processing to finish the processing of the original remote sensing image.
The image used in the embodiment is a satellite image in a certain area of a scientific and technological development area in north lake of Changchun city, Jilin province, the coverage range is 125.377-125.432 degrees from east longitude, 43.972-43.992 degrees from north latitude, the size of the image is 3584 pixels 10240 pixels, and the image has an obvious radiation non-uniform phenomenon due to the influence of one or more reasons such as illumination conditions, atmospheric conditions, sensor equipment, imaging conditions and the like, and a high brightness band exists in the middle area of the image. Fig. 2 is a luminance graph of an original remote sensing image (after stretching with standard deviation of n equal to 1), wherein the upper left black thick frame selected region is the ROI region of the present embodiment, and fig. 3 is a luminance graph after radiation correction by this method (after stretching with standard deviation of n equal to 1). As can be seen from the observation, the present embodiment can effectively suppress the phenomenon of radiation unevenness in the image.
The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and it is not intended that the invention be limited to these specific details. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.
Claims (3)
1. A universal correction method for remote sensing image radiation nonuniformity is characterized in that: the method is realized by the following steps:
firstly, carrying out multi-scale decomposition on an original image by adopting an NSCT (non-subsampled Contourlet transform) method;
dividing an original image according to wave bands, and performing multi-scale decomposition on each wave band of the original image by adopting an NSCT (non-subsampled Contourlet transform) method; obtaining a low-frequency sub-band image NSCT of each wave bandLowpassAnd a plurality of directional band-pass sub-band images NSCTBandpassSte；
Step two, performing NSCT on the low-frequency subband image obtained in the step oneLowpassPerforming Gaussian filtering to obtain an image background image, and removing the image background image;
thirdly, based on the local interesting area radiance information, carrying out linear restoration enhancement processing on the image radiance value to obtain corrected low-frequency image InfoFeatureLowCor；
Replacing the low-pass sub-band image, and performing NSCT reconstruction to obtain a corrected band image;
adopting the corrected low-frequency image Info obtained in the third stepFeatureLowCorReplacement of low frequency subband images NSCTLowpassCombining the multi-direction band-pass sub-band image NSCT of the step oneBandpassSteNSCT reconstruction is carried out to obtain each wave band of the corrected image, and each wave band is combined to realize general correction of the original image.
2. The universal correction method for remote sensing image radiation nonuniformity according to claim 1, characterized in that: the specific process of the second step is as follows:
NSCT of the low-frequency subband image by adopting Gaussian filteringLowpassProcessing to obtain non-uniform background image InfoBackgroundAnd the non-uniform background image information InfoBackgroundEliminating to obtain ground object radiation low-frequency image InfoFeatureLow；InfoFeatureLow＝NSCTLowpass-InfoBackground。
3. The universal correction method for remote sensing image radiation nonuniformity according to claim 1, characterized in that: the concrete process of the third step is as follows:
the ground object radiation low-frequency image Info obtained in the step twoFeatureLowSelecting a local region of interest, and fitting a low-frequency sub-band image NSCT in the local region of interest by a least square methodLowpassAnd ground object radiation low-frequency image InfoFeatureLowAnd the linear coefficient obtained by fitting is used for radiating the low-frequency image Info of the ground objectFeatureLowCorrecting to obtain corrected low-frequency image InfoFeatureLowCor。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011187801.5A CN112258430B (en) | 2020-10-30 | 2020-10-30 | Universal correction method for remote sensing image radiation nonuniformity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011187801.5A CN112258430B (en) | 2020-10-30 | 2020-10-30 | Universal correction method for remote sensing image radiation nonuniformity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112258430A true CN112258430A (en) | 2021-01-22 |
| CN112258430B CN112258430B (en) | 2022-06-21 |
Family
ID=74269035
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011187801.5A Active CN112258430B (en) | 2020-10-30 | 2020-10-30 | Universal correction method for remote sensing image radiation nonuniformity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112258430B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6108454A (en) * | 1998-04-27 | 2000-08-22 | The United States Of America As Represented By The Secretary Of The Navy | Line contrast difference effect correction for laser line scan data |
| CN102222322A (en) * | 2011-06-02 | 2011-10-19 | 西安电子科技大学 | Multiscale non-local mean-based method for inhibiting infrared image backgrounds |
| CN103226819A (en) * | 2013-03-29 | 2013-07-31 | 武汉大学 | Segmental counting-based relative radiation correction method |
| US20150078651A1 (en) * | 2011-04-14 | 2015-03-19 | Inb Vision Ag | Device and method for measuring surfaces |
| CN106295696A (en) * | 2016-08-09 | 2017-01-04 | 中国科学院遥感与数字地球研究所 | A kind of multi-source Remote Sensing Images radiation normalization method |
| CN107274365A (en) * | 2017-06-15 | 2017-10-20 | 中国矿业大学(北京) | A kind of mine image intensification method based on unsharp masking and NSCT algorithms |
| CN108830814A (en) * | 2018-06-15 | 2018-11-16 | 武汉大学 | A kind of relative radiometric correction method of remote sensing image |
| CN111343439A (en) * | 2018-12-19 | 2020-06-26 | 罗伯特·博世有限公司 | Method and apparatus for processing image signal of image sensor of vehicle |
-
2020
- 2020-10-30 CN CN202011187801.5A patent/CN112258430B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6108454A (en) * | 1998-04-27 | 2000-08-22 | The United States Of America As Represented By The Secretary Of The Navy | Line contrast difference effect correction for laser line scan data |
| US20150078651A1 (en) * | 2011-04-14 | 2015-03-19 | Inb Vision Ag | Device and method for measuring surfaces |
| CN102222322A (en) * | 2011-06-02 | 2011-10-19 | 西安电子科技大学 | Multiscale non-local mean-based method for inhibiting infrared image backgrounds |
| CN103226819A (en) * | 2013-03-29 | 2013-07-31 | 武汉大学 | Segmental counting-based relative radiation correction method |
| CN106295696A (en) * | 2016-08-09 | 2017-01-04 | 中国科学院遥感与数字地球研究所 | A kind of multi-source Remote Sensing Images radiation normalization method |
| CN107274365A (en) * | 2017-06-15 | 2017-10-20 | 中国矿业大学(北京) | A kind of mine image intensification method based on unsharp masking and NSCT algorithms |
| CN108830814A (en) * | 2018-06-15 | 2018-11-16 | 武汉大学 | A kind of relative radiometric correction method of remote sensing image |
| CN111343439A (en) * | 2018-12-19 | 2020-06-26 | 罗伯特·博世有限公司 | Method and apparatus for processing image signal of image sensor of vehicle |
Non-Patent Citations (3)
| Title |
|---|
| YONGXING JIA .ETC: "Research on the decomposition and fusion method for the infrared and visible images based on the guided image filtering and Gaussian filter", 《2017 3RD IEEE INTERNATIONAL CONFERENCE ON COMPUTER AND COMMUNICATIONS (ICCC)》 * |
| 李成范等: "《火山灰云遥感监测方法与应用》", 31 August 2016 * |
| 黄黎红: "单幅图像的快速去雾算法", 《光电子.激光》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112258430B (en) | 2022-06-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Singh et al. | Super-resolving noisy images | |
| CN108765347B (en) | Color enhancement method suitable for remote sensing image | |
| CN106934771B (en) | Infrared image stripe noise removing method based on local correlation | |
| Gupta et al. | Review of different local and global contrast enhancement techniques for a digital image | |
| CN111583123A (en) | Wavelet transform-based image enhancement algorithm for fusing high-frequency and low-frequency information | |
| Suganya et al. | Survey on image enhancement techniques | |
| Tang et al. | Single image dehazing via lightweight multi-scale networks | |
| DE102018119625A1 (en) | Reduction of structured IR patterns in stereoscopic depth sensor imaging | |
| CN110533600B (en) | Same/heterogeneous remote sensing image high-fidelity generalized space-spectrum fusion method | |
| CN104021532A (en) | Image detail enhancement method of infrared image | |
| CN109801233B (en) | Method for enhancing true color remote sensing image | |
| CN103268596A (en) | Method for reducing image noise and enabling colors to be close to standard | |
| CN108961325A (en) | Method for registering between more/high-spectrum remote sensing wave band | |
| Zhu et al. | Underwater image enhancement based on colour correction and fusion | |
| CN108537744A (en) | A kind of coloured image luminance component homomorphic filtering defogging method | |
| CN106981052B (en) | Adaptive uneven brightness variation correction method based on variation frame | |
| CN101889295B (en) | Method of generating a multiscale contrast enhanced image | |
| CN112258430B (en) | Universal correction method for remote sensing image radiation nonuniformity | |
| CN111553960B (en) | Ring artifact rapid correction method based on projection mean image | |
| CN110415188A (en) | A kind of HDR image tone mapping method based on Multiscale Morphological | |
| CN106780400B (en) | Image processing method and device | |
| CN107644226A (en) | Be advantageous to the image processing method of image recognition | |
| CN109886901B (en) | Night image enhancement method based on multi-channel decomposition | |
| CN115937302A (en) | Hyperspectral image sub-pixel positioning method combined with edge preservation | |
| CN115588060A (en) | Projection average image-based rapid CT (computed tomography) ring artifact correction method and system |
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 | ||
| CB02 | Change of applicant information |
Address after: No. 1299, Mingxi Road, Beihu science and Technology Development Zone, Changchun City, Jilin Province Applicant after: Changguang Satellite Technology Co.,Ltd. Address before: 130032 no.1299, Mingxi Road, Beihu science and Technology Development Zone, Changchun City, Jilin Province Applicant before: CHANG GUANG SATELLITE TECHNOLOGY Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A universal correction method for radiation non-uniformity in remote sensing images Effective date of registration: 20230724 Granted publication date: 20220621 Pledgee: Industrial and Commercial Bank of China Limited Changchun Southern Urban Economic Development Zone Branch Pledgor: Changguang Satellite Technology Co.,Ltd. Registration number: Y2023220000056 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right |