CN104537891A - Ship trajectory real-time predicting method - Google Patents
Ship trajectory real-time predicting method Download PDFInfo
- Publication number
- CN104537891A CN104537891A CN201410841564.8A CN201410841564A CN104537891A CN 104537891 A CN104537891 A CN 104537891A CN 201410841564 A CN201410841564 A CN 201410841564A CN 104537891 A CN104537891 A CN 104537891A
- Authority
- CN
- China
- Prior art keywords
- boats
- ships
- prime
- time
- real
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G3/00—Traffic control systems for marine craft
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G3/00—Traffic control systems for marine craft
- G08G3/02—Anti-collision systems
Abstract
The invention relates to a ship trajectory real-time predicting method. The ship trajectory real-time predicting method comprises the following steps that firstly, real-time and historical position information of a ship is obtained through sea surface radar and is primarily processed; then, ship trajectory data are preprocessed at each sampling time, then are clustered at each sampling time and are subjected to parameter training through a hidden Markov model at each sampling time, and then a hidden state q corresponding to the observation value at the current time is obtained by the adoption of a Viterbi algorithm according to parameters of the hidden Markov model at each sampling time; finally, by setting the prediction time domain W at each sampling time, on the basis of the hidden state q of the ship at the current time, the position prediction value O of the ship in the future is obtained, and therefore the ship trajectory in the future can be obtained in a rolling speculation mode at each sampling time. The ship trajectory real-time predicting method predicts the ship trajectory in real time in a rolling mode, is good in accuracy and provides powerful guarantee for subsequent ship conflict resolution.
Description
Technical field
The present invention relates to a kind of marine site traffic control method, particularly relate to a kind of boats and ships track real-time predicting method based on Rolling Planning strategy.
Background technology
Along with the fast development of global shipping business, the traffic in the busy marine site of part is further crowded.In the intensive complicated marine site of vessel traffic flow, still adopt sail plan can not adapt to the fast development of shipping business in conjunction with the regulation model that artificial interval is allocated for the collision scenario between boats and ships.For ensureing the personal distance between boats and ships, implement the emphasis that effective conflict allotment just becomes marine site traffic control work.Boats and ships conflict Resolution is a gordian technique in navigational field, frees scheme safely and efficiently for increasing marine site boats and ships flow and guaranteeing that sea-freight safety is significant.
In order to improve the efficiency of navigation of boats and ships, marine radar automatic plotter has been widely applied in ship monitor and collision prevention at present, and this equipment provides reference frame by extracting boats and ships relevant informations for the judgement of collision scenario between boats and ships.Although this kind equipment greatly reduces manual supervisory load, it does not have the automatic conflict Resolution function of boats and ships.And boats and ships conflict Resolution is on the basis based on the prediction to boats and ships track, in boats and ships real navigation, by the impact of the various factors such as meteorological condition, navigator and driver's operation, its running status often not exclusively belongs to a certain specific motion state, in boats and ships trajectory predictions process, need the impact considering various enchancement factor, by the up-to-date characteristic that obtains all kinds of enchancement factor rolling forecast implemented to its Future Trajectory and strengthen the robustness of its trajectory predictions.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of robustness good boats and ships track real-time predicting method, and the boats and ships trajectory predictions precision of the method is higher.
The technical scheme realizing the object of the invention is to provide a kind of boats and ships track real-time predicting method, comprises following several step:
1. obtain the real-time of boats and ships and historical position information by sea radar, the positional information of each boats and ships is discrete two-dimensional position sequence x '=[x
1', x
2' ..., x
n'] and y '=[y
1', y
2' ..., y
n'], by application wavelet transformation theory to original discrete two-dimensional position sequence x '=[x
1', x
2' ..., x
n'] and y '=[y
1', y
2' ..., y
n'] carry out rough handling, thus obtain the denoising discrete two-dimensional position sequence x=[x of boats and ships
1, x
2..., x
n] and y=[y
1, y
2..., y
n];
2. in each sampling instant to the pre-service of boats and ships track data, according to the boats and ships original discrete two-dimensional position sequence x=[x that obtains
1, x
2..., x
n] and y=[y
1, y
2..., y
n], adopt first order difference method to carry out processing new boats and ships discrete location sequence Δ x=[the Δ x of acquisition to it
1, Δ x
2..., Δ x
n-1] and Δ y=[Δ y
1, Δ y
2..., Δ y
n-1], wherein Δ x
i=x
i+1-x
i, Δ y
i=y
i+1-y
i(i=1,2 ..., n-1);
3. in each sampling instant to boats and ships track data cluster, to new boats and ships discrete two-dimensional position sequence Δ x and Δ y after process, by setting cluster number M ', adopt K-means clustering algorithm to carry out cluster to it respectively;
4. Hidden Markov Model (HMM) is utilized to carry out parameter training in each sampling instant to boats and ships track data, by the vessel motion track data Δ x after process and Δ y being considered as the aobvious observed reading of hidden Markov models, upgrade period τ ' by setting hidden state number N and parameter, according to nearest T ' individual position detection value and adopt B-W algorithm to roll to obtain up-to-date Hidden Markov Model (HMM) parameter lambda ';
5. in each sampling instant foundation Hidden Markov Model (HMM) parameter, the hidden state q corresponding to Viterbi algorithm acquisition current time observed reading is adopted;
6. in each sampling instant, by setting prediction time domain W, based on the hidden state q of boats and ships current time, obtain the position prediction value O of future time period boats and ships, thus infer the track to boats and ships in future time period in each sampling instant rolling.
Further, described step 1. in, by application wavelet transformation theory to original discrete two-dimensional position sequence x '=[x
1', x
2' ..., x
n'] and y '=[y
1', y
2' ..., y
n'] carry out rough handling, thus obtain the denoising discrete two-dimensional position sequence x=[x of boats and ships
1, x
2..., x
n] and y=[y
1, y
2..., y
n]: for given original two dimensional sequence data x '=[x
1', x
2' ..., x
n'], utilize the linear representation of following form to be similar to it respectively:
Wherein:
F ' (x ') represents the function expression obtained after data smoothing processing, and ψ (x ') represents female ripple, and δ, J and K are wavelet transformation constant, ψ
j, K(x ') represents the transition form of female ripple, c
j, Krepresent the function coefficients obtained by wavelet transform procedure, it embodies wavelet ψ
j, K(x '), to the weight size of whole approximation to function, if this coefficient is very little, so it means wavelet ψ
j, Kthe weight of (x ') is also less, thus can under the prerequisite of not influence function key property, by wavelet ψ from approximation to function process
j, K(x ') removes; In real data processing procedure, implemented " threshold transition " by setting threshold value χ, work as c
j, Kduring < χ, setting c
j, K=0; Choosing of threshold function table adopts the following two kinds mode:
For y '=[y
1', y
2' ..., y
n'], also adopt said method to carry out denoising.
Further, described step 4. in determine that flight path Hidden Markov Model (HMM) parameter lambda '=process of (π, A, B) is as follows:
4.1) variable initialize: application is uniformly distributed to variable π
i, a
ijand b
j(o
k) initialize
with
and make it meet constraint condition:
With
Obtain λ thus
0=(π
0, A
0, B
0), wherein o
krepresent a certain aobvious observed reading, π
0, A
0and B
0by element respectively
with
the matrix formed, makes parameter l=0, o=(o
t-T '+1..., o
t-1, o
t) be the individual historical position observed reading of the T ' before current time t;
4.2) E-M algorithm is performed:
4.2.1) E-step: by λ
1calculate ξ
e(i, j) and γ
e(s
i);
Variable
So
Wherein s represents a certain hidden state;
4.2.2) M-step: use
Estimate π respectively
i, a
ijand b
j(o
k) and obtain λ thus
l+1;
4.2.3) circulate: l=l+1, repeats E-step and M-step, until π
i, a
ijand b
j(o
k) convergence, namely
wherein parameter ε=0.00001, returns step 4.2.4);
4.2.4): make λ '=λ
l+1, algorithm terminates.
Further, described step 5. in determine that the iterative process of the best hidden status switch of ship track is as follows:
5.1) variable initialize: make g=2, β
t '(s
i)=1 (s
i∈ S), δ
1(s
i)=π
ib
i(o
1), ψ
1(s
i)=0, wherein,
, wherein variable ψ
g(s
j) represent make variable δ
g-1(s
i) a
ijget the hidden state s of ship track of maximal value
i, parameter S represents the set of hidden state;
5.2) recursive process:
5.3) moment upgrade: make g=g+1, if g≤T ', return step 5.2), otherwise iteration ends and forward step 5.4 to);
5.4)
Forward step 5.5 to);
5.5) optimum hidden status switch obtains:
5.5.1) variable initialize: make g=T '-1;
5.5.2) backward recursion:
5.5.3) moment upgrades: make g=g-1, if g >=1, return step 5.5.2), otherwise stop.
Further, described step 3. in, the value of cluster number M ' is 4.
Further, described step 4. in, the value of state number N is 3, parameter upgrade period τ ' be 30 seconds, T ' is 10.
Further, described step 6. in, prediction time domain W be 300 seconds.
The present invention has positive effect: (1) the present invention is in the process of boats and ships track real-time estimate, incorporate the impact of enchancement factor, the rolling track prediction scheme adopted can extract the changing condition of extraneous enchancement factor in time, improves the accuracy of boats and ships trajectory predictions.
(2) the present invention is based on different performance index, its boats and ships track real-time estimate result can provide for the multiple boats and ships that there is conflict frees trajectory planning scheme, improves the economy of vessel motion and the utilization factor of sea area resources.
Accompanying drawing explanation
Fig. 1 is the vessel motion short-term Track Pick-up schematic flow sheet in the present invention.
Embodiment
(embodiment 1)
See Fig. 1, a kind of boats and ships track real-time predicting method of the present embodiment comprises following several step:
1. obtain the real-time of boats and ships and historical position information by sea radar, the positional information of each boats and ships is discrete two-dimensional position sequence x '=[x
1', x
2' ..., x
n'] and y '=[y
1', y
2' ..., y
n'], by application wavelet transformation theory to original discrete two-dimensional position sequence x '=[x
1', x
2' ..., x
n'] and y '=[y
1', y
2' ..., y
n'] carry out rough handling, thus obtain the denoising discrete two-dimensional position sequence x=[x of boats and ships
1, x
2..., x
n] and y=[y
1, y
2..., y
n]: y=[y
1, y
2..., y
n]: for given original two dimensional sequence data x '=[x
1', x
2' ..., x
n'], utilize the linear representation of following form to be similar to it respectively:
Wherein:
F ' (x ') represents the function expression obtained after data smoothing processing, and ψ (x ') represents female ripple, and δ, J and K are wavelet transformation constant, ψ
j, K(x ') represents the transition form of female ripple, c
j, Krepresent the function coefficients obtained by wavelet transform procedure, it embodies wavelet ψ
j, K(x '), to the weight size of whole approximation to function, if this coefficient is very little, so it means wavelet ψ
j, Kthe weight of (x ') is also less, thus can under the prerequisite of not influence function key property, by wavelet ψ from approximation to function process
j, K(x ') removes; In real data processing procedure, implemented " threshold transition " by setting threshold value χ, work as c
j, Kduring < χ, setting c
j, K=0; Choosing of threshold function table adopts the following two kinds mode:
For y '=[y
1', y
2' ..., y
n'], also adopt said method to carry out denoising;
2. in each sampling instant to the pre-service of boats and ships track data, according to the boats and ships original discrete two-dimensional position sequence x=[x that obtains
1, x
2..., x
n] and y=[y
1, y
2..., y
n], adopt first order difference method to carry out processing new boats and ships discrete location sequence Δ x=[the Δ x of acquisition to it
1, Δ x
2..., Δ x
n-1] and Δ y=[Δ y
1, Δ y
2..., Δ y
n-1], wherein Δ x
i=x
i+1-x
i, Δ y
i=y
i+1-y
i(i=1,2 ..., n-1);
3. in each sampling instant to boats and ships track data cluster, to new boats and ships discrete two-dimensional position sequence Δ x and Δ y after process, by setting cluster number M ', adopt K-means clustering algorithm to carry out cluster to it respectively;
4. Hidden Markov Model (HMM) is utilized to carry out parameter training in each sampling instant to boats and ships track data, by the vessel motion track data Δ x after process and Δ y being considered as the aobvious observed reading of hidden Markov models, upgrade period τ ' by setting hidden state number N and parameter, according to nearest T ' individual position detection value and adopt B-W algorithm to roll to obtain up-to-date Hidden Markov Model (HMM) parameter lambda '; Determine that flight path Hidden Markov Model (HMM) parameter lambda '=process of (π, A, B) is as follows:
4.1) variable initialize: application is uniformly distributed to variable π
i, a
ijand b
j(o
k) initialize
with
and make it meet constraint condition:
With
Obtain λ thus
0=(π
0, A
0, B
0), wherein o
krepresent a certain aobvious observed reading, π
0, A
0and B
0by element respectively
with
the matrix formed, makes parameter l=0, o=(o
t-T '+1..., o
t-1, o
t) be the individual historical position observed reading of the T ' before current time t;
4.2) E-M algorithm is performed:
4.2.1) E-step: by λ
1calculate ξ
e(i, j) and γ
e(s
i);
Variable
So
Wherein s represents a certain hidden state;
4.2.2) M-step: use
Estimate π respectively
i, a
ijand b
j(o
k) and obtain λ thus
l+1;
4.2.3) circulate: l=l+1, repeats E-step and M-step, until π
i, a
ijand b
j(o
k) convergence, namely
wherein parameter ε=0.00001, returns step 4.2.4);
4.2.4): make λ '=λ
l+1, algorithm terminates.
5. in each sampling instant foundation Hidden Markov Model (HMM) parameter, the hidden state q corresponding to Viterbi algorithm acquisition current time observed reading is adopted:
5.1) variable initialize: make g=2, β
t '(s
i)=1 (s
i∈ S), δ
1(s
i)=π
ib
i(o
1), ψ
1(s
i)=0, wherein,
, wherein variable ψ
g(s
j) represent make variable δ
g-1(s
i) a
ijget the hidden state s of ship track of maximal value
i, parameter S represents the set of hidden state;
5.2) recursive process:
5.3) moment upgrade: make g=g+1, if g≤T ', return step 5.2), otherwise iteration ends and forward step 5.4 to);
5.4)
Forward step 5.5 to);
5.5) optimum hidden status switch obtains:
5.5.1) variable initialize: make g=T '-1;
5.5.2) backward recursion:
5.5.3) moment upgrades: make g=g-1, if g >=1, return step 5.5.2), otherwise stop.。
6. in each sampling instant, by setting prediction time domain W, based on the hidden state q of boats and ships current time, the position prediction value O of future time period boats and ships is obtained
The value of above-mentioned cluster number M ' is 4, the value of state number N is 3, and it is 30 seconds that parameter upgrades period τ ', and T ' is 10, and prediction time domain W is 300 seconds.
(application examples, navigation traffic control method)
The navigation traffic control method of the present embodiment comprises following several step:
Steps A, the boats and ships track real-time predicting method obtained according to embodiment 1 obtain the track of boats and ships boats and ships in the future time period that each sampling instant is inferred;
Step B, in each sampling instant, the running status current based on boats and ships and historical position observation sequence, obtain the numerical value of marine site wind field variable, its detailed process is as follows:
B.1) stop position setting boats and ships is track reference coordinate initial point;
B.2) when boats and ships are in straight running condition and at the uniform velocity turning running status, marine site wind field linear filtering model is built;
B.3) numerical value of wind field variable is obtained according to constructed Filtering Model.
Step C, in each sampling instant, the safety rule collection that need meet when running in marine site based on the running status of each boats and ships and the boats and ships of setting, when likely occurring violating the situation of safety rule when between boats and ships, provide warning information timely to its dynamic behaviour implementing monitoring and for maritime traffic control center;
Step D, when warning information occurs, under the prerequisite meeting boats and ships physical property and marine site traffic rules, by setting optimizing index function and incorporating wind field variable value, Model Predictive Control Theory method is adopted to carry out Rolling Planning to boats and ships collision avoidance track, and program results is transferred to the execution of each boats and ships, its detailed process is as follows:
D.1) termination reference point locations P, collision avoidance policy control time domain Θ, the trajectory predictions time domain γ of boats and ships collision avoidance trajectory planning is set;
D.2) under being set in the prerequisite of given optimizing index function, based on cooperative collision avoidance trajectory planning thought, give different weights by giving each boats and ships and incorporate real-time wind field variable filtering numerical value, obtain the collision avoidance track of each boats and ships and collision avoidance control strategy and program results is transferred to each boats and ships performing, and its first Optimal Control Strategy only implemented by each boats and ships in Rolling Planning interval;
D.3) in next sampling instant, step 5.2 is repeated) until each boats and ships all arrive it free terminal.
Above-mentioned termination reference point locations P is set as the next navigation channel point of vessel position conflict point, and collision avoidance policy control time domain Θ is 300 seconds; Trajectory predictions time domain γ is 300 seconds.
Obviously, above-described embodiment is only for example of the present invention is clearly described, and is not the restriction to embodiments of the present invention.For those of ordinary skill in the field, can also make other changes in different forms on the basis of the above description.Here exhaustive without the need to also giving all embodiments.And these belong to spirit institute's apparent change of extending out of the present invention or change and are still among protection scope of the present invention.
Claims (7)
1. a boats and ships track real-time predicting method, is characterized in that comprising following several step:
1. obtain the real-time of boats and ships and historical position information by sea radar, the positional information of each boats and ships is discrete two-dimensional position sequence x '=[x
1', x
2' ..., x
n'] and y '=[y
1', y
2' ..., y
n'], by application wavelet transformation theory to original discrete two-dimensional position sequence x '=[x
1', x
2' ..., x
n'] and y '=[y
1', y
2' ..., y
n'] carry out rough handling, thus obtain the denoising discrete two-dimensional position sequence x=[x of boats and ships
1, x
2..., x
n] and y=[y
1, y
2..., y
n];
2. in each sampling instant to the pre-service of boats and ships track data, according to the boats and ships original discrete two-dimensional position sequence x=[x that obtains
1, x
2..., x
n] and y=[y
1, y
2..., y
n], adopt first order difference method to carry out processing new boats and ships discrete location sequence Δ x=[the Δ x of acquisition to it
1, Δ x
2..., Δ x
n-1] and Δ y=[Δ y
1, Δ y
2..., Δ y
n-1], wherein Δ x
i=x
i+1-x
i, Δ y
i=y
i+1-y
i(i=1,2 ..., n-1);
3. in each sampling instant to boats and ships track data cluster, to new boats and ships discrete two-dimensional position sequence Δ x and Δ y after process, by setting cluster number M ', adopt K-means clustering algorithm to carry out cluster to it respectively;
4. Hidden Markov Model (HMM) is utilized to carry out parameter training in each sampling instant to boats and ships track data, by the vessel motion track data Δ x after process and Δ y being considered as the aobvious observed reading of hidden Markov models, upgrade period τ ' by setting hidden state number N and parameter, according to nearest T ' individual position detection value and adopt B-W algorithm to roll to obtain up-to-date Hidden Markov Model (HMM) parameter lambda ';
5. in each sampling instant foundation Hidden Markov Model (HMM) parameter, the hidden state q corresponding to Viterbi algorithm acquisition current time observed reading is adopted;
6. in each sampling instant, by setting prediction time domain W, based on the hidden state q of boats and ships current time, obtain the position prediction value O of future time period boats and ships, thus infer the track to boats and ships in future time period in each sampling instant rolling.
2. a kind of boats and ships track real-time predicting method according to claim 1, is characterized in that: described step 1. in, by application wavelet transformation theory to original discrete two-dimensional position sequence x '=[x
1', x
2' ..., x
n'] and y '=[y
1', y
2' ..., y
n'] carry out rough handling, thus obtain the denoising discrete two-dimensional position sequence x=[x of boats and ships
1, x
2..., x
n] and y=[y
1, y
2..., y
n]: for given original two dimensional sequence data x '=[x
1', x
2' ..., x
n'], utilize the linear representation of following form to be similar to it respectively:
Wherein:
ψ
J,K(x′)=δ·ψ(2
Jx′-K)
F ' (x ') represents the function expression obtained after data smoothing processing, and ψ (x ') represents female ripple, and δ, J and K are wavelet transformation constant, ψ
j, K(x ') represents the transition form of female ripple, c
j, Krepresent the function coefficients obtained by wavelet transform procedure, it embodies wavelet ψ
j, K(x '), to the weight size of whole approximation to function, if this coefficient is very little, so it means wavelet ψ
j, Kthe weight of (x ') is also less, thus can under the prerequisite of not influence function key property, by wavelet ψ from approximation to function process
j, K(x ') removes; In real data processing procedure, implemented " threshold transition " by setting threshold value χ, work as c
j, Kduring < χ, setting c
j, K=0; Choosing of threshold function table adopts the following two kinds mode:
For y '=[y
1', y
2' ..., y
n'], also adopt said method to carry out denoising.
3. a kind of boats and ships track real-time predicting method according to claim 1 and 2, is characterized in that: described step 4. in determine that flight path Hidden Markov Model (HMM) parameter lambda '=process of (π, A, B) is as follows:
4.1) variable initialize: application is uniformly distributed to variable π
i, a
ijand b
j(o
k) initialize
with
and make it meet constraint condition:
With
Obtain λ thus
0=(π
0, A
0, B
0), wherein o
krepresent a certain aobvious observed reading, π
0, A
0and B
0by element respectively
with
the matrix formed, makes parameter l=0, o=(o
t-T '+1..., o
t-1, o
t) be the individual historical position observed reading of the T ' before current time t;
4.2) E-M algorithm is performed:
4.2.1) E-step: by λ
lcalculate ξ
e(i, j) and γ
e(s
i);
Variable
So
Wherein s represents a certain hidden state;
4.2.2) M-step: use
Estimate π respectively
i, a
ijand b
j(o
k) and obtain λ thus
l+1;
4.2.3) circulate: l=l+1, repeats E-step and M-step, until π
i, a
ijand b
j(o
k) convergence, namely | P (o| λ
l+1)-P (o| λ
l) | < ε, wherein parameter ε=0.00001, return step 4.2.4);
4.2.4): make λ '=λ
l+1, algorithm terminates.
4., according to a kind of boats and ships track real-time predicting method one of claims 1 to 3 Suo Shu, it is characterized in that: described step 5. in determine that the iterative process of the best hidden status switch of ship track is as follows:
5.1) variable initialize: make g=2, β
t '(s
i)=1 (s
i∈ S), δ
1(s
i)=π
ib
i(o
1), ψ
1(s
i)=0, wherein,
Wherein variable ψ
g(s
j) represent make variable δ
g-1(s
i) a
ijget the hidden state s of ship track of maximal value
i, parameter S represents the set of hidden state;
5.2) recursive process:
5.3) moment upgrade: make g=g+1, if g≤T ', return step 5.2), otherwise iteration ends and forward step 5.4 to);
5.4)
Forward step 5.5 to);
5.5) optimum hidden status switch obtains:
5.5.1) variable initialize: make g=T '-1;
5.5.2) backward recursion:
5.5.3) moment upgrades: make g=g-1, if g >=1, return step 5.5.2), otherwise stop.
5., according to a kind of boats and ships track real-time predicting method one of Claims 1-4 Suo Shu, it is characterized in that: described step 3. in, the value of cluster number M ' is 4.
6., according to a kind of boats and ships track real-time predicting method one of claim 1 to 5 Suo Shu, it is characterized in that: described step 4. in, the value of state number N is 3, parameter upgrade period τ ' be 30 seconds, T ' is 10.
7., according to a kind of boats and ships track real-time predicting method one of claim 1 to 6 Suo Shu, it is characterized in that: described step 6. in, prediction time domain W be 300 seconds.
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610624893.6A CN106228850A (en) | 2014-12-30 | 2014-12-30 | Boats and ships track real-time predicting method based on Rolling Planning strategy |
CN201410841564.8A CN104537891B (en) | 2014-12-30 | 2014-12-30 | A kind of boats and ships track real-time predicting method |
CN201610620526.9A CN106251704A (en) | 2014-12-30 | 2014-12-30 | Boats and ships trajectory predictions method |
CN201610620437.4A CN106205213A (en) | 2014-12-30 | 2014-12-30 | A kind of boats and ships trajectory predictions method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410841564.8A CN104537891B (en) | 2014-12-30 | 2014-12-30 | A kind of boats and ships track real-time predicting method |
Related Child Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610620526.9A Division CN106251704A (en) | 2014-12-30 | 2014-12-30 | Boats and ships trajectory predictions method |
CN201610620437.4A Division CN106205213A (en) | 2014-12-30 | 2014-12-30 | A kind of boats and ships trajectory predictions method |
CN201610624893.6A Division CN106228850A (en) | 2014-12-30 | 2014-12-30 | Boats and ships track real-time predicting method based on Rolling Planning strategy |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104537891A true CN104537891A (en) | 2015-04-22 |
CN104537891B CN104537891B (en) | 2016-08-24 |
Family
ID=52853410
Family Applications (4)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610620526.9A Pending CN106251704A (en) | 2014-12-30 | 2014-12-30 | Boats and ships trajectory predictions method |
CN201610624893.6A Pending CN106228850A (en) | 2014-12-30 | 2014-12-30 | Boats and ships track real-time predicting method based on Rolling Planning strategy |
CN201410841564.8A Active CN104537891B (en) | 2014-12-30 | 2014-12-30 | A kind of boats and ships track real-time predicting method |
CN201610620437.4A Pending CN106205213A (en) | 2014-12-30 | 2014-12-30 | A kind of boats and ships trajectory predictions method |
Family Applications Before (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610620526.9A Pending CN106251704A (en) | 2014-12-30 | 2014-12-30 | Boats and ships trajectory predictions method |
CN201610624893.6A Pending CN106228850A (en) | 2014-12-30 | 2014-12-30 | Boats and ships track real-time predicting method based on Rolling Planning strategy |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610620437.4A Pending CN106205213A (en) | 2014-12-30 | 2014-12-30 | A kind of boats and ships trajectory predictions method |
Country Status (1)
Country | Link |
---|---|
CN (4) | CN106251704A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106154296A (en) * | 2015-06-26 | 2016-11-23 | 安徽华米信息科技有限公司 | The method of adjustment of a kind of path locus and device |
CN106251704A (en) * | 2014-12-30 | 2016-12-21 | 江苏理工学院 | Boats and ships trajectory predictions method |
CN106595665A (en) * | 2016-11-30 | 2017-04-26 | 耿生玲 | Prediction method for spatial-temporal trajectory of moving object in obstructed space |
CN112070312A (en) * | 2020-09-10 | 2020-12-11 | 龙马智芯(珠海横琴)科技有限公司 | Flight path prediction method, flight path prediction device, and computer-readable storage medium |
CN113627104A (en) * | 2021-08-12 | 2021-11-09 | 北京中安智能信息科技有限公司 | Underwater submarine track simulation method, device and equipment under multi-constraint condition |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10780968B2 (en) * | 2017-11-27 | 2020-09-22 | Applied Physical Sciences Corp. | Systems and methods for wave sensing and ship motion forecasting with operational period indicators |
CN110009936B (en) * | 2019-03-15 | 2021-01-05 | 北京海兰信数据科技股份有限公司 | Ship auxiliary collision prevention method for busy water area |
CN110706518B (en) * | 2019-08-19 | 2022-04-22 | 浙江海洋大学 | Ship port-hanging and berthing prediction technology |
CN111929715A (en) * | 2020-06-28 | 2020-11-13 | 杭州云起智慧校园科技有限公司 | Positioning method, device and equipment for school badge and storage medium |
CN113093735B (en) * | 2021-03-18 | 2022-07-26 | 华南理工大学 | Cross-section function control method, system and medium of under-actuated unmanned ship with limited output under random disturbance |
CN113221449B (en) * | 2021-04-27 | 2024-03-15 | 中国科学院国家空间科学中心 | Ship track real-time prediction method and system based on optimal strategy learning |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3725918A (en) * | 1970-11-18 | 1973-04-03 | Sperry Rand Corp | Collision avoidance display apparatus for maneuverable craft |
WO1995003598A1 (en) * | 1993-07-20 | 1995-02-02 | Philip Bernard Wesby | Locating/map dissemination system |
JP2005181077A (en) * | 2003-12-18 | 2005-07-07 | Tokimec Inc | Navigation supporting system for vessel |
CN101639988A (en) * | 2009-05-15 | 2010-02-03 | 天津七一二通信广播有限公司 | Method for preventing boats from colliding |
CN102831787A (en) * | 2012-07-30 | 2012-12-19 | 广东省公安边防总队 | Quay berth-based intelligent supervision method and system |
CN102890875A (en) * | 2012-10-15 | 2013-01-23 | 浙江大学 | Method for acquiring system collision risk of maritime intelligent transportation network |
CN103714718A (en) * | 2013-12-31 | 2014-04-09 | 武汉理工大学 | Inland river bridge area ship safe navigation precontrol system |
CN104010167A (en) * | 2014-06-16 | 2014-08-27 | 交通运输部天津水运工程科学研究所 | Real-time virtual ship video displaying method for AIS data |
CN104021700A (en) * | 2014-06-16 | 2014-09-03 | 交通运输部天津水运工程科学研究所 | Ship safety alarm device and method based on TETRA digital trunking system |
CN104091470A (en) * | 2014-07-15 | 2014-10-08 | 南京大学 | Channel traffic information prediction method and application based on multidata fusion |
KR20150064909A (en) * | 2013-12-04 | 2015-06-12 | 한국전자통신연구원 | Apparatus and method for tracking vessel using ais target information |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002312900A (en) * | 2001-04-13 | 2002-10-25 | Original Engineering Consultants Co Ltd | Three-dimensional position confirmation system |
JP4461450B2 (en) * | 2001-04-24 | 2010-05-12 | 日本無線株式会社 | Weather information provision system |
US7643989B2 (en) * | 2003-08-29 | 2010-01-05 | Microsoft Corporation | Method and apparatus for vocal tract resonance tracking using nonlinear predictor and target-guided temporal restraint |
US7983835B2 (en) * | 2004-11-03 | 2011-07-19 | Lagassey Paul J | Modular intelligent transportation system |
JP4709684B2 (en) * | 2006-04-13 | 2011-06-22 | 隼馬 今津 | Counterpart movement monitoring device |
JP2010033352A (en) * | 2008-07-29 | 2010-02-12 | Toyota Central R&D Labs Inc | Lane change alarm and program |
WO2012012550A2 (en) * | 2010-07-20 | 2012-01-26 | The University Of Memphis Research Foundation | Theft detection nodes and servers, methods of estimating an angle of a turn, methods of estimating a distance traveled between successive stops, and methods and servers for determining a path traveled by a node |
US8405531B2 (en) * | 2010-08-31 | 2013-03-26 | Mitsubishi Electric Research Laboratories, Inc. | Method for determining compressed state sequences |
CN102147981B (en) * | 2010-12-20 | 2014-01-15 | 成都天奥信息科技有限公司 | Method for warning of warning region of shipborne automatic identification system |
US8478711B2 (en) * | 2011-02-18 | 2013-07-02 | Larus Technologies Corporation | System and method for data fusion with adaptive learning |
US20130080365A1 (en) * | 2011-04-13 | 2013-03-28 | The Board Of Trustees Of The Leland Stanford Junior University | Phased Whole Genome Genetic Risk In A Family Quartet |
CN102568200B (en) * | 2011-12-21 | 2015-04-22 | 辽宁师范大学 | Method for judging vehicle driving states in real time |
CN103152253A (en) * | 2013-01-17 | 2013-06-12 | 王少夫 | Mining path prediction algorithm based on complex communication network channel data |
CN103235933B (en) * | 2013-04-15 | 2016-08-03 | 东南大学 | A kind of vehicle abnormality behavioral value method based on HMM |
CN103336863B (en) * | 2013-06-24 | 2016-06-01 | 北京航空航天大学 | The flight intent recognition methods of flight path observed data of flying based on radar |
CN103413443B (en) * | 2013-07-03 | 2015-05-20 | 太原理工大学 | Short-term traffic flow forecasting method based on hidden Markov model |
CN103473540B (en) * | 2013-09-11 | 2016-06-22 | 天津工业大学 | The modeling of intelligent transportation system track of vehicle increment type and online method for detecting abnormality |
CN103471589B (en) * | 2013-09-25 | 2015-10-21 | 武汉大学 | The method of the identification of a kind of indoor pedestrian's walking mode and trajectory track |
CN106251704A (en) * | 2014-12-30 | 2016-12-21 | 江苏理工学院 | Boats and ships trajectory predictions method |
-
2014
- 2014-12-30 CN CN201610620526.9A patent/CN106251704A/en active Pending
- 2014-12-30 CN CN201610624893.6A patent/CN106228850A/en active Pending
- 2014-12-30 CN CN201410841564.8A patent/CN104537891B/en active Active
- 2014-12-30 CN CN201610620437.4A patent/CN106205213A/en active Pending
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3725918A (en) * | 1970-11-18 | 1973-04-03 | Sperry Rand Corp | Collision avoidance display apparatus for maneuverable craft |
WO1995003598A1 (en) * | 1993-07-20 | 1995-02-02 | Philip Bernard Wesby | Locating/map dissemination system |
JP2005181077A (en) * | 2003-12-18 | 2005-07-07 | Tokimec Inc | Navigation supporting system for vessel |
CN101639988A (en) * | 2009-05-15 | 2010-02-03 | 天津七一二通信广播有限公司 | Method for preventing boats from colliding |
CN102831787A (en) * | 2012-07-30 | 2012-12-19 | 广东省公安边防总队 | Quay berth-based intelligent supervision method and system |
CN102890875A (en) * | 2012-10-15 | 2013-01-23 | 浙江大学 | Method for acquiring system collision risk of maritime intelligent transportation network |
KR20150064909A (en) * | 2013-12-04 | 2015-06-12 | 한국전자통신연구원 | Apparatus and method for tracking vessel using ais target information |
CN103714718A (en) * | 2013-12-31 | 2014-04-09 | 武汉理工大学 | Inland river bridge area ship safe navigation precontrol system |
CN104010167A (en) * | 2014-06-16 | 2014-08-27 | 交通运输部天津水运工程科学研究所 | Real-time virtual ship video displaying method for AIS data |
CN104021700A (en) * | 2014-06-16 | 2014-09-03 | 交通运输部天津水运工程科学研究所 | Ship safety alarm device and method based on TETRA digital trunking system |
CN104091470A (en) * | 2014-07-15 | 2014-10-08 | 南京大学 | Channel traffic information prediction method and application based on multidata fusion |
Non-Patent Citations (1)
Title |
---|
杨君兰: "基于复杂度建模的船舶碰撞预警研究", 《CNKI优秀硕士学位论文全文库》, 1 April 2013 (2013-04-01) * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106251704A (en) * | 2014-12-30 | 2016-12-21 | 江苏理工学院 | Boats and ships trajectory predictions method |
CN106154296A (en) * | 2015-06-26 | 2016-11-23 | 安徽华米信息科技有限公司 | The method of adjustment of a kind of path locus and device |
CN106595665A (en) * | 2016-11-30 | 2017-04-26 | 耿生玲 | Prediction method for spatial-temporal trajectory of moving object in obstructed space |
CN106595665B (en) * | 2016-11-30 | 2019-10-11 | 耿生玲 | The prediction technique of mobile object space-time trajectory in a kind of space with obstacle |
CN112070312A (en) * | 2020-09-10 | 2020-12-11 | 龙马智芯(珠海横琴)科技有限公司 | Flight path prediction method, flight path prediction device, and computer-readable storage medium |
CN112070312B (en) * | 2020-09-10 | 2021-11-02 | 龙马智芯(珠海横琴)科技有限公司 | Flight path prediction method, flight path prediction device, and computer-readable storage medium |
CN113627104A (en) * | 2021-08-12 | 2021-11-09 | 北京中安智能信息科技有限公司 | Underwater submarine track simulation method, device and equipment under multi-constraint condition |
CN113627104B (en) * | 2021-08-12 | 2024-02-06 | 北京中安智能信息科技有限公司 | Underwater submarine track simulation method, device and equipment under multiple constraint conditions |
Also Published As
Publication number | Publication date |
---|---|
CN106228850A (en) | 2016-12-14 |
CN104537891B (en) | 2016-08-24 |
CN106251704A (en) | 2016-12-21 |
CN106205213A (en) | 2016-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104537891B (en) | A kind of boats and ships track real-time predicting method | |
CN104484726A (en) | Ship track real-time prediction method | |
CN104462856A (en) | Ship conflict early warning method | |
CN104504277A (en) | Early-warning method for collision of ship | |
Liu et al. | STMGCN: Mobile edge computing-empowered vessel trajectory prediction using spatio-temporal multigraph convolutional network | |
Sheng et al. | Research on ship classification based on trajectory features | |
Perera et al. | Data analysis on marine engine operating regions in relation to ship navigation | |
CN111736595B (en) | Control method, system and device for ship running track and storage medium | |
CN105654139A (en) | Real-time online multi-target tracking method adopting temporal dynamic appearance model | |
CN114022847A (en) | Intelligent agent trajectory prediction method, system, equipment and storage medium | |
CN104504934A (en) | Marine traffic control method | |
CN104504935A (en) | Maritime traffic control method | |
CN109359661B (en) | Sentinel-1 radar image classification method based on convolutional neural network | |
CN112052802A (en) | Front vehicle behavior identification method based on machine vision | |
CN116383685A (en) | Vehicle lane change detection method based on space-time interaction diagram attention network | |
CN113706581A (en) | Target tracking method based on residual channel attention and multilevel classification regression | |
CN109711431A (en) | The method for tracking target of local block convolution, system and storage medium at one | |
Liu et al. | Ship route planning in the pirate area via hybrid probabilistic roadmap algorithm within the context of the Maritime Silk Road | |
Liu et al. | Intelligent data-driven vessel trajectory prediction in marine transportation cyber-physical system | |
CN115083199A (en) | Parking space information determination method and related equipment thereof | |
CN112350733B (en) | Rapid track compression method based on effective distance | |
Alajlan et al. | Automatic lane marking prediction using convolutional neural network and S-Shaped Binary Butterfly Optimization | |
US20220262097A1 (en) | Label estimation device, label estimation method, and label estimation program | |
CN105160677A (en) | Multi-feature image segmenting and positioning method suitable for spatial non-cooperative targets | |
Zhang et al. | Motion Prediction of Tugboats Using Hidden Markov Model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |