CN111126718A - Typhoon path prediction method - Google Patents
Typhoon path prediction method Download PDFInfo
- Publication number
- CN111126718A CN111126718A CN202010095852.9A CN202010095852A CN111126718A CN 111126718 A CN111126718 A CN 111126718A CN 202010095852 A CN202010095852 A CN 202010095852A CN 111126718 A CN111126718 A CN 111126718A
- Authority
- CN
- China
- Prior art keywords
- track
- point
- typhoon
- distance
- motion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000033001 locomotion Effects 0.000 claims abstract description 48
- 230000002596 correlated effect Effects 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000001914 filtration Methods 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 description 2
- 238000007418 data mining Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a typhoon path prediction method, which comprises the following steps: 1. importing motion path data of typhoon in the decade into a database, and taking a section of typhoon motion track A; 2. collecting the motion path data of the typhoon which is happening at present, and selecting a section of typhoon motion track B; 3. judging the correlation between the trajectories C and D of the two segments of typhoon motion, entering the step 4 if the trajectory C is correlated with the trajectory D, and returning to reselect the trajectory segment if the trajectory C is not correlated with the trajectory D; 4. and (3) extracting the overall motion trail of the typhoon belonging to the trail C, fitting the motion trail with the motion trail of the typhoon collected in the step (2), and extracting the motion trail of the typhoon in the future time period. The method can predict the trend of the typhoon path in advance, and has high accuracy.
Description
Technical Field
The invention belongs to the technical field of oceans, and particularly relates to a typhoon path prediction method.
Background
China is wide in territory, most of China is close to the sea, the sea area is wide, and many people living and working on the sea and surrounding cities are often attacked by typhoons, so that the economic development and social progress of China are hindered, and therefore, the typhoon forecasting research is very important.
Traditionally, a conceptual typhoon model based on a physical process is mostly used for describing a typhoon process, and the method is mature and can achieve a good prediction effect. However, the models are complex, the adaptability of the models for different regions is poor, and the difficulty in calibrating the model parameters is high. Therefore, data-driven typhoon process prediction methods are increasingly developed. In recent years, a relatively perfect information monitoring network is established in China, a large amount of typhoon related data is collected, the data contains the internal rules of the typhoon process, and how to improve the accuracy of hydrologic process prediction by using a data-driven prediction model established by a data mining technology is an important research direction.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the problems in the prior art, the typhoon path prediction method can predict the trend of the typhoon path in advance and has high accuracy.
The technical scheme is as follows: in order to solve the above technical problem, the present invention provides a typhoon path prediction method, which comprises the following steps:
(1) importing motion path data of typhoon in the decade into a database, and taking a section of typhoon motion track C;
(2) collecting the motion path data of the typhoon which is happening at present, and selecting a section of typhoon motion track D;
(3) judging the correlation between the tracks C and D of the two segments of typhoon motions, entering the step (4) if the track C is correlated with the track D, and returning to the step (1) if the track C is not correlated with the track D;
(4) and (3) extracting the overall motion track of the typhoon belonging to the track C, fitting the motion track with the motion track of the typhoon collected in the step (2), and extracting the motion track of the typhoon in the future time period.
Further, the specific steps of determining the correlation between the trajectories C and D of the two segments of typhoon motions in the step (3) are as follows:
(3.1) converting the longitude and latitude coordinates into a universal transverse-axis mercator projection coordinate system, and performing filtering processing on the track data by using a Kalman filtering algorithm;
(3.2) calculating the space distance between each point in the two sections of tracks C and D, and correcting;
(3.3) calculating the predicted distance between each point in the two sections of the tracks C and D;
(3.4) calculating the distance between each section of the two sections of tracks C and D;
(3.5) calculating the accumulated distance of the two tracks C and D;
(3.6) subjecting the result obtained in step (3.5) to a normalization process, i.e. normalizing the cumulative distance to [0,1], wherein 0 means that the two tracks are completely unrelated, and 1 means that the two tracks are completely similar, as long as the value is 0.5-1, there is a correlation.
Further, the specific step of calculating the spatial distance between each point in the two sections of trajectories C and D in the step (3.2) is as follows:
and selecting the ith track point in the track C and the jth track point in the track D, judging whether the distance between the point segment from the ith track point to the jth track point is equal to the distance between the point segment from the jth track point to the ith track point, and if not, correcting.
Further, the specific steps of calculating the predicted distance between each point in the two sections of trajectories C and D in the step (3.3) are as follows:
(3.3.1) taking the ith track point P in the track Ci(xi,yi) J-th track point SP in track Dj(xj,yj) (ii) a Comparison point PiAnd point SPjThe time sequence relation of (1) is that the point with large time is A and the point with small time is B; with time difference Δ t ═ tA-tB;
(3.3.2) calculating the predicted position B' of the point B, assuming that the point B is a moving point, traversing the time of each point, and searching tBAt + Δ t point B lies between which two trace points, i.e. where its predicted position is, assuming at tB+ Δ t point B is between the i-1 st and i-th points, the spatial coordinates (x) of the predicted position B' for point BB′,yB′) The calculation formula is as follows:
assuming that the motion between any two points on the track is uniform linear motion, the motion speed between the two points can be obtained, and the solving formula is as follows:
if time tB+ Δ t is absent on the B trace, then B' is calculated as follows:
wherein N is the total number of points of the track where the B point is located;
(3.3.3) calculating the predicted distance between A and B by the following formula:
distt(A,B)=dist(A,B′)
where dist (A, B ') is the Euclidean distance of A and B' in space coordinates.
Further, the specific steps of calculating the distance between each of the two sections of the trajectories C and D in the step (3.4) are as follows:
(3.4.1) wherein SiSection i of the track C, SSjIs the jth segment of the trajectory D; calculating SiAnd SSjAngle theta, S ofiTwo end points of (A) are Pi(xi,yi) And Pi+1(xi+1,yi+1),SSjAre SPj(xj,yj) And SPj+1(xj+1,yj+1) The calculation formula of the angle theta is as follows:
θ=|arctan2(yi+1-yi,xi+1-xi)-arctan2(yj+1-yj,xj+1-xj)|
(3.4.2) the spatio-temporal distance between segments is the sum of the spatio-temporal distances of the two end points of the segment, so SiAnd SSjInter-segment distance dist ofs(Si,SSj) The calculation formula of (2) is as follows:
dists(Si,SSj)=f(θ)(distst(Pi,SPj)+distst(Pi+1,SPj+1))
(3.4.3) solving f (theta), wherein the calculation formula is as follows:
wherein: ω is the negative factor of the adjustable parameter shape, with larger ω distances being less sensitive to the shape factor, where ω is 1, distsmid(Si,SSj) Is the midpoint SiAnd SSjSpace-time distance of midpoints, distmax(R, S) is the maximum spatiotemporal distance between any two points of the trajectory C and the trajectory D.
Compared with the prior art, the invention has the advantages that:
according to the method, the next step trend of the typhoon path can be predicted in advance better by improving the accuracy of the track similarity, and the accuracy is higher.
Drawings
FIG. 1 is a flow chart of the present invention.
Detailed Description
The invention is further elucidated with reference to the drawings and the detailed description. The described embodiments of the present invention are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, other embodiments obtained by a person of ordinary skill in the art without any creative effort belong to the protection scope of the present invention.
A typhoon path prediction method comprises the following steps:
(1) importing motion path data of typhoon in the decade into a database, and taking a section of typhoon motion track C;
(2) collecting the motion path data of the typhoon which is happening at present, and selecting a section of typhoon motion track D;
(3) judging the correlation between the tracks C and D of the two segments of typhoon motions, entering the step (4) if the track C is correlated with the track D, and returning to the step (1) if the track C is not correlated with the track D;
(4) and (3) extracting the overall motion track of the typhoon belonging to the track C, fitting the motion track with the motion track of the typhoon collected in the step (2), and extracting the motion track of the typhoon in the future time period.
Further, the specific steps of determining the correlation between the trajectories C and D of the two segments of typhoon motions in the step (3) are as follows:
(3.1) converting the longitude and latitude coordinates into a universal transverse-axis mercator projection coordinate system, and performing filtering processing on the track data by using a Kalman filtering algorithm;
(3.2) calculating the space distance between each point in the two sections of tracks C and D, and correcting;
(3.3) calculating the predicted distance between each point in the two sections of the tracks C and D;
(3.4) calculating the distance between each section of the two sections of tracks C and D;
(3.5) calculating the accumulated distance of the two tracks C and D;
(3.6) subjecting the result obtained in step (3.5) to a normalization process, i.e. normalizing the cumulative distance to [0,1], wherein 0 means that the two tracks are completely unrelated, and 1 means that the two tracks are completely similar, as long as the value is 0.5-1, there is a correlation.
Further, the specific step of calculating the spatial distance between each point in the two sections of trajectories C and D in the step (3.2) is as follows:
and selecting the ith track point in the track C and the jth track point in the track D, judging whether the distance between the point segment from the ith track point to the jth track point is equal to the distance between the point segment from the jth track point to the ith track point, and if not, correcting.
Further, the specific steps of calculating the predicted distance between each point in the two sections of trajectories C and D in the step (3.3) are as follows:
(3.3.1) taking the ith track point P in the track Ci(xi,yi) J-th track point SP in track Dj(xj,yj) (ii) a Comparison point PiAnd point SPjThe time sequence of (1) is defined as A, hourThe dots with small intervals are B; with time difference Δ t ═ tA-tB;
(3.3.2) calculating the predicted position B' of the point B, assuming that the point B is a moving point, traversing the time of each point, and searching tBAt + Δ t point B lies between which two trace points, i.e. where its predicted position is, assuming at tB+ Δ t point B is between the i-1 st and i-th points, the spatial coordinates (x) of the predicted position B' for point BB′,yB′) The calculation formula is as follows:
assuming that the motion between any two points on the track is uniform linear motion, the motion speed between the two points can be obtained, and the solving formula is as follows:
if time tE+ Δ t is absent on the B trace, then B' is calculated as follows:
wherein N is the total number of points of the track where the B point is located;
(3.3.3) calculating the predicted distance between A and B by the following formula:
distt(A,B)=dist(A,B′)
where dist (A, B ') is the Euclidean distance of A and B' in space coordinates.
Further, the specific steps of calculating the distance between each of the two sections of the trajectories C and D in the step (3.4) are as follows:
(3.4.1) wherein SiSection i of the track C, SSjIs the jth segment of the trajectory D; calculating SiAnd SSjAngle theta, S ofiTwo end points of (A) are Pi(xi,yi) And Pi+1(xi+1,yi+1),SSjTwo of (2)End point is SPj(xj,yj) And SPj+1(xj+1,yj+1) The calculation formula of the angle theta is as follows:
θ=|arctan2(yi+1-yi,xi+1-xi)-arctan2(yj+1-yj,xj+1-xj)|
(3.4.2) the spatio-temporal distance between segments is the sum of the spatio-temporal distances of the two end points of the segment, so SiAnd SSjInter-segment distance dist ofs(Si,SSj) The calculation formula of (2) is as follows:
dists(Si,SSj)=f(θ)(distst(Pi,SPj)+distst(Pi+1,SPj+1))
(3.4.3) solving f (theta), wherein the calculation formula is as follows:
wherein: ω is the negative factor of the adjustable parameter shape, with larger ω distances being less sensitive to the shape factor, where ω is 1, distsmid(Si,SSj) Is the midpoint SiAnd SSjSpace-time distance of midpoints, distmax(R, S) is the maximum spatiotemporal distance between any two points of the trajectory C and the trajectory D.
Claims (5)
1. A typhoon path prediction method is characterized by comprising the following steps:
(1) importing motion path data of typhoon in the decade into a database, and taking a section of typhoon motion track C;
(2) collecting the motion path data of the typhoon which is happening at present, and selecting a section of typhoon motion track D;
(3) judging the correlation between the tracks C and D of the two segments of typhoon motions, entering the step (4) if the track C is correlated with the track D, and returning to the step (1) if the track C is not correlated with the track D;
(4) and (3) extracting the overall motion track of the typhoon belonging to the track C, fitting the motion track with the motion track of the typhoon collected in the step (2), and extracting the motion track of the typhoon in the future time period.
2. The method for predicting the typhoon path according to claim 1, wherein the specific steps of judging the correlation between the trajectories C and D of the two segments of typhoon motion in the step (3) are as follows:
(3.1) converting the longitude and latitude coordinates into a universal transverse-axis mercator projection coordinate system, and performing filtering processing on the track data by using a Kalman filtering algorithm;
(3.2) calculating the space distance between each point in the two sections of tracks C and D, and correcting;
(3.3) calculating the predicted distance between each point in the two sections of the tracks C and D;
(3.4) calculating the distance between each section of the two sections of tracks C and D;
(3.5) calculating the accumulated distance of the two tracks C and D;
(3.6) subjecting the result obtained in step (3.5) to a normalization process, i.e. normalizing the cumulative distance to [0,1], wherein 0 means that the two tracks are completely unrelated, and 1 means that the two tracks are completely similar, as long as the value is 0.5-1, there is a correlation.
3. The typhoon path prediction method according to claim 1, characterized in that, the specific steps of calculating the spatial distance between each point in the two sections of trajectories C and D in the step (3.2) are as follows:
and selecting the ith track point in the track C and the jth track point in the track D, judging whether the distance between the point segment from the ith track point to the jth track point is equal to the distance between the point segment from the jth track point to the ith track point, and if not, correcting.
4. The typhoon path prediction method according to claim 1, characterized in that, the specific steps of calculating the predicted distance between each point in the two sections of trajectories C and D in the step (3.3) are as follows:
(3.3.1) taking the ith track in the tracks CPoint Pi(xi,yi) J-th track point SP in track Dj(xj,yj) (ii) a Comparison point PiAnd point SPjThe time sequence relation of (1) is that the point with large time is A and the point with small time is B; with time difference Δ t ═ tA-tB;
(3.3.2) calculating the predicted position B' of the point B, assuming that the point B is a moving point, traversing the time of each point, and searching tBAt + Δ t point B lies between which two trace points, i.e. where its predicted position is, assuming at tB+ Δ t point B is between the i-1 st and i-th points, the spatial coordinates (x) of the predicted position B' for point BB′,yB′) The calculation formula is as follows:
assuming that the motion between any two points on the track is uniform linear motion, the motion speed between the two points can be obtained, and the solving formula is as follows:
if time tB+ Δ t is absent on the B trace, then B' is calculated as follows:
wherein N is the total number of points of the track where the B point is located;
(3.3.3) calculating the predicted distance between A and B by the following formula:
distt(A,B)=dist(A,B′)
where dist (A, B ') is the Euclidean distance of A and B' in space coordinates.
5. The typhoon path prediction method according to claim 1, characterized in that, the specific steps of calculating the distance between each of the two sections of the trajectories C and D in the step (3.4) are as follows:
(3.4.1) wherein SiSection i of the track C, SSjIs the jth segment of the trajectory D; calculating SiAnd SSjAngle theta, S ofiTwo end points of (A) are Pi(xi,yi) And Pi+1(xi+1,yi+1),SSjAre SPj(xj,yj) And SPj+1(xj+1,yj+1) The calculation formula of the angle theta is as follows:
θ=|arctan2(yi+1-yi,xi+1-xi)-arctan2(yj+1-yj,xj+1-xj)|
(3.4.2) the spatio-temporal distance between segments is the sum of the spatio-temporal distances of the two end points of the segment, so SiAnd SSjInter-segment distance dist ofs(Si,SSj) The calculation formula of (2) is as follows:
dists(Si,SSj)=f(θ)(distst(Pi,SPj)+distst(Pi+1,SPj+1))
(3.4.3) solving f (theta), wherein the calculation formula is as follows:
wherein: ω is the negative factor of the adjustable parameter shape, with larger ω distances being less sensitive to the shape factor, where ω is 1, distsmid(Si,SSj) Is the midpoint SiAnd SSjSpace-time distance of midpoints, distmax(R, S) is the maximum spatiotemporal distance between any two points of the trajectory C and the trajectory D.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2019110830251 | 2019-11-07 | ||
CN201911083025 | 2019-11-07 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111126718A true CN111126718A (en) | 2020-05-08 |
Family
ID=70492110
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010095852.9A Withdrawn CN111126718A (en) | 2019-11-07 | 2020-02-17 | Typhoon path prediction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111126718A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111722306A (en) * | 2020-07-27 | 2020-09-29 | 郑州航空工业管理学院 | Typhoon landing intensity prediction method and system based on TCN network model |
CN111898815A (en) * | 2020-07-21 | 2020-11-06 | 中远海运科技(北京)有限公司 | Typhoon track prediction method and device, electronic equipment and computer readable medium |
CN116341287A (en) * | 2023-05-24 | 2023-06-27 | 自然资源部第二海洋研究所 | Track self-adaptive gridding processing method based on optimal typhoon track data |
-
2020
- 2020-02-17 CN CN202010095852.9A patent/CN111126718A/en not_active Withdrawn
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111898815A (en) * | 2020-07-21 | 2020-11-06 | 中远海运科技(北京)有限公司 | Typhoon track prediction method and device, electronic equipment and computer readable medium |
CN111722306A (en) * | 2020-07-27 | 2020-09-29 | 郑州航空工业管理学院 | Typhoon landing intensity prediction method and system based on TCN network model |
CN116341287A (en) * | 2023-05-24 | 2023-06-27 | 自然资源部第二海洋研究所 | Track self-adaptive gridding processing method based on optimal typhoon track data |
CN116341287B (en) * | 2023-05-24 | 2023-08-01 | 自然资源部第二海洋研究所 | Track self-adaptive gridding processing method based on optimal typhoon track data |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6759474B2 (en) | Vessel automatic tracking methods and systems based on deep learning networks and average shifts | |
CN111126718A (en) | Typhoon path prediction method | |
CN102567380A (en) | Method for searching vehicle information in video image | |
CN110751164B (en) | Old man travel abnormity detection method based on location service | |
CN107194408A (en) | A kind of method for tracking target of the sparse coordination model of mixed block | |
CN108986453A (en) | A kind of traffic movement prediction method based on contextual information, system and device | |
CN109543615B (en) | Double-learning-model target tracking method based on multi-level features | |
CN111008991B (en) | Background-aware related filtering target tracking method | |
CN108986143B (en) | Target detection tracking method in video | |
WO2014082480A1 (en) | Method and device for calculating number of pedestrians and crowd movement directions | |
CN104680559A (en) | Multi-view indoor pedestrian tracking method based on movement behavior mode | |
CN111242985B (en) | Video multi-pedestrian tracking method based on Markov model | |
CN107944628A (en) | A kind of accumulation mode under road network environment finds method and system | |
CN111722306A (en) | Typhoon landing intensity prediction method and system based on TCN network model | |
CN113989333A (en) | Pedestrian tracking method based on face and head and shoulder information | |
CN110531618A (en) | Closed loop based on effective key frame detects robot self-localization error cancelling method | |
CN113406659A (en) | Mobile robot position re-identification method based on laser radar information | |
CN115311617A (en) | Method and system for acquiring passenger flow information of urban rail station area | |
Ma et al. | A Method of Multiple Lane Detection Based on Constraints of Lane Information | |
CN112183221B (en) | Semantic-based dynamic object self-adaptive track prediction method | |
Chen et al. | Deep neural network based on feature fusion for indoor wireless localization | |
CN109350072B (en) | Step frequency detection method based on artificial neural network | |
CN112767438B (en) | Multi-target tracking method combining space-time motion | |
CN115880332A (en) | Target tracking method for low-altitude aircraft visual angle | |
CN115690157A (en) | Millimeter wave radar pedestrian trajectory prediction method based on Transformer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20200508 |
|
WW01 | Invention patent application withdrawn after publication |