Typhoon historical track similarity quick search and matching method

Abstract

The application provides a typhoon historical path similarity quick search matching method, comprising the following steps: step S1, obtaining historical typhoon data; step S2, typhoon generation source space search matching; step S3, establishing a typhoon path buffer zone; step S4, extracting a typhoon buffer zone intersection region; step S5, extracting a typhoon buffer intersection region center skeleton line; step S6, calculating typhoon similarity; and step S7, estimating typhoon path line density; the application can meet real-time and quick application and query, meanwhile, the buffer zone matching technology based on the typhoon moving path effectively overcomes the deficiency of the overall representation ability of the time point buffer zone, and realizes a typhoon moving path similarity quick and efficient search matching application.

Description

Technical Field

[0001] This invention belongs to the field of meteorological data processing technology, and in particular relates to a method for rapid retrieval and matching of historical typhoon paths. Background Technology

[0002] my country is one of the countries most severely affected by typhoons. Typhoons (including tropical storms, severe tropical storms, typhoons, strong typhoons, and super typhoons) are highly destructive weather systems that occur over tropical or subtropical oceans. Offshore oil and gas fields in the South my country Sea and East China Sea are significantly impacted by typhoon disasters. Statistics show that an average of 3-4 typhoons pose a major threat to offshore oil and gas platforms each year, causing significant economic losses and personnel safety hazards. Therefore, forecasting the movement paths of typhoons is particularly important for typhoon prevention and disaster mitigation for offshore oil and gas platforms. Historical typhoon similarity paths are a crucial reference factor, and typhoon path similarity retrieval and matching methods are a core technical component of this process.

[0003] Currently, similar typhoon matching algorithms primarily consider typhoon similarity indices, starting with geographical (path) similarity and further using season of occurrence, typhoon direction of movement, and center speed as similarity conditions for matching. A longer "matching duration" results in a larger matching path range, and vice versa. A series of buffer circles are drawn from the typhoon's formation point to its current position. If a typhoon's record point exists within any of these buffer circles, then that historical typhoon is considered a similar typhoon to the current real-time typhoon. However, this method, based on forecast time points, lacks sufficient ability to characterize the linear features of the overall tropical cyclone movement path. Summary of the Invention

[0004] In view of this, the present invention aims to propose a fast retrieval and matching method for the historical path similarity of typhoons, so as to overcome the shortcomings of existing methods in terms of computational complexity and path linearity representation ability.

[0005] To achieve the above objectives, the technical solution created by this invention is implemented as follows:

[0006] A fast typhoon historical path similarity retrieval and matching method includes the following steps:

[0007] Step S1: Obtain historical typhoon data;

[0008] Step S2: Spatial retrieval and matching of typhoon origin sites. The specific method is as follows: First, obtain the typhoon data to be matched. Using the geographic latitude and longitude coordinates of the current origin 3 of the typhoon to be matched as the center, set a radius of origin site to generate a circular origin site buffer zone 4. Then, use the GIS geographic point inclusion matching method to retrieve similar historical typhoons 2 that fall into the origin site buffer zone 4 from the historical typhoon data.

[0009] Step S3, establishing a typhoon path buffer zone, the specific method is: along the path of the similar historical typhoon is divided into equal distance to get the segmentation key point, with the segmentation key point as the center, set a path radius, draw a plurality of circular key point buffer zone, connect the tangent point outside each key point buffer zone, get the similar historical typhoon path buffer zone; similarly, along the path of the typhoon to be matched is divided into equal distance to get the segmentation key point, with the segmentation key point as the center, set a path radius, draw a plurality of circular key point buffer zone, connect the tangent point outside each key point buffer zone, get the typhoon to be matched path buffer zone;

[0010] Step S4, typhoon buffer zone intersection region extraction, the specific method is: the typhoon to be matched path buffer zone and the similar historical typhoon path buffer zone are intersected and cut one by one by GIS, the typhoon to be matched path buffer zone and the similar historical typhoon path buffer zone are superimposed, and the mutual overlapping common intersection region is obtained;

[0011] Step S5, typhoon buffer intersection region center skeleton line extraction, the specific method is: based on the common intersection region and geographic information image, the center skeleton line of the polygon region is calculated by using image processing technology;

[0012] Step S6, typhoon similarity calculation, the specific method is: according to the length of the center skeleton line and the length of the typhoon to be matched path, the typhoon similarity is obtained;

[0013] Step S7, typhoon path line density estimation, the specific method is: the geographic grid 12 is demarcated, the intersection curve of the typhoon path and the geographic grid 12 is obtained, and the line density probability distribution of all similar historical typhoon paths 13 is calculated by taking the typhoon similarity as a weight, and the center line of the maximum probability region is the typhoon maximum probability occurrence path.

[0014] Compared with the prior art, the typhoon historical path similarity quick retrieval matching method has the following advantages:

[0015] In order to overcome the deficiencies of the prior art method in calculation complexity and path linear representation ability, the application provides a similarity quick retrieval matching technology method based on typhoon path buffer zone intersection degree, which has smaller calculation complexity, can meet the real-time quick application and query, and effectively overcomes the deficiency of the overall representation ability of the time point buffer zone based on the buffer matching technology of the typhoon moving path, so that a typhoon moving path similarity quick and efficient retrieval matching application is realized. BRIEF DESCRIPTION OF DRAWINGS

[0016] The accompanying drawings, which form part of this invention, are used to provide a further understanding of the invention. The illustrative embodiments and descriptions of the invention are used to explain the invention and do not constitute an undue limitation of the invention. In the drawings:

[0017] Figure 1 A flowchart illustrating the method described in the embodiments of the present invention;

[0018] Figure 2 A schematic diagram of the matching and filtering process in step S2 of the present invention embodiment;

[0019] Figure 3 A schematic diagram of the typhoon path line density obtained in step S7 of the present invention.

[0020] Explanation of reference numerals in the attached figures:

[0021] 1-Similar historical typhoon origin; 2-Similar historical typhoon; 3-Current origin of the typhoon to be matched; 4-Origin buffer zone; 12-Geographic grid; 13-Similar historical typhoon path. Detailed Implementation

[0022] It should be noted that, unless otherwise specified, the embodiments and features described in the present invention can be combined with each other.

[0023] The invention will now be described in detail with reference to the accompanying drawings and embodiments.

[0024] like Figure 1 As shown, a fast typhoon historical path similarity retrieval and matching method includes the following steps:

[0025] Step S1: Obtain historical typhoon data;

[0026] Step S2: Spatial search and matching of typhoon origin locations, such as... Figure 2 As shown, the specific method is as follows: First, obtain the typhoon data to be matched, and take the geographic latitude and longitude coordinates of the current generation source 3 of the typhoon to be matched as the center, set a radius of origin, and generate a circular origin buffer zone 4; then, use the GIS geographic point inclusion matching method to retrieve similar historical typhoons 2 that fall into the origin buffer zone 4 from the historical typhoon data; further, during the retrieval, extract similar historical typhoons 2 that fall into the origin buffer zone 4 from the similar historical typhoon generation source 1.

[0027] Step S3, a typhoon path buffer zone is established, and the specific method is: the path of the similar historical typhoon is divided into segmentation key points at equal intervals, a path radius is set with the segmentation key point as the center, a plurality of circular key point buffer zones are drawn, the tangent points outside each key point buffer zone are connected to obtain the similar historical typhoon path buffer zone; similarly, the path of the to-be-matched typhoon is divided into segmentation key points at equal intervals, a path radius is set with the segmentation key point as the center, a plurality of circular key point buffer zones are drawn, the tangent points outside each key point buffer zone are connected to obtain the to-be-matched typhoon path buffer zone;

[0028] Step S4, typhoon buffer zone intersection region extraction, and the specific method is: the to-be-matched typhoon path buffer zone and the similar historical typhoon path buffer zone are subjected to GIS intersection clipping processing one by one, the to-be-matched typhoon path buffer zone and the similar historical typhoon path buffer zone are superimposed, and a common intersection region that overlaps each other is obtained;

[0029] Step S5, typhoon buffer intersection region center skeleton line extraction, and the specific method is: based on the common intersection region and geographic information image, an image processing technology is used to calculate the center skeleton line of the polygon region;

[0030] Step S6, typhoon similarity calculation, and the specific method is: according to the length of the center skeleton line and the length of the to-be-matched typhoon path, the typhoon similarity is obtained;

[0031] Step S7, typhoon path line density estimation, and the specific method is: as shown in Figure 3 , a geographic grid 12 is drawn, an intersection curve of the typhoon path and the geographic grid 12 is obtained, and the line density probability distribution of all similar historical typhoon paths 13 is calculated by taking the typhoon similarity as a weight, and the center line of the maximum probability region is the typhoon maximum probability occurrence path.

[0032] The method for calculating the center skeleton line of the polygon region in step S5 is:

[0033] First step, continuously erode to extract a pseudo skeleton, and the specific method is:

[0034] S5.1, extract the latest target contour and record the contour points;

[0035] S5.2, sequentially detect whether the 8-pixel neighborhood of the contour points contains only 3 connected pixels, if so, delete the point from the contour points and delete the corresponding point in the target image;

[0036] S5.3, sequentially detect whether the 8-pixel neighborhood of the remaining contour points in S5.2 contains only 3 or 4 connected pixels, if so, delete the point from the contour points and delete the corresponding point in the target image;

[0037] S5.4, sequentially detecting the 8-pixel neighborhood of the remaining contour points in S5.3, whether only contains 3 or 4 or 5 connected pixels, if so, delete this point from the contour points, and delete the corresponding point in the target image;

[0038] S5.5, sequentially detecting the 8-pixel neighborhood of the remaining contour points in S5.4, whether only contains 3 or 4 or 5 or 6 connected pixels, if so, delete this point from the contour points, and delete the corresponding point in the target image.

[0039] S5.6, sequentially detecting the 8-pixel neighborhood of the remaining contour points in S5.5, whether only contains 3 or 4 or 5 or 6 or 7 connected pixels, if so, delete this point from the contour points, and delete the corresponding point in the target image;

[0040] Second step, extracting the real skeleton from the pseudo skeleton, the pseudo skeleton has some areas of two-pixel width, while the target skeleton is single-pixel width; sequentially detecting the 8-pixel neighborhood of the pseudo skeleton in the target image, whether only contains 2 or 3 or 4 or 5 or 6 or 7 connected pixels, if so, delete this point from the pseudo skeleton, and get the center skeleton line.

[0041] In the step S6, the calculation formula of the typhoon similarity is:

[0042]

[0043] Wherein, L 中心骨架线 is the length of the center skeleton line, L 待匹配台风路径 is the length of the typhoon path to be matched, and AD is the typhoon similarity.

[0044] In the step S7, the calculation formula of the line density is:

[0045]

[0046] Wherein, Li represents the internal curve of each typhoon path intersecting with the geographical grid, i.e. the curve length falling into the geographical grid; Vi represents the similarity of each historical typhoon path falling into the geographical grid and the path to be matched, as the influence weight of the path pair line density; the line density LD is calculated by multiplying the length of each typhoon path falling into the geographical grid by the similarity weight coefficient and then dividing by the area of the geographical grid.

[0047] A typhoon historical path similarity fast retrieval matching method further comprises the following steps:

[0048] Step S8, typhoon path trend change analysis, such as Figure 3As shown, the specific method is: the larger the value of the line density of the geographical grid, the more typhoons in the historical similar typhoon path pass through the geographical grid at the same time, that is, the higher the probability of occurrence, according to the historical statistical characteristics, the place where the historical probability appears the highest is also the largest area through which the future typhoon path passes, which can be used as the largest possible path for typhoon path trend prediction.

[0049] The application provides a tropical cyclone historical path similarity quick retrieval matching technical method, mainly including the following steps: step S1, obtaining historical typhoon data; step S2, typhoon generation source space retrieval matching; step S3, establishing a typhoon path buffer area; step S4, extracting a typhoon buffer intersection area; step S5, extracting a typhoon buffer intersection area center skeleton line; step S6, calculating a typhoon similarity; and step S7, estimating a typhoon path line density. The method has a small calculation complexity, can meet real-time quick application and query, and effectively overcomes the deficiency of the overall representation ability of the time point buffer area based on the buffer area matching technology of the typhoon moving path, realizes a typhoon moving path similarity quick and efficient retrieval matching application, and overcomes the deficiency of the prior art method in calculation complexity and path linear representation ability.

[0050] The above description is only the preferred embodiment of the application, and is not intended to limit the application. Any modification, equivalent replacement, improvement, etc. within the spirit and principle of the application shall be included in the protection scope of the application.

Claims

1. A fast retrieval and matching method for historical typhoon paths, characterized in that, Includes the following steps: Step S1: Obtain historical typhoon data; Step S2: Spatial retrieval and matching of typhoon origin. The specific method is as follows: First, obtain the typhoon data to be matched. Take the geographic latitude and longitude coordinates of the current origin (3) of the typhoon to be matched as the center and set a radius of origin to generate a circular origin buffer (4). Then, use the GIS geographic point inclusion matching method to retrieve similar historical typhoons (2) that fall into the origin buffer (4) in the historical typhoon data. Step S3: Establish a typhoon path buffer zone. The specific method is as follows: Divide the path of the similar historical typhoon at equal distances to obtain segmentation key points. Using the segmentation key points as the center, set a path radius and draw multiple circular key point buffer zones. Connect the tangent points on the outside of each key point buffer zone to obtain the similar historical typhoon path buffer zone. Similarly, divide the path of the typhoon to be matched at equal distances to obtain segmentation key points. Using the segmentation key points as the center, set a path radius and draw multiple circular key point buffer zones. Connect the tangent points on the outside of each key point buffer zone to obtain the typhoon path buffer zone to be matched. Step S4: Extracting the intersection area of ​​the typhoon buffer zone. The specific method is as follows: Perform GIS intersection clipping processing on the typhoon path buffer zone to be matched and the similar historical typhoon path buffer zone one by one. Overlay the typhoon path buffer zone to be matched and the similar historical typhoon path buffer zone to obtain the common intersection area that overlaps with each other. Step S5: Extraction of the center skeleton line of the typhoon buffer intersection area. The specific method is as follows: Based on the common intersection area and the geographic information image, the center skeleton line of the polygon area is calculated using image processing technology. Step S6: Typhoon similarity calculation. The specific method is as follows: the typhoon similarity is obtained based on the length of the central skeleton line and the length of the path of the typhoon to be matched. In step S6, the formula for calculating typhoon similarity is: Among them, L 中心骨架线 L is the length of the central skeleton line. 待匹配台风路径 is the length of the typhoon path to be matched, and AD is the typhoon similarity score. Step S7: Typhoon path line density estimation. The specific method is as follows: delineate a geographic grid (12), obtain the intersection curve between the typhoon path and the geographic grid (12), and use the typhoon similarity as the weight to calculate the line density probability distribution of all similar historical typhoon paths (13). The center line of the region with the highest probability is the path of the typhoon with the highest probability.

2. The method for fast retrieval and matching of historical typhoon paths according to claim 1, characterized in that, The method for calculating the center skeleton line of the polygonal region in step S5 is as follows: The first step is to continuously erode the pseudo-skeleton, and the specific method is as follows: S5.1 Extract the latest target contour and record these contour points; S5.

2. Sequentially detect whether the 8-pixel neighborhood of these contour points contains only 3 connected pixels. If so, delete this point from the contour points and delete the corresponding point in the target image. S5.

3. Sequentially check the 8-pixel neighborhood of the remaining contour points in S5.2 to see if they contain only 3 or 4 connected pixels. If so, delete this point from the contour points and delete the corresponding point in the target image. S5.

4. Sequentially check the 8-pixel neighborhood of the remaining contour points in S5.3 to see if they contain only 3, 4, or 5 connected pixels. If so, delete this point from the contour points and delete the corresponding point in the target image. S5.

5. Sequentially check the 8-pixel neighborhood of the remaining contour points in S5.4 to see if they contain only 3, 4, 5, or 6 connected pixels. If so, delete this point from the contour points and delete the corresponding point in the target image. S5.

6. Sequentially check the 8-pixel neighborhood of the remaining contour points in S5.5 to see if they contain only 3 or 4 or 5 or 6 or 7 connected pixels. If so, delete this point from the contour points and delete the corresponding point in the target image. The second step is to extract the real skeleton from the pseudo skeleton. Some areas of the pseudo skeleton are two pixels wide, while the target skeleton is a single-layer pixel wide. The pseudo skeleton is then checked in the target image to see if its 8-pixel neighborhood contains only 2, 3, 4, 5, 6, or 7 connected pixels. If so, this point is removed from the pseudo skeleton to obtain the central skeleton line.

3. The method for fast retrieval and matching of historical typhoon paths according to claim 1, characterized in that: In step S7, the formula for calculating the linear density is: Wherein, Li represents the internal curve where each typhoon path intersects with the geographic grid (12), that is, the length of the curve that falls into the geographic grid (12), Vi represents the similarity between each historical typhoon path falling into the geographic grid (12) and the path to be matched, which is used as the influence weight of the path on the line density; the line density LD is calculated by multiplying the length of each typhoon path falling into the geographic grid (12) by the similarity weight coefficient and summing the results, and then dividing by the area of ​​the geographic grid.

4. The method for rapid retrieval and matching of historical typhoon paths according to any one of claims 1 to 3, characterized in that: It also includes the following steps: Step S8, Typhoon path trend change analysis, the specific method is: the larger the value of the line density of the geographic grid, the more typhoons that pass through the geographic grid (12) in the historical similar typhoon paths, that is, the higher the probability of occurrence. According to the historical statistical characteristics, the place with the highest historical probability is the largest area that the future typhoon path will pass through, which is the maximum possible path for typhoon path trend prediction.

Images