(57)【特許請求の範囲】
【請求項１】
少なくとも１つのＡＩＳチャンネルにて衛星によって受信された信号を示すデジタルデータを提供するステップと、
少なくとも１つのフラグメント状疑似メッセージがＡＩＳ送信機のデータベース内で識別された送信機に固有である１つ以上のスタティックフィールドからのビットから成る、少なくとも１つのフラグメント状疑似メッセージを含むステップと、
デジタルデータと少なくとも１つの前記フラグメント状疑似メッセージとを相関化することにより、デジタルデータ内の候補メッセージを識別するステップ（Ｅ１８）とを含む、ＡＩＳメッセージを検出するための方法において、
前記少なくとも１つのフラグメント状疑似メッセージが、送信時間および／または送信機の位置に依存すると共にＡＩＳ送信機のデータベース（２６）に基づいて決定される少なくとも１つのＡＩＳメッセージフラグメントを含むことを特徴とする、ＡＩＳメッセージを検出するための方法。
【請求項２】
前記フラグメント状疑似メッセージが、送信位置に応じて決まる少なくとも１つのＡＩＳメッセージフラグメントを含み、前記フラグメント状疑似メッセージを合成することは、
前記信号の受信時に衛星によるサービスを受けるゾーン内に位置するＡＩＳ送信機を、ＡＩＳ送信機のデータベースに基づいて識別するステップ（Ｅ１４）と、
前記識別された送信機の、前記データベースに記憶されていた位置に基づき、送信機の位置に応じた少なくとも１つのＡＩＳメッセージフラグメントを決定するステップとを含むことを特徴とする、請求項１に記載の方法。
【請求項３】
前記送信機の位置に応じて決まる少なくとも１つのＡＩＳメッセージフラグメントは、「経度」フィールドからの最大位ビットおよび／または「緯度」フィールドからの最大位ビットを含むことを特徴とする、請求項２に記載の方法。
【請求項４】
前記デジタルデータにはタイムスタンプが押されており、前記フラグメント状疑似メッセージは、送信時間に応じて決まる１つ以上のＡＩＳメッセージフラグメントを含み、前記フラグメント状疑似メッセージを合成するステップは、前記デジタルデータのタイムスタンプに基づき、送信時間に応じて決まる少なくとも１つのＡＩＳメッセージフラグメントを決定することを含む、請求項１〜３のうちのいずれか１項に記載の方法。
【請求項５】
送信時間に依存するする少なくとも１つのＡＩＳメッセージフラグメントは、「タイムスタンプ」フィールドからのビット、または「ＵＴＣ年」、「ＵＴＣ月」、「ＵＴＣ日」、「ＵＴＣ時」、「ＵＴＣ分」および／または「ＵＴＣ秒」フィールドからのビットを含むことを特徴とする、請求項４に記載の方法。
【請求項６】
前記フラグメント状疑似メッセージは、前記データベース内で識別された前記送信機に固有の１つ以上のスタティックフィールド、例えば「ユーザーＩＤ」フィールド、「ＡＩＳバージョンインジケータ」フィールド、「ＩＭＯ番号」フィールド、「コールサイン」フィールド、「ネーム」フィールド、「船舶のタイプおよび積荷のタイプ」フィールド、「位置に対する全寸法／基準」フィールド、「電子位置固定デバイスのタイプ」フィールドからのビットを含むことを特徴とする、請求項１〜５のうちのいずれか１項に記載の方法。
【請求項７】
前記フラグメント状疑似メッセージは、前記データベース内で識別された前記送信機に固有の１つ以上のアクティーブフィールド、例えば「ＲＡＩＭフラグ」フィールド、「特殊操縦インジケータ」フィールド、「推定到着時間」フィールド、「現在のドロート」フィールド、「目的地」フィールド、「ナビゲーションステータス」フィールド、「地上のコース」フィールド、「地上の速度」フィールドおよび／または「通信ステータス」フィールドからのビットを含むことを特徴とする、請求項１〜６のうちのいずれか１項に記載の方法。
【請求項８】
デジタルデータ内の候補メッセージが一旦識別された場合、前記候補メッセージをデコードすること（Ｅ２０）、および前記デコードされたメッセージ内に含まれる情報を使ってＡＩＳ送信機の前記データベースを更新すること（Ｅ２２）を特徴とする、請求項１〜７のうちのいずれか１項に記載の方法。
【請求項９】
前記デジタルデータと複数のフラグメント状疑似メッセージとを相関化することにより、前記候補メッセージ（Ｅ１８）を識別し、種々の相関値をシリアルおよび／またはパラレルに計算することを特徴とする、請求項１〜８のうちのいずれか１項に記載の方法。
【請求項１０】
前記フラグメント状疑似メッセージは、前記候補メッセージ内のスタッフィングビットが存在するか、存在しないかを検査するようになっていることを特徴とする、請求項９に記載の方法。
【請求項１１】
前記デジタルデータ内の１つ以上のＡＩＳメッセージを識別した送信機に対する三角測量により、前記信号受信時の前記衛星の位置を推定する（Ｅ１２”）ことを特徴とする、請求項１〜１０のうちのいずれか１項に記載の方法。
【請求項１２】
前記デジタルデータと前記フラグメント状疑似メッセージとを相関化するために、前記フラグメント状疑似メッセージのスタートと前記ＡＩＳ時間フレームのうちの時間インターバルのスタートとを同期化することを特徴とする、請求項１〜１１のうちのいずれか１項に記載の方法。
【請求項１３】
前記デジタルデータと前記フラグメント状疑似メッセージとを相関化するために、前記フラグメント状疑似メッセージに対する計算されたドップラーシフトを送信機位置の関数として補償することを特徴とする、請求項１〜１２のうちのいずれか１項に記載の方法。
【請求項１４】
前記ドップラーシフトを補償するために、前記データベースに記憶されている送信機のキャリア周波数の誤差を考慮する（２６）ことを特徴とする、請求項１３に記載の方法。
【請求項１５】
前記デジタルデータと前記フラグメント状疑似メッセージとを相関化するために、前記データベースに記憶されている前記送信機の変調指数を考慮する（２６）ことを特徴とする、請求項１〜１４のうちのいずれか１項に記載の方法。
【請求項１６】
コンピュータによる実行時に前記請求項１〜１５のうちのいずれか１項に記載の方法を、前記コンピュータに実行させるコンピュータで実行可能な命令を含む、コンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、船舶およびブイのための衛星監視システムに関する。
【背景技術】
【０００２】
自動識別システム（ＡＩＳ）とは、２０〜３０海里（３７〜５６ｋｍ）の距離内で船舶と沿岸局とナビゲーションサポートシステムとの間でショートメッセージを交換するためのシステムのことである。ＡＩＳは、ＶＨＦの船舶モバイルバンドにおいて自己組織型時間分割マルチアクセス（ＳＯＴＤＭＡ）を使用している。（１６２.９７５ＭＨｚを中心とする）ＡＩＳ  １チャンネルまたは（１６２.０２５ＭＨｚを中心とする）ＡＩＳ  ２チャンネルで通信が行われ、各チャンネルは、２５ｋＨｚの幅を有し、交換されるメッセージは、所定の長さを有し、船舶の識別（ＩＤ）、船舶の位置、船舶の航路および船舶の速度のようなデータを含む。
【０００３】
各通信チャンネルは、６０秒の長さのフレーム内で時間分割されており、各フレームは、２２５０個のスロットに分割されており、通常の運用では１つのＡＩＳメッセージ（またはレポート）が各スロットを占有できるようになっている。通常のモードでは、各船舶は２つの通信チャンネルで自己のＡＩＳメッセージを交互に送信する。所定の船舶からのＡＩＳメッセージの頻度は、その船舶の速度に応じて決まる。その船舶が航行中であれば、２〜１０秒ごとに１つのメッセージが送信され、船舶の速度が速くなればなるほど、連続する２つのメッセージの間のスペースは短くなる。船舶が停止している場合には、３分ごとにメッセージが送信される。このシステムに参加している異なるプラットフォームからのメッセージの間のインターバルは、例えばＩＴＵ−Ｒ  Ｍ.２０８４レポート「自動識別システムメッセージの衛星検出」に記載されている。
【０００４】
船舶および沿岸局のＡＩＳ受信機は、近くにあるすべての船舶からのメッセージを受信し、これら船舶がローカル環境の詳細マップを作成できるようにしている。
【０００５】
このＡＩＳは、船舶の時間分割アクセスを管理する中央ＴＤＭＡコントローラがなくても運用される。各ＡＩＳトランスポンダは、通信チャンネルをモニタし、各メッセージ内で次のメッセージを送信するために使用されるスロットを留保している。ローカルの地理的ゾーンでは、衝突解決プロトコルが同一スロットにおけるメッセージの衝突を防止している。
【０００６】
実情は、ＡＩＳメッセージを衛星で検出することではない。実際には低軌道を回る衛星のアンテナのフットプリントが、共にコーディネートされていない複数のローカルの地理的ゾーンを一般にカバーしている。従って、その衛星は同じスロットにおけるメッセージのランダムな衝突を観察するので、その結果、衝突するメッセージの１つ以上が失われることがあり得る。所定のＡＩＳメッセージの正しい検出がデコードされる確率は、衝突の数のレートが増加するにつれて低下する。この件に関するより詳細な情報は、ＩＴＵ−Ｒ  Ｍ.２０８４レポートに記載されている。本明細書も、衛星のキャパシティを改善するための可能な方法を示唆するものである。
〇衛星の動きの結果生じるドップラーシフトは、フットプリント内の送信機の位置に応じて決まるので、このドップラーシフトに基づき、複数のメッセージを弁別することができる。衛星の動きによって生じるドップラーシフトは、±４ｋＨｚに達し得る。それにもかかわらず、このアプローチではＡＩＳの規格は、船舶のＴＤＭＡ送信機に対し、キャリア周波数に対して±５００Ｈｚの誤差を許容していることを念頭に入れておかなければならない。
〇所定の船舶によって送信される連続する２つのＡＩＳメッセージの間には情報の重要な冗長性がある。ＩＴＵ−Ｒ  Ｍ.２０８４レポートは、地上をベースとする局で特殊な処理をすることを見込んで、特に衛星が１回の通過する時間長さに対応する時間スロットにおいて、先に受信された信号のスキャンされたコピーと着信中のＡＩＳ信号とを連続的に相関化するために、２つのＡＩＳチャンネルで受信されたデータをリアルタイムでダウンロードすることにより、連続送信することを示唆している。
【０００７】
米国特許第７，８７６，８６５号は、スペース内のＡＩＳメッセージを検出し、これらメッセージをデコードするためのシステムおよび方法を記載している。ＡＩＳ信号はデジタル化されており、異なるドップラーシフトに対するＡＩＳメッセージのコードシーケンスに対応する所定信号とデジタルデータとを相関化することにより、ＡＩＳ候補メッセージを識別している。相関化に使用されるメッセージコードシーケンスは、すべてのＡＩＳメッセージ、すなわち２４ビットのトレーニングシーケンスおよび８ビットのスタートフラグに共通するメッセージのスタートに対応している。ＡＩＳ候補メッセージの検出を改善するために、相関化に使用されるメッセージコードシーケンスは、既知メッセージの一部、例えば船舶識別子（この場合、ＭＭＳＩ番号１０）の一部によって拡張できる。従って、米国特許第７，８７６，８６５号は、直前に検出された送信機からのメッセージを再検出しようと試みるステップに述べているＩＴＵ−Ｒ  Ｍ.２０８４レポートからのアイデアについて記載している。このアプローチは、同じ船舶からの連続する２つのメッセージのうちのかなり多数のビットが同じとなっているという事実を活用するものである。
【０００８】
ＡＩＳの詳細については勧告書ＩＵＴ−Ｒ  Ｍ.１３７１−４（０４／２０１０）およびその付属書に記載されている。
【発明の概要】
【発明が解決しようとする課題】
【０００９】
本発明の目的は、ＡＩＳ送信機の地域マップまたは世界マップを作成するために、衛星によるＡＩＳメッセージの検出を改善することにある。
【課題を解決するための手段】
【００１０】
本発明によれば、ＡＩＳメッセージを検出するための方法は、少なくとも１つのＡＩＳチャンネルにて衛星によって受信された信号を示すデジタルデータを提供するステップと、
少なくとも１つのフラグメント状疑似メッセージがＡＩＳ送信機のデータベース内で識別された送信機に固有である１つ以上のスタティックフィールドからのビットから成る、少なくとも１つのフラグメント状疑似メッセージを含むステップと、デジタルデータと少なくとも１つの前記フラグメント状疑似メッセージとを相関化することにより、デジタルデータ内の候補メッセージを識別するステップとを含む。
前記少なくとも１つのフラグメント状疑似メッセージは、送信時間および／または送信機の位置に依存すると共にＡＩＳ送信機のデータベースに基づいて決定される少なくとも１つのＡＩＳメッセージフラグメントを含む。
【００１１】
従って、本発明に係わる方法は、ＡＩＳ送信機上の情報を含むデータベースに基づくものである。第１の場所において、この情報は衛星のオーバーパス(上空通過)中（すなわち当該デジタルデータに変換された信号の受信時）の衛星のフットプリント内に、どのＡＩＳ送信機を発見すべきかを決定するのに使用される。すべてのＡＩＳメッセージに共通するメッセージフラグメント（特にスタートフラグに続く可能性のあるトレーニングシーケンス）とデジタルデータとを相関化することにより、デジタルデータ内のＡＩＳメッセージをサーチする代わりに、本発明は、信号受信時に衛星によってカバーされているゾーン内において発見すべき（またはデータがリアルタイムで評価されない場合には、発見された）既知の特定の送信機からのＡＩＳメッセージをデータバンクからデジタルデータ内でサーチすることを提案するものである。従って、フラグメント状疑似メッセージは、データベースの助けにより、かつ衛星の軌跡から（例えば、時間を関数とする衛星の位置から）予想可能なメッセージフラグメントから成る。例えば各ＡＩＳメッセージは送信機を不明確でないように定義するフィールド「ユーザーＩＤ」（通常は船舶のＭＭＳＩ番号）を含む。トレーニングシーケンスとスタートシーケンスとデジタルデータと相関化するためのフィールド「ユーザーＩＤ」を含むフラグメント化されたメッセージを使用することにより、特定のＡＩＳ送信機から発信されたＡＩＳメッセージを明瞭にサーチできる。ノイズに対する相関化ピークの高さは、フラグメント化されたメッセージのサイズ（ビット数）が大きくなるにつれ高くなることが理解できよう。本発明によれば、候補メッセージのアクティーブ部分も予想する。ここで、「アクティーブ」部分またはフィールドなる用語は、例えば位置、ナビゲーションなどに関する情報がメッセージごとに変わる可能性が高いＡＩＳメッセージ部分を意味するものと理解する。これまで示したように、フラグメント状疑似メッセージのアクティーブ部分は、送信時間、（例えば、「タイムスタンプ」、「ＵＴＣ年」、「ＵＴＣ月」、「ＵＴＣ日」、「ＵＴＣ時」、「ＵＴＣ分」および／または「ＵＴＣ秒」フィールドからのビット）および／または送信機の位置（例えば、フィールド「経度」および／または「緯度」からのビット）に依存するＡＩＳメッセージの少なくとも１つのフラグメントを含む。種々のＡＩＳメッセージフィールドについては、勧告書ＩＴＵ−ＲＭ.１３７１−４（０４／２０１０）およびその付属書に記載されている。
【００１２】
本発明の方法の好ましい実施形態によれば、フラグメント状疑似メッセージは、送信機の位置に依存する少なくとも１つのＡＩＳメッセージフラグメントを含む。このメッセージフラグメントを合成することは、
〇前記信号の受信時に衛星によるサービスを受けるゾーン内に位置するＡＩＳ送信機を、ＡＩＳ送信機のデータベースに基づいて識別するステップと、
〇前記識別された送信機の、前記データベースに記憶されていた位置に基づき、送信機の位置に応じた少なくとも１つのＡＩＳメッセージフラグメントを決定するステップとを含む。
【００１３】
送信機の位置に応じて決まる少なくとも１つのＡＩＳメッセージフラグメントは、「経度」フィールドからの最大位ビットおよび／または「緯度」フィールドからの最大位ビットを含む。ここで、同一船舶の連続する２回の観察の間で、船舶の位置は少なくともほぼ同じであると仮定する。船舶の新しい位置は、その船舶の最後の既知の位置を中心とする円内に予想でき、この新しい位置は、船舶の最高速度および最終メッセージから経過した時間に対応する。その他のパラメータ、例えば船舶の目的地、その船舶の最終船首方位、その船舶の最終速度、その船舶の最終ナビゲーションステータスなどが既知である場合（例えばデータバンク内に記憶されている場合）、予想は更に正確になり得る。ある送信機が２〜１０秒ごとにその送信機のＡＩＳメッセージを送信すると仮定した場合、衛星の同じオーバーパス中に受信されるＡＩＳメッセージに対する位置の推定値は、極めて正確となる。衛星が２回オーバーパスする間のある時間（例えば１〜２時間）の間、送信機が観測されないことがあり得る。推定される位置がより正確になればなるほど、フィールド「経度」および「緯度」のより多くの最大位ビットが既知となる。
【００１４】
デジタルデータには（例えば衛星により）タイムスタンプを押すことが好ましい。こうして、デジタルデータのタイムスタンプに基づき、送信時間に依存する少なくとも１つのＡＩＳメッセージフラグメントを決定できる。送信時間に依存する少なくとも１つのＡＩＳメッセージフラグメントは、特に「タイムスタンプ」フィールドからのビット、すなわち「ＵＴＣ年」、「ＵＴＣ月」、「ＵＴＣ日」、「ＵＴＣ時」、「ＵＴＣ分」および／または「ＵＴＣ秒」フィールドからのビットを含む。
【００１５】
フラグメント状疑似メッセージは、データベース内で識別された送信機固有の１つ以上のスタティックフィールドからのビットを含むことが好ましい。ここで、「スタティック」フィールドなる用語は所定のＡＩＳ送信機に対してメッセージごとにコンテントが変化しないフィールドを意味するものと理解されよう。スタティックフィールドの例として、「ユーザーＩＤ」フィールド（ＭＭＳＩ数を含む）、「ＡＩＳバージョンインジケータ」フィールド、「ＩＭＯ番号」フィールド（国際船舶協会番号）、「コールサイン」フィールド、「ネーム」フィールド（船名）、「船舶のタイプおよび積荷のタイプ」フィールド、「位置に対する全寸法／基準」フィールド、「電子位置固定デバイスのタイプ」フィールドからのビットがある。
【００１６】
送信時間または送信機の位置に依存するフラグメントとは別に、フラグメント状疑似メッセージとして、データベース内で識別された送信機に固有の１つ以上のアクティーブフィールド、例えば「ＲＡＩＭフラグ」フィールド（送信機自律完全調節）、「特殊操縦インジケータ」フィールド、「推定到着時間」フィールド、「現在の喫水」フィールド、「目的地」フィールド、「ナビゲーションステータス」フィールド、「地上のコース」フィールド、「地上の速度」フィールドおよび／または「通信ステータス」フィールドからのビットがある。
【００１７】
デジタルデータ内の候補メッセージが一旦識別された場合、前記候補メッセージをデコードし、デコードされたメッセージ内に含まれる情報を使ってＡＩＳ送信機の前記データベースを更新することが好ましい。
【００１８】
デジタルデータと複数のフラグメント状疑似メッセージとを相関化することにより、候補メッセージを識別することが好ましい。実際に、特にメッセージ固有の信号が極めて微弱であるか、ノイズが多い場合には、候補メッセージを識別するためにはいくつかの想定を検討しなければならないことがある。その理由の１つは、候補メッセージの予想可能なフラグメントがギャップによって分離される得ることにある。基本的には、ＡＩＳ規格はフィールドの固定された公称長さを定義しているので、ギャップの長さは既知である。しかしながら、スタートフラグまたはエンドフラグと同じビットシーケンスの発生を防止するために、スタッフィングビットを挿入しなければならない場合、フィールド長さが、変わり得る。アプリオリに未知のように、ギャップ（候補メッセージの一部分）がスタッフィングビットを含む場合、候補メッセージ内にスタッフィングビットが存在するか否かを検査するように、複数のフラグメント状疑似メッセージを設計でき、異なるフラグメント状疑似メッセージに対応する相関性をシリアルにおよび／またはパラレルに計算できる。
【００１９】
信号受信時の衛星の位置は種々の方法で得ることができる。例えば、衛星には測位システム（ＧＰＳ、ドリス、ガリレオ、グロナス、コンパスまたはその他）を装備させることができ、こうして衛星はデジタルデータを衛星の軌跡データにリンクできる。別の可能な方法として、衛星コントロールセンター（ＴＴ＆Ｃセンター：テレメートリー、トラッキングおよびコマンド）または例えばレーザー処理測定システムを使用する監視ネットワークにより、衛星の軌跡データを得る方法がある。更に別の可能な方法として、デジタルデータ内の１つ以上のＡＩＳメッセージを識別する送信機に対する三角測量（距離を測定することにより）により、信号受信時の衛星の位置を推定する方法がある。従って、受信したＡＩＳメッセージに基づき、衛星の軌跡を再構成できる。
【００２０】
デジタルデータとフラグメント状疑似メッセージとを相関化するためには、フラグメント状疑似メッセージのスタート点とＡＩＳ時間フレームの時間インターバルのスタート点とを同期化することが好ましい。デジタルデータにはタイムスタンプが押されており、ＡＩＳフレームは、ＵＴＣ時間と同期化されているので、時間的な不確実性の範囲は、ＡＩＳの同期化ジッター（±３ビットに対応する±３１２μｓ）に船舶の位置の不確実性を補償するための所定のバッファを加えた値に低減される。全体として見れば、時間的な不確実性の範囲は、ＵＴＣ時間に対して±６ビットとなると想定できる。
【００２１】
デジタルデータとフラグメント状疑似メッセージとを相関化するために、フラグメント状疑似メッセージに対する計算されたドップラーシフトを送信機の位置の関数として補償することが好ましい。１つのＡＩＳメッセージのドップラーシフトは、衛星とメッセージ送信機との間の相対的速度に応じて決まる。地上基準フレームでは、ＡＩＳ送信機の速度を無視できるので、ドップラーシフトを衛星の基準フレーム内の送信機の位置の関数として表記できる。低軌道衛星の場合、ドップラーシフトは、±４.５ｋＨｚのインターバル内に含まれることになる。送信機の位置は、少なくとも最小１００ｋｍ（最後に検出された位置を中心とする約１００ｋｍの半径）であると判っており、衛星の軌跡も既知であると仮定した場合、ドップラーシフトの不確実性は、５００Ｈｚ未満となる。しかしながら衛星の高さが低い（ターゲット方向と水平線の間の角度が小さい）場合、ドップラーシフトの不確実性は、通常±１００Ｈｚのオーダーとなる。
【００２２】
しかしながらＡＩＳ規格は、±５００Ｈｚの送信機のキャリア周波数に対する誤差を許容しているので、周波数レンジ内の全不確実性は、±１０００Ｈｚ（ドップラーシフトの不確実性＋キャリア周波数に対する誤差）となる。所定のＡＩＳ送信機では、キャリア周波数に対する誤差は、通常時間に対して基本的には一定である。従ってこの誤差は、周波数レンジ内でのサーチを考慮できるようにデータベース内に記憶することが好ましい。従って、データベースからキャリア周波数に対する誤差を抽出できるＡＩＳ送信機の場合、周波数レンジ内の不確実性はドップラーシフト（±１００Ｈｚ）の不確実性だけとなる。データベース内に所定のＡＩＳ送信機のキャリア周波数に対する誤差がセーブされていない場合、相関性ピークを検出した周波数から、および（衛星並びに送信機の位置を使って得た）理論的ドップラーシフトからおよび将来使用するためにこの誤差をデータベース内に記憶することにより、この誤差を計算できる。
【００２３】
最後に、フラグメント状疑似メッセージとデジタルデータとを相関化するために、データベース内に記憶されているＡＩＳ送信機の変調指数も同じように考慮することが好ましい。特定のＡＩＳ送信機のメッセージをサーチするには、その変調指数（ＡＩＳに対する公称値：０.５）を知ることが望ましい。従って、この変調指数がアプリオリに未知である場合、変調指数（例えば０.５±０、１のレンジにある）のための複数の想定をするための相関化を行うことができ、相関性を最大にする変調指数の値はデータベース内に記憶される。
【００２４】
本発明の更なる特徴は、コンピュータによる実行時にコンピュータに上記方法を実行させるコンピュータで実行可能な命令を含むコンピュータプログラムにも関する。本発明の更に別の特徴は、プログラム命令を含むデータ記憶メディア（例えばランダムアクセスメモリまたはリードオンリーメモリ、ＣＤ−ＲＯＭ、ＵＳＢキー、ハードディスクなど）を含むコンピュータプログラム製品にも関する。
説明のために添付図面を参照して、以下に示す好ましい実施形態の詳細な説明から、本発明の上記以外の特性および特徴が明らかとなろう。
【図面の簡単な説明】
【００２５】
【図１】ＡＩＳメッセージのための衛星検出システムの全体のアウトライン図である。
【図２】ＡＩＳメッセージの構造の図である。
【図３】本発明の第１実施形態に係わるＡＩＳメッセージの検出方法の原理を示す略図である。
【図４】ＡＩＳメッセージを検出するためのフラグメント状疑似メッセージの図である。
【図５】船舶の新しい位置を予測するための原理の図である。
【図６】フラグメント状疑似メッセージのうちのメッセージフラグメント（既知のビット）およびギャップ（アプリオリに未知のビット）の図である。
【図７】ＬＥＯ衛星のフェーズドアレイアンテナのフットプリントの略図である。
【図８】本発明の第２実施形態に係わるＡＩＳメッセージの検出方法の原理を示す略図である。
【発明を実施するための形態】
【００２６】
図１は、ＡＩＳメッセージのための衛星検出システムの全体のアウトラインを示す。各々にＡＩＳトランシーバ１２が装備されている複数の船舶１０は、船舶の識別符号（ＭＭＳＩ番号）、ＷＧＳ８４座標における局所限定化（世界ジェオデティックシステム、１９８４年改訂）、船首方位、速度およびその他のステータス情報のようなデータを含むＡＩＳメッセージを送信し、近くにある船舶および沿岸局のＡＩＳ受信機によりこれらメッセージが受信される。こうして各船舶１０は、（約２０〜３０海里の半径内の自己の）ローカル環境におけるＡＩＳ送信機の詳細マップを作成でき、このマップはレーダーおよびその他のナビゲーション援助装置が供給する情報を補足する。
【００２７】
これらＡＩＳメッセージは、低地球軌道内、すなわち１６０〜２０００ｋｍの間の高度に存在し得る（ＬＥＯ）衛星１４によっても受信できる。１つの低軌道衛星のアンテナのフットプリントは、ＡＩＳがメッセージの衝突を防止するセルの大きさを大幅に超えるような広さとなっている地理的ゾーンを一般にカバーする。これらセルは一般に、セル間ではコーディネートされていないので、衛星が受信するＡＩＳ信号は時間および周波数がオーバーラップするメッセージを含み得る。特に船舶密度が高いゾーンで問題が生じる。
【００２８】
（１６２.９７５ＭＨｚを中心とする）ＡＩＳチャンネル１または（１６２.０２５ＭＨｚを中心とする）ＡＩＳチャンネル２（以下、これらチャンネルをＡＩＳ信号と称す）が衛星１４によって受信され、デジタル化される。このように得られたデジタル信号にはタイムスタンプが押され、この信号は、ダウンリンクを通して１つ以上の地上局１６へ送信される。船舶監視システムは、複数のＬＥＯ衛星と複数の地上局とを含むことができることに留意すべきである。これらＬＥＯ衛星は、１つ以上の衛星コントロールセンター１８によって制御されるが、明確にするために、以下、１つのＬＥＯ衛星について述べる。複数の衛星を使用する場合、複数の衛星によるデータ処理をパラレルに実行することが好ましい。
【００２９】
このＬＥＯ衛星１４には、地球測位システム、例えばＧＰＳ受信機を装備することが好ましい。こうすると、衛星はＧＮＳＳ衛星２０から発信された無線ナビゲーション信号２４によりその位置を決定できる。次に、ＡＩＳ信号を示すデジタルデータを、ＬＥＯ衛星の軌跡２２に、全単射的にリンクできる。各データシーケンスが、ＬＥＯ衛星の１つのインスタンスおよび１つの位置に明白に属すようにすることができる。
【００３０】
図２は、ＡＩＳメッセージの構造の図である。各ＡＩＳメッセージは２５６ビットを含み、ランプアップフィールド（８ビット）と、トレーニングシーケンス（２４ビット）と、メッセージスタートフラグ（８ビット）と、データシーケンス（１６８ビット）と、周期的冗長性チェック（１６ビット）と、メッセージエンドフラグ（８ビット）と、バッファ（２４ビット）とにサブ分割されている。勧告書ＩＴＵ−Ｒ  Ｍ.１３７１.４は、スタートバッファだけでなくエンドバッファについても述べている。本明細書では、これら２つのフィールドを１つのフィールドとして扱う。
【００３１】
ＡＩＳメッセージのデータシーケンスも所定の長さのいくつかのフィールドにサブ分割されている。データシーケンスおよび送信されるフィールドの構造は、メッセージのタイプに応じて決まる。各ＡＩＳメッセージは、すべてのＡＩＳメッセージに共通するビットシーケンス、例えばランプアップフィールドと、トレーニングシーケンスと、スタートフラグと、エンドフラグとを含むことに留意すべきである。あるタイプのＡＩＳメッセージしか検討しない場合、フィールド「ＩＤメッセージ」も既知である。その他のビットシーケンスは送信機に固有のものである。これらのうちで、送信機において、不定期間の間でプログラムされる送信機固有の「スタティック」なビットシーケンスとメッセージごとに変わり得る送信機固有の「ダイナミック」なビットシーケンスとを区別することもできる。
次の表は、送信機固有の「スタティック」なビットシーケンスを示す。
【表１】
0005903356.tif 000002
【００３２】
図３は、第１の好ましい実施形態に係わる本発明の方法のための原理の略図を示す。本発明は、メッセージを送信したＡＩＳ送信機が既知であり、かつ更にその送信機に関する補足情報、特にその位置が入手できる場合のＡＩＳメッセージのデータシーケンスの所定のビット数を予測できるということを活用するものである。この予測可能なビットは送信機固有の「スタティック」なフィールドのビットを特に含む。その理由は、ＡＩＳ送信機のデータベース内にこれらフィールドのコンテントをセーブするだけでよいからである。その他の予測可能なビットは、送信機固有の「ダイナミック」フィールドの所定のビットである。例えばフィールド「経度」および「緯度」の最大位ビットを含むが、視界から送信機が消えることがないよう十分高い頻度で当該送信機の位置が発見されることが条件となる。
【００３３】
上記のように、ＬＥＯ衛星のレベルでデータを補足する（ステップＥ１０）。タイムスタンプが押されたデジタルデータを衛星１４へのダイレクトリンクにより地上局１６内に位置し得るデータセンターへ送信するか、または地上局１６と通信状態にある別のロケーションへ送信する。
【００３４】
ＡＩＳメッセージを検出するために、ＬＥＯ衛星から（直接または間接的に）発信されたデジタルデータと通常トレーニングシーケンスとを相関化する。このシーケンスはかなり短いのでノイズの多い信号が存在している場合、相関化ピーク、従って候補メッセージを検出することはできない。デジタルデータ内でまだ検出されていないＡＩＳメッセージのうちのより多くの数のビットを予想するためには、オーバーパスの時間に衛星によってカバーされているゾーン内でＡＩＳメッセージを送信できたＡＩＳ送信機を識別する。
【００３５】
これを行うには、まず最初に衛星の軌跡を決定しなければならない（ステップＥ１２）。例えば測位システムを助けにして衛星のレベルで軌跡（時間を関数とする位置）を計算することが好ましい。この場合、ＡＩＳ信号データと同じ通信リンクを通して衛星の軌跡に関するデータをデータセンターへ送信することができる。これとは異なり、データセンターが任意の通信リンクを通して衛星制御センター１８（ＴＴ＆Ｃセンター）からのＬＥＯ軌跡データを受信することができる。
【００３６】
デジタルデータ内に含まれるＡＩＳメッセージを送信することができたＡＩＳ送信機を識別するために（ステップＥ１４）、ＡＩＳ信号の受信時に衛星の位置を決定する。次に、データベース２６にセーブされている最終位置は、衛星位置のまわりの所定ゾーン内にあるＡＩＳ送信機を示すＡＩＳ送信機２６のデータベースにコンサルトする。この所定のゾーンは、衛星のアンテナ（単数または複数）のフットプリントおよび可能な場合には周辺ゾーンも含むことができる。この周辺ゾーンは、データベースにセーブされた最終位置がフットプリントの外部にあり、かつ船舶位置が検出されてから経過した時間内に進入できる程度にフットプリントに十分近い船舶を考慮する。
【００３７】
デジタルデータ内に含まれるＡＩＳメッセージを送信できた各ＡＩＳ送信機に対し、データベース２６に含まれる情報に基づき予測できるビットを含む、１つまたは複数のフラグメント状疑似メッセージを合成する（ステップＥ１６）。次に、デジタルデータ内の候補メッセージとフラグメント状疑似メッセージとを相関化することにより、その候補メッセージのサーチを実行する（ステップＥ１８）。以下、これらステップについてより詳細に説明する。
【００３８】
デジタルデータ内の候補メッセージが一旦識別されると、その未知部分（すなわち、十分な確度で予想できないビットシーケンス）をデコードする（ステップＥ２０）。この未知の部分は、既知の技術、例えばビタルビアルゴリズムの助けによってもデコードしてもよい。次に、デコードされたメッセージ内に含まれる情報を用いてデータベース２６を更新する（ステップＥ２２）。
【００３９】
次のデジタルデータまでこの方法を続ける。衛星の位置が変化するにつれ、衛星フットプリントに進入したり、このフットプリントから離れる新しいＡＩＳ送信機を識別する。
【００４０】
デジタルデータへの加入者は、このデジタルデータにアクセスすることができ、特にＡＩＳ送信機のマップのディスプレイを得ることができる。（ステップＥ２４）。
【００４１】
ＡＩＳ送信機が衛星のフットプリントから離間するとき、ＡＩＳ送信機は、ＡＩＳメッセージの衛星検出システムの同じ衛星または別の衛星の次のオーバーパスまで、観察されない状態のままとなる。４台のＬＥＯ衛星が配置されたコンステレーションでは、ローカルの無観察時間は、２時間の長さとなり得る。このことは、衛星のフットプリントに進入する送信機に対してセーブされるデータが、平均２時間古いことを意味する。フラグメント状疑似メッセージを合成する際には、この時間を考慮することが好ましい。デジタルデータ内の候補メッセージを識別するには複数の仮想が必要となり得る。これら想定は、パラレルにおよび／またはシリアルに検討してよい。シリアルに検討することを選択するときには、最も可能性の高い確率の想定でスタートしなければならず、検出に失敗した場合には、確率の降順で次の想定を続けなければならない。
【００４２】
例えば送信機が移動しているとデータベースにセーブされている最終ナビゲーション状態が示しているような所定のＡＩＳ送信機に対しては、タイプ１のＡＩＳメッセージを検出しようとすることが好ましい。この場合、デジタルデータと相関化するためのフラグメント状疑似メッセージを正しく構成する。失敗した場合（例えば相関性ピークが検出されない場合）、時間および周波数の同じ想定を有する同じ送信機からのタイプ３のＡＩＳメッセージの検出を試みることができる。データベース内にセーブされている最終ナビゲーション状態が移動していないことを示すＡＩＳ送信機を探すときには、最初にタイプ３のＡＩＳメッセージを検出しようとする。失敗した場合、時間および周波数の同じ想定を有する同じ送信機からのタイプ１のＡＩＳメッセージを検出することを試みることができる。カバーされている地理的ゾーン内に存在すると考えられるすべての送信機から、タイプ５のメッセージをパラレルに検出しようと試みることが好ましい。
【００４３】
次の表は、タイプ１、２または３のＡＩＳメッセージに対する予測可能なビットを示す。一般的に、予測可能なビット数は予想に使用されるモデルの複雑さ、および予想時のデータベースの鮮度に応じて決まる。一般に送信機固有のスタティックフィールドのビットはすべて予想可能である。これと対照的に、メッセージごとに変化し得る送信機固有のダイナミックフィールドに対しては、データベースの鮮度は、これらフィールドのうちのどれだけ多くのビットを予想できるかを決定する。
【表２−１】
0005903356.tif 000003
【表２−２】
0005903356.tif 000004
【００４４】
図４は、ＡＩＳ送信機が通信ステートＳＯＴＤＭＡとなっている固定された（アンカー係留されているかまたはつながれている）船舶からのタイプ１のＡＩＳメッセージを検出するためのフラグメント状疑似メッセージの図である。予想可能なビットシーケンス（メッセージフラグメント）はハッチングされている。これら予想可能なビットシーケンスはいくつかのフィールドにわたって拡張できる。所定のフィールド、例えば経度および緯度のフィールドに対しては、ビットの一部を十分な確度で予想できる。
【００４５】
次の表は、（固定された基地局から送信された）タイプ４のＡＩＳメッセージに対する予想ビットを示す。
【表３】
0005903356.tif 000005
【００４６】
次の表は、（ＡＩＳフレーム内の連続する２つのスロットにわたって拡張された）タイプ５のＡＩＳメッセージに対する予想可能なビットを示す。
【表４】
0005903356.tif 000006
【００４７】
次の表は、タイプ２１のＡＩＳメッセージに対する予想可能なビットを示す（ナビゲーション支援局によって送信され：ＡＩＳフレーム内の連続する２つのスロットにわたって延びている）ことを示す。
【表５】
0005903356.tif 000007
【００４８】
従って、ＡＩＳ送信機が既知の場合、各ＡＩＳメッセージは予想可能なビットのシーケンスで開始する。フラグメント状疑似メッセージの第１フラグを示すこのビットシーケンスは、特に「トレーニングシーケンス」、「スタートフラグ」、「ＩＤメッセージ」、「リピートインジケータ」および「ＩＤユーザー」の各フィールド、すなわち少なくとも７０ビットを含む。識別すべきメッセージのタイプによっては、この第１メッセージフラグメントは大幅に大きくなり得る（タイプ４のメッセージの場合は１１０ビット、タイプ５のメッセージに対しては３０６ビット）。
【００４９】
上記表では、可能なスタッフィングビットは考慮していない。各ＡＩＳメッセージでは、エンドメッセージフラグと同じビットシーケンス（すなわちビットシーケンス［０１１１１１１０］）が発生することを防止するためにデータシーケンスおよびフレームチェックシーケンス内にスタッフィングビットが自動的に挿入される。ＡＩＳ送信機の側では、このことは出力ビットストリーム内に連続する５つの１が発見された場合、これらの１が５つ連続した後にゼロを挿入することを意味する。受信側では１が５つ連続した後の最初のゼロを除去しなければならない。これらスタッフィングビットは、検出すべきメッセージの既知の部分だけでなく、未知の部分も長くし得る。従って、異なる可能なケースを考慮するためにデジタルデータと相関化されたフラグメント状疑似メッセージ内で、ギャップ（未知の部分）の長さを変えなければならない。
【００５０】
フラグメント状疑似メッセージは、ＡＩＳ送信機の送信時間（例えばタイプ１、２、３または２１のメッセージのフィールド「タイムスタンプ」）および／または位置（例えばタイプ１、２、３、４および２１のメッセージのフィールド「経度」および「緯度」）に応じて決まる少なくとも１つのメッセージフラグメントを含む。送信機の位置に応じて決まるメッセージフラグメントのビットを推定（予想）するには、当該ＡＩＳ送信機の（データベースにセーブされている）既知の最終位置を考慮しなければならない。いくつかのシナリオが可能である。特定の送信機からの最後にデコードされたＡＩＳメッセージと同じ送信機からの識別すべきＡＩＳメッセージとの間で、ナビゲーションステートが変化しないと想定することが好ましい。このような想定は、すべてのケースで必ずしも真であるわけではないが、当該観察されていない時間の間（長くて数時間の間）では、この想定は極めてあり得ることである。従って、船舶がアンカー係留されているか、つながれていることをデータベースが表示した場合、船舶の位置は（１０メートル内の）同じ位置に留まるとの想定に基づき、フラグメント状疑似メッセージを合成する。相関ピークを検出できなければ、船舶が航行中であると仮定する。船舶が航行中であることを最後のナビゲーションステートが表示した場合、まず最初にこの想定を検討しなければならない。図５に示されるようなこのようなケースでは、最後の既知の位置３０のまわりの円の内側に、船舶２８の新しいアプリオリに未知の位置を見つけなければならない。ここで、円の半径Ｒは、最後にデコードされたＡＩＳメッセージからの時間ΔＴと船舶の最大速度ｖｍａｘとの積である。すなわちＲ＝ΔＴ×ｖｍａｘである。船舶の最大速度（この値がデータベースに記憶されている場合）をデータベース内で見つけることができる。最大速度が未知である場合、その速度は所定のスレッショルド値（例えば５０ｋｍ／ｈ）よりも低いと想定することができる。円の内部にある陸地３２の面積は無視できる。（しかしながら当該船舶が河川のボートである場合、明らかにナビゲート可能な水路を排除してはならない。）Ｒ＝１００ｋｍの場合の新しい位置の精度では９ビット、フィールド「経度」およびフィールド「緯度」では８ビットを予想できる。データベースで船舶に関する別の情報を入手できる場合、特にタイプ５のメッセージのフィールド「船舶のタイプおよび積み荷のタイプ」、「推定到着時間」および「目的地」が入手できる場合、予想精度を改善できる。例えば、当該船舶が過去のルートの知られているフェリーであると判っている場合、その船舶は恐らくこれらルートのうちの１つにあり、位置の不確実性ゾーンは、かなり低下する。到着推定時間および到着港が判っている場合、このことは、この船舶がとるルートに関する予想も可能にする。
【００５１】
次にデジタルデータ内の候補メッセージに対するサーチステップの好ましい実行（ステップＥ１８）について、より詳細に説明する。これらデジタルデータは、予想ビットのセグメント（メッセージフラグメント）と一義的に相関化されることに留意すべきである。未知のビット（異なるメッセージフラグメント間のギャップ）とＡＩＳ信号のデジタルデータとは相関化されない。デジタルデータとＵＴＣ時間との同期化およびＡＩＳ送信機の位置の最小近似知識により、各メッセージからのＡＩＳメッセージの伝搬遅延時間を計算できる。各ＡＩＳ送信機がそのメッセージをＵＴＣ時間に（直接または間接的に）同期化する際、伝搬遅延時間を考慮することにより、±６ビット（時間的な不確実性レンジ）の精度で所定の送信機からＡＩＳメッセージを開始できるインスタンスを知る。この時間的不確実性のレンジは、船舶の位置の不確実性を補償するためにＡＩＳの同期化ジッター（±３ビットに対応する±３１２μｓ）＋所定のバッファから構成される。相関化のために時間サーチレンジは、時間的不確実性レンジ（例えば１２ビット）に対応する。ＡＩＳフレームの各タイムスロットのスタート時に相関化が実行される。
【００５２】
送信機のキャリア周波数に対するドップラー効果および誤差を検討するために、キャリア周波数のまわりのサーチ周波数レンジにわたって相関化方法を実行する。データベース内にキャリア周波数に対する誤差が含まれているような送信機からのＡＩＳメッセージを探した場合、この誤差を補償することができ、周波数の残留ドメイン内における不確実性は、ドップラーシフトの不確実性（通常±１００Ｈｚ）だけとなる。従って、周波数サーチインターバルは、一般に２００Ｈｚの幅となる。
【００５３】
ハーフビットの各積分の間でのデジタルデータとフラグメント状疑似メッセージとの間の遅延時間を変えることにより、各相関化演算は、２４回の積分を含む。周波数サーチインターバルを（２５Ｈｚ幅の）８つのセルに分割することにより、（以下の変調指数のための異なる想定を含めない）時間／周波数サーチスペース内のセルごとに、総計１９２回の積分が行われる。周波数軸上の不明確性を生じることなく、主たる相関化ピークを検出できるように、周波数におけるセルの幅を選択することが望ましいことに留意すべきである。
【００５４】
キャリア周波数に対するシフトが一旦判ると、理論的なドップラーシフト量を減算でき、残留誤差を使ってデータベース内にセーブされているキャリア周波数に対する誤差を補正できる。
【００５５】
時間／周波数サーチスペースの各セルに対する変調指数の仮想のための相関化方法は、次のステップを含むことが好ましい。
ステップａ．２つの可能な位相（ＮＲＺＩコーディングされたデータ）を使用することにより、デジタルデータと（Ｋ１の既知のビット（図６参照）を有する）第１メッセージフラグメントとを相関化するステップ。ステップａから２つの相関化の結果が得られる。
ステップｂ．デジタルデータと（Ｋ２の既知のビットを有する）第２メッセージフラグメントを相関化するステップ。第１ギャップ（第１フラグメントと第２フラグメントとの間の未知のビット）内の可能なスタッフィングビットを考慮するためには、第１フラグメントから異なる距離に相関化のスタート点を位置させなければならない。第１ギャップの未知のビットの数は、既知である（Ｌ１）。従って、第１ギャップの長さは、Ｌ１＋スタッフィングビットの数となる。後者の数は未知であるので、異なる可能なギャップ長さを試みなければならない（Ｌ１、Ｌ１＋１、....、Ｌ１＋Ｂ１、ここでＢ１はＬ１ビットのシーケンス内のスタッフィングビットの最大数である）。ＮＲＺＩコーディングを考慮するために２つの可能な位相に対して各相関化を実行しなければならない。ステップｂからは２×Ｂ１の相関化結果が得られる。
ステップｃ．デジタルデータと（Ｋ３の既知のビットを有する）第３メッセージフラグメントを相関化するステップ。第２ギャップ（第２フラグメントと第３フラグメントとの間の未知のビット）内の可能なスタッフィングビットを考慮するためには第２フラグメントから異なる距離に相関化のスタート点を位置させなければならない。第２ギャップの未知のビットの数は、既知である（Ｌ２）。従って、第２ギャップの長さはＬ２＋スタッフィングビットの数となる。後者の数は未知であるので、異なる可能なギャップ長さを試みなければならない（Ｌ１、Ｌ１＋１、....、Ｌ１＋Ｂ２、ここでＢ２はＬ２ビットのシーケンス内のスタッフィングビットの最大数である）。ＮＲＺＩコーディングを考慮するように２つの可能な位相に対する各相関化も実行しなければならない。ステップｂからは２×Ｂ２の相関化結果が得られる。
ステップｄ．等々
ステップｅ．パラレルな処理技術を採用する場合、ステップａ、ｂ、ｃおよびｄの可能な組み合わせに従い、異なる相関化の結果を加算しなければならない。従って、テストされる想定の数は、２×（２×Ｂ１）（２×Ｂ２）....となる。
ステップｆ．検討された想定のうちで、（フラグメントの組に対し）最大の相関化ピークを有する結果のサーチを行う。更に確実な検出をするには、相関化ピークが、経験的に強化された（ノイズに対する）分析的にあらかじめ定められた所定のスレッショルドを超えていなければならない。
【００５６】
すべての想定（すべての可能な組み合わせ）を検討する代わりに、ステップａ、ｂ、ｃおよびｄに対し、各メッセージフラグメントに対して最大の相関化ピークを発生した想定だけを維持してもよい。このオプションには、フラグメントの正しい組み合わせに到達することができないという、より高いリスクが存在する。この欠点を緩和するために、ステップａ、ｂ、ｃおよびｄを通過する際に、最大の相関化ピークを有する限られた数の組み合わせを保持できる可能性があり、その他の組み合わせは廃棄する。
【００５７】
ＬＥＯ衛星は、各々がＡＩＳ信号を受信する一群の基本アンテナ（フェーズドアレイアンテナを構成する）を含むことが好ましい。次に、これらアンテナ群の意図する指向性が得られるように、異なるアンテナが受信した信号を変換することによって得られるデジタルデータを組み合わせることができる。この場合、アンテナ群によってカバーされる海抜におけるゾーンを人工的に縮小し、一度に検討すべきＡＩＳ送信機の数（および衝突メッセージの数）を減少できる。換言すれば、このゾーンの外の地理的ゾーンを縮小しながら、所定の地理ゾーンに対するアンテナの感度を高めることができる。次に、異なるデジタルデータフローの間の位相関係を変えることによって、アンテナの群のフットプリントをスイープできる。この周波数選択技術（狭周波数セル）によって、干渉する信号のパワーを大幅に低下させることが可能となっている。
【００５８】
図７は、衛星のアンテナの群のフットプリント３４を略図で示している。衛星は（矢印３６が示す）速度ｖｓａｔで移動し、これによって衛星が受信するＡＩＳ信号のドップラーシフトが生じる。参照番号３８により等ドップラーラインが示されている。ハッチングのつけられたバンド４０は、複数の周波数の１つのセルに対応している。フェーズドアレイアンテナにより、検出ゾーンをフットプリント３４よりも実質的に狭いローカルゾーン４２に縮小することが可能となっている。次に、データベース内に記録された最終位置の結果が、このローカルゾーン４２に存在し得る船舶からのＡＩＳメッセージを検出する試みを行うことができる。ローカルゾーン４２の内部では、ドップラーシフトの範囲はフットプリント３４に対してかなり狭く、従って各ローカルゾーン４２に対してテストすべき周波数セルの数は大幅に減少する。
【００５９】
図８は、第２の好ましい実施形態に係わる本発明の方法に対する原理の略図を示す。この図の方法のステップと同一または類似のステップを示すのに、同じ参照番号が使用されている。前の例と同じように、ＬＥＯ衛星の高さでデータを収集する（ステップＥ１０）。コンピュータセンターにはタイムスタンプが押されたデジタルデータを送信する。図３の方法と異なる主な違いは、衛星の位置を送信機の位置までの距離を三角測量で得ている点にあり、送信機のＡＩＳメッセージはデジタルデータで識別される。ＬＥＯ衛星に地球測位システムが装備されていない場合、このシステムが故障している場合、または制御センターにより衛星に関する軌道データが伝送されていない場合、本発明のこの変形例は重要である。
【００６０】
この方法を開始するには、デジタルデータ内に含まれたＡＩＳメッセージを送信できたＡＩＳ送信機を識別する（ステップＥ１４）。初期の可能性を低減できるようにするには、衛星位置の初期推定値を入手することが好ましい（ステップＥ１２’）。衛星位置の初期推定値が入手できない場合、（相関化のためのＡＩＳメッセージすべてに共通するスタティック部分だけを使用することにより）従来の方法でデジタルデータ内の最初のＡＩＳメッセージを検出することを試みることができる。最初のＡＩＳメッセージが一旦デコードされると、受信された最初のＡＩＳ送信機の位置が判り、この位置を衛星位置の大ざっぱな推定値として使用できる。（上記ステップＥ１６、Ｅ１８およびＥ２０を続けることにより）発見された第１送信機の近くにある船舶からのデジタルデータ内のメッセージを探す試みを行う。（異なる送信機からの）少なくとも４つのＡＩＳメッセージが一旦検出されると、次の式の系を解くことにより、衛星位置を推定できる（ステップＥ１２”）
【数１】
0005903356.tif 000008
ここで、Ｔｉ，ｉ＝Ａ，．．．，Ｄは、衛星の時間系（衛星のクロック）内の送信機ｉか
0005903356.tif 000009
メッセージの送信時間であり、δＳＡＴ-ＵＴＣ（ｔ）は、衛星の時間基準と基準ＵＴＣの項との間のバイアスである。ＡＩＳフレームは、ＵＴＣ時間に同期化されているので、
0005903356.tif 000010
されているので、位置ｉ＝Ａ，．．．．，Ｄは、既知である。衛星のクロックが４つの測定値に対して安定しており、従ってＵＴＣ時間に対する安定したバイアスを示すこと、および一般に１００ｍ未満の誤差を有する時間Ｔｉ、ｉ＝Ａ，．．．．，Ｄに対して衛星位置が同じであると仮定できる。Ｔｍ＝１／４（ＴＡ＋ＴＢ＋ＴＣ＋ＴＤ）と設定すると、次の式が得られる。
【数２】
0005903356.tif 000011
【００６１】
0005903356.tif 000012
考慮することにより、より高い精度に位置を決定できる。最小二乗方法を使用するか、または衛星軌道の数学的モデルに基づくカルマンフィルタを使用することが望ましい。
【００６２】
このように、衛星位置を見つけた後に、衛星が受信したメッセージ、従ってデジタルデータに含まれるメッセージを送信できたＡＩＳ送信機を、より高い精度で識別できる。既に処理したデータ内のメッセージをサーチし、次にデジタルデータ内のメッセージをサーチし続けることができる。
【００６３】
図３の例に示されるようにデータベースを更新する。
【符号の説明】
【００６４】
１０    船舶
１２    ＡＩＳトランシーバ
１４    ＬＥＯ衛星
１６    地上局
１８    コントロールセンター
２０    ＧＮＳＳ衛星
２２    ＬＥＯ衛星の軌跡
２４    無線ナビゲーション信号
２６    データベース
２８    船舶の新しい位置
３０    船舶の最終位置
３２    閉じた陸上エリア
３４    フットプリント
３６    ｖｓａｔを示す矢印
３８    等ドップラーライン
４０    周波数のセルに対応するハッチングのつけられたバンド
４２    ローカルゾーン
Ｅ１０    データを補足するステップ
Ｅ１２    衛星の軌跡を決定するステップ
Ｅ１２’  衛星位置を最初に推定するステップ
Ｅ１２”  衛星位置を推定するステップ
Ｅ１４    ＡＩＳ送信機を識別するステップ
Ｅ１６    フラグメント状疑似メッセージを合成するステップ
Ｅ１８    デジタルデータ内の候補メッセージをサーチするステップ
Ｅ２０    発見されたメッセージの未知の部分をデコードするステップ
Ｅ２２    データベースを更新するステップ
Ｅ２４    ＡＩＳ送信機のマップをディスプレイするステップ
Ｒ        半径
【図１】
0005903356.tif 000013
【図２】
0005903356.tif 000014
【図３】
0005903356.tif 000015
【図４】
0005903356.tif 000016
【図５】
0005903356.tif 000017
【図６】
0005903356.tif 000018
【図７】
0005903356.tif 000019
【図８】
0005903356.tif 000020