CN1710379A - Atmosphere correction method of airosol optical thickness of aeronautical high-spectrum remote-sensing inversion boundary layer - Google Patents
Atmosphere correction method of airosol optical thickness of aeronautical high-spectrum remote-sensing inversion boundary layer Download PDFInfo
- Publication number
- CN1710379A CN1710379A CN 200510027524 CN200510027524A CN1710379A CN 1710379 A CN1710379 A CN 1710379A CN 200510027524 CN200510027524 CN 200510027524 CN 200510027524 A CN200510027524 A CN 200510027524A CN 1710379 A CN1710379 A CN 1710379A
- Authority
- CN
- China
- Prior art keywords
- module
- aerosol
- optical depth
- atmospheric correction
- boundary layer
- 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
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
A method for inverting optical thickness of boundary layer aerosol includes picking up boundary layer aerosol optical thickness data from atmosphere transmission spectrum by utilizing transmission and radiation principle of atmosphere, obtaining expression reflectivity by carrying out reinterpolation for atmosphere and aerosol modes in boundary layer, thus realizing inversion of calculation for boundary layer aerosol optical thickness . The atmospheric calibration can be carried out based on result calculated out by above method.
Description
Technical field
The present invention relates to a kind of atmospheric correction method of aviation high-spectrum remote-sensing inversion boundary layer aerosol optical depth, the high spectrum image that the aviation high-spectrum remote-sensing device that this method is applied to develop Shanghai Institute of Technical Physics of Chinese Academy of Sciences airborne remote sensing research department obtains, by atmospheric correction method based on radiation delivery mechanism, realize the inverting of urban atmosphere boundary-layer aerosol optical depth, belong to atmospheric environment remote sensing application field.
Background technology
Atmospheric aerosol is meant the radius that suspends in the atmosphere less than tens microns solid-state or liquid particulate, and gasoloid is being played the part of important role in earth atmosphere radiation budget balance and global climate, be the important research object in the atmospheric physics.On the one hand, gasoloid by scattering and absorb solar radiation and terrestrial radiation directly influence land-the radiation budget balance of gas system; On the other hand, gasoloid also participates in a plurality of physical processes of atmosphere, manages mechanism and ozone equilibrium etc. as the speck that cloud and mist forms; With the signal that absorbs and scattering method is disturbing remote sensor to receive.Therefore, accurately the Measurement and analysis gasoloid for understanding climate change, is removed the atmospheric effect in the remotely-sensed data, and it is all significant to improve remote sensing quantitative Application level.
Aerosol optical depth is one of most important parameter of gasoloid, is the important physical amount that characterizes atmospheric turbidity, also is to determine a key factor of gasoloid climatic effect and an important parameter of Atmospheric models.Survey aerosol optical depth and can adopt the ground based detection method, as actinometer, corpuscular counter, radiation summary table etc.Though the ground based detection method can accurately provide local gasoloid information, can not obtain interior gasoloid spatial and temporal distributions on a large scale.The remote-sensing inversion aerosol optical depth can overcome the inadequate natural endowment of ground based detection method, provides possibility for the gasoloid in round-the-clock, the real-time understanding of people on a large scale changes.
In recent years, the remote-sensing inversion aerosol optical depth has become the means that obtain atmospheric aerosol information quickly and efficiently, has especially obtained the excellent research achievement aspect satellite remote sensing, and the inversion algorithm of comparative maturity has been arranged
[5-15], but mainly be to adopt dark pixel method (or dark goal method) to set up look-up table by 6S to realize aerocolloidal inverting, and these algorithms are mostly at satellite data.Dark pixel method utilizes most of lands face in red (0.6 ~ 0.68 μ m) and blue (0.40 ~ 0.48 μ m) characteristic that the wave band reflectivity is low, with vegetation index (NDVI) or near infrared passage (2.1 μ m) reflectivity forest being declared knowledge is dark pixel, is used for the inverting aerosol optical depth
[6]Chaetomers Thailand
[13]When test utilizes dark pixel method inverting Beijing, Hong Kong municipal pollution aerosol optical depth, think, determine that with fixed proportion coefficient relational expression the dark pixel of vegetation is red, the method for blue channel earth surface reflection rate exists than mistake by near infrared passage apparent reflectance in the Beijing area.This shows and utilizes dark pixel method inverting Beijing area aerosol optical depth to have certain difficulty.Control methods is to study the satellite remote sensing method that the land pollution gasoloid adopts in early days
[16]
Say that in principle U.S. NASA can provide the aerosol optical depth in most of area, the whole world with the MODIS image, but its spatial resolution is 10km only, and satellite altitude is at more than 700 kilometers, what obtain is aerosol optical depth in the whole troposphere, and gasoloid mainly concentrates on from ground to the urban boundary layer in such vertical range, therefore becomes the problem that we pay close attention to from airborne high spectrum image inverting urban boundary layer aerosol optical depth.
Practical modularization imaging spectrometer (OMIS-I) is developed by Shanghai Institute of Technical Physics of Chinese Academy of Sciences airborne remote sensing research department, have 128 continuous spectrum passages from visible, near infrared to thermal infrared (wavelength coverage is at 0.46-12.5 μ m), spectral resolution reaches the 10nm rank, and the detailed technology index sees Table-1.The OMIS average flight altitude is 2km, and ground instantaneous field of view is 3mrad, and total visual field is spent greater than 70, the about 6m*6m of the ground resolution of center pixel, and the spatial resolution of image border point is low.
Summary of the invention
In order to overcome above-mentioned weak point, fundamental purpose of the present invention aims to provide a kind of practical modularization imaging spectrometer (OMIS-I) at the development of Shanghai Institute of Technical Physics of Chinese Academy of Sciences airborne remote sensing research department, in conjunction with the atmosphere radiation transmission principle, from atmospheric transmission spectrum, extract the atmospheric correction method of the aviation high-spectrum remote-sensing inversion boundary layer aerosol optical depth of boundary layer aerosol optical depth information.
The technical problem to be solved in the present invention is: the radiation calibration problem that solve high-spectrum remote sensing; Solve and how on the practical modularization imaging spectrometer of aviation high-spectrum remote-sensing device 0MIS image, to choose dark target, how to obtain showing reflectivity, problems such as inversion boundary layer aerosol optical depth calculating; Solve and how from atmospheric transmission spectrum, to extract boundary layer aerosol optical depth information, carry out relevant technologies problems such as atmospheric correction.
The technical solution adopted for the present invention to solve the technical problems is: the atmospheric correction method of aviation high-spectrum remote-sensing inversion boundary layer aerosol optical depth is by the atmosphere radiation transmission principle, from atmospheric transmission spectrum, extract boundary layer aerosol optical depth information, by being carried out once more interpolation in the boundary layer, atmospherical model, gasoloid pattern obtain showing reflectivity, realize the calculating of inversion boundary layer aerosol optical depth, its concrete calculation procedure is:
Step 1: the radiation calibration of high-spectrum remote sensing
A), read high-spectrum remote sensing
Read airborne aviation high-spectrum remote sensing by standard format;
B), convert radiation value to
The output signal of high-spectrum remote sensing module is sent to the input end of radiation value module, slope and two coefficients of intercept according to each wave band correspondence in the radiation calibration file, according to formula radiation value=DN value * slope+intercept, convert the digital number DN value of high-spectrum remote sensing to radiation value;
Step 2, calculating apparent reflectance are chosen dark target
A), calculate
Radiation value with step 1 output is calculated as the apparent reflectance that obtains on the sensor by following formula:
R=π*L/(μ*f)
In the formula: R is an apparent reflectance;
L is a radiation value;
μ is the cosine of solar zenith angle;
F is an aeropause solar radiation flux density;
B), choose dark target
The output signal of radiation value module is sent to dark target apparent reflectance module, chooses the apparent reflectance on the near infrared 2.1um passage, and the pixel in 0.036~0.044 scope is dark target;
Step 3, atmospheric correction
The output signal of dark target apparent reflectance module is sent to the atmospheric correction module, and the dark target apparent reflectance that obtains according to step 2 carries out atmospheric correction; Its concrete job step is as follows:
A), definition geometric parameter
The input aircraft is with respect to elevation angle, the zenith angle of pixel in the geometric parameter module, and sun altitude, zenith angle and observation date, the output signal of its geometric parameter module is sent to the atmospheric correction module;
B), definition atmospherical model
The atmospherical model that defining mode adopted in the atmospherical model module, the output signal of its atmospherical model module is sent to the atmospheric correction module;
C), definition aerosol type pattern
Delustring and meteorological range type in aerosol type in the used boundary layer gasoloid pattern of defining mode, if defined visibility simultaneously, then replace the meteorological range of default value definition in the aerosol type, the output signal of aerosol type module is sent to the atmospheric correction module;
D), definition aerosol concentration pattern
The used aerosol concentration pattern of defining mode in aerosol concentration, aerosol optical depth or the horizontal visibility of input 550nm, the output signal of aerosol concentration module is sent to the atmospheric correction module;
E), input atural object elevation
The output signal of input atural object elevation module is sent to the atmospheric correction module;
F), input aircraft flight height
The output signal of input aircraft flight altitude module is sent to the atmospheric correction module;
Step 4, atmospheric correction parameter extraction
A), signal transmits
The output signal of atmospheric correction module is sent to aerosol scattering phase function module and gasoloid single scattering reflectivity module respectively;
B), carry out atmospheric correction
According to each input parameter in the step 3, the dark target on the aviation high-spectrum remote sensing of selecting is carried out atmospheric correction;
C), correction parameter
From the output result, extract correction parameter aerosol scattering phase function and gasoloid single scattering reflectivity;
Step 5, computation bound layer aerosol optical depth
A), signal transmits
The output signal of aerosol scattering phase function module and gasoloid single scattering reflectivity module is sent to the aerosol optical depth module;
B), calculate
Apparent reflectance by following formula computation bound atmosphere:
C), computation bound layer aerosol optical depth
Aerosol scattering phase function and gasoloid single scattering reflectivity by step 4 draws calculate the boundary layer aerosol optical depth;
Step 6, judgement
A), signal transmits
The output signal of aerosol optical depth module is sent to comparison module;
B), calculate the ground sight distance
With synchronous ground sight distance data according to formula
V=3.91·H·1/τ
In the formula: V is the ground sight distance;
H is the gasoloid absolute altitude of Various Seasonal;
τ is an aerosol optical depth;
C), relatively
The aerosol optical depth that the aerosol optical depth value that calculates and step 5 are obtained compares;
When both difference<0.1, the output signal of comparison module is sent to output module, be output as aerosol optical depth;
Otherwise feed back to the atmospheric correction module, get back to step 3 and carry out atmospheric correction again, re-enter the needed parameters of atmospheric correction, circulation repeats.
The atmospherical model of the atmospheric correction method of described aviation high-spectrum remote-sensing inversion boundary layer aerosol optical depth comprises: tropical atmosphere, middle latitude summer, middle latitude winter, subpolar zone summer, subpolar zone winter and 1976 United States standard atmospheres etc.
The aerosol type of the atmospheric correction method of described aviation high-spectrum remote-sensing inversion boundary layer aerosol optical depth comprises: or be no gasoloid; Or be rural extinction coefficient, default meteorological range=23km; Or be the city extinction coefficient, default meteorological range=5km; Or be the troposphere extinction coefficient, default meteorological range=50km.
The invention has the beneficial effects as follows: do not have so far at airborne high light spectrum image-forming data, especially OMIS imaging spectrometer, carry out the method for the detailed complete of aerosol optical depth inverting, the method with the detailed complete of carrying out the aerosol optical depth inverting of the present invention; And,, from atmospheric transmission spectrum, extract aerosol optical depth information and carry out atmospheric correction in conjunction with boundary layer atmosphere radiation transmission principle at the OMIS imaging spectrometer.
Description of drawings
The present invention is further described below in conjunction with drawings and Examples.
Accompanying drawing 1 is an inversion boundary layer aerosol optical depth calculation process synoptic diagram of the present invention;
Accompanying drawing 2 is position views shown in the dark surface in the embodiment of the invention;
Accompanying drawing 3 be in the embodiment of the invention dark surface respectively at 13 and 14 s' apparent reflectance curve synoptic diagram;
Accompanying drawing 4 is that the atmospheric transmittance that calculates of the embodiment of the invention is with the wavelength change curve synoptic diagram;
Accompanying drawing 5 is aerosol optical depth synoptic diagram that embodiment of the invention atmospheric correction method according to the present invention is finally inversed by from airborne OMIS high-spectrum remote sensing;
Accompanying drawing 6 atmospheric transmittance that to be embodiment of the invention calculate based on LOWTRAN and comparison synoptic diagram based on the atmospheric transmittance of OMIS;
Label declaration in the accompanying drawing:
The 10-high-spectrum remote sensing;
The 20-radiation value;
The dark target apparent reflectance of 30-;
The 40-atmospheric correction;
The 41-geometric parameter;
The 42-atmospherical model;
The 43-aerosol type;
The 44-aerosol concentration;
45-atural object elevation;
46-aircraft flight height;
51-aerosol scattering phase function;
52-gasoloid single scattering reflectivity;
The 60-aerosol optical depth;
70-relatively;
80-output;
Embodiment
See also shown in the accompanying drawing 1, the present invention is a kind of atmospheric correction method of aviation high-spectrum remote-sensing inversion boundary layer aerosol optical depth, this method is by the atmosphere radiation transmission principle, from atmospheric transmission spectrum, extract boundary layer aerosol optical depth information, by being carried out once more interpolation in the boundary layer, atmospherical model, gasoloid pattern obtain showing reflectivity, realize the calculating of inversion boundary layer aerosol optical depth, its concrete calculation procedure is:
Step 1: the radiation calibration of high-spectrum remote sensing
A), read high-spectrum remote sensing (10)
Read airborne aviation high-spectrum remote sensing (10) by standard format;
B), convert radiation value (20) to
The output signal of high-spectrum remote sensing (10) module is sent to the input end of radiation value (20) module, slope and two coefficients of intercept according to each wave band correspondence in the radiation calibration file, according to formula radiation value=DN value * slope+intercept, convert the digital number DN value of high-spectrum remote sensing (10) to radiation value (20);
Step 2, calculating apparent reflectance are chosen dark target
A), calculate
Radiation value (20) with step 1 output is calculated as the apparent reflectance that obtains on the sensor by following formula:
R=π*L/(μ*f)
In the formula: R is an apparent reflectance;
L is a radiation value;
μ is the cosine of solar zenith angle;
F is an aeropause solar radiation flux density;
B), choose dark target
The output signal of radiation value (20) module is sent to dark target apparent reflectance (30) module, chooses the apparent reflectance on the near infrared 2.1um passage, and the pixel in 0.036~0.044 scope is dark target;
Step 3, atmospheric correction (40)
The output signal of dark target apparent reflectance (30) module is sent to atmospheric correction (40) module, and the dark target apparent reflectance (30) that obtains according to step 2 carries out atmospheric correction (40); Its concrete job step is as follows:
A), definition geometric parameter (41)
The input aircraft is with respect to elevation angle, the zenith angle of pixel in geometric parameter (41) module, and sun altitude, zenith angle and observation date, the output signal of its geometric parameter (41) module is sent to atmospheric correction (40) module;
B), definition atmospherical model (42)
The atmospherical model that defining mode adopted in atmospherical model (42) module, the output signal of its atmospherical model (42) module is sent to atmospheric correction (40) module;
C), definition aerosol type (43) pattern
Delustring and meteorological range type in aerosol type (43) in the used boundary layer gasoloid pattern of defining mode, if defined visibility simultaneously, then replace the meteorological range of default value definition in the aerosol type, the output signal of aerosol type (43) module is sent to atmospheric correction (40) module;
D), definition aerosol concentration (44) pattern
The used aerosol concentration pattern of defining mode in aerosol concentration (44), aerosol optical depth or the horizontal visibility of input 550nm, the output signal of aerosol concentration (44) module is sent to atmospheric correction (40) module;
E), input atural object elevation (45)
The output signal that to import atural object elevation (45) module is sent to atmospheric correction (40) module;
F), input aircraft flight height (46)
The output signal that to import aircraft flight height (46) module is sent to atmospheric correction (40) module;
Step 4, atmospheric correction parameter extraction
A), signal transmits
The output signal of atmospheric correction (40) module is sent to aerosol scattering phase function (51) module and gasoloid single scattering reflectivity (52) module respectively;
B), carry out atmospheric correction
According to each input parameter in the step 3, the dark target on the aviation high-spectrum remote sensing of selecting is carried out atmospheric correction;
C), correction parameter
From the output result, extract correction parameter aerosol scattering phase function (51) and gasoloid single scattering reflectivity (52);
Step 5, computation bound layer aerosol optical depth (60)
A), signal transmits
The output signal of aerosol scattering phase function (51) module and gasoloid single scattering reflectivity (52) module is sent to aerosol optical depth (60) module;
B), calculate
Apparent reflectance by following formula computation bound atmosphere:
C), computation bound layer aerosol optical depth (60)
Aerosol scattering phase function (51) and gasoloid single scattering reflectivity (52) by step 4 draws calculate boundary layer aerosol optical depth (60);
Step 6, judgement
A), signal transmits
The output signal of aerosol optical depth (60) module is sent to comparison (70) module;
B), calculate the ground sight distance
With synchronous ground sight distance data according to formula
V=3.91·H·1/τ
In the formula: V is the ground sight distance;
H is the gasoloid absolute altitude of Various Seasonal;
τ is an aerosol optical depth;
C), relatively
The aerosol optical depth (60) that the aerosol optical depth value that calculates and step 5 are obtained compares (70);
When both difference<0.1, the output signal that will compare (70) module is sent to output (80) module, is output as aerosol optical depth (60);
Otherwise feed back to atmospheric correction (40) module, get back to step 3 and carry out atmospheric correction again, re-enter the needed parameters of atmospheric correction, circulation repeats.
The atmospherical model (42) of the atmospheric correction method of described aviation high-spectrum remote-sensing inversion boundary layer aerosol optical depth comprises: tropical atmosphere, middle latitude summer, middle latitude winter, subpolar zone summer, subpolar zone winter and 1976 United States standard atmospheres etc.
The aerosol type (43) of the atmospheric correction method of described aviation high-spectrum remote-sensing inversion boundary layer aerosol optical depth comprises: or be no gasoloid; Or be rural extinction coefficient, default meteorological range=23km; Or be the city extinction coefficient, default meteorological range=5km; Or be the troposphere extinction coefficient, default meteorological range=50km.
Principle of work of the present invention is: the airborne high-spectrum remote-sensing platform of practical modularization imaging spectrometer OMIS is an aircraft, its flying height is roughly on the 2km height, basically belong to the boundary layer scope, therefore, can be according to airborne high-spectrum remote sensing data, by the atmosphere radiation transporting mechanism, the inversion boundary layer aerosol optical depth.
At the apparent reflectance of boundary layer atmosphere, can be written as
Because aircraft is in the boundary layer, therefore, shortened the path of clutter reflections radiation through atmosphere arrival sensor, so up radiation path need be made amendment to remove the atmospheric action factor more than the aircraft altitude in the up radiative process, this can obtain by atmospherical model, gasoloid pattern are carried out once more interpolation in the boundary layer.
Because the boundary layer intensive amount is very little, therefore can not consider; But the moisture content conversion is obvious in the boundary layer, and is responsive if the observation passage changes steam, then must consider the influence of steam real-time change to apparent reflectance, and need carry out the simultaneous observation of boundary layer steam content this moment aboard.
The concrete job step of the embodiment of the invention:
Use said method and choose dark target from the aviation high-spectrum remote-sensing device OMIS image of Shanghai Institute of Technical Physics of Chinese Academy of Sciences airborne remote sensing research department development, inversion boundary layer aerosol optical depth calculation process as shown in Figure 1:
Step 1: the radiation calibration of high-spectrum remote sensing
The airborne aviation high-spectrum remote sensing (10) of the standard format that provides by Shanghai Institute of Technical Physics of the Chinese Academy of Sciences is provided, slope and two coefficients of intercept according to each wave band correspondence in the radiation calibration file, according to formula radiation value=DN value * slope+intercept, convert digital number DN value to radiation value (20);
Step 2, calculating apparent reflectance are chosen dark target
With step 1 output radiation value (20), according to formula
R=π*L/(μ*f)
Wherein, R is an apparent reflectance; L is a radiation value; μ is the cosine of solar zenith angle; F is an aeropause solar radiation flux density; Be calculated to be the apparent reflectance that obtains on the sensor, the pixel of apparent reflectance in 0.036~0.044 scope of choosing then on the near infrared 2.1um passage is dark target.
Step 3, atmospheric correction
The apparent reflectance of the dark target that obtains according to step 2 carries out atmospheric correction.Step is as follows:
Step 3-1, definition geometric parameter
In geometric parameter (41) module, import elevation angle, the zenith angle of aircraft, sun altitude, zenith angle and observation date (MMDD) with respect to pixel;
Step 3-2, definition atmospherical model
The atmospherical model that defining mode adopted in atmospherical model (42) module comprises 6 kinds of atmospherical models such as tropical atmosphere, middle latitude summer, middle latitude winter, subpolar zone summer, subpolar zone winter, 1976 United States standard atmospheres;
Step 3-3, definition aerosol type pattern
Delustring and meteorological range type in the used boundary layer gasoloid of defining mode in aerosol type pattern (43) module (0~2km the is highly interval) pattern if defined visibility simultaneously, then replace the meteorological range of default value definition in the aerosol type.Aerosol type mainly comprises following several selection: (1) no gasoloid; (2) rural extinction coefficient, default meteorological range=23km; (3) city extinction coefficient, default meteorological range=5km; (4) troposphere extinction coefficient, default meteorological range=50km.
Step 3-4, definition aerosol concentration pattern
The used aerosol concentration pattern of defining mode in aerosol concentration (44), aerosol optical depth or the horizontal visibility (km) of input 550nm;
Step 3-5, input atural object elevation (km)
Step 3-6, input aircraft flight height (km)
Step 4, atmospheric correction parameter extraction
According to each input parameter in the step 3, the dark target on the aviation high-spectrum remote sensing of selecting is carried out atmospheric correction, from the output result, extract correction parameter aerosol scattering phase function (51) and gasoloid single scattering reflectivity (52);
Step 5, computation bound layer aerosol optical depth
According to formula 1., aerosol scattering phase function (51) and gasoloid single scattering reflectivity (052) by step 4 draws calculate boundary layer aerosol optical depth (60);
Step 6,
With synchronous ground sight distance data according to formula
V=3.91·H·1/τ
Note: V is ground sight distance (m); H is the gasoloid absolute altitude (Shanghai winter get 776.4m) of Various Seasonal; τ is an aerosol optical depth.The aerosol optical depth (06) that aerosol optical depth value that calculates and step 5 obtain compares (07), when both difference<0.1, just carries out 7, output (06) aerosol optical depth; Otherwise carry out 8, get back to step 3 and carry out atmospheric correction again, re-enter the needed parameters of atmospheric correction.
The specific embodiment of atmospheric correction of the present invention is as follows:
See also shown in the accompanying drawing 2,3,4,5,6, Shanghai practical modularization imaging spectrometer OMIS high spectrum image on October 7th, 2002 is chosen dark target, according to techniqueflow chart of the present invention, the inverting aerosol optical depth, at first the DN value of practical modularization imaging spectrometer OMIS high spectrum image is converted to radiation value through step 1, therefrom select dark surface (corresponding to polluted-water on Huangpu River) according to the size of radiation value, Fig. 2 is a sketch shown in the dark surface among the embodiment; Fig. 3 be dark surface respectively at 13 and 14 s' apparent reflectance curve, the x axle is wavelength (nanometer), the y axle is an apparent reflectance; Then, from ground Object Spectra DataBase, select corresponding object spectrum reflectivity, progressively import each parameter (3-1: solar zenith angle 41.8-44.8 degree, position angle 212 degree-219.6 degree according to step 3; 3-2: middle latitude atmospherical model in winter; 3-3: aerosols from major cities pattern; 3-4: visibility is 8km; 3-5: atural object elevation 0.004km; 3-6: aircraft flight height 2km), calculate atmospheric transmittance, Fig. 4 is that the atmospheric transmittance that calculates is with the wavelength change curve, dotted line is 13 an atmospheric transmittance curve, solid line is 14 an atmospheric transmittance curve, and the x axle is wavelength (nanometer), and the y axle is an atmospheric transmittance; Then inquire into according to formula 3 according to step 6 again and aerosol optical depth, Fig. 5 is exactly the aerosol optical depth that atmospheric correction method according to the present invention is finally inversed by from airborne practical modularization imaging spectrometer OMIS high-spectrum remote sensing, black line is represented 13 aerosol optical depth value, red line is represented 14 aerosol optical depth, the x axle is wavelength (nanometer), and the y axle is the aerosol optical depth value.
Table-1: practical modularization imaging spectrometer 0MIS spectrometer important technological parameters Table-1 Main parameters of the OMIS instrument
Type | Spectral range (μ m) | Sample interval/port number | The wave band number | Total visual field | Instantaneous field of view | Row pixel number | Sweep speed (line/second) | Digital coding |
??OMIS-I | ??0.46~1. ??1 ??1.06~1. ??7 ??2.0~2.5 ??3~5 ??8~12.5 | ??10nm/64 ? ??40nm/16 ? ??15nm/32 ??250nm/8 ??500nm/8 | ??64 ? ??16 ? ??32 ??8 ??8 | ??>70° ? ? ? | ??3mrad | ??512 | ??5 ??10 ??15 ??20 | ??12b ??it |
GPS bearing accuracy gyrostabilized platform speed height ratio system signal | 20 meters lasting accuracies are better than ± 4 '≤0.216 radian per second (speed unit: km/hour, height unit: rice≤300 |
(continued on next page)
Find that according to data-searching the method for carrying out the detailed complete of aerosol optical depth inverting at airborne high light spectrum image-forming data (especially OMIS imaging spectrometer) is not arranged so far.For this reason, we in conjunction with boundary layer atmosphere radiation transmission principle, extract aerosol optical depth information at practical modularization imaging spectrometer OMIS imaging spectrometer from atmospheric transmission spectrum.
Table-2: practical modularization imaging spectrometer OMIS inverting aerosol optical depth:
According to Shanghai Environment Monitoring Center's air quality historical summary query display: Shanghai City, 2002-10-713 point left and right sides air quality situation is totally better, sulphuric dioxide (SO
2) concentration is 0.051mg/m
3, nitrogen dioxide (NO
2) concentration is 0.044mg/m
3, oxides of nitrogen (NOx) concentration is 0.051mg/m
3, pellet (PM
10) concentration is 0.124mg/m
3, air pollution index (API) is in good level, and air quality is equivalent to that " ambient air quality (GB3095-1996) II level level, major pollutants are pellet PM
10, a little less than the pollution.The aerosol optical depth that inverting obtains numerically is an acceptable.
Practical modularization imaging spectrometer OMIS high spectrum image data on October 7th, 2002, in conjunction with the atmosphere radiation transmission equation, by computation bound atmosphere transmitance, from the atmospheric transmission spectrum of properties, inquire into aerosol optical depth, provided preliminary inversion result, in the aerosol optical depth value at 502nm-590nm wave band place between 0.175-0.314.Atmospheric transmittance under the rural gasoloid pattern (23km) that inversion result and the aerosol optical depth result (table-3) who calculates according to meteorological optical range and Fig. 6 are calculated by the LOWTRAN radiation delivery, aerosols from major cities (5km), the aerosols from major cities (8km) compares, and consistance is preferably arranged.
The aerosol optical depth that table 3 calculates according to the visibility data
Table?3?Calculated?AOD?from?visibility
Constantly | Pudong's instantaneous sight distance (m) of standing | Pudong's ten minutes average sight distances (m) in station | The instantaneous sight distance in station, Xuhui (m) | The ten minutes average sight distances (m) in station, Xuhui | The instantaneous sight distance in station, park, Huangpu (m) | The ten minutes average sight distances (m) in station, park, Huangpu | Fuel corporation's instantaneous sight distance (m) of standing | Fuel corporation's ten minutes average sight distances (m) in station | Aod mean value (550 nm) |
??13 ??:00 | ??1105 ??0 | ??9897 | ??8536 | ??9573 | ??14917 | ??15242 | ??10809 | ??10368 | ??0.27 ??9 |
??13 ??:10 | ??9648 | ??9397 | ??1028 ??3 | ??9273 | ??12899 | ??15504 | ??12267 | ??11102 | ??0.27 ??8 |
??13 ??:20 | ??9253 | ??9678 | ??8217 | ??8865 | ??14058 | ??16321 | ??7145 | ??10122 | ??0.31 ??0 |
??13 ??:30 | ??1232 ??0 | ??9585 | ??9610 | ??8751 | ??14471 | ??13295 | ??10306 | ??9892 | ??0.28 ??3 |
??13 ??:40 | ??9735 | ??9442 | ??8313 | ??8376 | ??13846 | ??15522 | ??9743 | ??9909 | ??0.29 ??9 |
??13 ??:50 | ??8556 | ??9775 | ??8663 | ??8485 | ??11368 | ??13140 | ??9096 | ??10231 | ??0.31 ??3 |
??14 ??:00 | ??8693 | ??8684 | ??7390 | ??7752 | ??11925 | ??12961 | ??9931 | ??9850 | ??0.32 ??5 |
??0.29 ??8 |
Claims (3)
1, a kind of atmospheric correction method of aviation high-spectrum remote-sensing inversion boundary layer aerosol optical depth, it is characterized in that: this method is by the atmosphere radiation transmission principle, from atmospheric transmission spectrum, extract boundary layer aerosol optical depth information, by being carried out once more interpolation in the boundary layer, atmospherical model, gasoloid pattern obtain showing reflectivity, realize the calculating of inversion boundary layer aerosol optical depth, its concrete calculation procedure is:
Step 1: the radiation calibration of high-spectrum remote sensing
A), read high-spectrum remote sensing (10)
Read airborne aviation high-spectrum remote sensing (10) by standard format;
B), convert radiation value (20) to
The output signal of high-spectrum remote sensing (10) module is sent to the input end of radiation value (20) module, slope and two coefficients of intercept according to each wave band correspondence in the radiation calibration file, according to formula radiation value=DN value * slope+intercept, convert the digital number DN value of high-spectrum remote sensing (10) to radiation value (20);
Step 2, calculating apparent reflectance are chosen dark target
A), calculate
Radiation value (20) with step 1 output is calculated as the apparent reflectance that obtains on the sensor by following formula:
R=π*L/(μ*f)
In the formula: R is an apparent reflectance;
L is a radiation value;
μ is the cosine of solar zenith angle;
F is an aeropause solar radiation flux density;
B), choose dark target
The output signal of radiation value (20) module is sent to dark target apparent reflectance (30) module, chooses the apparent reflectance on the near infrared 2.1um passage, and the pixel in 0.036~0.044 scope is dark target;
Step 3, atmospheric correction (40)
The output signal of dark target apparent reflectance (30) module is sent to atmospheric correction (40) module, and the dark target apparent reflectance (30) that obtains according to step 2 carries out atmospheric correction (40); Its concrete job step is as follows:
A), definition geometric parameter (41)
The input aircraft is with respect to elevation angle, the zenith angle of pixel in geometric parameter (41) module, and sun altitude, zenith angle and observation date, the output signal of its geometric parameter (41) module is sent to atmospheric correction (40) module;
B), definition atmospherical model (42)
The atmospherical model that defining mode adopted in atmospherical model (42) module, the output signal of its atmospherical model (42) module is sent to atmospheric correction (40) module;
C), definition aerosol type (43) pattern
Delustring and meteorological range type in aerosol type (43) in the used boundary layer gasoloid pattern of defining mode, if defined visibility simultaneously, then replace the meteorological range of default value definition in the aerosol type, the output signal of aerosol type (43) module is sent to atmospheric correction (40) module;
D), definition aerosol concentration (44) pattern
The used aerosol concentration pattern of defining mode in aerosol concentration (44), aerosol optical depth or the horizontal visibility of input 550nm, the output signal of aerosol concentration (44) module is sent to atmospheric correction (40) module;
E), input atural object elevation (45)
The output signal that to import atural object elevation (45) module is sent to atmospheric correction (40) module;
F), input aircraft flight height (46)
The output signal that to import aircraft flight height (46) module is sent to atmospheric correction (40) module;
Step 4, atmospheric correction parameter extraction
A), signal transmits
The output signal of atmospheric correction (40) module is sent to aerosol scattering phase function (51) module and gasoloid single scattering reflectivity (52) module respectively;
B), carry out atmospheric correction
According to each input parameter in the step 3, the dark target on the aviation high-spectrum remote sensing of selecting is carried out atmospheric correction;
C), correction parameter
From the output result, extract correction parameter aerosol scattering phase function (51) and gasoloid single scattering reflectivity (52);
Step 5, computation bound layer aerosol optical depth (60)
A), signal transmits
The output signal of aerosol scattering phase function (51) module and gasoloid single scattering reflectivity (52) module is sent to aerosol optical depth (60) module;
B), calculate
Apparent reflectance by following formula computation bound atmosphere:
C), computation bound layer aerosol optical depth (60)
Aerosol scattering phase function (51) and gasoloid single scattering reflectivity (52) by step 4 draws calculate boundary layer aerosol optical depth (60);
Step 6, judgement
A), signal transmits
The output signal of aerosol optical depth (60) module is sent to comparison (70) module;
B), calculate the ground sight distance
With synchronous ground sight distance data according to formula
V=3.91·H·1/τ
In the formula: V is the ground sight distance;
H is the gasoloid absolute altitude of Various Seasonal;
τ is an aerosol optical depth;
C), relatively
The aerosol optical depth (60) that the aerosol optical depth value that calculates and step 5 are obtained compares (70);
When both difference<0.1, the output signal that will compare (70) module is sent to output (80) module, is output as aerosol optical depth (60);
Otherwise feed back to atmospheric correction (40) module, get back to step 3 and carry out atmospheric correction again, re-enter the needed parameters of atmospheric correction, circulation repeats.
2, the atmospheric correction method of aviation high-spectrum remote-sensing inversion boundary layer aerosol optical depth according to claim 1, it is characterized in that: described atmospherical model (42) comprising: tropical atmosphere, middle latitude summer, middle latitude winter, subpolar zone summer, subpolar zone winter and 1976 United States standard atmospheres.
3, the atmospheric correction method of aviation high-spectrum remote-sensing inversion boundary layer aerosol optical depth according to claim 1, it is characterized in that: described aerosol type (43) comprising: or be no gasoloid; Or be rural extinction coefficient, default meteorological range=23km; Or be the city extinction coefficient, default meteorological range=5km; Or be the troposphere extinction coefficient, default meteorological range=50km.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100275240A CN100362318C (en) | 2005-07-05 | 2005-07-05 | Atmosphere correction method of airosol optical thickness of aeronautical high-spectrum remote-sensing inversion boundary layer |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100275240A CN100362318C (en) | 2005-07-05 | 2005-07-05 | Atmosphere correction method of airosol optical thickness of aeronautical high-spectrum remote-sensing inversion boundary layer |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1710379A true CN1710379A (en) | 2005-12-21 |
CN100362318C CN100362318C (en) | 2008-01-16 |
Family
ID=35706652
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005100275240A Expired - Fee Related CN100362318C (en) | 2005-07-05 | 2005-07-05 | Atmosphere correction method of airosol optical thickness of aeronautical high-spectrum remote-sensing inversion boundary layer |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100362318C (en) |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329173B (en) * | 2008-07-07 | 2010-06-02 | 武汉大学 | Method for correcting feculent water body atmosphere |
CN101936881A (en) * | 2010-09-15 | 2011-01-05 | 吉林大学 | Tomography method of using limb sounding data for inverting atmospheric ozone profiles |
CN102288956A (en) * | 2011-05-10 | 2011-12-21 | 中国资源卫星应用中心 | Atmospheric correction method for multispectral data of remote sensing satellite |
CN102346070A (en) * | 2010-07-30 | 2012-02-08 | 中国科学院遥感应用研究所 | HJ-1A satellite hyper-spectral imaging instrument orbit radiation scaling method |
CN101515987B (en) * | 2008-12-30 | 2012-04-25 | 中国资源卫星应用中心 | Method for radiometric correction of remote sensing image taken by rotary scan multiple parallel-scan infrared camera |
CN101477036B (en) * | 2009-01-13 | 2012-06-06 | 中国科学院遥感应用研究所 | Inland water chlorophyll a concentration remote-sensing monitoring method based on segmenting cooperation model |
CN102495405A (en) * | 2011-11-30 | 2012-06-13 | 武汉大学 | Evaluation method of TM/ETM (thematic mapper/enhanced thematic mapper) and image-based atmospheric correction product quality |
CN102628940A (en) * | 2012-04-20 | 2012-08-08 | 中国科学院遥感应用研究所 | Remote sensing image atmospheric correction method |
CN101699509B (en) * | 2009-11-11 | 2012-10-03 | 耿则勋 | Method for recovering atmosphere fuzzy remote image with meteorological data |
CN102955154A (en) * | 2012-10-16 | 2013-03-06 | 中国科学院遥感应用研究所 | High-resolution remote sensing data atmospheric correction method |
CN102955878A (en) * | 2012-09-05 | 2013-03-06 | 环境保护部卫星环境应用中心 | Inland water body optical classification method based on medium-resolution imaging spectrometer (MERIS) full-resolution image data |
CN103135113A (en) * | 2013-02-05 | 2013-06-05 | 中国科学院大气物理研究所 | Method for measuring height of atmospheric boundary layer |
CN103267533A (en) * | 2013-05-08 | 2013-08-28 | 中国科学院光电研究院 | Practical high spectrum remote sensing image atmosphere self-correction method and apparatus |
CN103605123A (en) * | 2013-12-04 | 2014-02-26 | 中国科学院遥感与数字地球研究所 | Parameterization remote sensing method based on oxygen A channel aerosol scattering effect |
CN103674815A (en) * | 2013-12-12 | 2014-03-26 | 中国科学院遥感与数字地球研究所 | Aerosol MTF (modulation transfer function) monitoring method based on sensor characteristics |
CN103674794A (en) * | 2013-12-16 | 2014-03-26 | 中国科学院遥感与数字地球研究所 | Multivariable regression method for remote sensing monitoring of near-surface fine particle matter PM2.5 mass concentration |
CN105136746A (en) * | 2015-07-17 | 2015-12-09 | 南京农业大学 | Multispectral crop growth sensor spectral reflectivity calibration method |
CN105953921A (en) * | 2016-04-15 | 2016-09-21 | 北京航空航天大学 | Quick simulation method for earth observation radiation image on condition of different aerosol parameters |
CN106096247A (en) * | 2016-06-06 | 2016-11-09 | 华北电力大学(保定) | Determination of Aerosol Optical method of estimation based on multiple-factor model of fit |
CN106096246A (en) * | 2016-06-06 | 2016-11-09 | 华北电力大学(保定) | Aerosol optical depth method of estimation based on PM2.5 and PM10 |
CN106407487A (en) * | 2015-07-27 | 2017-02-15 | 中国科学院遥感与数字地球研究所 | Method and system for assessing influence of aerosol scattering on CO2 remote sensing detection precision |
CN106569226A (en) * | 2016-10-28 | 2017-04-19 | 兰州大学 | Method using laser radar data to invert boundary layer height under cloudy condition |
CN107132545A (en) * | 2017-04-28 | 2017-09-05 | 北京理工大学 | The analysis method that particulate matter in aerosol influences on the quality of remotely sensed image |
CN108387547A (en) * | 2017-12-29 | 2018-08-10 | 华中科技大学 | A kind of survey spectrum compensation method of infared spectrum association remote sensing equipment |
CN111579504A (en) * | 2020-06-29 | 2020-08-25 | 中国科学技术大学 | Atmospheric pollution component vertical distribution inversion method based on optical remote sensing |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5471299A (en) * | 1992-02-21 | 1995-11-28 | The Secretary Of State For Defence In Her Britannic Majesty's Government Of The United Kingdom Of Great Britain And Northern Ireland | Apparatus and method for the analysis of particle characteristics using monotonically scattered light |
JP2004133302A (en) * | 2002-10-11 | 2004-04-30 | Toshiba Corp | Method for adhering photomask and apparatus therefor |
-
2005
- 2005-07-05 CN CNB2005100275240A patent/CN100362318C/en not_active Expired - Fee Related
Cited By (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101329173B (en) * | 2008-07-07 | 2010-06-02 | 武汉大学 | Method for correcting feculent water body atmosphere |
CN101515987B (en) * | 2008-12-30 | 2012-04-25 | 中国资源卫星应用中心 | Method for radiometric correction of remote sensing image taken by rotary scan multiple parallel-scan infrared camera |
CN101477036B (en) * | 2009-01-13 | 2012-06-06 | 中国科学院遥感应用研究所 | Inland water chlorophyll a concentration remote-sensing monitoring method based on segmenting cooperation model |
CN101699509B (en) * | 2009-11-11 | 2012-10-03 | 耿则勋 | Method for recovering atmosphere fuzzy remote image with meteorological data |
CN102346070A (en) * | 2010-07-30 | 2012-02-08 | 中国科学院遥感应用研究所 | HJ-1A satellite hyper-spectral imaging instrument orbit radiation scaling method |
CN101936881B (en) * | 2010-09-15 | 2012-10-03 | 吉林大学 | Tomography method of using limb sounding data for inverting atmospheric ozone profiles |
CN101936881A (en) * | 2010-09-15 | 2011-01-05 | 吉林大学 | Tomography method of using limb sounding data for inverting atmospheric ozone profiles |
CN102288956A (en) * | 2011-05-10 | 2011-12-21 | 中国资源卫星应用中心 | Atmospheric correction method for multispectral data of remote sensing satellite |
CN102495405A (en) * | 2011-11-30 | 2012-06-13 | 武汉大学 | Evaluation method of TM/ETM (thematic mapper/enhanced thematic mapper) and image-based atmospheric correction product quality |
CN102628940B (en) * | 2012-04-20 | 2014-04-16 | 中国科学院遥感应用研究所 | Remote sensing image atmospheric correction method |
CN102628940A (en) * | 2012-04-20 | 2012-08-08 | 中国科学院遥感应用研究所 | Remote sensing image atmospheric correction method |
CN102955878A (en) * | 2012-09-05 | 2013-03-06 | 环境保护部卫星环境应用中心 | Inland water body optical classification method based on medium-resolution imaging spectrometer (MERIS) full-resolution image data |
CN102955878B (en) * | 2012-09-05 | 2015-07-29 | 环境保护部卫星环境应用中心 | Based on the Inland Water optics sorting technique of MERIS full resolution image data |
CN102955154A (en) * | 2012-10-16 | 2013-03-06 | 中国科学院遥感应用研究所 | High-resolution remote sensing data atmospheric correction method |
CN102955154B (en) * | 2012-10-16 | 2014-04-16 | 中国科学院遥感应用研究所 | High-resolution remote sensing data atmospheric correction method |
CN103135113A (en) * | 2013-02-05 | 2013-06-05 | 中国科学院大气物理研究所 | Method for measuring height of atmospheric boundary layer |
CN103135113B (en) * | 2013-02-05 | 2014-08-06 | 中国科学院大气物理研究所 | Method for measuring height of atmospheric boundary layer |
CN103267533A (en) * | 2013-05-08 | 2013-08-28 | 中国科学院光电研究院 | Practical high spectrum remote sensing image atmosphere self-correction method and apparatus |
CN103267533B (en) * | 2013-05-08 | 2015-09-30 | 中国科学院光电研究院 | A kind of practical high-spectrum remote sensing air automatic correcting method |
CN103605123B (en) * | 2013-12-04 | 2016-08-31 | 中国科学院遥感与数字地球研究所 | Parametrization remote sensing technique based on oxygen A channel aerosol scattering effect |
CN103605123A (en) * | 2013-12-04 | 2014-02-26 | 中国科学院遥感与数字地球研究所 | Parameterization remote sensing method based on oxygen A channel aerosol scattering effect |
CN103674815A (en) * | 2013-12-12 | 2014-03-26 | 中国科学院遥感与数字地球研究所 | Aerosol MTF (modulation transfer function) monitoring method based on sensor characteristics |
CN103674815B (en) * | 2013-12-12 | 2016-03-23 | 中国科学院遥感与数字地球研究所 | A kind of gasoloid MTF monitoring method based on sensor characteristics |
CN103674794A (en) * | 2013-12-16 | 2014-03-26 | 中国科学院遥感与数字地球研究所 | Multivariable regression method for remote sensing monitoring of near-surface fine particle matter PM2.5 mass concentration |
CN103674794B (en) * | 2013-12-16 | 2016-06-01 | 中国科学院遥感与数字地球研究所 | Remote sensing monitoring near surface fine particle quality concentration PM2.5Multiple regression procedure |
CN105136746A (en) * | 2015-07-17 | 2015-12-09 | 南京农业大学 | Multispectral crop growth sensor spectral reflectivity calibration method |
CN105136746B (en) * | 2015-07-17 | 2018-11-30 | 南京农业大学 | A kind of multispectral plant growth sensor spectrum reflectivity correction method |
CN106407487A (en) * | 2015-07-27 | 2017-02-15 | 中国科学院遥感与数字地球研究所 | Method and system for assessing influence of aerosol scattering on CO2 remote sensing detection precision |
CN106407487B (en) * | 2015-07-27 | 2019-08-13 | 中国科学院遥感与数字地球研究所 | Aerosol scattering is assessed to CO2The method and system that remote sensing precision influences |
CN105953921B (en) * | 2016-04-15 | 2018-11-06 | 北京航空航天大学 | The rapid simulation method of earth observation radiation image under the conditions of aerosol parameters difference |
CN105953921A (en) * | 2016-04-15 | 2016-09-21 | 北京航空航天大学 | Quick simulation method for earth observation radiation image on condition of different aerosol parameters |
CN106096247A (en) * | 2016-06-06 | 2016-11-09 | 华北电力大学(保定) | Determination of Aerosol Optical method of estimation based on multiple-factor model of fit |
CN106096247B (en) * | 2016-06-06 | 2018-08-07 | 华北电力大学(保定) | Determination of Aerosol Optical method of estimation based on multiple-factor model of fit |
CN106096246B (en) * | 2016-06-06 | 2018-08-31 | 华北电力大学(保定) | Aerosol optical depth method of estimation based on PM2.5 and PM10 |
CN106096246A (en) * | 2016-06-06 | 2016-11-09 | 华北电力大学(保定) | Aerosol optical depth method of estimation based on PM2.5 and PM10 |
CN106569226A (en) * | 2016-10-28 | 2017-04-19 | 兰州大学 | Method using laser radar data to invert boundary layer height under cloudy condition |
CN106569226B (en) * | 2016-10-28 | 2019-02-05 | 兰州大学 | A method of have and utilizes laser radar Data Inversion Boundary Layer Height in the case of cloud |
CN107132545B (en) * | 2017-04-28 | 2019-05-28 | 北京理工大学 | The analysis method that particulate matter in aerosol influences the quality of remotely sensed image |
CN107132545A (en) * | 2017-04-28 | 2017-09-05 | 北京理工大学 | The analysis method that particulate matter in aerosol influences on the quality of remotely sensed image |
CN108387547A (en) * | 2017-12-29 | 2018-08-10 | 华中科技大学 | A kind of survey spectrum compensation method of infared spectrum association remote sensing equipment |
CN108387547B (en) * | 2017-12-29 | 2019-09-27 | 华中科技大学 | A kind of survey spectrum compensation method of infared spectrum association remote sensing equipment |
CN111579504A (en) * | 2020-06-29 | 2020-08-25 | 中国科学技术大学 | Atmospheric pollution component vertical distribution inversion method based on optical remote sensing |
CN111579504B (en) * | 2020-06-29 | 2021-10-01 | 中国科学技术大学 | Atmospheric pollution component vertical distribution inversion method based on optical remote sensing |
Also Published As
Publication number | Publication date |
---|---|
CN100362318C (en) | 2008-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1710379A (en) | Atmosphere correction method of airosol optical thickness of aeronautical high-spectrum remote-sensing inversion boundary layer | |
Ouyang et al. | Albedo changes caused by future urbanization contribute to global warming | |
Huang et al. | Characteristics of carbonaceous aerosol in PM2. 5: Pearl Delta River region, China | |
CN1690727A (en) | Method of measuring ultraviolet radiation and ultraviolet measuring device | |
Wei et al. | Evaluation and uncertainty estimate of next-generation geostationary meteorological Himawari-8/AHI aerosol products | |
Yin et al. | Exploring the effects of crop residue burning on local haze pollution in Northeast China using ground and satellite data | |
Kim et al. | An analysis on the impact of large-scale transports of dust pollution on air quality in East Asia as observed in central Korea in 2014 | |
CN1774645A (en) | Correction of troposphere induced errors in global positioning systems | |
Saha et al. | One year measurements of aerosol optical properties over an urban coastal site: Effect on local direct radiative forcing | |
CN102539336A (en) | Method and system for estimating inhalable particles based on HJ-1 satellite | |
Golitsyn et al. | Studying the pollution of Moscow and Beijing atmospheres with carbon monoxide and aerosol | |
CN101034475A (en) | Method for generating non-shadow satellite remote sensing orthographic digital image by computer | |
Dumka et al. | Altitude variation of aerosol properties over the Himalayan range inferred from spatial measurements | |
Qi et al. | Seasonal variation of aerosol optical properties in an urban site of the Yangtze Delta Region of China | |
CN101034472A (en) | Landform transformation of satellite remote sensing digital image supported by GIS | |
CN1556381A (en) | Aviation high spectrum remote sensing flight ground synchronous scaling and reflectivity conversion method | |
CN112798482B (en) | PM2.5 and PM10 estimation method based on satellite remote sensing | |
Du et al. | Photometric measurements of spring aerosol optical properties in dust and non-dust periods in China | |
Babu et al. | Seasonal changes in aerosol characteristics over Arabian Sea and their consequence on aerosol short-wave radiative forcing: Results from ARMEX field campaign | |
Guleria et al. | Validation of space-born Moderate Resolution Imaging Spectroradiometer remote sensors aerosol products using application of ground-based Multi-wavelength Radiometer | |
Geng et al. | Preliminary study of vertical distribution of aerosols during dry haze periods around Shanghai based on CALIPSO | |
Das et al. | Fog-induced changes in optical and physical properties of transported aerosols over Sundarban, India | |
LI et al. | Diagnosis of nitrogen nutrition of winter wheat and summer corn using images from digital camera equipped on unmanned aerial vehicle | |
Yuan-Jian et al. | Impacts of agricultural fire on aerosol distribution over East China during summer harvest time | |
Liao et al. | Mixing-layer-height-referenced ozone vertical distribution in the lower troposphere of Chinese megacities: stratification, classification, and meteorological and photochemical mechanisms |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080116 Termination date: 20100705 |