CN110705089B - Fine-mode aerosol parameter inversion method - Google Patents
Fine-mode aerosol parameter inversion method Download PDFInfo
- Publication number
- CN110705089B CN110705089B CN201910921707.9A CN201910921707A CN110705089B CN 110705089 B CN110705089 B CN 110705089B CN 201910921707 A CN201910921707 A CN 201910921707A CN 110705089 B CN110705089 B CN 110705089B
- Authority
- CN
- China
- Prior art keywords
- aerosol
- fine
- data
- inversion
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000000443 aerosol Substances 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000003287 optical effect Effects 0.000 claims abstract description 42
- 238000002310 reflectometry Methods 0.000 claims abstract description 34
- 230000010287 polarization Effects 0.000 claims abstract description 27
- 230000005540 biological transmission Effects 0.000 claims description 9
- 230000005855 radiation Effects 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 8
- 239000010419 fine particle Substances 0.000 claims description 6
- 238000002834 transmittance Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 4
- 238000012795 verification Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 241000282414 Homo sapiens Species 0.000 description 2
- 239000005427 atmospheric aerosol Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000013618 particulate matter Substances 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
-
- G01N15/075—
Abstract
The invention discloses a fine-mode aerosol parameter inversion method, which utilizes the characteristic that the reflectivity of a polarized earth surface is basically not changed along with a wave band from visible light to near infrared wave band; the characteristic of small change of the reflectivity of the polarized earth surface is combined; and (3) inverting the polarized remote sensing data of the POLDER sensor to obtain the optical thickness of the fine-mode aerosol. The method can effectively invert the fine-mode aerosol parameters and the polarization earth surface reflectivity at the same time, and the inverted fine-mode aerosol optical thickness has higher spatial coverage, higher spatial resolution and higher precision. The method can be used for small-scale fine-mode aerosol research, and the coverage and precision of bright earth surfaces (cities, deserts and the like) are increased.
Description
Technical Field
The invention relates to the field of atmospheric pollution treatment, in particular to an atmospheric aerosol observation technology.
Background
Atmospheric aerosols, which typically have a particle size of from a few nanometers to several tens of micrometers, are generally large in size for better evaluation of their global radiation balance, climate change, atmospheric environment, and human healthSmall size separates aerosols into fine mode and coarse mode aerosols. The radius of the fine mode aerosol is generally 0.1-0.2 microns, and the main source is artificial discharge. The radius of the coarse mode aerosol is generally larger than 1 micron, and mainly comes from natural emission. According to the fifth assessment report by the inter-government commission on climate change (IPCC), fine mode aerosols have a greater impact on climate change due to their large absorption properties. In addition, the fine mode aerosol is used for estimating the surface fine particulate matter PM2.5Is an important parameter. However, there is a great uncertainty in its quantitative estimation due to the high temporal and spatial variability characteristics of the aerosol. Particularly, as the number and activities of human beings are continuously increased, the emission of fine-modal aerosol is increased and the influence is increased. Therefore, it is increasingly important to obtain fine-modal aerosol products with high temporal resolution, high spatial resolution, and high precision, to reduce their own uncertainty and impact on other factors. However, currently, the inversion surface of the fine-mode aerosol still faces numerous challenges.
The fine-mode aerosol inversion mainly comprises two types of ground and satellite observation, the ground inversion has the characteristic of high precision, is generally regarded as a true value and can be used for verifying products inverted by the satellite, but the space coverage is insufficient, and the space distribution characteristics of the fine-mode aerosol cannot be comprehensively reflected. On the contrary, the method for inverting the fine-mode aerosol by using the satellite remote sensing data has the characteristics of wide observation range, less terrain constraint, short acquisition time and the like. However, due to the complexity of the land surface type, the space coverage, the space resolution and the precision of the land space fine mode aerosol product obtained by satellite remote sensing data have certain problems. For example, a single angle intensity sensor mode Resolution Imaging Spectroradiometer (MODIS) is used for inverting the acquired fine mode aerosol product, and related verification studies show that the fine mode aerosol product has relatively low precision and the official party does not recommend the fine mode aerosol product to be used. The reason is that the intensity signal at a single angle is not sufficient to distinguish fine mode aerosol from coarse mode aerosol information, and the surface signal accounts for a relatively large proportion of the total signal. In contrast, the polarized signal is only sensitive to the fine-mode aerosol within a certain scattering angle range, and the proportion of the polarized surface signal is smaller than or equal to that of the atmospheric signal, which are very beneficial for inverting the fine-mode aerosol by using the polarized signal. However, at present, only the Earth reflection POLarization measurement instrument POLDER (POLARIZATION and Direction of the Earth's reflexes) can provide long-time-series multi-angle polarized satellite remote sensing data, and Deuze et al utilize a surface POLarization reflectivity model (BPDF) to invert fine-mode aerosol products, and the algorithm is also a POLDER official algorithm. The released fine particle aerosol product is widely applied to various related fields, but the product has low spatial resolution and cannot be applied to small-scale aerosol research, and secondly, the product has serious deficiency conditions in cities, deserts, ice and snow and other bright earth surfaces and water bodies and other areas, and has relatively low precision. A large number of scholars have inverted terrestrial space-fine modal aerosols based on the polar data, but are basically similar to the official algorithm, i.e. using BPDF to estimate the polarized earth surface reflectance. Since the BPDF is a statistical model and the polarization earth surface reflectivity cannot be well estimated on some specific earth surfaces, the problems of low spatial resolution, product loss on specific earth surfaces, low precision and the like of the currently inverted fine-mode aerosol exist, and therefore, a new technical scheme for inverting the fine-mode aerosol is urgently needed to be developed.
Disclosure of Invention
The invention aims to provide a fine mode aerosol parameter inversion method which has higher spatial resolution, wider spatial coverage and higher precision of the obtained fine mode aerosol optical thickness than a PLODER official secondary product; the specific technical scheme is as follows:
the method comprises the following steps:
1) importing polarization remote sensing data of a POLDER sensor, selecting data of apparent polarization reflectivity of at least two wavelengths, and extracting a solar zenith angle, a solar azimuth angle, an observation zenith angle, a relative azimuth angle, a cloud mask, an amphibious mask and an elevation from a data file to be used as data to be processed;
cosΘ=-cosθscosθv-sinθssinθvcosφ (1)
2) shielding data of a cloud layer and a region covered by a water surface by using cloud mask data and land and water mask data, and calculating a scattering angle of each inversion data unit through a solar zenith angle, an observation zenith angle and a relative azimuth angle according to a formula (1); taking only data with scattering angles ranging from 80 ° to 120 ° as inversion data;
wherein theta is the scattering angle thetasAt the zenith angle of the sun, thetavPhi is a relative azimuth angle for observing a zenith angle;
3) processing the inversion data by using a lookup table for forward simulation of a 6SV radiation transmission model by taking an image element in an wxw window as an inversion unit to obtain a preprocessing lookup table;
4) performing linear interpolation in the preprocessing lookup table according to the observation geometric angle corresponding to each effective pixel in the inversion unit to obtain process radiation, upward atmospheric transmittance and downward atmospheric transmittance;
5) carrying out atmospheric correction on each effective pixel by using a formula (2) in combination with the polarization apparent reflectivity of each effective pixel to obtain the polarization earth surface reflectivity of different wave bands;
in the formula, lambda is the wave band,in order to polarize the apparent reflectance of light,in order to be the range radiation,for the downward atmospheric transmission rate,the upward air transmission rate is higher than that of the air,is the polarized surface reflectance;
6) calculating the ratio of the polarized earth surface reflectivities of the two wave bands according to a formula (3);
in the formula, λ1Is a first band, λ2In the second wavelength band, the first wavelength band,the ratio of the polarized earth surface reflectivity of the first wave band and the second wave band is obtained;
7) calculating the error of the ratio of the polarized earth surface reflectivity and the square root of 1 according to a formula (4), and selecting the effective fine-mode aerosol optical thickness and the polarized earth surface reflectivity when each aerosol model under an effective scattering angle enables the difference value of the polarized earth surface reflectivities of the two wave bands to be minimum;
8) respectively calculating the average values of the effective fine-mode aerosol optical thickness and the polarized earth surface reflectivity of the effective pixel according to formulas (5) and (6);
in formula (5), N1 is the effective fine mode aerosol optical thicknessNumber of (d), δmean_wThe average value of the optical thickness of the effective fine particle aerosol of the effective pixel;
in the formula (6), N2 represents the number of effective polarized surface reflectances,the effective polarization earth surface reflectivity average value of the second wave band of the effective pixel;
9) calculating a cost function eta of the inversion window unit according to formula (7), and selecting the value of eta which is the minimumAnd its corresponding aerosol optical thickness delta under aerosol modelmean_wRadius REObtaining an inversion result;
where STD represents the variance, w x w represents the selected window size,is the polarized earth surface reflectivity of 865nm wave band in the effective pixel.
Further, the wavelengths of the first and second bands are two of 490nm, 670nm and 865nm, respectively.
Further, the wavelength of the first waveband is 670 nm; the wavelength of the second waveband is 865 nm.
Further, the spatial resolution of the inversion is made to be consistent with the sensor resolution.
Further, the value of w is 3.
The method can effectively invert the fine-mode aerosol parameters and the polarization earth surface reflectivity at the same time, and the inverted fine-mode aerosol optical thickness has higher spatial coverage, higher spatial resolution and higher precision. The method can be used for small-scale fine-mode aerosol research, and the coverage and precision of bright earth surfaces (cities, deserts and the like) are increased.
Drawings
FIG. 1 is a comparison graph of the fine-mode aerosol parameter inversion method of embodiment 1 of the present invention and the effect of the prior art;
FIG. 2 is a comparison graph of the fine mode aerosol parameter inversion method of embodiment 2 of the present invention and the effect of the prior art;
FIG. 3 is a comparison of the fine mode aerosol parameter inversion method of embodiment 3 of the present invention and the prior art effect in FIG. 1;
fig. 4 is a comparison between the fine-mode aerosol parameter inversion method of embodiment 3 of the present invention and the prior art, and fig. 2.
In FIGS. 1-2: (a) the method comprises the following steps Satellite image true color picture; (b) the method comprises the following steps The method of the invention produces a result of the optical thickness of the fine-mode aerosol; (c) the method comprises the following steps Official fine mode aerosol optical thickness results.
Detailed Description
The present invention will now be more fully described with reference to the following examples. This invention may, however, be embodied in many different forms and should not be construed as limited to the exemplary embodiments set forth herein.
For ease of description, spatially relative terms, such as "upper," "lower," "left," "right," and the like, may be used herein to describe one element or feature's relationship to another element or feature as illustrated in the figures. It will be understood that the spatial terms are intended to encompass different orientations of the device in use or operation in addition to the orientation depicted in the figures. For example, if the device in the figures is turned over, elements described as "below" other elements or features would then be oriented "above" the other elements or features. Thus, the exemplary term "lower" can encompass both an upper and a lower orientation. The device may be otherwise oriented (rotated 90 degrees or at other orientations) and the spatially relative descriptors used herein interpreted accordingly.
In order to realize the invention, the invention discloses a high-spatial-resolution fine-modal aerosol parameter inversion method based on polarization data, which is applied to polarization remote sensing data of a POLDER sensor, and the specific implementation scheme is as follows:
in order to realize the invention, the method is applied to polarization remote sensing data of a POLDER sensor, wherein the POLDER sensor comprises 9 wave bands (443, 490, 565, 670, 763, 765, 865, 910 and 1020nm), and the 490, 670 and 865nm 3 wave bands comprise polarization channels. The POLDER L1 level data is a binary file and comprises a header file and a data file, wherein the header file comprises information of sensors and data processing, and the data file comprises sensor observation and auxiliary data information. The sensor observation information comprises multi-angle intensity observation data of 9 wave bands and multi-angle polarization observation data of 3 wave bands. The auxiliary data comprise data of a solar zenith angle, an observation zenith angle, a relative azimuth angle, a cloud mask, a land-water mask, an elevation and the like. Aiming at POLDER sensor data, the data used by the invention comprise 670 and 865nm multi-angle polarization observation data, and solar zenith angle, observation zenith angle, relative azimuth angle, cloud mask, land-water mask and elevation auxiliary data. The specific embodiment of the invention is as follows:
(1) simulating different wave bands, different observation geometries, different elevations, different aerosol models and different aerosol optical thicknesses (AOD) according to a forward simulation lookup table of a 6SV radiation transmission modelf) The specific parameters of the lookup table of parameter combinations are shown in table one.
Table one, lookup table parameter settings
(2) The method takes a 3 multiplied by 3 window as a unit, and input data comprise 670 and 865nm effective multi-angle apparent polarization reflectivity observed by each pixel in the window, and a corresponding solar zenith angle, an observation zenith angle, a relative azimuth angle, a cloud mask, an elevation and an amphibious mask. Data in the 470nm band may also be combined with data in the 670nm or 865nm band.
(3) The input data is quality controlled using a sea-land mask, a cloud mask and a scattering angle Θ, which is calculated by formula (1). Only image elements on land and clear sky and satisfying Θ between 80 ° and 120 ° can be used for fine mode aerosol parametric inversion.
cosΘ=-cosθscosθv-sinθssinθvcosφ (1)
Wherein theta is the scattering angle thetasIs a Chinese character ofZenith angle of the sun, thetavTo observe zenith angles, phi is the relative azimuth angle.
(4) Carrying out linear interpolation on different observation geometric angles of the effective pixel in a 3 multiplied by 3 window in a lookup table to obtain the scattering angle theta of the effective pixel_wCorresponding different aerosol models, different fine particle aerosol optical thicknesses delta_wAtmospheric parameters including range radiationUpward atmospheric transmittanceAnd downward atmospheric transmittanceAnd combined with apparent reflectivity of polarization for each effective pixel elementAtmosphere correction is carried out according to the formula (2), and further different aerosol optical thicknesses delta under different aerosol types of the effective pixel w are obtained_wDifferent scattering angle theta_w670 and 865nm band of polarized earth surface reflectivityand
In the formula, lambda is the wave band,in order to polarize the apparent reflectance of light,in order to be the range radiation,for the downward atmospheric transmission rate,the upward air transmission rate is higher than that of the air,is the polarized surface reflectance.
(5) Calculating different scattering angles theta under the effective pixel w according to the formula (3)_wDifferent aerosol types, different aerosol optical thicknesses δ_wPolarized surface reflectance ratio R of lower 670 and 865nm wave bandsratio。
In the formula, λ1Is 670nm band, lambda2Is a wave band of 865nm,the surface reflectance ratio is 670nm and 865nm wave band polarization.
(6) According to the characteristic that the reflectivity of the polarized earth surface does not change along with the wave band from visible light to near infrared wave band, different scattering angles theta under different effective pixels are calculated according to the formula (4)_wDifferent aerosol types, different aerosol optical thicknesses δ_wLower 670 and 865nm polarized earth surface reflectivity ratio Rratio_wAnd 1 square root error.
(7) For each effective aerosol model, for each effective scattering angle Θ_wBy selecting the fine mode aerosol optical thickness delta_wLet DiffR _wOf minimum value, i.e. polarized surface reflectanceAndthe closest; i.e. the optical thickness delta of the fine mode aerosol_wA fine mode aerosol optical thickness as a corresponding effective aerosol type at an effective scattering angle for the pixel. Therefore, the optimal fine mode aerosol parameters and the optimal polarization earth surface reflectivity of each effective aerosol model under each effective scattering angle of each effective pixel are obtained.
(8) Respectively calculating the fine mode aerosol optical thickness delta of each effective aerosol model under each effective pixel according to the formulas (5) and (6) aiming at each effective aerosol model of each effective pixel_wAverage value of (2) and polarized surface reflectanceIncluding the mean value delta of the optical thickness of the fine-particle aerosolmean_wAnd 865nm band of polarized earth surface reflectivity
Wherein N1 is the effective fine mode aerosol optical thicknessNumber, deltamean_wThe average value of the optical thickness of the aerosol of the fine particles of the effective pixel is shown.
Wherein N2 is effective polarization earth surface reflectivity of 865nm wave bandThe number of the first and second groups is,is the average value of the surface reflectivity of the effective pixel 865nm wave band.
(9) Calculating a cost function eta according to the formula (7), selecting the one which minimizes the value of etaAnd its corresponding aerosol optical thickness delta under aerosol modelmean_wRadius REThe space resolution can reach about 6km per pixel, namely the inversion result; with higher resolution sensors (e.g., DPC sensors, resolution of 3km per pixel), the resolution of the inversion can also be achieved for single pixel inversion.
In the formula, STD represents variance, 3 × 3 represents selected window size, and 4x4 or 5x5 window size can be selected;is the polarized earth surface reflectivity of 865nm wave band of the effective pixel. The inversion efficiency is best compared to the inversion results with a 3x3 window.
In order to qualitatively evaluate the spatial distribution characteristics of the inverted fine-mode aerosol optical thickness, the spatial distribution characteristics are compared with the POLDER official algorithm result.
Example 1
Fig. 1(a) is a true color image of the polar satellite image at 11/15/2011 in eastern china, (b) is the optical thickness result of the fine mode aerosol obtained by the method of the present invention, and (c) is the optical thickness result of the official fine mode aerosol. The land surface type of the area is mainly vegetation, cultivated land and city, and as can be seen from the figure, the official inversion resolution is about 18km per pixel; the optical thickness of the fine-mode aerosol inverted by the method is consistent with the spatial distribution of official results on the whole, but the method has higher spatial coverage, and particularly compensates the condition that a large amount of products are lost in the area of an official algorithm in coastal areas with gathered population and developed economy. In addition, the optical thickness spatial resolution of the fine-mode aerosol is higher, the detail change is more obvious, and the research of the small and medium-scale fine-mode aerosol is more facilitated.
Example 2
Fig. 2(a) is a true color image of the polar satellite image at 11/18/2011 in northern west of china and india, (b) is the optical thickness result of the fine mode aerosol of the present invention, and (c) is the optical thickness result of the official fine mode aerosol. The land surface type of the region mainly comprises desert, mountain land and cultivated land, and as is obvious from the figure, the official inversion resolution ratio is about 18km per pixel; the space coverage of the inverted fine-modal aerosol optical thickness in desert, mountain and coastal areas is obviously higher than that of an official algorithm, the condition that the official algorithm has a large number of product defects in the surface type area is made up, and particularly in Qinghai-Tibet plateau areas with complex terrain and high aerosol inversion difficulty is solved. In addition, the fine mode aerosol optical thickness product of the algorithm successfully captures the condition of heavy pollution on the plain of the indian constant river, which is low in the Qinghai-Tibet plateau, and is more obvious in detail change.
Example 3
In order to quantitatively evaluate the precision of the inverted fine-mode aerosol optical thickness, foundation AERONET site results respectively located in different types of tables are selected for verification. Fig. 3 is a result of verification between the fine mode aerosol optical thickness inverted by the method and the official algorithm in 2010-2012 based on the POLDER data and the ground-based site, respectively, and it can be seen from the result that the fine mode aerosol optical thickness inverted by the method is better than the official algorithm result no matter the inverted number N, the correlation R, the root mean square error RMSE or the ratio Gfrac falling within the error range, which indicates that the algorithm is obviously better than the official algorithm. Fig. 4 is a verification result of the method and the official algorithm of the invention based on the POLDER data in 2010-2012 and 3 respectively with an SACOL foundation site, the SACOL site is located in the northwest semiarid area in China, the vegetation coverage is low, the inversion difficulty of the fine-mode aerosol is high, from the verification result, the optical thickness result of the fine-mode aerosol inverted by the method is obviously more than the official result, and the inversion accuracy is closer to the observation result of the foundation, which indicates that the algorithm of the invention has higher applicability in the semiarid area.
In general, whether the qualitative comparison with the result of the official algorithm or the quantitative verification of the result of the foundation AERONET site shows that the algorithm has higher spatial coverage, higher spatial resolution and higher precision than the fine-modal aerosol optical thickness inverted by the official algorithm, and the defects of the official algorithm product are overcome.
The above examples are only for illustrating the present invention, and besides, there are many different embodiments, which can be conceived by those skilled in the art after understanding the idea of the present invention, and therefore, they are not listed here.
Claims (4)
1. A method for inverting parameters of a fine-mode aerosol is characterized by comprising the following steps:
1) importing polarization remote sensing data of a POLDER sensor, selecting data of at least two apparent polarization reflectances with wavelengths in a range from visible light to near infrared, and extracting a solar zenith angle, an observation zenith angle, a relative azimuth angle, a cloud mask, a land-water mask and an elevation from a data file as data to be processed;
2) shielding data of a cloud layer and a region covered by a water surface by using cloud mask data and land and water mask data, and calculating a scattering angle of each inversion data unit through a solar zenith angle, an observation zenith angle and a relative azimuth angle according to a formula (1); taking only data with scattering angles ranging between 80 and 120 ° as inversion data;
cosΘ=-cosθscosθv-sinθssinθvcosφ (1)
wherein theta is the scattering angle thetasAt the zenith angle of the sun, thetavPhi is a relative azimuth angle for observing a zenith angle;
3) processing the inversion data by using a forward simulation lookup table of a 6SV radiation transmission model by taking an wxw pixel as an inversion unit to obtain a preprocessing lookup table;
4) performing linear interpolation in the preprocessing lookup table according to the observation geometric angle corresponding to each effective pixel in the inversion unit to obtain process radiation, upward atmospheric transmittance and downward atmospheric transmittance;
5) carrying out atmospheric correction on each effective pixel according to a formula (2) by combining the polarization apparent reflectivity of each effective pixel to obtain the polarization earth surface reflectivity of different wave bands;
in the formula, lambda is the wave band,in order to polarize the apparent reflectance of light,in order to be the range radiation,for the downward atmospheric transmission rate,the upward air transmission rate is higher than that of the air,is the polarized surface reflectance;
6) calculating the ratio of the polarized earth surface reflectivities of the two wave bands according to a formula (3);
in the formula, λ1Is a first band, λ2In the second wavelength band, Rratioδ_wPolarising the earth surface for a first band and a second bandA ratio of refractive indices;
7) calculating the error of the square root of the ratio of the reflectivity of the polarized earth surface and 1 according to the formula (4);
8) selecting the optical thickness of the fine mode aerosol with the minimum difference of the polarized earth surface reflectivities of the two wave bands from each aerosol model corresponding to each effective scattering angle according to the scattering angle
In formula (5), N1 is the effective fine mode aerosol optical thicknessNumber, deltamean_wThe average value of the optical thickness of the effective fine particle aerosol of the effective pixel;
in the formula (6), N2 represents the number of effective polarized surface reflectances,the effective polarization earth surface reflectivity average value of the second wave band of the effective pixel;
9) calculating a cost function eta of the window according to the formula (7), and selecting the one which minimizes the eta valueAnd its corresponding aerosol optical thickness delta under aerosol modelmean_wRadius REAs a result of the inversion;
2. The method of fine mode aerosol parametric inversion of claim 1, wherein the wavelengths of the first and second bands are two of 490nm, 670nm and 865nm, respectively.
3. The method of fine mode aerosol parametric inversion of claim 2, wherein the wavelengths of the first and second bands are 670nm and 865nm, respectively.
4. The method of fine mode aerosol parametric inversion of claim 1, wherein w has a value of 3.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910921707.9A CN110705089B (en) | 2019-09-27 | 2019-09-27 | Fine-mode aerosol parameter inversion method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910921707.9A CN110705089B (en) | 2019-09-27 | 2019-09-27 | Fine-mode aerosol parameter inversion method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110705089A CN110705089A (en) | 2020-01-17 |
CN110705089B true CN110705089B (en) | 2021-04-16 |
Family
ID=69197687
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910921707.9A Active CN110705089B (en) | 2019-09-27 | 2019-09-27 | Fine-mode aerosol parameter inversion method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110705089B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116660106B (en) * | 2023-07-21 | 2023-10-17 | 中国科学院空天信息创新研究院 | Aerosol parameter iterative inversion method for collaborative satellite-borne scalar and polarization observation data |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7337065B2 (en) * | 2001-01-23 | 2008-02-26 | Spectral Sciences, Inc. | Methods for atmospheric correction of solar-wavelength hyperspectral imagery over land |
WO2009108795A1 (en) * | 2008-02-26 | 2009-09-03 | Battelle Memorial Institute | Biological and chemical microscopic targeting |
CN103927455A (en) * | 2014-04-24 | 2014-07-16 | 中国科学院遥感与数字地球研究所 | Land aerosol optical property retrieval method based on Gaofen-1 satellite |
CN106225693A (en) * | 2016-08-29 | 2016-12-14 | 中国科学院遥感与数字地球研究所 | A kind of fine particle aerosol optical thickness and aerosol type Simultaneous Inversion method |
CN109781593A (en) * | 2019-01-17 | 2019-05-21 | 南京泛在地理信息产业研究院有限公司 | A kind of aerosol quadratic inversion method based on PARASOL multi-angle polarization data |
CN110186823A (en) * | 2019-06-26 | 2019-08-30 | 中国科学院遥感与数字地球研究所 | A kind of aerosol optical depth inversion method |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8345252B2 (en) * | 2005-04-25 | 2013-01-01 | X-Rite, Inc. | Method and system for enhanced formulation and visualization rendering |
US8558884B2 (en) * | 2011-05-04 | 2013-10-15 | Raytheon Company | In-scene determination of aerosol parameters from imagery |
CN110186822B (en) * | 2019-05-13 | 2020-09-11 | 中国科学院遥感与数字地球研究所 | Aerosol optical thickness remote sensing inversion method |
-
2019
- 2019-09-27 CN CN201910921707.9A patent/CN110705089B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7337065B2 (en) * | 2001-01-23 | 2008-02-26 | Spectral Sciences, Inc. | Methods for atmospheric correction of solar-wavelength hyperspectral imagery over land |
WO2009108795A1 (en) * | 2008-02-26 | 2009-09-03 | Battelle Memorial Institute | Biological and chemical microscopic targeting |
CN103927455A (en) * | 2014-04-24 | 2014-07-16 | 中国科学院遥感与数字地球研究所 | Land aerosol optical property retrieval method based on Gaofen-1 satellite |
CN106225693A (en) * | 2016-08-29 | 2016-12-14 | 中国科学院遥感与数字地球研究所 | A kind of fine particle aerosol optical thickness and aerosol type Simultaneous Inversion method |
CN109781593A (en) * | 2019-01-17 | 2019-05-21 | 南京泛在地理信息产业研究院有限公司 | A kind of aerosol quadratic inversion method based on PARASOL multi-angle polarization data |
CN110186823A (en) * | 2019-06-26 | 2019-08-30 | 中国科学院遥感与数字地球研究所 | A kind of aerosol optical depth inversion method |
Non-Patent Citations (4)
Title |
---|
Directional Polarimetric Camera (DPC): Monitoring aerosol spectral optical properties over land from satellite observation;Zhengqiang Li,et al.;《Journal of Quantitative Spectroscopy & Radiative Transfer》;20180707;21-37 * |
Retrieval of aerosol microphysical and optical properties over land using a multimode approach;Guangliang Fu,et al.;《Atmospheric Measurement Techniques》;20181217;6627-6650 * |
基于多角度偏振信息反演海洋上空卷云云顶高度;李树等;《红外与毫米波学报》;20180831;第37卷(第4期);445-458 * |
基于多角度标量和偏振卫星数据的气溶胶光学参数反演研究;张洋;《中国博士学位论文全文数据库 工程科技I辑》;20190415;B027-29 * |
Also Published As
Publication number | Publication date |
---|---|
CN110705089A (en) | 2020-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ghonima et al. | A method for cloud detection and opacity classification based on ground based sky imagery | |
Sirguey et al. | Subpixel monitoring of the seasonal snow cover with MODIS at 250 m spatial resolution in the Southern Alps of New Zealand: Methodology and accuracy assessment | |
CN101963664B (en) | Microwave remote sensing pixel element decomposing method based on land and water living beings classifying information | |
Li et al. | An evaluation of the use of atmospheric and BRDF correction to standardize Landsat data | |
Liang et al. | Global LAnd Surface Satellite (GLASS) products: algorithms, validation and analysis | |
Zhang et al. | A Landsat 8 OLI-based, semianalytical model for estimating the total suspended matter concentration in the slightly turbid Xin’anjiang Reservoir (China) | |
Guo et al. | Atmospheric correction comparison of SPOT-5 image based on model FLAASH and model QUAC | |
König et al. | Application of sentinel-2 MSI in Arctic research: Evaluating the performance of atmospheric correction approaches over Arctic sea ice | |
Cui et al. | Assessment of atmospheric correction methods for historical Landsat TM images in the coastal zone: A case study in Jiangsu, China | |
Mishra et al. | Assessment of different topographic corrections in AWiFS satellite imagery of Himalaya terrain | |
Qin et al. | A geometry-dependent surface Lambertian-equivalent reflectivity product for UV–Vis retrievals–Part 1: Evaluation over land surfaces using measurements from OMI at 466 nm | |
Wang et al. | An adaptive atmospheric correction algorithm for the effective adjacency effect correction of submeter-scale spatial resolution optical satellite images: Application to a WorldView-3 panchromatic image | |
Goyens et al. | High angular resolution measurements of the anisotropy of reflectance of sea ice and snow | |
CN108764326A (en) | Urban impervious surface extracting method based on depth confidence network | |
Zhang et al. | A multiband model with successive projections algorithm for bathymetry estimation based on remotely sensed hyperspectral data in Qinghai lake | |
Bair et al. | Divergence of apparent and intrinsic snow albedo over a season at a sub-alpine site with implications for remote sensing | |
CN110705089B (en) | Fine-mode aerosol parameter inversion method | |
CN111650128B (en) | High-resolution atmospheric aerosol inversion method based on surface reflectivity library | |
CN116664947A (en) | Blue algae bloom monitoring method and system based on satellite observation data | |
CN105259145A (en) | Method for simultaneous remote sensing of underwater terrain and features of island | |
Sartika et al. | Determining the Precision of Spectral Patterns Arising from Atmospheric Correction Utilizing MODTRAN-FLAASH and 6S Approaches on High-Resolution SPOT-6 Imagery | |
Gao et al. | Estimation of surface sediment moisture content in muddy tidal flats using analytical radiative transfer model | |
Yin et al. | Quantitative typical land cover remote sensing and its application in earthquake evaluation | |
Chen et al. | Atmospheric correction of remote sensing imagery based on the surface spectrum’s vector space | |
Balick et al. | LANL MTI science team experience |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |