KR102411506B1 - Apparatus and method for deriving deriving tropospheric ozone motion vector - Google Patents

Abstract

The present invention relates to an apparatus and method for calculating a tropospheric ozone movement vector, wherein the apparatus for calculating a tropospheric ozone movement vector according to an aspect of the present invention includes ozone concentration values previously matched to a plurality of pixels corresponding to adjacent latitudes and longitudes. Ozone vertical data, an acquisition unit for acquiring an ozone observation time point corresponding to the ozone vertical data, a target area setting unit for setting a target area from the ozone vertical data corresponding to the first time point, and a preset point before or after the first time point Based on the ozone vertical data corresponding to each time point, the matching region, which is the region with the highest cross-correlation coefficient with the target region, is detected, and displacement values for each time point are calculated based on the position difference between the target region and the matching region, respectively. and an ozone movement vector calculation module for calculating an ozone movement vector based on the displacement values and matching the calculated ozone movement vector with a target area.

Description

Apparatus and method for calculating tropospheric ozone movement vector

The present invention relates to an apparatus and method for calculating a tropospheric ozone movement vector, and more particularly, to a tropospheric ozone movement vector calculating apparatus capable of tracking tropospheric ozone movement using vertical ozone data calculated from ultraviolet observation values observed from satellites; it's about how

In the stratosphere, ozone plays a beneficial role in protecting the ecosystem from UV rays, but in the troposphere, long-term exposure to air pollutants harmful to the human body can cause skin cancer and lung damage. In particular, tropospheric ozone can act more fatally for children or the elderly with weakened immune systems.

Tropospheric ozone is currently increasing globally, and although the increase was remarkable in Europe in the past, the ozone concentration value in East Asia is increasing rapidly in recent years. This is related to an international trend such as the recent rapid industrialization in East Asia. In particular, as China became industrialized and major industrial complexes in China were built on the windward side of the Korean Peninsula, damage from air pollutants including fine dust and ozone from China has occurred. Accordingly, in order to efficiently manage air quality, the importance of technology that clearly identifies the cause is increasing.

In this regard, technologies for tracking the movement of air pollutants using satellite, ground, or aircraft observation data have been disclosed. There is this.

However, the conventional technology for tracking the movement of air pollutants using satellite observation data tracks the movement based on an image image in which an object such as a cloud is displayed, and it was difficult to track the movement of ozone by applying invisible ozone data .

In addition, even using the visualized ozone data, the conventional techniques only disclose two-dimensional data such as the direction and speed of movement by analyzing the movement based on the image image. There was a problem that it couldn't.

Korean Patent Application No. 10-2008-0072809 Republic of Korea Patent Application No. 10-2016-0150303

An object of the present invention is to solve the above problems, and a tropospheric ozone movement vector capable of tracking the movement of ozone in the tropospheric space using ozone data obtained from ozone vertical data calculated from ultraviolet observation values observed from artificial satellites. It is to provide a calculation apparatus and method.

Another object of the present invention is to provide an apparatus and method for calculating a tropospheric ozone movement vector capable of calculating the amount of ozone change in an arbitrary region based on ozone vertical data calculated through ultraviolet observation values observed from artificial satellites.

The object of the present invention is not limited to the object mentioned above, and other objects not mentioned will be clearly understood from the description below.

To achieve the above object, an apparatus for calculating a tropospheric ozone movement vector according to an aspect of the present invention includes ozone vertical data including ozone concentration values previously matched to a plurality of pixels corresponding to adjacent latitudes and longitudes, and ozone vertical data An acquisition unit for acquiring an ozone observation point corresponding to the first time point, a target area setting unit for setting a target area in the ozone vertical data corresponding to the first time point, and ozone vertical data corresponding to each preset time point before or after the first time point Detects the matching region, which is the region with the highest cross-correlation coefficient with the target region, based on and an ozone movement vector calculation module that calculates the vector and matches the calculated ozone movement vector with the target area.

A tropospheric ozone movement vector calculation method according to another aspect of the present invention includes ozone vertical data including ozone concentration values pre-matched to a plurality of pixels corresponding to adjacent latitudes and longitudes, and an ozone observation point corresponding to the ozone vertical data. a step of obtaining a, setting a target area in the ozone vertical data corresponding to the first time point, a cross-correlation coefficient with the target area based on the ozone vertical data corresponding to each of the preset time points before or after the first time point Detecting the matching region, which is the region with the highest , calculating displacement values for each time based on the position difference between the target region and the matching region, respectively, and calculating the ozone movement vector based on the displacement values for each time. and matching the vector with the target region.

According to the present invention, there is an effect of providing an apparatus and method for calculating a tropospheric ozone movement vector capable of tracking the movement of ozone in the troposphere using ozone vertical data calculated through ultraviolet observation values observed from artificial satellites.

In addition, according to the present invention, it is possible to calculate the ozone movement vector and calculate the ozone change amount or the ozone discharge amount for an arbitrary area using the ozone movement vector and the ozone concentration value.

The effects of the present invention are not limited to the above-mentioned effects, and other effects not mentioned will be clearly understood by those skilled in the art from the description of the claims.

1 is a block diagram of an apparatus for calculating a tropospheric ozone movement vector according to an embodiment of the present invention.
2 is a flowchart of a method for calculating a tropospheric ozone movement vector according to another embodiment of the present invention.
3A is a diagram illustrating an ozone movement vector calculated according to embodiments of the present invention.
3B is a scatter diagram comparing the ozone movement vector and the wind field vector calculated in the embodiments of the present invention.
4A is a diagram illustrating an ozone movement vector calculated from embodiments of the present invention.
4B is a diagram illustrating an ozone movement vector in which a quality factor of 0.6 or more is calculated according to embodiments of the present invention.
4C is a diagram illustrating an ozone movement vector in which a quality factor of 0.7 or higher is calculated according to embodiments of the present invention.
4D is a diagram illustrating an ozone movement vector in which a quality factor of 0.8 or more is calculated according to embodiments of the present invention.
5A is a graph showing a speed bias according to a threshold value.
5B is a graph illustrating a root mean square error according to a threshold value.

Advantages and features of the present invention and methods of achieving them will become apparent with reference to the embodiments described below in detail in conjunction with the accompanying drawings. However, the present invention is not limited to the embodiments disclosed below, but will be implemented in various different forms, and only these embodiments allow the disclosure of the present invention to be complete, and those of ordinary skill in the art to which the present invention pertains. It is provided to fully inform the person of the scope of the invention, and the present invention is only defined by the description of the claims. On the other hand, the terms used in this specification are for describing the embodiments, and are not intended to limit the present invention. In this specification, the singular also includes the plural unless specifically stated otherwise in the phrase.

1 is a block diagram of an apparatus for calculating a tropospheric ozone movement vector according to an embodiment of the present invention. Referring to FIG. 1 , a tropospheric ozone movement vector calculating device 10 according to an embodiment of the present invention includes an acquisition unit 101 , a target area setting unit 110 , an ozone movement vector calculation module 120 , and a wind field input unit ( 130), quality factor calculation unit 140, threshold value setting unit 150, ozone movement vector output unit 160, ozone flux calculation unit 170, monitoring area setting unit 180, ozone change amount calculation unit 190 ) may be included.

The acquisition unit 101 may acquire ozone vertical data including ozone concentration values previously matched to a plurality of pixels corresponding to adjacent latitudes and longitudes, and an ozone observation point corresponding to the ozone vertical data.

The ozone vertical data may include a vertical ozone concentration value for each predetermined vertical layer calculated from observation values of ultraviolet rays observed for a plurality of pixels corresponding to latitudes and longitudes adjacent to each other in the artificial satellite.

The acquisition unit 101 may have a spatial resolution corresponding to the latitude and longitude intervals according to the preset resolution and a temporal resolution of a preset time interval (eg, 1 hour) according to the observation period of the artificial satellite.

In addition, the ozone vertical data may include a vertical ozone concentration value in units of Dobson Unit [DU].

The target area setting unit 110 sets the target area in the ozone vertical data corresponding to the first time point.

The size of the target area may be preset according to a user's input.

The target area setting unit 110 selects a plurality of reference areas by dividing the ozone vertical data corresponding to the first time point at a predetermined interval, and is located within a predetermined radius based on each of a plurality of pixels included in the reference area. The standard deviation of the ozone concentration value matched to the pixel may be calculated, and the target area may be set by changing the central position of the reference area so as to correspond to the pixel for which the largest standard deviation is calculated.

Accordingly, the target area setting unit 110 may set the target area to include many characteristic ozone concentration values.

The ozone movement vector calculation module 120 detects a matching region, which is a region having the highest cross-correlation coefficient with the target region, based on the ozone vertical data corresponding to each of the preset time points before or after the first time point, and the target region The process may include calculating displacement values for each time point based on the difference in position of the matching area, calculating an ozone movement vector based on the displacement values for each time point, and matching the calculated ozone movement vector with the target area.

More specifically, the ozone movement vector calculation module 120 is configured to calculate the region having the highest cross-correlation coefficient with the target region in the ozone vertical data corresponding to each of the second time point before the first time point and the third time point after the first time point. A matching area is detected, and a plurality of displacement values for each of the second and third time points are calculated according to the difference between the center position of the target area and the center position of the matching area, and a plurality of displacement values calculated for the second time point are calculated. 1 vector obtained by dividing the value by the time between the first time point and the second time point, and a second vector obtained by dividing the plurality of displacement values calculated for the third time point by the time between the first time point and the third time point and calculating the ozone movement vector by averaging the first vector and the second vector to match the target area.

The ozone movement vector calculation module 120 includes a tracking area setting unit 121 , a matching area detection unit 123 , a displacement value calculation unit 125 , a vector generation unit 127 , and an ozone movement vector calculation unit 129 . it could be

The tracking area setting unit 121 has a size larger than the target area in the ozone vertical data corresponding to each of the second time point before the first time point, and the third time point after the first time point, and sets the same tracking area as the center position and the target area. may be doing

For example, when the target area has a size of 6X6 pixels, the tracking area may be set to have a size of 14X14 pixels, but is not limited thereto.

The matching area detection unit 123 is a part of the tracking area set from the ozone vertical data corresponding to each of the second and third time points for the target area, and a cross-correlation coefficient between a plurality of partial areas having the same size as the target area and the target area is calculated, and a partial region having the largest calculated cross-correlation coefficient among the plurality of partial regions may be detected as the matching region.

The matching region detection unit 123 may calculate a cross-correlation coefficient (CC) between the target region and a partial region having the same size as the target region as a part of the tracking region using the following equation.

는 추적영역과 표적영역에 속하는 화소에 매칭된 오존 농도값들의 공분산을 의미하며,

는 표적영역 내 화소에 매칭된 오존 농도값들의 표준편차를 의미하고,

는 추적영역 내에서 표적영역과 대응하는 화소에 대응되는 오존농도값들의 표준편차를 의미한다.where (m, n) means the relative position of the pixel in the tracking area,

is the covariance of ozone concentration values matched to pixels belonging to the tracking area and the target area,

is the standard deviation of ozone concentration values matched to pixels in the target area,

is the standard deviation of ozone concentration values corresponding to the pixel corresponding to the target area in the tracking area.

)는 아래의 수학식을 통해 산출되는 것일 수 있다.The standard deviation of ozone concentration values matched to pixels in the target area (

) may be calculated through the following equation.

는 표적영역 내에 i행 j열의 화소에 매칭된 오존 농도값이며,

는 표적영역 내 화소에 매칭된 오존 농도값들의 평균이다.Here, N is the pixel length of one side of the square target area,

is the ozone concentration value matched to the pixels in row i and column j in the target area,

is the average of ozone concentration values matched to pixels in the target area.

)는 아래의 수학식을 통해 산출되는 것일 수 있다.The standard deviation (

) may be calculated through the following equation.

는 추적영역 내에서 표적영역의 어느 한 화소와 대응되는 화소와 매칭된 오존 농도값이고,

는 추적영역에서 표적영역과 대응되는 화소들에 매칭된 오존 농도값들의 평균이다.here,

is the ozone concentration value matched with the pixel corresponding to any one pixel of the target area in the tracking area,

is the average of ozone concentration values matched to pixels corresponding to the target area in the tracking area.

The displacement value calculator 125 calculates a plurality of displacement values for each of the second and third time points according to the difference between the center position of the matching region and the center position of the target region detected for each of the second and third time points. It may be to calculate

The displacement value calculator 125 uses the latitude and longitude corresponding to the center position of the matching area detected for each of the second and third time points, the latitude and longitude corresponding to the center position of the target area, and the pre-stored earth radius. Thus, a plurality of displacement values for each of the second and third time points may be calculated.

,

)을 산출하는 것일 수 있다.The displacement value calculator 125 calculates a plurality of displacement values (

,

) can be calculated.

는 지구 반지름이고,

와

은 각각 표적영역 중심에 대응되는 화소에 매칭된 위도와 경도이며,

와

는 각각 정합영역의 중심에 대응되는 화소에 매칭된 위도와 경도이다.here,

is the radius of the earth,

Wow

is the latitude and longitude matched to the pixel corresponding to the center of each target area,

Wow

is the latitude and longitude matched to the pixel corresponding to the center of the matching area, respectively.

The vector generator 127 divides the plurality of displacement values calculated for the second time point by the time between the first time point and the second time point, respectively. The second vector may be generated based on values obtained by dividing the plurality of displacement values calculated for the third time point by the time between the first time point and the third time point, respectively.

The ozone movement vector calculator 129 may calculate the ozone movement vector by simply averaging the first vector and the second vector and match the ozone movement vector with the target area.

3A is a diagram illustrating the ozone movement vector calculated in the embodiments of the present invention, and FIG. 3B is a scatter diagram comparing the ozone movement vector and the wind field vector calculated in the embodiments of the present invention.

Referring to FIG. 3 , it can be confirmed that the ozone movement vector calculated according to the embodiments of the present invention is similar to the actual wind field, so that the ozone movement prediction is effectively performed.

The wind field input unit 130 may receive wind field data including wind field vectors previously matched to a plurality of pixels corresponding to adjacent latitudes and longitudes from the outside.

The quality factor calculating unit 140 is a quality factor for the ozone movement vector based on the difference between the wind field vector matched to the pixel corresponding to the target zero in the wind field data input from the outside, and the ozone movement vector matched to the target area. may be to calculate

More specifically, the quality factor calculator 140 calculates a first quality value according to an angle difference between a first vector and a second vector generated with respect to the target region, and a difference in magnitude between the absolute values of the first vector and the second vector. The second quality value, the third quality value according to the vector difference between the first vector and the second vector, and the fourth quality value according to the vector difference between the ozone movement vector and the wind field vector are respectively calculated, and the first quality value and the second quality value are calculated. The quality factor may be calculated by multiplying and summing each of the quality value, the third quality value, and the fourth quality value by a preset weight.

The quality factor calculation unit 140 calculates the quality factor based on the first quality value, the second quality value, the third quality value, and the fourth quality value, so that the vectors calculated at each time point before and after the first time point are It is possible to calculate the quality factor reflecting the consistency of wind direction, wind speed, vector, and the consistency of ozone movement vector and wind field.

In addition, the quality factor calculator 140 further calculates a fifth quality value according to the vector difference between the ozone movement vector for the target target area, which is one target area, and the ozone movement vector matched to another target area adjacent to the target target area, and , the first to fourth quality values and the fifth quality values may be multiplied and summed by a preset weight to calculate the quality factor.

Accordingly, the quality factor calculator 140 may calculate a quality factor that further reflects the spatial homogeneity of the ozone movement vector.

The quality factor calculator 140 may calculate the first quality value using the following equation.

Here, Diff is the angle difference between the first vector and the second vector, Avg is the average angle of the first vector and the second vector, and coeff.A, coeff.B, coeff.C, and coeff.D are preset constant values. .

The quality factor calculator 140 may calculate the second to fourth quality values using the following equation, respectively.

Here, Diff is the absolute difference between the first vector and the second vector in the calculation of the second quality value, the vector difference between the first vector and the second vector in the calculation of the third quality value, and the ozone movement vector and the wind field in the calculation of the fourth quality value The vector difference with the vector is the vector difference between the ozone movement vector and the ozone movement vector matched to another target area adjacent to the target area where the ozone movement vector and the ozone movement vector are matched in the calculation of the fifth quality value, and Avg is the first vector in the calculation of the second quality value and the absolute value average of the second vector, the vector average of the first and second vectors in the third quality value calculation, the vector average of the ozone movement vector and the wind field vector in the fourth quality value calculation, and the fifth quality value calculation The ozone movement vector and the ozone movement vector are the vector averages of the ozone movement vectors matched to the matched target area and other adjacent target areas, and coeff.A, coeff.B, coeff.C, and coeff.D are preset constant values.

Table 1 below shows constant values (coeff.A, coeff.B, coeff.C, and coeff.D) applied to the calculation of the first to fifth quality values.

coeff.A oeff.B coeff.C coeff.D first quality value 20 10 10 4*10 ⁷ second quality value 0.2 0.01 5 4 third quality value 0.3 0.01 5 40 4th quality value 0.4 0.01 5 15 5th quality value One 0.01 One 2

The threshold value setting unit 150 may set a threshold value according to a user input.

The threshold value setting unit 150 may set a preset initial value (eg, 0) as a threshold value, and set a value input according to a user's input as the threshold value.

The ozone movement vector output unit 160 may output an ozone movement vector matched to a target area in which a quality factor equal to or greater than a preset threshold value is calculated among a plurality of ozone movement vectors matched to each of the plurality of target areas.

Referring to FIG. 4 , it can be seen that although the number of output ozone movement vectors decreases as the threshold value increases, the overall distribution of ozone movement vectors does not change.

Figure 5a is a graph showing the speed bias (Speed-bias) according to the threshold value, Figure 5b is a graph showing the root mean square error (RMSE) according to the threshold value, referring to Figure 5, the threshold value is As it increases from 0 to 0.9, both the speed bias and the root mean square error become smaller. It can be seen that by applying the quality factor, the ozone movement vector with large error can be effectively separated from the ozone movement vector with small error.

According to the present invention, only the ozone movement vector for which the quality factor equal to or greater than the preset threshold is calculated is output, so that the high quality ozone movement vector is output.

The ozone flux calculator 170 may calculate the ozone flux for the target area by multiplying the ozone movement vector matched to the target area by the ozone concentration value previously matched to the pixels included in the target area.

The ozone flux calculator 170 calculates the first flux value obtained by multiplying the ozone concentration value by the cosine component value (for example, the x-axis change amount) of the ozone movement vector matched to the target region, and the sine component value of the ozone movement vector ( For example, the ozone flux including the second flux value obtained by multiplying the y-axis change amount by the ozone concentration value may be calculated.

The monitoring area setting unit 180 may set a monitoring area including a plurality of pixels corresponding to latitude and longitude adjacent to each other according to a user's input.

The ozone change calculation unit 190 extracts target areas located within a predetermined radius of the monitoring area from among the plurality of target areas set in the target area setting unit, and a predetermined first direction adjacent to the monitoring area among the plurality of extracted target areas. Ozone fluxes for target areas corresponding to (for example, the left or lower side of the monitoring area) are summed, and a preset second direction different from the first direction adjacent to the monitoring area among the plurality of extracted target areas is added. The amount of change in ozone for the monitoring area may be calculated by subtracting the ozone flux for the target areas located corresponding to the direction (eg, to the right or above the monitoring area).

The preset radius in the monitoring area means within a preset pixel length range based on the edge of the monitoring area.

The ozone change calculation unit 190 extracts target areas located within a preset radius of the monitoring area from among the plurality of target areas set in the target area setting unit, and a preset first area adjacent to the monitoring area among the plurality of extracted target areas. The first component values of the ozone flux for target areas located at (eg, the left side with respect to the monitoring area) are summed, and a preset second area different from the first area (eg, the right side with respect to the monitoring area) is added. ) subtracting the first component value of the ozone flux for the target areas located in Amount of ozone change for the monitoring area by summing the values and subtracting the value of the second component of the ozone flux for target areas located in a preset fourth area different from the third area (eg, above the monitoring area) may be to yield

The example described above is described based on the ozone inflow, and contrary to the above example (for example, summing as subtraction and subtraction as summation), the ozone outflow for the monitoring area can be calculated as the calculation is carried out.

That is, in embodiments of the present invention, the ozone change calculation unit 190 includes an ozone movement vector matched to a target area located in the vicinity of a monitoring area set according to a user's input, and ozone matched to pixels corresponding to the target area. Using the concentration value, it is possible to calculate the amount of ozone inflow or outflow into the monitoring area, that is, the change in ozone in the monitoring area.

For example, when a monitoring area surrounding the Korean peninsula is set, the ozone change amount calculating unit 190 can quantitatively calculate the ozone change amount in the Korean peninsula by multiplying the ozone movement vector corresponding to the edge of the monitoring area and the ozone concentration value.

2 is a flowchart of a method for calculating a tropospheric ozone movement vector according to another embodiment of the present invention.

The method for calculating a tropospheric ozone movement vector according to another embodiment of the present invention may be performed through the tropospheric ozone movement vector calculating apparatus 10 according to an embodiment of the present invention. Hereinafter, the tropospheric ozone movement vector calculating apparatus described above Contents and configurations overlapping with (10) coincide with reference numerals and detailed descriptions are omitted for convenience of description.

The tropospheric ozone movement vector calculating device 10 acquires ozone vertical data including ozone concentration values previously matched to a plurality of pixels corresponding to latitude and longitude adjacent to each other from the outside, and an ozone observation point corresponding to the ozone vertical data. and (S103), a target area is set in the ozone vertical data corresponding to the first time point (S103).

Thereafter, the tropospheric ozone movement vector calculating device 10 sets a tracking area in the ozone vertical data corresponding to each of the second time point before the first time point and the third time point after the first time point (S105), and the target in the tracking area A matching region having the highest cross-correlation coefficient with the region is detected (S107), and a plurality of displacement values for each of the second and third time points are calculated based on the central position of the matching region and the center position of the target region (S109). ), calculates an ozone movement vector for the target area based on a plurality of displacement values for each of the second and third time points and matches them with the target area (S111).

Here, the tropospheric ozone movement vector calculating device 10 converts the plurality of displacement values for each of the second time point and the third time point into the time between the first time point and the second time point and the time between the first time point and the third time point, respectively. The ozone movement vector may be calculated by generating the first vector and the second vector, respectively, based on the divided values, and simply averaging the first vector and the second vector.

Thereafter, the tropospheric ozone movement vector calculating device 10 receives wind field data including a wind field vector previously matched to a plurality of pixels corresponding to latitude and longitude adjacent to each other from the outside (S113), A quality factor is calculated based on the first vector, the second vector, the ozone movement vector, the wind field vector matched to the pixel corresponding to the target area, and the ozone movement vector matched to another target area adjacent to the target area ( S115), an ozone movement vector matched to a target area for which a quality factor equal to or greater than a preset threshold is calculated according to a user's input is output (S117).

The tropospheric ozone movement vector calculating device 10 calculates the ozone flux by multiplying the output ozone movement vector by the ozone concentration value of the target area matched to the ozone movement vector (S119), and a monitoring area preset according to a user's input The ozone flux of the target area located nearby is added up or subtracted depending on whether the target area is located in a predetermined first direction or a second direction with respect to the monitoring area to calculate the ozone change amount in the monitoring area (S121).

According to embodiments of the present invention, an ozone movement vector may be calculated using ozone concentration values for a predetermined vertical layer corresponding to latitude and longitude calculated from ultraviolet observation values observed from artificial satellites at regular time intervals.

Accordingly, the movement of ozone is intuitively grasped through the visualization of the ozone movement vector.

In addition, according to the embodiments of the present invention, there is an advantage that the amount of ozone change according to the inflow and outflow of ozone to an arbitrary area can be quantitatively traced using the ozone movement vector and the ozone concentration value.

Those of ordinary skill in the art to which the present invention pertains will understand that the present invention may be embodied in other specific forms without changing the technical spirit or essential features thereof. Therefore, it should be understood that the embodiments described above are illustrative in all respects and not restrictive. The scope of the present invention is indicated by the claims to be described later rather than the above detailed description, and all changes or modifications derived from the claims and their equivalent concepts should be construed as being included in the scope of the present invention.

10: tropospheric ozone movement vector calculation device
101: Acquisition
110: target area setting unit
120: ozone movement vector calculation module
130: wind field input unit
140: quality factor calculation unit
150: threshold setting unit
160: ozone movement vector output unit
170: ozone flux calculator
180: monitoring area setting unit
190: ozone change amount calculator

Claims (12)

an acquisition unit configured to acquire ozone vertical data including ozone concentration values previously matched to a plurality of pixels corresponding to adjacent latitudes and longitudes, and an ozone observation point corresponding to the ozone vertical data;
a target area setting unit for setting a target area in the ozone vertical data corresponding to the first time point;
A matching region, which is a region having the highest cross-correlation coefficient with the target region, is detected based on ozone vertical data corresponding to each of the preset time points before or after the first time point, and the positions of the target region and the matching region an ozone movement vector calculating module that calculates displacement values for each time point based on the difference, calculates an ozone movement vector based on the displacement values for each time point, and matches the calculated ozone movement vector with the target area;
an ozone flux calculator configured to calculate an ozone flux for the target region by multiplying the ozone movement vector matched to the target region by an ozone concentration value previously matched to a pixel included in the target region;
a monitoring area setting unit configured to set a monitoring area including a plurality of pixels corresponding to latitude and longitude adjacent to each other according to a user's input; and
Extracting target areas located within a predetermined radius in the monitoring area from among the plurality of target areas set in the target area setting unit, and located corresponding to a predetermined first direction adjacent to the monitoring area among the plurality of extracted target areas By summing the ozone flux for the target areas, and subtracting the ozone flux for the target areas located corresponding to a second predetermined direction different from the first direction adjacent to the monitoring area among the plurality of extracted target areas A tropospheric ozone movement vector calculating device comprising a; an ozone change amount calculating unit that calculates the ozone change amount for the monitoring area accordingly. delete delete an acquisition unit configured to acquire ozone vertical data including ozone concentration values previously matched to a plurality of pixels corresponding to adjacent latitudes and longitudes, and an ozone observation point corresponding to the ozone vertical data;
a target area setting unit for setting a target area in the ozone vertical data corresponding to the first time point; and
A matching region, which is a region having the highest cross-correlation coefficient with the target region, is detected based on ozone vertical data corresponding to each of the preset time points before or after the first time point, and the positions of the target region and the matching region an ozone movement vector calculation module that calculates displacement values for each time point based on the difference, calculates an ozone movement vector based on the displacement values for each time point, and matches the calculated ozone movement vector with the target area;
The target area setting unit
A plurality of reference areas are selected by dividing the ozone vertical data corresponding to the first time point at a predetermined interval, and ozone concentrations matched to pixels located within a predetermined radius based on each of the plurality of pixels included in the reference area Calculating the standard deviation of the values, and setting the target area by changing the central position of the reference area so as to correspond to the pixel for which the largest standard deviation is calculated
Phosphorus Tropospheric Ozone Transport Vector Calculator. an acquisition unit configured to acquire ozone vertical data including ozone concentration values previously matched to a plurality of pixels corresponding to adjacent latitudes and longitudes, and an ozone observation point corresponding to the ozone vertical data;
a target area setting unit for setting a target area in the ozone vertical data corresponding to the first time point;
A matching region, which is a region having the highest cross-correlation coefficient with the target region, is detected based on ozone vertical data corresponding to each of the preset time points before or after the first time point, and the positions of the target region and the matching region an ozone movement vector calculating module that calculates displacement values for each time point based on the difference, calculates an ozone movement vector based on the displacement values for each time point, and matches the calculated ozone movement vector with the target area;
Wind field data including a wind field vector previously matched to a plurality of pixels corresponding to adjacent latitudes and longitudes is received from the outside, and the wind field vector matched to the pixel corresponding to the target zero and the target area are applied to the target area. a quality factor calculator for calculating a quality factor for the ozone motion vector based on a difference between the matched ozone motion vectors; and
Tropospheric ozone movement vector calculation comprising a; Device. 6. The method of claim 5,
The ozone movement vector calculation module is
In the ozone vertical data corresponding to each of the second time point before the first time point and the third time point after the first time point, a matching area, which is an area having the highest cross-correlation coefficient with the target area, is detected, and A plurality of displacement values for each of the second time point and the third time point are calculated according to a difference between the center position and the center position of the matching region, and the plurality of displacement values calculated for the second time point are used as the first time point. and a first vector obtained by dividing by the time between the second time point and a second vector obtained by dividing a plurality of displacement values calculated for the third time point by the time between the first time point and the third time point generating and calculating an ozone movement vector by averaging the first vector and the second vector to match the target area,
The quality factor calculation unit
A first quality value according to an angle difference between the first vector and the second vector, a second quality value according to a difference in magnitudes of the absolute values of the first vector and the second vector, the first vector and the second vector A third quality value according to a vector difference between vectors and a fourth quality value according to a vector difference between the ozone movement vector and the wind field vector are calculated, respectively, and the first quality value, the second quality value, and the third quality Calculating a quality factor by multiplying each value and the fourth quality value by a predetermined weight and adding them together
Phosphorus Tropospheric Ozone Transport Vector Calculator. acquiring vertical ozone data including ozone concentration values previously matched to a plurality of pixels corresponding to adjacent latitudes and longitudes, and an ozone observation time point corresponding to the ozone vertical data;
setting a target area in the ozone vertical data corresponding to the first time point;
detecting a matching region, which is a region having the highest cross-correlation coefficient with the target region, based on the ozone vertical data corresponding to each preset time point before or after the first time point;
calculating displacement values for each time point based on a position difference between the target region and the matching region;
calculating an ozone movement vector based on the displacement values for each time point and matching the calculated ozone movement vector with the target area;
calculating an ozone flux for the target region by multiplying the ozone movement vector matched to the target region by an ozone concentration value previously matched to a pixel included in the target region; and
Among the plurality of target areas set in the step of setting the target area, target areas located within a preset radius are extracted from the monitoring area including a plurality of pixels corresponding to latitude and longitude adjacent to each other and preset according to a user's input, , summing ozone fluxes for target areas positioned corresponding to a predetermined first direction adjacent to the monitoring area among the plurality of extracted target areas, and in the first direction adjacent to the monitoring area among the plurality of extracted target areas Calculating the amount of change in ozone for the monitoring area by subtracting ozone flux for target areas located corresponding to a second preset direction different from that of the tropospheric ozone movement vector. delete delete acquiring vertical ozone data including ozone concentration values previously matched to a plurality of pixels corresponding to adjacent latitudes and longitudes, and an ozone observation time point corresponding to the ozone vertical data;
setting a target area in the ozone vertical data corresponding to the first time point;
detecting a matching region, which is a region having the highest cross-correlation coefficient with the target region, based on the ozone vertical data corresponding to each preset time point before or after the first time point;
calculating displacement values for each time point based on a position difference between the target region and the matching region; and
Calculating an ozone movement vector based on the displacement values for each time point and matching the calculated ozone movement vector with the target area;
The step of setting the target area is
A plurality of reference areas are selected by dividing the ozone vertical data corresponding to the first time point at a predetermined interval, and ozone concentrations matched to pixels located within a predetermined radius based on each of the plurality of pixels included in the reference area Calculating the standard deviation of the values, and setting the target area by changing the central position of the reference area so as to correspond to the pixel for which the largest standard deviation is calculated
A method for calculating the phosphorus tropospheric ozone transport vector. acquiring vertical ozone data including ozone concentration values previously matched to a plurality of pixels corresponding to adjacent latitudes and longitudes, and an ozone observation time point corresponding to the ozone vertical data;
setting a target area in the ozone vertical data corresponding to the first time point;
detecting a matching region, which is a region having the highest cross-correlation coefficient with the target region, based on the ozone vertical data corresponding to each preset time point before or after the first time point;
calculating displacement values for each time point based on a position difference between the target region and the matching region;
calculating an ozone movement vector based on the displacement values for each time point and matching the calculated ozone movement vector with the target area;
receiving wind field data including wind field vectors previously matched to a plurality of pixels corresponding to adjacent latitudes and longitudes;
calculating a quality factor for the ozone movement vector based on a difference between a wind field vector matched to a pixel corresponding to the target area and an ozone movement vector matched to the target area; and
and outputting an ozone movement vector matched to a target area in which a quality factor equal to or greater than a predetermined threshold value is calculated among a plurality of ozone movement vectors matched to each of the plurality of target areas. 12. The method of claim 11,
The step of detecting the matching region is
Detecting a matching region, which is a region having the highest cross-correlation coefficient with the target region, in ozone vertical data corresponding to each of the second time point before the first time point and the third time point after the first time point,
The step of calculating the displacement values for each time point is
Calculating a plurality of displacement values for each of the second time point and the third time point according to the difference between the center position of the matching area and the center position of the target area detected for each of the second time point and the third time point ,
The step of calculating the ozone movement vector and matching the calculated ozone movement vector with the target area includes:
A first vector obtained by dividing a plurality of displacement values calculated with respect to the second time point by the time between the first time point and the second time point, and a plurality of displacement values calculated for the third time point, are generating a second vector according to dividing by the time between the time point and the third time point, calculating an ozone movement vector by averaging the first vector and the second vector to match the target area;
The step of calculating the quality factor is
A first quality value according to an angle difference between the first vector and the second vector, a second quality value according to a difference in magnitudes of the absolute values of the first vector and the second vector, the first vector and the second vector A third quality value according to a vector difference between vectors and a fourth quality value according to a vector difference between the ozone movement vector and the wind field vector are calculated, respectively, and the first quality value, the second quality value, and the third quality Calculating a quality factor by multiplying and summing each of the value and the fourth quality value by a preset weight
A method for calculating the phosphorus tropospheric ozone transport vector.

Images