NovAtel PwrPak7 静止時性能分析

今回は山中湖で停止した状態時の性能を評価します。

BESTGNSSPOSA, BESTPOSA, INSPVAXA三者の状態をプロットしたものを以下に示します。BESTGNSSPOSAのデータがばらつきが大きく、他の両者の比較的に小さく、形状も似ているものになっています。

第1章 静的データの特性とシステムサンプリング構造の分析

NovAtel SPAN（Synchronized Position Attitude Navigation）システムの性能を評価する上で、データフローのサンプリングロジックと各ログ間の物理的な関連性を理解することは不可欠です。提供された3組のデータサンプルは、サンプル数において顕著な差を示しており、これはOEM7プラットフォームがリアルタイムナビゲーション解を処理する際の優先順位を反映しています。

1.1 サンプリング周波数とデータ規模の統計的相関

サンプル規模を確認すると、INSPVAXA ログは 28,604 サンプルであるのに対し、BESTPOSA と BESTGNSSPOSA はいずれも 2,861 サンプルとなっています [ユーザー入力]。この約 10:1 の比率は、NovAtel SPAN システムの標準的な出力設定と正確に一致しています。通常、慣性航法装置（INS）の複合解は、自動運転や精密制御に必要な高頻度な姿勢・位置情報を提供するため 100 Hz 以上で出力され、一方でGNSS関連の測位解はベースバンド処理ユニットの更新周期に同期した 10 Hz 程度に設定されるのが一般的です 。

静的テストにおいて、100 Hz という高頻度サンプリングは大量の冗長データを生成しますが、その核心的価値は、ジャイロスコープのランダムウォークや加速度計のバイアス不安定性といったIMUセンサーのノイズ特性を膨大なサンプル数で完全に捉えられる点にあります。以下に、静止状態における3つのログの主要な統計指標をまとめます。

統計指標	INSPVAXA (INS/GNSS 複合解)	BESTPOSA (複合/最適解)	BESTGNSSPOSA (純GNSS解)
サンプル数 (Samples)	28,604	2,861	2,861
CEP50 (m)	0.3279	0.3270	0.3101
CEP95 (m)	0.6918	0.6910	0.8366
RMS (m)	0.4221	0.4218	0.4221
1$\sigma$ 標準偏差 (m)	0.1778	0.1778	0.2168
2$\sigma$ 標準偏差 (m)	0.3556	0.3557	0.4336

1.2 静的測位性能の差異に関する予備的考察

上表の横断的な比較から、注目すべき統計学的現象が見て取れます。50%確率円誤差（CEP50）の指標では、純GNSS測位（BESTGNSSPOSA）の数値が 0.3101 m であり、慣導補強された複合解（0.3279 m）よりも僅かに優れています [ユーザー入力]。しかし、システムの堅牢性を表す 95%信頼区間（CEP95）では、複合解（0.6918 m）が純GNSS解（0.8366 m）よりも約 17.3% 向上しています [ユーザー入力]。

この「中心部はGNSS、周辺部はSPANが安定」という特性は、SPANフィルターの核心的な作用メカニズムを示しています。GNSS単独測位は、理想的な開口環境下では中心値付近に密集しますが、電離層のシンチレーション、衛星星座の切り替え、マルチパスなどの影響を瞬時に受けやすく、大きな「誤差外れ値（Outliers）」を発生させることがあります。SPANシステムはカルマンフィルターを通じてIMUの物理的制約を導入することで、短期的にはフィルターの収束過程やIMUのバイアスドリフトにより中心点が微細にオフセットされる可能性があるものの、GNSS特有の突発的な飛び（スパイク）を効果的に抑制し、軌跡全体を滑らかで物理的法則に即したものにします 。

第2章 NovAtel SPAN測位モードとアルゴリズムロジックの解析データ

差異の物理的な原因を深く探るためには、ログに記録された測位タイプ（Position Type）とその内部処理プロセスを解剖する必要があります。

2.1 INS_PSRSP モードの定義と動作原理提供されたログの生データにおいて、INSPVAXA と BESTPOSA はいずれも測位タイプを INS_PSRSP と報告しています 6。NovAtel OEM7 の技術定義によれば、INS_PSRSP は「最後に適用された位置更新がGNSS単独測位（SINGLE）解に基づいている慣導計算位置」を意味します 7。このモードでは、RTK基準局やPPP（精密単独測位）補正サービスの支援がないため、絶対位置精度は民用GNSS信号の固有誤差に制限されます。しかし、SPANの深結合アーキテクチャにより、受信機は生の慣性データと衛星の疑似距離（Pseudorange）観測値をフィルター層で直接融合させることができます。GNSS 観測層: 受信機は L1/L2 搬送波周波数で衛星信号を捕捉し、コード位相（Code Phase）を利用して衛星との距離を算出します。この時点の誤差源は主に約1～3mの電離層遅延、0.5mの対流層遅延、および衛星軌道誤差です 1。INS 推算層: IMUは極めて高い周波数（例：100 Hz）で加速度と角速度を計測します。システムはニュートンの運動法則に基づいて積分推算を行い、瞬時の位置・速度・姿勢（PVA）を算出します 9。カルマンフィルタ融合: GNSSが有効な場合、フィルターはGNSS位置を観測値として使用してINSのドリフトを修正し、IMUのバイアスパラメータをリアルタイムで推定します。静止状態が識別されると（ZUPT技術など）、システムはさらにゼロ速度制約を導入し、位置のドリフトを強制的に制限します 4。

2.2 静止下での平滑化メカニズム：ノイズ抑制と外れ値の除去BESTGNSSPOSA は純GNSSログであり、その出力位置は衛星観測値の幾何学的な解のみに基づいています。静止環境下では、衛星の高度角の変化や大気条件の微細な変動により、算出される緯度・経度が座標軸上で高周波のジッタ（振れ）を発生させます。対照的に、INSPVAXA が示すRMSと標準偏差の関係は、より高い統計的一貫性を体現しています。低精度測位モード（SINGLEやWAASなど）において、慣導が報告する標準偏差は、実際のGNSS精度よりも「楽観的」に見えることがあります。これは精度が実際にセンチメートル級に達したわけではなく、フィルターが「静止している物体がGNSS解のような高周波のランダム変動を起こすはずがない」という物理的判断に基づき、強力な平滑化（スムージング）を行っているためです 3。ユーザーのデータでは、INSPVAXA の 1$\sigma$ が 0.1778 m であるのに対し、BESTGNSSPOSA の 1$\sigma$ は 0.2168 m に達しています [ユーザー入力]。この差は、静止環境下で慣性デバイスが「ローパスフィルター」として機能し、GNSS解のランダムなホワイトノイズを除去して、短期間の測位結果をより安定させていることを明確に証明しています。

第3章 生ログフィールドの深度デコードと診断分析ASCII形式の生ログ行をビット単位でデコードすることで、テスト期間中の受信機の具体的な動作環境と状態を洞察できます。

3.1 ログヘッダーと時間システムの分析例：#INSPVAXA,FILE,0,62.0,FINESTEERING,2396,57620.300,02004000,471d,15341;... 6FINESTEERING: このフィールドは、受信機の内部時計が衛星信号を通じて高精度に同期され、時計のオフセットとドリフトが完全に制御下にあることを示しています。これは高精度な測位解を得るための大前提です 10。2396, 57620.300: GNSS週数が2396、週秒が57620.300であることを表します。この精密なタイムスタンプにより、IMUのサンプリングデータと衛星エフェメリスデータがナノ秒レベルで厳密にアライメントされます 2。

3.2 受信機ステータスワード（Receiver Status）：02004000これは32ビットの16進数マスクであり、ハードウェアの健康状態をリアルタイムで監視するために使用されます。Bit 14 (0x00004000) - Antenna Gain State: OEM7アーキテクチャにおいて、このビットが立っていることは通常、アンテナゲインの状態が正常であることを示します。受信機内蔵の自動利得制御（AGC）回路が、RFフロントエンドの信号強度が想定範囲内にあり、深刻な帯域外干渉やアンテナの断線・短絡が発生していないことをモニターしています 11。Bit 25 (0x02000000) - Version Bits: NovAtel受信機の系譜において、Bit 25とBit 26の特定の組み合わせはハードウェアプラットフォームの識別に使われます。Bit 25が1でBit 26が0であることは、OEM7シリーズ受信機の標準的な識別記号であり、コマンドセットやログフォーマットの下位互換性を保証します 13。総合判定: ステータスワード内の Bit 0（致命的なエラー）がセットされていないことから、受信機全体の動作ロジックは非常に堅牢であり、ハードウェアレベルのダウングレード警告は発生していないと判断できます 10。

3.3 拡張解決ステータス（Extended Solution Status）：0b0030c7このフィールドは、SPANフィルター内部の意思決定を深く理解するための鍵となります。Bit 0-3 (0x7): 通常、現在の計算収束フェーズや特定の測位ソースを示します。Bit 12-15 (0x3): 静的テストにおいて、この値はシステムが特定のキネマティック制約（静止制約やゼロ速度更新など）を適用していることを反映することがよくあります。Bit 22 (0x00400000): このビットが変動する場合、主/予備フィルター間の切り替えを示唆している可能性があります。SPANシステムは主フィルターとGNSS観測値の間に乖離が生じた際、PVA出力の連続性を確保するために自動的にバックアップフィルターを有効にする機能を備えています 4。全記録期間を通じてこのフィールドが高度に一致していることから、SPANフィルターは静止環境下で「安定巡航」状態に入っていることがわかります。

第4章 静的測位誤差分布の数学的検証と統計的推論

測位精度指標の深い分析には、確率分布理論に基づいた考察が必要です。一般に、水平面上の静的測位誤差はレイリー分布（Rayleigh Distribution）に従うと考えられています。

4.1 RMS と標準偏差（$\sigma$）の内在的ロジックユーザーの報告によると、INSPVAXA の RMS は 0.4221 m、1$\sigma$ は 0.1778 m です [ユーザー入力]。標準的な2次元測位誤差分析において、RMS（具体的には DRMS: Distance Root Mean Square）の計算式は以下の通りです。

南北方向と東西方向の標準偏差がほぼ等しいと仮定すると、$DRMS \approx 1.414 \times \sigma$ となります。計算すると、$1.414 \times 0.1778 = 0.251$ m です。問題点の分析: 実測の RMS (0.4221 m) は、標準偏差から推算される理論値 (0.251 m) よりも顕著に高くなっています。これは静的な単独測位の核心的な特徴である系統的偏差（バイアス）の存在を明らかにしています。標準偏差（精度/Precision）: データ点自体の平均（Mean Position）に対するバラツキの程度を反映します。RMS（正確度/Accuracy）: データ点が真値（True Position）からどれだけ離れているかを反映します 16。自律測位（SINGLE）モードでは、差分補正がないため、電離層の残差誤差や衛星軌道誤差によって測位の中心自体が真の座標から数十～数百センチメートル系統的にオフセットされます。RMSが高くなっていることは、システムが測位の「再現性」は保っているものの、その測位中心自体に固定的な偏差があることを示しています。

4.2 確率信頼区間（CEP）の収束性に関する論証正規分布とレイリー分布の変換理論に基づくと、異なるパーセンタイルの誤差半径は特定の倍数関係を満たす必要があります 2。CEP95 / CEP50 の理論比率 $\approx 2.08$2$\sigma$ / 1$\sigma$ の理論比率 $= 2.0$INSPVAXA の場合:実測 CEP95 / CEP50 = $0.6918 / 0.3279 \approx 2.109$ [ユーザー入力]実測 2$\sigma$ / 1$\sigma$ = $0.3556 / 0.1778 = 2.000$ [ユーザー入力]結論: 複合解の誤差分布は、理論上のガウス・レイリー分布モデルに完璧に適合しています。これはSPANフィルターのカルマンゲインが最適に調整されており、システムが高度に制御された物理的実体として動作していることを示しています。BESTGNSSPOSA（純GNSS）の場合:実測 CEP95 / CEP50 = $0.8366 / 0.3101 \approx 2.698$ [ユーザー入力]深度インサイト: 純GNSSの実測比率（2.7）は理論値（2.08）を大きく上回っています。これは数学的に、GNSS単独測位が静止下において深刻な「ファットテール効果（Fat-tail effect）」を持つことを証明しています。つまり、50%の点は集中しているものの、残りの5%の点が信号マルチパスや星座配置の変動により、大きな位置の飛散を起こしているということです。SPANシステムはこの比率を2.7から2.1に押し下げることで、これら「非物理的」な飛散点を合理的な範囲内に拘束することに成功し、測位の予測可能性を大幅に強化しています。

第5章 慣導姿勢観測と静的方位角維持能力

位置精度以外に、INSPVAXA が BESTGNSSPOSA に対して持つ最大の利点は、完全な姿勢（Attitude）情報を提供できる点であり、これは静的キャリブレーションのシナリオにおいて極めて重要です。

5.1 姿勢計算の精度と安定性生ログ 6 のPVAデータセグメントには、リアルタイムの姿勢値とその推定標準偏差が記録されています。ロール角 (Roll): -0.812357°、標準偏差 0.0248° 6ピッチ角 (Pitch): 6.662142°、標準偏差 0.0253° 6方位角 (Azimuth): 203.506416°、標準偏差 0.1170° 6ロールとピッチの精度が高い（約0.02°）のは、静止下において重力ベクトル（Gravity Vector）が極めて強力な垂直基準の制約を提供するためです。方位角の標準偏差（0.117°）はこれに比べてやや大きいものの、静的自律モードとしては極めて優秀なパフォーマンスです。運動による向心加速度がないため、方位角の解は主に地球自転分量に対するジャイロ感応（ジャイロコンパス）やGNSS航向補正に依存します 4。

5.2 多センサーシステムにおける姿勢データの価値LiDARやカメラ等の視覚センサーを搭載したシステムにとって、静止状態での INSPVAXA の精度は外参キャリブレーション（Extrinsic Calibration）の成否を直結させます。0.11°の方位角不確定性は、10メートル先での投影誤差がわずか約2センチメートルであることを意味します。純GNSS測位（BESTGNSSPOSA）では姿勢リファレンスを提供できないため、静止状態におけるSPANのPVA出力は、センサーの空間座標系マッピングを構築するための唯一の信頼できる情報源となります 19。

第6章 NovAtel 独自技術による静的性能への潜在的貢献

現在のログには明示されていませんが、OEM7プラットフォームの標準コンポーネントとして、GLIDEおよびSTEADYLINE技術が静的安定性の維持に重要な役割を果たしています。

6.1 GLIDE アルゴリズムの平滑化ロジックNovAtelのGLIDE技術は、搬送波位相（Carrier Phase）と疑似距離観測値を融合させることで、測位結果の「ジャンプ」を排除することを目指しています。メカニズム: GLIDEは搬送波位相を用いて位置の微細な変化（相対変位）を測定し、それを絶対的な疑似距離解に重ね合わせます 21。静的性能: 静止状態において、GLIDEはエポック間（Epoch-to-epoch）の位置変動を1センチメートル以内に抑えることができます 23。これが、ユーザーの INSPVAXA 散布図が短期間において極めて高い局部的な緻密さを見せる理由の一つです。

6.2 STEADYLINE モード遷移衛星信号の遮蔽や干渉が発生した際、STEADYLINEは異なる測位精度レベル間（単独測位から慣導推算のみへの転落など）の滑らかな遷移を担当します 24。"Maintain" または "Transition" モードを通じて、衛星星座の瞬時の入れ替わりによる座標系のオフセットを防ぎます。これは静的な長期モニタリングにおいて、堅牢性を保証するための重要な機能です。

第7章 静的テスト環境下におけるハードウェア健康状態の診断

6 に含まれる衛星可視性ログ SATVIS2A に基づき、現在の観測環境の「健康診断」を行うことができます。

7.1 衛星星座の利用率ログによれば、システムは最大32機の衛星を追尾・利用しており、GPS、GLONASS、QZSS、およびBeiDouを網羅しています 25。GPS: 追尾数は6～12機の間で変動し、PDOP値は極めて良好なレベルにあります 6。マルチシステム融合: 受信機はQZSS（準天頂衛星システム）と多数のBeiDou衛星を同時に使用しています 6。マルチ星座の導入によりDOP（幾何学的精度低下率）が顕著に改善されており、これが BESTGNSSPOSA のCEP50が0.31メートルという高い水準に達している要因の一つです 1。

7.2 信号強度と干渉の評価受信機ステータスワード 02004000 と FINESTEERING の時刻状態から、アンテナ端において深刻な電磁干渉（Jamming）やなりすまし（Spoofing）の不確実性はないと判断できます 12。Undulation（ジオイド高）が 40.70 メートルで固定されていることは、座標系変換パラメータが正しくロードされていることを示し、楕円体高が 988.4～988.9 メートルの範囲に収まっていることは、単独測位の垂直精度が水平精度の約1.5倍であるという一般的な期待値に合致しています 1。

第8章 SPAN システム静的性能の限界と最適化

パスSPANシステムは平滑度と堅牢性において単独GNSSをはるかに凌駕していますが、「自律単独測位」という技術的前提の下では、依然として絶対精度の「天井」に直面しています。

8.1 絶対精度と相対精度のパラドックスINS_PSRSP モードにおいて、システムは慣導を通じて極めて高い相対精度（例：10分以内に観察される位置の変動がデシメートル級に留まるなど）を維持できますが、その絶対精度は依然として地球座標系に対するアライメント精度に制限されます 3。もしアプリケーションのシナリオがセンチメートル級の絶対位置（地殻変動モニタリングなど）を要求する場合、現在の 0.42 メートルの RMS では不十分です。

8.2 推奨される最適化戦略搬送波位相差分（RTK）の導入: NTRIP差分ソースに接続することで、測位タイプを INS_PSRSP から INS_RTKFIXED へアップグレードできます。これにより、CEP50 は 0.3 メートルオーダーから 0.01 メートルオーダーへと飛躍的に向上します 1。静的予熱と運動アライメントの実施: 慣導システムは起動直後、状態が不安定です。正式なデータ収集の前に15～20分程度の静的な予熱を行い、さらに典型的な「8の字」や「S字」運動を行ってIMUの全自由度をアクティブにすることを推奨します。これにより姿勢角（特にAzimuth）の標準偏差を顕著に低減できます 4。IMU の剛性固定: 静的テストであっても、微小な構造振動がIMUに捉えられ積分される可能性があります。IMUブラケットと基準点の剛性接続を確実にし、測位以外の要因による誤差項を排除する必要があります 5。