CN109085100B - Method and device for determining concentration of pollutant - Google Patents

Method and device for determining concentration of pollutant Download PDF

Info

Publication number
CN109085100B
CN109085100B CN201811083073.6A CN201811083073A CN109085100B CN 109085100 B CN109085100 B CN 109085100B CN 201811083073 A CN201811083073 A CN 201811083073A CN 109085100 B CN109085100 B CN 109085100B
Authority
CN
China
Prior art keywords
grid
data
sub
pollutant concentration
concentration data
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
CN201811083073.6A
Other languages
Chinese (zh)
Other versions
CN109085100A (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.)
Beijing Yingshi Ruida Technology Co ltd
Original Assignee
Beijing Yingshi Ruida 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 Beijing Yingshi Ruida Technology Co ltd filed Critical Beijing Yingshi Ruida Technology Co ltd
Priority to CN201811083073.6A priority Critical patent/CN109085100B/en
Publication of CN109085100A publication Critical patent/CN109085100A/en
Application granted granted Critical
Publication of CN109085100B publication Critical patent/CN109085100B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Measurement Of Radiation (AREA)

Abstract

The invention provides a method for determining the concentration of pollutants, which comprises the following steps: acquiring first to third point location information and first to third pollutant concentration data of first to third devices; acquiring all sub-grid information of a first grid; calculating first to third distances of each sub-grid and the first to third devices according to the point location information and the position information of each sub-grid; calculating original pollutant concentration data of each sub-grid according to the first data of the first area and the position information of each grid; taking the third correction data of each sub-grid as the pollutant concentration data of the sub-grid according to the first to third distances; determining original pollutant concentration data of the first grid according to a preset relation between the first grid and all the sub-grid information and the pollutant concentration data of each sub-grid; repeating the steps, and taking the sixth correction data as the pollutant concentration data of the first grid. Thereby, the accuracy of the intra-grid contaminant concentration data is improved.

Description

