PNT Moni 月次QCレポート
GEONET 観測品質 — Free Tier(無料公開版)
方法論バージョン: 1.1.0 · QC ツール: teqc 2019Feb25 · convbin: RTKLIB 2.4.2 · 電子基準点数: 1,302 · 対象日数: 30 · https://pntmoni.com/methodology/qc/v1.1.0
1 改訂履歴
| バージョン | 日付 | 内容 |
|---|---|---|
| 1.0 | 2026-07-08 | 初版。 |
| 1.1 | 2026-07-10 | 視覚デザイン刷新(ブランド統一、ADR 0017)。数値変更なし。PDF は新表紙(Typst)に移行 |
| 1.2 | 2026-07-18 | 低 n セル表示改訂: 3 局未満セルの空白表示を撤廃し、点線リング付き全表示へ。方法論 qc v1.1.0(数値変更なし) |
2 はじめに
本レポートは、国土地理院(GSI)が運用する GEONET(GNSS 連続観測網) 電子基準点網の 観測品質 の状態を、2026年6月 について まとめたものです。PNT Moni モニタリングサービスの Free tier 提供物 であり、GNSS 補強サービスの事業者・研究者・利用者が、基盤となる観測 インフラの稼働状況を把握できるよう、全国の集約的な視点を提供します。
本レポートは対象範囲を意図的に限定しています。以下は 含みません:
- 電子基準点単位の詳細分析や電子基準点名リスト(Pro tier)。
- 測位領域の精度や TTFF 結果 — これらは完全版の月次レポートに掲載 されます。
- インテグリティ指標、電離層・対流圏補正の診断、エンジン間比較 (Phase 1 以降)。
本レポートの生成に用いたデータと手法は公開文書化されており、公開されている GEONET 入力から再現可能です。§評価手法 に要約を示します。 完全な方法論文書は https://pntmoni.com/methodology/qc/v1.1.0 を参照してください。
3 評価手法
本節は方法論の要約です。完全な方法論文書は https://pntmoni.com/methodology/qc/v1.1.0 を参照してください。
データソース(方法論 §2): GEONET RINEX 3.02 の観測ファイルおよび 放送暦ファイル(30 秒サンプリング、24 時間日次ファイル)。GSI の アーカイブから Public Data License 1.0 のもとで取得。2026年6月 に稼働していた連続観測中の GEONET 電子基準点 1,302 局すべてを 対象とします。
QC パイプライン(方法論 §3):
RINEX 3 (GEONET)
↓ convbin (RTKLIB 2.4.2) v3 → v2.11 conversion
↓ teqc 2019Feb25 +qc -R +L2C_L2 +L5
→ .{yy}S summary per station/day
↓ pntmoni-pipeline qc summarize parse + wide Parquet
→ data/processed/qc_summary/... 1 row × 639 cols per station-day
QZSS の LNAV 放送は IS-QZSS-PNT-006 §4.1.2.3(4) の 2024 年以降のヘルスフラグ符号化を使用します。本パイプラインは QNAV を teqc に渡す前に L1C/B ヘルスビットをマスクし、L1C/A 放送中の QZS 衛星を 正しく QC します(方法論 §2.1)。
指標(方法論 §4): Shiono & Kubo (2025)(ION GNSS+ 2025)で確立し、 Shiono & Kubo (2026)(NAVIGATION, doi:10.33012/navi.762)で拡張した 4 つのコア指標ファミリ — 可視率・SNR・マルチパス・サイクルスリップ — にデータアベイラビリティを加えたものです。Galileo E5b 系(SN7、MP17、 MP71)は参考情報としてのみ報告し、絶対評価ベースラインには含みません。 記号ごとの定義は 付録の指標用語集 を参照してください。
集計(方法論 §5): マルチパス・SNR・スリップの値はすべて teqc の 5° 仰角ビンから、仰角マスク 15° を適用して算出します (19 ビン中 15 ビンを採用)。マルチパス・SNR は採用ビン間で観測数重み付け 平均、スリップは合算です。これは Shiono & Kubo (2025) の手法と、測位 エンジン(CLASLIB、RTKLIB)の仰角マスクの双方に一致します。
3.1 評価ネットワーク
図 1: 本レポートが対象とする GEONET 電子基準点。
4 エグゼクティブサマリ
| 指標 | 2026年6月 |
|---|---|
| 全国データアベイラビリティ(エポック被覆率・日次平均) | 99.97 % |
| 全国可視率(仰角10°超のトラック完全性・日次平均) | 95.60 % |
| L1 C/A マルチパス — 基準局・月間の中央値 | 0.201 m |
| L1 C/A SNR — 基準局・月間の中央値 | 43.8 dB-Hz |
| 報告電子基準点数 | 1,302 |
| 対象日数 | 30 |
要点。 GEONET 網は 2026年6月 を通じて高い観測忠実度で 稼働し、主要指標は健全な基準局網に期待される水準で推移しました。これら サマリ数値の背後にある分布と傾向は §全国概況 を参照してください。
5 全国概況 — 観測品質
5.1 データアベイラビリティと可視率
5.1.1 データアベイラビリティ
エポックレベルの被覆率: 各局は 30 秒ごとにデータを記録していたか?
図 2: データアベイラビリティ — 1日あたり公称 2,880 エポックのうち観測が記録されたエポックの割合。
| 統計量 | 基準局・月間のデータアベイラビリティ |
|---|---|
| 5 パーセンタイル | 100.00 % |
| 中央値 | 100.00 % |
| 95 パーセンタイル | 100.00 % |
| 最良局 | 100.00 % |
| 最悪局 | 12.50 % |
5.1.2 可視率
仰角10°超の衛星トラック完全性: 局が稼働していたとき、本来見えるはず の衛星を実際に追尾できていたか?
図 3: 可視率 — teqc による仰角10°超の Complete obs / Possible obs。
| 統計量 | 基準局・月間の可視率 |
|---|---|
| 5 パーセンタイル | 86.61 % |
| 中央値 | 97.10 % |
| 95 パーセンタイル | 99.71 % |
| 最良局 | 100.00 % |
| 最悪局 | 75.28 % |
5.2 指標の分布
5.2.1 信号対雑音比(SNR)
図 4: SNR 分布 — L1 C/A と L5(Shiono & Kubo 2025 のコアベースライン)。
図 5: SNR 分布 — L2(コア)と Galileo E5b(参考。絶対評価ベースラインには含まれない)。
5.2.2 マルチパス
図 6: マルチパス RMS 分布 — L1/L2 ペア(MP12 = L2 をパートナーとする L1 マルチパス、MP21 = L1 をパートナーとする L2 マルチパス)。Shiono & Kubo 2025 のコアベースライン。
図 7: マルチパス RMS 分布 — L1/L5 ペア(MP15 = L5 をパートナーとする L1 マルチパス、MP51 = L1 をパートナーとする L5 マルチパス)。コアベースライン。
図 8: 参考マルチパス RMS 分布 — L1/Galileo E5b ペア(MP17/MP71)。L1/E5b 周波数ペアの係数が大きいため数値スケールも大きい。絶対評価ベースラインには含まれない(E5b は局間ばらつきが大きく、Shiono & Kubo 2025 では評価指標として採用されなかった)。
5.2.3 サイクルスリップ
図 9: 基準局・月間のサイクルスリップ(ION とコア MP 組合せのスリップを仰角マスク15°のビンで合算)。論文準拠の定義(Shiono & Kubo 2025 Table 1)。対数ビン。
5.3 空間分布
以下の地図は、1,302 局の基準局・月間値を 全国の六角形グリッドに集約したものです。各六角形は 対辺間 約60 km で、 セル内に含まれる局の基準局・月間値の 中央値 で着色しています。グリッドは Albers 正積図法上に描画され、全セルが同一の地表面積を覆うため、全国間・月間で 色を直接比較できます。報告局数が 3 局未満の セルは点線の内側リング付きで表示します — これは信頼性の表示であり、 局数の少ないセルの中央値は統計的に薄いため、参考値として読んでください (方法論 §5.3、2026-07-18 改訂)。
5.3.1 信号対雑音比(SNR)
図 10: 基準局・月間の SNR — L1 C/A と L5(コアベースライン)。信号ごとに固有のカラースケール。
図 11: 基準局・月間の SNR — L2(コア)と Galileo E5b(参考。絶対評価ベースラインには含まれない)。信号ごとに固有のカラースケール。
5.3.2 マルチパス
図 12: 基準局・月間のマルチパス RMS — L1/L2 標準ペア(MP12 = L2 をパートナーとする L1 マルチパス、MP21 = L1 をパートナーとする L2 マルチパス)。
図 13: 基準局・月間のマルチパス RMS — L1/L5 ペア(MP15 = L5 をパートナーとする L1 マルチパス、MP51 = L1 をパートナーとする L5 マルチパス)。MP51 は MP12 より低い傾向 — 近年の L5 信号は GEONET 局でマルチパスが少ない。
図 14: 参考マルチパス — L1/Galileo E5b ペア(MP17/MP71)。L1/E5b 周波数ペアの係数が大きいため数値スケールも大きい。モニタリング目的で表示。絶対評価ベースラインには含まれない。
5.3.3 サイクルスリップ
図 15: 基準局・日あたりのコアサイクルスリップ(ION + MP12/21/15/51 のスリップ、仰角マスク15°、月内の日次平均)。論文準拠の定義。対数カラースケール。
Free tier の範囲として、六角形集約(セル中央値)の値のみを表示します — 個々の局識別子や局単位の値は Pro tier レポートで提供します。カラースケールは 信号ごとに固定(付録の指標用語集 を参照)のため、月をまたいで 一目で比較できます。
5.4 受信機・アンテナ構成
図 17: ネットワーク全体での受信機×アンテナの組合せ上位10。
受信機×アンテナの組合せ上位10で、報告のあった 1,302 局のうち 74.2 % を占めます。組合せは全体で 46 種類存在します。
前月比の機材変更 は集約カウントのみで把握します(例: 「N 局が受信機を 変更、M 局がアンテナを変更」)。個々の局の同定は Pro tier レポートで提供 します。前月差はベースライン月が確立し次第表示します。
5.5 絶対局クオリフィケーション
Coming Soon。 絶対局クオリフィケーション — すなわち各 GEONET 局を Shiono & Kubo (2025) で定義された 4 つのコア QC 基準に対して合否判定 する処理 — は、移動窓インフラ(日次粒度の 3 ヶ月ローリング窓、 クオリフィケーション周期あたり局ごとに n ≈ 91 の日次値)が整い次第ここに 掲載します。クオリフィケーション基準は、3 ヶ月ローリング窓における日次 指標の 99.73 パーセンタイル(3σ)で、対象は以下の各項目です:
- 可視率 —
Complete obs / Possible obs > 10° - SNR — L1・L2・L5 の平均 SNR(
SN1、SN2、SN5) - マルチパス —
MP12、MP21、MP15、MP51の RMS - サイクルスリップ —
cycle_slips_core(ION + MP12/21/15/51 スリップ)
日次×3 ヶ月ローリング窓は、公開済みの Shiono & Kubo (2025) の週次× ローリング年手法からの意図的な運用上の逸脱であり、運用上の即時性 (直近3ヶ月の報告周期との整合)、環境適応、統計的頑健性の向上 (n ≈ 91 対 53)のために採用しています。
原著に従い、Galileo E5b ファミリ(SN7、MP17、MP71)は意図的に 除外 します:4 年ベースライン(2021–2024)の時点で L1/E5b 値は L1/L2/L5 ペアよりも局間ばらつきが大きく、合否しきい値として信頼できない ためです。これらは 参考・補助モニタリング としてレポートに残します。
本節の出力は全国の合格率(例: 「1,298 局中 1,055 局が 4 つのコア基準を すべて満たした」)で、集約カウントのみ です — 個々の局名は Pro tier レポートで提供します。
6 付録
6.1 GNSS システムイベント: NAGU / NANU / NAQU
2026年6月 中に GNSS コンステレーション運用に対して発出された 通知。観測異常や信号可用性との相互参照のために記録します。公式の NAGU(Galileo、EUSPA/GSC)・NANU(GPS、USCG-NAVCEN)・ NAQU(QZSS、QSS)の通知フィードより取得。
| 日付 | 衛星 | 通知 | イベント種別 | 重大度 |
|---|---|---|---|---|
| 2026-06-04 02:59 UTC | QZS PRN200 | NAQU 2026174, NAQU 2026184 | scheduled_maintenance | total |
| 2026-06-04 06:00 UTC | GPS PRN17 | NANU 2026045, NANU 2026046 | scheduled_maintenance | total |
| 2026-06-08 04:16 UTC | QZS PRN205 | NAQU 2026186, NAQU 2026188 | other | partial |
| 2026-06-08 04:17 UTC | QZS PRN195 | NAQU 2026185, NAQU 2026187 | other | partial |
| 2026-06-15 03:59 UTC | QZS PRN199 | NAQU 2026190, NAQU 2026191 | scheduled_maintenance | total |
| 2026-06-16 08:01 UTC | QZS | NAQU 2026193, NAQU 2026194 | other | partial |
| 2026-06-17 08:00 UTC | GAL | NAGU 2026035, NAGU 2026037 | scheduled_maintenance | informational |
| 2026-06-17 08:00 UTC | GAL | NAGU 2026036, NAGU 2026038 | scheduled_maintenance | informational |
| 2026-06-18 06:00 UTC | GAL | NAGU 2026034, NAGU 2026039 | scheduled_maintenance | informational |
| 2026-06-18 19:59 UTC | QZS PRN195 | NAQU 2026189, NAQU 2026195 | scheduled_maintenance | total |
| 2026-06-20 15:01 UTC | QZS | NAQU 2026196, NAQU 2026197 | other | partial |
| 2026-06-22 08:05 UTC | GAL | NAGU 2026040 | other | total |
| 2026-06-23 02:59 UTC | QZS PRN200 | NAQU 2026192, NAQU 2026198 | scheduled_maintenance | total |
| 2026-06-26 16:12 UTC | GPS PRN11 | NANU 2026048, NANU 2026049 | other | total |
| 2026-06-28 17:48 UTC | QZS PRN200 | NAQU 2026200, NAQU 2026202 | other | partial |
| 2026-06-28 17:48 UTC | QZS PRN210 | NAQU 2026199, NAQU 2026201 | other | partial |
2026年6月 に 16 件の通知を検出しました(公式の NAGU(Galileo)・NANU(GPS)・NAQU(QZSS)フィードより収集)。
6.2 地震一覧
Coming Soon。 GEONET 観測異常の解釈に関連しうる、対象期間中に 発生した地震をここに掲載します(気象庁(JMA)より取得)。ソースフィード 連携の実装待ちです。
6.3 衛星一覧
取得時点でのコンステレーション運用者の 公式ステータスページの スナップショット(期間サマリではなく、運用者ページ自体がライブ スナップショット)。出典:GPS = USCG NAVCEN、QZSS = QSS、Galileo = GSC Europa。GLONASS は本レポートの QC パイプラインから除外しています (方法論による)。ただし PNT Moni は GEONET 局で受信される Galileo 信号は 評価します。
取得時刻: 2026-07-07 11:13 UTC
6.3.1 GPS
| PRN | SVN | Block | Slot | Clock | Status | Notice |
|---|---|---|---|---|---|---|
| G01 | 80 | III | D2 | RB | Unusable | UNUSABLE |
| G02 | 61 | IIR | D1 | RB | Operational | FCSTSUMM |
| G03 | 69 | IIF | E1 | RB | Operational | FCSTDV |
| G04 | 74 | III | F4 | RB | Operational | FCSTDV |
| G05 | 50 | IIR-M | E3 | RB | Operational | FCSTSUMM |
| G06 | 67 | IIF | D4 | RB | Operational | FCSTDV |
| G07 | 48 | IIR-M | A4 | RB | Operational | FCSTSUMM |
| G08 | 72 | IIF | C3 | RB | Operational | FCSTSUMM |
| G09 | 68 | IIF | F3 | RB | Operational | FCSTSUMM |
| G10 | 73 | IIF | E2 | RB | Operational | FCSTSUMM |
| G11 | 78 | III | D5 | RB | Unusable | UNUSABLE |
| G12 | 58 | IIR-M | B4 | RB | Operational | FCSTSUMM |
| G13 | 43 | IIR | A5 | RB | Decommissioned | DECOM |
| G14 | 77 | III | B6 | RB | Operational | FCSTSUMM |
| G15 | 55 | IIR-M | F2 | RB | Operational | FCSTSUMM |
| G16 | 56 | IIR | B1 | RB | Operational | FCSTSUMM |
| G17 | 53 | IIR-M | C4 | RB | Operational | FCSTSUMM |
| G18 | 75 | III | D6 | RB | Operational | FCSTSUMM |
| G19 | 59 | IIR | C5 | RB | Operational | FCSTSUMM |
| G20 | 82 | III | F5 | RB | Operational | FCSTSUMM |
| G21 | 81 | III | E6 | RB | Operational | FCSTSUMM |
| G22 | 44 | IIR | B3 | RB | Operational | FCSTSUMM |
| G23 | 76 | III | E5 | RB | Operational | FCSTSUMM |
| G24 | 65 | IIF | A1 | RB | Unusable | UNUSABLE |
| G25 | 62 | IIF | B2 | RB | Operational | FCSTSUMM |
| G26 | 71 | IIF | B5 | RB | Operational | FCSTSUMM |
| G27 | 66 | IIF | C2 | RB | Operational | FCSTSUMM |
| G28 | 79 | III | A6 | RB | Unusable | UNUSABLE |
| G29 | 57 | IIR-M | C1 | RB | Operational | FCSTSUMM |
| G30 | 64 | IIF | A3 | RB | Operational | FCSTSUMM |
| G31 | 52 | IIR-M | A2 | RB | Operational | FCSTSUMM |
| G32 | 70 | IIF | F1 | RB | Operational | FCSTDV |
6.3.2 QZSS (Michibiki)
| PRN | SVN | Block | Signals | Status | Notice |
|---|---|---|---|---|---|
| J194 | 002 | Ⅱ-Q | L1C/A, L1C, L2C, L5, L1S, L5S, L6D, L6E | Operational | |
| J195 | 004 | Ⅱ-Q | L1C/A, L1C, L2C, L5, L1S, L5S, L6D, L6E | Operational | |
| J196 | 005 | ⅡA-Q | L1C, L2C, L5, L1C/B, L1S, L5S, L6D, L6E | Operational | |
| J199 | 003 | Ⅱ-G | L1C/A, L1C, L2C, L5, L1S, L5S, L1Sb, L6D, L6E, Sr/Sf | Outage | NAQU 2026203 |
| J200 | 007 | Ⅲ-G | L1C, L5, L1C/B, L5S, L1Sb, L6D, L6E | Operational |
6.3.3 Galileo
| E# | GSAT | Clock | Status | Notice |
|---|---|---|---|---|
| E01 | GSAT0210 | RAFS | Unusable | GENERAL 2023048 |
| E02 | GSAT0211 | PHM | Operational | |
| E03 | GSAT0212 | PHM | Operational | |
| E04 | GSAT0213 | RAFS | Operational | |
| E05 | GSAT0214 | PHM | Operational | |
| E06 | GSAT0227 | PHM | Operational | |
| E07 | GSAT0207 | PHM | Operational | |
| E08 | GSAT0208 | PHM | Operational | |
| E09 | GSAT0209 | RAFS | Operational | |
| E10 | GSAT0224 | PHM | Operational | |
| E11 | GSAT0101 | RAFS | Operational | |
| E12 | GSAT0102 | RAFS | Operational | |
| E13 | GSAT0220 | PHM | Operational | |
| E14 | GSAT0202 | RAFS | Unusable | GENERAL 2021008 |
| E15 | GSAT0221 | PHM | Operational | |
| E16 | GSAT0232 | PHM | Operational | |
| E18 | GSAT0201 | PHM | Unusable | GENERAL 2021008 |
| E19 | GSAT0103 | RAFS | Operational | |
| E21 | GSAT0215 | PHM | Operational | |
| E22 | GSAT0204 | RAFS | Unusable | GENERAL 2017045 |
| E23 | GSAT0226 | PHM | Operational | |
| E25 | GSAT0216 | PHM | Operational | |
| E26 | GSAT0203 | PHM | Operational | |
| E27 | GSAT0217 | PHM | Operational | |
| E28 | GSAT0233 | PHM | Operational | |
| E29 | GSAT0225 | PHM | Operational | |
| E30 | GSAT0206 | PHM | Operational | |
| E31 | GSAT0218 | PHM | Operational | |
| E32 | GSAT0234 | Commissioning | LAUNCH 2025061 | |
| E33 | GSAT0222 | PHM | Operational | |
| E34 | GSAT0223 | PHM | Operational | |
| E36 | GSAT0219 | PHM | Operational |
再現性。 本レポートはオープンソースの pntmoni-pipeline ツールチェーンにより、公開されている GEONET アーカイブから生成して います。方法論 1.1.0 · QC ツール teqc 2019Feb25 · convbin RTKLIB 2.4.2 — レポートと併せてバージョン管理されています。
6.4 指標用語集
本レポートで使用する指標の記号ごとの定義です(Shiono & Kubo 2025, 2026 で確立)。完全な方法論文書は https://pntmoni.com/methodology/qc/v1.1.0 を参照してください。
6.4.1 コア指標(絶対評価ベースライン)
論文の Table 1 の局選定に用いられた 4 ファミリであり、計画中の 絶対局クオリフィケーション の 基礎 です (§5.6 を参照)。
| 指標 | 記号 | 説明 |
|---|---|---|
| データアベイラビリティ | data_availability |
epochs_w_obs / expected_epochs — 1日の公称 2,880 個の 30 秒エポックのうち、局が少なくとも 1 つの観測を記録したエポックの割合。局が稼働・記録していたかを測ります。 |
| 可視率 | visibility |
Complete obs >10° / Possible obs >10° — teqc の仰角10°マスク超で本来得られるべき衛星エポックトラックのうち、実際に取得された割合。局が稼働していたときのトラック完全性を測ります。値が低い場合、毎エポック記録している局でもアンテナ遮蔽・RF 干渉・受信機の問題を示唆します。 |
| L1 C/A SNR | SN1 |
L1 C/A の平均 SNR(dB-Hz)。 |
| L2 SNR | SN2 |
L2(P コードまたは L2C)の平均 SNR(dB-Hz)。 |
| L5 SNR | SN5 |
L5 の平均 SNR(dB-Hz)。 |
| マルチパス(L1、L2 パートナー) | MP12 |
L2 をパートナーとする L1 マルチパス組合せの RMS(m)。 |
| マルチパス(L2、L1 パートナー) | MP21 |
L1 をパートナーとする L2 マルチパス組合せの RMS(m)。 |
| マルチパス(L1、L5 パートナー) | MP15 |
L5 をパートナーとする L1 マルチパス組合せの RMS(m)。 |
| マルチパス(L5、L1 パートナー) | MP51 |
L1 をパートナーとする L5 マルチパス組合せの RMS(m)。 |
| サイクルスリップ(コア) | cycle_slips_core |
teqc が報告する ION スリップと 4 つのコア MP 組合せスリップ(MP12 + MP21 + MP15 + MP51)の合計。論文準拠の定義。 |
6.4.2 参考・補助モニタリング(Galileo E5b)
これらは 参考情報 として報告します。L1/E5b 周波数ペアは L1/L2/L5 ペアよりも局間ばらつきが大きく、Shiono & Kubo (2025) では 評価指標として採用されなかったため、絶対評価ベースラインには 含みません。
| 指標 | 記号 | 説明 |
|---|---|---|
| Galileo E5b SNR | SN7 |
E5b(1207.140 MHz)の平均 SNR(dB-Hz)。 |
| マルチパス(L1、E5b パートナー) | MP17 |
E5b をパートナーとする L1 マルチパス組合せの RMS(m)。マルチパスの式における 2/(α-1) 係数が L1/E5b 周波数ペアで大きいため、数値も大きくなります。 |
| マルチパス(E5b、L1 パートナー) | MP71 |
L1 をパートナーとする E5b マルチパス組合せの RMS(m)。 |
| サイクルスリップ(拡張) | cycle_slips_extended |
cycle_slips_core に MP17 と MP71 のスリップを加えたもの。Galileo E5b を追尾する受信機の評価に有用ですが、クオリフィケーションには使用しません。 |