CN110310249B - Visual enhancement method for remote sensing image - Google Patents
Visual enhancement method for remote sensing image Download PDFInfo
- Publication number
- CN110310249B CN110310249B CN201910417005.7A CN201910417005A CN110310249B CN 110310249 B CN110310249 B CN 110310249B CN 201910417005 A CN201910417005 A CN 201910417005A CN 110310249 B CN110310249 B CN 110310249B
- Authority
- CN
- China
- Prior art keywords
- remote sensing
- image
- sensing image
- resolution
- low
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 230000000007 visual effect Effects 0.000 title claims abstract description 17
- 238000012800 visualization Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 239000011248 coating agent Substances 0.000 claims description 2
- 238000000576 coating method Methods 0.000 claims description 2
- 238000002360 preparation method Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 7
- 230000002708 enhancing effect Effects 0.000 description 4
- 230000001105 regulatory effect Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4053—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution
-
- 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 using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a visual enhancement method for a remote sensing image, which provides a low-resolution multiband remote sensing image; and a high resolution remote sensing image having the same scene and the same size as the low resolution multiband remote sensing image; and (3) utilizing the high-resolution remote sensing image to participate in adjusting the overall dynamic display range of the low-resolution multiband remote sensing image so as to enhance the spatial detail information of each waveband, and finally obtaining the visually enhanced multiband remote sensing image. The problem of current remote sensing image visual effect not good is solved.
Description
[ technical field ] A method for producing a semiconductor device
The invention relates to a method for enhancing a remote sensing image, in particular to a visual enhancement method for the remote sensing image.
[ background of the invention ]
Some satellite-borne/airborne sensors, such as WorldView-3, AVIRIS and the like, can provide multi-band/high-spectrum remote sensing images with multiple bands, and due to the restriction of the existing hardware level, the image data can not record the real information of the scene observed by human eyes, and the brightness, the contrast, the color and the like of the image data are limited. Therefore, it is necessary to adjust and enhance the color saturation and detail information of such images, so as to improve the visualization effect and provide a more realistic sensory experience for the viewer.
[ summary of the invention ]
The invention aims to provide a visualization enhancing method for a remote sensing image, which aims to solve the problem that the existing remote sensing image is poor in visualization effect.
The invention adopts the following technical scheme: a visualization enhancing method for remote sensing images provides,
a low resolution multiband remote sensing image; and a process for the preparation of a coating,
the high-resolution remote sensing image and the low-resolution multiband remote sensing image have the same scene and the same size;
and (3) utilizing the high-resolution remote sensing image to participate in adjusting the overall dynamic display range of the low-resolution multiband remote sensing image so as to enhance the spatial detail information of each waveband, and finally obtaining the visually enhanced multiband remote sensing image.
Further, the high-resolution remote sensing image is recorded as H, and the processing method of the high-resolution remote sensing image H comprises the following steps:
firstly, carrying out logarithmic transformation on the high-resolution remote sensing image H, and outputting and recording as HLOG，
HLOG＝LOG(H) (1)；
Second, image H is filtered using a low pass filterLOGFiltering to obtain a low-frequency image HL；
Thirdly, calculating a detail image HDThe following were used:
HD＝HLOG-HL(2)；
fourthly, adjusting the detail image H by using an adjusting coefficient C1DAdjusting the low-frequency image H by using the adjustment coefficient C2LThen obtaining an enhanced image H by calculation according to the formula (3)E，
HE＝C1×HD+C2×HL(3)；
The fifth step is to enhance the image HEInverse logarithmic transformation is carried out, and the obtained enhanced image is recorded as HN，
HN＝LOG-1(HE) (4)。
Further, it is assumed that the low-resolution multiband remote sensing image comprises N bands, which are marked as L1,L2,…LNThe processing method of the low-resolution multiband remote sensing image comprises the following steps:
first step, calculating low resolution multiband remote sensing image (L)1,…,LN) The intensity component of (a) is denoted as I;
second, enhance the image HNDivided by the intensity component I to obtain a ratio image R,
R＝HN/I (5)；
thirdly, a low-resolution multiband remote sensing image (L) with N wave bands1,…,LN) Respectively multiplied by the ratio images R to obtain a visual enhanced N-waveband remote sensing image (V)1,…,VN)
Vi＝R×Li,i＝1,...,N (6)。
The method has the advantages that the auxiliary high-resolution remote sensing image is introduced to participate in adjusting the overall dynamic display range of the low-resolution multiband remote sensing image, so that the spatial detail information of each waveband is enhanced; the method of the invention enhances each wave band of the low-resolution multiband remote sensing image, breaks through the enhanced display only in the red, green and blue wave band range in the prior art, can realize the visual enhancement of pseudo color, is an effective visual enhancement method suitable for the multiband remote sensing image, well displays the remote sensing image with high dynamic range on the common display with low dynamic range with richer detail information, and achieves the shocking and exciting effect.
[ description of the drawings ]
FIG. 1 is a flow chart of a visualization enhancing method for remote sensing images of the present invention.
[ detailed description ] embodiments
The technical solution of the present invention is further described in detail by the accompanying drawings and embodiments.
The invention provides a visual enhancement method for a remote sensing image, which provides two images, comprising the following steps:
the system comprises a low-resolution multiband remote sensing image and a high-resolution multiband remote sensing image, wherein the high-resolution multiband remote sensing image and the low-resolution multiband remote sensing image firstly have the same scene, and secondly are strictly registered to ensure that the high-resolution multiband remote sensing image and the low-resolution multiband remote sensing image have the same line number and column number, namely the image sizes are the same;
and then, the high-resolution remote sensing image is used for participating in adjusting the overall dynamic display range of the low-resolution multiband remote sensing image so as to enhance the spatial detail information of each waveband of the low-resolution multiband remote sensing image, and finally the multiband visual enhanced multiband remote sensing image is obtained.
The invention provides a visual enhancement method for a remote sensing image, which specifically comprises the following steps:
suppose that the low-resolution multiband remote sensing image comprises N wave bands, which are marked as L1,L2,…LNAnd the auxiliary image is a high-resolution remote sensing image with the same scene and is marked as H. Image L1,L2,…LNHas been registered exactly with image H and has the same number of rows and columns, i.e. the image size.
Firstly, carrying out logarithmic transformation on the high-resolution remote sensing image H, and outputting and recording as HLOGIs calculated as follows
HLOG＝LOG(H) (1)；
Second, image H is filtered using a low pass filterLOGFiltering to obtain a low-frequency image HL；
Thirdly, calculating a detail image HDAs follows
HD＝HLOG-HL(2)；
Fourthly, adjusting the detail image H by using an adjusting coefficient C1DAdjusting the low-frequency image H by using the adjustment coefficient C2LThen obtaining an enhanced image H by calculation according to the formula (3)E；
HE＝C1×HD+C2×HL(3)；
Wherein, the value ranges of the regulating coefficients C1 and C2 are respectively that C1 is more than or equal to 1.1 and less than or equal to 1.5 and C2 is more than or equal to 0.90 and less than or equal to 0.99;
the fifth step is to enhance the image HEInverse logarithmic transformation is carried out, and the obtained enhanced image is recorded as HN
HN＝LOG-1(HE) (4)；
Sixthly, calculating a low-resolution multiband remote sensing image (L) with N wave bands1,…,LN) The intensity component of (a) is denoted as I;
seventh step, enhance the image HNDividing the intensity component I to obtain a ratio image R
R＝HN/I (5)；
Eighthly, the N wave band low-resolution multiband remote sensing image (L)1,…,LN) Respectively multiplied by the ratio images R to obtain a visual enhanced N-waveband remote sensing image (V)1,…,VN)
Vi＝R×Li,i＝1,...,N (6)，
V1,…,VNI.e. the output.
Examples
The experimental data was selected from the low resolution eight band images and high resolution images provided by the WorldView-2 satellite. The spatial resolution of the low-resolution eight-waveband image is 2 meters, the spatial resolution of the high-resolution image is 0.5 meter, the image size is 1600 pixels × 1600 pixels, and the quantization bit number is 13 bits. The image scene is a partial area of the italian roman city, and the content comprises automobiles, buildings, roads, trees, sports fields and the like. It can be easily seen from the RGB composite image of the low-resolution eight-band image that the overall definition is not high and the colors show unnatural phenomena, especially in the shadow areas around buildings and trees and light-colored roof parts with high reflection characteristics.
Selecting a gaussian low-pass filter as the low-pass filter used in the second step, the window radius of the filter being selected to be 2 and the standard deviation being selected to be 3; the regulating coefficient C1 in the fourth step is 1.2, and C2 is 0.98. In the output obtained visualization enhancement result, the overall definition of the image is greatly improved, the outlines and edges of small ground objects such as roads, automobiles and the like are clearly displayed, the content in a shadow area can be clearly displayed, the highly-reflective white roof has layering, and the visualization effect of the image is obviously enhanced. Closing device
The value ranges of the C1 and C2 adjustment coefficients are respectively defined as 1.1-1.25C 1 and 0.90-0.24C 2 and 0.99, and the following four experimental results of the value ranges of the C1 and the C2 adjustment coefficients exceeding the suggested range are listed:
1. if C1 is 2 and C2 is 2, then:
both C1 and C2 are larger than the suggested value range, and in the output enhancement result, the white roof is displayed brighter, the shadow area is darker, the green vegetation is displayed darker, and the visual effect is poor;
2. if C1 is 0.5 and C2 is 0.5, then:
both C1 and C2 are smaller than the suggested value range, the overall color of the output image is displayed too light, and obvious spectral distortion phenomena appear, such as yellow green on green lawn, orange red on sports ground, and very unnatural display effect.
3. If C1 is 2 and C2 is 0.98, then:
only C1 is larger than the suggested range, the output image is excessively sharpened, resulting in a distortion of the spatial information, such as a white roof that has changed from original white to off-white and milky-white.
4. If C1 is 1.2 and C2 is 0.5, then:
only C2 is smaller than the suggested value range, the output image has obvious spectral distortion, the color distortion is concentrated on a green vegetation part and a white roof part, in addition, the image also has obvious noise points, and the visual effect is poor.
From the above analysis, if the adjustment coefficient is not selected within the setting range of the present invention, the effect is poor.
According to the invention, the auxiliary high-resolution remote sensing image is introduced to participate in adjusting the overall dynamic display range of the low-resolution multiband remote sensing image, so that the spatial detail information of each waveband is enhanced; the method of the invention enhances each wave band of the low-resolution multiband remote sensing image, breaks through the enhanced display only in the red, green and blue wave band range in the prior art, can realize the visual enhancement of pseudo color, is an effective visual enhancement method suitable for the multiband remote sensing image, well displays the remote sensing image with high dynamic range on the common display with low dynamic range with richer detail information, and achieves the shocking and exciting effect.
Claims (1)
1. A visualization enhancement method for remotely sensed images, characterized by providing,
a low resolution multiband remote sensing image; and a process for the preparation of a coating,
a high-resolution remote sensing image having the same scene and the same size as the low-resolution multiband remote sensing image;
utilizing the high-resolution remote sensing image to participate in adjusting the overall dynamic display range of the low-resolution multiband remote sensing image so as to enhance the spatial detail information of each waveband of the low-resolution multiband remote sensing image and finally obtain a visualized enhanced multiband remote sensing image;
recording the high-resolution remote sensing image as H, wherein the processing method of the high-resolution remote sensing image H comprises the following steps:
firstly, carrying out logarithmic transformation on the high-resolution remote sensing image H, and outputting and recording as HLOG，
HLOG＝LOG(H) (1)；
Second, image H is filtered using a low pass filterLOGFiltering to obtain a low-frequency image HL；
Thirdly, calculating a detail image HDThe following were used:
HD＝HLOG-HL(2)；
fourthly, adjusting the detail image H by using an adjusting coefficient C1DAdjusting the low-frequency image H by using the adjustment coefficient C2LThen obtaining an enhanced image H by calculation according to the formula (3)E，
HE＝C1×HD+C2×HL(3)；
The fifth step is to enhance the image HEInverse logarithmic transformation is carried out, and the obtained enhanced image is recorded as HN，
HN＝LOG-1(HE) (4)；
Sixthly, calculating a low-resolution multiband remote sensing image L with N wavebands1,…,LNThe intensity component of (a) is denoted as I;
seventh step, enhance the image HNDividing the intensity component I to obtain a ratio image R
R＝HN/I (5)；
Eighthly, the N-waveband low-resolution multiband remote sensing image L1,…,LNRespectively multiplied by the ratio images R to obtain a visual enhanced N-waveband remote sensing image V1,…,VN
Vi＝R×Li,i＝1,...,N (6)，
V1,…,VNI.e. the output.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910417005.7A CN110310249B (en) | 2019-05-20 | 2019-05-20 | Visual enhancement method for remote sensing image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910417005.7A CN110310249B (en) | 2019-05-20 | 2019-05-20 | Visual enhancement method for remote sensing image |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110310249A CN110310249A (en) | 2019-10-08 |
| CN110310249B true CN110310249B (en) | 2020-09-08 |
Family
ID=68074707
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910417005.7A Active CN110310249B (en) | 2019-05-20 | 2019-05-20 | Visual enhancement method for remote sensing image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110310249B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111242882A (en) * | 2020-01-08 | 2020-06-05 | 长春工程学院 | Hyperspectral remote sensing geological survey control system, method and application |
| CN113436205A (en) * | 2021-06-16 | 2021-09-24 | 中国电子科技集团公司第五十四研究所 | Remote sensing image rapid interpretation method based on sight tracking |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102521815A (en) * | 2011-11-02 | 2012-06-27 | 薛笑荣 | Fast fusion system and fast fusion method for images |
| CN103530860A (en) * | 2013-09-26 | 2014-01-22 | 天津大学 | Adaptive autoregressive model-based hyper-spectral imagery super-resolution method |
| CN107220957A (en) * | 2017-04-25 | 2017-09-29 | 西北工业大学 | It is a kind of to utilize the remote sensing image fusion method for rolling Steerable filter |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9020256B2 (en) * | 2011-04-22 | 2015-04-28 | Exelis Inc. | System and method for combining color information with spatial information in multispectral images |
| US10147170B2 (en) * | 2015-09-17 | 2018-12-04 | Raytheon Company | Systems and methods for sharpening multi-spectral imagery |
-
2019
- 2019-05-20 CN CN201910417005.7A patent/CN110310249B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102521815A (en) * | 2011-11-02 | 2012-06-27 | 薛笑荣 | Fast fusion system and fast fusion method for images |
| CN103530860A (en) * | 2013-09-26 | 2014-01-22 | 天津大学 | Adaptive autoregressive model-based hyper-spectral imagery super-resolution method |
| CN107220957A (en) * | 2017-04-25 | 2017-09-29 | 西北工业大学 | It is a kind of to utilize the remote sensing image fusion method for rolling Steerable filter |
Non-Patent Citations (3)
| Title |
|---|
| Spatial resolution enhancement of ASTER thermal bands;Aiazzi, B.et al.;《Proceedings of the SPIE - The International Society for Optical Engineering》;20050920;第5982卷;1-10 * |
| 一种多光谱图像和全色图像融合算法;边境 等;《微电子学与计算机》;20080720;第22卷(第6期);5-11 * |
| 基于Curvelet变换的多光谱图像与全色波段图像融合;张强 等;《系统工程与电子技术》;20061231;第28卷(第12期);1786-1789 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110310249A (en) | 2019-10-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101681613B (en) | Image color balance adjustment for display panels with 2d subpixel layouts | |
| DE60106197T2 (en) | COLOR SIGNAL PROCESSING | |
| CN110310249B (en) | Visual enhancement method for remote sensing image | |
| CN102509262B (en) | Method for removing thin cloud of remote sensing image | |
| Crippen et al. | Unveiling the lithology of vegetated terrains in remotely sensed imagery | |
| Suganya et al. | Survey on image enhancement techniques | |
| US6912307B2 (en) | Method for automatic color and intensity contrast adjustment of still and video images | |
| CN102436640A (en) | Foggy day image sharpening method based on HIS space multi-scale Retinex model | |
| CN113222054A (en) | Remote sensing image fusion method, system, equipment and medium based on characteristic ratio index | |
| CN101533599A (en) | Method for increasing gamma accuracy in quantized display systems | |
| CN102063700A (en) | Satellite remote sensing image generating method and system | |
| CN106384332A (en) | Method for fusing unmanned aerial vehicle image and multispectral image based on Gram-Schmidt | |
| DE102009001122B4 (en) | 1 - 15Camera arrangement and method for determining image signals with color values | |
| CN103268596A (en) | Method for reducing image noise and enabling colors to be close to standard | |
| CN107967668A (en) | A kind of image processing method and device | |
| Murata et al. | True color imagery rendering for Himawari-8 with a color reproduction approach based on the CIE XYZ color system | |
| CN107169946B (en) | Image fusion method based on nonnegative sparse matrix and hypersphere color transformation | |
| Toet | Colorizing single band intensified nightvision images | |
| CN112992039B (en) | Real-time acquisition method for optical radiation visual health information of display screen | |
| CN106023130A (en) | Gradient filtering and PCA-based unmanned aerial vehicle (UAV) image and multispectral image fusion method | |
| CN106204508A (en) | WorldView 2 remote sensing PAN and multi-spectral image interfusion method based on non-negative sparse matrix | |
| CN111563866A (en) | Multi-source remote sensing image fusion method | |
| CN102183244B (en) | Proportional light homogenizing method for aerial remote sensing image | |
| CN103824250B (en) | image tone mapping method based on GPU | |
| CN117935079A (en) | Remote sensing image fusion method, system and readable storage medium |
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 |