CN101493525A - Radiation correcting method for satellite remote sensing data - Google Patents
Radiation correcting method for satellite remote sensing data Download PDFInfo
- Publication number
- CN101493525A CN101493525A CNA2008102198420A CN200810219842A CN101493525A CN 101493525 A CN101493525 A CN 101493525A CN A2008102198420 A CNA2008102198420 A CN A2008102198420A CN 200810219842 A CN200810219842 A CN 200810219842A CN 101493525 A CN101493525 A CN 101493525A
- Authority
- CN
- China
- Prior art keywords
- water body
- remote sensing
- mountain area
- correcting method
- cleaning water
- 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.)
- Pending
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 37
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 42
- 238000001228 spectrum Methods 0.000 claims abstract description 10
- 238000004140 cleaning Methods 0.000 claims description 24
- 230000003287 optical effect Effects 0.000 abstract description 3
- 238000002310 reflectometry Methods 0.000 description 10
- 238000005259 measurement Methods 0.000 description 3
- 239000000443 aerosol Substances 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000000149 penetrating effect Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000005427 atmospheric aerosol Substances 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 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
- 238000005316 response function Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a radiation correction method used for satellite remote sensing data. The radiation correction method comprises the step of carrying out radiation correction on the original image by mountainous shade and/or clean water body spectrum. The radiation correction method of the satellite data can improve the radiation precision of the satellite remote sensing data by more than one magnitude; furthermore, the radiation correction method can be applicable to any optical sensors; and according to the test results, the method has high precision and achieves the requirement of practicability.
Description
Technical field
The invention belongs to the satellite remote sensing technology field, particularly a kind of radiation correcting method of satellite remote sensing date.
Background technology
After remote sensing satellite is launched successfully, reasons such as change have taken place owing to environmental factors such as temperature, pressure, the distribution of the response function of sensor changes, and has caused radiation value that satellite remote sensing obtains and the error between the actual emanations value.The method that addresses this problem at present has two, and the one, on satellite, measure direct sunlight and standard sources and come remotely-sensed data is carried out radiation correcting as the calibration parameter.The 2nd, by at ground and satellite synchro measure object spectrum, according to the object spectrum that records remotely-sensed data is carried out radiation correcting then.But owing to multiple reason, there is certain fluctuation inevitably in the emittance of standard lamp on the star, still has certain error after causing radiation correcting.Above-mentioned ground calibration method is because also there is tangible error in the variation of atmospheric aerosol.At present, after the calibration parameter that external satellite remote sensing date raw data provides was carried out radiation correcting, the reflectivity error still reached 1%~2%; The radiation error of homemade satellite more reaches 2%~4%.This error has produced CR Critical influence to the quantitative Analysis of aspects such as water quality remote sensing, becomes a big obstacle of water quality quantitative remote sensing.
Summary of the invention
The objective of the invention is to by utilizing atmospheric effect less, and stable mountain area shade and/or the cleaning water body spectrum of influence carries out radiation correcting, the radiation precision of satellite remote sensing date can be improved one more than the order of magnitude, and this kind radiation correcting method is applicable to any optical sensor.
The technical solution adopted in the present invention: a kind of radiation correcting method of satellite remote sensing date, described radiation correcting method are with mountain area shade and/or cleaning water body spectrum original image to be carried out radiation correcting.
Above-mentioned radiation correcting method may further comprise the steps:
(1) measures and calculates the pixel spoke brightness that obtains mountain area shade and/or cleaning water body via satellite;
(2) calculate the calibration parameter that obtains mountain area shade and/or cleaning water body;
(3) with the mountain area shade and/or the cleaning water body the calibration parameter to the original image spoke brightness carry out radiation correcting.
The pixel spoke brightness L of the mountain area shade in the above-mentioned steps (1)
sComputing formula be: L
s=L
g+ L
p+ L
e, in the formula: L
g, L
pAnd L
eBe respectively path scattering spoke brightness and the surrounding environment reflected light of ground return to sky reflection of light, atmosphere.
The pixel spoke brightness L of the cleaning water body in the above-mentioned steps (1)
wComputing formula be: L
w=L '
w+ L
p+ L
g+ L
Sky, in the formula: L '
w, L
pBe respectively the path scattering spoke brightness of the reflected light and the atmosphere of water body, L
gBe water surface sun credit light, L
SkyAs be the direct reflection spoke brightness of the water surface to skylight.
The calibration parameter b of the mountain area shade in the above-mentioned steps (2)
λComputing formula be: b
λ=L
S λ-L
λ, in the formula: L
S λBe the pixel spoke brightness of the mountain area shade of corresponding wave band, L
λBe the brightness of original image spoke.
The calibration parameter b of the cleaning water body in the above-mentioned steps (2)
λComputing formula be: b
λ=L
W λ-L
λ, in the formula: L
W λBe the pixel spoke brightness of the cleaning water body of corresponding wave band, L
λBe the brightness of original image spoke.
Error≤0.2% of the calibration parameter of above-mentioned mountain area shade and/or cleaning water body.
The radiation correcting method of satellite data of the present invention can improve the radiation precision of satellite remote sensing date one more than the order of magnitude, and this kind radiation correcting method is applicable to any optical sensor.Test findings shows that this method has higher precision, reaches the requirement of practicability.
Description of drawings
Fig. 1 is the image of the original mini-bus satellite in Delta of the Pearl River area;
Fig. 2 is the image through the mini-bus satellite in the Delta of the Pearl River area behind the radiation correcting of the present invention.
Embodiment
Below in conjunction with accompanying drawing concrete structure of the present invention is further described.
1. the principle of method
1.1 mountain area shade and cleaning water body spectrum characteristic
1.1.1 the spectrum characteristic of mountain area shade
In the atmospheric envelope, the composition of normal atmosphere molecule is stable, reduces but its concentration is negative exponent with the elevation increase, and its scattering nature is a Rayleigh scattering.With the normal atmosphere molecular mimicry, gasoloid is counted concentration also to be increased with elevation and is negative exponent and reduces, and the speed that reduces with height is faster.According to our measurement result, general area is more than 600 meters, and aerocolloidal concentration is just extremely low, and is mainly the stable background gas colloidal sol of horizontal distribution, and its scattering nature also is a Rayleigh scattering.In the shade of mountain area, the forest shade is modal shadow region on the remote sensing images.The shadow region direct sunlight is 0, and the spoke brightness of pixel is:
L
s=L
g+L
p+L
e
L
g, L
pAnd L
eBe respectively path scattering spoke brightness and the surrounding environment reflected light of ground return to sky reflection of light, atmosphere.
E
0Be solar irradiance; ω=ω
m+ ω
aThe whole atmosphere scattering coefficient equals gas and has divided scattering coefficient and aerosol dispersion to penetrate the coefficient sum; P (θ) for Rayleigh with penetrating phase function, the mountain area aerosol particle that is higher than 600 meters for height above sea level is superfine, is Rayleigh scattering also, has
E
SkyBe the skylight irradiance, and
R
vBe mountain area forest cover reflectivity, can be by measuring.E
eFor surrounding environment reflection irradiance, depend on exclosure type and clutter reflections rate, T=e
-ωBe atmospheric transmittance, T was less when height above sea level was higher, and T=1-ω is arranged.
For visible light wave range, the E of forest shade
eMinimum, can ignore.
λ is a wavelength; H
m=8.44km is atmosphere folding letter height; Z is the sea level elevation of pixel.
β, H
aRespectively background gas colloidal sol turbidity coefficient and folding letter height, the two is to certain geographic position and be season than steadiness parameter, can be by measuring.
1.1.2 mountain area cleaning water body spectrum characteristic
The spoke brightness of mountain area clean water body image unit is:
L
w=L′
w+L
p+L
g+L
sky
L '
w, L
pBe respectively the path scattering spoke brightness of the reflected light and the atmosphere of water body, L
gBe water surface sun credit light.
Be the direct reflection spoke brightness of the water surface to skylight, R
Mw=0.025 is the specular reflectance of the water surface to skylight.For calm water surface, sun credit light can be ignored:
R
wBe the clean water volume reflectivity, can be by measuring.Other CALCULATION OF PARAMETERS is the same.
1.2 based on the radiation correcting of mountain area shade with cleaning water body spectrum
1.2.1 visible light wave range radiation correcting method based on the mountain area shade
Use above formula, can more accurately calculate each wave band spoke brightness L of visible light wave range forest shade pixel
S λCan obtain each wave band calibration parameter b then
λ
b
λ=L
sλ-L
λ
L
λBe the brightness of original image spoke.
1.2.1 radiation correcting method based on mountain area cleaning water body
Visible light wave range radiation correcting method based on the mountain area shade is high at the reflectivity of near-infrared band owing to vegetation, and the surrounding environment reflected light can not be applicable to near-infrared band by force.The available following radiation correcting method of the correction of near-infrared band based on mountain area cleaning water body.For the cleaning water body, near-infrared band is strong the absorption, and its reflectivity is 0, so the pixel spoke brightness of water body only is the air path scattering:
The calibration parameter is: b
λ=L
W λ-L
λ
Regional widely in the waters, can b
λAll sweep traces are corrected.
2. the estimation of error of method
2.1 visible light wave range radiation correcting estimation of error based on the mountain area shade
This method is mainly from the error of shadow region vegetation albedo measurement value and actual reflectance.The Different Forest shadow region is because vegetation type or coverage rate are variant, its reflectivity has certain error, but at visible light wave range since its to penetrate reflectivity very low, add that measurement of reflectivity measures in actual shadow region, the reflected light of atural object also obtains considering to a certain extent on every side, all obtains under bright day gas condition so the error of visible light wave range reflectivity in 2%, adds qualified image, skylight only is 1/10 of a direct sunlight, so error≤0.2% of calibration parameter.
2.2 near-infrared band radiation correcting estimation of error based on mountain area cleaning water body
This method is mainly from the influence of steam on the water body.The steam scattered light can be calibrated CALCULATION OF PARAMETERS to this method and cause certain error.But under the sunny weather, steam mainly exists with the gas pattern, and precision is thin, and with penetrating, the scattering strength of infrared band is minimum based on Rayleigh, and its amplitude of variation is littler, so the error of calibration parameter also≤0.2%.
3. application example
We utilize the radiation calibration to TM of Guangdong Province and mini-bus satellite image data, and the result shows that this method has higher precision, reaches the requirement of practicability.Fig. 1 is the image of the original mini-bus satellite in Delta of the Pearl River area, and Fig. 2 utilizes correcting method that above-mentioned Delta of the Pearl River mini-bus satellite image is carried out radiation correcting example afterwards.
Can find out obviously that from figure original image is deposited significantly band of ground, on the same image, wave band 1 and wave band 4 errors are all above 3%, and behind the radiation correcting with the carrying out of this method, the radiation error of reflectivity reduces in 0.2%.
Claims (7)
1. the radiation correcting method of a satellite remote sensing date is characterized in that, described radiation correcting method is with mountain area shade and/or cleaning water body spectrum original image to be carried out radiation correcting.
2. according to the radiation correcting method of the described satellite remote sensing date of claim 1, it is characterized in that, may further comprise the steps:
(1) measures and calculates the pixel spoke brightness that obtains mountain area shade and/or cleaning water body via satellite;
(2) calculate the calibration parameter that obtains mountain area shade and/or cleaning water body;
(3) with the mountain area shade and/or the cleaning water body the calibration parameter to the original image spoke brightness carry out radiation correcting.
3. according to the radiation correcting method of the described satellite remote sensing date of claim 2, it is characterized in that the pixel spoke brightness I of the mountain area shade in the described step (1)
sComputing formula be: L
s=L
g+ L
p+ L
e
In the formula: L
g, L
pAnd L
eBe respectively path scattering spoke brightness and the surrounding environment reflected light of ground return to sky reflection of light, atmosphere.
4. according to the radiation correcting method of the described satellite remote sensing date of claim 2, it is characterized in that the pixel spoke brightness L of the cleaning water body in the described step (1)
wComputing formula be: L
w=L '
w+ I
p+ L
g+ L
Sky
In the formula: L '
w, L
pBe respectively the path scattering spoke brightness of the reflected light and the atmosphere of water body, L
gBe water surface sun credit light, L
SkyBe the direct reflection spoke brightness of the water surface to skylight.
5. according to the radiation correcting method of the described satellite remote sensing date of claim 2, it is characterized in that the calibration parameter b of the mountain area shade in the described step (2)
λComputing formula be: b
λ=L
S λ-L
λ
In the formula: L
S λBe the pixel spoke brightness of the mountain area shade of corresponding wave band, L
λBe the brightness of original image spoke.
6. according to the radiation correcting method of the described satellite remote sensing date of claim 2, it is characterized in that the calibration parameter b of the cleaning water body in the described step (2)
λComputing formula be: b
λ=L
W λ-L
λ
In the formula: L
W λBe the pixel spoke brightness of the cleaning water body of corresponding wave band, L
λBe the brightness of original image spoke.
7. according to the radiation correcting method of claim 5 or 6 described satellite remote sensing dates, it is characterized in that error≤0.2% of the calibration parameter of described mountain area shade and/or cleaning water body.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2008102198420A CN101493525A (en) | 2008-12-11 | 2008-12-11 | Radiation correcting method for satellite remote sensing data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2008102198420A CN101493525A (en) | 2008-12-11 | 2008-12-11 | Radiation correcting method for satellite remote sensing data |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101493525A true CN101493525A (en) | 2009-07-29 |
Family
ID=40924215
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2008102198420A Pending CN101493525A (en) | 2008-12-11 | 2008-12-11 | Radiation correcting method for satellite remote sensing data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101493525A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103363962A (en) * | 2013-07-23 | 2013-10-23 | 中国科学院遥感与数字地球研究所 | Remote sensing evaluation method of lake water reserves based on multispectral images |
CN105092476A (en) * | 2015-08-20 | 2015-11-25 | 中山大学 | Method for simultaneously inverting turbidity, COD and chlorophyll concentration of inland water |
CN105205789A (en) * | 2015-08-20 | 2015-12-30 | 中山大学 | Method for eliminating specular reflection influence of water area remote sensing data |
CN105259145A (en) * | 2015-08-20 | 2016-01-20 | 中山大学 | Method for simultaneous remote sensing of underwater terrain and features of island |
CN110006463A (en) * | 2019-05-23 | 2019-07-12 | 中国科学院合肥物质科学研究院 | A kind of in-orbit absolute radiation calibration method and system of Optical remote satellite |
-
2008
- 2008-12-11 CN CNA2008102198420A patent/CN101493525A/en active Pending
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103363962A (en) * | 2013-07-23 | 2013-10-23 | 中国科学院遥感与数字地球研究所 | Remote sensing evaluation method of lake water reserves based on multispectral images |
CN103363962B (en) * | 2013-07-23 | 2015-04-08 | 中国科学院遥感与数字地球研究所 | Remote sensing evaluation method of lake water reserves based on multispectral images |
CN105092476A (en) * | 2015-08-20 | 2015-11-25 | 中山大学 | Method for simultaneously inverting turbidity, COD and chlorophyll concentration of inland water |
CN105205789A (en) * | 2015-08-20 | 2015-12-30 | 中山大学 | Method for eliminating specular reflection influence of water area remote sensing data |
CN105259145A (en) * | 2015-08-20 | 2016-01-20 | 中山大学 | Method for simultaneous remote sensing of underwater terrain and features of island |
CN105205789B (en) * | 2015-08-20 | 2017-12-26 | 中山大学 | A kind of method for eliminating water area remote sensing data specular reflections effect |
CN105259145B (en) * | 2015-08-20 | 2018-06-12 | 中山大学 | A kind of while remote sensing islands and reefs underwater topography and the method for atural object |
CN110006463A (en) * | 2019-05-23 | 2019-07-12 | 中国科学院合肥物质科学研究院 | A kind of in-orbit absolute radiation calibration method and system of Optical remote satellite |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101598543B (en) | Practical atmospheric correction method for remote sensing images | |
Gui et al. | Evaluation of radiosonde, MODIS-NIR-Clear, and AERONET precipitable water vapor using IGS ground-based GPS measurements over China | |
He et al. | Using geostationary satellite ocean color data to map the diurnal dynamics of suspended particulate matter in coastal waters | |
CN110006463B (en) | On-orbit absolute radiation calibration method and system for optical remote sensing satellite | |
Mueller et al. | The CM-SAF operational scheme for the satellite based retrieval of solar surface irradiance—A LUT based eigenvector hybrid approach | |
Higurashi et al. | A study of global aerosol optical climatology with two-channel AVHRR remote sensing | |
CN103198314B (en) | Remote sensing images radiation correction method | |
CN104406686B (en) | Sun shortwave incident radiation evaluation method under MODEL OVER COMPLEX TOPOGRAPHY | |
CN102338869B (en) | Inversion method and system of downlink shortwave radiation and photosynthetically active radiation data | |
Xie et al. | Calculating NDVI for Landsat7-ETM data after atmospheric correction using 6S model: A case study in Zhangye city, China | |
CN101915914B (en) | Lookup table based pixel-by-pixel atmospheric correction method of remote sensing images | |
CN103018736A (en) | Satellite-borne remote sensor radiation calibration method based on atmospheric parameter remote sensing retrieval | |
CN104156567B (en) | Technique for acquiring surface reflectance by coupling satellite remote-sensing image atmospheric correction and topographical correction processes | |
CN1948953B (en) | Based on laser rediffusion type atmospheric visibility testing method and testing instrument thereof | |
CN102288956A (en) | Atmospheric correction method for multispectral data of remote sensing satellite | |
CN103499815A (en) | Method for conducting inland water body atmospheric correction based on oxygen and water-vapor absorption wavebands | |
CN101493525A (en) | Radiation correcting method for satellite remote sensing data | |
Soufflet et al. | Atmospheric effects on NOAA AVHRR data over Sahelian regions | |
CN114218786B (en) | On-orbit polarized radiation characteristic inversion method for unpolarized satellite sensor | |
CN116664947A (en) | Blue algae bloom monitoring method and system based on satellite observation data | |
CN102073792A (en) | Method for performing inversion on optical property of aerosol of coastal zone by using MODIS image | |
CN106198460B (en) | Random fluctuation sea surface albedo full-spectrum segment calculation method applied to marine environment | |
Jakobson et al. | Diurnal variability of precipitable water in the Baltic region, impact on transmittance of the direct solar radiation | |
CN113740263A (en) | Aerosol optical thickness inversion method and atmospheric particulate matter remote sensing inversion method | |
Nahmani et al. | Sensitivity of GPS tropospheric estimates to mesoscale convective systems in West Africa |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20090729 |