CN107202750A - A kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates - Google Patents
A kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates Download PDFInfo
- Publication number
- CN107202750A CN107202750A CN201710349064.6A CN201710349064A CN107202750A CN 107202750 A CN107202750 A CN 107202750A CN 201710349064 A CN201710349064 A CN 201710349064A CN 107202750 A CN107202750 A CN 107202750A
- Authority
- CN
- China
- Prior art keywords
- gray haze
- ground
- mrow
- remote sensing
- satellite
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000012545 processing Methods 0.000 title claims abstract description 23
- 238000005267 amalgamation Methods 0.000 title claims abstract description 18
- 239000000443 aerosol Substances 0.000 claims abstract description 51
- 239000002245 particle Substances 0.000 claims abstract description 39
- 238000009826 distribution Methods 0.000 claims abstract description 28
- 230000003287 optical effect Effects 0.000 claims abstract description 23
- 239000000926 atmospheric chemistry Substances 0.000 claims abstract description 6
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 6
- 231100000719 pollutant Toxicity 0.000 claims abstract description 6
- 230000008033 biological extinction Effects 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 claims description 9
- 239000005427 atmospheric aerosol Substances 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 239000012528 membrane Substances 0.000 claims description 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 230000008901 benefit Effects 0.000 abstract description 3
- 230000007613 environmental effect Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 239000000203 mixture Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000011109 contamination Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000007499 fusion processing Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
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
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates, mainly include step a:Obtain gray haze aerosol particle subtype or component parsing;Step b:Remote sensing estimation model is improved with priori based on Beijing-tianjin-hebei Region gray haze aerosol particle sub-feature observation data, obtained gray haze optical thickness inverting PM2.5 concentration distributions near the ground;Step c:Gray haze optical thickness inverting PM2.5 concentration distributions near the ground, pollutant ground observation data and meteorological, environment auxiliary information, and gray haze aerosol particle subtype or component parsing are combined, ground gray haze integrated data is formed;Step d:Gray haze integrated data in ground in step c and Atmospheric Chemistry or air quality model, and satellite remote sensing, ground observation multi-source data are combined, the short-period forecast of gray haze pollution distribution is realized.A kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates of the present invention, makes the gray haze data more high advantage of accurate, reliability of acquisition.
Description
Technical field
The present invention relates to environmental monitoring technology field, in particular it relates to which a kind of satellite-ground integrated monitoring of Atmospheric particulates is quantitative
RS fusion processing method.
Background technology
As science and technology is developed rapidly, our environment is also in increasingly exhaustion, and the health and environment of the mankind is closely bound up, protection
Environment is very urgent.Then the acquisition of China's environmental monitoring index, main effective acquisition using sensor utilizes various algorithms
Obtain the thematic product of fuse information of the high spatial resolution, high time resolution and high detection accuracy of environmental information.But, it is right
In the integration technology of environmental monitoring big data, also in developing stage, technology is also immature.In face of substantial amounts of data, such as where
The valid data that reason could obtain environmental information turn into the key of research.
At present, the acquisition for environmental information there is shortage of data, imperfect and data accuracy it is not high, insecure
Defect.
The content of the invention
It is an object of the present invention in view of the above-mentioned problems, propose a kind of satellite-ground integrated monitoring quantitative remote sensing of Atmospheric particulates
Method for amalgamation processing, to realize the total data of fast acquiring gray haze, make the last related data for obtaining gray haze it is more accurate, can
By the high advantage of property.
To achieve the above object, the technical solution adopted by the present invention is:A kind of satellite-ground integrated monitoring of Atmospheric particulates is quantitative
RS fusion processing method, methods described is specifically included:
Step a:Gray haze remote sensing recognition is realized based on multi-source, polymorphic type satellite remote sensing date, to obtain gray haze gas
Sol particles type or component parsing;
Step b:Based on Beijing-tianjin-hebei Region gray haze aerosol particle sub-feature observation data and heavily contaminated aerosol properties priori
Knowledge improves remote sensing estimation model, so as to realize the quantitative inversion of this area's gray haze optical thickness, obtained gray haze optics is thick
Spend inverting PM2.5 concentration distributions near the ground;
Step c:By gray haze optical thickness inverting PM2.5 concentration distributions near the ground, pollutant ground observation data with it is meteorological,
Environment auxiliary information, and the gray haze aerosol particle subtype or component parsing are combined, and form ground gray haze integrated data;
Step d:By gray haze integrated data in ground described in step c with being based on Atmospheric Chemistry or air quality model, Yi Jiwei
Star remote sensing, ground observation multi-source data are combined, and realize the short-period forecast of gray haze pollution distribution.
Further, gray haze aerosol particle subtype is obtained in the step a or component parsing is specifically included:
Assuming that the gray haze to some region of the Beijing-Tianjin wing carries out comprehensive monitoring, provided with n sensor, including domestic and international
High Resolution Remote Sensing Satellites and ground transaucer synchronize monitoring, and monitoring sample is 1,2,3 ... n,
Step a1:The gray haze in some region of Jing-jin-ji region is monitored using n sensor;
Step a2:The big data that monitoring is obtained is analyzed and processed using PCA, relevant information is obtained;
Step a3:The relevant information of acquisition is subjected to computer disposal, the contribution rate collection of illustrative plates of gray haze is obtained.
Further, the PCA methods in the step a2 are specifically included:
Step 1:Conversion is standardized to the sample that monitoring is obtained
Standardized acquisition p dimensions random vector x=(x1, X2 ..., Xp) T of original index data) n sample xi=
(xi1, xi2 ..., xip) T, i=1,2 ..., n, n > p, construction sample battle array, following standardized transformation is carried out to sample array element:
WhereinA gust Z must be standardized;
Step 2:To seeking correlation matrix to standardization battle array Z
Correlation matrix is asked to standardization battle array Z
Wherein,
Step 3:Solve sample correlation matrix R characteristic equation | R- λ Ip|=0 p characteristic root, determine principal component byM values are determined, make the utilization rate of information up to more than 85%, to each λ j, j=1,2 ..., m, solving equations
Rb=λ jb obtain unit character vector
Step 4:Target variable after standardization is converted into principal component
U1 is referred to as first principal component, and U2 is referred to as Second principal component, ..., Up is referred to as pth principal component;
Step 5:Overall merit is carried out to m principal component
Summation is weighted to m principal component, final evaluation of estimate is produced, flexible strategy are the variance contribution ratio of each principal component.
Further, computer disposal is carried out to relevant information in the step a3, it is main to use MATLAB software programmings
Relative program simultaneously carries out the processing of related data to obtain contribution rate collection of illustrative plates.
Further, the quantitative inversion model of this area's gray haze optical thickness is in the step b:
Assuming that aerosol is that concentration is different under any height or thickness, the ratio between component is identical, therefore disappears
Backscatter extinction logarithmic ratio simply changes with height or thickness.So, different height or the aerosol of concentration are exactly by this single
Equivalent particle composition, simply particle book is different.Equivalent particle number at height z is n (z), if the delustring of equivalent particle
Coefficient is q, and mass concentration is p, then is respectively for the extinction coefficient q (z) and mass concentration m (z) of aerosol at z in height:
Q (z)=n (z) * q,
M (z)=n (z) * p;
It can be obtained with reference to above formula,
M (z)=n (z) * p=q (z)/q*p=x*q;
Wherein x=p/q is proportionality coefficient, relevant with atmospheric aerosol yardstick Spectral structure, component and light refractive index;
So, aerosol is just directly proportional in the mass concentration of different height or thickness to the extinction coefficient at this.
Further, the inverse model of gray haze optical thickness inverting PM2.5 concentration distributions near the ground is obtained in the step b
Flow is:
Step b1:Obtain the mass concentration value at ground heavily contaminated Aerosol Extinction value moment corresponding with same place.
Step b2:The extinction coefficient under the same terms with mass concentration to putting together, pass through iterative method computation model
Parameter.
Step b3:The model parameter and Beijing-tianjin-hebei Region gray haze aerosol particle sub-feature obtained by calculating observes disappearing for data
Backscatter extinction logarithmic ratio, inverting PM2.5 concentration distributions near the ground.
Further, Aerosol Extinction can be obtained by laser radar in the step b1, and mass concentration can pass through filter
Film sampling is weighed or oscillating balance measurement is obtained.
A kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates of the present invention, method is specifically included:
Step a:Gray haze remote sensing recognition is realized based on multi-source, polymorphic type satellite remote sensing date, with obtain gray haze aerosol particle subtype or
Component is parsed;Step b:Based on Beijing-tianjin-hebei Region gray haze aerosol particle sub-feature observation data and heavily contaminated aerosol properties priori
Knowledge improves remote sensing estimation model, so as to realize the quantitative inversion of this area's gray haze optical thickness, obtained gray haze optics is thick
Spend inverting PM2.5 concentration distributions near the ground;Step c:By gray haze optical thickness inverting PM2.5 concentration distributions near the ground, pollutant
Ground observation data and meteorological, environment auxiliary information, and the gray haze aerosol particle subtype or component parsing are combined, shape
Into ground gray haze integrated data;Step d:By gray haze integrated data in ground described in step c with being based on Atmospheric Chemistry or air quality
Pattern, and satellite remote sensing, ground observation multi-source data are combined, and realize the short-period forecast of gray haze pollution distribution.
Mainly obtain following technique effect:
(1) the multi-source remote sensing automatic identification of gray haze
Lack the remote-sensing monitoring method for China's gray haze distribution in the world at present, gray haze is realized using multi- source Remote Sensing Data data
The blank in the domestic and international field will be effectively filled up in the accurate extraction of distribution, is the important wound to China's atmosphere pollution satellite monitoring
Newly.
(2) inverting of gray haze optical thickness
Current international mainstream aerosol optical inversion method is primarily directed to lighter clear sky condition is polluted, to east China
Weight haze pollution is inapplicable or can not inverting;This research will do important development, innovative realization to existing aerosol inversion method
To the quantitative inversion of Beijing-tianjin-hebei Region gray haze heavily contaminated.
(3) PM2.5 retrieving concentrations near the ground under the conditions of gray haze
Estimate that particle concentration near the ground is always international forward position and difficulties based on satellite remote sensing, estimate both at home and abroad at present
Calculation method, which is also primarily adapted for use in, pollutes lighter situation, and PM2.5 near the ground is realized for the extremely strong gray haze heavily contaminated condition of delustring
The estimation of concentration, not only makes important innovations to the field inversion method, while being also that Contamination Assessment and gray haze prediction are provided
Support.
(4) gray haze aerosol particle subtype or component parsing
Using n sensor of satellite-ground integrated monitoring, all data of Atmospheric particulates can be obtained, it is to avoid the something lost of data
Leakage;Data are analyzed using PCA, the maximal correlation data of gray haze can be effectively obtained, point of gray haze is more quickly obtained
Butut;The processing of contribution rate collection of illustrative plates is carried out using MATLAB softwares, makes data more precisely, reliably.
Other features and advantages of the present invention will be illustrated in the following description, also, partly becomes from specification
Obtain it is clear that or being understood by implementing the present invention.
Below by drawings and examples, technical scheme is described in further detail.
Brief description of the drawings
Accompanying drawing is used for providing a further understanding of the present invention, and constitutes a part for specification, the reality with the present invention
Applying example is used to explain the present invention together, is not construed as limiting the invention.In the accompanying drawings:
Fig. 1 is a kind of flow of the satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates of the present invention
Figure;
Fig. 2 is at a kind of PCA of the satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates of the present invention
Manage flow chart;
Fig. 3 is a kind of contribution rate of the satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates of the present invention
Collection of illustrative plates.
Embodiment
The preferred embodiments of the present invention are illustrated below in conjunction with accompanying drawing, it will be appreciated that preferred reality described herein
Apply example to be merely to illustrate and explain the present invention, be not intended to limit the present invention.
As shown in figure 1, a kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates, methods described is specific
Including:
Step a:Gray haze remote sensing recognition is realized based on multi-source, polymorphic type satellite remote sensing date, to obtain gray haze gas
Sol particles type or component parsing;
Step b:Based on Beijing-tianjin-hebei Region gray haze aerosol particle sub-feature observation data and heavily contaminated aerosol properties priori
Knowledge improves remote sensing estimation model, so as to realize the quantitative inversion of this area's gray haze optical thickness, obtained gray haze optics is thick
Spend inverting PM2.5 concentration distributions near the ground;
Step c:By gray haze optical thickness inverting PM2.5 concentration distributions near the ground, pollutant ground observation data with it is meteorological,
Environment auxiliary information, and the gray haze aerosol particle subtype or component parsing are combined, and form ground gray haze integrated data;
Step d:By gray haze integrated data in ground described in step c with being based on Atmospheric Chemistry or air quality model, Yi Jiwei
Star remote sensing, ground observation multi-source data are combined, and realize the short-period forecast of gray haze pollution distribution.
Gray haze aerosol particle subtype is obtained in the step a or component parsing is specifically included:
Assuming that the gray haze to some region of the Beijing-Tianjin wing carries out comprehensive monitoring, provided with n sensor, including domestic and international
High Resolution Remote Sensing Satellites and ground transaucer synchronize monitoring, and monitoring sample is 1,2,3 ... n,
Step a1:The gray haze in some region of Jing-jin-ji region is monitored using n sensor;
Step a2:The big data that monitoring is obtained is analyzed and processed using PCA, relevant information is obtained;
Step a3:The relevant information of acquisition is subjected to computer disposal, the contribution rate collection of illustrative plates of gray haze is obtained.
PCA methods in the step a2 are specifically included:
Step 1:Conversion is standardized to the sample that monitoring is obtained
Standardized acquisition p dimensions random vector x=(x1, X2 ..., Xp) T of original index data) n sample xi=
(xi1, xi2 ..., xip) T, i=1,2 ..., n, n > p, construction sample battle array, following standardized transformation is carried out to sample array element:
WhereinA gust Z must be standardized;
Step 2:To seeking correlation matrix to standardization battle array Z
Correlation matrix is asked to standardization battle array Z
Wherein,
Step 3:Solve sample correlation matrix R characteristic equation | R- λ Ip|=0 obtains p characteristic root, determines principal component
PressDetermine m values, make the utilization rate of information up to more than 85%, to each λ j, j=1,2 ...,
m,
Solving equations Rb=λ jb obtain unit character vector
Step 4:Target variable after standardization is converted into principal component
U1 is referred to as first principal component, and U2 is referred to as Second principal component, ..., Up is referred to as pth principal component;
Step 5:Overall merit is carried out to m principal component
Summation is weighted to m principal component, final evaluation of estimate is produced, flexible strategy are the variance contribution ratio of each principal component.
Computer disposal is carried out to relevant information in the step a3, it is main to use MATLAB software programmings relative program simultaneously
The processing of related data is carried out to obtain contribution rate collection of illustrative plates.
As shown in Fig. 2 PCA handling processes include:1st, data are read in;2nd, pre-process, standardization;3rd, principal component contribution is calculated
Rate;4th, contribution rate sorts, and chooses principal component;5th, each principal component load is calculated;6th, output pattern.
As shown in figure 3, in the Monitoring Data of two groups of gray hazes, the contribution rate of first principal component is significantly greater than other masters
The contribution rate of composition.
The quantitative inversion model of this area's gray haze optical thickness is in the step b:
Assuming that aerosol is that concentration is different under any height or thickness, the ratio between component is identical, therefore disappears
Backscatter extinction logarithmic ratio simply changes with height or thickness.So, different height or the aerosol of concentration are exactly by this single
Equivalent particle composition, simply particle book is different.Equivalent particle number at height z is n (z), if the delustring of equivalent particle
Coefficient is q, and mass concentration is p, then is respectively for the extinction coefficient q (z) and mass concentration m (z) of aerosol at z in height:
Q (z)=n (z) * q,
M (z)=n (z) * p;
It can be obtained with reference to above formula,
M (z)=n (z) * p=q (z)/q*p=x*q;
Wherein x=p/q is proportionality coefficient, relevant with atmospheric aerosol yardstick Spectral structure, component and light refractive index;
So, aerosol is just directly proportional in the mass concentration of different height or thickness to the extinction coefficient at this.
Further, the inverse model of gray haze optical thickness inverting PM2.5 concentration distributions near the ground is obtained in the step b
Flow is:
Step b1:Obtain the mass concentration value at ground heavily contaminated Aerosol Extinction value moment corresponding with same place.
Step b2:The extinction coefficient under the same terms with mass concentration to putting together, pass through iterative method computation model
Parameter.
Step b3:The model parameter and Beijing-tianjin-hebei Region gray haze aerosol particle sub-feature obtained by calculating observes disappearing for data
Backscatter extinction logarithmic ratio, inverting PM2.5 concentration distributions near the ground.
Aerosol Extinction can be obtained by laser radar in the step b1, and mass concentration can be weighed by filter membrane sampling
Or oscillating balance measurement is obtained.
Following beneficial effect can at least be reached:
A kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates of the present invention, method is specifically included:
Step a:Gray haze remote sensing recognition is realized based on multi-source, polymorphic type satellite remote sensing date, with obtain gray haze aerosol particle subtype or
Component is parsed;Step b:Based on Beijing-tianjin-hebei Region gray haze aerosol particle sub-feature observation data and heavily contaminated aerosol properties priori
Knowledge improves remote sensing estimation model, so as to realize the quantitative inversion of this area's gray haze optical thickness, obtained gray haze optics is thick
Spend inverting PM2.5 concentration distributions near the ground;Step c:By gray haze optical thickness inverting PM2.5 concentration distributions near the ground, pollutant
Ground observation data and meteorological, environment auxiliary information, and the gray haze aerosol particle subtype or component parsing are combined, shape
Into ground gray haze integrated data;Step d:By gray haze integrated data in ground described in step c with being based on Atmospheric Chemistry or air quality
Pattern, and satellite remote sensing, ground observation multi-source data are combined, and realize the short-period forecast of gray haze pollution distribution.
Mainly obtain following technique effect:
(1) the multi-source remote sensing automatic identification of gray haze
Lack the remote-sensing monitoring method for China's gray haze distribution in the world at present, gray haze is realized using multi- source Remote Sensing Data data
The blank in the domestic and international field will be effectively filled up in the accurate extraction of distribution, is the important wound to China's atmosphere pollution satellite monitoring
Newly.
(2) inverting of gray haze optical thickness
Current international mainstream aerosol optical inversion method is primarily directed to lighter clear sky condition is polluted, to east China
Weight haze pollution is inapplicable or can not inverting;This research will do important development, innovative realization to existing aerosol inversion method
To the quantitative inversion of Beijing-tianjin-hebei Region gray haze heavily contaminated.
(3) PM2.5 retrieving concentrations near the ground under the conditions of gray haze
Estimate that particle concentration near the ground is always international forward position and difficulties based on satellite remote sensing, estimate both at home and abroad at present
Calculation method, which is also primarily adapted for use in, pollutes lighter situation, and PM2.5 near the ground is realized for the extremely strong gray haze heavily contaminated condition of delustring
The estimation of concentration, not only makes important innovations to the field inversion method, while being also that Contamination Assessment and gray haze prediction are provided
Support.
(4) gray haze aerosol particle subtype or component parsing
Using n sensor of satellite-ground integrated monitoring, all data of Atmospheric particulates can be obtained, it is to avoid the something lost of data
Leakage;Data are analyzed using PCA, the maximal correlation data of gray haze can be effectively obtained, point of gray haze is more quickly obtained
Butut;The processing of contribution rate collection of illustrative plates is carried out using MATLAB softwares, makes data more precisely, reliably.
Finally it should be noted that:The preferred embodiments of the present invention are the foregoing is only, are not intended to limit the invention,
Although the present invention is described in detail with reference to the foregoing embodiments, for those skilled in the art, it still may be used
To be modified to the technical scheme described in foregoing embodiments, or equivalent substitution is carried out to which part technical characteristic.
Within the spirit and principles of the invention, any modification, equivalent substitution and improvements made etc., should be included in the present invention's
Within protection domain.
Claims (7)
1. a kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates, it is characterised in that methods described is specific
Including:
Step a:Gray haze remote sensing recognition is realized based on multi-source, polymorphic type satellite remote sensing date, to obtain gray haze aerosol particle subclass
Type or component parsing;
Step b:Based on Beijing-tianjin-hebei Region gray haze aerosol particle sub-feature observation data and heavily contaminated aerosol properties priori
To improve remote sensing estimation model, so as to realize the quantitative inversion of this area's gray haze optical thickness, obtained gray haze optical thickness is anti-
Drill PM2.5 concentration distributions near the ground;
Step c:By gray haze optical thickness inverting PM2.5 concentration distributions near the ground, pollutant ground observation data and meteorological, environment
Auxiliary information, and the gray haze aerosol particle subtype or component parsing are combined, and form ground gray haze integrated data;
Step d:It is by gray haze integrated data in ground described in step c and distant based on Atmospheric Chemistry or air quality model, and satellite
Sense, ground observation multi-source data are combined, and realize the short-period forecast of gray haze pollution distribution.
2. the satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of a kind of Atmospheric particulates according to claim 1, it is special
Levy and be, gray haze aerosol particle subtype is obtained in the step a or component parsing is specifically included:
Assuming that the gray haze to some region of the Beijing-Tianjin wing carries out comprehensive monitoring, provided with n sensor, including domestic and international high score
Resolution remote sensing satellite and ground transaucer synchronize monitoring, and monitoring sample is 1,2,3 ... n,
Step a1:The gray haze in some region of Jing-jin-ji region is monitored using n sensor;
Step a2:The big data that monitoring is obtained is analyzed and processed using PCA, relevant information is obtained;
Step a3:The relevant information of acquisition is subjected to computer disposal, the contribution rate collection of illustrative plates of gray haze is obtained.
3. the satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of a kind of Atmospheric particulates according to claim 2, it is special
Levy and be, the PCA methods in the step a2 are specifically included:
Step 1:Conversion is standardized to the sample that monitoring is obtained
Standardized acquisition p dimensions random vector x=(x1, X2 ..., Xp) T of original index data) n sample xi=(xi1,
Xi2 ..., xip) T, i=1,2 ..., n, n > p, construction sample battle array, to the following standardized transformation of sample array element progress:
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mi>j</mi>
</msub>
</mrow>
<msub>
<mi>s</mi>
<mi>j</mi>
</msub>
</mfrac>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>n</mi>
<mo>;</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>p</mi>
</mrow>
WhereinA gust Z must be standardized;
Step 2:To seeking correlation matrix to standardization battle array Z
Correlation matrix is asked to standardization battle array Z
<mrow>
<mi>R</mi>
<mo>=</mo>
<msub>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mi>p</mi>
</msub>
<mi>x</mi>
<mi>p</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>Z</mi>
<mi>T</mi>
</msup>
<mi>Z</mi>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
</mrow>
Wherein,
Step 3:Solve sample correlation matrix R characteristic equation | R- λ Ip|=0 obtains p characteristic root, determines principal component
PressM values are determined, make the utilization rate of information up to more than 85%, to each λ j, j=1,2 ..., m, solution
Equation group Rb=λ jb obtain unit character vector
Step 4:Target variable after standardization is converted into principal component
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>z</mi>
<mi>i</mi>
<mi>T</mi>
</msubsup>
<msubsup>
<mi>b</mi>
<mi>j</mi>
<mi>o</mi>
</msubsup>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>m</mi>
</mrow>
U1 is referred to as first principal component, and U2 is referred to as Second principal component, ..., Up is referred to as pth principal component;
Step 5:Overall merit is carried out to m principal component
Summation is weighted to m principal component, final evaluation of estimate is produced, flexible strategy are the variance contribution ratio of each principal component.
4. the satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of a kind of Atmospheric particulates according to claim 2, it is special
Levy and be, computer disposal is carried out to relevant information in the step a3, it is main to use MATLAB software programmings relative program simultaneously
The processing of related data is carried out to obtain contribution rate collection of illustrative plates.
5. the satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of a kind of Atmospheric particulates according to claim 1, it is special
Levy and be, the quantitative inversion model of this area's gray haze optical thickness is in the step b:
Assuming that aerosol is that concentration is different under any height or thickness, the ratio between component is identical, therefore delustring system
Number simply changes with height or thickness.So, different height or the aerosol of concentration are exactly by this single etc.
Imitate grain molecular, simply particle book is different.Equivalent particle number at height z is n (z), if the extinction coefficient of equivalent particle
For q, mass concentration is p, then is respectively for the extinction coefficient q (z) and mass concentration m (z) of aerosol at z in height:
Q (z)=n (z) * q,
M (z)=n (z) * p;
It can be obtained with reference to above formula,
M (z)=n (z) * p=q (z)/q*p=x*q;
Wherein x=p/q is proportionality coefficient, relevant with atmospheric aerosol yardstick Spectral structure, component and light refractive index;
So, aerosol is just directly proportional in the mass concentration of different height or thickness to the extinction coefficient at this.
6. the satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of a kind of Atmospheric particulates according to claim 1, it is special
Levy and be, the inverse model flow that gray haze optical thickness inverting PM2.5 concentration distributions near the ground are obtained in the step b is:
Step b1:Obtain the mass concentration value at ground heavily contaminated Aerosol Extinction value moment corresponding with same place.
Step b2:The extinction coefficient under the same terms with mass concentration to putting together, pass through the ginseng of iterative method computation model
Number.
Step b3:The model parameter and Beijing-tianjin-hebei Region gray haze aerosol particle sub-feature obtained by calculating observes the delustring system of data
Number, inverting PM2.5 concentration distributions near the ground.
7. the satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of a kind of Atmospheric particulates according to claim 6, it is special
Levy and be, Aerosol Extinction can be obtained by laser radar in the step b1, and mass concentration can be weighed by filter membrane sampling
Or oscillating balance measurement is obtained.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710349064.6A CN107202750A (en) | 2017-05-17 | 2017-05-17 | A kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710349064.6A CN107202750A (en) | 2017-05-17 | 2017-05-17 | A kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107202750A true CN107202750A (en) | 2017-09-26 |
Family
ID=59906549
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710349064.6A Pending CN107202750A (en) | 2017-05-17 | 2017-05-17 | A kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107202750A (en) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107766799A (en) * | 2017-09-28 | 2018-03-06 | 中国地质大学(武汉) | The analysis method and system of multi- source Remote Sensing Data data source remittance landscape based on scale effect |
CN108426815A (en) * | 2018-04-20 | 2018-08-21 | 中国科学院遥感与数字地球研究所 | A kind of fine particle concentration of component evaluation method near the ground |
CN108805368A (en) * | 2018-09-06 | 2018-11-13 | 天津珞雍空间信息研究院有限公司 | A kind of atmospheric environment intellectualized management system based on space-air-ground integration |
CN109213839A (en) * | 2018-09-12 | 2019-01-15 | 北京英视睿达科技有限公司 | Hot spot grid pollutant data capture method based on the study of gridding depths of features |
CN109543935A (en) * | 2018-10-09 | 2019-03-29 | 北京英视睿达科技有限公司 | Environmental data processing method based on hot spot grid |
CN109934377A (en) * | 2018-12-29 | 2019-06-25 | 段文旭 | A kind of the interaction factor control analysis of industrial pollution source synthesis and prediction technique of fine particle |
CN110261272A (en) * | 2019-07-05 | 2019-09-20 | 西南交通大学 | Based on geographical detection with PCA to the Key Influential Factors screening technique of PM2.5 concentration distribution |
CN110411918A (en) * | 2019-08-02 | 2019-11-05 | 中国科学院遥感与数字地球研究所 | A kind of PM2.5 concentration remote-sensing evaluation method based on satellite polarization technology |
CN110673229A (en) * | 2019-10-23 | 2020-01-10 | 新亚优华科技有限公司 | Atmospheric pollutant diffusion track tracking method based on hotspot grid technology |
CN110909309A (en) * | 2019-11-21 | 2020-03-24 | 中国科学院遥感与数字地球研究所 | Method for acquiring hourly high-resolution PM2.5 data |
CN110989038A (en) * | 2019-11-22 | 2020-04-10 | 北京中科锐景科技有限公司 | Near-surface particle concentration estimation method and system |
CN111220517A (en) * | 2019-12-06 | 2020-06-02 | 南京大学 | PM based on satellite remote sensing2.5Inversion method for components with different particle sizes |
CN111323352A (en) * | 2020-04-09 | 2020-06-23 | 中南大学 | Regional PM2.5 remote sensing inversion model fusing fine particulate matter concentration data |
CN111709180A (en) * | 2020-06-24 | 2020-09-25 | 中国科学院空天信息创新研究院 | Atmospheric vector radiation transmission simulation method of coupling absorptive aerosol model |
CN112068158A (en) * | 2020-09-08 | 2020-12-11 | 中国科学院空天信息创新研究院 | Black carbon aerosol aging state estimation method based on back scattering observation |
CN112819026A (en) * | 2020-12-11 | 2021-05-18 | 中科三清科技有限公司 | Fusion method and device for optical thickness data of multisource satellite aerosol |
CN113447623A (en) * | 2021-09-02 | 2021-09-28 | 航天宏图信息技术股份有限公司 | Atmospheric environment monitoring method and system |
CN114397276A (en) * | 2022-01-19 | 2022-04-26 | 北华航天工业学院 | Regional soil humidity monitoring method based on equivalent rainfall estimation method |
CN115201074A (en) * | 2022-06-30 | 2022-10-18 | 中国科学院大气物理研究所 | Method, system, device and computer-readable storage medium for remote sensing inversion of aerosol component distribution |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102944503A (en) * | 2012-10-11 | 2013-02-27 | 中国科学院安徽光学精密机械研究所 | PM2.5 mass concentration value automatic inversion algorithm based on sun photometer and lidar |
CN103674815A (en) * | 2013-12-12 | 2014-03-26 | 中国科学院遥感与数字地球研究所 | Aerosol MTF (modulation transfer function) monitoring method based on sensor characteristics |
CN103927454A (en) * | 2014-04-24 | 2014-07-16 | 中国科学院遥感与数字地球研究所 | Dust-haze pollution monitoring method based on environmental satellite |
KR101609740B1 (en) * | 2015-08-24 | 2016-04-08 | 대한민국 | Vertical distribution calculation method of Aerosol mass concentration |
CN106124374A (en) * | 2016-07-22 | 2016-11-16 | 中科宇图科技股份有限公司 | Atmospheric particulates remote-sensing monitoring method based on data fusion |
CN106407656A (en) * | 2016-08-29 | 2017-02-15 | 中国科学院遥感与数字地球研究所 | Retrieval method for aerosol optical thickness based on high resolution satellite image data |
CN106446307A (en) * | 2015-08-05 | 2017-02-22 | 中国科学院遥感与数字地球研究所 | Aerosol foundation data-based AOD vertical correction effect evaluation method and system |
-
2017
- 2017-05-17 CN CN201710349064.6A patent/CN107202750A/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102944503A (en) * | 2012-10-11 | 2013-02-27 | 中国科学院安徽光学精密机械研究所 | PM2.5 mass concentration value automatic inversion algorithm based on sun photometer and lidar |
CN103674815A (en) * | 2013-12-12 | 2014-03-26 | 中国科学院遥感与数字地球研究所 | Aerosol MTF (modulation transfer function) monitoring method based on sensor characteristics |
CN103927454A (en) * | 2014-04-24 | 2014-07-16 | 中国科学院遥感与数字地球研究所 | Dust-haze pollution monitoring method based on environmental satellite |
CN106446307A (en) * | 2015-08-05 | 2017-02-22 | 中国科学院遥感与数字地球研究所 | Aerosol foundation data-based AOD vertical correction effect evaluation method and system |
KR101609740B1 (en) * | 2015-08-24 | 2016-04-08 | 대한민국 | Vertical distribution calculation method of Aerosol mass concentration |
CN106124374A (en) * | 2016-07-22 | 2016-11-16 | 中科宇图科技股份有限公司 | Atmospheric particulates remote-sensing monitoring method based on data fusion |
CN106407656A (en) * | 2016-08-29 | 2017-02-15 | 中国科学院遥感与数字地球研究所 | Retrieval method for aerosol optical thickness based on high resolution satellite image data |
Non-Patent Citations (1)
Title |
---|
亢红霞等: "基于卫星遥感数据(AOD)估算PM2.5的研究进展", 《环境科学与管理》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107766799A (en) * | 2017-09-28 | 2018-03-06 | 中国地质大学(武汉) | The analysis method and system of multi- source Remote Sensing Data data source remittance landscape based on scale effect |
CN108426815A (en) * | 2018-04-20 | 2018-08-21 | 中国科学院遥感与数字地球研究所 | A kind of fine particle concentration of component evaluation method near the ground |
CN108426815B (en) * | 2018-04-20 | 2021-04-27 | 中国科学院遥感与数字地球研究所 | Method for estimating concentration of components of near-surface fine particulate matters |
CN108805368A (en) * | 2018-09-06 | 2018-11-13 | 天津珞雍空间信息研究院有限公司 | A kind of atmospheric environment intellectualized management system based on space-air-ground integration |
CN109213839A (en) * | 2018-09-12 | 2019-01-15 | 北京英视睿达科技有限公司 | Hot spot grid pollutant data capture method based on the study of gridding depths of features |
CN109543935A (en) * | 2018-10-09 | 2019-03-29 | 北京英视睿达科技有限公司 | Environmental data processing method based on hot spot grid |
CN109934377A (en) * | 2018-12-29 | 2019-06-25 | 段文旭 | A kind of the interaction factor control analysis of industrial pollution source synthesis and prediction technique of fine particle |
CN110261272A (en) * | 2019-07-05 | 2019-09-20 | 西南交通大学 | Based on geographical detection with PCA to the Key Influential Factors screening technique of PM2.5 concentration distribution |
CN110261272B (en) * | 2019-07-05 | 2020-08-18 | 西南交通大学 | Method for screening key influence factors on PM2.5 concentration distribution based on geographic detection and PCA (principal component analysis) |
CN110411918A (en) * | 2019-08-02 | 2019-11-05 | 中国科学院遥感与数字地球研究所 | A kind of PM2.5 concentration remote-sensing evaluation method based on satellite polarization technology |
CN110411918B (en) * | 2019-08-02 | 2020-09-08 | 中国科学院遥感与数字地球研究所 | PM2.5 concentration remote sensing estimation method based on satellite polarization technology |
CN110673229A (en) * | 2019-10-23 | 2020-01-10 | 新亚优华科技有限公司 | Atmospheric pollutant diffusion track tracking method based on hotspot grid technology |
CN110909309A (en) * | 2019-11-21 | 2020-03-24 | 中国科学院遥感与数字地球研究所 | Method for acquiring hourly high-resolution PM2.5 data |
CN110989038A (en) * | 2019-11-22 | 2020-04-10 | 北京中科锐景科技有限公司 | Near-surface particle concentration estimation method and system |
CN110989038B (en) * | 2019-11-22 | 2022-03-04 | 北京中科锐景科技有限公司 | Near-surface particle concentration estimation method and system |
CN111220517A (en) * | 2019-12-06 | 2020-06-02 | 南京大学 | PM based on satellite remote sensing2.5Inversion method for components with different particle sizes |
CN111323352A (en) * | 2020-04-09 | 2020-06-23 | 中南大学 | Regional PM2.5 remote sensing inversion model fusing fine particulate matter concentration data |
CN111323352B (en) * | 2020-04-09 | 2021-08-24 | 中南大学 | Regional PM2.5 remote sensing inversion model fusing fine particulate matter concentration data |
CN111709180A (en) * | 2020-06-24 | 2020-09-25 | 中国科学院空天信息创新研究院 | Atmospheric vector radiation transmission simulation method of coupling absorptive aerosol model |
CN111709180B (en) * | 2020-06-24 | 2023-10-13 | 中国科学院空天信息创新研究院 | Atmospheric vector radiation transmission simulation method of coupling absorptive aerosol model |
CN112068158A (en) * | 2020-09-08 | 2020-12-11 | 中国科学院空天信息创新研究院 | Black carbon aerosol aging state estimation method based on back scattering observation |
CN112819026A (en) * | 2020-12-11 | 2021-05-18 | 中科三清科技有限公司 | Fusion method and device for optical thickness data of multisource satellite aerosol |
CN112819026B (en) * | 2020-12-11 | 2022-02-15 | 中科三清科技有限公司 | Fusion method and device for optical thickness data of multisource satellite aerosol |
CN113447623A (en) * | 2021-09-02 | 2021-09-28 | 航天宏图信息技术股份有限公司 | Atmospheric environment monitoring method and system |
CN114397276A (en) * | 2022-01-19 | 2022-04-26 | 北华航天工业学院 | Regional soil humidity monitoring method based on equivalent rainfall estimation method |
CN114397276B (en) * | 2022-01-19 | 2023-12-01 | 北华航天工业学院 | Regional soil humidity monitoring method based on equivalent precipitation estimation method |
CN115201074A (en) * | 2022-06-30 | 2022-10-18 | 中国科学院大气物理研究所 | Method, system, device and computer-readable storage medium for remote sensing inversion of aerosol component distribution |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107202750A (en) | A kind of satellite-ground integrated monitoring quantitative remote sensing method for amalgamation processing of Atmospheric particulates | |
Smith et al. | Regional estimation of pollutant gas dry deposition in the UK: model description, sensitivity analyses and outputs | |
Chen et al. | MICS-Asia III: multi-model comparison and evaluation of aerosol over East Asia | |
Vignati et al. | Global scale emission and distribution of sea-spray aerosol: Sea-salt and organic enrichment | |
Bakwin et al. | Measurements of carbon dioxide on very tall towers: results of the NOAA/CMDL program | |
Pan et al. | Wet and dry deposition of atmospheric nitrogen at ten sites in Northern China | |
Velasco et al. | Cities as net sources of CO2: Review of atmospheric CO2 exchange in urban environments measured by eddy covariance technique | |
Yuan et al. | Identification and spatiotemporal variations of dominant PM10 sources over Hong Kong | |
Sarmiento et al. | A US carbon cycle science plan | |
Baldini et al. | An investigation of the controls on Irish precipitation δ 18 O values on monthly and event timescales | |
CN110058258B (en) | Atmospheric boundary layer detection method based on hybrid laser radar | |
Erisman et al. | Monitoring and modelling of biosphere/atmosphere exchange of gases and aerosols in Europe | |
CN106404620A (en) | Method and system for inverting ground PM 2.5 through combination of geostatistical interpolation and satellite remote sensing | |
Droste et al. | Assessing the potential and application of crowdsourced urban wind data | |
Goto et al. | Application of a global nonhydrostatic model with a stretched-grid system to regional aerosol simulations around Japan | |
CN110889196B (en) | Water environment bearing capacity assessment method and device based on water quality model and storage medium | |
Bai et al. | Comparison of satellite-based PM2. 5 estimation from aerosol optical depth and top-of-atmosphere reflectance | |
CN110595968A (en) | PM2.5 concentration estimation method based on geostationary orbit satellite | |
Erisman et al. | Deposition monitoring in Europe | |
Begum et al. | Identification of haze-creating sources from fine particulate matter in Dhaka aerosol using carbon fractions | |
CN113435068A (en) | Radionuclide assimilation prediction method based on logarithmic variational assimilation | |
Chen et al. | Significant influence of carbonates on determining organic carbon and black carbon: a case study in Tajikistan, central Asia | |
CN110907318B (en) | Near-ground atmospheric total suspended particulate matter mass concentration remote sensing physical estimation method | |
CN108614803A (en) | A kind of meteorological data method of quality control and system | |
Gan et al. | Prediction of air pollutant concentration based on luong attention mechanism Seq2Seq model |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170926 |
|
WD01 | Invention patent application deemed withdrawn after publication |