Method and device for determining concentration of pollutant
Technical Field
The present invention relates to the field of data processing, and in particular, to a method and apparatus for determining a concentration of a contaminant.
Background
With the rapid development of various industries, a large amount of harmful substances such as smoke, sulfur dioxide, nitrogen oxides, carbon monoxide, hydrocarbons, etc. are generated. The harmful substances are continuously discharged into the atmosphere, and when the content exceeds the limit of the environment, natural physical, chemical and ecological balance is destroyed, so that atmospheric pollution is formed, and life, work and health of people are endangered. With the advent of nationwide and wide range of haze weather, the term PM 2.5 is coming into public view. PM 2.5 refers to particulates having an ambient aerodynamic equivalent diameter of 2.5 microns or less. It can be suspended in air for a longer time, and the higher the content concentration of the suspension in the air is, the more serious the air pollution is.
With the rapid development of the economic society, environmental problems become one of important obstacle factors for the development of the society, and solving good environmental problems becomes an urgent problem for various countries.
One of the important bases for solving the environmental problems is to accurately grasp the current environmental situation, including which specific environmental problems exist, etc., and the environmental monitoring work is also the key for solving the environmental problems and knowing the current environmental situation in time, wherein the accuracy of the environmental monitoring data becomes the key point and the key link of the environmental monitoring work.
The environment monitoring data is the basis for formulating environment protection policies and measures, and is also the basis for environment management, law enforcement, statistics, information release and environment protection target responsibility system assessment. Therefore, whether the quality of the environment detection data is positive for the environment protection work.
The atmospheric pollution monitoring is to measure the type and concentration of pollutants in the atmospheric environment and observe the time-space distribution and change rule. The atmospheric pollution monitoring aims to identify pollutant in the atmosphere, master the distribution and diffusion rule of the pollutant, and monitor the emission and control conditions of an atmospheric pollution source. Because the monitoring area is large in range, manpower and material resources are limited, and difficulty is brought to atmospheric pollution monitoring.
However, in the prior art, the monitoring area is generally divided according to an administrative unit, such as a city, so that the obtained pollutant concentration data is specific to the administrative unit, for example, the pollutant concentration data of a certain city, and since the interior of the administrative unit further comprises small administrative units with different pollution degrees, the pollutant concentration data cannot reliably reflect the pollution degrees of the administrative units.
Disclosure of Invention
The embodiment of the invention aims at overcoming the defects in the prior art and provides a method and a device for determining the concentration of pollutants.
To solve the above problems, in a first aspect, the present invention provides a method of determining a concentration of a contaminant, the method comprising:
acquiring first point position information and first pollutant concentration data of at least three first devices in a first grid;
acquiring second point location information and second pollutant concentration data of second equipment;
acquiring third point location information and third pollutant concentration data of third equipment;
acquiring all pieces of sub-grid information of a first grid, wherein each piece of sub-grid information in all pieces of sub-grid information comprises the position information of the sub-grid;
acquiring first data of a first area, wherein the first data comprises position information of the first area and pollutant concentration data of the first area;
calculating original pollutant concentration data of each sub-grid according to the position information of the first area, the pollutant concentration data of the first area and the position information of each sub-grid;
according to the first point position information and the position information of each sub-grid, respectively calculating a first distance between each sub-grid and the first device;
correcting the original pollutant concentration data of each sub-grid according to the first distance and the original pollutant concentration data of each sub-grid to obtain first correction data;
Respectively calculating a second distance between each sub-grid and the second device according to the second point location information and the position information of each sub-grid; correcting the first correction data of each sub-grid according to the second distance and the second pollutant concentration data to obtain second correction data;
according to the third point location information and the position information of each sub-grid, calculating a third distance between each sub-grid and the third device;
correcting the second correction data of each sub-grid according to the third distance and the third pollutant concentration data to obtain third correction data, and taking the third correction data of each sub-grid as pollutant concentration data;
determining original pollutant concentration data of the first grid according to the preset relation between the first grid and all the sub-grid information and the pollutant concentration data of each sub-grid;
acquiring position information of the first grid;
calculating a fourth distance between the first grid and the first device according to the first point position information and the position information of the first grid;
correcting the original pollutant concentration data of the first grid according to the fourth distance and the first pollutant concentration data to obtain fourth correction data;
Calculating a fifth distance between the first grid and the second device according to the second point location information and the position information of the first grid;
correcting the fourth correction data according to the fifth distance and the second pollutant concentration data to obtain fifth correction data;
calculating a sixth distance between the first grid and the third device according to the third point location information and the position information of the first grid;
and correcting the fifth correction data according to the sixth distance and the third pollutant concentration data to obtain sixth correction data, and taking the sixth correction data as the pollutant concentration data of the first grid.
In one possible implementation manner, the calculating the first distance between each sub-grid and the first device according to the first point location information and the position information of each sub-grid specifically includes:
acquiring longitude and latitude data of a first point bit;
acquiring longitude and latitude data of each sub-grid, wherein the longitude and latitude data of each sub-grid comprises upper left longitude and latitude data, lower left longitude and latitude data, upper right longitude and latitude data, lower right longitude and latitude data and middle longitude and latitude data;
And calculating a first distance between each sub-grid and the first device according to the longitude and latitude data of the first point location, the upper left longitude and latitude data, the lower left longitude and latitude data, the upper right longitude and latitude data, the lower right longitude and latitude data and the middle longitude and latitude data.
In one possible implementation manner, the calculating the original pollutant concentration data of each sub-grid according to the position information of the first area, the pollutant concentration data of the first area and the position information of each sub-grid specifically includes:
determining an overlapping area of the first area and the sub-grids according to the position information of the first area and the position information of each sub-grid;
determining contaminant concentration data of the sub-grid of the overlapping region according to the contaminant concentration data of the first region;
and obtaining pollutant concentration data of the subgrid of the non-overlapping area through interpolation, wherein the pollutant concentration data of the subgrid of the overlapping area and the pollutant concentration data of the subgrid of the non-overlapping area form original pollutant concentration data of all subgrids.
In one possible implementation manner, the correcting the original pollutant concentration data of each sub-grid according to the first distance and the original pollutant concentration data of each sub-grid to obtain first corrected data specifically includes:
Determining a correlation coefficient according to a preset radius constant and the position information of the sub-grid;
and obtaining first correction data according to the original pollutant concentration data of each sub-grid, the correlation coefficient and the first distance.
In one possible implementation manner, the determining the original pollutant concentration data of the first grid according to the preset relation between the first grid and the all sub-grid information and the pollutant concentration data of each sub-grid specifically includes:
according to the preset relation between the first grid and all the sub-grid information, determining the number of sub-grids;
determining a sum of the contaminant concentration data from the contaminant concentration data for each sub-grid;
and determining the original pollutant concentration data of the first grid according to the sum of the number of the sub-grids and the pollutant concentration data.
In a second aspect, the present invention provides a pollutant concentration determining apparatus, comprising:
an acquisition unit for acquiring first point location information and first contaminant concentration data of at least three first devices within a first grid;
the acquisition unit is further used for acquiring second point location information and second pollutant concentration data of second equipment;
The acquisition unit is further used for acquiring third point location information and third pollutant concentration data of third equipment;
the acquiring unit is further configured to acquire all pieces of sub-grid information of the first grid, where each piece of sub-grid information in the all pieces of sub-grid information includes location information of the sub-grid;
the acquisition unit is further used for acquiring first data of the first area, wherein the first data comprises position information of the first area and pollutant concentration data of the first area;
a calculation unit, said calculation unit being adapted to,
calculating original pollutant concentration data of each sub-grid according to the position information of the first area, the pollutant concentration data of the first area and the position information of each sub-grid;
the computing unit is further configured to respectively compute a first distance between each sub-grid and the first device according to the first point position information and the position information of each sub-grid;
the correction unit is used for correcting the original pollutant concentration data of each sub-grid according to the first distance and the original pollutant concentration data of each sub-grid to obtain first correction data;
The calculating unit is further used for calculating a second distance between each sub-grid and the second device according to the second point location information and the position information of each sub-grid;
the correction unit is further used for correcting the first correction data of each sub-grid according to the second distance and the second pollutant concentration data to obtain second correction data;
the calculating unit is further configured to calculate a third distance between each sub-grid and the third device according to the third point location information and the position information of each sub-grid;
the correction unit is further configured to correct the second correction data of each sub-grid according to the third distance and the third pollutant concentration data, to obtain third correction data, and use the third correction data of each sub-grid as pollutant concentration data thereof;
the determining unit is used for determining original pollutant concentration data of the first grid according to the preset relation between the first grid and all the sub-grid information and the pollutant concentration data of each sub-grid;
the acquisition unit is further used for acquiring the position information of the first grid;
The computing unit is further configured to calculate a fourth distance between the first grid and the first device according to the first point position information and the position information of the first grid;
the correction unit is further used for correcting the original pollutant concentration data of the first grid according to the fourth distance and the first pollutant concentration data to obtain fourth correction data;
the calculating unit is further used for calculating a fifth distance between the first grid and the second device according to the second point location information and the position information of the first grid;
the correction unit is further configured to correct the fourth correction data according to the fifth distance and the second pollutant concentration data, so as to obtain fifth correction data;
the calculating unit is further configured to calculate a sixth distance between the first grid and the third device according to the third point location information and the position information of the first grid;
the correction unit is further configured to correct the fifth correction data according to the sixth distance and the third pollutant concentration data, obtain sixth correction data, and use the sixth correction data as the pollutant concentration data of the first grid.
In one possible implementation, the computing unit is specifically configured to:
acquiring longitude and latitude data of a first point bit;
acquiring longitude and latitude data of each sub-grid, wherein the longitude and latitude data of each sub-grid comprises upper left longitude and latitude data, lower left longitude and latitude data, upper right longitude and latitude data, lower right longitude and latitude data and middle longitude and latitude data;
and calculating a first distance between each sub-grid and the first device according to the longitude and latitude data of the first point location, the upper left longitude and latitude data, the lower left longitude and latitude data, the upper right longitude and latitude data, the lower right longitude and latitude data and the middle longitude and latitude data.
In one possible implementation, the computing unit is specifically configured to:
determining an overlapping area of the first area and the sub-grids according to the position information of the first area and the position information of each sub-grid;
determining contaminant concentration data of the sub-grid of the overlapping region according to the contaminant concentration data of the first region;
and obtaining pollutant concentration data of the subgrid of the non-overlapping area through interpolation, wherein the pollutant concentration data of the subgrid of the overlapping area and the pollutant concentration data of the subgrid of the non-overlapping area form original pollutant concentration data of all subgrids.
In a possible implementation manner, the correction unit is specifically configured to:
determining a correlation coefficient according to a preset radius constant and the position information of the sub-grid;
and obtaining first correction data according to the original pollutant concentration data of each sub-grid, the correlation coefficient and the first distance.
In a possible implementation manner, the determining unit is specifically configured to:
determining the number of the sub-grids according to the preset relation between the first grid and the whole sub-grid information;
determining a sum of the contaminant concentration data from the contaminant concentration data for each sub-grid;
and determining the original pollutant concentration data of the first grid according to the sum of the number of the sub-grids and the pollutant concentration data.
In a third aspect, the present invention provides an apparatus comprising a memory for storing a program and a processor for performing the method of any of the first aspects.
In a fourth aspect, the present invention provides a computer program product comprising instructions which, when run on a computer, cause the computer to perform the method of any of the first aspects.
In a fifth aspect, the present invention provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements a method according to any of the first aspects.
By applying the method for determining the pollutant concentration, provided by the invention, the interpolation method is utilized to calculate the original pollutant concentration data of the small grids, then the original pollutant concentration data of each sub-grid is sequentially controlled with the first pollutant concentration data of the first equipment, the second pollutant concentration data of the second equipment and the third pollutant concentration data of the third equipment to obtain the pollutant concentration data of the small grids, and then the pollutant concentration data of each small grid is sequentially controlled with the first pollutant concentration data of the first equipment, the second pollutant concentration data of the second equipment and the third pollutant concentration data of the third equipment to obtain the pollutant concentration data of the first grid, so that the area is divided into grids, the pollutant concentration data of the sub-grids in the grids are obtained, a novel method for obtaining the pollutant concentration data is provided, the accuracy of the pollutant concentration data is improved, and the accuracy and the effectiveness of the whole environment monitoring work are improved.
Drawings
FIG. 1 is a flow chart of a method for determining a concentration of a contaminant according to a first embodiment of the present application;
fig. 2 is a schematic diagram of a determining device for determining a concentration of a contaminant according to a second embodiment of the present application.
Detailed Description
The application is described in further detail below with reference to the drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the application and are not limiting of the application. It should be noted that, for convenience of description, only the portions related to the present application are shown in the drawings.
It should be noted that, without conflict, the embodiments of the present application and features of the embodiments may be combined with each other. The application will be described in detail below with reference to the drawings in connection with embodiments.
In order to facilitate a better description of the method according to the application, a description of the "grid" is first provided below.
In order to achieve the purpose of fine control and management of regional atmospheric pollution, a target region is divided into different grids according to different monitoring requirements and environmental characteristics to carry out point location arrangement, and the concentration of relevant pollutants in each grid is monitored in real time, which is called gridding monitoring. The key pollution area divided by the urban meshing supervision work is called as a hot spot mesh. High density meshed monitoring network
Various functional monitoring points are reasonably distributed in the area, so that the air quality change of the key polluted area can be reflected, the requirement of area environment air monitoring is met, and the air quality of the key polluted area is objectively evaluated.
The distribution condition of the pollutants can be evaluated according to the diffusion, migration and conversion rules of the local pollutants, and reasonable monitoring points can be determined by combining the feasibility of resources and economy, so that the obtained data is representative.
At the determined reasonable monitoring point location, a gridding monitoring device can be set. The gridding monitoring equipment is a detection method adopting light scattering, has small volume and light weight, and is used for continuously and automatically monitoring the pollutant condition in the ambient air.
At a site, there is typically a standard monitoring device (also referred to as a national control device or a provincial control device), and at least 3 meshed monitoring devices may be installed in a certain range of the site, where the 3 meshed monitoring devices are referred to as quality control devices, and the three quality control devices form a quality control point.
Besides the quality control point, a plurality of point positions can be provided with grid monitoring equipment. After dividing a relatively large area into grids, the grids can be divided into a plurality of sub-grids, and the sub-grids can be divided into a plurality of small sub-grids so as to improve the accuracy of pollutant monitoring.
It will be appreciated that the concept of a grid is relative. For example, the X market may be divided into 36 grids, and each of the 36 grids may be divided into 36 small grids. The grids involved in the present application may be the 36 grids or the 36 small grids, and in particular, should be determined in practical applications, which is not limited by the present application.
The following first and second etc. numbers are merely for distinction and are not actually meant.
Fig. 1 is a schematic flow chart of a method for determining a concentration of a contaminant according to a first embodiment of the present application. The application scenario of the method is a meshed monitoring network, the execution subject of the method may be a device with a calculation function, for example, a computer, a mobile phone or a pollutant concentration determining device, etc., where the computer, the mobile phone or the pollutant concentration determining device may be connected to the meshed monitoring device, and the connection may be performed by a wireless or wired communication manner, which is not limited in this aspect of the application. As shown in fig. 1, the method comprises the steps of:
step 101, acquiring first point location information and first pollutant concentration data of at least three first devices in a first grid.
In order to obtain more accurate pollutant concentration data, the area may be gridded, and at least three first devices are included in the first grid, where the first devices may be, by way of example and not limitation, gridded monitoring devices (also referred to as quality control devices) of the quality control point, and the number of the first devices is 4.
After the grid monitoring equipment is put into the quality control point, the grid monitoring equipment can acquire pollutant concentration data of the point in real time, and the pollutant concentration data can comprise the type of pollutants and the concentration value of the pollutants under the type. In this case, the "real-time" may be set on the grid-type monitoring device, and may be set as needed, for example, but not limited to, 60 pollutant concentration data may be acquired in one minute, and the more data are acquired in one minute, the more pollution is generated at the subsequent first point location
The more accurate the dye concentration data.
By way of example and not limitation, the contaminants may be any of fine particulate matter (PM 2.5), inhalable particulate matter (PM 10), nitrogen dioxide (NO 2), sulfur dioxide (SO 2), carbon monoxide (CO), ozone (O3), and total volatile organic compounds (Total Volatile Organic Compounds, TVOC).
It will be appreciated that in subsequent studies, the contaminants may be any combination of the above contaminants, and that units of different contaminants may be processed by normalization to obtain normalized contaminant concentration data, which may be used to calculate the contaminant concentration data of the grid after comprehensive determination of the normalized contaminant concentration data.
It will be appreciated that contaminant concentration data herein and hereinafter contaminant concentration data, refer to valid data, for which the default is-999.
Step 102, obtaining second point location information and second pollutant concentration data of a second device.
And step 103, acquiring third point location information and third pollutant concentration data of third equipment.
The second equipment is provincial control equipment, and the third equipment is national control equipment. The second device and the third device may be located in the first mesh or may be located outside the first mesh, and the positional relationship between the second device and the third device and the first mesh is not limited in the present application.
Step 104, acquiring all pieces of sub-grid information of the first grid, wherein each piece of sub-grid information in all pieces of sub-grid information comprises the position information of the sub-grid.
In particular, the first grid may include a number of sub-grids, by way of example and not limitation, and the first grid includes 36 sub-grids.
Step 105, acquiring first data of a first area, wherein the first data comprises position information of the first area and contaminant concentration data of the first area.
The first data may be nc data of WRF, which is data obtained by training using a preset model, and the data is stored in a database. The first region is a region that is larger than each small grid, but smaller than the first grid.
The nc data of the WRF includes positional information of the first region and contaminant concentration data of the first region. The location information may be latitude and longitude data and the contaminant concentration data of the first region may be a plurality.
And 106, calculating the original pollutant concentration data of each sub-grid according to the position information of the first area, the pollutant concentration data of the first area and the position information of each sub-grid.
Step 106 specifically includes: firstly, determining an overlapping area of a first area and each sub-grid according to the position information of the first area and the position information of each sub-grid; then, determining contaminant concentration data of the sub-grid of the overlapping region from the contaminant concentration data of the first region; and finally, obtaining the pollutant concentration data of the subgrid of the non-overlapping area through interpolation, wherein the pollutant concentration data of the subgrid of the overlapping area and the pollutant concentration data of the subgrid of the non-overlapping area form the original pollutant concentration data of all subgrids.
Step 107, calculating the first distance between each sub-grid and the first device according to the first point position information and the position information of each sub-grid.
Step 107 specifically includes:
firstly, acquiring longitude and latitude data of a first point bit; then, acquiring longitude and latitude data of each sub-grid, wherein the longitude and latitude data of each sub-grid comprises upper left longitude and latitude data, lower left longitude and latitude data, upper right longitude and latitude data, lower right longitude and latitude data and middle longitude and latitude data; and finally, calculating the first distance between each sub-grid and the first device according to the longitude and latitude data of the first point location, the upper left longitude and latitude data, the lower left longitude and latitude data, the upper right longitude and latitude data, the lower right longitude and latitude data and the middle longitude and latitude data.
And step 108, correcting the original pollutant concentration data of each sub-grid according to the first distance and the original pollutant concentration data of each sub-grid to obtain first correction data.
Wherein, the first correction data can be obtained by using an inverse distance weight interpolation method.
Specifically, firstly, determining a correlation coefficient according to a preset radius constant and position information of a sub-grid; and finally, obtaining first correction data according to the original pollutant concentration data, the correlation coefficient and the first distance of each sub-grid.
Specifically, dis, dis=sqrt (dx+dy) may be calculated according to the longitude and latitude of each sub-grid.
Then, the correlation coefficient can be calculated from the influence radius constants R and dis.
Correlation coefficient= (R-dis) 2/(r+dis) 2×0.8.
Finally, first correction data = raw contaminant concentration data + correlation coefficient for each sub-grid (first contaminant concentration data-first distance).
Step 109, calculating the second distance between each sub-grid and the second device according to the second point location information and the position information of each sub-grid.
Step 110, the first correction data of each sub-grid is corrected according to the second distance and the second pollutant concentration data, so as to obtain second correction data.
Step 111, calculating a third distance between each sub-grid and the third device according to the third point location information and the position information of each sub-grid.
And step 112, correcting the second correction data of each sub-grid according to the third distance and the third pollutant concentration data to obtain third correction data, and taking the third correction data of each sub-grid as the pollutant concentration data.
Step 113, determining the original pollutant concentration data of the first grid according to the preset relation between the first grid and all the sub-grids and the pollutant concentration data of each sub-grid.
Step 113 specifically includes:
determining the number of the sub-grids according to the preset relation between the first grid and all the sub-grid information;
determining a sum of the contaminant concentration data from the contaminant concentration data for each sub-grid;
the raw contaminant concentration data for the first grid is determined based on the sum of the number of sub-grids and the contaminant concentration data. Thus, the average value of the pollutant concentration data of the small grids is used as the original pollutant concentration data of the first grid, and the accuracy of the data is improved.
Step 114, obtaining location information of the first grid.
The position information of the first grid comprises upper left longitude and latitude data, lower left longitude and latitude data, upper right longitude and latitude data, lower right longitude and latitude data and middle longitude and latitude data of the first grid.
Step 115, calculating a fourth distance between the first grid and the first device based on the first point location information and the location information of the first grid.
And step 116, correcting the original pollutant concentration data of the first grid according to the fourth distance and the first pollutant concentration data to obtain fourth correction data.
Step 117, calculating a fifth distance between the first grid and the second device according to the second point location information and the position information of the first grid.
Step 118, correcting the fourth correction data according to the fifth distance and the second pollutant concentration data to obtain fifth correction data.
Step 119, calculating a sixth distance between the first grid and the third device according to the third point location information and the position information of the first grid.
And step 120, correcting the fifth correction data according to the sixth distance and the third pollutant concentration data to obtain sixth correction data, and taking the sixth correction data as the pollutant concentration data of the first grid.
It will be appreciated that, when the second to sixth correction data are obtained later, an inverse distance weight interpolation method may be also used, and the specific steps are similar to those for obtaining the first correction data, and will not be repeated here.
By applying the method for determining the pollutant concentration, provided by the invention, the interpolation method is utilized to calculate the original pollutant concentration data of the small grids, then the original pollutant concentration data of each sub-grid is sequentially controlled with the first pollutant concentration data of the first equipment, the second pollutant concentration data of the second equipment and the third pollutant concentration data of the third equipment to obtain the pollutant concentration data of the small grids, and then the pollutant concentration data of each small grid is sequentially controlled with the first pollutant concentration data of the first equipment, the second pollutant concentration data of the second equipment and the third pollutant concentration data of the third equipment to obtain the pollutant concentration data of the first grid, so that the area is divided into grids, the pollutant concentration data of the sub-grids in the grids are obtained, a novel method for obtaining the pollutant concentration data is provided, the accuracy of the pollutant concentration data is improved, and the accuracy and the effectiveness of the whole environment monitoring work are improved.
Fig. 2 is a schematic diagram of a determining device for determining a concentration of a contaminant according to a second embodiment of the present invention. The pollutant concentration determining device can be applied to a pollutant concentration determining method. As shown in fig. 2, the contaminant concentration determination device 200 includes: an acquisition unit 210, a calculation unit 220, a correction unit 230, and a determination unit 240.
The acquiring unit 210 is configured to acquire first point location information and first contaminant concentration data of at least three first devices within the first grid.
The obtaining unit 210 is further configured to obtain second point location information and second pollutant concentration data of the second device.
The obtaining unit 210 is further configured to obtain third point location information and third pollutant concentration data of a third device.
The obtaining unit 210 is further configured to obtain all pieces of sub-grid information of the first grid, where each piece of sub-grid information in the all pieces of sub-grid information includes location information of the sub-grid.
The acquisition unit is further configured to acquire first data of the first region, wherein the first data includes positional information of the first region and contaminant concentration data of the first region.
The calculation unit 220 is configured to calculate raw contaminant concentration data of each sub-grid based on the position information of the first region, the contaminant concentration data of the first region, and the position information of each sub-grid.
The calculating unit 220 is further configured to calculate a first distance between each sub-grid and the first device according to the first point location information and the position information of each sub-grid.
The calculating unit 220 is further configured to calculate a second distance between each sub-grid and the second device according to the second point location information and the position information of each sub-grid.
The correction unit 230 is further configured to correct the first correction data of each sub-grid according to the second distance and the second contaminant concentration data, so as to obtain second correction data.
The calculating unit 220 is further configured to calculate a third distance between each sub-grid and the third device according to the third point location information and the position information of each sub-grid.
The correction unit 230 is further configured to correct the second correction data of each sub-grid according to the third distance and the third pollutant concentration data, to obtain third correction data, and to use the third correction data of each sub-grid as the pollutant concentration data thereof.
The determining unit 240 is configured to determine the original contaminant concentration data of the first grid according to the preset relationship between the first grid and the information of all the sub-grids and the contaminant concentration data of each sub-grid.
The obtaining unit 210 is further configured to obtain location information of the first grid.
The calculating unit 220 is further configured to calculate a fourth distance between the first grid and the first device based on the first point location information and the location information of the first grid.
The correction unit 230 is further configured to correct the original contaminant concentration data of the first grid according to the fourth distance and the first contaminant concentration data, so as to obtain fourth corrected data.
The calculating unit 220 is further configured to calculate a fifth distance between the first mesh and the second device according to the second point location information and the position information of the first mesh.
The correction unit 230 is further configured to correct the fourth correction data according to the fifth distance and the second contaminant concentration data, so as to obtain fifth correction data.
The calculating unit 220 is further configured to calculate a sixth distance between the first mesh and the third device according to the third point location information and the position information of the first mesh.
The correction unit 230 is further configured to correct the fifth correction data according to the sixth distance and the third pollutant concentration data, obtain sixth correction data, and use the sixth correction data as the pollutant concentration data of the first grid.
Further, the calculating unit 220 is specifically configured to:
acquiring longitude and latitude data of a first point bit;
acquiring longitude and latitude data of each sub-grid, wherein the longitude and latitude data of each sub-grid comprises upper left longitude and latitude data, lower left longitude and latitude data, upper right longitude and latitude data, lower right longitude and latitude data and middle longitude and latitude data;
And calculating the first distance between each sub-grid and the first device according to the longitude and latitude data of the first point location, the upper left longitude and latitude data, the lower left longitude and latitude data, the upper right longitude and latitude data, the lower right longitude and latitude data and the middle longitude and latitude data.
Further, the calculating unit 220 is specifically configured to:
determining an overlapping region of the first region and each sub-grid according to the position information of the first region and the position information of each sub-grid;
determining contaminant concentration data for the sub-grid of the overlapping region based on the contaminant concentration data for the first region;
and obtaining the pollutant concentration data of the subgrid of the non-overlapping area through interpolation, wherein the pollutant concentration data of the subgrid of the overlapping area and the pollutant concentration data of the subgrid of the non-overlapping area form the original pollutant concentration data of all subgrids. Further, the correction unit 230 is specifically configured to:
determining a correlation coefficient according to a preset radius constant and position information of the sub-grid;
and obtaining first correction data according to the original pollutant concentration data of each sub-grid, the correlation coefficient and the first distance.
Further, the determining unit 240 is specifically configured to:
determining the number of the sub-grids according to the preset relation between the first grid and all the sub-grid information;
Determining a sum of the contaminant concentration data from the contaminant concentration data for each sub-grid;
the raw contaminant concentration data for the first grid is determined based on the sum of the number of sub-grids and the contaminant concentration data.
An embodiment of the present invention provides an apparatus, including a memory and a processor, where the memory is configured to store a program, and the memory may be connected to the processor through a bus. The memory may be non-volatile memory, such as a hard disk drive and flash memory, in which software programs and device drivers are stored. The software program can execute various functions of the above method provided in the first embodiment of the present invention; the device driver may be a network and interface driver. The processor is configured to execute a software program, where the software program is executed to implement the method provided in the first embodiment of the present invention.
A fourth embodiment of the present invention provides a computer program product containing instructions, which when executed on a computer, cause the computer to perform the method provided by the first embodiment of the present invention.
The fifth embodiment of the present invention provides a computer readable storage medium, where a computer program is stored, and when the computer program is executed by a processor, the method provided in the first embodiment of the present invention is implemented.
Those of skill would further appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the various illustrative elements and steps are described above generally in terms of function in order to clearly illustrate the interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
The steps of a method or algorithm described in connection with the embodiments disclosed herein may be embodied in hardware, in a software module executed by a processor, or in a combination of the two. The software modules may be disposed in Random Access Memory (RAM), memory, read Only Memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.
The foregoing detailed description of the invention has been presented for purposes of illustration and description, and it should be understood that the invention is not limited to the particular embodiments disclosed, but is intended to cover all modifications, equivalents, alternatives, and improvements within the spirit and principles of the invention.

Claims (10)

1. A method of determining a concentration of a contaminant, the method comprising:
acquiring first point position information and first pollutant concentration data of at least three first devices in a first grid, wherein the first devices are quality control devices;
acquiring second point location information and second pollutant concentration data of second equipment, wherein the second equipment is provincial control equipment;
acquiring third point location information and third pollutant concentration data of third equipment, wherein the third equipment is state control equipment;
acquiring all pieces of sub-grid information of a first grid, wherein each piece of sub-grid information in all pieces of sub-grid information comprises the position information of the sub-grid;
acquiring first data of a first area, wherein the first data comprises position information of the first area and pollutant concentration data of the first area;
calculating original pollutant concentration data of each sub-grid according to the position information of the first area, the pollutant concentration data of the first area and the position information of each sub-grid;
according to the first point position information and the position information of each sub-grid, respectively calculating a first distance between each sub-grid and the first device;
Correcting the original pollutant concentration data of each sub-grid according to the first distance and the original pollutant concentration data of each sub-grid to obtain first correction data;
respectively calculating a second distance between each sub-grid and the second device according to the second point location information and the position information of each sub-grid;
correcting the first correction data of each sub-grid according to the second distance and the second pollutant concentration data to obtain second correction data;
according to the third point location information and the position information of each sub-grid, calculating a third distance between each sub-grid and the third device;
correcting the second correction data of each sub-grid according to the third distance and the third pollutant concentration data to obtain third correction data, and taking the third correction data of each sub-grid as pollutant concentration data;
determining original pollutant concentration data of the first grid according to the preset relation between the first grid and all the sub-grid information and the pollutant concentration data of each sub-grid;
acquiring position information of the first grid;
calculating a fourth distance between the first grid and the first device according to the first point position information and the position information of the first grid;
Correcting the original pollutant concentration data of the first grid according to the fourth distance and the first pollutant concentration data to obtain fourth correction data;
calculating a fifth distance between the first grid and the second device according to the second point location information and the position information of the first grid;
correcting the fourth correction data according to the fifth distance and the second pollutant concentration data to obtain fifth correction data; calculating a sixth distance between the first grid and the third device according to the third point location information and the position information of the first grid;
and correcting the fifth correction data according to the sixth distance and the third pollutant concentration data to obtain sixth correction data, and taking the sixth correction data as the pollutant concentration data of the first grid.
2. The method for determining a concentration of a contaminant according to claim 1, wherein said calculating a first distance between each of said sub-grids and said first device based on said first point location information and location information of each of said sub-grids, respectively, comprises: acquiring longitude and latitude data of a first point bit;
acquiring longitude and latitude data of each sub-grid, wherein the longitude and latitude data of each sub-grid comprises upper left longitude and latitude data, lower left longitude and latitude data, upper right longitude and latitude data, lower right longitude and latitude data and middle longitude and latitude data;
And calculating a first distance between each sub-grid and the first device according to the longitude and latitude data of the first point location, the upper left longitude and latitude data, the lower left longitude and latitude data, the upper right longitude and latitude data, the lower right longitude and latitude data and the middle longitude and latitude data.
3. The method for determining a concentration of a contaminant according to claim 1, wherein said calculating the original concentration data of each sub-grid based on the position information of the first area, the concentration data of the first area, and the position information of each sub-grid specifically comprises:
determining an overlapping area of the first area and the sub-grids according to the position information of the first area and the position information of each sub-grid;
determining contaminant concentration data of the sub-grid of the overlapping region according to the contaminant concentration data of the first region;
and obtaining pollutant concentration data of the subgrid of the non-overlapping area through interpolation, wherein the pollutant concentration data of the subgrid of the overlapping area and the pollutant concentration data of the subgrid of the non-overlapping area form original pollutant concentration data of all subgrids.
4. The method for determining a concentration of a contaminant according to claim 1, wherein said correcting the original contaminant concentration data of each sub-grid based on said first distance and said original contaminant concentration data of each sub-grid to obtain first corrected data comprises:
Determining a correlation coefficient according to a preset radius constant and the position information of the sub-grid;
and obtaining first correction data according to the original pollutant concentration data of each sub-grid, the correlation coefficient and the first distance.
5. The method for determining a concentration of a contaminant according to claim 1, wherein said determining the original concentration data of the first grid according to the predetermined relation between the first grid and the information of all the sub-grids and the concentration data of each sub-grid specifically comprises:
determining the number of the sub-grids according to the preset relation between the first grid and the whole sub-grid information;
determining a sum of the contaminant concentration data from the contaminant concentration data for each sub-grid;
and determining the original pollutant concentration data of the first grid according to the sum of the number of the sub-grids and the pollutant concentration data.
6. A pollutant concentration determining apparatus, characterized in that the pollutant concentration determining apparatus comprises:
an acquisition unit for acquiring first point location information and first contaminant concentration data of at least three first devices within a first grid;
The acquisition unit is further used for acquiring second point location information and second pollutant concentration data of second equipment;
the acquisition unit is further used for acquiring third point location information and third pollutant concentration data of third equipment;
the acquiring unit is further configured to acquire all pieces of sub-grid information of the first grid, where each piece of sub-grid information in the all pieces of sub-grid information includes location information of the sub-grid;
the acquisition unit is further used for acquiring first data of the first area, wherein the first data comprises position information of the first area and pollutant concentration data of the first area;
a calculating unit, configured to calculate original contaminant concentration data of each sub-grid according to the position information of the first area, the contaminant concentration data of the first area, and the position information of each sub-grid;
the computing unit is further configured to respectively compute a first distance between each sub-grid and the first device according to the first point position information and the position information of each sub-grid;
the correction unit is used for correcting the original pollutant concentration data of each sub-grid according to the first distance and the original pollutant concentration data of each sub-grid to obtain first correction data;
The calculating unit is further used for calculating a second distance between each sub-grid and the second device according to the second point location information and the position information of each sub-grid;
the correction unit is further used for correcting the first correction data of each sub-grid according to the second distance and the second pollutant concentration data to obtain second correction data;
the calculating unit is further configured to calculate a third distance between each sub-grid and the third device according to the third point location information and the position information of each sub-grid;
the correction unit is further configured to correct the second correction data of each sub-grid according to the third distance and the third pollutant concentration data, to obtain third correction data, and use the third correction data of each sub-grid as pollutant concentration data thereof;
the determining unit is used for determining original pollutant concentration data of the first grid according to the preset relation between the first grid and all the sub-grid information and the pollutant concentration data of each sub-grid;
the acquisition unit is further used for acquiring the position information of the first grid;
The computing unit is further configured to calculate a fourth distance between the first grid and the first device according to the first point position information and the position information of the first grid;
the correction unit is further used for correcting the original pollutant concentration data of the first grid according to the fourth distance and the first pollutant concentration data to obtain fourth correction data;
the calculating unit is further used for calculating a fifth distance between the first grid and the second device according to the second point location information and the position information of the first grid;
the correction unit is further configured to correct the fourth correction data according to the fifth distance and the second pollutant concentration data, so as to obtain fifth correction data;
the calculating unit is further configured to calculate a sixth distance between the first grid and the third device according to the third point location information and the position information of the first grid;
the correction unit is further configured to correct the fifth correction data according to the sixth distance and the third pollutant concentration data, obtain sixth correction data, and use the sixth correction data as the pollutant concentration data of the first grid.
7. The pollutant concentration determination device of claim 6, wherein the computing unit is specifically configured to: acquiring longitude and latitude data of a first point bit;
acquiring longitude and latitude data of each sub-grid, wherein the longitude and latitude data of each sub-grid comprises upper left longitude and latitude data, lower left longitude and latitude data, upper right longitude and latitude data, lower right longitude and latitude data and middle longitude and latitude data;
and calculating a first distance between each sub-grid and the first device according to the longitude and latitude data of the first point location, the upper left longitude and latitude data, the lower left longitude and latitude data, the upper right longitude and latitude data, the lower right longitude and latitude data and the middle longitude and latitude data.
8. The pollutant concentration determination device of claim 6, wherein the computing unit is specifically configured to: determining an overlapping area of the first area and the sub-grids according to the position information of the first area and the position information of each sub-grid;
determining contaminant concentration data of the sub-grid of the overlapping region according to the contaminant concentration data of the first region;
and obtaining pollutant concentration data of the subgrid of the non-overlapping area through interpolation, wherein the pollutant concentration data of the subgrid of the overlapping area and the pollutant concentration data of the subgrid of the non-overlapping area form original pollutant concentration data of all subgrids.
9. The pollutant concentration determination device of claim 6, wherein the correction unit is specifically configured to: determining a correlation coefficient according to a preset radius constant and the position information of the sub-grid;
and obtaining first correction data according to the original pollutant concentration data of each sub-grid, the correlation coefficient and the first distance.
10. The pollutant concentration determination device of claim 6, wherein the determination unit is specifically configured to: determining the number of the sub-grids according to the preset relation between the first grid and the whole sub-grid information;
determining a sum of the contaminant concentration data from the contaminant concentration data for each sub-grid;
and determining the original pollutant concentration data of the first grid according to the sum of the number of the sub-grids and the pollutant concentration data.
CN201811083073.6A 2018-09-17 2018-09-17 Method and device for determining concentration of pollutant Active CN109085100B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811083073.6A CN109085100B (en) 2018-09-17 2018-09-17 Method and device for determining concentration of pollutant

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811083073.6A CN109085100B (en) 2018-09-17 2018-09-17 Method and device for determining concentration of pollutant

Publications (2)

Publication Number Publication Date
CN109085100A CN109085100A (en) 2018-12-25
CN109085100B true CN109085100B (en) 2023-10-31

Family

ID=64841798

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811083073.6A Active CN109085100B (en) 2018-09-17 2018-09-17 Method and device for determining concentration of pollutant

Country Status (1)

Country Link
CN (1) CN109085100B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109765149B (en) * 2018-12-27 2021-12-03 北京英视睿达科技有限公司 Method and device for determining dust reduction amount based on total suspended particulate matters
CN109814549B (en) * 2018-12-29 2023-11-10 北京英视睿达科技股份有限公司 Method, device and system for determining pollution source position
CN116738539B (en) * 2023-06-09 2024-01-09 西安交通大学 Pollutant diffusion prediction system and method based on hydrodynamic model

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103336093A (en) * 2013-06-26 2013-10-02 中山大学 Regional spatial quality analysis method
RU2503042C1 (en) * 2012-08-21 2013-12-27 Федеральное бюджетное учреждение науки "Федеральный научный центр медико-профилактических технологий управления рисками здоровью населения" (ФБУН "ФНЦ медико-профилактических технологий управления рисками здоровью населения") Method for spatial quantitative evaluation of ambient air contamination level
KR20150031577A (en) * 2013-09-16 2015-03-25 주식회사 에니텍 Method of Inverse Calculation for The Amount of Air Pollutants Emissions
CN106407714A (en) * 2016-10-14 2017-02-15 珠海富鸿科技有限公司 Air pollution assessment method and device based on CALPUFF system
CN107607450A (en) * 2017-09-15 2018-01-19 深圳市卡普瑞环境科技有限公司 A kind of air quality surveillance method and relevant device

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2503042C1 (en) * 2012-08-21 2013-12-27 Федеральное бюджетное учреждение науки "Федеральный научный центр медико-профилактических технологий управления рисками здоровью населения" (ФБУН "ФНЦ медико-профилактических технологий управления рисками здоровью населения") Method for spatial quantitative evaluation of ambient air contamination level
CN103336093A (en) * 2013-06-26 2013-10-02 中山大学 Regional spatial quality analysis method
KR20150031577A (en) * 2013-09-16 2015-03-25 주식회사 에니텍 Method of Inverse Calculation for The Amount of Air Pollutants Emissions
CN106407714A (en) * 2016-10-14 2017-02-15 珠海富鸿科技有限公司 Air pollution assessment method and device based on CALPUFF system
CN107607450A (en) * 2017-09-15 2018-01-19 深圳市卡普瑞环境科技有限公司 A kind of air quality surveillance method and relevant device

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
ArcGIS在大气污染源排放清单建立中的应用研究;张晓郁;陈振飞;徐文哲;;环境科学与管理(第01期);第135-139页 *

Also Published As

Publication number Publication date
CN109085100A (en) 2018-12-25

Similar Documents

Publication Publication Date Title
Ravindra et al. Air pollution in India: bridging the gap between science and policy
Super et al. Uncertainty analysis of a European high-resolution emission inventory of CO 2 and CO to support inverse modelling and network design
CN109085100B (en) Method and device for determining concentration of pollutant
Elbir A GIS based decision support system for estimation, visualization and analysis of air pollution for large Turkish cities
CN111125641B (en) Atmospheric pollution anomaly identification and cause analysis method, device and storage medium
Pesch et al. Mosses as bioindicators for metal accumulation: Statistical aggregation of measurement data to exposure indices
CN112991132B (en) PM2.5Method and device for identifying area cooperatively controlled by ozone
CN109524130A (en) Determination method and device of the vehicle to pollutant concentration contribution rate
Smith et al. Use of GIS and ancillary variables to predict volatile organic compound and nitrogen dioxide levels at unmonitored locations
CN115436573B (en) Intelligent monitoring method and device for atmospheric pollution source
CN114218976B (en) Method and system for tracing and positioning refined emission source in high-value area of ozone precursor
Kumar et al. Air quality interventions and spatial dynamics of air pollution in Delhi and its surroundings
ApSimon et al. The UK Integrated Assessment Model for source apportionment and air pollution policy applications to PM2. 5
Carnevale et al. Applying the delta tool to support the air quality directive: evaluation of the TCAM chemical transport model
Puyana-Romero et al. A 3D GIS tool for the detection of noise hot-spots from major roads
CN109613179B (en) Method for determining cumulative high value area
CN109856321A (en) The determination method of abnormal high level point
Dios et al. A mixed top-down and bottom-up methodology in spatial segregation of emissions based on GIS tools
Halimi et al. Modeling spatial distribution of Tehran air pollutants using geostatistical methods incorporate uncertainty maps
Sartini et al. From emissions to source allocation: Synergies and trade-offs between top-down and bottom-up information
CN109507367A (en) Determine the method and device of atmosphere pollution fining distribution
CN111461163B (en) Urban interior PM2.5 concentration simulation and population exposure evaluation method and device
CN109521155B (en) Quality control method and device
Zhu et al. A data-driven approach for optimal design of integrated air quality monitoring network in a chemical cluster
CN115705510A (en) Factory gas pollution distribution prediction method and system, electronic equipment and storage medium

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
CB02 Change of applicant information

Address after: 100070 Beijing Fengtai Auto Museum West Road No. 8 Courtyard 1 Building 6 Floor 606

Applicant after: Beijing Yingshi Ruida Technology Co.,Ltd.

Address before: 1202, No. 18, Yumin Middle Road, Xicheng District, Beijing 100029

Applicant before: BEIJING INSIGHTS VALUE TECHNOLOGY Co.,Ltd.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant