CN106950574B - The remote sensing measuring method and device of gray haze total amount in a kind of atmosphere - Google Patents
The remote sensing measuring method and device of gray haze total amount in a kind of atmosphere Download PDFInfo
- Publication number
- CN106950574B CN106950574B CN201710244393.4A CN201710244393A CN106950574B CN 106950574 B CN106950574 B CN 106950574B CN 201710244393 A CN201710244393 A CN 201710244393A CN 106950574 B CN106950574 B CN 106950574B
- Authority
- CN
- China
- Prior art keywords
- gray haze
- haze
- gray
- measured target
- region
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/95—Lidar systems specially adapted for specific applications for meteorological use
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The present invention provides a kind of remote sensing measuring method of gray haze total amount in atmosphere and devices, comprising: obtains the moonscope image of selected gray haze Polluted area；According to the corresponding measured target identical criterion of the different-waveband of satellite, the gray haze region in the moonscope image is identified；The measured target identical criterion includes at least gray haze identical criterion；According in the gray haze region gray haze intensity, gray haze region area, gray haze thickness and gray haze vertical distribution profiles, determine the gray haze total amount in selected gray haze Polluted area；It can be according to satellite in the corresponding measured target identical criterion of different-waveband, gray haze Polluted area is identified in real time and further calculates gray haze total amount, realize gray haze pollution layer vertical space structure detection, haze pollution flood detecting thickness and haze pollution spatial distribution monitoring, can Quantitative Monitoring take place frequently regional gray haze pollution and monitoring result accuracy it is higher.
Description
Technical field
The present invention relates to atmospheric remote sensing pollution detection technical fields, in particular to gray haze total amount in a kind of atmosphere
Remote sensing measuring method and device.
Background technique
Gray haze is that a large amount of particulate pollutants gather in an atmosphere, causes the stuffy phenomenon of On The Deterioration of Visibility Over.In recent years,
During sustained and rapid development of economy, the factors such as sharply increase of industrial activity and motor vehicle constantly aggravate China's atmosphere
A large amount of discharges of the deterioration of quality, high concentration Atmospheric particulates and its precursor, cause serious regional pollution.In recent years,
Intensive North China, city village, gray haze pollution take place frequently, and the hundreds of thousands square kilometre that is usually continuous during gray haze occurs, strong to the public
Health and life bring great threat and inconvenience.
Gray haze is polluted and to be formed and the observation of correlation properties and research are the bases for effectively preventing atmosphere pollution, is also simultaneously
Uncertain very big problem in Atmospheric Chemistry and weather Changeement.Fine particle in troposphere aerosol, the service life is only one
Or so week, time-space resalved technique are very big.Industrial fossil-fueled, motor-vehicle tail-gas, Sand-dust Transmission, the Local anesthesias such as biomass combustion
With regional pollution under DIFFERENT METEOROLOGICAL CONDITIONS driving, air quality and gray haze distribution are affected, therefore, big region gray haze is carried out
Fast slowdown monitoring and quantitative analysis has important practical significance and scientific value.
Due to the aerosol properties and discharge, diffusion characteristic of China's complexity, traditional ground station monitoring exists a lot not true
It is qualitative.Since gray haze pollution layer is vertically and horizontally all being likely to non-uniform Distribution, the ground station of limited quantity
The pollution distribution situation in big region cannot be characterized well, and single-point sampling is likely to over-evaluate the contribution in neighbouring source, monitors generation
Table is poor.Only the method by ground station monitoring is difficult the distribution for the regional pollution that monitoring takes place frequently, change procedure, dirt
Contaminate the large-scale characteristics of total amount etc..
Summary of the invention
In view of this, the embodiment of the present invention be designed to provide in a kind of atmosphere the remote sensing measuring method of gray haze total amount and
Device can identify the gray haze region in gray haze Polluted area in real time and determine the gray haze total amount in the gray haze region, and monitoring is tied
Fruit accuracy is higher and can monitor the regional pollution to take place frequently.
In a first aspect, the embodiment of the invention provides a kind of remote sensing measuring methods of gray haze total amount in atmosphere, comprising:
Obtain the moonscope image of selected gray haze Polluted area；
According to the corresponding measured target identical criterion of the different-waveband of satellite, the gray haze in the moonscope image is identified
Region；The measured target identical criterion includes at least gray haze identical criterion；
According in the gray haze region gray haze intensity, gray haze region area, gray haze thickness and gray haze vertical distribution
Profile determines the gray haze total amount in selected gray haze Polluted area.
With reference to first aspect, the embodiment of the invention provides the first possible embodiments of first aspect, wherein defends
The acquisition methods of the corresponding measured target identical criterion of the different-waveband of star, comprising:
The moonscope image that each wave band of satellite obtains is counted, observed image sample set is set up；
The observed image sample set is classified according to different height layer and ground mulching type, obtains multiple samples
Area；
The selection target sample area from multiple sample areas, and respectively to the measured target in the target sample area into
Row identification and statistics obtain the measured target identical criterion of each wave band corresponding to the satellite.
The possible embodiment of with reference to first aspect the first, the embodiment of the invention provides second of first aspect
Possible embodiment, wherein according to the corresponding measured target identical criterion of the different-waveband of satellite, identify the moonscope
Gray haze region in image, comprising:
According to the corresponding each measured target identical criterion of the different-waveband of satellite, identify in the moonscope image
Each measured target, the measured target include at least: gray haze, cloud, bright earth's surface and snow；
According to the type of each measured target, each measured target obtained respectively to identification is identified；
All gray haze marks are counted, and all gray hazes are identified into corresponding region as gray haze region.
The possible embodiment of second with reference to first aspect, the embodiment of the invention provides the third of first aspect
Possible embodiment, wherein according to the corresponding each measured target identical criterion of the different-waveband of satellite, identify the satellite
Each measured target in observed image, comprising:
Moonscope image is identified according to each wave band of satellite corresponding measured target identical criterion respectively；
When the measured target recognition result obtained according to the measured target identical criterion only includes gray haze result or only
When including non-gray haze result, determine that the measured target is gray haze or non-gray haze；
When the measured target recognition result obtained according to the measured target identical criterion includes gray haze result and at least another
When a kind of non-gray haze result, determine that the measured target is non-gray haze.
The possible embodiment of second with reference to first aspect, the embodiment of the invention provides the 4th kind of first aspect
Possible embodiment, wherein the determination method of gray haze region area, comprising:
Extract the gray haze pixel in the gray haze region；
Projection transform processing is carried out to the gray haze pixel, obtains the gray haze vector data with projection information；
According to the gray haze vector data with projection information, the gray haze region area is calculated.
The possible embodiment of second with reference to first aspect, the embodiment of the invention provides the 5th kind of first aspect
Possible embodiment, wherein the determination method of gray haze thickness, comprising:
Extract the edge gray haze pixel at gray haze edge in the gray haze region；
Whether the height above sea level where judging the edge gray haze pixel is greater than the corresponding height above sea level in region of no relief；
If so, the edge gray haze pixel is matched with digital complex demodulation data, the gray haze edge is obtained
Corresponding longitude and latitude and DEM numerical value；The DEM numerical value, that is, gray haze thickness.
The possible embodiment of second with reference to first aspect, the embodiment of the invention provides the 6th kind of first aspect
Possible embodiment, wherein the determination method of the vertical distribution profiles of gray haze, comprising:
The extinction coefficient profile of the gray haze of acquisition laser radar acquisition and again analytical model meteorological data；It is described to analyze again
Mode meteorological data includes: wind field data and humidity data；
Humidity is carried out to the extinction coefficient profile according to the humidity data in the meteorological data of analytical model again to correct,
Obtain Dry aerosol extinction coefficient profile；
According to the Dry aerosol extinction coefficient profile and the Dry aerosol extinction coefficient profile and ground monitoring ash
The statistical model relationship of haze calculates the mass concentration profile of the ground monitoring gray haze；The mass concentration profile of the gray haze is logical
PM2.5 is crossed to be characterized；
According to the mass concentration profile of the particulate matter feature of different height, meteorologic factor and the ground monitoring gray haze,
Calculate the normalization profile of the gray haze mass concentration；Wherein, the particulate matter feature includes at least: particulate matter component and partial size
Size.
The 6th kind of possible embodiment with reference to first aspect, the embodiment of the invention provides the 7th kind of first aspect
Possible embodiment, wherein according to gray haze intensity, gray haze region area, gray haze thickness and the ash in the gray haze region
The vertical distribution profiles of haze determine the gray haze total amount in selected gray haze Polluted area, comprising:
Interpolation calculation is carried out according to mass concentration of the gray haze region area to the ground monitoring gray haze, is selected
Gray haze Polluted area in gray haze mass concentration；
According to the normalization profile of the gray haze mass concentration, to the quality of gray haze described in selected gray haze Polluted area
Concentration carries out horizontal integral and vertical integral, obtains the gray haze total amount in selected gray haze Polluted area.
With reference to first aspect, the possible embodiment of the first of first aspect to first aspect the 6th kind of possible reality
The possible embodiment of any one in mode is applied, the embodiment of the invention provides the 8th kind of possible embodiment party of first aspect
Formula, wherein according to the corresponding measured target identical criterion of the different-waveband of satellite, identify in the moonscope image
Each measured target, comprising:
When the measured target in the moonscope image in any region meets the gray haze identical criterion, determination is worked as
The measured target of forefoot area is gray haze；The gray haze identical criterion includes at least the following conditions: R0.47Greater than Thr1 and (R0.64-
R2.25)/(R0.64+R2.25) it is greater than Thr2 less than Thr3 and R2.25Greater than Thr4 and BT11Greater than 275 and DEM is greater than 1200 or IGBP
Equal to 16 or 7 or 10；Wherein, R0.47、R0.64、R2.25The moonscope image is respectively indicated in visible light wave range 1, visible light
The reflectivity of wave band 3 and visible light wave range 6；BT11Indicate the moonscope image in the bright temperature of infrared band 14；Thr1,
Thr2, Thr3, Thr4 are the threshold value of corresponding criterion, different according to season different specific threshold values；DEM is altitude data；IGBP is
Ground mulching categorical data.
Second aspect, the embodiment of the invention also provides a kind of remote sensing of gray haze total amount in atmosphere to calculate device, comprising:
Module is obtained, for obtaining the moonscope image of selected gray haze Polluted area；
Identification module identifies that the satellite is seen for the corresponding measured target identical criterion of different-waveband according to satellite
Gray haze region in altimetric image；The measured target identical criterion includes at least gray haze identical criterion；
First determining module, for according in the gray haze region gray haze intensity, gray haze region area, gray haze thickness with
And the vertical distribution profiles of gray haze, determine the gray haze total amount in selected gray haze Polluted area.
The remote sensing measuring method and device of gray haze total amount in a kind of atmosphere provided in an embodiment of the present invention, comprising: obtain choosing
The moonscope image of fixed gray haze Polluted area；According to the corresponding measured target identical criterion of the different-waveband of satellite, identification
Gray haze region in moonscope image；Measured target identical criterion includes at least gray haze identical criterion；According in gray haze region
Gray haze intensity, gray haze region area, gray haze thickness and gray haze vertical distribution profiles, determine selected gray haze Polluted area
In gray haze total amount, it is poor with the method monitoring precision of traditional ground station monitoring gray haze in the prior art and be difficult to monitor
The regional pollution to take place frequently is compared, and can be identified in real time according to the corresponding measured target identical criterion of different-waveband of satellite
Gray haze region in selected gray haze Polluted area simultaneously monitors the gray haze total amount in the gray haze region, realizes gray haze pollution layer
Vertical distribution structure detection, the detection of haze pollution flood thickness and haze polluted space distribution monitoring, can monitor the area to take place frequently
Domain contact scar and monitoring result accuracy is higher.
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, preferred embodiment is cited below particularly, and cooperate
Appended attached drawing, is described in detail below.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached
Figure is briefly described, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not construed as pair
The restriction of range for those of ordinary skill in the art without creative efforts, can also be according to this
A little attached drawings obtain other relevant attached drawings.
Fig. 1 shows the process of the remote sensing measuring method of gray haze total amount in a kind of atmosphere provided by the embodiment of the present invention
Figure；
Fig. 2 shows the processes of the remote sensing measuring method of gray haze total amount in another kind atmosphere provided by the embodiment of the present invention
Figure；
Fig. 3 shows the process of the remote sensing measuring method of gray haze total amount in another kind atmosphere provided by the embodiment of the present invention
Figure；
Fig. 4 shows the process of the remote sensing measuring method of gray haze total amount in another kind atmosphere provided by the embodiment of the present invention
Figure；
Fig. 5 shows the process of the remote sensing measuring method of gray haze total amount in another kind atmosphere provided by the embodiment of the present invention
Figure；
Fig. 6 shows the process of the remote sensing measuring method of gray haze total amount in another kind atmosphere provided by the embodiment of the present invention
Figure；
Fig. 7 shows the process of the remote sensing measuring method of gray haze total amount in another kind atmosphere provided by the embodiment of the present invention
Figure；
The structure that Fig. 8 shows the remote sensing measuring and calculating device of gray haze total amount in a kind of atmosphere provided by the embodiment of the present invention is shown
It is intended to；
Fig. 9 shows the structure of the remote sensing measuring and calculating device of gray haze total amount in another kind atmosphere provided by the embodiment of the present invention
Schematic diagram.
Major Symbol explanation: 11, module is obtained；12, identification module；13, the first computing module；14, the second computing module；
15, statistical module；16, categorization module；17, statistical module is identified.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
Middle attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is only
It is a part of the embodiment of the present invention, instead of all the embodiments.The present invention being usually described and illustrated herein in the accompanying drawings is real
The component for applying example can be arranged and be designed with a variety of different configurations.Therefore, of the invention to what is provided in the accompanying drawings below
The detailed description of embodiment is not intended to limit the range of claimed invention, but is merely representative of selected reality of the invention
Apply example.Based on the embodiment of the present invention, those skilled in the art institute obtained without making creative work
There are other embodiments, shall fall within the protection scope of the present invention.
Since gray haze pollution layer is vertically and horizontally all being likely to non-uniform Distribution, the earth station of limited quantity
Point cannot mark the pollution distribution situation for levying big region well, and single-point sampling is likely to over-evaluate the contribution in neighbouring source, monitor generation
Table is poor.Only the method by ground station monitoring is difficult the distribution for the regional pollution that monitoring takes place frequently, change procedure, dirt
Contaminate the large-scale characteristics of total amount etc..
The fast development of satellite atmosphere remote sensing technology and laser radar atmospheric sounding techniques is INTEGRATED SIGHT and parsing area
The pollution of domain property gray haze provides new opportunity.The detection of gray haze pollution layer vertical space may be implemented in laser radar, and satellite can be real
The detection and spatial distribution monitoring of existing haze pollution flood provide safeguard for realization gray haze pollution real time monitoring.The embodiment of the present invention
Design is that further the information such as the formation of research gray haze and pollution transportation lay base for the cloud, mist, gray haze recognizer in China
Plinth.
Based on this, the embodiment of the invention provides a kind of remote sensing measuring method of gray haze total amount in atmosphere and device, below
It is described by embodiment.
With reference to flow chart shown in FIG. 1, the embodiment of the invention provides a kind of remote sensing measuring and calculating sides of gray haze total amount in atmosphere
Method, which comprises
S101, the moonscope image for obtaining selected gray haze Polluted area.
S102, the corresponding measured target identical criterion of different-waveband according to satellite, identify in the moonscope image
Gray haze region；The measured target identical criterion includes at least gray haze identical criterion.
In the embodiment of the present invention, in advance for the gray haze Polluted area selected in step 101, i.e. monitoring region, statistics is defended
The pixel value for the historical data (the moonscope image in the monitoring region i.e. in preset time period) that star is acquired by each wave band,
These historical datas are set up as observed image sample set, each tested mesh for then including to the observed image sample set of establishment
Mark distinguishes, wherein above-mentioned measured target may include gray haze, cloud, bright earth's surface, clear sky and snow etc., and the present invention is implemented certainly
It is not limited to above-mentioned measured target in example, then respectively to above-mentioned each measured target, such as gray haze, cloud, bright earth's surface, clear sky and snow
Sampling statistics are carried out, determine that above-mentioned each measured target corresponds to the measured target identical criterion of each wave band of satellite, then in advance
First by the storage of above-mentioned measured target identical criterion into database.
In the embodiment of the present invention, satellite mainly passes through the historical data of visible light wave range and infrared band acquisition, wherein
In visible light wave range, the corresponding reflectivity of different measured targets is different, and in infrared band, different measured targets are corresponding
Bright Wen Butong, and be finally calculated according to two parameters of reflectivity and bright temperature and correspond respectively to visible light wave range and infrared waves
The measured target identical criterion of section；Measured target identification data include at least: gray haze identical criterion, cloud identical criterion, clear sky
Identical criterion, bright earth's surface identical criterion and snow identical criterion.
Specifically, satellite is after acquiring moonscope image by visible light wave range and infrared band, likewise, directly
The measured target identical criterion for calling from database and corresponding to visible light wave range and infrared band is connect, above-mentioned quilt is then passed through
Survey target identification criterion respectively in moonscope image reflectivity and bright temperature be compared, distinguish and identify that the satellite is seen
Altimetric image acquisition selection area in measured target, determine each measured targets such as cloud, gray haze, bright earth's surface, clear sky, snow and
Domain of the existence of the above-mentioned each measured target in selection area, and area existing for all gray hazes is counted in above-mentioned recognition result
Domain obtains all gray haze regions.
S103, the gray haze intensity according in the gray haze region, gray haze region area, gray haze thickness and gray haze it is vertical
Distribution profiles determine the gray haze total amount in selected gray haze Polluted area.
In the embodiment of the present invention, the calculating process of each steric parameter is specifically included in above-mentioned gray haze region:
First, the calculation method of gray haze intensity: the intensity of above-mentioned gray haze can be by the PM2.5 data characterization of ground survey station；
Second, the calculation method of gray haze region area: after identification obtains gray haze region, to the gray haze picture in gray haze region
Member carries out projection transform processing, and calculates above-mentioned gray haze region area according to projection transform processing result.
The calculation method of gray haze thickness: third after identification obtains gray haze region, determines the gray haze side in gray haze region
Edge, and the edge gray haze pixel at gray haze edge is obtained, according to DEM (Digital Elevation Model, digital elevation model)
Matching treatment is carried out to edge gray haze pixel, calculates the height of the corresponding longitude and latitude in gray haze edge and gray haze at corresponding longitude and latitude
Degree, the space characteristics of above-mentioned altitude information quantitatively characterizing gray haze obtained by calculation；Above-mentioned altitude information, that is, gray haze thickness.
4th, the calculation method of the vertical distribution profiles of gray haze: after identification obtains gray haze region, surveyed based on laser radar
The extinction coefficient profile (i.e. Aerosol Extinction profile profile) of the gray haze obtained, combined ground and meteorological data (are analyzed again
Mode meteorological data), calculate the mass concentration profile of ground monitoring gray haze, the i.e. vertical distribution of gray haze (PM2.5 mass concentration)
Function f (x), x are highly unit m, and dense divided by the quality of the PM2.5 on minimum one layer (i.e. ground) with the f (x) being calculated
Degree obtains the normalization profile P (x) of gray haze mass concentration；Wherein, P (x) is characteristic, and x is height, unit m.
Then wide according to the vertical distribution of obtained above-mentioned gray haze intensity, gray haze region area, gray haze thickness and gray haze
Line determines the gray haze total amount in selected gray haze Polluted area.It is small that specific determining method, which includes: selected gray haze Polluted area,
Region obtains a laser radar in the selection area and obtains the normalization exterior feature of gray haze for the selection area of above-mentioned zonule
Line P0(x), P is utilized0(x) in the normalization gray haze profile of gray haze height H laser radar below, obtained result is expressed as
Gray haze is distributed in selection area, then, is carried out according to gray haze distribution in selection area to the mass concentration of gray haze in selection area
Horizontal and vertical integral obtains gray haze total amount.
The remote sensing measuring method of gray haze total amount in a kind of atmosphere provided in an embodiment of the present invention is monitored with traditional ground station
The limitation that gray haze technology is difficult to monitor the features such as the range of regional pollution, change procedure, vertical property is compared, being capable of root
According to the corresponding measured target identical criterion of the different-waveband of satellite, the gray haze region in selected gray haze Polluted area is identified in real time
And the gray haze total amount in the gray haze region is monitored, realize the Atmospheric Survey of gray haze pollution layer vertical space, haze pollutes flood
Detection and haze polluted space distribution monitoring, can monitor the regional pollution to take place frequently and monitoring result accuracy is higher.
Further, it with reference to Fig. 2, in the embodiment of the present invention, is identified according to the corresponding measured target of the different-waveband of satellite
Criterion, before identifying the gray haze region in the moonscope image, the different-waveband for obtaining the satellite first is corresponding tested
Target identification criterion, specific acquisition methods, comprising:
The moonscope image that S201, statistics satellite are obtained with each wave band, sets up observed image sample set.
Specifically, being directed to one or more corresponding satellite of survey region, the history of each wave band acquisition of satellite is obtained
Data count these historical data (the moonscope image in the monitoring region i.e. in preset time period) wave band pixel values, set up
Observed image sample set obtains the corresponding quilt of satellite different-waveband to be trained according to obtained observed image sample set
Survey target identification criterion.
In the embodiment of the present invention, satellite mainly passes through the historical data of visible light wave range and infrared band acquisition.
S202, the observed image sample set is classified according to different height layer and ground mulching type, is obtained more
A sample area.
Specifically, different height layer is usually indicated with altitude data, covered in conjunction with the different height layer and earth's surface of altitude data
Lid type classifies the observed image sample set, obtains multiple sample areas, from this multiple sample area select one or
The multiple target sample areas of person obtain satellite not so that the historical data in the target sample area according to selection is trained identification
The corresponding measured target identical criterion with wave band.
S203, the selection target sample area from multiple sample areas, and respectively to tested in the target sample area
Target carries out identification and statistics, obtains the measured target identical criterion of each wave band corresponding to the satellite.
Specifically, the historical data in acquisition target sample area is (i.e. after having selected one or more target sample area
Moonscope image), then to the measured target (such as gray haze, cloud, clear sky, bright earth's surface and snow) in these moonscope images
It is identified and distinguished between, it is each corresponding to satellite to exclude cloud, clear sky, bright earth's surface and avenge the interference to gray haze, and train to obtain
The identical criterion of each measured target of wave band.
In the embodiment of the present invention, satellite mainly passes through the historical data of visible light wave range and infrared band acquisition, wherein
In visible light wave range, the corresponding reflectivity of different measured targets is different, and in infrared band, different measured targets are corresponding
Bright Wen Butong, and be finally calculated according to two parameters of reflectivity and bright temperature and correspond respectively to visible light wave range and infrared waves
The measured target identical criterion of section；Measured target identification data include at least: gray haze identical criterion, cloud identical criterion, clear sky
Identical criterion, bright earth's surface identical criterion and snow identical criterion.
Further, with reference to Fig. 3, in the embodiment of the present invention, in step 102, according to the corresponding quilt of the different-waveband of satellite
Target identification criterion is surveyed, identifies the gray haze region in the moonscope image, comprising:
S1021, the corresponding each measured target identical criterion of different-waveband according to satellite, identify the moonscope figure
Each measured target as in.
Specifically, when same atural object (i.e. same measured target) had not only been identified as gray haze but also had been identified as non-gray haze, then
The atural object is ultimately determined to non-gray haze.
S1022, according to the type of each measured target, the obtained each measured target of identification is identified respectively.
All gray haze marks of S1023, statistics, and all gray hazes are identified into corresponding region as gray haze region.
In conjunction with step 1021- step 1023, above-mentioned identical criterion is included at least: gray haze identical criterion, cloud identical criterion, fine
Empty identical criterion, bright earth's surface identical criterion and snow identical criterion；
Specifically, being based on above-mentioned each measured target identical criterion, gray haze is tentatively identified with gray haze identical criterion first
Range, and be 1 by the gray haze area identification of identification；Cloud is identified with cloud identical criterion, and is 2 by the cloud sector domain identifier of identification；With
Identical criterion identification snow is avenged, and is 3 by the snow area identification of identification；Then, bright picture dot earth's surface is identified with bright earth's surface identical criterion,
And the bright picture dot earth surface area of identification is set as 0, the bright earth's surface area of removal gray haze misrecognition.It is carried out to moonscope image
After identification, statistics is identified as 1 all areas, as gray haze region.
When the measured target in the moonscope image in any region meets the gray haze identical criterion, determination is worked as
The measured target of forefoot area is gray haze；The gray haze identical criterion includes at least the following conditions: R0.47Greater than Thr1 and (R0.64-
R2.25)/(R0.64+R2.25) it is greater than Thr2 less than Thr3 and R2.25Greater than Thr4 and BT11Greater than 275 and DEM is greater than 1200 or IGBP
Equal to 16 or 7 or 10；Wherein, R0.47、R0.64、R2.25The moonscope image is respectively indicated in visible light wave range 1, visible light
The reflectivity of wave band 3 and visible light wave range 6；BT11Indicate the moonscope image in the bright temperature of infrared band 14；Thr1,
Thr2, Thr3, Thr4 are the threshold value of corresponding criterion, different according to season different specific threshold values；DEM is altitude data；IGBP is
Ground mulching categorical data.
When the moonscope image meets cloud identical criterion, the quilt of matching area in the moonscope image is determined
Survey target is cloud；The identification of cloud includes following criterion:
1)R0.47The absolute standard of 3*3 window it is poor (standard deviation * average reflectance/) it is greater than Thrc1, standard deviation is greater than
Thrc2, then all wave band 3*3 windows are all judged as cloud.
2)R0.47Greater than Thrc3, criterion identification is spissatus.
3)R0.87-R0.64Greater than Thrc5 and it is less than Thrc6 or is less than Thrc7, which identifies spissatus and water body.
4)R0.87/R0.64Greater than Thrc4, cloud is further identified.
5)BT3.9-BT10.4Greater than Thrc9, which identifies the thin cloud of low clouds and fritter cumulus.
6)BT11-BT13.3Less than Thrc11, which identifies cloud edge.
7)BT13.3Less than Thrc9, which identifies high cloud and Yun Bianyuan.
Wherein, R0.47、R0.64、R0.87Moonscope image is respectively indicated in visible light wave range 1, visible light wave range 3 and visible
The reflectivity of optical band 4；BT3.9、BT10.4、BT11、BT13.3, respectively indicate satellite infrared band 7, infrared band 13,
The bright temperature of infrared band 14 and infrared band 16；Thrc is the threshold value that cloud identifies corresponding criterion.When the moonscope figure
As when meeting bright earth's surface identical criterion, determining that the measured target of matching area in the moonscope image is bright earth's surface；Brightly
The identification of table includes following criterion:
1) (R0.64-R2.25)/(R0.64+R2.25) is less than Thrb1, while R2.25 is greater than Thrb2.
2) when pixel vegetation account for it is leading, then (R0.64-R2.25)/(R0.64+R2.25) be greater than Thrb3.
Wherein, R1.6、R2.25 respectively indicates moonscope image in visible light wave range 5, the reflectivity of visible light wave range 6；
Thrb is the threshold value that cloud identifies corresponding criterion, different and different according to season, as Various Seasonal be corresponding with the Thrb1 of different value with
Thrb2。
When the moonscope image meets snow identical criterion, the quilt of matching area in the moonscope image is determined
Surveying target is snow；In the embodiment of the present invention, snow identical criterion is R0.51-R0.64Greater than 0.4；Wherein, R0.51、R0.64Indicate that satellite is seen
Altimetric image is in visible light wave range 2, the reflectivity of visible light wave range 3.
Specifically, by the reflectivity and bright temperature of visible light wave range each in moonscope image and each infrared band point
(specifically include: gray haze identical criterion, cloud identical criterion, clear sky identical criterion, bright earth's surface do not identify with measured target identical criterion
Criterion and snow identical criterion) in reflectivity and bright temperature condition be compared, each region that identification satellite observed image includes
Measured target information.
Further, in atmosphere provided in an embodiment of the present invention gray haze total amount remote sensing measuring method, step 1021, according to
The corresponding each measured target identical criterion of the different-waveband of satellite, identifies each tested mesh in the moonscope image
Mark, comprising:
1021a, respectively according to the corresponding measured target identical criterion of each wave band of satellite to moonscope image into
Row identification.
1021b, when the measured target recognition result obtained according to the measured target identical criterion only includes gray haze result
Determine that the measured target is gray haze or non-gray haze when perhaps only including non-gray haze result.
1021c, when the measured target recognition result obtained according to the measured target identical criterion include gray haze result and
At least another non-gray haze result when, determine the measured target be non-gray haze.
Specifically, the priority of cloud identical criterion is higher than gray haze identical criterion in the embodiment of the present invention, identical criterion is avenged
Priority is higher than gray haze identical criterion, and the priority of bright earth's surface identical criterion is higher than gray haze identical criterion, i.e., when moonscope figure
Some region part as in is not only identified as gray haze, but also is identified as cloud or snow or bright earth's surface, then finally determines the area
Domain part is cloud or snow or bright earth's surface, rather than gray haze.
Further, with reference to Fig. 4, in the embodiment of the present invention, the remote sensing measuring method of gray haze total amount, gray haze region in atmosphere
The determination method of area, comprising:
Gray haze pixel in S301, the extraction gray haze region.
S302, projection transform processing is carried out to the gray haze pixel, obtains the gray haze vector data with projection information.
S303, the gray haze vector data that projection information is had according to, calculate the gray haze region area.
In conjunction with step 301- step 303, specifically, in above-mentioned steps, the embodiment of the present invention is by the region of all marks 1
As gray haze region, then take out be identified as 1 gray haze region, all pixels in gray haze region are obtained in this step, will be taken
All pixels out are converted to shapefile file, then set the projection information of shapefile file to be suitble to China
Albert authalic projection is finally calculated the area of shapefile file using arcgis software, obtains gray haze region area S.
Further, with reference to Fig. 5, in the embodiment of the present invention, in atmosphere in the remote sensing measuring method of gray haze total amount, gray haze is thick
The determination method of degree, comprising:
S401, the edge gray haze pixel for extracting gray haze edge in the gray haze region.
S402, judge the edge gray haze pixel where height above sea level whether be greater than the corresponding height above sea level in region of no relief.
S403, if so, the edge gray haze pixel is matched with digital complex demodulation data, obtain the ash
The corresponding longitude and latitude in haze edge and DEM numerical value；The DEM numerical value, that is, gray haze thickness.
In conjunction with step 401- step 403, in the embodiment of the present invention, based on landform (such as Yanshan Mountain in North China region, Taihang Mountain with
And the middle areas in Shandong Province mountainous region in Shandong) method of measuring and calculating gray haze thickness is illustrated: it is homogeneous within a certain range that atmosphere, which is assumed to be,
, based on the gray haze region identified in above-mentioned steps, gray haze edge is determined from the gray haze region, and extract the picture at gray haze edge
Then member then ignore by pixel progress height differentiation if it is determined that pixel is in region of no relief to the pixel at gray haze edge；
It then proceeds as follows, specifically includes if it is determined that pixel is in mountain area: by the pixel at gray haze edge and dem data progress
Match, obtains the corresponding longitude and latitude in gray haze edge and DEM numerical value；Wherein, above-mentioned DEM numerical value is gray haze at corresponding longitude and latitude
Highly (i.e. gray haze thickness)；It is assumed that it can be by height approximation characterization one of the gray haze at corresponding longitude and latitude based on above-mentioned homogenieity
Determine the height H of gray haze in geographic coverage, the space characteristics of quantitatively characterizing gray haze.
Further, with reference to Fig. 6, in the embodiment of the present invention, the calculation method of the vertical distribution profiles of the gray haze, packet
It includes:
The extinction coefficient profile for the gray haze that S501, acquisition laser radar acquire and again analytical model meteorological data；It is described
Analytical model meteorological data includes: wind field data and humidity data again.
In the embodiment of the present invention, the particle contamination in the preset height layer of laser radar website longitude and latitude corresponding position is obtained
The extinction coefficient of object；Wherein, above-mentioned preset height layer such as 1000hpa-700hpa corresponds to sea level and plays 0m-3000m's or so
Highly, which may be considered Boundary Layer Height, and general pollutant all concentrates in this layer.
Above-mentioned pattern simulation analyzes meteorological field data (i.e. analytical model meteorological data again) again, is by actual observation gas
As the revised meteorological field data of data, there is higher reliability than forecast fields data, having become in the world at present can quilt
Direct applied meteorological data resource.Analyzing a data (i.e. analytical model meteorological data again) again can not only provide by height layer
Wind field data, and can equally obtain the relative humidity data by height layer, wind field data and relative humidity data all have
There is higher precision, the spatial resolution of mesoscale model output data (such as WRF, MM5 mode) is up to 5km ﹡ 5km, completely may be used
To meet accuracy requirement and space requirement to wind field and humidity field data.And analyzing field data again is that pattern simulation assimilation obtains
, acquisition is relatively easy, and cost is very cheap.
S502, humidity is carried out to the extinction coefficient profile according to the humidity data in the analytical model meteorological data again
It corrects, obtains Dry aerosol extinction coefficient profile.
Since in natural situation, particulate absorbs the moisture in surrounding air, moisture absorption occurs and increases.Laser radar
The extinction coefficient for measuring particulate in surrounding air is the extinction coefficient of wet particle, and the particulate pollutant of ground observation
Mass concentration be then it is dry after dried particle (i.e. Dry aerosol particle) mass concentration, so needing according to analytical model again
Humidity data in meteorological data carries out humidity to the extinction coefficient that laser radar measures and corrects, to obtain Dry aerosol delustring system
Number.
S503, it is supervised according to the Dry aerosol extinction coefficient profile and the Dry aerosol extinction coefficient profile and ground
The statistical model relationship for surveying gray haze, calculates the mass concentration profile of the ground monitoring gray haze；The mass concentration of the gray haze is wide
Line is characterized by PM2.5.
Specifically, to calculate laser radar transport fluxes, it is necessary first to establish Dry aerosol extinction coefficient KEX_d with it is right
The conversion relationship of the PM of the observation at moment is answered, particulate matter is then calculated according to above-mentioned conversion relationship and Dry aerosol extinction coefficient and is hung down
Straight distribution profiles, wherein the mass concentration PM of the particulate matter in the particulate matter vertical distribution profiles including flood.
S504, according to the mass concentration of the particulate matter feature of different height, meteorologic factor and the ground monitoring gray haze
Profile calculates the normalization profile of the gray haze mass concentration；Wherein, the particulate matter feature includes at least: particulate matter component
And particle size.
Specifically, the mass concentration of ground monitoring gray haze is obtained divided by the mass concentration of the gray haze on minimum one layer (i.e. ground)
Obtain the normalization profile of gray haze mass concentration.
In conjunction with step 501- step 504, based on the Aerosol Extinction profile that laser radar measures, combined ground is gentle
As data, the vertical distribution function f (x) for obtaining gray haze (PM2.5 mass concentration) is calculated, using obtained f (x) divided by minimum one
The mass concentration of the PM2.5 of layer (ground), obtains the normalization profile P (x) of gray haze；Above-mentioned P (x) is characteristic, wherein x
For height, unit m.
Specifically, above-mentioned P (x) is the PM2.5 mass concentration normalization coefficient array on one group of different height, difference is represented
Multiple of the gray haze concentration of height layer relative to ground gray haze concentration.Assumed based on atmosphere homogenieity, it is believed that in this laser radar
The gray haze profile P of detection0, represent the vertical distribution of gray haze in a certain range.For zonule, a laser radar number can be used
According to representing gray haze profile；For extensive area, laser radar quantity can be increased, to obtain bigger representativeness, Mei Gelei
The gray haze profile reached is Pi, i is any one radar.Layouting in combination with meteorological field and contamination characteristics for big band radar, carries out
There is selection to layout.
Further, with reference to Fig. 7, in atmosphere provided in an embodiment of the present invention in the remote sensing measuring method of gray haze total amount, step
In rapid 103, according in the gray haze region gray haze intensity, gray haze region area, gray haze thickness and gray haze vertical distribution
Profile determines the gray haze total amount in selected gray haze Polluted area, comprising:
S1041, interpolation calculation is carried out according to mass concentration of the gray haze region area to the ground monitoring gray haze, obtained
The mass concentration of gray haze into selected gray haze Polluted area.
S1042, according to the normalization profile of the gray haze mass concentration, to gray haze described in selected gray haze Polluted area
Mass concentration carry out horizontal integral and vertical integral, obtain the gray haze total amount in selected gray haze Polluted area.
In conjunction with step 1041- step 1042, the intensity of gray haze by ground survey station PM2.5 data characterization, gray haze range and
Area is provided by satellite, the PM2.5 of survey station within the scope of gray haze can be carried out space interpolation, be obtained the mass concentration of region gray haze
PMi, unit μ g/m3.For zonule, the gray haze profile P of a laser radar is utilized0(x), interception gray haze height H is below swashs
The gray haze profile of optical radar, it is believed that gray haze is distributed all as shown in profile, by following formula to gray haze mass concentration in zonule
Horizontal integral and vertical integral are carried out, gray haze total amount is obtained:Wherein, formula (1)
Middle S is gray haze area, and H is gray haze thickness；For big region, respectively a certain range of gray haze profile is represented with several radars
Pi(x), it is assumed that used N platform radar altogether, then gray haze total amount isI.e. respectively in each thunder
Horizontal and vertical integral is carried out in the representative region reached.
The embodiment of the present invention provides a kind of remote sensing measuring method of gray haze total amount in atmosphere, monitors ash with traditional ground station
The limitation that haze technology is difficult to monitor the features such as the range of regional pollution, change procedure, vertical property is compared, being capable of basis
The corresponding measured target identical criterion of the different-waveband of satellite, the gray haze region in gray haze Polluted area that identification is selected in real time is simultaneously
The gray haze total amount in the gray haze region is monitored, the spy of the Atmospheric Survey, haze pollution flood of gray haze pollution layer vertical space is realized
Survey and haze polluted space distribution monitoring, can monitor the regional pollution to take place frequently and monitoring result accuracy is higher.
The embodiment of the invention also provides a kind of remote sensing of gray haze total amount in atmosphere to calculate device, and described device is for executing
The remote sensing measuring method of gray haze total amount in above-mentioned atmosphere, with reference to Fig. 8, described device is specifically included:
Module 11 is obtained, for obtaining the moonscope image of selected gray haze Polluted area；
Identification module 12 identifies the satellite for the corresponding measured target identical criterion of different-waveband according to satellite
Gray haze region in observed image；The measured target identical criterion includes at least gray haze identical criterion；
First determining module 13, for according to the gray haze intensity in the gray haze region, gray haze region area, gray haze thickness
And the vertical distribution profiles of gray haze, determine the gray haze total amount in selected gray haze Polluted area.
Further, the remote sensing measuring and calculating device of gray haze total amount is also wrapped with reference to Fig. 9 in atmosphere provided in an embodiment of the present invention
It includes:
Statistical module 15, the moonscope image obtained for counting each wave band of satellite, sets up observed image sample set；
Categorization module 16, for being divided the observed image sample set according to different height layer and ground mulching type
Class obtains multiple sample areas；
Statistical module 17 is identified, for the selection target sample area from multiple sample areas, and respectively to the target
Measured target in sample area carries out identification and statistics, and the measured target identification for obtaining each wave band corresponding to the satellite is sentenced
According to.
Further, device, identification module 12, packet are calculated in the remote sensing of gray haze total amount in atmosphere provided in an embodiment of the present invention
It includes:
Recognition unit is defended described in identification for the corresponding each measured target identical criterion of the different-waveband according to satellite
Each measured target in star observed image；
Unit is identified, for the type according to each measured target, each measured target that identification obtains is carried out respectively
Mark；
Statistic unit identifies corresponding region as gray haze area for counting all gray haze marks, and using all gray hazes
Domain.
Further, device is calculated in the remote sensing of gray haze total amount in atmosphere provided in an embodiment of the present invention, and recognition unit is specific
Moonscope image is identified for each wave band corresponding measured target identical criterion according to satellite；When according to institute
State the measured target recognition result that measured target identical criterion obtains only include gray haze result or only include non-gray haze result when,
Determine that the measured target is gray haze or non-gray haze；When the measured target identification obtained according to the measured target identical criterion
As a result include gray haze result and at least another non-gray haze result when, determine the measured target for non-gray haze.
Further, device is calculated in the remote sensing of gray haze total amount in atmosphere provided in an embodiment of the present invention, further includes:
First extraction unit, for extracting the gray haze pixel in the gray haze region；
Projection transform processing unit is obtained for carrying out projection transform processing to the gray haze pixel with projection information
Gray haze vector data；
First computing unit, for calculating the gray haze region according to the gray haze vector data for having projection information
Area.
Further, device is calculated in the remote sensing of gray haze total amount in atmosphere provided in an embodiment of the present invention, further includes:
Second extraction unit, for extracting the edge gray haze pixel at gray haze edge in the gray haze region；
Whether judging unit is greater than the corresponding sea in region of no relief for the height above sea level where judging the edge gray haze pixel
Degree of lifting, and under conditions of being, the edge gray haze pixel is matched with digital complex demodulation data, is obtained described
The corresponding longitude and latitude in gray haze edge and DEM numerical value；The DEM numerical value, that is, gray haze thickness.
Further, device is calculated in the remote sensing of gray haze total amount in atmosphere provided in an embodiment of the present invention, further includes:
First acquisition unit, the extinction coefficient profile and analytical model gas again of the gray haze for obtaining laser radar acquisition
As data；The meteorological data of analytical model again includes: wind field data and humidity data；
Humidity corrects unit, for the humidity data in the basis analytical model meteorological data again to the extinction coefficient
Profile carries out humidity and corrects, and obtains Dry aerosol extinction coefficient profile；
Second computing unit, for according to the Dry aerosol extinction coefficient profile and the Dry aerosol extinction coefficient
The statistical model relationship of profile and ground monitoring gray haze calculates the mass concentration profile of the ground monitoring gray haze；The gray haze
Mass concentration profile characterized by PM2.5；
Third computing unit, for particulate matter feature, meteorologic factor and the ground monitoring ash according to different height
The mass concentration profile of haze calculates the normalization profile of the gray haze mass concentration；Wherein, the particulate matter feature is at least wrapped
It includes: particulate matter component and particle size.
Further, in atmosphere provided in an embodiment of the present invention in the remote sensing measuring and calculating device of gray haze total amount, first determines mould
Block 13, comprising:
Interpolation calculation unit, for being carried out according to mass concentration of the gray haze region area to the ground monitoring gray haze
Interpolation calculation obtains the mass concentration of gray haze in selected gray haze Polluted area；
Integral calculation unit, for the normalization profile according to the gray haze mass concentration, to selected gray haze contaminated area
The mass concentration of gray haze described in domain carries out horizontal integral and vertical integral, and the gray haze obtained in selected gray haze Polluted area is total
Amount.
Further, in atmosphere provided in an embodiment of the present invention in the remote sensing measuring and calculating device of gray haze total amount, recognition unit tool
Body is used for, and when the measured target in the moonscope image in any region meets the gray haze identical criterion, determination is worked as
The measured target of forefoot area is gray haze；The gray haze identical criterion includes at least the following conditions: R0.47Greater than Thr1 and (R0.64-
R2.25)/(R0.64+R2.25) it is greater than Thr2 less than Thr3 and R2.25Greater than Thr4 and BT11Greater than 275 and DEM is greater than 1200 or IGBP
Equal to 16 or 7 or 10；Wherein, R0.47、R0.64、R2.25The moonscope image is respectively indicated in visible light wave range 1, visible light
The reflectivity of wave band 3 and visible light wave range 6；BT11Indicate the moonscope image in the bright temperature of infrared band 14；Thr1,
Thr2, Thr3, Thr4 are the threshold value of corresponding criterion, different according to season different specific threshold values；DEM is altitude data；IGBP is
Ground mulching categorical data.
The remote sensing measuring method of gray haze total amount in a kind of atmosphere provided in an embodiment of the present invention, with tradition in the prior art
The method monitoring precision of ground station monitoring gray haze is poor and is difficult the regional pollution that takes place frequently of monitoring and compares, being capable of basis
The corresponding measured target identical criterion of the different-waveband of satellite, the gray haze region in gray haze Polluted area that identification is selected in real time is simultaneously
The gray haze total amount in the gray haze region is monitored, the spy of the Atmospheric Survey, haze pollution flood of gray haze pollution layer vertical space is realized
Survey and haze polluted space distribution monitoring, can monitor the regional pollution to take place frequently and monitoring result accuracy is higher.
The device of the remote sensing measuring and calculating of gray haze total amount can be in equipment in a kind of atmosphere provided by the embodiment of the present invention
Specific hardware or the software being installed in equipment or firmware etc..Device provided by the embodiment of the present invention, realization principle and
The technical effect of generation is identical with preceding method embodiment, and to briefly describe, Installation practice part does not refer to place, can refer to
Corresponding contents in preceding method embodiment.It is apparent to those skilled in the art that for the convenience and letter of description
It is clean, system, the specific work process of device and unit of foregoing description, corresponding to during reference can be made to the above method embodiment
Journey, details are not described herein.
In embodiment provided by the present invention, it should be understood that disclosed device and method, it can be by others side
Formula is realized.The apparatus embodiments described above are merely exemplary, for example, the division of the unit, only one kind are patrolled
Function division is collected, there may be another division manner in actual implementation, in another example, multiple units or components can combine or can
To be integrated into another system, or some features can be ignored or not executed.Another point, shown or discussed is mutual
Coupling, direct-coupling or communication connection can be INDIRECT COUPLING or communication link by some communication interfaces, device or unit
It connects, can be electrical property, mechanical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit
The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple
In network unit.It can select some or all of unit therein according to the actual needs to realize the mesh of this embodiment scheme
's.
In addition, each functional unit in embodiment provided by the invention can integrate in one processing unit, it can also
To be that each unit physically exists alone, can also be integrated in one unit with two or more units.
It, can be with if the function is realized in the form of SFU software functional unit and when sold or used as an independent product
It is stored in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be embodied in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be a
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes: that USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic or disk.
It should also be noted that similar label and letter indicate similar terms in following attached drawing, therefore, once a certain Xiang Yi
It is defined in a attached drawing, does not then need that it is further defined and explained in subsequent attached drawing, in addition, term " the
One ", " second ", " third " etc. are only used for distinguishing description, are not understood to indicate or imply relative importance.
Finally, it should be noted that embodiment described above, only a specific embodiment of the invention, to illustrate the present invention
Technical solution, rather than its limitations, scope of protection of the present invention is not limited thereto, although with reference to the foregoing embodiments to this hair
It is bright to be described in detail, those skilled in the art should understand that: anyone skilled in the art
In the technical scope disclosed by the present invention, it can still modify to technical solution documented by previous embodiment or can be light
It is readily conceivable that variation or equivalent replacement of some of the technical features；And these modifications, variation or replacement, do not make
The essence of corresponding technical solution is detached from the spirit and scope of technical solution of the embodiment of the present invention.Should all it cover in protection of the invention
Within the scope of.Therefore, protection scope of the present invention should be based on the protection scope of the described claims.
Claims (9)
1. the remote sensing measuring method of gray haze total amount in a kind of atmosphere characterized by comprising
Obtain the moonscope image of selected gray haze Polluted area；
According to the corresponding measured target identical criterion of the different-waveband of satellite, the gray haze area in the moonscope image is identified
Domain；The measured target identical criterion includes at least gray haze identical criterion；
According in the gray haze region gray haze intensity, gray haze region area, gray haze thickness and gray haze vertical distribution profiles,
Determine the gray haze total amount in selected gray haze Polluted area；
According to the corresponding measured target identical criterion of the different-waveband of satellite, identify each in the moonscope image
Measured target, comprising:
When the measured target in the moonscope image in any region meets the gray haze identical criterion, then primarily determine
The measured target of current region is gray haze；
The gray haze identical criterion includes at least the following conditions:
R0.47Greater than Thr1；
(R0.64-R2.25)/(R0.64+R2.25) it is greater than Thr2 less than Thr3 and R2.25Greater than Thr4；
BT11Greater than 275 and DEM is greater than 1200；Alternatively, BT11Greater than 275 and IGBP is equal to 16；Alternatively, BT11Greater than 275 and
IGBP is equal to 10；
Wherein, R0.47、R0.64、R2.25The moonscope image is respectively indicated in visible light wave range 1, visible light wave range 3 and visible
The reflectivity of optical band 6；BT11Indicate the moonscope image in the bright temperature of infrared band 14；Thr1,Thr2,Thr3,
Thr4 is the threshold value of corresponding criterion, different according to season different specific threshold values；DEM is altitude data；IGBP is ground mulching class
Type data；
Identify cloud, snow, bright earth's surface simultaneously, when the gray haze primarily determined is identified as cloud or snow or bright earth's surface, then the region
It is ultimately determined to non-gray haze, when the gray haze primarily determined and is not recognized as cloud or snow or bright earth's surface, then the region is final
It is determined as gray haze.
2. the remote sensing measuring method of gray haze total amount in atmosphere according to claim 1, which is characterized in that the different waves of satellite
The acquisition methods of the corresponding measured target identical criterion of section, comprising:
The moonscope image that each wave band of satellite obtains is counted, observed image sample set is set up；
The observed image sample set is classified according to different height layer and ground mulching type, obtains multiple sample areas；
The selection target sample area from multiple sample areas, and the measured target in the target sample area is known respectively
It not and counts, obtains the measured target identical criterion of each wave band corresponding to the satellite.
3. the remote sensing measuring method of gray haze total amount in atmosphere according to claim 2, which is characterized in that not according to satellite
The corresponding measured target identical criterion with wave band, identifies the gray haze region in the moonscope image, comprising:
According to the corresponding each measured target identical criterion of the different-waveband of satellite, identify each in the moonscope image
Measured target；
According to the type of each measured target, each measured target obtained respectively to identification is identified；
All gray haze marks are counted, and all gray hazes are identified into corresponding region as gray haze region.
4. the remote sensing measuring method of gray haze total amount in atmosphere according to claim 3, which is characterized in that not according to satellite
With wave band, corresponding each measured target identical criterion, identifies each measured target in the moonscope image, comprising:
Moonscope image is identified according to each wave band of satellite corresponding measured target identical criterion respectively；
When the measured target recognition result obtained according to the measured target identical criterion only includes gray haze result or only includes
When non-gray haze result, determine that the measured target is gray haze or non-gray haze；
When the measured target recognition result obtained according to the measured target identical criterion includes gray haze result and at least another
When non-gray haze result, determine that the measured target is non-gray haze.
5. the remote sensing measuring method of gray haze total amount in atmosphere according to claim 3, which is characterized in that gray haze region area
Determination method, comprising:
Extract the gray haze pixel in the gray haze region；
Projection transform processing is carried out to the gray haze pixel, obtains the gray haze vector data with projection information；
According to the gray haze vector data with projection information, the gray haze region area is calculated.
6. the remote sensing measuring method of gray haze total amount in atmosphere according to claim 3, which is characterized in that gray haze thickness is really
Determine method, comprising:
Extract the edge gray haze pixel at gray haze edge in the gray haze region；
Whether the height above sea level where judging the edge gray haze pixel is greater than the corresponding height above sea level in region of no relief；
If so, the edge gray haze pixel is matched with digital complex demodulation data, it is corresponding to obtain the gray haze edge
Longitude and latitude and DEM numerical value；The DEM numerical value, that is, gray haze thickness.
7. the remote sensing measuring method of gray haze total amount in atmosphere according to claim 3, which is characterized in that vertical point of gray haze
The determination method of cloth profile, comprising:
The extinction coefficient profile of the gray haze of acquisition laser radar acquisition and again analytical model meteorological data；The analytical model again
Meteorological data includes: wind field data and humidity data；
Humidity is carried out to the extinction coefficient profile according to the humidity data in the meteorological data of analytical model again to correct, and is obtained
Dry aerosol extinction coefficient profile；
According to the Dry aerosol extinction coefficient profile and the Dry aerosol extinction coefficient profile and ground monitoring gray haze
Statistical model relationship calculates the mass concentration profile of the ground monitoring gray haze；The mass concentration profile of the gray haze passes through
PM2.5 is characterized；
According to the mass concentration profile of the particulate matter feature of different height, meteorologic factor and the ground monitoring gray haze, calculate
The normalization profile of the gray haze mass concentration；Wherein, the particulate matter feature includes at least: particulate matter component and partial size are big
It is small.
8. the remote sensing measuring method of gray haze total amount in atmosphere according to claim 7, which is characterized in that according to the gray haze
The vertical distribution profiles of gray haze intensity, gray haze region area, gray haze thickness and gray haze in region determine that selected gray haze is dirty
Contaminate the gray haze total amount in region, comprising:
Interpolation calculation is carried out according to mass concentration of the gray haze region area to the ground monitoring gray haze, obtains selected ash
The mass concentration of gray haze in haze Polluted area；
According to the normalization profile of the gray haze mass concentration, to the mass concentration of gray haze described in selected gray haze Polluted area
Horizontal integral and vertical integral are carried out, the gray haze total amount in selected gray haze Polluted area is obtained.
9. device is calculated in the remote sensing of gray haze total amount in a kind of atmosphere characterized by comprising
Module is obtained, for obtaining the moonscope image of selected gray haze Polluted area；
Identification module identifies the moonscope figure for the corresponding measured target identical criterion of different-waveband according to satellite
Gray haze region as in；The measured target identical criterion includes at least gray haze identical criterion；
First determining module, for according to gray haze intensity, gray haze region area, gray haze thickness and the ash in the gray haze region
The vertical distribution profiles of haze determine the gray haze total amount in selected gray haze Polluted area；
The recognition unit is specifically used for, when the measured target in the moonscope image in any region meets the gray haze
When identical criterion, then primarily determine that the measured target of current region is gray haze；
The gray haze identical criterion includes at least the following conditions:
R0.47Greater than Thr1；
(R0.64-R2.25)/(R0.64+R2.25) it is greater than Thr2 less than Thr3 and R2.25Greater than Thr4；
BT11Greater than 275 and DEM is greater than 1200；Alternatively, BT11Greater than 275 and IGBP is equal to 16；Alternatively, BT11Greater than 275 and
IGBP is equal to 10；
Wherein, R0.47、R0.64、R2.25The moonscope image is respectively indicated in visible light wave range 1, visible light wave range 3 and visible
The reflectivity of optical band 6；BT11Indicate the moonscope image in the bright temperature of infrared band 14；Thr1,Thr2,Thr3,
Thr4 is the threshold value of corresponding criterion, different according to season different specific threshold values；DEM is altitude data；IGBP is ground mulching class
Type data；
Identify cloud, snow, bright earth's surface simultaneously, when the gray haze primarily determined is identified as cloud or snow or bright earth's surface, then the region
It is ultimately determined to non-gray haze, when the gray haze primarily determined and is not recognized as cloud or snow or bright earth's surface, then the region is final
It is determined as gray haze.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710244393.4A CN106950574B (en) | 2017-04-14 | 2017-04-14 | The remote sensing measuring method and device of gray haze total amount in a kind of atmosphere |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710244393.4A CN106950574B (en) | 2017-04-14 | 2017-04-14 | The remote sensing measuring method and device of gray haze total amount in a kind of atmosphere |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106950574A CN106950574A (en) | 2017-07-14 |
| CN106950574B true CN106950574B (en) | 2019-11-08 |
Family
ID=59476171
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710244393.4A Active CN106950574B (en) | 2017-04-14 | 2017-04-14 | The remote sensing measuring method and device of gray haze total amount in a kind of atmosphere |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106950574B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107972887B (en) * | 2017-11-22 | 2019-11-08 | 北京航空航天大学 | The monitoring method and device of PM2.5 based on Cubesat satellite |
| CN108427905A (en) * | 2017-11-27 | 2018-08-21 | 易瓦特科技股份公司 | The method, apparatus and system warned to gray haze source based on remote server |
| CN108427906A (en) * | 2017-11-27 | 2018-08-21 | 易瓦特科技股份公司 | The method, apparatus and system warned to gray haze source based on earth station |
| CN112102433B (en) * | 2020-09-17 | 2021-06-01 | 中科三清科技有限公司 | Method and device for drawing vertical distribution map of air pollutants and storage medium |
| CN114005049B (en) * | 2022-01-05 | 2022-04-12 | 国家卫星气象中心(国家空间天气监测预警中心) | Haze area detection method, device, equipment and medium based on multiband remote sensing information |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102176073A (en) * | 2011-01-24 | 2011-09-07 | 环境保护部卫星环境应用中心 | Ambient air quality comprehensive monitoring method based on first ambient satellite |
| CN102539336A (en) * | 2011-02-01 | 2012-07-04 | 环境保护部卫星环境应用中心 | Method and system for estimating inhalable particles based on HJ-1 satellite |
| CN106124374A (en) * | 2016-07-22 | 2016-11-16 | 中科宇图科技股份有限公司 | Atmospheric particulates remote-sensing monitoring method based on data fusion |
| CN106126903A (en) * | 2016-06-21 | 2016-11-16 | 北京市环境保护监测中心 | The computational methods of PM2.5 three-dimensional transport fluxes and device |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8340461B2 (en) * | 2010-02-01 | 2012-12-25 | Microsoft Corporation | Single image haze removal using dark channel priors |
-
2017
- 2017-04-14 CN CN201710244393.4A patent/CN106950574B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102176073A (en) * | 2011-01-24 | 2011-09-07 | 环境保护部卫星环境应用中心 | Ambient air quality comprehensive monitoring method based on first ambient satellite |
| CN102539336A (en) * | 2011-02-01 | 2012-07-04 | 环境保护部卫星环境应用中心 | Method and system for estimating inhalable particles based on HJ-1 satellite |
| CN106126903A (en) * | 2016-06-21 | 2016-11-16 | 北京市环境保护监测中心 | The computational methods of PM2.5 three-dimensional transport fluxes and device |
| CN106124374A (en) * | 2016-07-22 | 2016-11-16 | 中科宇图科技股份有限公司 | Atmospheric particulates remote-sensing monitoring method based on data fusion |
Non-Patent Citations (1)
| Title |
|---|
| 霾光谱特性分析与卫星遥感识别算法;葛巍 等;《光谱学与光谱分析》;20161231;第36卷(第12期);第3818-3823页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106950574A (en) | 2017-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106950574B (en) | The remote sensing measuring method and device of gray haze total amount in a kind of atmosphere | |
| Banks et al. | Performance evaluation of the boundary-layer height from lidar and the Weather Research and Forecasting model at an urban coastal site in the north-east Iberian Peninsula | |
| Moseley et al. | Probing the precipitation life cycle by iterative rain cell tracking | |
| CA2792185C (en) | System and method for identifying individual trees in lidar data using local view | |
| Qi et al. | A real‐time automated convective and stratiform precipitation segregation algorithm in native radar coordinates | |
| Gossart et al. | Blowing snow detection from ground-based ceilometers: application to East Antarctica | |
| Su et al. | Improving MODIS sea ice detectability using gray level co-occurrence matrix texture analysis method: A case study in the Bohai Sea | |
| Anderson et al. | Long-term statistical assessment of Aqua-MODIS aerosol optical depth over coastal regions: bias characteristics and uncertainty sources | |
| CN105678085A (en) | PM2.5 concentration estimation method and system | |
| Torbick et al. | Spatiotemporal lake skin summer temperature trends in the Northeast United States | |
| Lei et al. | Observed characteristics of dust storm events over the western United States using meteorological, satellite, and air quality measurements | |
| CN111191673B (en) | Ground surface temperature downscaling method and system | |
| CN111209533B (en) | Road surface fine particulate matter inversion method and device, computing equipment and storage medium | |
| US12032659B2 (en) | Method for identifying dry salt flat based on sentinel-1 data | |
| Huo et al. | Cloud determination of all-sky images under low-visibility conditions | |
| CN109407177A (en) | Dense fog identifying system and methods for using them based on machine learning and conventional meteorological observation | |
| CN104573662B (en) | A kind of cloud sentences method and system | |
| US20120004847A1 (en) | System and method for analyzing trees in lidar data using views | |
| CN114384534A (en) | Tree barrier growth prediction analysis method for unmanned aerial vehicle line patrol | |
| Banks et al. | Retrieval of boundary layer height from lidar using extended Kalman filter approach, classic methods, and backtrajectory cluster analysis | |
| CN113313273B (en) | Public facility detection method, system and storage medium based on big data environment | |
| CN114111705B (en) | Method for calculating cloud base height by using meteorological satellite data | |
| Cavayas et al. | Mapping urban vegetation cover using WorldView-2 imagery | |
| Yang et al. | Radar‐Derived Quantitative Precipitation Estimation Based on Precipitation Classification | |
| Ismanto et al. | Classification tree analysis (CTA) of smoke detection using Himawari_8 satellite data over Sumatera–Borneo Island, Indonesia |
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 |