CN108154478A - A kind of remote sensing image processing method - Google Patents
A kind of remote sensing image processing method Download PDFInfo
- Publication number
- CN108154478A CN108154478A CN201611101877.5A CN201611101877A CN108154478A CN 108154478 A CN108154478 A CN 108154478A CN 201611101877 A CN201611101877 A CN 201611101877A CN 108154478 A CN108154478 A CN 108154478A
- Authority
- CN
- China
- Prior art keywords
- remote sensing
- sensing images
- value
- image
- pixel points
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 10
- 230000002401 inhibitory effect Effects 0.000 claims abstract description 22
- 230000004044 response Effects 0.000 claims abstract description 17
- 239000000284 extract Substances 0.000 claims abstract description 4
- 238000012545 processing Methods 0.000 claims description 22
- 230000000694 effects Effects 0.000 claims description 12
- 238000000034 method Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 3
- 230000005764 inhibitory process Effects 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 description 6
- 230000007850 degeneration Effects 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000004043 responsiveness Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G06T5/73—
-
- G06T5/80—
Abstract
The present invention provides a kind of remote sensing image processing method, including obtaining original remote sensing images I0And its image parameter；Obtain the satellite and camera parameter for shooting the remote sensing images；To original remote sensing images I0It is pre-processed, extracts remote sensing images I0On pixel value catastrophe point, using catastrophe point as target pixel points；Original remote sensing images I0In, the original pixel value based on target pixel points calculates the center response of the target pixel points, and the center response is more than the original pixel value；First total inhibiting value of the target pixel points is calculated, the center response is subtracted described first total inhibition is worth to difference, and using the difference as the pixel value of target pixel points described after image enhancement；Using 0 value as the pixel value of target pixel points described after image enhancement；Obtain the remote sensing images I after image procossing1.The present invention provides it is a kind of efficiently, quickly and accurately carry out remote sensing image processing method, picture quality can be efficiently modified.
Description
Technical field
The present invention relates to field of remote sensing image processing, and in particular, to a kind of remote sensing image processing method.
Background technology
With the rapid development of space technology and sensor technology, remote sensing image data has become people and obtains information
Important means plays more and more important in fields such as military surveillance, environmental monitoring, resource investigation, land use and urban plannings
Effect.However, the influence of the factors such as condition, atmospheric condition, sensor is illuminated by the light, the remote sensing figure acquired in imaging device
Picture and the atural object that is taken often there are certain deviations.Image fault can inevitably influence the correct of subsequent analysis and interpretation result
Property and stability.In face of the magnanimity remote sensing images exponentially increased at present, how automatically, remote sensing figure is quickly and stably handled
Picture becomes the hot issue of people's care and research and the difficulties with height challenge.Therefore, research and development are accurate, real
Remote sensing image processing method becomes particularly urgent and necessary.
In terms of image rectification, domestic and international researcher and technical staff have been carried out certain explore and study, and obtain
Preliminary achievement.Wherein more representational remote sensing software includes：The ERDAS IMAGINE and moral of ERDAS companies of the U.S.
The OrthoVista of INPHO companies of state, these softwares can improve the situation of remote sensing image；General commercial image processing software,
Mainly there is the Photoshop of Adobe companies.However, these softwares need a large amount of man-machine interactively and confirmation operation, including parameter
Setting, method choice, selection for whether subsequently being enhanced etc., handling result depend on the Heuristics or needs of operating personnel
It is operated repeatedly by multiple adjusting parameter, processing procedure is cumbersome, inefficient.With multi-platform, more spatial resolutions, multidate
The appearance of remote sensing images, daily remote sensing image data amount to be treated sharply increase, there is an urgent need to it is automatic, efficiently handle school
Positive technology, and existing technology is difficult in adapt to the needs of social development in terms of the degree of automation, processing speed and stability.
Invention content
Technological deficiency based on this field, the present invention provides a kind of remote sensing image processing methods, it is characterised in that：
Step 1, original remote sensing images I is obtained0And its image parameter；Obtain the satellite and camera for shooting the remote sensing images
Parameter；
Step 2, to original remote sensing images I0It is pre-processed, extracts the pixel value catastrophe point on remote sensing images I0, it will
Catastrophe point is as target pixel points；
Step 3, original remote sensing images I0In, the original pixel value based on target pixel points is calculated in the target pixel points
Heart response, the center response are more than the original pixel value；
Step 4, first total inhibiting value of the target pixel points is calculated, first total inhibiting value is the object pixel
For all neighborhood territory pixel point of point to the summation of the inhibiting value of the target pixel points, the neighborhood territory pixel point is that preassign be institute
The neighborhood territory pixel point of target pixel points is stated, the neighborhood territory pixel point refers to the neighborhood territory pixel to the inhibiting value of the target pixel points
Point plays target pixel points the value represented by inhibiting effect；
Step 5, when the center response is more than first total inhibiting value, the center response is subtracted described
First total inhibition is worth to difference, and using the difference as the pixel value of target pixel points described after image enhancement；When described
When center response is less than first total inhibiting value, using 0 value as the pixel value of target pixel points described after image enhancement.
Step 6, the remote sensing images I after image procossing is obtained1, and obtain its image parameter；
Step 7, image processing effect testing model is established, is tested to remote sensing image processing effect.
Preferably, wherein, step 2, to original remote sensing images I0It is pre-processed, if coloured image, pretreatment is also
Including image gray processing.
Preferably, wherein, step 2, to original remote sensing images I0It is pre-processed, including：Defogging processing, geometric correction.
Preferably, wherein, step 1, original remote sensing images I is obtained0And its image parameter, image parameter include resolution ratio.
Preferably, wherein, step 6, the remote sensing images I after image procossing is obtained1, and its image parameter is obtained, image ginseng
Number includes resolution ratio.
The present invention provides it is a kind of efficiently, quickly and accurately carry out remote sensing image processing method, image can be efficiently modified
Quality, and complete, accurate, objective inspection can be carried out to image processing effect.
Description of the drawings
Method flow diagram proposed by the invention Fig. 1.
Specific embodiment
For a better understanding of the present invention, with reference to the description of the embodiment of the accompanying drawings, the method for the present invention is carried out
Further instruction.
In order to fully understand the present invention, numerous details are referred in the following detailed description.But art technology
Personnel are it should be understood that the present invention may not need these details and realize.In embodiment, it is not described in detail well known side
Method, process, component, circuit, in order to avoid unnecessarily make embodiment cumbersome.
A kind of remote sensing image processing method shown in Figure 1, proposed by the invention, it is characterised in that：
Step 1, original remote sensing images I is obtained0And its image parameter；Obtain the satellite and camera for shooting the remote sensing images
Parameter；
Step 2, to original remote sensing images I0It is pre-processed, extracts the pixel value catastrophe point on remote sensing images I0, it will
Catastrophe point is as target pixel points；
Step 3, original remote sensing images I0In, the original pixel value based on target pixel points is calculated in the target pixel points
Heart response, the center response are more than the original pixel value；
Step 4, first total inhibiting value of the target pixel points is calculated, first total inhibiting value is the object pixel
For all neighborhood territory pixel point of point to the summation of the inhibiting value of the target pixel points, the neighborhood territory pixel point is that preassign be institute
The neighborhood territory pixel point of target pixel points is stated, the neighborhood territory pixel point refers to the neighborhood territory pixel to the inhibiting value of the target pixel points
Point plays target pixel points the value represented by inhibiting effect；
Step 5, when the center response is more than first total inhibiting value, the center response is subtracted described
First total inhibition is worth to difference, and using the difference as the pixel value of target pixel points described after image enhancement；When described
When center response is less than first total inhibiting value, using 0 value as the pixel value of target pixel points described after image enhancement.
Step 6, the remote sensing images I after image procossing is obtained1, and obtain its image parameter；
Step 7, image processing effect testing model is established, is tested to remote sensing image processing effect.
Preferably, wherein, step 2, to original remote sensing images I0It is pre-processed, if coloured image, pretreatment is also
Including image gray processing.
Preferably, wherein, step 2, to original remote sensing images I0It is pre-processed, including：Defogging processing, geometric correction.
Preferably, wherein, step 1, original remote sensing images I is obtained0And its image parameter, image parameter include resolution ratio.
Preferably, wherein, step 6, the remote sensing images I after image procossing is obtained1, and its image parameter is obtained, image ginseng
Number includes resolution ratio.
Preferably, wherein, the step 1 obtains the satellite and camera parameter for shooting the remote sensing images,
The satellite parametric reduction includes：Satellite designation, pointing accuracy, orbit altitude, panchromatic resolution ratio, multispectral resolution rate,
It is imaged breadth；
The camera parameter includes：Camera title, focal length, f-number, pixel, items MTF parameters.
Preferably, wherein, original remote sensing images is carried out with processing and includes image gray processing.
Preferably, wherein, the step 7 establishes image processing effect testing model, further includes calculating image processing effect
Test rating, the Testing index includes：Signal-to-noise ratio index, MTF indexs, information content index, dependent degeneration degree index.
Preferably, wherein, the signal-to-noise ratio index S/N is specially：
Wherein, f is camera focus, and F is camera aperture value, and B is the radiance of camera inlet, τopticsFor optical system
Unite transmitance, η be optical system cover bar ratio, NtdiSeries, t are integrated for TDICCDintFor TDICCD row integration periods, RCCDFor CCD
Responsiveness, PefFor camera valid pixel, Nshot、Ndark、Nread、NFPNThere are shots in respectively CCD and signal processing circuit to make an uproar
Sound, reads noise, fixed pattern noise at dark current noise.
Preferably, wherein, MTF is system modulation transmission function, and the MTF indexs are specially：
MTF=MTFIt is static·MTFDynamically
MTFIt is static=MTFOptical design·MTFOptical manufacturing·MTFOptics adjustment·MTFCCD device·MTFCircuit
Wherein, N is TDICCD integration series, and Δ p is improper as shifting amount, fNormalizationIt is normalized spatial frequency.
Preferably, wherein, described information content's index C is specially：
Wherein, the gray scale value set of each pixel i is { S in imagei, i=1,2 ..., n }, correspond to the probability occurred
For { Ai, i=1,2 ..., n }, W is imaging breadth；
Preferably, wherein, after dependent degeneration degree refers to image procossing, compared with original remote sensing images, picture quality is moved back
The degree of change, the dependent degeneration degree index RL are specially：
Wherein, S0、S1The gray value of pixel in respectively original remote sensing images and treated remote sensing images,
The mean value of pixel gray value in respectively original remote sensing images and treated remote sensing images, the ash of each pixel i in image
Angle value collection is combined into { Si, i=1,2 ..., n }, R0And R1The resolution ratio of respectively original remote sensing images and treated remote sensing images.
Preferably, wherein, the step 7, establishing image processing effect testing model IQT is specially：
Wherein, wiFor validity check parameters weighting input by user, Σ wi=1,For validity check parameters weighting wi
Fitting function, as the coefficient of validity check model indices, IQT scoring output areas are 0-100 points, wherein 100 points
Most satisfied for picture quality, 0 point least satisfied for picture quality, preferably using substar high quality graphic as 100 points with reference to mark
It is accurate.
As it can be seen that by the way that picture quality can be efficiently modified, and image processing effect can be carried out complete, accurate, objective
The inspection of sight.
Here the preferred embodiment of the present invention is only illustrated, but its meaning is not intended to limit the scope of the invention, applicability and is matched
It puts.On the contrary, detailed explanation of the embodiments can be implemented by those skilled in the art.It will be understood that without departing from appended power
In the case of the spirit and scope of the invention that sharp claim determines, changes and modifications may be made to details.
Claims (5)
1. a kind of remote sensing image processing method, it is characterised in that：
Step 1, original remote sensing images I is obtained0And its image parameter；Obtain the satellite and camera parameter for shooting the remote sensing images；
Step 2, to original remote sensing images I0It is pre-processed, extracts remote sensing images I0On pixel value catastrophe point, by catastrophe point
As target pixel points；
Step 3, original remote sensing images I0In, the center that the original pixel value based on target pixel points calculates the target pixel points is rung
It should be worth, the center response is more than the original pixel value；
Step 4, first total inhibiting value of the target pixel points is calculated, first total inhibiting value is the target pixel points institute
For some neighborhood territory pixel points to the summation of the inhibiting value of the target pixel points, the neighborhood territory pixel point is that preassign be the mesh
The neighborhood territory pixel point of pixel is marked, the neighborhood territory pixel point refers to the inhibiting value of the target pixel points neighborhood territory pixel point pair
Target pixel points play the value represented by inhibiting effect；
Step 5, when the center response is more than first total inhibiting value, the center response is subtracted described first
It is total to inhibit to be worth to difference, and using the difference as the pixel value of target pixel points described after image enhancement；When the center
When response is less than first total inhibiting value, using 0 value as the pixel value of target pixel points described after image enhancement.
Step 6, the remote sensing images I after image procossing is obtained1, and obtain its image parameter；
Step 7, image processing effect testing model is established, is tested to remote sensing image processing effect.
2. the method for claim 1, wherein step 2, to original remote sensing images I0It is pre-processed, if cromogram
Picture, pretreatment further include image gray processing.
3. the method for claim 1, wherein step 2, to original remote sensing images I0It is pre-processed, including：At defogging
Reason, geometric correction.
4. the method for claim 1, wherein step 1 obtains original remote sensing images I0And its image parameter, image parameter
Including resolution ratio.
5. the method for claim 1, wherein step 6 obtains the remote sensing images I after image procossing1, and obtain its image
Parameter, image parameter include resolution ratio.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611101877.5A CN108154478B (en) | 2016-12-02 | 2016-12-02 | A kind of remote sensing image processing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611101877.5A CN108154478B (en) | 2016-12-02 | 2016-12-02 | A kind of remote sensing image processing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108154478A true CN108154478A (en) | 2018-06-12 |
| CN108154478B CN108154478B (en) | 2019-10-08 |
Family
ID=62470835
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611101877.5A Active CN108154478B (en) | 2016-12-02 | 2016-12-02 | A kind of remote sensing image processing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108154478B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112950490A (en) * | 2021-01-25 | 2021-06-11 | 宁波市鄞州区测绘院 | Unmanned aerial vehicle remote sensing mapping image enhancement processing method |
| CN116664449A (en) * | 2023-07-26 | 2023-08-29 | 中色蓝图科技股份有限公司 | Satellite image processing method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104123697A (en) * | 2013-04-23 | 2014-10-29 | 华为技术有限公司 | Image enhancement method and device |
| CN104318526A (en) * | 2014-10-20 | 2015-01-28 | 南京理工大学 | MTF based remote sensor on-track automatic on-track parameter optimization method |
| CN105430378A (en) * | 2015-11-26 | 2016-03-23 | 航天恒星科技有限公司 | Image quality evaluation system and method |
-
2016
- 2016-12-02 CN CN201611101877.5A patent/CN108154478B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104123697A (en) * | 2013-04-23 | 2014-10-29 | 华为技术有限公司 | Image enhancement method and device |
| CN104318526A (en) * | 2014-10-20 | 2015-01-28 | 南京理工大学 | MTF based remote sensor on-track automatic on-track parameter optimization method |
| CN105430378A (en) * | 2015-11-26 | 2016-03-23 | 航天恒星科技有限公司 | Image quality evaluation system and method |
Non-Patent Citations (3)
| Title |
|---|
| 莫毅君，王一，明朝辉: "基于天绘一号卫星三线阵影像在轨几何检校与精度分析", 《数字通信世界》 * |
| 薛旭成 等: "空间遥感相机TDI CCD积分级数和增益的优化设置", 《光学精密工程》 * |
| 黄哲: "光学遥感卫星在轨参数与图像质量参数", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112950490A (en) * | 2021-01-25 | 2021-06-11 | 宁波市鄞州区测绘院 | Unmanned aerial vehicle remote sensing mapping image enhancement processing method |
| CN112950490B (en) * | 2021-01-25 | 2022-07-19 | 宁波市鄞州区测绘院 | Unmanned aerial vehicle remote sensing mapping image enhancement processing method |
| CN116664449A (en) * | 2023-07-26 | 2023-08-29 | 中色蓝图科技股份有限公司 | Satellite image processing method |
| CN116664449B (en) * | 2023-07-26 | 2023-10-13 | 中色蓝图科技股份有限公司 | Satellite image processing method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108154478B (en) | 2019-10-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gao et al. | Gradient-based interpolation method for division-of-focal-plane polarimeters | |
| CN108154479A (en) | A kind of method that remote sensing images are carried out with image rectification | |
| Chianucci | An overview of in situ digital canopy photography in forestry | |
| Pekin et al. | Measurement of crown cover and leaf area index using digital cover photography and its application to remote sensing | |
| Stritzinger et al. | The Carnegie Supernova Project I-Photometry data release of low-redshift stripped-envelope supernovae | |
| CN105389591B (en) | A kind of verification method of infrared imaging sensor typical effect emulation | |
| Li et al. | Using sub-resolution features for self-compensation of the modulation transfer function in remote sensing | |
| Minařík et al. | Radiometric and atmospheric corrections of multispectral μMCA camera for UAV spectroscopy | |
| Young et al. | Signal processing and performance analysis for imaging systems | |
| CN104200457A (en) | Wide-angle camera shooting based discrete type canopy leaf area index detection system and method | |
| CN104008543A (en) | Image fusion quality evaluation method | |
| CN104412298B (en) | Method and apparatus for changing image | |
| Song et al. | Improving accuracy of canopy hemispherical photography by a constant threshold value derived from an unobscured overcast sky | |
| CN108154478B (en) | A kind of remote sensing image processing method | |
| CN104813217B (en) | For the method for the passive single channel imager for being designed to estimate the depth of field | |
| CN110827375B (en) | Infrared image true color coloring method and system based on low-light-level image | |
| Jabari et al. | Improving UAV imaging quality by optical sensor fusion: an initial study | |
| CN108154470A (en) | A kind of remote sensing image processing method | |
| CN106705942B (en) | A method of examine remote sensing image to handle quality | |
| Carrasco et al. | The dwarf galaxy population in nearby groups: The data | |
| Popowicz | PSF photometry for BRITE nano-satellite mission | |
| CN106249220B (en) | Hot spot enhances processing method, device and laser range finder | |
| Wu et al. | Recovering sensor spectral sensitivity from raw data | |
| Manakov | Evaluation of computational radiometric and spectral sensor calibration techniques | |
| CN111504912A (en) | Air pollution observation system and method based on image recognition |
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 | ||
| CB02 | Change of applicant information |
Address after: 101399 No. 2 East Airport Road, Shunyi Airport Economic Core Area, Beijing (1st, 5th and 7th floors of Industrial Park 1A-4) Applicant after: Zhongke Star Map Co., Ltd. Address before: 101399 Building 1A-4, National Geographic Information Technology Industrial Park, Guomen Business District, Shunyi District, Beijing Applicant before: Space Star Technology (Beijing) Co., Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |