JP6243308B2 - 特定装置および特定方法 - Google Patents

Description

本発明は、情報を特定する特定装置および特定方法に関する。

船舶の安全航行のため、船舶自動識別装置（Ａｕｔｏｍａｔｉｃ Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ Ｓｙｓｔｅｍ，ＡＩＳ）が搭載されている船舶がある。船舶の位置や速度などを含む航跡データ（以降、ＡＩＳデータとする）が、ＡＩＳから定期的に近隣の港湾や船舶に無線で発信される。近年、衛星を通してＡＩＳデータを収集する仕組みが導入され、全世界のＡＩＳデータをより簡単に入手可能である。そこで、地図上にリアルタイムで船舶の位置情報を表示する既存システムに対し、ＡＩＳデータを用いて海運業界の経営戦略策定や船舶運行管理を支援する付加機能が要望されている。

例えば、船舶の到着遅延により貨物ターミナルのゲート入口付近で渋滞を回避するためのＡＩＳデータを用いて船舶の予定到着時刻をリアルタイムに提供する船舶動静予測技術（たとえば、下記特許文献１を参照。）や、燃費を削減するための船舶の就航データと予測された気象データを用いて最適航路を計算する技術（たとえば、下記特許文献２を参照）などがある。

特開２００９−０２５８６０号公報 特開２００９−２８６２３０号公報

荷主は、海運会社や用船仲介者に個別に問い合わせることで、指定された時刻から一定期間内に、指定された海域（以降では、対象海域）に到達する船舶を列挙して、用船交渉を行っている。

荷主にとって、用船交渉をする際に交渉可能な船舶や海運会社が増えるほど、選択肢が増え、用船料の削減交渉に有利になる。また、荷主は、積荷を運びたい目的港に実際に到着する船舶以外の他の船舶も考慮することができるのであれば、さらに選択肢を増やすことができる。たとえば、他の船舶とは、目的港に到着すべく対象海域に向かっており、要望があれば一定期間内に目的港とは異なる他の港に到着できる船舶である。

また、荷主（発荷主または着荷主のいずれもでよい）にとって、将来どのタイミングで用船契約するのが最も良いかを検討するとき、将来の各タイミングに目的港に到着する船舶を予測するだけでは、用船可能の船舶や海運会社が限定されてしまい、用船契約の交渉を有利に運べないという問題がある。また、上述した特許文献１，２の従来技術では、対象海域への到達以前の指定された時刻から一定期間内に対象海域に到達する船舶を確認することができないため、用船可能の船舶や海運会社が限定されてしまう。

本発明は、対象海域への到達以前の指定された時刻から一定期間内に対象海域に到着する船舶を予測することを目的とする。

本願において開示される発明の一側面となる特定装置および特定方法では、プロセッサは、対象海域内の目的港へ航行する船舶の前記目的港に到着するまでの時系列な位置および時刻を示す一連の時系列データを取得する取得処理と、前記取得処理によって取得された一連の時系列データの中の連続する２つの時系列データと、前記一連の時系列データのうち前記船舶が前記目的港に位置するときの時系列データが示す到着時刻と、予測期間と、に基づいて、前記一連の時系列データのうち前記船舶が前記予測期間以内で前記対象海域に進入できる時系列データの中の最古の時系列データの位置を、前記対象海域からの前記船舶の最遠位置として特定する特定処理と、を実行することを特徴とする。

本発明の代表的な実施の形態によれば、対象海域への到達以前の指定された時刻から一定期間内に対象海域に到着する船舶を予測することができる。前述した以外の課題、構成及び効果は、以下の実施例の説明により明らかにされる。

海域を含む地図データを示す説明図である。 特定装置の構成例を示すブロック図である。 データメモリに含まれるＡＩＳデータのデータ構造例を示す説明図である。 時系列データのデータ構造例を示す説明図である。 データメモリに含まれる港到着履歴データのデータ構造例を示す説明図である。 データメモリに含まれる海域データのデータ構造例１を示す説明図である。 データメモリに含まれる海域データのデータ構造例２を示す説明図である。 最遠位置データのデータ構造例を示す説明図である。 データメモリに含まれる予測船舶データのデータ構造を示す説明図である 特定装置による特定処理手順例を示すフローチャートである。 航海データ取得処理（ステップＳ１００４）の詳細な処理手順例を示すフローチャートである。 最遠位置特定処理（ステップＳ１００５）の具体例を示す説明図である。 最遠位置特定処理（ステップＳ１００５）の詳細な処理手順例を示すフローチャートである。 例外除去処理（ステップＳ１００６）の詳細な処理手順例を示すフローチャートである。 航行枠生成処理（ステップＳ１００７）の詳細な処理手順例を示すフローチャートである。 到着船舶予測処理（ステップＳ１０１３）の詳細な処理手順例を示すフローチャートである。 選択条件での到着船舶予測処理（ステップＳ１６０３）の詳細な処理手順例を示すフローチャートである。 予測船舶データ生成処理（ステップＳ１７０６）における予測期間ごとの海域範囲に含まれる船舶を特定する例を示す説明図である。 予測船舶データ生成処理（ステップＳ１７０６）の詳細な処理手順例を示すフローチャートである。 図１０の予測表示処理（Ｓ１０１４）における航行枠表示処理部により表示される表示画面例１を示す説明図である。 図１０の予測表示処理（Ｓ１０１４）における航行枠表示処理部により表示される表示画面例２を示す説明図である。 図１０の予測表示処理（Ｓ１０１４）における予測表示処理部による表示画面例１を示す説明図である。 図１０の予測表示処理（Ｓ１０１４）における予測表示処理部による表示画面例２を示す説明図である。 図１０の予測表示処理（Ｓ１０１４）における予測表示処理部による表示画面例３を示す説明図である。

以下、添付図面を参照しながら、本発明の特定装置を実施するための最良の実施形態を詳細に説明する。図１〜図２４は、本発明の実施例を示す説明図である。これらの図において、同一の符号を付した部分は同一物を表し、基本的な構成および動作は同様であるものとする。尚、本発明の実施形態において、使用される機器、手法等は一例であり、本発明はこれらに限定されるものではないことは勿論である。

更に、本発明の実施形態は、後述されるように、汎用コンピュータ上で稼動するソフトウェアで実装しても良いし専用ハードウェア又はソフトウェアとハードウェアの組み合わせで実装しても良い。

なお、以後の説明では「テーブル」形式によって本発明の各情報について説明するが、これら情報は必ずしもテーブルによるデータ構造で表現されていなくても良く、リスト、ＤＢ、キュー等のデータ構造やそれ以外で表現されていても良い。そのため、データ構造に依存しないことを示すために「テーブル」、「リスト」、「ＤＢ」、「キュー」等について単に「情報」と呼ぶことがある。

また、各情報の内容を説明する際に、「識別情報」、「識別子」、「名」、「名前」、「ＩＤ」という表現を用いることが可能であり、これらについてはお互いに置換が可能である。

以下では「プログラム」を主語（動作主体）として本発明の実施形態における各処理について説明を行うが、プログラムはプロセッサによって実行されることで定められた処理をメモリおよび通信ポート（通信制御装置）を用いながら行うため、プロセッサを主語とした説明としてもよい。また、プログラムを主語として開示された処理は管理サーバ等の計算機、情報処理装置が行う処理としてもよい。プログラムの一部またはすべては専用ハードウェアで実現してもよく、また、モジュール化されていても良い。各種プログラムはプログラム配布サーバや記憶メディアによって各計算機にインストールされてもよい。

＜対象海域と航行枠の一例＞
図１は、海域を含む地図データを示す説明図である。特定装置１００は、地図データ１を記憶する。地図データ１には、対象海域２（ここでは、例として東南アジア海域）が表示される。米印である符号３は、船舶の位置を示す位置データを示す。位置データ３を中心とする円４は、船舶位置の誤差を許容する許容範囲である。また、円４を内接する太線の多角形は、特定装置１００により作成された航行枠５である。

本実施例の特定装置１００は、たとえば、荷主（発荷主または着荷主のいずれもでよい）に利用される。本実施例の特定装置１００は、過去の航行履歴を基に、船舶ごとに海域に到達（または進入）する時点から所定期間前の船舶位置を航行枠５として設定して表示する。すなわち、位置データ３で特定される航行枠の境界に位置する船舶は、所定期間経過後に対象海域２に到達（または進入）することを意味する。これにより、荷主は、対象海域２内の船舶のみならず航行枠５と対象海域２との間に位置する船舶を把握することができる。したがって、用船交渉先となる可能な船舶や海運会社といった選択肢が増え、用船料の削減交渉が荷主にとって有利になる。

＜特定装置１００の構成例＞
図２は、特定装置１００の構成例を示すブロック図である。特定装置１００は、ＡＩＳデータＤＢ２００と、港到着履歴データＤＢ２０１と、海域データＤＢ２０２と、予測船舶データＤＢ２０３と、入力装置２０４と、表示装置２０５と、中央処理装置２０６と、プログラムメモリ２０７と、中央処理装置（プロセッサ）２０６の処理に必要なデータを格納するデータメモリ２０８と、を有する。特定装置１００は、１台以上のコンピュータにより構成される。

ＡＩＳデータＤＢ２００は、ＡＩＳデータ２１８を記憶するデータベースである（図３、図４を参照。）。港到着履歴データＤＢ２０１は、港到着履歴データ２１９を記憶するデータベースである（図５を参照。）。海域データＤＢ２０２は、海域データ２２０を記憶するデータベースである（図６〜図８を参照。）。予測船舶データＤＢ２０３は、予測船舶データを記憶するデータベースである（図９を参照。）。

中央処理装置２０６は、予測処理部２０９と、予測表示処理部２１０と、を有する。予測処理部２０９および予測表示処理部２１０は、コンピュータ上で実行されるプログラムの機能の一部として実現される。当該プログラムは、プログラムメモリ２０７に格納され、処理実行の際に中央処理装置２０６が内部メモリに読み込むようにしている。

予測処理部２０９は、航海データ取得処理部２１１と、最遠位置特定処理部２１２と、例外除去処理部２１３と、航行枠生成処理部２１４と、到着船舶予測処理部２１５と、を有する。また、予測表示処理部２１０は、航行枠表示処理部２１６と、到着船舶予測表示処理部２１７と、を有する。

データメモリ２０８は、ＡＩＳデータＤＢ２００から読み込まれたＡＩＳデータ２１８と、港到着履歴データＤＢ２０１から読み込まれた港到着履歴データ２１９と、海域データＤＢ２０２から読み込まれた、または、航行枠生成処理部２１４より生成された海域データ２２０と、到着船舶予測処理部２１５より予測された予測船舶データ２２１と、を格納する。

＜データ構造例＞
つぎに、上述したＡＩＳデータ２１８、港到着履歴データ２１９、海域データ２２０、および予測船舶データ２２１のデータ構造例について、図３〜図９を用いて説明する。

（ＡＩＳデータ２１８）
図３は、データメモリ２０８に含まれるＡＩＳデータ２１８のデータ構造例を示す説明図である。ＡＩＳデータ２１８は、ＡＩＳデータＤＢ２００に記憶されており、中央処理装置２０６により読み出されてデータメモリ２０８に格納される。

ＡＩＳデータ２１８は、船舶番号３００と、ＩＭＯ（Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｍａｒｉｔｉｍｅ Ｏｒｇａｎｉｚａｔｉｏｎ、国際海事機関）番号（船舶識別番号）３０１と、船舶名３０２と、船種３０３と、サイズ３０４と、最大積載量３０５と、海運会社３０６と、建船日付３０７と、時系列データ３０８と、を含み、たとえば配列の形式でこれらの情報が保有される。

船舶番号３００には、ＡＩＳデータ２１８に対して一意に振られた値が保持される。ＩＭＯ番号３０１には、値として、ＩＭＯにおいて大型船に割り振られた船舶の識別番号が保持される。船舶名３０２には、値として、ＩＭＯに登録した船の名前が保持される。船種３０３には、値として、ＩＭＯに登録した船舶の種類名称が保持される。サイズ３０４には、値として、ＩＭＯに登録した船舶の長さや幅などが保持される。最大積載量３０５には、値として、船舶別に積載可能の最大容量が保持される。海運会社３０６には、値として、船舶のオペレーター会社の名称が保持される。建船日付３０７には、値として、船舶が造られた日付が保持される。時系列データ３０８には、時刻ごとの船舶位置や速度などの情報が保持される。

（時系列データ３０８）
図４は、時系列データ３０８のデータ構造例を示す説明図である。時系列データ３０８は、時系列番号４００と、受信時刻４０１と、経度４０２と、緯度４０３と、速度４０４と、喫水４０５と、を含み、たとえば配列の形式でこれらの情報が保有される。時系列番号４００には、時系列データ３０８に対して一意に振られた値が保持される。受信時刻４０１には、値として、該当時系列データが受信された時刻が保持される。経度４０２および緯度４０３には、値として、該当時系列データ３０８が発信されたときの船舶の座標情報が保持される。速度４０４には、値として、該当時系列データ３０８が発信されたときの船舶の速度が保持される。喫水４０５には、値として、該当時系列データ３０８が発信されたときに船舶が水上にある際に船体が沈む深さが保持される。

（港到着履歴データ２１９）
図５は、データメモリ２０８に含まれる港到着履歴データ２１９のデータ構造例を示す説明図である。港到着履歴データ２１９は、港到着履歴データＤＢ２０１に記憶されており、中央処理装置２０６により読み出されてデータメモリ２０８に格納される。港到着履歴データ２１９は、履歴番号５００と、港番号５０１と、港名５０２と、国５０３と、港経度５０４と、港緯度５０５と、船舶番号５０６と、到着時刻５０７と、を含み、たとえば配列の形式でこれらの情報が保有される。

履歴番号５００には、港到着履歴データ２１９に対して一意に振られた値が保持される。港番号５０１には、値として、世界中の港に対して振られた唯一の値が保持される。港名５０２には、値として、世界共通の港名称が保持される。国５０３には、値として、港の所属国の名称が保持される。港経度５０４および港緯度５０５には、値として、港の座標情報が保持される。船舶番号５０６には、該当港に到着した該当船舶における船舶番号３００が保持される。到着時刻５０７には、値として、船舶が該当港に到着したときの時刻が保持される。

（海域データ２２０）
図６および図７は、データメモリ２０８に含まれる海域データ２２０のデータ構造を示す説明図である。海域データ２２０は、海域データＤＢ２０２に記憶されており、中央処理装置２０６により読み出されてまたは算出されてデータメモリ２０８に格納される。海域データ２２０は、種類６０１によって値の構成が異なるため、種類６０１ごとに説明する。図６の海域データ２２０は、種類６０１の値が「対象海域」である場合のデータ構造例であり、図７の海域データ２２０は、種類６０１の値が「航行枠」である場合のデータ構造例である。なお、データ構造自体は、図６も図７も同一である。

図６および図７において、海域データ２２０は、海域番号６００と、種類６０１と、対象海域番号６０２と、海域名６０３と、船種６０４と、海域範囲６０５と、季節６０６と、予測期間６０７と、最遠位置データ６０８と、を含み、たとえば配列の形式でこれらの情報が保有される。海域番号６００には、海域データ２２０に対して一意に振られた値が保持される。

種類６０１には、値として、海域データ２２０の種類を識別するための値が保持される。上述したように、海域データ２２０には、「対象海域」（図６）と「航行枠」（図７）の２種類である。種類が「対象海域」である海域データ２２０は、海域番号６００および海域名６０３の各値により特定される海域を示す海域データ２２０であり、あらかじめ用意される。一方、種類６０１が「航行枠」である海域データ２２０は、対象海域番号６０２の値である対象海域２に対して航行枠５を規定する海域データ２２０であり、本実施例の特定装置１００により生成される。

対象海域番号６０２には、値として、「航行枠」種類の海域データ２２０において対応する「対象海域」種類の海域データ２２０の海域番号６００が保持される。なお、「対象海域」種類の海域データ２２０の場合、例として、対象海域番号６０２に「−」が登録される。

海域名６０３には、値として、海域の名称が保持される。なお、「航行枠」種類の海域データ２２０の場合、航行枠５は海域外であるため、例として、海域名６０３に「−」が登録される。

船種６０４には、値として、予測対象船舶の種類名称が保持される。「対象海域」種類の海域データの場合、予測対象船舶が特定されていないため、例として、船種６０４に「−」が登録される。

海域範囲６０５には、種類６０１によって特定される領域の範囲を示す図形情報が保持される。たとえば、種類６０１が「対象海域」である場合、海域番号６００および海域名６０３で特定される対象海域２を特定する図形情報（たとえば、円）が保持される。一方、種類６０１が「航行枠」である場合、航行枠５を特定する図形情報（たとえば、多角形）が保持される。

季節６０６には、値として、季節の名称が保持される。なお、「対象海域」種類の海域データ２２０の場合、例として、季節６０６に「−」が登録される。

予測期間６０７には、値として、船舶が、「航行枠」種類の海域データ２２０が示す航行枠５から、「対象海域」種類の海域データ２２０が示す対象海域２に進入するまで必要な期間が保持される。具体的には、図１０で説明する設定処理が実行される場合に選択される予測期間ΔＴが格納される。なお、「対象海域」種類の海域データの場合、例として、予測期間６０７に「０」が登録される。

最遠位置データ６０８には、対象海域２に進入した直後から予測期間６０７の値である一定期間前の船舶位置などの情報が保持される。

（最遠位置データ６０８）
図８は、最遠位置データ６０８のデータ構造例を示す説明図である。最遠位置データ６０８は、最遠位置番号８００と、受信時刻８０１と、経度８０２と、緯度８０３と、例外フラグ８０４と、を含み、たとえば、配列の形式でこれらの情報が保有される。なお、最遠位置データ６０８は、時系列データ３０８を基に生成されるデータである。

最遠位置番号８００には、最遠位置データ６０８に対して一意に振られた値が保持される。受信時刻８０１には、該当時系列データ３０８における受信時刻４０１が保持される。経度８０２および緯度８０３には、該当時系列データ３０８における経度４０２および緯度４０３が保持される。例外フラグ８０４には、該当「航行枠」種類の海域データ２２０の海域範囲６０５を生成するときに、最遠位置データ６０８が使われている可否かを識別するための値が保持される。なお、デフォルトの場合、例として、例外フラグ８０４に「Ｆａｌｓｅ」が保持され、海域範囲６０５を生成するときに使われることを明記する。

（予測船舶データ２２１）
図９は、データメモリ２０８に含まれる予測船舶データ２２１のデータ構造を示す説明図である。予測船舶データ２２１は、中央処理装置２０６により生成されて、予測船舶データＤＢ２０３に記憶される。また、予測船舶データ２２１は、中央処理装置２０６により読み出されてデータメモリ２０８に格納される。

予測船舶データ２２１は、予測船舶番号９００と、対象海域番号９０１と、船舶番号９０２と、ＩＭＯ番号９０３と、船舶名９０４と、船種９０５と、最大積載量９０６と、海運会社９０７と、受信時刻９０８と、経度９０９と、緯度９１０と、速度９１１と、喫水９１２と、必要期間９１３と、予測日付９１４と、を含み、たとえば、配列の形式でこれらの情報が保有される。なお、予測船舶データ２２１は、ＡＩＳデータ２１８を基に生成される。

予測船舶番号９００には、予測船舶データ２２１に対して一意に振られた値が保持される。対象海域番号９０１には、該当する「対象海域」種類の海域データ２２０における海域番号６００が保持される。

船舶番号９０２には、該当船舶における船舶番号３００が保持される。ＩＭＯ番号９０３には、該当船舶におけるＩＭＯ番号３０１が保持される。船舶名９０４には、該当船舶における船舶名３０２が保持される。船種９０５には、該当船舶における船種３０３が保持される。最大積載量９０６には、該当船舶における最大積載量３０５が保持される。

海運会社９０７には、該当船舶における海運会社３０６が保持される。受信時刻９０８には、予測時点に該当船舶の最新の時系列データ３０８における受信時刻４０１が保持される。経度５０９および緯度５１０には、予測時点に該当船舶の最新の時系列データ３０８における経度４０２および緯度４０３が保持される。

速度５１２には、予測時点での該当船舶の最新の時系列データ３０８における速度４０４が保持される。喫水９１２には、予測時点での該当船舶の最新の時系列データ３０８における喫水４０５が保持される。必要期間９１３には、予測された船舶が該当対象海域に入るまで必要な航行期間が保持される。予測日付９１４には、該当対象海域において到着船舶を予測するときの日付が保持される。

＜特定処理手順例＞
つぎに、特定装置１００による特定処理手順例について図１０〜図２４を用いて説明する。

図１０は、特定装置１００による特定処理手順例を示すフローチャートである。図１０のフローチャートは、指定された船種において、受信されたＡＩＳデータ２１８と港到着履歴データ２１９に基づいて一定期間（以降、予測期間と略す）以内に対象海域２に到着する可能の航行枠５を生成し、航行枠５とＡＩＳデータ２１８の座標値における位置関係を求め、予測期間以内に対象海域２に到着する船舶を予測する処理例である。これにより、各船種、各海域および各予測期間において、荷主は、対象海域２に到着可能な船舶を確認することができる。

なお、ここでは、図１０のステップＳ１００１〜ステップＳ１０１４までの処理を、たとえば、１日１回のバッチ処理で実行する例について説明するが、バッチ処理の間隔は任意に変更することができる。また、バッチ処理に限られず、任意のタイミングでも実行可能である。

図１０において、まず、特定装置１００は、海域データＤＢ２０２から「対象海域」種類の海域データ２２０をデータメモリ２０８に格納する（ステップＳ１０００）。つぎに、特定装置１００は、複数の季節（たとえば、春夏秋冬）から未選択の季節Ｘを選択する（ステップＳ１００１）。また、特定装置１００は、複数の予測期間（たとえば、１週間、２週間、３週間、４週間）から未選択の予測期間ΔＴを選択する（ステップＳ１００２）。本例の場合、季節Ｘが４種類、予測期間ΔＴも４種類であるため、１６通りの組み合わせについて以降の処理が実行される。

つぎに、特定装置１００は、データメモリ２０８から１件の「対象海域」種類の海域データ２２０である海域データＡ１を取得する（ステップＳ１００３）。そして、特定装置１００は、航海データ取得処理部２１１を用いて、航海データ取得処理を実行する（ステップＳ１００４）。航海データ取得処理（ステップＳ１００４）は、港到着履歴データ２１９に基づいて、季節Ｘに海域データＡ１内の港へ航行した船舶の航海データを抽出する処理である。航海データは、ＡＩＳデータ２１８の時系列データ３０８である。

例えば、対象海域２である東南アジア海域の夏は６月から８月までとすると、選択季節Ｘ＝夏に東南アジア海域内の港に到着した航海データを抽出するため、特定装置１００は、港到着履歴データ２１９を基づいて、６月から８月に東南アジア海域内の港に到着した船舶の航海データを当該船舶のＡＩＳデータ２１８から抽出する。なお、航海データ取得処理（ステップＳ１００４）についての処理は、図１１を用いて詳細に後述する。

つぎに、特定装置１００は、航海データ取得処理部２１１より得られた航海データにおいて、最遠位置特定処理部２１２を用いて、最遠位置特定処理を実行する（ステップＳ１００５）。最遠位置特定処理（ステップＳ１００５）は、海域データＡ１に入った直後から予測期間分遡った船舶位置を特定する処理である。例えば、特定装置１００は、対象海域２である東南アジア海域に入った直後から１週間前の船舶位置を特定する。なお、ステップＳ１００４についての処理は、図１３を用いて詳細に後述する。

つぎに、特定装置１００は、最遠位置特定処理部２１２より得られた最遠位置データ４０８について、例外除去処理部２１３を用いて、例外除去処理を実行する（ステップＳ１００５）。例外除去処理（ステップＳ１００６）は、他の最遠位置データ４０８から遠く離れた最遠位置データ４０８を表示対象から除去できるよう設定する処理である。なお、例外除去処理（ステップＳ１００６）は、図１４を用いて詳細に後述する。

つぎに、特定装置１００は、航行枠生成処理部２１４を用いて、例外を除外した後の最遠位置データ４０８を基に、航行枠生成処理を実行する（ステップＳ１００７）。航行枠生成処理（ステップＳ１００７）は、「航行枠」種類の海域データＡ１を作成する処理である。なお、航行枠生成処理（ステップＳ１００７）の詳細は、図１５を用いて詳細に後述する。

つぎに、特定装置１００は、航行枠生成処理（ステップＳ１００６）によって新しく生成された「航行枠」種類の海域データ２２０を海域データＤＢ２０２に書込む（ステップＳ１００８）。具体的には、図７および図８に示したような海域データ２２０および最遠位置データ６０８が海域データＤＢ２０２に書き込まれる。このあと、特定装置１００は、データメモリ２０８にある「航行枠」種類の海域データ２２０を削除する（ステップＳ１００９）。

つぎに、特定装置１００は、データメモリ２０８上に、未取得の「対象海域」種類の海域データ２２０があるか否かを判断する（ステップＳ１０１０）。未取得の「対象海域」種類の海域データ２２０が有る場合（ステップＳ１０１０：Ｙｅｓ）、ステップＳ１００３に戻り、海域データＡ１の再取得を試行する。

一方、未取得の海域データ２２０がない場合（ステップＳ１０１０：Ｎｏ）、未処理のすべての「対象海域」種類の海域データ２２０が処理されたこととなる。この場合、特定装置１００は、未選択の予測期間ΔＴがあるか否かを判断する（ステップＳ１０１１）。未選択の予測期間ΔＴがある場合（ステップＳ１０１１：Ｙｅｓ）、ステップＳ１００２に戻り、特定装置１００は、予測期間を再選択する（ステップＳ１００２）。

一方、未選択の予測期間ΔＴがない場合（ステップＳ１１１１：Ｎｏ）、特定装置１００は、到着船舶予測処理部２１５を用いて、対象海域２において、到着船舶予測処理を実行する（ステップＳ１０１３）。

到着船舶予測処理（ステップＳ１０１３）は、予測期間以内に対象海域２に到着する船舶を予測する処理である。例えば、特定装置１００は、１週間（指定期間）以内に東南アジア海域（対象海域２）に到着する船舶を予測する。なお、到着船舶予測処理（ステップＳ１０１３）の詳細は、図１６〜図１８を用いて詳細に後述する。

このあと、ユーザ（たとえば、荷主）からの操作により、特定装置１００は、航行枠表示処理部２１６と到着船舶予測表示処理部２１７を用いて、予測結果を表示する（ステップＳ１０１４）。なお、航行枠表示処理部２１６と到着船舶予測表示処理部２１７（ステップＳ１０１２）の表示内容例と処理は、図２０〜図２４を用いて詳細に後述する。

（航海データ取得処理（ステップＳ１００４））
図１１は、航海データ取得処理（ステップＳ１００４）の詳細な処理手順例を示すフローチャートである。航海データ取得処理（ステップＳ１００４）は、対象海域２内の目的港へ航行する船舶の目的港に到着するまでの時系列な位置および時刻を示す一連の時系列データを取得する。具体的には、航海データ取得処理（ステップＳ１００４）は、指定された船種において、受信されたＡＩＳデータと港到着履歴データに基づいて、該当季節に対象海域にある港へ航行した航海データを抽出する処理である。

図１１において、まず、特定装置１００は、港到着履歴データＤＢ２０１から港到着履歴データ２１９をデータメモリに格納する（ステップＳ１１００）。つぎに、特定装置１００は、データメモリから１件の港到着履歴データＰ１を取得し（ステップＳ１１０１）、港到着履歴データＰ１の港経度５０４および港緯度５０５が「対象海域」種類の海域データＡ１の海域範囲６０５にあるか否かを判断する（ステップＳ１１０２）。

港到着履歴データＰ１の港経度５０４または港緯度５０５が海域データＡ１の海域範囲６０５の外にあった場合（ステップＳ１１０２：Ｎｏ）、ステップＳ１１０４に移行する。一方、Ｐ１の港経度５０４および港緯度５０５がＡ１の海域範囲６０５にあった場合（ステップＳ１１０２：Ｙｅｓ）、特定装置１００は、港到着履歴データＰ１の到着時刻５０７が季節Ｘに属するか否かを判断する（ステップＳ１１０３）。

港到着履歴データＰ１の到着時刻５０７が季節Ｘに属しない場合（ステップＳ１１０３：Ｎｏ）、特定装置１００は、港到着履歴データＰ１をデータメモリから削除する（ステップＳ１１０４）。一方、季節Ｘに属する場合（ステップＳ１１０３：Ｙｅｓ）、特定装置１００は、ＡＩＳデータＤＢ２００から港到着履歴データＰ１の船舶番号と一致する船舶のＡＩＳデータ２１８をデータメモリ２０８に格納する（ステップＳ１１０５）。

そして、特定装置１００は、港到着履歴データＰ１の到着時刻５０７以降に受信されたＡＩＳデータ２１８の時系列データ３０８をデータメモリ２０８から削除して（ステップＳ１１０６）、ステップＳ１１０７に移行する。例えば、ある船舶のＡＩＳデータ２１８には、受信時刻４０１が「２０１２−０２−０１」から「２０１４−０７−２１」までの時系列データ３０８が含まれており、港到着履歴データＰ１の到着時刻５０７が「２０１２−０５−０１」とする。

この場合、航海データにおける時系列データ３０８は、受信時刻４０１が「２０１２−０２−０１」から「２０１４−０７−２１」までの時系列データ３０８から、受信時刻４０１が到着時刻５０７「２０１２−０５−０１」以降の時系列データ３０８が削除された時系列データ３０８である。すなわち、航海データにおける時系列データ３０８の受信時刻４０１は、「２０１２−０２−０１」から「２０１２−０５−０１」までである。

このあと、特定装置１００は、未取得の港到着履歴データが存在するか否かを判断する（ステップＳ１１０７）。未取得の港到着履歴データ２１９が有る場合（ステップＳ１１０７：Ｙｅｓ）、ステップＳ１１０１に戻り、特定装置１００は、データメモリ２０８から未取得の１件の港到着履歴データＰ１を取得する（ステップＳ１１０１）。一方、未取得の港到着履歴データ２１９がない場合（ステップＳ１１０７：Ｎｏ）、特定装置１００は、データメモリ２０８にある港到着履歴データ２１９を削除し（ステップＳ１１０８）、航海データ取得処理（ステップＳ１００４）を終了する。

（最遠位置特定処理（ステップＳ１００５））
つぎに、最遠位置特定処理（ステップＳ１００５）について説明する。最遠位置特定処理（ステップＳ１００５）は、船舶が目的港に位置するときの時系列データ３０８が示す到着時刻と、一定期間（すなわち、予測期間ΔＴ）と、に基づいて、船舶が予測期間ΔＴ以内で対象海域２に進入できる時系列データのうち最古の時系列データの位置を、対象海域２からの船舶の最遠位置として特定する。具体的には、最遠位置特定処理（ステップＳ１００５）は、指定された船種において、航海データ取得処理（ステップＳ１００４）より得られた航海データについて、時系列データ３０８における受信時刻４０１の降順で対象海域に入った直後に受信された時系列データ３０８の受信時刻を特定し、その受信時刻より更に予測期間ΔＴ前の船舶位置を特定する処理である。

図１２は、最遠位置特定処理（ステップＳ１００５）の具体例を示す説明図である。×印は船舶の位置を示す時系列データ３０８である。（Ａ）は、航海データ取得処理（ステップＳ１００４）より得られた航海データ、すなわち、時系列データ３０８の集合である。ｔ１〜ｔ１０は受信時刻４０１であり、ｔ１０が港１２００の到着時刻である。

特定装置１００は、時刻ｔ１〜ｔ１０のうち最も直近となる時刻ｔ１０を基準として、時刻ｔ１〜ｔ１０から未選択の直近の２つの時刻の組み合わせ（ａ，ｂ）を選択する。ｂはａの１つ前の受信時刻４０１である。特定装置１００は、（ａ，ｂ）＝（ｔ１０，ｔ９）、（ｔ９，ｔ８）、（ｔ８，ｔ７）、…の順に選択し、（ａ，ｂ）が対象海域２の境界Ａｂをまたぐ組を選択する。

（Ａ）の例では、特定装置１００は、（ａ，ｂ）＝（ｔ７，ｔ６）の場合に、対象海域２の境界Ａｂをまたぐことになる。したがって、ｔ７の時系列データ３０８が対象海域２に最初に進入した船舶の位置を示すことになる。特定装置１００は、最初に対象海域２への進入が特定された時刻ｔ７での船舶位置を基準として、対象海域２からの船舶の最遠位置を特定することになる。

（Ｂ）および（Ｃ）は、（Ａ）の時刻ｔ７での船舶位置を基準時刻Ｔとして最遠位置特定例を示す。横軸は時間軸である。ここでは、予測期間ΔＴ＝１週間として説明する。特定装置１００は、時刻ｔ７以前の時刻ｔ１〜ｔ６から直近の２つの時刻の組み合わせ（ａ，ｂ）を選択する。ｂはａの次の時刻である。（Ｂ）では、（ａ，ｂ）＝（ｔ６，ｔ５）が選択される。特定装置１００は、基準時刻Ｔ（＝ｔ７）から予測期間前の時刻Ｔ−ΔＴを特定する。時刻Ｔ−ΔＴが（ａ，ｂ）の間の時刻となれば、ａの船舶位置を示す時系列データ３０８が最遠位置を示す最遠位置データとなる。

（Ｂ）に示した（ａ，ｂ）＝（ｔ６，ｔ５）の場合、ｔ５、ｔ６のいずれも時刻Ｔ−ΔＴ以降の時刻となるため、最遠位置が特定されない。この場合、特定装置１００は、時系列データ３０８の時刻を１つ過去にシフトして、（ａ，ｂ）＝（ｔ５，ｔ４）を選択する。このように、時刻Ｔ−ΔＴが（ａ，ｂ）の間の時刻となるまで、特定装置１００はシフトを続ける。

（Ｃ）は、（ａ，ｂ）＝（ｔ３，ｔ２）の場合の例である。この場合、ｔ３が時刻Ｔ−ΔＴ以前の時刻であり、ｔ２が時刻Ｔ−ΔＴ以降の時刻となるため、時刻ｔ３の時系列データ３０８が、予測期間ΔＴ内において対象海域２からの船舶の最遠位置を示す最遠位置データとなる。

図１３は、最遠位置特定処理（ステップＳ１００５）の詳細な処理手順例を示すフローチャートである。最遠位置特定処理（ステップＳ１００５）は、図１２に示したように、予測期間ΔＴでの対象海域２からの船舶の最遠位置を特定する処理である。

まず、特定装置１００は、１件の新しい「航行枠」種類の海域データｎｅｗＡ１を作成する（ステップＳ１３００）。「航行枠」種類の海域データｎｅｗＡ１は、海域番号６００〜予測期間６０７が設定された海域データ２２０である。種類６０１は「航行枠」である。

つぎに、特定装置１００は、データメモリ２０８から１件の未取得のＡＩＳデータａｉｓを取得し（ステップＳ１３０１）、ＡＩＳデータａｉｓの時系列データ３０８を受信時刻４０１の降順でソートする（ステップＳ１３０２）。そして、特定装置１００は、図１０のステップＳ１００３で取得された「対象海域」種類の海域データＡ１に進入した直後の時刻Ｔを時刻Ｔ０で初期化する（ステップＳ１３０３）。たとえば、特定装置１００は、時刻Ｔを時刻Ｔ０＝「１９００−０１−０１」に設定する。

つぎに、特定装置１００は、ＡＩＳデータａｉｓにおけるｉ件目の時系列データａを取得する（ステップＳ１３０４）、たとえば、特定装置１００は、未取得の時系列データ３０８のうち、受信時刻４０１が最新の時刻を、ＡＩＳデータａｉｓにおけるｉ件目の時系列データａとして取得する。

また、特定装置１００は、時系列データａの１つ前の時刻の時系列データ３０８ｂを取得する（ステップＳ１３０５）。これにより、（ａ，ｂ）の組が特定される。つぎに、特定装置１００は、海域データＡ１に進入した直後の時刻ＴがＴ＝Ｔ０であるか否かを判断する（ステップＳ１３０６）。すなわち、特定装置１００は、時刻Ｔが初期化した時刻Ｔ０のままであるか否かを判断する。

時刻Ｔが更新されず、Ｔ＝Ｔ０である場合（ステップＳ１３０６：Ｙｅｓ）、特定装置１００は、ａの経度４０２および緯度４０３が海域データＡ１の海域範囲６０５内であり、かつ、ｂの経度４０２および緯度４０３が海域データＡ１の海域範囲６０５外にあるか否かを判断する（ステップＳ１３０７）。すなわち、特定装置１００は、図１２に示した（Ａ）の状態であるか否かを判断する。

ａの経緯度が海域データＡ１の海域範囲６０５内でかつｂの経緯度が海域データＡ１の海域範囲６０５外である場合（ステップＳ１３０７：Ｙｅｓ）、特定装置１００は、時刻Ｔにａの受信時刻４０１を設定する（ステップＳ１３０８）。一方、ａの経緯度が海域データＡ１の海域範囲６０５外またはｂの経緯度が海域データＡ１の海域範囲６０５内である場合（ステップＳ１３０７：Ｎｏ）、ステップＳ１３１１に移行する。

また、ステップＳ１３０６において、Ｔ≠Ｔ０である場合（ステップＳ１３０６：Ｎｏ）、すなわち、時刻Ｔが時系列データａの受信時刻４０１に更新された場合、特定装置１００は、ａの受信時刻４０１が（Ｔ−ΔＴ）以降で、かつ、ｂの受信時刻４０１が（Ｔ−ΔＴ）以前であるか否かを判断する（ステップＳ１３０９）。すなわち、特定装置１００は、図１２の（Ｃ）の状態になるか否かを判断する。

ａが（Ｔ−ΔＴ）以降で、かつ、ｂが（Ｔ−ΔＴ）以前に受信された時系列データ３０８である場合（ステップＳ１３０９：Ｙｅｓ）、図１２の（Ｃ）の状態となる。したがって、特定装置１００は、ａを基に海域データｎｅｗＡ１に１件の新しい最遠位置データｆ１を追加し（ステップＳ１３１０）、ステップＳ１３１２に移行する。

なお、１件の新しい最遠位置データｆ１では、ａの受信時刻４０１が最遠位置データｆ１の受信時刻８０１、ａの経度４０２が最遠位置データｆ１の経度８０２、ａの緯度４０３が最遠位置データｆ１の緯度８０３に設定される。例外フラグ８０４は未設定である。

また、ａが（Ｔ−ΔＴ）以前、または、ｂが（Ｔ−ΔＴ）以降に受信された時系列データ３０８である場合（ステップＳ１３０９：Ｎｏ）、図１２の（Ｂ）の状態であるため、ステップＳ１３１１に移行する。

ステップＳ１３１１では、特定装置１００は、ＡＩＳデータａｉｓの未取得の時系列データ３０８が２件以上あるか否かを判断する（ステップＳ１３１１）。２件以上ある場合（ステップＳ１３１１：Ｙｅｓ）、（ａ、ｂ）の組が選択可能であるため、ステップＳ１３０４に戻る。一方、２件以上残存していない場合（ステップＳ１３１１：Ｎｏ）、ステップＳ１３１２に移行する。

ステップＳ１３１２では、特定装置１００は、未取得の時系列データ３０８があるか否かを判断する（ステップＳ１３１２）。未取得の時系列データ３０８がある場合（ステップＳ１３１２：Ｙｅｓ）、ステップＳ１３０１に戻る。一方、未取得の時系列データ３０８がない場合（ステップＳ１３１２：Ｎｏ）、特定装置１００は、データメモリ２０８にあるＡＩＳデータを削除して（ステップＳ１３１３）、最遠位置特定処理（ステップＳ１００５）を終了する。

（例外除去処理（ステップＳ１００６））
例外除去処理（ステップＳ１００６）は、複数の最遠位置データのうち、最近傍の最遠位置データとの距離が所定距離以上となる特定の最遠位置データを、航行枠５の海域範囲となる図形の作成元からの除外対象に設定する。具体的には、例外除去処理（ステップＳ１００６）は、指定された船種において、最遠位置特定処理（ステップＳ１００６）より得られた最遠位置データにおいて、各最遠位置データおよびその最近傍である最遠位置データにおける距離を算出し、距離が一定閾値以上の最遠位置データを例外データとして除去できるように設定する処理である。ここでは、例外データとなる最遠位置データ６０８には、「Ｔｒｕｅ」というフラグが設定される。例外データとなった最遠位置データ６０８は、当該最遠位置データ６０８をユーザ操作により画面から除去可能である。

図１４は、例外除去処理（ステップＳ１００６）の詳細な処理手順例を示すフローチャートである。まず、特定装置１００は、最遠位置特定処理（ステップＳ１００５）で作成された「航行枠」種類の海域データｎｅｗＡ１における未取得の最遠位置データｆ１を取得する（ステップＳ１４００）。

つぎに、特定装置１００は、既存の最近傍手法を用いて、「航行枠」種類の海域データｎｅｗＡ１の最遠位置データ４０８の中から、最遠位置データｆ１との空間距離が最も近い最遠位置データｎｎｆ１を特定する（ステップＳ１４０１）。

そして、特定装置１００は、最遠位置データｆ１と最遠位置データｎｎｆ１との空間距離をｄｉｓｔとし、空間距離ｄｉｓｔが一定閾値ΔＤＩＳＴ以上であるか否かを判断する（ステップＳ１４０２）。空間距離ｄｉｓｔがΔＤＩＳＴ以上である場合（ステップＳ１４０２：Ｙｅｓ）、特定装置１００は、最遠位置データｆ１の例外フラグ８０５に「Ｔｒｕｅ」を設定して（ステップＳ１４０４）、ステップＳ１４０５に移行する。「Ｔｒｕｅ」は、当該最遠位置データが、航行枠５から除去してもよい例外的な最遠位置データであることを示す。

例えば、ΔＤＩＳＴ＝２０００［ｋｍ］に設定されたとすると、最近傍である最遠位置データｎｎｆ１との空間距離ｄｉｓｔが２０００［ｋｍ］以上の最遠位置データｆ１は例外データとして除去設定する。

一方、空間距離ｄｉｓｔがΔＤＩＳＴ未満である場合（ステップＳ１４０２：Ｎｏ）、特定装置１００は、最遠位置データｆ１の例外フラグ８０５に「Ｆａｌｓｅ」を設定して（ステップＳ１４０５）、ステップＳ１４０５に移行する。「Ｆａｌｓｅ」は、当該最遠位置データが、例外データではない最遠位置データであることを示す。

つぎに、ステップＳ１４０５において、特定装置１００は、「航行枠」種類の海域データｎｅｗＡ１に未取得の最遠位置データ４０８があるか否かを判断する（ステップＳ１４０５）。未取得の最遠位置データ４０８が有る場合（ステップＳ１４０６：Ｙｅｓ）、ステップＳ１４００に戻り、「航行枠」種類の海域データｎｅｗＡ１から未取得の１件の最遠位置データｆ１を取得する（ステップＳ１４００）。一方、未取得の最遠位置データ４０８がない場合（ステップＳ１４０６：Ｎｏ）、例外除去処理（ステップＳ１００６）を終了する。

（航行枠生成処理（ステップＳ１００７））
航行枠生成処理（ステップＳ１００７）は、船舶ごとの最遠位置を示す複数の最遠位置データにより対象海域２を包含する図形を生成し、複数の最遠位置データと図形と対象海域２を指定する情報と予測期間ΔＴを含む海域データ（図７の「航行枠」種類の海域データ２２０）を生成する。具体的には、航行枠生成処理（ステップＳ１００７）は、指定された船種において、各最遠位置データにおいて一定の空間範囲（例えば、円範囲）に拡張し、それらの空間範囲をマージした新しい空間範囲を基に航行枠５の海域範囲を作成する処理である。一定の空間範囲は、最遠位置の誤差を許容する範囲である。また、例外除去処理（ステップＳ１００７）より例外データが設定される場合には、当該例外データを除外した最遠位置データが航行枠５の海域範囲の作成対象となる。

図１５は、航行枠生成処理（ステップＳ１００７）の詳細な処理手順例を示すフローチャートである。まず、特定装置１００は、「航行枠」種類の海域データｎｅｗＡ１の海域範囲６０５に海域範囲ＡｒｅａとしてＮＵＬＬを設定して初期化する（ステップＳ１５００）。つぎに、特定装置１００は、「航行枠」種類の海域データｎｅｗＡ１において航行枠生成処理（ステップＳ１００７）で未取得の最遠位置データｆ１を取得する（ステップＳ１５０１）。

そして、特定装置１００は、最遠位置データｆ１の例外フラグ８０５が「Ｆａｌｓｅ」であるか否かを判断する（ステップＳ１５０２）。最遠位置データｆ１の例外フラグ８０５に「Ｔｒｕｅ」が設定されている場合（ステップＳ１５０２：Ｎｏ）、ステップＳ１５０５に移行する。一方、最遠位置データｆ１の例外フラグ８０５に「Ｆａｌｓｅ」が設定されている場合（ステップＳ１５０２：Ｙｅｓ）、特定装置１００は、最遠位置データｆ１の経度８０２および緯度８０３で特定される位置を中心にして一定半径ｒの円範囲ｃｉｒｃｌｅを作成する（ステップＳ１５０３）。なお、円に限らず多角形でもよい。また、ステップＳ１５０３を省略してもよい。

つぎに、特定装置１００は、円範囲ｃｉｒｃｌｅと海域範囲Ａｒｅａをマージした新しい海域範囲Ａｒｅａを海域範囲６０５に設定する（ステップＳ１５０４）。つぎに、特定装置１００は、「航行枠」種類の海域データｎｅｗＡ１に未取得の最遠位置データ４０８があるか否かを判断する（ステップＳ１５０５）。未取得の最遠位置データ４０８が有る場合（ステップＳ１５０５：Ｙｅｓ）、ステップＳ１５０１に戻り、特定装置１００は、「航行枠」種類の海域データｎｅｗＡ１から未取得の１件の最遠位置データｆ１を取得する（ステップＳ１５０１）。

一方、未取得の最遠位置データ４０８がない場合（ステップＳ１５０５：Ｎｏ）、特定装置１００は、既存の多角形生成方法を用いて海域範囲Ａｒｅａを包含する図形を作成する（ステップＳ１５０６）。たとえば、特定装置１００は、海域範囲Ａｒｅａが内接するような図形を作成する。作成した図形が航行枠５である。そして、特定装置１００は、作成した図形（航行枠５）を「航行枠」種類の海域データｎｅｗＡ１の海域範囲６０５に設定し（ステップＳ１５０７）、航行枠生成処理を終了する。

（到着船舶予測処理（ステップＳ１０１３））
到着船舶予測処理（ステップＳ１０１３）は、指定された船種において、現時点まで受信された最新の時系列データ３０８を基づいて、各時系列データ３０８において航行枠生成処理（ステップＳ１００７）より得られた航行枠５との位置関係を確認し、予測期間以内に対象海域に到着する船舶を予測する処理である。時系列データ３０８の座標値は、該当航行枠５の中にあり、また、予測期間ΔＴがより短い航行枠５および対象海域２外にある船舶も出力対象となる。

図１６は、到着船舶予測処理（ステップＳ１０１３）の詳細な処理手順例を示すフローチャートである。まず、特定装置１００は、予測日付ｅＤａｙを設定する（ステップＳ１６００）。本フローチャートがバッチ処理により実行される場合、予測日付ｅＤａｙは、バッチ処理の実行日となる。なお、バッチ処理のタイミングに限らず、あらかじめ設定されたタイミングやユーザが操作したタイミングで到着船舶予測処理（ステップＳ１０１３）が実行されてもよい。この場合、当該タイミングの日付が予測日付ｅＤａｙに設定される。

つぎに、特定装置１００は、図１０のステップＳ１０００においてデータメモリ２０８に格納された「対象海域」種類の海域データ群から、未選択の海域データＡ１を選択する（ステップＳ１６０１）。そして、特定装置１００は、未選択の船舶種類ＳＫを選択する（ステップＳ１６０２）。ステップＳ１６０１〜ステップＳ１６０３で設定または選択された情報を「選択条件」とする。このあと、特定装置１００は、選択条件での到着船舶予測処理を実行する（ステップＳ１６０３）。これにより、到着船舶予測処理（ステップＳ１０１３）を終了する。

（選択条件での到着船舶予測処理（ステップＳ１６０３））
図１７は、選択条件での到着船舶予測処理（ステップＳ１６０３）の詳細な処理手順例を示すフローチャートである。まず、特定装置１００は、海域データＤＢ２０２から、選択条件に該当する「航行枠」種類の海域データ２２０をデータメモリ２０８に格納する（ステップＳ１７００）。選択条件に該当する「航行枠」種類の海域データ２２０とは、船種３０３および予測日付ｅＤａｙが属する季節４０６が、海域データＡ１の船種３０３および季節４０６と一致し、かつ、対象海域番号６０２が海域データＡ１の海域番号６００と一致する海域データである。

たとえば、ＶＬＣＣ船において、予測日付ｅＤａｙ＝「２０１４−０７−２１」に設定されたとすると、予測期間ΔＴ（たとえば、１週間）以内に東南アジア海域に到着する船舶を予測する必要がある。このため、特定装置１００は、船種３０３が「ＶＬＣＣ船」、季節４０６が「夏」、対象海域番号６０２が「東南アジア海域」の海域番号６００である「航行枠」種類の海域データ２２０を、海域データＤＢ３０２から取得することになる。

つぎに、特定装置１００は、ＡＩＳデータＤＢ２００から、全ＡＩＳデータ２１８をデータメモリ２０８に格納する（ステップＳ１７０１）。そして、特定装置１００は、ステップＳ１７００においてデータメモリ２０８に格納された「航行枠」種類の海域データ群から、１件の「航行枠」種類の海域データｎｅｗＡ１を取得する（ステップＳ１７０２）。

つぎに、特定装置１００は、データメモリ２０８から１件の未取得のＡＩＳデータａｉｓを取得し（ステップＳ１７０３）、ＡＩＳデータａｉｓにおける最新に受信された時系列データａを取得する（ステップＳ１７０４）。

そして、特定装置１００は、時系列データａの受信時刻４０１が予測日付ｅＤａｙと同じ日であるか否かを判断する（ステップＳ１７０５）。時系列データａの受信時刻４０１が予測日付ｅＤａｙと異なる日である場合（ステップＳ１７０５：Ｎｏ）、ステップＳ１７０７に移行する。一方、同じ日である場合（ステップＳ１７０５：Ｙｅｓ）、特定装置１００は、予測船舶データ生成処理を実行して（ステップＳ１７０６）、ステップＳ１７０７に移行する。予測船舶データ生成処理（ステップＳ１７０６）の詳細については、図１８および図１９で後述する。

つぎに、特定装置１００は、データメモリ２０８上に未取得のＡＩＳデータ２１８があるか否かを判断する（ステップＳ１７０７）。未処理のＡＩＳデータ２１８が有る場合（ステップＳ１７０７：Ｙｅｓ）、ステップＳ１７０４に戻り、データメモリ２０８から未取得の１件のＡＩＳデータ２１８を取得する（ステップＳ１７０３）。

一方、未取得のＡＩＳデータ２１８がない場合（ステップＳ１７０７：Ｎｏ）、特定装置１００は、データメモリ２０８上に未取得の「航行枠」種類の海域データ２２０があるか否かを判断する（ステップＳ１７０８）。未取得の「航行枠」種類の海域データ２２０がある場合（ステップＳ１７０８：Ｙｅｓ）、ステプＳ１７０２に戻り、特定装置１００は、データメモリ２０８から未取得の１件の「航行枠」種類の海域データ２２０を取得する（ステップＳ１７０３）。

一方、未取得の「航行枠」種類の海域データ２２０がない場合（ステップＳ１７０８：Ｎｏ）、特定装置１００は、予測船舶データ生成処理（ステップＳ１７０６）によって生成された予測船舶データ２２１を予測船舶データＤＢ２０３に書き込む（ステップＳ１７０９）。そして、特定装置１００は、データメモリ２０８にあるＡＩＳデータ２１８を削除し（ステップＳ１７１０）、到着船舶予測処理を終了する。

（予測船舶データ生成処理（ステップＳ１７０６））
予測船舶データ生成処理（ステップＳ１７０６）は、予測期間ΔＴ以内に対象海域２に到着する船舶を予測する処理である。予測期間ΔＴが複数存在する場合は、予測期間ΔＴが重複しないように、予測期間ΔＴごとに予測期間ΔＴ以内に対象海域に到着する船舶を予測する。たとえば、予測期間ΔＴが１週間、２週間、３週間の３種類である場合、予測船舶データ生成処理（ステップＳ１７０６）は、予測日付から１週間以内で対象海域２に到着する船舶、予測日付から１週間経過後２週間以内で対象海域２に到着する船舶、予測日付から２週間経過後３週間以内で対象海域２に到着する船舶を予測する。

図１８は、予測船舶データ生成処理（ステップＳ１７０６）における予測期間ΔＴごとの海域範囲に含まれる船舶を特定する例を示す説明図である。（Ａ）は予測期間ΔＴ＝１週間、（Ｂ）は予測期間ΔＴ＝２週間、（Ｃ）は予測期間ΔＴ＝３週間の場合の例である。

（Ａ）において、時系列データ１８００は、対象海域２内の船舶位置を示し、時系列データ１８０１は、予測期間ΔＴ＝１週間の海域範囲５１（対象海域２は含まない）内の船舶位置を示す。したがって、予測期間ΔＴ＝１週間の場合、時系列データ１８００は、生成される予測船舶データには含まれないが、時系列データ１８０１は、生成される予測船舶データに含まれる。

（Ｂ）において、時系列データ１８０２は、予測期間ΔＴ＝２週間の海域範囲５２（対象海域２および海域範囲６１は含まない）内の船舶位置を示す。したがって、予測期間ΔＴ＝２週間の場合、時系列データ１８００，１８０１は、生成される予測船舶データには含まれないが、時系列データ１８０２は、生成される予測船舶データに含まれる。

（Ｃ）において、時系列データ１８０３は、予測期間ΔＴ＝３週間の海域範囲５３（対象海域２、海域範囲５１、５２は含まない）内の船舶位置を示す。したがって、予測期間ΔＴ＝３週間の場合、時系列データ１８００，１８０１，１８０２は、生成される予測船舶データには含まれないが、時系列データ１８０３は、生成される予測船舶データに含まれる。

図１９は、予測船舶データ生成処理（ステップＳ１７０６）の詳細な処理手順例を示すフローチャートである。まず、特定装置１００は、図１７のステップＳ１７０４で取得された時系列データａの経度４０２および緯度４０３が、図１６のステップＳ１６０１で選択された「対象海域」種類の海域データＡ１の海域範囲６０５内にあるか否かを判断する（ステップＳ１９００）。

時系列データａの経度４０２および緯度４０３が海域データＡ１の海域範囲６０５内にある場合（ステップＳ１９００：Ｙｅｓ）、時系列データａは、図１８の（Ａ）の時系列データ１８００と同じ状態であるため、予測船舶データ生成処理（ステップＳ１７０６）を終了する。一方、時系列データａの経度４０２または緯度４０３が、海域データＡ１の海域範囲６０５外にある場合（ステップＳ１９００：Ｎｏ）、時系列データａの経度４０２および緯度４０３が、図１７のステップＳ１７０２で取得された「航行枠」種類の海域データｎｅｗＡ１の海域範囲６０５内にあるか否かを判断する（ステップＳ１９０１）。

時系列データａの経度４０２または緯度４０３が「航行枠」種類の海域データｎｅｗＡ１の海域範囲６０５外にあった場合（ステップＳ１９０１：Ｎｏ）、予測船舶データ生成処理を終了する。たとえば、時系列データａが時系列データ１８００であり、海域範囲６０５が海域範囲５１である場合は、時系列データ１８００は海域範囲外となる。

一方、時系列データａの経度４０２および緯度４０３が「航行枠」種類の海域データｎｅｗＡ１の海域範囲６０５内にある場合（ステップＳ１９０１：Ｙｅｓ）、特定装置１００は、データメモリ２０８から予測期間６０７が「航行枠」種類の海域データｎｅｗＡ１の予測期間６０７より短い「航行枠」種類の海域データ２２０の一覧情報Ｌｉｓｔを取得する（ステップＳ１９０２）。

つぎに、特定装置１００は、取得した一覧情報Ｌｉｓｔから、１件の未取得の「航行枠」種類の海域データｎｅｗＡ２を取得する（ステップＳ１９０３）。そして、特定装置１００は、時系列データａの経度４０２および緯度４０３が海域データｎｅｗＡ２の海域範囲６０５内にあるか否かを判断する（ステップＳ１９０４）。

時系列データａの経度４０２および緯度４０３が海域データｎｅｗＡ２の海域範囲６０５内にある場合（ステップＳ１９０４：Ｙｅｓ）、予測船舶データ生成処理（ステップＳ１７０６）を終了する。たとえば、海域データｎｅｗＡ１の海域範囲６０５が海域範囲５１、海域データｎｅｗＡ２の海域範囲６０５が海域範囲５２であるとする。時系列データａが時系列データ１８０２である場合、時系列データ３０８は海域範囲５１外でかつ海域範囲５２内であるため、予測船舶データ生成処理（ステップＳ１７０６）を終了する。

一方、時系列データａの経度４０２および緯度４０３が海域データｎｅｗＡ２の海域範囲６０５外にある場合（ステップＳ１９０４：Ｎｏ）、一覧情報Ｌｉｓｔに未取得の「航行枠」種類の海域データ２２０があるか否かを判断する（ステップＳ１９０５）。未取得の「航行枠」種類の海域データ２２０がある場合（ステップＳ１９０５：Ｎｏ）、ステップＳ１９０３に戻り、特定装置１００は、一覧情報Ｌｉｓｔから未取得の１件の「航行枠」種類の海域データｎｅｗＡ２を取得する（ステップＳ１９０３）。

一方、未取得の「航行枠」種類の海域データ２２０がない場合（ステップ１２０５：Ｎｏ）、特定装置１００は、図１７のステップＳ１７０３で取得されたＡＩＳデータａｉｓと時系列データａとに基づいて、新規予測船舶データｓｈｉｐを生成し（ステップＳ１９０６）、予測船舶データ生成処理を終了する。

具体的には、たとえば、特定装置１００は、図１６のステップＳ１６０１で選択された海域データＡ１の海域番号６００を、新規予測船舶データｓｈｉｐの対象海域番号９０１に設定する。また、特定装置１００は、図１７のステップＳ１７０３で取得されたＡＩＳデータａｉｓの船舶番号３００を、新規予測船舶データｓｈｉｐの船舶番号９０２に設定する。

また、特定装置１００は、ＡＩＳデータａｉｓのＩＭＯ番号３０１を、新規予測船舶データｓｈｉｐのＩＭＯ番号９０３に設定する。また、特定装置１００は、ＡＩＳデータａｉｓの船舶名３０２を、新規予測船舶データｓｈｉｐの船舶名９０４に設定する。また、特定装置１００は、ＡＩＳデータａｉｓの船種３０３を、新規予測船舶データｓｈｉｐの船種９０５に設定する。

また、特定装置１００は、ＡＩＳデータａｉｓの最大積載量３０５を、新規予測船舶データｓｈｉｐの最大積載量９０６に設定する。また、特定装置１００は、ＡＩＳデータａｉｓの海運会社３０６を、新規予測船舶データｓｈｉｐの海運会社９０７に設定する。また、特定装置１００は、時系列データａの受信時刻４０１を、新規予測船舶データｓｈｉｐの受信時刻９０８に設定する。

また、特定装置１００は、時系列データａの経度４０２を、新規予測船舶データｓｈｉｐの経度９０９に設定する。また、特定装置１００は、時系列データａの緯度４０３を、新規予測船舶データｓｈｉｐの緯度９１０に設定する。また、特定装置１００は、時系列データａの速度４０４を、新規予測船舶データｓｈｉｐの速度９１１に設定する。

また、特定装置１００は、時系列データａの喫水４０５を、新規予測船舶データｓｈｉｐの喫水９１２に設定する。また、特定装置１００は、一覧情報Ｌｉｓｔにある予測期間６０７の最大値から「航行枠」種類の海域データｎｅｗＡ１の予測期間６０７までの期間を新規予測船舶データｓｈｉｐの必要期間９１３に設定する。また、特定装置１００は、図１６のステップＳ１６００で設定された予測日付ｅＤａｙを、新規予測船舶データｓｈｉｐの予測日付９１４に設定する。

＜航行枠表示処理部２１６による表示処理例＞
予測表示処理（Ｓ１０１４）は、航行枠生成処理（ステップＳ１００７）によって生成された「航行枠」種類の海域データ２２０を出力する。

図２０および図２１は、図１０の予測表示処理（Ｓ１０１４）における航行枠表示処理部２１６により表示される表示画面例を示す説明図である。ユーザは、航行枠生成処理部２１４より生成された「航行枠」種類の海域データ２２０が航行枠表示処理部２１６により表示された表示画面２０００から、確認または修正することができる。

図２０の表示画面２０００では、特定装置１００は、ユーザ操作により、プルダウンメニューである船種指定メニュー２００１から船種、予測期間２００２から予測期間、対象海域２００３から対象海域を選択する。例えば、それぞれ「ＶＬＣＣ船」と「１週間」と「東南アジア海域」が選択されたとする。そして、特定装置１００は、ユーザ操作により、航行枠季節２００４の表示季節を指定し、航行枠表示ボタン２００５を押下する。

これにより、特定装置１００は、海域データＤＢ２０２から東南アジア海域および東南アジア海域を対象海域とする「航行枠」種類の海域データ２２０を取得し、航行枠情報２００６を表示する。航行枠情報２００６には、例えば、対象海域２の海域範囲２００７と、最遠位置データ６０８より拡張した円範囲２００８と、航行枠５の海域範囲２００９と、例外データとして除外した最遠位置データ２０１０と、メニューボックス「例外に設定」または「対象に設定」を選択可能なメニュー２０１１が表示される。

図２１の表示画面２０００は、図２０において、ユーザ操作により、例えば、メニュー２０１１を用いて、ある最遠位置データ２０１２について「例外に設定」が選択された場合、特定装置１００は、航行枠５の海域範囲を再生成する。これにより、図２０の航行枠５の海域範囲２００９が航行枠５の海域範囲２１００に更新され、例外設定された最遠位置データ２０１２は、航行枠５の海域範囲２１００外となる。

＜予測表示処理部２１０による表示処理例＞
図２２〜図２４は、図１０の予測表示処理（Ｓ１０１４）における予測表示処理部２１０による表示画面例を示す説明図である。ユーザは、到着船舶予測処理部２１５より予測された予測船舶データ２２１が到着船舶予測表示処理部２１７により表示された表示画面２２００から、海運会社別に予測された到着船舶の数や最大積載量、喫水および船舶位置などを確認することができる。

図２２に示す表示画面２２００では、特定装置１００は、ユーザ操作により、プルダウンメニューである船種指定メニュー２２０１より船種、予測日付２２０２より予測日付、対象海域２２０３より対象海域を選択する。例えば、それぞれ「ＶＬＣＣ船」と「２０１４−０７−１１」と「東南アジア海域」が選択されたとする。そして、特定装置１００は、ユーザ操作により、航行枠季節２２０４の表示季節を指定し、予測ボタン２２０５を押下する。

これにより、特定装置１００は、予測船舶データＤＢから、「２０１４―０７−１１」に予測された東南アジア海域に到着可能なＶＬＣＣ船を示す予測船舶データ２２１を取得して、予測船舶結果２２０６を表示する。予測船舶結果２２０６には、例えば、東南アジア海域に到着するまでの必要期間２２０７と、対象海運会社を選択するチェックボックス２２０８と、海運会社２２０９と、海運会社別に必要期間２２０７に東南アジアに到着する合計船舶数２２１０と、喫水１２［ｍ］未満の船舶数２２１１と、喫水１２［ｍ］以上の船舶数２２１２と、詳細船舶一覧ボタン２２１３と、位置表示ボタン２２１４と、が表示される。

特定装置１００は、ユーザ操作により、チェックボックス２２０８の複数の海運会社を選択する。例えば、Ａ社とＢ社とＣ社が選択されたとする。特定装置１００は、ユーザ操作により、詳細船舶一覧ボタン２２１３を押下する。これにより、図２３に示すように、表示画面２２００に、選択された海運会社Ａ〜Ｃの予測船舶一覧２３００が表示される。予測船舶一覧２３００には、例えば、東南アジア海域に到着するまでの必要期間２２０７と、対象船舶を選択するチェックボックス２３０１と、海運会社２３０２と、ＩＭＯ番号２３０３と、船舶名２３０４と、最大積載量２３０５と、喫水２３０６と、速度２３０７と、位置表示ボタン２２１４と、が表示される。

特定装置１００は、ユーザ操作により、チェックボックス２３０１の複数の船舶を選択する。例えば、Ａ社の上位５隻の船舶が選択されたとする。特定装置１００は、ユーザ操作により、位置表示ボタン２２１４を押下する。これにより、図２４に示すように、表示画面２２００に、選択された船舶の現在位置が地図上にマッピングされた予測船舶位置情報２４００が表示される。予測船舶位置情報２４００には、例えば、東南アジア海域に到着するまでの必要期間２２０７と、対象海域２の海域範囲２００７と、航行枠５の海域範囲２１００と、船舶位置を記す船舶マック２４０１と、マウスが指定した船舶において船舶情報を表示するボックス２４０２と、が表示される。

本実施例によれば、対象海域２からの船舶の最遠位置を示す最遠位置データ６０８を特定することができるため、ユーザ（たとえば、荷主）は、対象海域への到達以前の指定された時刻から一定期間内に対象海域に到着する船舶を予測することができる。したがって、ユーザにとっては、交渉可能な船舶や海運会社といった選択肢が多くなり、用船契約の交渉に役立てることができる。

また、過去の情報を用いることにより、対象海域２までの単純な直線距離ではなく、水陸分布や船舶サイズによって船舶の現在地から対象海域２に行くまでの迂回航路を考慮した航行枠５を生成することができる。すなわち、船舶の中には、対象海域２までの直線距離が短いとしても、実際に航行する距離が長くなる船舶もある。したがって、本実施例によれば、予測期間以内に対象海域２に到着可能な船舶が列挙されるため、対象海域２への到着船舶を高精度に予測することができる。

また、例外除去処理をすることにより、特定の最遠位置データを航行枠５の生成元から除外することができるため、航行枠５の生成の自由度の向上を図ることができる。

なお、上述した実施例は、図２０〜図２４に示したように、予測表示処理部２１０によって表示画面２０００、２２００を表示する例について説明したが、本実施例の特定装置１００が、表示画面２０００、２２００を表示させるための情報を外部装置に送信し、外部装置が、表示画面２０００、２２００を表示することとしてもよい。

なお、本実施例では、ＡＩＳデータ２１８と港到着履歴データ２１９とを基づいて作成された航行枠５、および、航行枠５を基に予測された到着船舶の特徴についての予測結果を確認する図面例を示した。この他にも、特定装置１００は、船舶の航行目的地と組み合わせて表示してもよい。また、特定装置１００は、前日（または一定期間だけ前の）からの船舶の位置変化および位置変化を基に対象海域に向かっているかどうかを判別した結果と組み合わせて表示してもよい。また、特定装置１００は、用船仲介者が提供する成約情報や予定到着船舶情報や需要情報などと組み合わせて表示したり、静的データや動的データなどと組み合わせて表示してもよい。例えば、将来の１週間に対象海域２に到着するフリー船を確認するため、まず１週間後に到着する船舶を予測し、予測された船舶は用船仲介者が提供する成約情報と照合することで、予測された到着船舶中の成約済船およびフリー船を判別することができる。

以上説明したように、本実施例によれば、特定装置１００は、ＡＩＳデータ２１８と港到着履歴データ２１９を基づいて、予測期間ΔＴ以内に対象海域２へ到着する航行枠５を作成し、航行枠５および予測時点の船舶位置における関係を確認し、一定期間以内に対象海域に到着する船舶を予測し、算出された航行枠および予測された到着船舶を可視化する。したがって、用船契約の効率化を図ることができる。

なお、本発明は前述した実施例に限定されるものではなく、添付した特許請求の範囲の趣旨内における様々な変形例及び同等の構成が含まれる。例えば、前述した実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに本発明は限定されない。また、ある実施例の構成の一部を他の実施例の構成に置き換えてもよい。また、ある実施例の構成に他の実施例の構成を加えてもよい。また、各実施例の構成の一部について、他の構成の追加・削除・置換をしてもよい。

また、前述した各構成、機能、処理部、処理手段等は、それらの一部又は全部を、例えば集積回路で設計する等により、ハードウェアで実現してもよく、プロセッサがそれぞれの機能を実現するプログラムを解釈し実行することにより、ソフトウェアで実現してもよい。

各機能を実現するプログラム、テーブル、ファイル等の情報は、メモリ、ハードディスク、ＳＳＤ（Ｓｏｌｉｄ Ｓｔａｔｅ Ｄｒｉｖｅ）等の記憶装置、又は、ＩＣカード、ＳＤカード、ＤＶＤ等の記録媒体に格納することができる。

また、制御線や情報線は説明上必要と考えられるものを示しており、実装上必要な全ての制御線や情報線を示しているとは限らない。実際には、ほとんど全ての構成が相互に接続されていると考えてよい。

１００ 特定装置
２０９ 予測処理部
２１０ 予測表示処理部
２１１ 航海データ取得処理部
２１２ 最遠位置特定処理部
２１３ 例外除去処理部
２１４ 航行枠生成処理部
２１５ 到着船舶予測処理部
２１６ 航行枠表示処理部
２１７ 到着船舶予測表示処理部

Claims (7)

1. プログラムを実行するプロセッサと、前記プロセッサが実行するプログラムを格納するメモリと、を備える特定装置であって、
   前記プロセッサは、
   対象海域内の目的港へ航行する船舶の前記目的港に到着するまでの時系列な位置および時刻を示す一連の時系列データを取得する取得処理と、
   前記取得処理によって取得された一連の時系列データの中の連続する２つの時系列データと、前記一連の時系列データのうち前記船舶が前記目的港に位置するときの時系列データが示す到着時刻と、予測期間と、に基づいて、前記一連の時系列データのうち前記船舶が前記予測期間以内で前記対象海域に進入できる時系列データの中の最古の時系列データの位置を、前記対象海域からの前記船舶の最遠位置として特定する特定処理と、
   を実行することを特徴とする特定装置。
2. 前記取得処理では、前記プロセッサは、船舶ごとに、前記一連の時系列データを取得し、
   前記特定処理では、前記プロセッサは、前記船舶ごとに、前記船舶の最遠位置を特定し、
   前記プロセッサは、
   前記船舶ごとの最遠位置を示す複数の最遠位置データにより前記対象海域を包含する図形を生成し、前記複数の最遠位置データと前記図形と前記対象海域を示す情報と前記予測期間とを含む海域データを生成する海域データ生成処理を実行することを特徴とする請求項１に記載の特定装置。
3. 前記プロセッサは、
   前記複数の最遠位置データのうち、最近傍の最遠位置データとの距離が所定距離以上となる特定の最遠位置データを、前記図形の生成元からの除外対象に設定する設定処理を実行することを特徴とする請求項２に記載の特定装置。
4. 前記プロセッサは、
   前記海域データ生成処理によって生成された海域データを出力する出力処理を実行することを特徴とする請求項２または３に記載の特定装置。
5. 前記プロセッサは、
   前記海域データ生成処理によって生成された海域データを出力する出力処理を実行し、
   前記海域データ生成処理では、前記プロセッサは、前記出力処理による前記海域データの出力後に、外部から前記特定の最遠位置データの除外指示が与えられた場合、前記複数の最遠位置データのうち前記特定の最遠位置データを除いた最遠位置データ群を用いて前記図形を更新し、
   前記出力処理では、前記プロセッサは、更新後の図形を出力することを特徴とする請求項３に記載の特定装置。
6. 前記特定装置は、
   前記海域データ生成処理によって生成された海域データを記憶するデータベースを有し、
   前記プロセッサは、
   時系列な位置および時刻を示す予測対象船舶に関する一連の時系列データの中のある予測日付に該当する特定の時系列データの位置を包含する図形を有する第１の海域データを前記データベースから特定し、前記第１の海域データに含まれる第１の予測期間よりも短い第２の予測期間を含む第２の海域データを前記データベースから特定し、前記特定の時系列データの位置が前記第２の海域データに包含されていない場合、前記予測対象船舶と特定する情報と、前記特定の時系列データと、前記第１の予測期間と、を含む予測船舶データを生成する予測船舶データ生成処理を実行することを特徴とする請求項２に記載の特定装置。
7. プログラムを実行するプロセッサと、前記プロセッサが実行するプログラムを格納するメモリと、を備える特定装置による特定方法であって、
   前記プロセッサは、
   対象海域内の目的港へ航行する船舶の前記目的港に到着するまでの時系列な位置および時刻を示す一連の時系列データを取得する取得処理と、
   前記取得処理によって取得された一連の時系列データの中の連続する２つの時系列データと、前記一連の時系列データのうち前記船舶が前記目的港に位置するときの時系列データが示す到着時刻と、予測期間と、に基づいて、前記一連の時系列データのうち前記船舶が前記予測期間以内で前記対象海域に進入できる時系列データの中の最古の時系列データの位置を、前記対象海域からの前記船舶の最遠位置として特定する特定処理と、
   を実行することを特徴とする特定方法。