CN112378828B - Method and device for inverting concentration of atmospheric fine particulate matters based on satellite remote sensing data - Google Patents

Method and device for inverting concentration of atmospheric fine particulate matters based on satellite remote sensing data Download PDF

Info

Publication number
CN112378828B
CN112378828B CN202011453310.0A CN202011453310A CN112378828B CN 112378828 B CN112378828 B CN 112378828B CN 202011453310 A CN202011453310 A CN 202011453310A CN 112378828 B CN112378828 B CN 112378828B
Authority
CN
China
Prior art keywords
data
concentration
atmospheric
remote sensing
fine
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
Application number
CN202011453310.0A
Other languages
Chinese (zh)
Other versions
CN112378828A (en
Inventor
吴剑斌
梁倩
肖林鸿
陈焕盛
秦东明
王文丁
张稳定
陈婷婷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
3Clear Technology Co Ltd
Original Assignee
3Clear Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 3Clear Technology Co Ltd filed Critical 3Clear Technology Co Ltd
Priority to CN202011453310.0A priority Critical patent/CN112378828B/en
Publication of CN112378828A publication Critical patent/CN112378828A/en
Application granted granted Critical
Publication of CN112378828B publication Critical patent/CN112378828B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/06Investigating concentration of particle suspensions
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/17Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/06Investigating concentration of particle suspensions
    • G01N15/075Investigating concentration of particle suspensions by optical means

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Algebra (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Image Processing (AREA)

Abstract

The invention discloses a method and a device for inverting PM2.5 concentration based on satellite remote sensing data. The method includes acquiring raw data sets of a plurality of satellites; preprocessing the acquired original data set to obtain a variable data set of a unified geographical block; performing data fusion of the same type of variable data on the preprocessed variable data sets by using a data fusion model; and respectively matching the geographical position information of the variable data with the geographical position information of the land cover type aiming at various land cover types, monitoring the atmospheric fine particulate matter concentration according to the satellite remote sensing spatial variation coefficient regression model corresponding to the land cover type, and obtaining the atmospheric fine particulate matter concentration of the geographical area where the land cover type is located according to the variable data matched with the geographical position information. The method and the device consider that the PM2.5 concentration is greatly influenced by the underlying surface, and respectively establish an inversion model aiming at different underlying surfaces, so that a PM2.5 data set with high space-time distribution can be obtained.

Description

Method and device for inverting concentration of atmospheric fine particulate matters based on satellite remote sensing data
Technical Field
The invention relates to the technical field of atmospheric satellite remote sensing, in particular to a method and a device for inverting the concentration of atmospheric fine particles based on satellite remote sensing data.
Background
The atmospheric fine particulate matter (PM2.5) refers to particles with aerodynamic diameters of less than or equal to 2.5 μm, is an important component of atmospheric pollutants, and is also a main cause of haze. The PM2.5 has small particle size, wide distribution range and strong particle activity, is easy to attach toxic and harmful substances such as heavy metals and microorganisms, is easy to stay and gather in the atmosphere for a long time, and has great influence on the health of a human body and the quality of the atmospheric environment.
The PM2.5 concentration monitoring mode mainly comprises ground monitoring and satellite remote sensing monitoring. Recently, China has built a plurality of PM2.5 concentration monitoring networks to monitor fine particulate matters in real time. Although ground monitoring is accurate, due to the fact that observation stations are discrete and uneven in distribution, the data coverage degree is still low, and enough data cannot be acquired to study the PM2.5 pollution condition of the whole area.
At present, how to acquire PM2.5 observation data with high resolution, consistent quality and wide coverage range becomes a research focus. The high time resolution characteristic of the satellite remote sensing data can effectively make up the defects of the ground station. However, a method for systematically and comprehensively monitoring the concentration of PM2.5 particles by using a remote sensing means is not popularized yet. The most studied method in recent years is to invert the distribution situation of the ground PM2.5 through the atmospheric Aerosol Optical Depth (AOD) of a satellite remote sensing data product. The data for inverting PM2.5 is divided into the following steps: obtaining AOD data, establishing an inversion ground PM2.5 concentration model and optimizing the model. The most critical of which is the selection of AOD source data.
Research shows that the optical thickness of the aerosol obtained by satellite inversion is closely related to the concentration of particulate matters with the particle size of 0.1-2 microns, and most of PM2.5 is in the particle size range, so that the satellite AOD product is an effective tool for inverting the concentration of the PM2.5 theoretically.
Considering that a satellite sensor can obtain an Optical characteristic parameter of the atmospheric Aerosol, namely the Optical Aerosol thickness (AOD), remote sensing monitoring of air quality in a large range by using a satellite becomes a feasible technical means and also becomes a hotspot for research of broad students. In order to estimate the aerosol property under air pollution by using the optical thickness parameter measured by the satellite, the university scholars of NASA and Maryland (Kaufman and Fraser, 1983) constructs a relation between the mass concentration of aerosol particles and the satellite remote sensing AOD through theoretical analysis and observation experiments.
Successful emission from moderate resolution spectral imager MODIS sensors (mounted on Terra, emitted in 1999, and Aqua satellites, emitted in 2002) provides a wealth of data for the application of satellite telemetry to aerosol monitoring. The scholars of the university of yalaba, usa (wang and Christopher, 2003) have constructed a relation between the optical thickness of a satellite-measured aerosol and the mass concentration of near-surface aerosol particles on the basis of a study of the optical properties of aerosol particles:
τ=f(RH)Qdext(0)Mdaer(0)Heff
Figure BDA0002832319640000021
wherein τ represents the optical thickness of the whole layer of the aerosol; (rh) denotes an aerosol hygroscopy growth factor; RH represents ambient air relative humidity; qdext(0) Denotes the extinction coefficient in m of the near-surface unit mass concentration dry (RH is less than or equal to 40%) aerosol particles2/μg;Mdaer(0) Denotes the mass concentration of the near-surface "dry" aerosol particles in μ g/m3;HeffThe equivalent aerosol elevation is expressed in m; beta is aext(z) represents the extinction coefficient of the aerosol at the height z in m-1(ii) a TOA represents the height of the Top of the Atmosphere (Top of Atmosphere) in m.
The formula provides a theoretical basis for satellite remote sensing monitoring of near-ground aerosol particle concentration, and a subsequent near-ground particle mass concentration estimation method is developed on the basis. On the basis of the theory, three parameters of moisture absorption growth factors, extinction coefficients of aerosol particles with unit mass concentration and aerosol elevation must be obtained for quantitatively obtaining the mass concentration of the aerosol particles near the ground from data of satellite remote sensing spectrum.
The statistical model mainly comprises a univariate simple linear model, a multivariate linear regression model, a multivariate nonlinear statistical model, a multivariate spatial statistical model and the like. The empirical physical mechanism model mainly comprises an empirical double-correction model and a numerical simulation correction model.
Currently used satellite products mainly include medium resolution imaging spectrometers (MODIS) AODs, multi-angle imaging spectrometers (MISR) AODs, Visible Infrared Imaging Radiometers (VIIRS) AODs, and Hiwari-8 AODs.
The PM2.5 remote sensing inversion statistical model mainly uses three parameters of aerosol optical thickness (AOD), boundary layer height HPBL and specific humidity RH to obtain the near-ground PM2.5 concentration. However, historical data of boundary layer height HPBL and specific humidity RH are mostly acquired by ground station observation, and the method has the defect of low data coverage degree and limits PM2.5 remote sensing inversion in a large-range area.
In the prior art, the most common AOD data product of the method for inverting the PM2.5 concentration is AOD data monitored by a MODIS (model Resolution Imaging spectrometer) sensor carried on a Terra & Aqua polar orbit satellite, and a stationary satellite is rarely researched as a data source. Through the development of the last 10 years, the advantages of a geostationary satellite that the same target area is observed at high frequency and a rapidly changing weather system is captured can be utilized, the geostationary satellite AOD data product is used as source data, polar orbit satellite AOD data is updated and supplemented, and the method for efficiently and accurately monitoring the PM2.5 concentration in the specified monitoring area can be provided. When the ground PM2.5 concentration is inverted, satellite AOD data with various resolutions and sources can be selected for use, AOD data can be selected from space and time scales or various AOD data can be combined for use.
Disclosure of Invention
The invention innovatively provides a method and a device for inverting the concentration of atmospheric fine particulate matters based on satellite remote sensing data, and solves the problem that the conventional satellite PM2.5 inversion needs to depend on a large amount of ground station observation data.
In order to achieve the technical purpose, on one hand, the invention discloses a method for inverting the concentration of atmospheric fine particles based on satellite remote sensing data. The method for inverting the concentration of the atmospheric fine particulate matters based on the satellite remote sensing data comprises the following steps: acquiring an original data set of a plurality of satellites; preprocessing the acquired original data set to obtain a variable data set of a unified geographical block; performing data fusion of the same type of variable data on the preprocessed variable data set by using a data fusion model of each type of variable data; respectively matching the geographical position information of variable data with the geographical position information of the land cover type aiming at various land cover types, monitoring a space change coefficient regression model of the atmosphere fine particulate matter concentration according to satellite remote sensing corresponding to the land cover types, and obtaining the atmosphere fine particulate matter concentration of a geographical area where the land cover type is located according to the variable data matched with the geographical position information.
Further, for the method for inverting the concentration of the atmospheric fine particulate matter based on the satellite remote sensing data, the regression model of the spatial variation coefficient of the satellite remote sensing monitoring atmospheric fine particulate matter concentration is obtained through the following processes: collecting hourly data of the concentration of the fine atmospheric particulates, data of all factors influencing the concentration of the fine atmospheric particulates and land coverage type data, and matching geographical position information of the concentration data of the fine atmospheric particulates, geographical position information of the data of all factors influencing the concentration of the fine atmospheric particulates and geographical position information of the land coverage type; and aiming at each land cover type, establishing a training data set on the basis of historical data of the concentration of the fine air particles and historical data of various factors influencing the concentration of the fine air particles, and training by adopting a neural network algorithm to obtain a spatial variation coefficient regression model applicable to different land cover types.
Further, for the method for inverting the concentration of the atmospheric fine particulate matter based on the satellite remote sensing data, the preprocessing is performed on the obtained original data set, and the method comprises the following steps: and matching the variable data to a unified geographical block by at least one interpolation method of linear interpolation, polynomial interpolation and spline curve interpolation.
Further, for the method for inverting the concentration of the atmospheric fine particulate matters based on the satellite remote sensing data, the geographic area is a longitude and latitude pixel grid, and the geographic position information is longitude and latitude information.
Further, for the method for inverting the concentration of the atmospheric fine particulate matter based on the satellite remote sensing data, the preprocessing is performed on the obtained original data set, and the method comprises the following steps: and performing quality control on the variable data according to at least one quality control index in the uncertainty and the reliability, and eliminating data which do not meet the quality control standard.
Further, for the method for inverting the concentration of the atmospheric fine particulate matters based on the satellite remote sensing data, the data fusion model is obtained through the following processes: and distributing proportional weights to the variable data of the same type of different satellites, and establishing a fusion relation model between the variable data of the same type of different satellites.
Further, for the method for inverting the concentration of the atmospheric fine particulate matter based on the satellite remote sensing data, the raw data set comprises aerosol optical thickness data, cloud cover data, atmospheric pollutant data and land cover type data.
In order to achieve the technical purpose, the invention discloses a device for inverting the concentration of atmospheric fine particles based on satellite remote sensing data. The device for inverting the concentration of the atmospheric fine particulate matters based on the satellite remote sensing data comprises: a data acquisition unit for acquiring raw data sets of a plurality of satellites; the preprocessing unit is used for preprocessing the acquired original data set to obtain a variable data set of a unified geographic area; the data fusion unit is used for performing data fusion of the same type of variable data on the preprocessed variable data sets by using the data fusion model of each type of variable data; and the atmosphere fine particle concentration calculation unit is used for matching the geographical position information of variable data with the geographical position information of the land cover type respectively aiming at various land cover types, monitoring a space change coefficient regression model of the atmosphere fine particle concentration according to satellite remote sensing corresponding to the land cover type, and obtaining the atmosphere fine particle concentration of a geographical area where the land cover type is located according to the variable data matched with the geographical position information.
To achieve the above technical object, in yet another aspect, the present invention discloses a computing device. The computing device includes: one or more processors, and a memory coupled with the one or more processors, the memory storing instructions that, when executed by the one or more processors, cause the one or more processors to perform the above-described method.
To achieve the above technical objects, in yet another aspect, the present invention discloses a machine-readable storage medium. The machine-readable storage medium stores executable instructions that, when executed, cause the machine to perform the above-described method.
The invention has the beneficial effects that:
the method and the device for inverting the concentration of the atmospheric fine particulate matters based on the satellite remote sensing data provided by the embodiment of the invention adopt observation data of a plurality of satellites, and an algorithm uses a neural network to establish a correlation relation between influencing factors and considers distinguishing different land cover types. Independent variable factors are comprehensive and easy to obtain, and in order to identify regional clouds, fog and high-reflectivity ground objects, the embodiment of the invention considers the variables which are not considered in the past and can be obtained by satellite remote sensing, such as clouds and sulfur dioxide (SO)2) And/or nitrogen dioxide (NO)2) And adding PM2.5 inversion model. The problem that the conventional satellite PM2.5 inversion needs to depend on a large amount of ground station observation data is solved.
The method and the device for inverting the concentration of the atmospheric fine particulate matters based on the satellite remote sensing data provided by the embodiment of the invention provide a new method for establishing a PM2.5 inversion model by using satellite data, and a PM2.5 data set with high space-time distribution can be obtained. Considering that the PM2.5 concentration is greatly influenced by the underlying surface, the method respectively establishes an inversion model for different underlying surfaces.
Drawings
In the figure, the position of the upper end of the main shaft,
FIG. 1 is a flowchart of a method for inverting the concentration of atmospheric fine particulate matter based on satellite remote sensing data according to an embodiment of the present invention;
FIG. 2 is a flow chart of a multi-source satellite remote sensing monitoring PM2.5 concentration calculation according to an example of the present invention;
FIG. 3 is a flowchart of the preprocessing method of step S120 in FIG. 1 according to an exemplary embodiment of the present invention;
fig. 4 is a schematic structural diagram of an apparatus for inverting the concentration of atmospheric fine particulate matter based on satellite remote sensing data according to embodiment 3 of the present invention;
fig. 5 is a block diagram of a computing device for inversion processing of atmospheric fine particulate matter concentration based on satellite remote sensing data according to an embodiment of the invention.
Detailed Description
The method and the device for inverting the concentration of the atmospheric fine particulate matters based on the satellite remote sensing data provided by the invention are explained and explained in detail below with reference to the attached drawings of the specification.
In order to solve the problem that the conventional satellite PM2.5 inversion needs a large amount of observation data depending on ground stations, the embodiment of the invention integrates the advantages of all satellite detectors and an algorithm based on multi-source satellite data information to obtain satellite AOD data, satellite cloud cover (cloud cover) data and atmospheric pollutant data (which can include SO)2Concentration and/or NO2Concentration), water vapor and/or land cover type and the like are used as source data, modeling is carried out by utilizing observation station historical data and using a neural network algorithm, and satellite remote sensing monitoring PM2.5 concentration spatial data are obtained through calculation.
Fig. 1 is a flowchart of a method for inverting the concentration of atmospheric fine particulate matters based on satellite remote sensing data according to an embodiment of the present invention. Fig. 2 is a flow chart of calculating the concentration of PM2.5 by remote sensing monitoring using a multi-source satellite according to an example of the present invention.
As shown in fig. 1 and 2, raw data sets of a plurality of satellites are acquired at step S110. The raw data set may include, among other things, aerosol optical thickness data, cloud cover data, atmospheric pollutants data, and land cover type data. Atmospheric pollutantsThe data may include sulfur dioxide (SO)2) Concentration and/or nitrogen dioxide (NO)2) Concentration data. As an alternative, the plurality of satellites herein may be a plurality of different types of satellites.
As a more specific example, aerosol optical thickness (AOD) data and cloud (clouded cover) data may be obtained from Himapari-8/9 geostationary satellite, sunflower satellite, SuomiNPP-VIIRS polar orbit, and atmospheric pollutants (which may include SO) data from Sentinel-5P TROPOMI polar orbit2Concentration and/or NO2Concentration, etc.) from Terra&The Aqua/MODIS polar orbit satellite acquires aerosol optical thickness (AOD) data, cloud cover data, land cover type data, and the like, in order to identify regional clouds, fog, and high reflectivity terrain.
In step S120, the obtained original data set is preprocessed to obtain a variable data set of a unified geographic region. The geographic zone may be a grid of latitude and longitude pixels.
Specifically, as shown in fig. 3, step S120 may include the steps of:
and step S122, performing quality control on the variable data according to at least one quality control index in the uncertainty and the credibility, and removing the data which do not accord with the quality control standard. Qualifying satellite data may include: screening quality control factors according to quality control indexes such as Uncertainty (uncertaintiy) and/or credibility (QA flag), and the like, and carrying out aerosol optical thickness (AOD) data, cloud cover (cloud cover) data and nitrogen dioxide (NO) on a plurality of satellites2) Concentration, sulfur dioxide (SO)2) And (3) performing quality control on the variable data such as the concentration and the like, and rejecting the data with poor quality interfered by factors such as cloud and/or surface topography and the like.
And step S124, matching the variable data to a unified geographical block by at least one interpolation method of linear interpolation, polynomial interpolation and spline curve interpolation. And (5) unifying longitude and latitude projections. As a more specific example, the plurality of satellites obtained in step S110 or step S122 can be interpolated by linear interpolation, polynomial interpolation, and/or spline curve interpolationAerosol optical thickness (AOD) data, cloud cover (clouded cover) data, nitrogen dioxide (NO)2) Concentration, sulfur dioxide (SO)2) Matching the variable data such as concentration and the like to a uniform longitude and latitude pixel grid to obtain AOD, cloud cover and NO in the research area2Concentration and SO2And (4) concentration.
In step S130, the data fusion model of each type of variable data is used to perform data fusion of the same type of variable data on the preprocessed variable data set. The data fusion model can be obtained through the following processes: and distributing proportional weights to the variable data of the same type of different satellites, and establishing a fusion relation model between the variable data of the same type of different satellites. In addition, the land cover type may be three-dimensional data composed of land categories and longitude and latitude information, and different land categories may be represented by different values, such as 0, 1, 2, 3, etc., which respectively represent different land categories. When a fusion relation model between variable data of the same type of different satellites is established, areas belonging to the same land coverage type are extracted, and a fusion relation model corresponding to the land coverage type is established in the area.
As a more specific example, for the aerosol optical thickness (AOD) data, the AOD data of a sunflower Hiwari-8/9 satellite, a SuomiNPP-VIIRS satellite and a Terra & Aqua/MODIS polar orbit satellite can be subjected to fusion processing, proportional weights are distributed, a relation model between the three types of satellite aerosol optical thicknesses is established, and the fusion AOD data of uniform longitude and latitude pixel grid points is obtained through inversion. The formula is as follows:
AOD=f(AOD1,AOD2,AOD3)=a1AOD1+a2AOD2+a3AOD3
a1+a2+a3=1
wherein, AOD1Aerosol optical thickness (AOD) for sunflower Hiwari-8/9 satellite, AOD2AOD for SuomiNPP-VIIRS satellites, AOD3Is Terra&AOD of Aqua/MODIS polar orbit satellite, a1、a2And a3Are weight coefficients.
As a more specific example, for the cloud amount data, the weight proportion can be distributed to the cloud amount data of a sunflower Hiwari-8/9 satellite and a SuomiNPP-VIIRS satellite, a relation model between the two satellite cloud amounts is established, and the unified longitude and latitude pixel grid point fusion cloud amount is obtained through inversion. The formula is as follows:
CLOUD=f(CLOUD1,CLOUD2)=b1CLOUD1+b2CLOUD2
b1+b2=1
wherein CLOUD1Is CLOUD, CLOUD of Hiwari-8/9 satellite2Cloud cover for SuomiNPP-VIIRS satellites, b1And b2Are weight coefficients.
In step S140, the geographical position information of the variable data is matched with the geographical position information of the land cover type for each land cover type, the atmospheric fine particle concentration of the geographical area in which the land cover type is located is obtained from the variable data matched with the geographical position information according to the spatial variation coefficient regression model for monitoring the atmospheric fine particle concentration by satellite remote sensing corresponding to the land cover type. The geographic location information may be latitude and longitude information.
The space variation coefficient regression model for monitoring the concentration of the atmospheric fine particulate matters through satellite remote sensing can be obtained through the following processes: collecting the time-by-time data of the concentration of the fine atmospheric particulates, the data of each factor influencing the concentration of the fine atmospheric particulates and the data of the land coverage type, and matching the geographical position information of the concentration data of the fine atmospheric particulates, the geographical position information of the data of each factor influencing the concentration of the fine atmospheric particulates and the geographical position information of the land coverage type; and aiming at each land cover type, establishing a training data set on the basis of historical data of the concentration of the atmospheric fine particulate matters and historical data of various factors influencing the concentration of the atmospheric fine particulate matters, and training by adopting a neural network algorithm to obtain a spatial variation coefficient regression model applicable to different land cover types. After a preset time, the space change coefficient regression model for monitoring the concentration of the fine atmospheric particulates by satellite remote sensing can be updated according to newly acquired concentration data of the fine atmospheric particulates, data of various factors influencing the concentration of the fine atmospheric particulates and/or land cover type data.
As a more specific example, the method comprises the steps of collecting the PM2.5 hourly data of the ground stations of the China region, the data of various factors influencing the concentration of the atmospheric fine particulate matters and the Terra of the China region published by the central environmental monitoring station in the last 3 years of 2017 and 2019&And matching the data of the land coverage type of the Aqua/MODIS polar orbit satellite with the longitude and latitude information of the observation station, the longitude and latitude information of the data of all factors influencing the concentration of the atmospheric fine particulate matters and the longitude and latitude information of the land coverage type, and classifying the land coverage type of the position of the observation station. The MODIS land cover types can be classified using a variety of criteria. For example, the schemes of extracting leaf area index/photosynthetically active radiation component (LAI/fPAR) according to satellite remote sensing can be divided into 9 types, which are water area, cereal and herbaceous crops, shrubs, broadleaf crops, spars grassland, broadleaf forest, coniferous forest, vegetation-free coverage area and cities respectively. The types of land cover may also be classified into 17 categories as defined by the international terrestris biosphere (IGBP) program. Aiming at each type of land cover type, a training data set is established on the basis of historical data of an observation station and data of various factors influencing the concentration of the atmospheric fine particulate matter, and a neural network algorithm is adopted to obtain a function relation applicable to different land cover types through training. As an alternative embodiment, the data for various factors affecting the concentration of atmospheric fine particulate matter may include aerosol optical thickness (AOD) data, cloud cover (clouded cover) data, NO2Concentration, and SO2Satellite data of concentration. The AOD data, cloud cover data, and NO can be obtained by Principal Component Analysis (PCA)2Concentration and SO2The weight coefficients of the influence of the concentration on PM2.5 are respectively c1、c2、c3And c4. Specifically, the formula of the functional relationship may be as follows:
PM2.5L=fL(AOD,CLOUD,NO2,SO2)=c1f(AOD)+c2f(CLOUD)+c3f(NO2)+c4f(SO2)
c1+c2+c3+c4=1
wherein L is a land cover type, AOD and CLOOD are respectively AOD concentration and CLOUD cover amount obtained by fusing in step S130, NO2 and SO2 are respectively NO2 concentration and SO2 concentration obtained by preprocessing in step S120, and c1、c2、c3And c4Are weight coefficients.
And matching longitude and latitude geographic position information of the data according to each land coverage type, and fitting the PM2.5 concentration at each longitude and latitude pixel grid point according to the spatial variation coefficient regression model of the PM2.5 concentration monitored by satellite remote sensing obtained in the step S140.
Fig. 4 is a schematic structural diagram of an apparatus for inverting the concentration of atmospheric fine particulate matters based on satellite remote sensing data according to another embodiment of the present invention. As shown in fig. 4, the apparatus 400 for inverting the concentration of atmospheric fine particulate matter based on satellite remote sensing data provided by this embodiment includes a data acquisition unit 410, a preprocessing unit 420, a data fusion unit 430, and an atmospheric fine particulate matter concentration calculation unit 440.
The data acquisition unit 410 is used to acquire raw data sets for a plurality of satellites. The operation of the data acquisition unit 410 may refer to the operation of step S110 described above with reference to fig. 1.
The preprocessing unit 420 is configured to preprocess the acquired original data set to obtain a variable data set of a unified geographic area. The operation of the preprocessing unit 420 may refer to the operation of step S120 described above with reference to fig. 1.
The data fusion unit 430 is configured to perform data fusion of the same type of variable data on the preprocessed variable data sets by using a data fusion model of each type of variable data. The operation of the data fusion unit 430 may refer to the operation of step S130 described above with reference to fig. 1.
The atmosphere fine particle concentration calculation unit 440 is configured to match geographical position information of the variable data with geographical position information of the land cover type for each land cover type, obtain the atmosphere fine particle concentration of a geographical area where the land cover type is located from the variable data matched with the geographical position information according to a spatial variation coefficient regression model for monitoring the atmosphere fine particle concentration by satellite remote sensing corresponding to the land cover type. The operation of the atmospheric fine particulate matter concentration calculation unit 440 may refer to the operation of step S140 described above with reference to fig. 1.
Fig. 5 is a block diagram of a computing device for inversion processing of atmospheric fine particulate matter concentration based on satellite remote sensing data according to an embodiment of the invention.
As shown in fig. 5, computing device 500 may include at least one processor 510, memory 520, memory 530, communication interface 540, and internal bus 550, and at least one processor 510, memory 520, memory 530, and communication interface 540 are connected together via bus 550. The at least one processor 510 executes at least one computer-readable instruction (i.e., an element described above as being implemented in software) stored or encoded in a computer-readable storage medium (i.e., memory 520).
In one embodiment, stored in the memory 520 are computer-executable instructions that, when executed, cause the at least one processor 510 to: acquiring an original data set of a plurality of satellites; preprocessing the acquired original data set to obtain a variable data set of a unified geographical block; performing data fusion of the same type of variable data on the preprocessed variable data set by using a data fusion model of each type of variable data; respectively matching the geographical position information of variable data with the geographical position information of the land cover type aiming at various land cover types, monitoring a space change coefficient regression model of the atmosphere fine particulate matter concentration according to satellite remote sensing corresponding to the land cover types, and obtaining the atmosphere fine particulate matter concentration of a geographical area where the land cover type is located according to the variable data matched with the geographical position information.
It should be understood that the computer-executable instructions stored in the memory 520, when executed, cause the at least one processor 510 to perform the various operations and functions described above in connection with fig. 1-4 in the various embodiments of the present disclosure.
In the present disclosure, computing device 500 may include, but is not limited to: personal computers, server computers, workstations, desktop computers, laptop computers, notebook computers, mobile computing devices, smart phones, tablet computers, cellular phones, Personal Digital Assistants (PDAs), handheld devices, messaging devices, wearable computing devices, consumer electronics, and so forth.
According to one embodiment, a program product, such as a non-transitory machine-readable medium, is provided. A non-transitory machine-readable medium may have instructions (i.e., elements described above as being implemented in software) that, when executed by a machine, cause the machine to perform various operations and functions described above in connection with fig. 1-4 in various embodiments of the present disclosure.
Specifically, a system or apparatus may be provided which is provided with a readable storage medium on which software program code implementing the functions of any of the above embodiments is stored, and causes a computer or processor of the system or apparatus to read out and execute instructions stored in the readable storage medium.
In this case, the program code itself read from the readable medium can realize the functions of any of the above-described embodiments, and thus the machine-readable code and the readable storage medium storing the machine-readable code form part of the present invention.
Examples of the readable storage medium include floppy disks, hard disks, magneto-optical disks, optical disks (e.g., CD-ROMs, CD-R, CD-RWs, DVD-ROMs, DVD-RAMs, DVD-RWs), magnetic tapes, nonvolatile memory cards, and ROMs. Alternatively, the program code may be downloaded from a server computer or from the cloud via a communications network.
The above description is only an embodiment of the present invention, and not intended to limit the scope of the claims, and all equivalent structures or equivalent processes that are transformed by the content of the specification and the drawings, or directly or indirectly applied to other related technical fields are included in the scope of the claims.

Claims (9)

1. A method for inverting the concentration of atmospheric fine particles based on satellite remote sensing data is characterized by comprising the following steps:
acquiring an original data set of a plurality of satellites;
preprocessing the acquired original data set to obtain a variable data set of a unified geographical block;
performing data fusion of the same type of variable data on the preprocessed variable data set by using a data fusion model of each type of variable data;
respectively matching the geographical position information of variable data with the geographical position information of the land cover type aiming at various land cover types, monitoring a space change coefficient regression model of the concentration of the atmospheric fine particulate matters according to satellite remote sensing corresponding to the land cover types, and obtaining the concentration of the atmospheric fine particulate matters in a geographical area where the land cover type is located according to the variable data matched with the geographical position information; the space variation coefficient regression model for monitoring the concentration of the atmospheric fine particulate matters through satellite remote sensing is obtained through the following processes:
collecting hourly data of the concentration of the fine atmospheric particulates, data of all factors influencing the concentration of the fine atmospheric particulates and land coverage type data, and matching geographical position information of the concentration data of the fine atmospheric particulates, geographical position information of the data of all factors influencing the concentration of the fine atmospheric particulates and geographical position information of the land coverage type;
and aiming at each land cover type, establishing a training data set on the basis of historical data of the concentration of the fine air particles and historical data of various factors influencing the concentration of the fine air particles, and training by adopting a neural network algorithm to obtain a spatial variation coefficient regression model applicable to different land cover types.
2. The method for inverting the concentration of the atmospheric fine particulate matter based on the satellite remote sensing data according to claim 1, wherein preprocessing is performed on the acquired raw data set, and comprises:
and matching the variable data to a unified geographical block by at least one interpolation method of linear interpolation, polynomial interpolation and spline curve interpolation.
3. The method for inverting the concentration of the atmospheric fine particulate matter based on the satellite remote sensing data according to any one of claims 1-2, wherein the geographic area is a latitude and longitude pixel grid, and the geographic position information is latitude and longitude information.
4. The method for inverting the concentration of atmospheric fine particulate matter based on satellite remote sensing data according to claim 1 or 2, wherein preprocessing is performed on the acquired raw data set, and comprises:
and performing quality control on the variable data according to at least one quality control index in the uncertainty and the reliability, and eliminating data which do not meet the quality control standard.
5. The method for inverting the concentration of the atmospheric fine particulate matter based on the satellite remote sensing data according to claim 1, wherein the data fusion model is obtained through the following processes:
and distributing proportional weights to the variable data of the same type of different satellites, and establishing a fusion relation model between the variable data of the same type of different satellites.
6. The method for inverting the concentration of atmospheric fine particulate matter based on satellite remote sensing data according to claim 1, wherein the raw data set comprises aerosol optical thickness data, cloud cover data, atmospheric pollutant data, and land cover type data.
7. A device based on thin particulate matter concentration of satellite remote sensing data retrieval, its characterized in that includes:
a data acquisition unit for acquiring raw data sets of a plurality of satellites;
the preprocessing unit is used for preprocessing the acquired original data set to obtain a variable data set of a unified geographic area;
the data fusion unit is used for performing data fusion of the same type of variable data on the preprocessed variable data sets by using the data fusion model of each type of variable data;
the atmosphere fine particle concentration calculation unit is used for matching geographical position information of variable data with geographical position information of various land cover types respectively according to various land cover types, monitoring a space change coefficient regression model of the atmosphere fine particle concentration according to satellite remote sensing corresponding to the land cover types, and obtaining the atmosphere fine particle concentration of a geographical area where the land cover type is located according to the variable data matched with the geographical position information;
the space variation coefficient regression model for monitoring the concentration of the atmospheric fine particulate matters through satellite remote sensing is obtained through the following processes:
collecting hourly data of the concentration of the fine atmospheric particulates, data of all factors influencing the concentration of the fine atmospheric particulates and land coverage type data, and matching geographical position information of the concentration data of the fine atmospheric particulates, geographical position information of the data of all factors influencing the concentration of the fine atmospheric particulates and geographical position information of the land coverage type;
and aiming at each land cover type, establishing a training data set on the basis of historical data of the concentration of the fine air particles and historical data of various factors influencing the concentration of the fine air particles, and training by adopting a neural network algorithm to obtain a spatial variation coefficient regression model applicable to different land cover types.
8. A computing device, comprising:
one or more processors, and
a memory coupled with the one or more processors, the memory storing instructions that, when executed by the one or more processors, cause the one or more processors to perform the method of any of claims 1-6.
9. A machine-readable storage medium having stored thereon executable instructions that, when executed, cause the machine to perform the method of any one of claims 1 to 6.
CN202011453310.0A 2020-12-11 2020-12-11 Method and device for inverting concentration of atmospheric fine particulate matters based on satellite remote sensing data Active CN112378828B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011453310.0A CN112378828B (en) 2020-12-11 2020-12-11 Method and device for inverting concentration of atmospheric fine particulate matters based on satellite remote sensing data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011453310.0A CN112378828B (en) 2020-12-11 2020-12-11 Method and device for inverting concentration of atmospheric fine particulate matters based on satellite remote sensing data

Publications (2)

Publication Number Publication Date
CN112378828A CN112378828A (en) 2021-02-19
CN112378828B true CN112378828B (en) 2021-09-17

Family

ID=74591002

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011453310.0A Active CN112378828B (en) 2020-12-11 2020-12-11 Method and device for inverting concentration of atmospheric fine particulate matters based on satellite remote sensing data

Country Status (1)

Country Link
CN (1) CN112378828B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113554305B (en) * 2021-07-22 2022-04-22 生态环境部卫星环境应用中心 Method for evaluating influence of straw burning on concentration of atmospheric fine particulate matters
CN113654959B (en) * 2021-07-29 2023-11-17 中国科学院合肥物质科学研究院 Rapid inversion method and system for smoke cloud concentration space-time distribution
CN113804595B (en) * 2021-09-15 2024-04-05 汉威科技集团股份有限公司 Multi-parameter air quality monitoring system
CN114091605A (en) * 2021-11-23 2022-02-25 云南电网有限责任公司电力科学研究院 Fusion method and device based on current dirty area map and remote sensing data
CN114707396B (en) * 2021-12-28 2024-08-20 华东师范大学 Near real-time production method of seamless grid point data of PM2.5 concentration in whole day
CN114414446B (en) * 2021-12-29 2023-10-20 中国科学院空天信息创新研究院 Near-ground particulate matter concentration estimation method, device and equipment based on satellite remote sensing
CN116698688B (en) * 2023-04-20 2024-08-20 兰州大学 Method for estimating concentration of atmospheric particulates based on double-star of cloud number 4

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000131218A (en) * 1998-10-26 2000-05-12 Ntt Data Corp Method and device for estimating concentration of substance in water, and record medium
CN103674794A (en) * 2013-12-16 2014-03-26 中国科学院遥感与数字地球研究所 Multivariable regression method for remote sensing monitoring of near-surface fine particle matter PM2.5 mass concentration
CN106442236A (en) * 2015-07-30 2017-02-22 中国科学院遥感与数字地球研究所 Ground PM2.5 inversion method and system based on satellite remote sensing
CN108491667A (en) * 2018-04-13 2018-09-04 武汉大学 A kind of PM2.5 remote sensing monitoring inversion methods based on Himawar-8 AOD
CN108761484A (en) * 2018-04-26 2018-11-06 江苏省气象台 A kind of sea fog monitoring method based on Multi-sensor satellite remote sensing
CN109583516A (en) * 2018-12-24 2019-04-05 天津珞雍空间信息研究院有限公司 A kind of space and time continuous PM2.5 inversion method based on ground and moonscope
KR101967250B1 (en) * 2017-12-29 2019-04-09 한국교원대학교 산학협력단 System for detecting the large-scale transport of haze and method
CN110108609A (en) * 2019-05-07 2019-08-09 国网湖南省电力有限公司 A kind of PM2.5 concentration value distribution situation acquisition methods and system based on satellite multichannel
WO2019173381A1 (en) * 2018-03-06 2019-09-12 Applied Particle Technology, Inc. Wireless exposure monitor
CN110398446A (en) * 2019-08-01 2019-11-01 成都信息工程大学 A kind of PM of combination AOD and meteorologic parameter2.5Concentration approximating method
CN110726653A (en) * 2019-09-25 2020-01-24 中国电子科技集团公司第二十七研究所 PM based on heaven and earth integration information2.5Concentration monitoring method
CN110954482A (en) * 2019-12-02 2020-04-03 生态环境部卫星环境应用中心 Atmospheric pollution gridding monitoring method based on static satellite and polar orbit satellite
WO2020078560A1 (en) * 2018-10-18 2020-04-23 Cellfacts Analytics Limited A method and apparatus for monitoring microbial contaminants in an industrial process
CN111079835A (en) * 2019-12-16 2020-04-28 安徽大学 Himapari-8 atmospheric aerosol inversion method based on deep full-connection network
CN111257241A (en) * 2020-01-20 2020-06-09 中国科学院地理科学与资源研究所 Atmospheric carbon dioxide concentration inversion algorithm based on DEEI (DeEI)
WO2020113275A1 (en) * 2018-12-06 2020-06-11 Newcastle City Council A housing
CN111538940A (en) * 2020-05-12 2020-08-14 自然资源部第二海洋研究所 Suspended matter concentration inversion model determination method and suspended matter concentration determination method
CN111912754A (en) * 2020-07-23 2020-11-10 安徽省气象科学研究所 Remote sensing inversion method for near-surface particulate matter concentration

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106124374A (en) * 2016-07-22 2016-11-16 中科宇图科技股份有限公司 Atmospheric particulates remote-sensing monitoring method based on data fusion
CN109030723A (en) * 2018-06-14 2018-12-18 北京国电高科科技有限公司 atmospheric monitoring system and method

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000131218A (en) * 1998-10-26 2000-05-12 Ntt Data Corp Method and device for estimating concentration of substance in water, and record medium
CN103674794A (en) * 2013-12-16 2014-03-26 中国科学院遥感与数字地球研究所 Multivariable regression method for remote sensing monitoring of near-surface fine particle matter PM2.5 mass concentration
CN106442236A (en) * 2015-07-30 2017-02-22 中国科学院遥感与数字地球研究所 Ground PM2.5 inversion method and system based on satellite remote sensing
KR101967250B1 (en) * 2017-12-29 2019-04-09 한국교원대학교 산학협력단 System for detecting the large-scale transport of haze and method
WO2019173381A1 (en) * 2018-03-06 2019-09-12 Applied Particle Technology, Inc. Wireless exposure monitor
CN108491667A (en) * 2018-04-13 2018-09-04 武汉大学 A kind of PM2.5 remote sensing monitoring inversion methods based on Himawar-8 AOD
CN108761484A (en) * 2018-04-26 2018-11-06 江苏省气象台 A kind of sea fog monitoring method based on Multi-sensor satellite remote sensing
WO2020078560A1 (en) * 2018-10-18 2020-04-23 Cellfacts Analytics Limited A method and apparatus for monitoring microbial contaminants in an industrial process
WO2020113275A1 (en) * 2018-12-06 2020-06-11 Newcastle City Council A housing
CN109583516A (en) * 2018-12-24 2019-04-05 天津珞雍空间信息研究院有限公司 A kind of space and time continuous PM2.5 inversion method based on ground and moonscope
CN110108609A (en) * 2019-05-07 2019-08-09 国网湖南省电力有限公司 A kind of PM2.5 concentration value distribution situation acquisition methods and system based on satellite multichannel
CN110398446A (en) * 2019-08-01 2019-11-01 成都信息工程大学 A kind of PM of combination AOD and meteorologic parameter2.5Concentration approximating method
CN110726653A (en) * 2019-09-25 2020-01-24 中国电子科技集团公司第二十七研究所 PM based on heaven and earth integration information2.5Concentration monitoring method
CN110954482A (en) * 2019-12-02 2020-04-03 生态环境部卫星环境应用中心 Atmospheric pollution gridding monitoring method based on static satellite and polar orbit satellite
CN111079835A (en) * 2019-12-16 2020-04-28 安徽大学 Himapari-8 atmospheric aerosol inversion method based on deep full-connection network
CN111257241A (en) * 2020-01-20 2020-06-09 中国科学院地理科学与资源研究所 Atmospheric carbon dioxide concentration inversion algorithm based on DEEI (DeEI)
CN111538940A (en) * 2020-05-12 2020-08-14 自然资源部第二海洋研究所 Suspended matter concentration inversion model determination method and suspended matter concentration determination method
CN111912754A (en) * 2020-07-23 2020-11-10 安徽省气象科学研究所 Remote sensing inversion method for near-surface particulate matter concentration

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Study on the temporal and spatial distribution of chlorophyll a in Erhai Lake based on multispectral data from environmental satellites";Xingmin Wang et al;《Ecological Informatics》;20201112;第61卷;第1-13页 *
"高分五号卫星多角度偏振成像仪细粒子气溶胶光学厚度遥感反演";谢一凇 等;《上海航天》;20191231;第36卷;第219-226页 *

Also Published As

Publication number Publication date
CN112378828A (en) 2021-02-19

Similar Documents

Publication Publication Date Title
CN112378828B (en) Method and device for inverting concentration of atmospheric fine particulate matters based on satellite remote sensing data
Zhao et al. Reconstruction of daytime land surface temperatures under cloud-covered conditions using integrated MODIS/Terra land products and MSG geostationary satellite data
Schirrmann et al. Monitoring agronomic parameters of winter wheat crops with low-cost UAV imagery
Román et al. Use of in situ and airborne multiangle data to assess MODIS-and Landsat-based estimates of directional reflectance and albedo
US20230316691A1 (en) Method And System For Detecting, Quantifying, And Attributing Gas Emissions Of Industrial Assets
Rosas et al. Sensitivity of Landsat 8 surface temperature estimates to atmospheric profile data: A study using MODTRAN in dryland irrigated systems
Ma et al. A global long-term (1981–2000) land surface temperature product for NOAA AVHRR
CN111738175A (en) Agricultural drought monitoring system based on remote sensing image and convolutional neural network
Carrer et al. Land surface albedo derived on a ten daily basis from Meteosat Second Generation Observations: The NRT and climate data record collections from the EUMETSAT LSA SAF
Dadashazar et al. Cloud drop number concentrations over the western North Atlantic Ocean: seasonal cycle, aerosol interrelationships, and other influential factors
Yu et al. A physical-based atmospheric correction algorithm of unmanned aerial vehicles images and its utility analysis
Bai et al. Comparison of satellite-based PM2. 5 estimation from aerosol optical depth and top-of-atmosphere reflectance
Xu et al. A new algorithm for Himawari-8 aerosol optical depth retrieval by integrating regional PM₂. ₅ concentrations
Mei et al. XBAER-derived aerosol optical thickness from OLCI/Sentinel-3 observation
CN112819026A (en) Fusion method and device for optical thickness data of multisource satellite aerosol
Mishra et al. Predicting ground level PM2. 5 concentration over Delhi using Landsat 8 satellite data
Bao et al. High-spatial-resolution aerosol optical properties retrieval algorithm using Chinese high-resolution earth observation satellite I
Shaohua et al. Advances of ozone satellite remote sensing in 60 years
Dong et al. Retrieval of aerosol single scattering albedo using joint satellite and surface visibility measurements
Bingcheng et al. BRDF feature observation method and modeling of desert site based on UAV platform
Li et al. DeepSAT4D: Deep learning empowers four-dimensional atmospheric chemical concentration and emission retrieval from satellite
Alonso-Montesinos et al. Impact of a Saharan dust intrusion over southern Spain on DNI estimation with sky cameras
Bagheri Using deep ensemble forest for high-resolution mapping of PM2. 5 from MODIS MAIAC AOD in Tehran, Iran
Grzegorski et al. Multi-Sensor Retrieval of Aerosol Optical Properties for Near-Real-Time Applications Using the Metop Series of Satellites: Concept, Detailed Description, and First Validation
Fréville et al. Analysis of the Transport of Aerosols over the North Tropical Atlantic Ocean Using Time Series of POLDER/PARASOL Satellite Data

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20210219

Assignee: Beijing Zhongke Sanqing Environmental Technology Co.,Ltd.

Assignor: 3CLEAR TECHNOLOGY Co.,Ltd.

Contract record no.: X2022980012305

Denomination of invention: Method and device for retrieving atmospheric fine particle concentration based on satellite remote sensing data

Granted publication date: 20210917

License type: Common License

Record date: 20220815

EE01 Entry into force of recordation of patent licensing contract