KR102492322B1 - System For Classifying Types of Effects of Through Environmental Satellite And Its Method - Google Patents

Abstract

The present invention relates to a system and a method for classifying the types of aerosol effects through environmental satellites, which determine whether a high-concentration aerosol is generated, and when a high-concentration aerosol is identified, determine whether the aerosol originates from abroad or domestically, and plan and establish countermeasures against domestic effects of the aerosol. The classification system of the present invention includes a map generation module, an aerosol detection module, a region setting module, a region division module, a variable calculation module, an initial classification module, a reclassification module, and a type setting module.

Description

System For Classifying Types of Effects of Through Environmental Satellite And Its Method

The present invention relates to an aerosol detection and classification system. More specifically, it is a technology related to a system that determines and classifies high concentrations of aerosols detected using environmental satellites and whether or not they are affected at home or abroad.

As the industry developed, the environment was polluted and destroyed. Many people are making great efforts to purify the polluted environment in hopes of living in a clean environment. For example, eco-friendly policies are being announced and developed in various fields, such as the prohibition of fossil fuel-based power generation and the development of new and renewable energy, mainly in the United States and Europe. However, fine dust is still generated from coal power generation in India and China, and exhaust gas is generated from the operation of ships transporting cargo between continents.

Currently, various researches and developments are being conducted to purify air pollution caused by fine dust and exhaust gas. In that study, studies on the causes of environmental pollution and the detection of polluted air are being actively conducted. For example, research and development on a technology that can determine the aerosol concentration in a specific area through geostationary orbit and a technology that can check the movement of pollutants in real time over a wide area in the case of a geostationary satellite is being conducted.

However, in addition to the above-mentioned research and development, research and development on causes that can identify and solve the causes of air pollution, that is, detection of aerosols, inflow of aerosols to a specific area, and movement of the inflowed aerosols, are insufficient.

Republic of Korea Patent Registration No. 10-2201084 (Announcement date: 2021.01.11)

The present invention is intended to solve the problem of not being able to properly respond to the inflow of polluted air because the inflow and outflow paths of the polluted atmosphere and the concentration of the polluted atmosphere cannot be grasped.

The problem to be solved by the present invention is not limited to the problems mentioned above, and other technical problems not mentioned will be understood by those skilled in the art from the following description.

In order to achieve the object to be solved, the system for classifying the impact type of aerosol through environmental satellites of the present invention includes a communication module for receiving location data and aerosol optical depth (AOD) data from the satellite, and a communication module. A map generation module that generates a grid map in which a grid area is generated with the received location data, and an aerosol detection module that detects aerosols from the aerosol optical thickness data received by the communication module and calculates the number of aerosol pixels by calculating the detected aerosols according to pixels. module, a region setting module for generating at least one set region in the grid map set in the map generation module, a region division module that divides the set region into a domestic region and a plurality of gateway regions, and the information received from the aerosol detection module during the set period A variable value calculation module that calculates the median value of the aerosol optical thickness ( m ), the standard deviation of the aerosol optical thickness ( σ ), and the average value of the aerosol optical thickness ( μ ) within the set area through the number of aerosol pixels, and during the set period Comparing the number of aerosol pixels in the setting area (E) with the standard number, if the number of aerosol pixels in the setting area (E) exceeds the reference number during the setting period, the first classification process is determined, and the aerosol pixels in the setting area (E) during the setting period If the number is less than the reference number, the initial classification module determines the second classification process, and if the initial classification module determines the first classification process, the average value ( μ2 ) is calculated, and the median value ( m2 ) of the aerosol optical thickness within the plurality of gateway areas calculated during the set period through the variable value calculation module and the standard deviation value (σ2) of the aerosol optical thickness are calculated by the second formula ( m2+3σ2 ) ) to calculate the gateway area boundary value, and compare the calculated average value ( μ 2 ) of the aerosol optical thickness in the plurality of gateway areas and the gateway area boundary value, and the average value ( μ 2 ) of the aerosol optical thickness in the plurality of gateway areas If the gateway area boundary value is exceeded, the 1-1 classification process is determined, and if the average value ( μ 2 ) of the aerosol optical thickness in a plurality of gateway areas is less than the gateway area boundary value, the first classification process is determined. If the reclassification module that determines the -2 classification process and the reclassification module determines the 1-1 classification process, the aerosol's moving direction is reversely tracked, and when the aerosol passes through the gateway area within the preset time, the aerosol is transported over a long distance. When it is determined as the movement (LRT) type and the aerosol does not pass through the gateway area within the preset time, the aerosol is determined as the first non-LRT type and the reclassification module determines the first-second classification process. Then, the average value ( μ3 ) of the aerosol optical thickness within the domestic area during the setting period is calculated through the variable value calculation module, and the median value ( m3 ) of the aerosol optical thickness within the domestic area calculated during the setting period through the variable value calculation module and the aerosol optical thickness The standard deviation of the thickness ( σ3 ) is substituted into the third formula ( m3 + 3σ3 ) to calculate the domestic boundary value, and the calculated average value of the aerosol optical thickness ( μ3 ) in the domestic region is compared with the domestic boundary value Therefore, if the average value ( μ3 ) of the aerosol optical thickness in the domestic area exceeds the domestic area boundary value, the aerosol is determined as the second non-long-distance travel (Non-LRT) type, and the average value ( μ3 ) of the aerosol optical thickness in the domestic area is If it is below the boundary value of the area, the aerosol is determined as a clean type that does not affect the domestic area, and if the initial classification module determines the second classification process, the aerosol is determined as a cloud in the airspace of the domestic area. Contains the type setting module.

Here, the re-classification module checks the number of aerosol pixels in the plurality of gateway areas calculated by the variable value calculation module when the first classification process is determined, and determines the first-third classification process when the number of aerosol pixels identified is zero. In addition, the type setting module determines the number of aerosols to be unknown when the reclassification module determines the first to third classification processes.

In addition, the aerosol influence type classification system through the environmental satellite calculates the median value ( m ) of the aerosol optical thickness and the standard deviation value ( σ ) of the aerosol optical thickness calculated in the variable value calculation module by the first formula ( m - 3σ ) ) to calculate a first formula value, and may further include a background density setting module for setting an average value of the first formula value calculated during the setting period as the background density of the grid area.

In addition, in the system for classifying the impact type of aerosols through environmental satellites of the present invention, when the vector aerosol generation module and the initial classification module that generate wind vector aerosols to which wind power is applied to the aerosols detected by the aerosol detection module determine the second classification process, The median value ( mX ) of the aerosol optical thickness of the wind vector in the grid area calculated by the variable value calculation module, the standard deviation value of the aerosol optical thickness ( σX ), and the average value of the aerosol optical thickness ( μX ) are calculated by the second or third equation By substituting, a value larger than the value calculated by substituting, that is, the difference between the excess value of the second formula and the background concentration value set in the grid area through the background concentration setting module, is substituted into the preset vertical variance equation in the map generation module It may further include a movement amount calculation module that calculates the inflow and outflow per generated grid area.

In addition, it is connected to the map generation module, area setting module, area division module, variable value calculation module, background concentration setting module, initial classification module, reclassification module, type setting module, displacement calculation module, and vector aerosol generation module, and is output from each module. It may further include an expression module that expresses the result to the screen.

Another method for classifying the impact type of aerosols through environmental satellites of the present invention to achieve the above object is the step (a) of the communication module receiving location data and aerosol optical thickness data from the satellite, map generation module (b) step (S20) of generating a grid map in which a grid area is generated with the location data received from the communication module, the aerosol detection module detects aerosol from the aerosol optical thickness data received from the communication module, and converts the detected aerosol into pixels Step (c) of calculating the number of aerosol pixels by calculating according to (S30), step (d) of generating an area of a certain size as a set area in the grid map set by the area setting module in the map generation module (S40), area division (e) step (S50) of dividing the setting area into a domestic area and a plurality of gateway areas (F) by the module, and calculating the aerosol optical thickness within the setting area (E) from the aerosol optical thickness data received during the setting period by the variable value calculation module. (f) step (S60) of calculating the median value ( m ), the standard deviation value ( σ ) of the aerosol optical thickness, and the average value ( μ ) of the aerosol optical thickness, respectively, the number of aerosol pixels in the setting area during the setting period by the initial classification module In comparison with the reference number, determining the first classification process when the number of aerosol pixels in the setting area exceeds the reference number during the setting period, and determining the second classification process when the number of aerosol pixels in the setting area is less than the reference number during the setting period ( g) In step (S70), when the reclassification module determines the first classification process in the initial classification module, the average value ( μ2 ) of the aerosol optical thickness in the plurality of gateway areas is calculated through the variable value calculation module during the set period, and the variable The median value ( m2 ) of the aerosol optical thickness within the plurality of gateway areas (F) and the standard deviation (σ2) of the aerosol optical thickness calculated during the set period through the value calculation module are calculated by the second formula ( m2+3σ2 ) The gateway area boundary value is calculated by substitution, and the average value of the aerosol optical thickness ( μ2 ) in the plurality of gateway areas is compared with the gate area boundary value, and the average value ( μ2 ) of the aerosol optical thickness in the plurality of gateway areas is If the boundary value is exceeded, the 1-1 classification (h) step (S80) of determining the first and second classification process when the average value ( μ2 ) of the aerosol optical thickness in the plurality of gateway areas is less than the gateway area boundary value and the type setting module in the reclassification module When the 1-1 classification process is determined, the aerosol's movement direction is reversely tracked, and when the aerosol passes through the gateway area within a preset time, the aerosol is determined as a long-distance movement (LRT) type, and the aerosol passes through the barrier within a preset time. When it does not pass through the area, the aerosol is determined as the first non-LRT type, and if the first-second classification process is determined in the reclassification module, the variable value calculation module determines the aerosol within the domestic area during the set period. The average value ( μ3 ) of the aerosol optical thickness is calculated, and the median ( m3 ) and aerosol optical thickness of the aerosol optical thickness within the domestic area calculated during the set period through the variable value calculation module are calculated in the third formula ( m3 + 3σ3 ). Substitute the standard deviation value ( σ3 ) of to calculate the domestic area boundary value, and compare the calculated average value ( μ3 ) of the domestic area aerosol optical thickness and the domestic area boundary value, and compare the average value ( μ3 ) of the aerosol optical thickness of the domestic area. ) exceeds the domestic boundary value, the aerosol is determined as the second non-long-distance travel (Non-LRT) type, and if the average value ( μ3 ) of the aerosol optical thickness in the domestic region (G) is less than the domestic boundary value, the aerosol is If it is determined as a clean type that does not affect the domestic area and the second classification process is determined in the initial classification module, (i) step (S90) of determining the type of aerosol as a cloud in the airspace of the domestic area is included. do.

Here, in step (h), the reclassification module checks the number of aerosol pixels in the plurality of gateway areas calculated by the variable value calculation module, and when the checked number is 0, further steps of determining the 1-3 classification process. Step (i) may further include determining that the aerosol is in an unknown state when the step of determining the first-third classification process in step (h) by the type setting module proceeds.

The present invention determines whether a high concentration of aerosol is generated, and when it is determined as a high concentration aerosol, it is possible to determine whether the aerosol is from overseas or from Korea, so that countermeasures against aerosol affecting Korea can be planned and established.

In other words, the present invention enables those involved in improving air pollution to constantly monitor the concentration of air quality and to plan correct policies for air pollutant generation.

1 is a diagram showing an operating state of a system for classifying impact types of aerosols through environmental satellites according to the present invention.
2 is a diagram showing location data and aerosol optical thickness data generated by a satellite.
3 is a block diagram of components constituting a system for classifying the impact types of aerosols through environmental satellites according to an embodiment of the present invention.
FIG. 4 is a diagram showing an operating state of the map generating module of FIG. 3 .
5 is a view showing an operating state of the aerosol detection module of FIG. 3;
FIG. 6 is a diagram showing an operating state of the area setting module and the area division module of FIG. 3 .
7 is a diagram showing an operating state of the variable value calculation module of FIG. 3 .
FIG. 8 is a diagram showing an operating state of the background concentration setting module of FIG. 3 .
9 is a diagram showing an operating state of the initial classification module of FIG. 3 .
10 and 11 are views showing an operating state of the reclassification module of FIG. 3 .
12 is a diagram showing an operating state of the type setting module of FIG. 3 .
13 and 14 are views showing an operating state of the expression module of FIG. 3 .
15 is a flow chart of a method for classifying impact types of aerosols through environmental satellites according to an embodiment of the present invention.
16 is a block diagram of components constituting a system for classifying impact types of aerosols through environmental satellites according to another embodiment of the present invention.
17 and 18 are views showing an operating state of the vector aerosol generating module of FIG. 16 .
FIG. 19 is a diagram showing an operating state of the movement amount calculation module of FIG. 3 .

Advantages and features of the present invention and a system for achieving them may be shown through the embodiments described below in detail in conjunction with the accompanying drawings. However, the present invention is not limited to the embodiments disclosed below and may be implemented in a variety of different forms.

However, the embodiments disclosed in this specification are one example that allows those skilled in the art to fully understand the scope of the invention while ensuring the disclosure of the present invention is complete. Accordingly, the same reference numerals described throughout the specification denote the same components so that the disclosure of the present invention is complete and those skilled in the art can fully understand the present invention through the posted examples.

Even if the description of one embodiment of the present invention is specifically described in the specific content for carrying out the invention, the claims of the present invention are not limited to one embodiment. The claims of the present invention are defined solely by the claims.

Hereinafter, with reference to FIGS. 1 to 18, a system for classifying impact types of aerosols through environmental satellites and a method for classifying impact types of aerosols through environmental satellites will be described in detail. However, in order to concisely understand the description of the aerosol impact type classification system through the environmental satellite and the aerosol impact type classification method through the environmental satellite of the present invention, with reference to FIGS. 1 to 3, the system of the present invention and the present invention The overall content covering the method of the invention will be described.

Here, FIG. 1 is a diagram showing the operating state of the aerosol impact type classification system through environmental satellites of the present invention, and FIG. 2 is a diagram showing location data and aerosol optical thickness data generated by the satellite. 3 is a block diagram of components constituting a system for classifying the impact type of aerosol through environmental satellites according to an embodiment of the present invention.

After explaining FIGS. 1 to 3 , components constituting the system of the present invention will be described in detail with reference to FIGS. 4 to 14 . In addition, based on the description of the aerosol impact type classification system through the above-described environmental satellite and the method of classifying the aerosol impact type through the environmental satellite of the present invention based on FIG. 15, a detailed description will be given.

The system (1) for classifying the impact type of aerosols through environmental satellites of the present invention and the method receive data including meteorological information generated while the satellite (A) photographs the earth (B), and determines whether aerosols exist around the Korean Peninsula from the data. to be able to detect In addition, it is possible to grasp the inflow path and current status of the detected aerosol. Here, the satellite A may generate and transmit location data and aerosol optical depth (AOD) data. The aerosol influence type classification system (1) and its method through environmental satellites observe a wide area every hour, detect whether high-concentration aerosols have foreign influences, and quantify the detected aerosols.

Therefore, the classification system (1) of impact types of aerosols through environmental satellites and its method enable those involved in improving air pollution to monitor the concentration of air quality at all times, and to establish correct policies for the generation of air pollutants. allow this to be planned.

The aerosol impact type classification system 1 through the environmental satellites includes a communication module 10 capable of wireless communication with the satellite A, processes data received from the satellite A, and outputs the processed information. It can be a computer with a program installed. Aerosol impact type classification system (1) through environmental satellites, that is, the computer includes a communication module (10), a map generation module (20), an aerosol detection module (30), a region setting module (40), and a region division module (50) , a variable value calculation module 60, a background density setting module 65, an initial classification module 70, a reclassification module 80, and a type setting module 90 may be included as components. And it may further include the expression module 100 as a component. Here, the communication module 10 is a module that wirelessly communicates with the satellite A, that is, the geostationary satellite, and receives location data and aerosol optical thickness data, that is, environmental data. In addition, the communication module 10 communicates with the Korea Meteorological Administration and receives wind data from the Korea Meteorological Administration.

Hereinafter, with reference to FIGS. 4 to 14, components constituting the system for classifying the impact type of aerosol through environmental satellites according to the present invention will be described in detail.

4 is a view showing an operating state of the map generation module of FIG. 3, FIG. 5 is a diagram showing an operating state of the aerosol detection module of FIG. 3, and FIG. 6 is an operating state of the region setting module and region segmentation module of FIG. 3 is a drawing showing 7 is a diagram showing the operating state of the variable value calculation module of FIG. 3, FIG. 8 is a diagram showing the operating state of the background concentration setting module of FIG. 3, and FIG. 9 is a diagram showing the operating state of the initial classification module of FIG. is the drawing shown. 10 and 11 are diagrams showing an operating state of the reclassification module of FIG. 3 , and FIG. 12 is a diagram showing an operating state of the type setting module of FIG. 3 . 13 and 14 are diagrams showing an operating state of the expression module of FIG. 3 .

The map generation module 20 extracts location information from the location data received by the communication module 10, and then creates a grid map C in which a grid area D is created in the extracted location information. For example, as shown in FIG. 4, the map generation module 20 creates a world map, forms a plurality of horizontal lines and a plurality of vertical lines at regular intervals on the generated world map, and generates a plurality of grid areas D. Create a grid map (C). An aerosol may be generated on the grid area D.

The aerosol detection module 30 detects aerosols from the aerosol optical thickness data received by the communication module 10 and calculates the number of aerosol pixels by calculating the detected aerosols according to pixels. As an example, the aerosol detection module 30, as shown in FIGS. 3 and 4, aerosol ( I) to detect Here, the aerosol optical depth (AOD) data is a value obtained by vertically integrating an aerosol dissipation coefficient, that is, a constant representing the degree of light absorption by the aerosol. Accordingly, there is no unit in the value of such aerosol optical thickness. The aerosol detection module 30 may calculate the number of aerosol pixels by calculating the number of aerosols based on the aerosol optical thickness data located in pixel units.

As shown in FIG. 6, the area setting module 40 sets the set area (E, blue solid line area) and the overseas emission source area (H, red dotted line area) in the grid map (C) set in the map generation module 20. create As an example, as shown in FIG. 6, the region setting module 40 sets the region of latitude 32 ° to 40 ° and longitude 122 ° to 130 ° in the grid map (C) as an area E, latitude 28 ° to 130 °. An area of 48°, longitude of 115° to 122.5°, latitude of 42° to 48°, and longitude of 122.5° to 128° may be created as the overseas emission source area (H). At this time, the created setting area E may be an area including a part of the territory of the Republic of Korea, and the overseas emission source area H may be an area including a territory outside the territory of the Republic of Korea.

In addition, the area division module 50 divides the set area E into a plurality of gateway areas F and at least one domestic area G. That is, the region division module 50 creates a plurality of gateway regions F and at least one domestic region G. For example, as shown in FIG. 6, the area division module 50 divides the first gateway area F1, the second gateway area F2, and the third gateway area F3 to the west of the setting area E. and create a fourth gateway area (F4) and a fifth gateway area (F5) in the north of the setting area (E), and the first gateway area (F1) to the fifth gateway area (F5) in the setting area (E) The domestic area (G) is created in the area except for .

The variable value calculation module 60 calculates the median value of the aerosol optical thickness ( m ), the standard deviation of the aerosol optical thickness ( σ ), and the average value of the aerosol optical thickness within the setting area (E) from the aerosol optical thickness data received during the setting period. Calculate ( μ ), respectively. That is, the variable value calculation module 60 calculates the median value ( m ), standard deviation value ( σ ), and average value ( μ ) of the aerosol optical thickness. Then, the variable value calculation module 60 generates and stores a plurality of equations using the calculated values as variables. As an example, the variable value calculation module 60, as shown in FIG. 7, has a first formula ( m -3σ ) or a median value ( m ) + 3 * standard deviation value having an expression of median value - 3 * standard deviation value After creating the 2nd formula (m+3σ ) and 3rd formula (m+3σ ) with the formula of ( σ ), it can be saved.

The background concentration setting module 65 calculates the median value ( m ) of the aerosol optical thickness calculated by the variable value calculation module 60 in the first preset formula ( m-3σ ) generated by the variable value calculation module 60 and The first formula is calculated by substituting the standard deviation ( σ ) of the aerosol optical thickness. And, as shown in FIG. 8 , the background concentration setting module 65 calculates an average value μ of the first formula values calculated during the setting period. Then, the average value ( μ ) of the calculated first formula is set as the background density of the grid area (D). Here, the setting period may be the current time or a period up to 15 days from the current time, that is, 15 days.

The initial classification module 70 compares the number of aerosol pixels in the setting area E with the reference number during the setting period. Thereafter, when the number of aerosol pixels in the setting area E exceeds the reference number during the setting period, a first classification process is determined so that the aerosol is classified through the first classification process. And, if the number of aerosol pixels in the setting area E is less than the reference number during the setting period, a second classification process is determined so that the aerosol is classified through the second classification process. For example, when the standard number is 50 and the number of aerosol pixels is 70, as shown in FIG. 8 , the initial classification module 70 classifies the aerosol into the first classification process. And when the number of aerosol pixels is 30, the initial classification module 70 classifies the aerosol into the second classification process.

When the initial classification module 70 determines the first classification process, the re-classification module 80 calculates the average value μ2 of the aerosol optical thickness within the plurality of gateway areas F during the set period through the variable value calculation module 60. yield As an example, as shown in FIGS. 10 and 11 , the average value ( μ 2 ) of the aerosol optical thicknesses of the first to fifth gateway regions F for 15 days is calculated. In addition, the reclassification module 80 calculates the intermediate value of the aerosol optical thickness within the plurality of gateway areas F1 to F5 during the set period according to the second formula ( m2+3σ2 ) generated and set through the variable value calculation module 60. ( m2 ) And the standard deviation ( σ2 ) of the aerosol optical thickness is substituted to calculate the boundary value of the gateway region. Then, the average value ( μ2 ) of the aerosol optical thickness within the calculated plurality of gateway areas (F1 to F5) is compared with the gateway area boundary value. At this time, when the average value ( μ 2 ) of the aerosol optical thickness in the plurality of gateway areas (F1 to F5) exceeds the gateway area boundary value, the 1-1 classification process is determined.

On the other hand, when the average value ( μ 2 ) of the aerosol optical thickness within the plurality of gateway regions F is less than or equal to the boundary value of the gateway region, the first-second classification process is determined. In addition, when the first classification process is determined, the re-classification module 80 checks the number of aerosol pixels in the plurality of gateway areas F calculated by the variable value calculation module 60, and if the number is 0, the re-classification module 80 determines the number of aerosol pixels. 1-3 Determine the classification process.

As shown in FIG. 12, the type setting module 90 may classify the detected aerosol into long-distance travel (LRT) type, non-long-distance travel (Non-LRT) type, clean type, undetermined type, and cloud type. . More specifically, the type setting module 90 reversely tracks the moving direction of the aerosol when the reclassification module 80 determines the 1-1 classification process, so that the aerosol is transferred to the overseas emission source region H and the gateway within a preset time. When passing through the area (F), the aerosol is determined as a long-distance travel (LRT) type. In addition, when the aerosol does not pass through the overseas emission source region (H) and the gateway region (F) within the preset time, the aerosol is determined as the first non-LRT type. In addition, the type setting module 90, when the reclassification module 80 determines the first-second classification process, uses the variable value calculation module 60 to set the average value of the aerosol optical thickness within the domestic area G during the setting period ( μ ) is calculated. And the median value ( m3 ) of the aerosol optical thickness within the domestic area (G) and the standard deviation of the aerosol optical thickness ( Substitute σ3 ) to calculate the domestic territory boundary value.

In addition, the type setting module 90 compares the calculated average value ( μ 3 ) of the aerosol optical thickness of the domestic region (G) and the boundary value of the domestic region, so that the average value ( μ 3 ) of the aerosol optical thickness of the domestic region (G) is the domestic region If the threshold value is exceeded, the aerosol is determined as the second non-LRT type. In addition, the type setting module 90 determines the aerosol as a clean type that does not affect the domestic area G when the average value ( μ 3 ) of the aerosol optical thickness of the domestic area G is equal to or less than the domestic area boundary value. In addition, the type setting module 90 determines the aerosol type as a cloud in the airspace of the domestic area G when the initial classification module 70 determines the second classification process. In addition, the type setting module 90 may determine the aerosol to be in an unknown state when the reclassification module 80 determines the first to third classification processes. At this time, the long-distance travel (LRT) type means a type of aerosol that is generated abroad and affects the country. And the non-LRT type refers to a type of aerosol that is generated in Korea and affects the country. And the clean type means that there are no aerosols in the airspace of the domestic area (G), and the cloud type means that clouds in the airspace of the domestic area (G) are detected as aerosols. . And the undetermined type means that the aerosol cannot be detected beyond the detection area.

The expression module 100 includes a map generation module 20, an aerosol detection module 30, a region setting module 40, a region division module 50, a variable value calculation module 60, a background concentration setting module 65, It is connected to the initial classification module 70, the reclassification module 80, and the type setting module 90, and the result output from each module and the type determined by the type setting module 90 can be displayed at the bottom of the screen. For example, if the type setting module 90 determines the long-distance movement (LRT) type, the display module 100 can output 'LRT' to the lower right of the screen as shown in (a) of FIG. . And when the non-long-distance movement (Non-LRT) type is determined in the type setting module 90, the expression module 100 outputs 'Non-LRT' at the bottom right of the screen as shown in FIG. 13(b). can do.

When the type setting module 90 determines the clean type, the display module 100 can output 'Clear' to the lower right corner of the screen as shown in (c) of FIG. 14 . In addition, when the cloud type is determined by the type setting module 90, the expression module 120 can output 'Cloud' at the bottom right of the screen as shown in FIG. 14(d). According to an embodiment of the present invention, the system 1 for classifying the impact type of aerosol through environmental satellites processes data by organically connecting a plurality of modules as described above, and the satellite A photographs the earth B, From the generated meteorological data, it detects the presence of aerosols around the Korean Peninsula and extracts path information of the inflowing aerosols. And it can respond to the inflowing aerosol.

Hereinafter, with reference to FIG. 15, a method for classifying impact types of aerosols through environmental satellites according to an embodiment of the present invention will be described in detail.

15 is a flow chart of a method for classifying impact types of aerosols through environmental satellites according to an embodiment of the present invention.

In the method for classifying impact types of aerosols through environmental satellites according to an embodiment of the present invention, components constituting the system for classifying impact types of aerosols through environmental satellites according to an embodiment of the present invention are subject to a series of orders. It goes on. Accordingly, the main body performing each step of the method for classifying the impact type of aerosol through environmental satellites according to an embodiment of the present invention and the components of the system for the impact type of aerosol through environmental satellites according to an embodiment of the present invention are the same. can do.

Therefore, the description of the aerosol impact type system through environmental satellites according to an embodiment of the present invention can be applied as it is to the aerosol impact type method description through environmental satellites according to an embodiment of the present invention.

According to an embodiment of the present invention, the aerosol impact type method through the environmental satellite includes (a) receiving location data and aerosol optical thickness data (S10), (b) generating a grid map (S20), and aerosol After detecting , (c) step of calculating the number of aerosol pixels (S30), (d) step of generating a set area (S40), (e) step of dividing the gateway area and domestic area (S50), and aerosol optics (f) step (S60) of calculating the median value of the thickness, the standard deviation of the aerosol optical thickness, and the average value of the aerosol optical thickness, and the type of aerosol is initially classified according to the number of aerosol pixels in the first or second classification process. (g) step (S70), (h) step (S80) of reclassifying the initially classified type, and (i) step (S90) of setting the reclassified aerosol type.

Hereinafter, steps (a) to (i) will be described in more detail. Step (a) is a step (S10) in which the communication module 10 receives location data and aerosol optical thickness data from the satellite A. In step (b), the map generation module 20 extracts location information from the location data received from the communication module 10, and then generates a grid map C in which a grid area D is created in the extracted location information. It becomes a step (S20) of doing. And in step (c), the aerosol detection module 30 detects aerosols from the aerosol optical thickness data received from the communication module 10, and calculates the number of aerosol pixels by calculating the detected aerosols according to pixels (S30). Step (d) is a step (S40) in which the region setting module 40 creates a region of a certain size as a setting region E in the grid map C set in the map generation module 20. Step (e) is a step (S50) in which the area division module 50 divides the set area (E) into a domestic area (G) and a plurality of gateway areas (F). And step (f) calculates the median value ( m ) of the aerosol optical thickness within the setting area (E) and the standard deviation ( σ ) of the aerosol optical thickness from the aerosol optical thickness data received by the variable value calculation module 60 during the setting period. And a step (S60) of calculating the average value ( μ ) of the aerosol optical thickness. And in step (g), the initial classification module 70 compares the number of aerosol pixels in the setting area (E) with the reference number during the setting period, and if the number of aerosol pixels in the setting area (E) exceeds the reference number during the setting period, the number of aerosol pixels in the setting area (E) is exceeded. When the first classification process is determined and the number of aerosol pixels in the setting area (E) is less than the reference number during the setting period, the second classification process is determined (S70), and step (h) is the initial stage when the reclassification module 80 When the first classification process is determined in the classification module 70, the average value ( μ 2 ) of the aerosol optical thickness in the plurality of gateway areas F is calculated through the variable value calculation module 60 during the set period, and the variable value calculation module The median value ( m2 ) of the aerosol optical thickness within the plurality of gateway areas (F) calculated during the set period through (60) and the standard deviation value (σ2) of the aerosol optical thickness are calculated by the second formula ( m2+3σ2 ) The average value of the aerosol optical thickness within the plurality of gateway regions F is calculated by substituting the gateway region boundary value, and the calculated average value ( μ2 ) of the aerosol optical thickness within the plurality of gate regions F is compared with the gateway region boundary value. If ( μ2 ) exceeds the gateway area boundary value, the 1-1 classification process is determined, and if the average value ( μ 2 ) of the aerosol optical thickness within the plurality of gateway areas (F) is less than the gateway area boundary value, the 1-2 classification process is determined. It becomes the step (S80) of determining a process.

And in step (i), when the type setting module 90 determines the 1-1 classification process in the re-classification module 80, the aerosol moves in a reverse direction to track the aerosol to the gateway area F within a preset time. When passing through, the aerosol is determined as a long-distance movement (LRT) type, and when the aerosol does not pass through the gateway area (F) within a preset time, the aerosol is determined as the first non-long-distance movement (Non-LRT) type, , When the 1-2 classification process is determined in the re-classification module 80, the average value ( μ 3 ) of the aerosol optical thickness within the domestic area (G) is calculated through the variable value calculation module 60 during the set period, and the preset In the third formula ( m3 + 3σ3 ), the median value ( m3 ) of the aerosol optical thickness within the domestic area (G) calculated during the set period through the variable value calculation module 60 and the standard deviation value ( σ3 ) of the aerosol optical thickness The domestic area boundary value is calculated by substitution, and the average value ( μ3 ) of the aerosol optical thickness of the domestic area (G) is compared with the calculated average value ( μ3 ) of the aerosol optical thickness of the domestic area (G). If the domestic area boundary value is exceeded, the aerosol is determined as the second non-LRT type, and if the average value ( μ3 ) of the aerosol optical thickness in the domestic area (G) is below the domestic area boundary value, the aerosol is classified as the domestic area. If it is determined as a clean type that does not affect (G) and the second classification process is determined in the initial classification module 70, determining the type of aerosol as a cloud in the airspace of the domestic area (G) (S90). At this time, in step (h), the reclassification module 80 checks the number of aerosol pixels in the plurality of gateway areas F calculated by the variable value calculation module 60, and when the number of confirmed aerosol pixels is 0, the first - The step (i) further includes the step of determining the 3 classification process, and in step (i), when the step of determining the 1-3 classification process in step (h) by the type setting module 90 proceeds, the aerosol is in an unknown state. It further includes the step of determining.

The method for classifying the impact type of aerosols through environmental satellites according to an embodiment of the present invention proceeds from steps (a) to (i) as a series of steps, and data is obtained through organic interlocking of subjects in each step. It processes quickly, detects the existence of aerosols around the Korean Peninsula from meteorological data through satellite (A), and extracts route information of aerosols flowing into the Korean Peninsula.

Hereinafter, with reference to FIGS. 15 to 18, a system and method for classifying the impact type of aerosol through environmental satellites according to another embodiment of the present invention will be described in detail.

The contents of the aerosol impact type classification system and method through environmental satellites according to one embodiment of the present invention are substantially the same as those of the aerosol impact type classification system and method through environmental satellites of another embodiment of the present invention described above. same.

Therefore, in order to simplify the description of the aerosol impact type classification system and method through environmental satellites of another embodiment of the present invention, the system and method of aerosol impact type classification through environmental satellites of one embodiment of the present invention are duplicated. For the part to be replaced with the description of the system and method for classifying the impact type of aerosol through environmental satellites according to an embodiment of the present invention, only the parts with differences will be described.

The differences between the aerosol impact type classification system and method through environmental satellites according to another embodiment of the present invention and the aerosol impact type classification system and method through environmental satellites according to an embodiment of the present invention are the vector aerosol generation module and the method. It is located in the movement calculation module.

Accordingly, the vector aerosol generating module 110 and the displacement calculation module 120 will be intensively described.

The vector aerosol generating module 110 generates wind vector aerosol J by applying wind power to the aerosol optical thickness data received by the communication module 10 and the aerosol detected by the aerosol detection module 30. As an example, the vector aerosol generating module 110, as shown in FIG. 17, has a size area of 0.2 ° (20 km) X 0.2 ° (20 km) of the grid area D in the aerosol detection module 30 and is shown in FIG. A wind vector aerosol (J) can be generated by applying wind force to the aerosol detected from the detected aerosol optical thickness (AOD) data. In addition, as shown in FIG. 18, the vector aerosol generating module 110 synthesizes a plurality of wind vector aerosols from the ground surface with a size of 20 km X 20 km to a height of 5.6 km (500 hPa) to form one synthesized wind vector aerosol (JA). ) can be created.

The movement amount calculation module 120 quantitatively indicates the amount of movement of the aerosol moving per time. For example, the movement amount calculation module 120 may indicate the movement amount of the wind vector aerosol generated in the vector aerosol generation module 110 .

The movement amount calculation module 120 calculates the median value ( mX ) of the aerosol optical thickness of the wind vector within the grid area (D) calculated by the variable value calculation module 60, the standard deviation value ( σX ) of the aerosol optical thickness, and the aerosol optics By substituting the average value of the thickness ( μX ) into the second formula or the third formula, a value greater than the value calculated by substituting, that is, a value exceeding the second formula value and the background concentration setting module 65 determines the grid area (D) Calculate the inflow and outflow per grid area generated by the map generation module using the difference value between the values of the background concentration set within the preset vertical variance equation. In this way, the movement amount calculation module 120 calculates the movement amount of the aerosol. At this time, the aerosol can be fine dust, and the Mass Extinction Efficiency (MEE) of this aerosol is a value determined by the aerosol type, size distribution, relative humidity, etc., which is not directly observed. It can be obtained from chemical and physicochemical properties. Therefore, it is known that the aerosol mass dissipation efficiency can be obtained through theoretical calculations or statistical analysis of observational data (Hand and Malm, 2007). The aerosol mass dissipation efficiency can be calculated by the following formula.

가 될 수 있다. 여기서, m(z): 질량수직분포 [g/m3]

가 되고, U(z): 고도별 바람 정보 [m/s]가 된다. 그리고 L: 격자 길이 [m] 가 되고, h: 단면적 높이[m]가 된다.In addition, the vertical integral variance equation is as shown in FIG.

can be where m(z): mass vertical distribution [g/m3]

, and U(z): wind information by altitude [m/s]. And L: lattice length [m], h: cross-section height [m].

The display module 100 is connected to the vector aerosol generation module 110 and the movement amount calculation module 120, and displays the direction of the wind and the movement of the aerosol on the screen and the concentration of the aerosol in the grid area D and further around the setting area. Allows determination of the amount of aerosol transported in an area.

In addition, the aerosol impact type classification system 1-1 through environmental satellites according to another embodiment of the present invention moves direction and movement of aerosols flowing in and out per grid area D through a vector aerosol generating module and a movement amount calculation module. can indicate the amount of

A method for classifying impact types of aerosols through environmental satellites according to another embodiment of the present invention includes steps (c-1) performed after step (c) of the method for classifying impact types of aerosols through environmental satellites according to an embodiment of the present invention; ) step (j) performed after the step is further included.

Here, step (c-1) is the wind vector aerosol (J ). In step (j), the movement amount calculation module 120 calculates the wind vector within the grid area D calculated by the variable value calculation module 60 ( mX ) and the standard deviation ( σX ) of the aerosol optical thickness. ) and the average value of the aerosol optical thickness ( μX ) into the second formula or the third formula, and a value greater than the calculated value and the background concentration value set within the grid area (D) through the background concentration setting module 65 It is a step of calculating the movement amount of the aerosol by substituting the difference value of the pre-set vertical variance equation.

A method for classifying impact types of aerosols through environmental satellites according to another embodiment of the present invention includes steps included in the method for classifying impact types of aerosols through environmental satellites according to an embodiment of the present invention, steps (c-1), and steps (j). Including, it is possible to check the direction of wind and the movement direction of aerosol on the screen, the concentration of aerosol in the grid area D, and the amount of aerosol moved in the surrounding area of the setting area.

Although the embodiments of the present invention have been described with reference to the accompanying drawings, those skilled in the art to which the present invention pertains can be implemented in other specific forms without changing the technical spirit or essential features of the present invention. you will be able to understand Therefore, it should be understood that the embodiments described above are illustrative in all respects and not restrictive.

1, 1-1: Aerosol impact type classification system through environmental satellites
10: communication module 20: map generation module
30: aerosol detection module 40: area setting module
50: area division module 60: variable value calculation module
65: background density setting module 70: initial classification module
80: reclassification module 90: type setting module
100: expression module 110: vector aerosol generating module
120: movement amount calculation module
A: Satellite B: Earth
C: grid map D: grid domain
E: setting area F: gateway area
G: Domestic area H: Overseas emission source area
I: aerosol J: wind vector aerosol

Claims (5)

A communication module 10 for receiving location data and aerosol optical depth (AOD) data from the satellite A;
a map generating module 20 for generating a grid map C in which a grid area D is generated based on the location data received by the communication module 10;
an aerosol detection module 30 that detects aerosols from the aerosol optical thickness data received by the communication module 10 and calculates the number of aerosol pixels by calculating the detected aerosols according to pixels;
an area setting module 40 for generating a setting area E in the grid map C set in the map generating module 20;
Area division module 50 that divides the set area (E) into a domestic area (G) and a plurality of gateway areas (F);
Through the number of aerosol pixels received from the aerosol detection module 30 during the setting period, the median value ( m ) of the aerosol optical thickness in the setting area (E), the standard deviation value ( σ ) of the aerosol optical thickness, and the average value of the aerosol optical thickness ( A variable value calculation module 60 that calculates μ );
Comparing the number of aerosol pixels in the setting area (E) with the reference number during the setting period, if the number of aerosol pixels in the setting area (E) exceeds the reference number during the setting period, the first classification process is determined, and the setting area (E) is determined during the setting period. An initial classification module 70 for determining a second classification process when the number of aerosol pixels in ) is less than the reference number;
When the initial classification module 70 determines the first classification process, the average value ( μ 2 ) of the aerosol optical thickness in the plurality of gate areas F is calculated through the variable value calculation module 60 during the set period, and the variable value is calculated. The median value ( m2 ) of the aerosol optical thickness within the plurality of gateway areas (F) calculated during the set period through the module 60 and the standard deviation value (σ2) of the aerosol optical thickness are calculated by the second formula ( m2+3σ2 ) Substituting in to calculate the gateway area boundary value, and comparing the calculated average value ( μ 2 ) of the aerosol optical thickness in the plurality of gateway areas (F) and the gateway area boundary value, the aerosol optical thickness in the plurality of gateway areas (F) If the average value ( μ2 ) exceeds the gateway area boundary value, the 1-1 classification process is determined, and if the average value ( μ 2 ) of the aerosol optical thickness within the plurality of gateway areas (F) is less than the gateway area boundary value, the 1-2 classification process is determined. a reclassification module 80 for determining a process; and
When the reclassification module 80 determines the 1-1 classification process, the moving direction of the aerosol is reversely tracked, and when the aerosol passes through the gateway area F within the preset time, the aerosol is classified as a long-distance movement (LRT) type. and when the aerosol does not pass through the gateway area F within the preset time, the aerosol is determined as the first non-LRT type, and the reclassification module 80 performs the first-second classification process. When is determined, the average value ( μ 3 ) of the aerosol optical thickness in the domestic area (G) is calculated through the variable value calculation module 60 during the setting period, and the domestic area calculated during the setting period through the variable value calculation module 60 ( G) The median value of the aerosol optical thickness ( m3 ) and the standard deviation of the aerosol optical thickness ( σ3 ) are substituted into the third formula ( m3 + 3σ3 ) to calculate the domestic area boundary value, Comparing the average value of the aerosol optical thickness ( μ3 ) of G) and the boundary value of the domestic region, if the average value of the optical thickness of the aerosol ( μ3 ) of the domestic region (G) exceeds the boundary value of the domestic region, the aerosol moves to the second non-long distance ( Non-LRT) type, and if the average value ( μ3 ) of the aerosol optical thickness in the domestic area (G) is less than the domestic area boundary value, the aerosol is determined as a clean type that does not affect the domestic area (G), and the initial When the classification module 70 determines the second classification process, the aerosol through the environmental satellite includes a type setting module 90 that determines the type of aerosol as a cloud in the airspace of the domestic area (G). of the impact type classification system. According to claim 1,
The reclassification module 80,
When the first classification process is determined, the number of aerosol pixels in the plurality of gateway areas (F) calculated by the variable value calculation module 60 is checked, and when the confirmed number is 0, the first-third classification process is determined,
The type setting module 90,
When the reclassification module 80 determines the first to third classification processes, the aerosol influence type classification system through the environmental satellite determines the number of aerosols as unknown. According to claim 1,
The first formula is calculated by substituting the median value ( m ) of the aerosol optical thickness and the standard deviation ( σ ) of the aerosol optical thickness calculated in the variable value calculation module 60 into the first formula ( m - 3σ ). and a background concentration setting module 65 for setting the average value of the first formula calculated during the setting period as the background concentration of the grid area D. (a) receiving, by the communication module 10, location data and aerosol optical thickness data from the satellite A (S10);
(b) generating, by the map generating module 20, a grid map C in which the grid area D is created with the location data received from the communication module 10 (S20);
(c) the aerosol detection module 30 detects aerosols from the aerosol optical thickness data received from the communication module 10 and calculates the number of aerosol pixels by calculating the detected aerosols according to pixels (S30);
(d) generating, by the area setting module 40, an area of a certain size as a setting area E in the grid map C set in the map generating module 20 (S40);
(e) dividing the set area (E) into a domestic area (G) and a plurality of gateway areas (F) by the area division module 50 (S50);
(f) From the aerosol optical thickness data received by the variable value calculation module 60 during the setting period, the median value ( m ) of the aerosol optical thickness within the setting area (E), the standard deviation value ( σ ) of the aerosol optical thickness, and the aerosol optical thickness Calculating an average value ( μ ) of thickness (S60);
(g) The initial classification module 70 compares the number of aerosol pixels in the setting area (E) with the reference number during the setting period, and if the number of aerosol pixels in the setting area (E) exceeds the reference number during the setting period, the first classification process Determining a second classification process when the number of aerosol pixels in the setting area (E) is less than the reference number during the setting period (S70);
(h) When the re-classification module 80 determines the first classification process in the initial classification module 70, the average value of the aerosol optical thickness within the plurality of gateway areas F during the set period through the variable value calculation module 60 ( μ2 ) is calculated, and the median value ( m2 ) of the aerosol optical thickness within the plurality of gateway areas (F) calculated during the setting period through the variable value calculation module 60 and the standard deviation value ( σ2 ) of the aerosol optical thickness are calculated. The gateway area boundary value is calculated by substituting it into the previously set second formula ( m2+3σ2 ) , and the average value (μ2) of the aerosol optical thickness in the calculated plurality of gateway areas (F) is compared with the gateway area boundary value, If the average value ( μ2 ) of the aerosol optical thickness in the gateway area (F) exceeds the gateway area boundary value, the 1-1 classification process is determined, and the average value ( μ2 ) of the aerosol optical thickness in the plurality of gateway areas (F) is determined as the gateway. determining a 1-2 classification process when the value is less than the region boundary value (S80); and
(i) When the type setting module 90 determines the 1-1 classification process in the reclassification module 80, the aerosol's movement direction is tracked in reverse and the aerosol passes through the gateway area F within a preset time. When the aerosol does not pass through the gateway area (F) within the preset time, the aerosol is determined as the first non-LRT type and reclassified. When the 1-2 classification process is determined in the module 80, the average value ( μ 3 ) of the aerosol optical thickness within the domestic area (G) during the set period is calculated through the variable value calculation module 60, and the third formula is set. ( m3 + 3σ3 ) by substituting the median value ( m3 ) of the aerosol optical thickness within the domestic area (G) calculated during the setting period through the variable value calculation module 60 and the standard deviation value ( σ3 ) of the aerosol optical thickness The area boundary value is calculated, and the calculated average value ( μ3 ) of the aerosol optical thickness of the domestic area (G) and the domestic area boundary value are compared, and the average value ( μ3 ) of the aerosol optical thickness of the domestic area (G) is the domestic area boundary value. If the value is exceeded, the aerosol is determined as the second non-LRT type, and if the average value ( μ3 ) of the aerosol optical thickness in the domestic area (G) is below the domestic area boundary value, the aerosol is classified as the domestic area (G). (S90) A method for classifying impact types of aerosols through environmental satellites, including a. According to claim 4,
In step (h),
The re-classification module 80 checks the number of aerosol pixels in the plurality of gateway areas F calculated by the variable value calculation module 60, and determines the 1-3 classification process when the checked number is 0 Including more,
In step (i),
When the step of determining the 1-3 classification process in step (h) by the type setting module 90 proceeds, further comprising the step of determining the aerosol as unknown, classifying the impact type of the aerosol through the environmental satellites. Way.

Images