[go: up one dir, main page]

JP7038281B2 - 測位方法および測位端末 - Google Patents

測位方法および測位端末 Download PDF

Info

Publication number
JP7038281B2
JP7038281B2 JP2019523397A JP2019523397A JP7038281B2 JP 7038281 B2 JP7038281 B2 JP 7038281B2 JP 2019523397 A JP2019523397 A JP 2019523397A JP 2019523397 A JP2019523397 A JP 2019523397A JP 7038281 B2 JP7038281 B2 JP 7038281B2
Authority
JP
Japan
Prior art keywords
solution
positioning
coordinates
fixed
fixed solution
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.)
Active
Application number
JP2019523397A
Other languages
English (en)
Other versions
JPWO2018225421A1 (ja
Inventor
靖久 山崎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Panasonic Intellectual Property Management Co Ltd
Original Assignee
Panasonic Intellectual Property Management Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Panasonic Intellectual Property Management Co Ltd filed Critical Panasonic Intellectual Property Management Co Ltd
Publication of JPWO2018225421A1 publication Critical patent/JPWO2018225421A1/ja
Application granted granted Critical
Publication of JP7038281B2 publication Critical patent/JP7038281B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/396Determining accuracy or reliability of position or pseudorange measurements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/35Constructional details or hardware or software details of the signal processing chain
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/421Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • G01S19/426Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/45Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/48Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system
    • G01S19/49Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system whereby the further system is an inertial position system, e.g. loosely-coupled
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/07Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/14Receivers specially adapted for specific applications
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/43Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Description

本開示は、測位衛星(以下、測位に利用できる人工衛星を総称して「衛星」とする)からの信号を利用して干渉測位を行う場合の測位方法および測位端末に関する。
従来、静止状態の対象物を高精度に測量するために、RTK(Real Time Kinematic)法による干渉測位(RTK演算)が利用されている。RTK法とは、衛星が送信する測位信号の搬送波位相積算値を用いて所定の地点の測位を行うものである。このRTK法による干渉測位を、移動体の測位に適用することにより、移動体の高精度な測位を実現することが期待されている。
車輌等の移動体に取り付けられた測位端末は、RTK演算を行う際、GNSS(GlobalNavigation Satellite System)の衛星(図示せず)からの測位信号を受信する。なお、GNSSとは、GPS(Global Positioning System)、BeiDou、GLONASS等の民間航空航法に使用可能な性能(精度・信頼性)を持つ衛星航法システムの総称である。測位信号には、GPS衛星から送信されるL1信号(1575.42MHz)、L2信号(1227.60MHz)等がある。
特許文献1には、衛星からの測位信号に基づく自車輌の位置(GPS位置)と、自車輌の速度に基づく自車輌の位置(推測航法位置)の両方を測定し、GPS位置の測位精度が、設定レベルを満足する場合にはGPS位置を出力し、設定レベルを満足しない場合には推測航法位置を出力する測位端末が開示されている。
このように、従来から、RTK演算等の測位演算において、高精度な測位解(フィックス解)が算出された場合には該測位解(フィックス解)を出力し、低精度な測位解(フロート解)が算出された場合には推測航法により推測された座標を出力する技術が開示されている。
特開平7-104847号公報
都市部など、遮蔽が多い環境では、RTK演算において、実際の移動体の座標と大きく異なるフィックス解(以下、「誤フィックス」という)が算出されることが散発する。
本開示の一態様は、誤フィックスを出力することを低減し、測位の精度を向上させることができる測位方法および測位端末を開示する。
本開示の一態様に係る測位方法は、測位端末が複数の衛星から送信される情報に基づいて測位演算を行うことで、移動体の座標を決定する測位方法であって、前記測位端末は、前記測位演算によって得られる解であるフロート解または前記フロート解より精度の高いフィックス解を算出し、前記フィックス解の示す座標を基準として、前記移動体の座標の推測値であるDR(Dead Reckoning)解を算出し、前記フィックス解が第1時間連続して得られている状態において、前記測位演算によって新たにフィックス解が得られた場合、前記新たに得られたフィックス解と前回のフィックス解との差が閾値以上の場合には前記DR解を前記移動体の座標として出力する。
本開示の一態様に係る測位端末は、複数の衛星から送信される測位信号を受信する受信部と、前記測位信号に含まれる情報に基づいて測位演算を行うことで、移動体の座標を決定するプロセッサと、を具備し、前記プロセッサは、前記測位演算によって得られる解であるフロート解または前記フロート解より精度の高いフィックス解を算出し、前記フィックス解の示す座標を基準として、前記移動体の座標の推測値であるDR(Dead Reckoning)解を算出し、前記フィックス解が第1時間連続して得られている状態において、前記測位演算によって新たにフィックス解が得られた場合、前記新たに得られたフィックス解と前回のフィックス解との差が閾値以上の場合には前記DR解を前記移動体の座標として出力する。
本開示の一態様によれば、誤フィックスを出力することを低減し、測位の精度を向上させることができる。
図1は、一実施の形態に係る測位システムの構成を示す図である。 図2は、一実施の形態に係る基準局の構成を示すブロック図である。 図3は、一実施の形態に係る測位端末の構成を示すブロック図である。 図4は、一実施の形態に係る測位端末の座標出力の一例を示す図である。 図5Aは、一実施の形態に係る測位処理を示すフロー図である。 図5Bは、一実施の形態に係る測位処理を示すフロー図である。 図5Cは、一実施の形態に係る測位処理を示すフロー図である。
以下、図面を適宜参照して、本開示の実施の形態について、詳細に説明する。但し、必要以上に詳細な説明は省略する場合がある。例えば、既によく知られた事項の詳細説明や実質的に同一の構成に対する重複説明を省略する場合がある。これは、以下の説明が不必要に冗長になるのを避け、当業者の理解を容易にするためである。
なお、添付図面および以下の説明は、当業者が本開示を十分に理解するために、提供されるのであって、これらにより請求の範囲に記載の主題を限定することは意図されていない。
<測位システムの構成>
まず、本実施の形態に係る測位システム1の構成について図1を用いて説明する。図1に示すように、測位システム1は、基準局10と、測位端末20と、から構成される。基準局10は、地球上の座標が既知である箇所に設置される。測位端末20は、座標を求める対象である移動体(例えば車輌など)に設置される。
測位システム1は、測位端末20の位置を計測し、測位端末20の地球上の座標を求める。座標は、例えば、緯度・経度・高度の三次元座標が一般的であるが、緯度・経度などの二次元座標であってもよい。
基準局10は、GNSS衛星から受信した測位信号に基づいて基準局10の測位データ(以下、「基準局測位データ」という)を生成し、測位端末20に送信する。なお、測位データの詳細については後述する。
測位端末20は、GNSS衛星から受信した測位信号に基づいて測位端末20の測位データ(以下、「測位端末測位データ」という)を生成し、基準局測位データ及び測位端末測位データを用いて1エポック毎にRTK法による干渉測位処理を行い、移動体の座標を出力する。エポックは、データ取得時刻のことであり、エポック間隔はデータ取得時刻の時間間隔(周期)を表す時間単位である。例えば、測位端末20が5Hzで動作する場合、1秒間に5つのデータが取得されるため、エポック間隔は0.2秒となる。なお、測位端末20には、測位用の専用端末、測位機能を有するパーソナルコンピュータ、スマートフォン、タブレット、測位サービスを行うサーバ等が含まれる。
<基準局の構成>
次に、本実施の形態に係る基準局10の構成について図2を用いて説明する。図2に示すように、基準局10は、プロセッサ101と、記憶部102と、入力部103と、出力部104と、通信部105と、受信部106と、バス110と、を有している。
プロセッサ101は、バス110を介して基準局10の他の要素を制御する。プロセッサ101として、例えば、汎用CPU(Central Processing Unit)が用いられる。また、プロセッサ101は、所定のプログラムを実行することにより、測位信号に基づいて基準局測位データを生成する。
記憶部102は、他の要素から様々な情報を取得し、一時的あるいは恒久的にその情報を保持する。記憶部102は、いわゆる一次記憶装置と二次記憶装置の総称である。記憶部102は、物理的に複数配置されても良い。記憶部102として、例えば、DRAM(Direct Random Access Memory)、HDD(Hard Disk Drive)、SSD(Solid State Drive)が用いられる。
入力部103は、外部からの情報を受け付ける。入力部103が受け付ける外部からの情報には、基準局10の操作者からの入力に関する情報などが含まれる。一例としてキーボード等の入力インターフェースを用いることで入力部103を構成することができる。
出力部104は、外部へ情報を提示する。出力部104が提示する情報には、測位に関する情報などが含まれる。一例としてディスプレイ等の既存の出力インターフェースを用いることで出力部104を構成することができる。
通信部105は、通信路を介して外部の機器と通信を行う。通信部105が通信する対象(通信対象)の機器には、測位端末20が含まれる。一例として無線LAN通信網、3G通信網など既存の通信網と通信可能な通信インターフェースを用いることで通信部105を構成することができる。
受信部106は、衛星からの測位信号を受信し、バス110を介して測位信号をプロセッサ101に出力する。
なお、上記の基準局10の構成は一例である。基準局10の各構成要素の一部を統合して構成することもできる。基準局10の各構成要素の一部を複数の要素に分割して構成することもできる。基準局10の各構成要素の一部を省略することもできる。基準局10に他の要素を付加して構成することもできる。
<測位端末の構成>
次に、本実施の形態に係る測位端末20の構成について図3を用いて説明する。図3に示すように、測位端末20は、プロセッサ201と、記憶部202と、入力部203と、出力部204と、通信部205と、受信部206と、バス210と、を備えている。
プロセッサ201は、バス210を介して測位端末20の他の要素を制御する。プロセッサ201として、例えば、汎用CPUが用いられる。また、プロセッサ201は、所定のプログラムを実行することにより、測位信号に基づいて測位端末測位データを生成する。また、本実施の形態では、プロセッサ201が、移動体の座標を出力する機能を備えている。このプロセッサ201の機能の詳細については後述する。
記憶部202は、他の要素から様々な情報を取得し、一時的あるいは恒久的にその情報を保持する。記憶部202は、いわゆる一次記憶装置と二次記憶装置の総称である。記憶部202は、物理的に複数配置されても良い。記憶部202として、例えば、DRAM、HDD、SSDが用いられる。
入力部203は、外部からの情報を受け付ける。入力部203が受け付ける外部からの情報には、測位端末20の操作者からの入力に関する情報などが含まれる。一例としてキーボード等の入力インターフェースを用いることで入力部203を構成することができる。
出力部204は、外部へ情報を提示する。出力部204が提示する情報には、測位に関する情報などが含まれる。一例としてディスプレイ等の既存の出力インターフェースを用いることで出力部204を構成することができる。
通信部205は、通信路を介して外部の機器と通信を行う。通信部205が通信する対象(通信対象)の機器には、基準局10が含まれる。一例として無線LAN通信網、3G通信網など既存の通信網と通信可能な通信インターフェースを用いることで通信部205を構成することができる。
受信部206は、衛星からの測位信号を受信し、バス210を介して測位信号をプロセッサ201に出力する。
なお、上記の測位端末20の構成は一例である。測位端末20の各構成要素の一部を統合して構成することもできる。測位端末20の各構成要素の一部を複数の要素に分割して構成することもできる。測位端末20の各構成要素の一部を省略することもできる。測位端末20に他の要素を付加して構成することもできる。
<測位端末のプロセッサの移動体座標出力機能>
次に、測位端末20のプロセッサ201の移動体の座標を出力する機能について詳細に説明する。
プロセッサ201は、基準局測位データと測位端末測位データとに基づいて1エポック毎にRTK法による干渉測位(RTK演算)を実行し、測位解(フィックス解またはフロート解)を算出する。以下、RTK演算によって得られる測位解を「RTK測位解」という。プロセッサ201は、RTK演算によって得られるAR(Ambiguity Ratio)値を用いて品質チェックを行い、AR値が所定の閾値(例えば3.0)以上の場合には、RTK測位解がフィックス解であると判定し、AR値が所定の閾値(例えば3.0)未満の場合には、RTK測位解がフロート解であると判定する。
また、プロセッサ201は、移動体の以前の座標と移動体の速度に関する情報とに基づいて、移動体の現在の座標を推測する。移動体の以前の座標とは、例えば、前回(すなわち、1エポック前に)出力された座標である。プロセッサ201は、1エポック毎に移動体の座標を出力する。なお本開示において移動体の以前の座標と対比する意味で「移動体の現在の座標」という語を用いるが、「移動体の現在の座標とは」「移動体の以前の座標」より後に移動体が存在すると推測される座標を意味し、「現在」は必ずしも座標出力装置が動作している時刻を意味するものではない。移動体の速度に関する情報は、例えば、測位データに含まれるドップラー周波数から算出した移動体の速度、または、移動体に備えられた速度検出モジュール(図示せず)から出力される移動体の速度である。速度の情報は、例えば「南方向に秒速Xメートル、東方向に秒速Yメートル、高さ方向に秒速Zメートル」など所定の次元に対する単位時間当たりの移動量(移動方向ベクトル)の対で定義される情報である。プロセッサ201は、ドップラー周波数情報またはセンサ情報に基づいて移動方向ベクトル(1エポック分の移動量)を算出し、移動体の以前の座標(例えば、1エポック前の座標)に、移動方向ベクトルを加算して、移動体の現在の座標を推測する。このようにして推測された座標は、推測航法座標あるいはDR解(Dead Reckoning解)とも呼ばれる。
また、プロセッサ201は、フィックス解が得られた場合、フィックス解が連続して得られている時間(以下、「第1時間」という)を計時し、第1時間が所定時間T1(例えば20秒)に達したか否かを判定する。また、プロセッサ201は、DR解が連続して出力されている時間(以下、「第2時間」という)を計時し、第2時間が所定時間T2(例えば60秒)に達したか否かを判定する。なお、プロセッサ201は、第1時間、第2時間に相当する測位処理(データ取得)の回数をカウントしても良い。例えば、エポック間隔が0.2秒である場合、フィックス解が100回連続して得られると第1時間が所定時間T1(=20秒)に達したことになり、DR解が300回連続して出力されると第2時間が所定時間T2(=60秒)に達したことになる。
そして、プロセッサ201は、所定の条件により、RTK測位解(フィックス解またはフロート解)あるいはDR解を移動体の現在の座標とし、1エポック毎に移動体の現在の座標を出力部204に出力する。なお、プロセッサ201が測位を行い、移動体の現在の座標としてRTK測位解あるいはDR解を選択して出力するまでの測位処理のフローについては後述する。
<測位データ>
次に、測位データについて説明する。本実施の形態において測位データには擬似距離情報、搬送波位相情報およびドップラー周波数情報が含まれる。
擬似距離情報とは、衛星と自局(基準局10あるいは測位端末20)との距離に関する情報である。プロセッサ(プロセッサ101あるいはプロセッサ201)は、測位信号を解析することにより衛星と自局との距離を算出することができる。具体的には、プロセッサは、まず、(1)測位信号が搬送したコードのパターンと自局が生成したコードのパターンとの相違、および、(2)測位信号に含まれるメッセージ(NAVDATA)に含まれる衛星の信号生成時刻と自局の信号受信時刻、の2つの情報に基づいて測位信号の到達時間を求める。そして、プロセッサは、当該到達時間に光速を乗ずることにより衛星と自局との擬似距離を求める。この距離には衛星のクロックと自局のクロックとの相違等に起因する誤差が含まれる。通常、この誤差を軽減させるために4つ以上の衛星に対して擬似距離情報が生成される。
搬送波位相情報とは、自局が受信した測位信号の位相である。測位信号は所定の正弦波である。プロセッサは、受信した測位信号を解析することにより測位信号の位相を算出することができる。
ドップラー周波数情報とは、衛星と自局との相対的な速度に関する情報である。プロセッサは、測位信号を解析することによりドップラー周波数情報を生成することができる。
以上のようにして、基準局10のプロセッサ101および測位端末20のプロセッサ201によって、それぞれ測位データが生成される。
<RTK演算>
RTK演算について説明する。RTK演算は干渉測位の一つであるRTK法を実行する演算である。
RTK法とは、衛星が送信する測位信号の搬送波位相積算値を用いて所定の地点の測位を行うものである。搬送波位相積算値とは、衛星から所定の地点までの(1)測位信号の波の数と(2)位相との和である。搬送波位相積算値が求まれば、測位信号の周波数(および波長)が既知であるので、衛星から所定の地点までの距離を求めることができる。測位信号の波の数は、未知数であるので整数アンビギュイティまたは整数値バイアスと呼ばれる。
RTK法を実行するにあたって重要なことはノイズの除去と、整数アンビギュイティの推定(決定)である。
RTK法では、二重差と呼ばれる差を演算することにより、ノイズの除去を行うことができる。二重差とは2つの衛星に対する1つの受信機の搬送波位相積算値の差(一重差)を2つの受信機(本実施の形態においては基準局10と測位端末20)の間でそれぞれ算出した値の差である。本実施の形態においてはRTK法を用いた測位のために4つ以上の衛星を使用する。従って、4つ以上の衛星の組み合わせの数だけ二重差を演算することになる。この演算では、基準局測位データおよび測位端末測位データが用いられる。
RTK法では、整数アンビギュイティの推定を様々な方法で行うことができる。例えば、(1)最小二乗法によるフロート解の推定、および、(2)フロート解に基づくフィックス解の検定という手順を実行することにより整数アンビギュイティの推定を行うことができる。
最小二乗法によるフロート解の推定は、時間単位毎に生成した二重差の組み合わせを用いて連立方程式を作成し、作成した連立方程式を最小二乗法によって解くことにより実行される。連立方程式はエポック毎に生成される。この演算では、基準局測位データ、測位端末測位データおよび基準局10の既知の座標が用いられる。このようにして求められた整数アンビギュイティの実数推定値をフロート解(推測解)と呼ぶ。
以上のようにして求められたフロート解は実数であるのに対して、整数アンビギュイティの真の値は整数である。よって、フロート解を丸めることにより整数値にする作業が必要になる。しかし、フロート解を丸める組み合わせには複数通りの候補が考えられる。従って、候補の中から正しい整数値を検定する必要がある。検定によって整数値バイアスとしてある程度確からしいとされた解をフィックス解(精密測位解)と呼ぶ。本実施の形態ではRTK演算によって得られるAR値を用いて品質チェックを行い、品質チェックの結果に基づいて正しい整数値を検定する。なお、整数値の候補の絞込みを効率化するために基準局測位データが用いられる。
<座標出力の一例>
次に、本実施の形態に係る測位端末20の座標出力の一例について、図4を用いて説明する。なお、図4では、フィックス解は四角印、フロート解は三角印、DR解は丸印で図示される。また、出力される解は黒塗り、出力されない解は白塗りで図示される。
測位端末20は、フィックス解が第1時間(図4の例では5個)連続して得られるまではRTK測位解(フィックス解あるいはフロート解)を出力する(領域401)。
また、測位端末20は、フィックス解が第1時間以上連続して得られている状態で、RTK演算においてフロート解が算出された場合、DR解を算出して出力する(タイミング402)。その後、フィックス解が第1時間連続して得られるまではDR解を出力する(領域403)。
また、測位端末20は、フィックス解が第1時間以上連続して得られている状態で、新たに得られたフィックス解に基づく現在の速度ベクトルと前回の速度ベクトルとの差分である速度差ベクトルの大きさが閾値より小さい場合にはフィックス解を出力する(領域404、406)。
一方、測位端末20は、フィックス解が第1時間以上連続して得られている状態で、速度差ベクトルの大きさが閾値以上の場合にはDR解を出力する(領域405)。
また、測位端末20は、DR解が第2時間(図4の例では10個)連続して出力された場合には(領域407)、その後、RTK測位解(フィックス解あるいはフロート解)を出力する(領域408)。
<測位処理のフロー>
次に、本実施の形態に係る測位処理のフローについて図5A、図5B、図5Cを用いて説明する。なお、本実施の形態では、測位端末20が測位処理を行う例を説明する。ただし、本開示の測位処理は、測位端末20によって行われるものに限定されず、例えば、測位システム1の内部に追加された汎用コンピュータによって実行されても良い。測位処理を開始するタイミングについては特に限定は無い。例えば、測位端末20の電源が投入された際に、測位処理を開始しても良い。また、測位端末20の入力部203によって測位処理を開始するコマンドが入力された際に、測位処理を開始しても良い。
まず、ST501において、プロセッサ201は、動作モードを「通常モード」に設定する。なお、「通常モード」とは、DR解を出力せず、常に、RTK測位解(フィックス解またはフロート解)を出力するモードである。
次に、ST502において、受信部206が、受信可能な全ての衛星のそれぞれから測位信号を受信する。そして、ST503において、プロセッサ201が、測位信号を用いて測位端末測位データを生成する。また、ST504において、通信部205が基準局10から基準局測位データを受信する。
次に、ST505において、プロセッサ201が、基準局測位データ、測位端末測位データを用いてRTK演算を実行し、現在のRTK測位解Xを算出する(nは現在を基準とする時刻の変数である)。
次に、ST506において、プロセッサ201が、現在のRTK測位解Xを記憶部202に記憶させる。
ST506の後、「通常モード」では(ST507:NO)、ST508において、プロセッサ201が、フィックス解が連続して得られている第1時間が所定時間T1に達しているか否か、すなわち、現在までの所定時間T1におけるRTK測位解Xn-T1~Xの全てがフィックス解か否かを確認する。
第1時間が所定時間T1に達している場合、すなわち、RTK測位解Xn-T1~Xの全てがフィックス解である場合(ST509:YES)、ST510において、プロセッサ201が、現在のRTK測位解Xと前回のRTK測位解Xn-1との差分により現在の速度ベクトルVを算出する。
さらに、ST511において、プロセッサ201が、前回のRTK測位解Xn-1と前々回のRTK測位解Xn-2との差分により前回の速度ベクトルVn-1を算出する。
そして、ST512において、プロセッサ201が、現在の速度ベクトルVと前回の速度ベクトルVn-1との差分である速度差ベクトルAを算出する。
次に、ST513において、プロセッサ201が、速度差ベクトルの大きさ|A|と所定の閾値sとを比較する。
速度差ベクトルの大きさ|A|が閾値sより小さい場合(ST514:YES)、ST515において、プロセッサ201が、動作モードを「DRモード」に設定する。なお、「DRモード」とは、DR解またはRTK測位解(フィックス解)を選択的に出力するモードである。そして、ST516において、プロセッサ201が、移動体の現在の座標として、現在のRTK測位解X(フィックス解)を出力する。
一方、速度差ベクトルの大きさ|A|が閾値s以上の場合(ST514:NO)、ST516において、プロセッサ201が、移動体の現在の座標として、現在のRTK測位解X(フィックス解)を出力する。
また、RTK測位解Xn-T1~Xのいずれかがフィックス解ではない(フロート解である)場合も(ST509:NO)、ST516において、プロセッサ201が、移動体の現在の座標として、現在のRTK測位解X(フィックス解あるいはフロート解)を出力する。
ST516の後、測位を継続しない場合(ST517:NO)、測位処理は終了となる。一方、測位を継続する場合(ST517:YES)、ST518において、プロセッサ201が、DR解が連続して出力されている第2時間が所定時間T2に達しているか否かを確認する。
第2時間が所定時間T2に達していない場合(ST519:NO)、フローはST502に戻る。また、第2時間が所定時間T2に達している場合(ST519:YES)、フローはST501に戻る。
ST506の後、「DRモード」では(ST507:YES)、ST520において、プロセッサ201が、フィックス解が連続して得られている第1時間が所定時間T1に達しているか否か、すなわち、現在までの所定時間T1におけるRTK測位解Xn-T1~Xの全てがフィックス解か否かを確認する。
第1時間が所定時間T1に達している場合、すなわち、RTK測位解Xn-T1~Xの全てがフィックス解である場合(ST521:YES)、ST522において、プロセッサ201が、現在のRTK測位解Xと前回のRTK測位解Xn-1との差分により現在の速度ベクトルVを算出する。
さらに、ST523において、プロセッサ201が、前回のRTK測位解Xn-1と前々回のRTK測位解Xn-2との差分により前回の速度ベクトルVn-1を算出する。
そして、ST524において、プロセッサ201が、現在の速度ベクトルVと前回の速度ベクトルVn-1との差分である速度差ベクトルAを算出する。
次に、ST525において、プロセッサ201が、速度差ベクトルの大きさ|A|と所定の閾値sとを比較する。
速度差ベクトルの大きさ|A|が閾値sより小さい場合(ST526:YES)、ST527において、プロセッサ201が、移動体の現在の座標として、現在のRTK測位解X(フィックス解)を出力する。その後、フローはST517に進む。
一方、速度差ベクトルの大きさ|A|が閾値s以上の場合(ST526:NO)、ST528において、プロセッサ201が、前回の出力値に移動方向ベクトルを加算してDR解を算出する。
また、RTK測位解Xn-T1~Xのいずれかがフィックス解ではない(フロート解である)場合も(ST521:NO)、ST528において、プロセッサ201が、前回の出力値に移動方向ベクトルを加算してDR解を算出する。
そして、ST528の後、ST529において、プロセッサ201が、移動体の現在の座標として、DR解を出力する。その後、フローはST517に進む。
<効果>
このように、本実施の形態では、所定時間連続してフィックス解が得られている状態において、フィックス解が得られた場合であっても、該フィックス解の速度差ベクトルの大きさ|A|が閾値s以上の場合には、該フィックス解が誤フィックスであると推測し、DR解を出力する。
これにより、誤フィックスを出力することを低減することができる。また、この場合に出力されるDR解は、所定時間連続してフィックス解が得られた後の最新のフィックス解を基準にして求められたものであり、この最新のフィックス解は、正しい座標を示している可能性が高い。そのため、このような状況では、反射などにより誤った座標にフィックスしている可能性が高いフィックス解(すなわち、速度差ベクトルが大きいフィックス解)よりも、過去の正しい座標を基準にしたDR解の方が、精度が高まることが期待できる。したがって、本実施の形態は、測位の精度を向上させることができる。
また、本実施の形態では、所定時間連続してフィックス解が得られている状態において、フロート解が得られた場合には、DR解を出力する。また、この場合に出力されるDR解は、所定時間連続してフィックス解が得られた後の最新のフィックス解を基準にして求められたものであり、この最新のフィックス解は、正しい座標を示している可能性が高い。そのため、このような状況では、フロート解よりもDR解の方が、精度が高まることが期待できる。したがって、本実施の形態は、高精度のDR解を出力することができるので、測位の精度を向上させることができる。
また、本実施の形態では、DR解が第2時間連続して出力された場合にはRTK測位解を出力する。DR解は、基準となる座標に対して速度を加算していくこと等により算出されるので、時間の経過と共に誤差が蓄積し、精度が劣化していく。そのため、長時間経過した後は、たとえRTK測位解が十分に信頼できない状態であっても、DR解の精度の方がRTK測位解以上に、劣化している可能性が高い。したがって、このような場合には出力する解をRTK測位解に切り替えることにより、測位精度の劣化を低減することができる。
なお、本開示は、部材の種類、配置、個数等は前述の実施の形態に限定されるものではなく、その構成要素を同等の作用効果を奏するものに適宜置換する等、発明の要旨を逸脱しない範囲で適宜変更することができる。
例えば、上記の実施の形態では、第1時間連続してフィックス解が得られている状態において、フィックス解の速度差ベクトルの大きさ|A|が閾値sより小さい場合には常に該フィックス解を出力する場合について説明したが、本開示はこれに限られず、第1時間連続してフィックス解が得られている状態においてDR解を出力した際に、第3時間連続して、フィックス解の速度差ベクトルの大きさ|A|が閾値sより小さくなった場合に該フィックス解を出力し、それまではDR解を出力するようにしても良い。なお、DR解の出力をやめ、フィックス解の出力に切り替える条件は、これに限られるものではない。例えば、ユーザが自らの入力により出力を切り替えられるようにしてもよい。このような構成は、例えばDR解で算出される座標がユーザの目から見ても明らかに劣化している場合に有用である。
また、上記の実施の形態では、誤フィックスを検知するために速度差ベクトルの大きさを用いる場合について説明したが、本開示はこれに限られるものではなく、誤フィックスの発生を推測できる他の情報を用いても良い。例えば、速度ベクトルの角度の変化を調べることでも誤フィックスの発生を推測できる。角度を使用する場合は、速度ベクトルの方向が所定の角度以上(例えば180度)回転したことを検知することなどが考えられる。また、誤フィックスの判定の精度を向上させるため、速度差と角度を併用するなど複数の基準を複合的に用いても構わない。上記の変形例において、第3時間連続して速度差ベクトルの大きさの小さいフィックス解が得られた場合に、DR解の出力をやめてフィックス解の出力に戻す場合も、速度差ベクトルに替えて角度の変化などを用いてもよい。
また、上記の実施の形態において、現在出力されている解が、RTK演算による解とDR解のどちらであるのかを測位端末20がユーザに通知するようにしてもよい。通知の方法としては、例えば、移動体の表示部の画面(図示せず)に表示したり、音声で通知したりするなど、様々なものが考えられる。
また、上記の実施の形態では、測位演算の一例としてRTK演算を行う場合について説明したが、本開示はこれに限られずRTK演算以外の測位演算を行っても良い。RTK演算以外の測位演算を行う場合、衛星航法システム以外の測位演算を用いても良い。具体的には、衛星航法を利用しにくい屋内などにおいては壁などに設置されたマーカー等をステレオカメラによって撮影して視差から位置を逆算することなどが考えられる。すなわち、本開示では、誤りが発生していなければ高精度に位置を推定可能な測位方法であれば、どのような方法を用いても良い。
上記の実施の形態の説明に用いた各機能ブロックは、典型的には集積回路であるLSIとして実現される。これらは個別に1チップ化されてもよいし、一部または全てを含むように1チップ化されてもよい。ここでは、LSIとしたが、集積度の違いにより、IC、システムLSI、スーパーLSI、ウルトラLSIと呼称されることもある。
また、集積回路化の手法はLSIに限るものではなく、専用回路または汎用プロセッサで実現してもよい。LSI製造後に、プログラムすることが可能なFPGA(Field Programmable Gate Array)、又は、LSI内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサーを利用してもよい。
さらには、半導体技術の進歩または派生する別技術によりLSIに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行ってもよい。バイオ技術の適用等が可能性としてありえる。
本開示は、衛星からの信号を利用して干渉測位を行う場合に用いるに好適である。
1 測位システム
10 基準局
20 測位端末
101、201 プロセッサ
102、202 記憶部
103、203 入力部
104、204 出力部
105、205 通信部
106、206 受信部
110、210 バス

Claims (7)

  1. 測位端末が複数の衛星から送信される情報に基づいて測位演算を行うことで、移動体の座標を決定する測位方法であって、
    前記測位端末は、
    前記測位演算によって得られる解であるフロート解または前記フロート解より精度の高いフィックス解を算出し、
    前記フィックス解の示す座標を基準として、前記移動体の座標の推測値であるDR(Dead Reckoning)解を算出し、
    前記フィックス解が第1時間連続して得られている状態において、前記測位演算によって新たにフィックス解が得られた場合、前記新たに得られたフィックス解と前回のフィックス解との差が閾値以上の場合には前記DR解を前記移動体の座標として出力する、
    測位方法。
  2. 前記フィックス解が前記第1時間連続して得られていない場合、前記測位演算によって新たに得られた解を、当該解がフロート解であるかフィックス解であるかにかかわらず、前記移動体の座標として出力する、
    請求項1に記載の測位方法。
  3. 前記差は、前記新たに得られたフィックス解に基づく現在の速度ベクトルと前記前回のフィックス解に基づく速度ベクトルとの差分である速度差ベクトルの大きさである、
    請求項1または2に記載の測位方法。
  4. 前記フィックス解が前記第1時間連続して得られている状態において、前記測位演算によって前記フィックス解が新たに得られず、前記フロート解が得られた場合には、前記フロート解と前記前回のフィックス解との差にかかわらず、前記DR解を前記移動体の座標として出力する、
    請求項1から3のいずれか一項に記載の測位方法。
  5. 前記DR解が第2時間連続して出力された場合には前記測位演算によって得られる解を、当該解がフロート解であるかフィックス解であるかにかかわらず、前記移動体の座標として出力する、
    請求項1から4のいずれか一項に記載の測位方法。
  6. 前記フィックス解が第1時間連続して得られている状態においてDR解が出力された際に、第3時間連続して前記新たに得られたフィックス解と前記前回のフィックス解との差が前記閾値より小さい場合には前記フィックス解を前記移動体の座標として出力する、
    請求項1から5のいずれか一項に記載の測位方法。
  7. 複数の衛星から送信される測位信号を受信する受信部と、
    前記測位信号に含まれる情報に基づいて測位演算を行うことで、移動体の座標を決定するプロセッサと、
    を具備し、
    前記プロセッサは、
    前記測位演算によって得られる解であるフロート解または前記フロート解より精度の高いフィックス解を算出し、
    前記フィックス解の示す座標を基準として、前記移動体の座標の推測値であるDR(Dead Reckoning)解を算出し、
    前記フィックス解が第1時間連続して得られている状態において、前記測位演算によって新たにフィックス解が得られた場合、前記新たに得られたフィックス解と前回のフィックス解との差が閾値以上の場合には前記DR解を前記移動体の座標として出力する、
    測位端末。
JP2019523397A 2017-06-09 2018-04-27 測位方法および測位端末 Active JP7038281B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2017114627 2017-06-09
JP2017114627 2017-06-09
PCT/JP2018/017139 WO2018225421A1 (ja) 2017-06-09 2018-04-27 測位方法および測位端末

Publications (2)

Publication Number Publication Date
JPWO2018225421A1 JPWO2018225421A1 (ja) 2020-05-21
JP7038281B2 true JP7038281B2 (ja) 2022-03-18

Family

ID=64603938

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019523397A Active JP7038281B2 (ja) 2017-06-09 2018-04-27 測位方法および測位端末

Country Status (5)

Country Link
US (1) US11047992B2 (ja)
EP (1) EP3637146A4 (ja)
JP (1) JP7038281B2 (ja)
CN (1) CN110678781B (ja)
WO (1) WO2018225421A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10830902B2 (en) * 2018-06-21 2020-11-10 Unicore Communications, Inc. Real-time kinematic using standalone global navigation satellite system receiver
CN111077546A (zh) * 2019-12-26 2020-04-28 广州比逊电子科技有限公司 应用于无人机的定位跳点处理方法以及接收机设备
CN114675313B (zh) * 2020-12-24 2025-09-26 嘉兴国电通新能源科技有限公司 修正rtk移动站差分方程的浮点解坐标的解算方法
CN115808702B (zh) * 2022-11-25 2025-07-25 重庆宝图科技发展有限公司 一种基于rtk的自动测量方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000029520A (ja) 1998-07-08 2000-01-28 Fuji Heavy Ind Ltd 自律走行方法及び自律走行車
JP2002004261A (ja) 2000-04-20 2002-01-09 Shimizu Corp 盛土締固め管理システム
JP2008128793A (ja) 2006-11-20 2008-06-05 Toyota Motor Corp 移動体位置測位装置
JP2013122406A (ja) 2011-12-09 2013-06-20 Tokyo Keiki Inc 衛星航法システムを用いた位置出力装置
US20140375497A1 (en) 2013-06-25 2014-12-25 Caterpillar Inc. Positioning Error Detection and Mitigation System and Method
WO2015168460A1 (en) 2014-05-02 2015-11-05 Trimble Navigation Limited Dead reckoning system based on locally measured movement
WO2015186845A1 (ja) 2015-06-29 2015-12-10 株式会社小松製作所 作業機械の制御システム及び作業機械の制御方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07104847A (ja) 1993-09-30 1995-04-21 Fuji Heavy Ind Ltd 自律走行車の走行制御装置
JP3338304B2 (ja) * 1996-10-01 2002-10-28 横河電子機器株式会社 航法装置
JP5586278B2 (ja) 2010-03-12 2014-09-10 三菱電機株式会社 測位データ評価装置、測位データ評価方法および測位データ評価プログラム

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000029520A (ja) 1998-07-08 2000-01-28 Fuji Heavy Ind Ltd 自律走行方法及び自律走行車
JP2002004261A (ja) 2000-04-20 2002-01-09 Shimizu Corp 盛土締固め管理システム
JP2008128793A (ja) 2006-11-20 2008-06-05 Toyota Motor Corp 移動体位置測位装置
JP2013122406A (ja) 2011-12-09 2013-06-20 Tokyo Keiki Inc 衛星航法システムを用いた位置出力装置
US20140375497A1 (en) 2013-06-25 2014-12-25 Caterpillar Inc. Positioning Error Detection and Mitigation System and Method
WO2015168460A1 (en) 2014-05-02 2015-11-05 Trimble Navigation Limited Dead reckoning system based on locally measured movement
WO2015186845A1 (ja) 2015-06-29 2015-12-10 株式会社小松製作所 作業機械の制御システム及び作業機械の制御方法

Also Published As

Publication number Publication date
CN110678781A (zh) 2020-01-10
WO2018225421A1 (ja) 2018-12-13
US11047992B2 (en) 2021-06-29
CN110678781B (zh) 2023-07-21
EP3637146A1 (en) 2020-04-15
EP3637146A4 (en) 2020-05-27
US20200116873A1 (en) 2020-04-16
JPWO2018225421A1 (ja) 2020-05-21

Similar Documents

Publication Publication Date Title
JP7153842B2 (ja) 測位方法および測位端末
JP6952250B2 (ja) 測位方法および測位端末
JP6523526B2 (ja) 測位装置
JP7065277B2 (ja) 測位方法および測位端末
US10816675B2 (en) Coordinate output method and coordinate output device
JP7038281B2 (ja) 測位方法および測位端末
JPWO2015145718A1 (ja) 測位装置
US11294072B2 (en) Method, device and server for estimation of IFB calibration value
WO2018110025A1 (ja) 測位方法および測位端末
WO2020004538A1 (ja) 構造物監視サーバおよび構造物監視システム
WO2018110026A1 (ja) 測位方法、配信方法、測位端末および測位システム

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210121

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210121

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210907

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211104

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220118

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220131

R151 Written notification of patent or utility model registration

Ref document number: 7038281

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151