観測所の結果:30 日間のバルト海の過去分と感度の研究
バージョン kshana-observatory/adsb/v1、kshana-interference-map/adsb/v2 上(パラメーターは基本手法で事前登録済みで、変更なし)。MODELLED:外部のオラクルも、真値も、原因の特定もありません。劣化セルとは、その日にそのセルでサンプリングされた航空機のうち、航法精度が低いと報告した機体の割合が高かったことを示すだけで、その理由は示しません。観測されていない、または存在しないセルは評価されていないことを意味し、問題がないこととは異なります。
実行した内容#
- ソース:adsb.lol の readsb 履歴(ODbL 1.0)、GitHub リリース
v<YYYY>.<MM>.<DD>-planes-readsb-prod-0(adsblol/globe_history_2026のもの)。scripts/fetch_adsblol_day.shで一パートずつ取得します(一日あたり約 3.4 GB、2 パート)。公開元はチェックサムを示していないため、各パートの SHA-256 をback-catalogue/fetch/<tag>.parts.jsonに記録しています。 - 地域と日数:範囲は 53.5、9.0、66.0、31.0(南、西、北、東の順。バルト海とその沿岸)、2026-09-15 から 2026-10-14 までの UTC 30 日分、セルは 0.5 度、
scripts/observatory_backcatalogue.shを使用(一日ずつ処理し、日ごとに処理後に生ファイルを削除)。 - 公開ルール:最新の日(2026-10-14)は、生成時に日の終了から 24 時間を超えて経過していました。そうでなければコマンドは実行を拒否します。地域単位の実行では、機器ベースラインとその日のバックグラウンドを、範囲内の報告のみから計算します(
method.scope)。 - ファイル:120 個の成果物ファイル(GeoJSON、PNG、インシデント JSON、日ごとのマニフェスト)は git に含まれていません。各ファイルのハッシュは
back-catalogue/MANIFEST.jsonにあり、アーカイブから一日分を再構築して比較できます。
一日分の読み方の例:2026-09-15 のアーカイブには 79,826 個のトレースファイルに 110,619,797 行の位置データがありました。うち 103.0 百万行は範囲外、5.6 百万行は ADS-B ではなく、1.2 百万行は地上または 5,000 ft 未満、0.6 百万行は精度フィールドを持たず、212,426 行が集計されました(manifest.json の入力カウンター)。
再現性#
三日分(2026-09-18、2026-09-26、2026-10-03)はアーカイブから再取得して 二回目の実行を行いました。GeoJSON、PNG、インシデントファイルは 一回目の実行とバイト単位で同一(SHA-256 が一致)であり、識別子のハッシュソルトが実行ごとにランダムであっても、成果物は決定的です。
バックカタログの概要#
| UTC 日 | 公開セル数 | 劣化 | 劣化なし | 未観測(公開) | インシデント |
|---|---|---|---|---|---|
| 2026-09-15 | 643 | 7 | 461 | 175 | 3 |
| 2026-09-16 | 686 | 6 | 472 | 208 | 5 |
| 2026-09-17 | 691 | 5 | 460 | 226 | 2 |
| 2026-09-18 | 694 | 8 | 473 | 213 | 5 |
| 2026-09-19 | 663 | 2 | 416 | 245 | 2 |
| 2026-09-20 | 673 | 1 | 464 | 208 | 1 |
| 2026-09-21 | 674 | 3 | 465 | 206 | 3 |
| 2026-09-22 | 674 | 4 | 455 | 215 | 3 |
| 2026-09-23 | 668 | 2 | 456 | 210 | 2 |
| 2026-09-24 | 683 | 3 | 468 | 212 | 2 |
| 2026-09-25 | 701 | 3 | 469 | 229 | 3 |
| 2026-09-26 | 656 | 1 | 423 | 232 | 1 |
| 2026-09-27 | 682 | 2 | 459 | 221 | 1 |
| 2026-09-28 | 683 | 2 | 472 | 209 | 2 |
| 2026-09-29 | 666 | 3 | 454 | 209 | 2 |
| 2026-09-30 | 687 | 6 | 465 | 216 | 4 |
| 2026-10-01 | 691 | 1 | 476 | 214 | 1 |
| 2026-10-02 | 684 | 4 | 462 | 218 | 3 |
| 2026-10-03 | 650 | 0 | 397 | 253 | 0 |
| 2026-10-04 | 675 | 5 | 458 | 212 | 4 |
| 2026-10-05 | 689 | 2 | 466 | 221 | 1 |
| 2026-10-06 | 662 | 4 | 448 | 210 | 3 |
| 2026-10-07 | 671 | 3 | 458 | 210 | 2 |
| 2026-10-08 | 664 | 4 | 447 | 213 | 4 |
| 2026-10-09 | 687 | 4 | 465 | 218 | 2 |
| 2026-10-10 | 640 | 3 | 402 | 235 | 2 |
| 2026-10-11 | 678 | 3 | 459 | 216 | 2 |
| 2026-10-12 | 688 | 4 | 467 | 217 | 4 |
| 2026-10-13 | 672 | 0 | 457 | 215 | 0 |
| 2026-10-14 | 669 | 1 | 474 | 194 | 1 |
30 日間:劣化セルは合計 96 個、少なくとも一個あった日は 30 日中 28 日、インシデント(接する劣化セルのまとまり)は 70 件で、日全体に影響が及んだ日はありません。影響を受けた航空機の割合のセル中央値(その日のバックグラウンド)は 29 日で 0、一日で 0.006 だったため、バックグラウンドガード(セルはバックグラウンドを 0.15 上回る必要がある)がこの地域で判定を変えたことはありません。一日あたりの劣化セル数は 0 から 8 の範囲です(平均 3.2)。公開セルの約三分の一は毎日 not_observed です(サンプリングされた航空機が少なすぎる)。この成果物の限界は空ではなく、カバレッジです。
これが示さないこと:劣化セルが多い日は、物理的な意味で「悪い」わけではありません。持続的に影響を受ける領域は過小評価されます(終日その内側だけを飛ぶ航空機は、その外側にベースラインがなく、カウントされません)。また、ここからは、誰が、または何が精度低下を引き起こしたかは分かりません。
感度調査(三日分)#
kshana interference-map sensitivity は、ある一日分の記録済み報告を、登録済みパラメーターを 一つずつ変えて再生します(各二つの値と 1.0 度のグリッド)。対象の日:2026-09-18(劣化セルが最多の 8 個)、2026-09-26(劣化セル 1 個)、2026-10-03(0 個)。完全なレポート:sensitivity/<day>.json(各 23 通りのバリアント。再生した行数はそれぞれ 2.2、1.7、1.7 百万行)。各バリアントでの劣化セル数(括弧内は登録値。登録済みの結果は劣化セル 8 個、1 個、0 個):
| パラメーター(登録値) | 値 | 09-18 (8) | 09-26 (1) | 10-03 (0) |
|---|---|---|---|---|
low_nacp_max (6) |
5 / 7 | 6 / 8 | 1 / 3 | 0 / 2 |
low_nic_max (5) |
4 / 6 | 3 / 16 | 0 / 12 | 0 / 8 |
good_nacp_min (8) |
7 / 9 | 8 / 6 | 1 / 0 | 0 / 1 |
min_alt_ft (5000) |
3000 / 8000 | 6 / 6 | 1 / 1 | 0 / 1 |
affected_report_share (0.5) |
0.4 / 0.6 | 8 / 4 | 2 / 1 | 2 / 0 |
degraded_share (0.30) |
0.2 / 0.4 | 14 / 1 | 5 / 1 | 4 / 0 |
min_aircraft_sampled (10) |
8 / 15 | 8 / 8 | 2 / 0 | 2 / 0 |
min_reports_per_aircraft_in_cell (3) |
2 / 5 | 11 / 2 | 12 / 0 | 6 / 1 |
background_margin (0.15) |
0.10 / 0.20 | 8 / 8 | 1 / 1 | 0 / 0 |
min_baseline_good_reports (5) |
3 / 10 | 8 / 8 | 1 / 1 | 0 / 0 |
day_confound_background (0.15) |
0.10 / 0.25 | 8 / 8 | 1 / 1 | 0 / 0 |
| グリッド 0.5 度 | 1.0 | 2(粗いセル) | 3 | 1 |
勝者を決めずに読み取ると:
- 判定は精度コードのしきい値とサンプリング規則に敏感で、件数が少ないため大きく変動します。
low_nic_max、degraded_share、min_reports_per_aircraft_in_cellは、いずれの日でも劣化セル数を 二倍以上変化させます。low_nic_maxを 6 に緩めると、09-26 では劣化セル 1 個が 12 個になります(登録済みの集合との Jaccard は 0.08)。09-18 でdegraded_shareを 0.4 に上げると、8 個のうち 1 個が残ります(Jaccard 0.125)。しきい値付近のセル判定は、データについてであると同時に、しきい値についての記述でもあります。 min_aircraft_sampledは主に、セルを「劣化なし」と「未観測」の間で移します:09-18 では劣化セル数が 8 のままで、21 から 56 個のセルの状態が変わります。- 三つのガードはここでは作動していません:
background_margin、min_baseline_good_reports、day_confound_backgroundは、三日のいずれでも何も変えません(バックグラウンドの割合は 0)。したがって、この地域と期間からは、これらがどの程度うまく機能するかは分かりません。地域的な影響でバックグラウンドが上がる場合には、重要になります。 - より粗いグリッド(1.0 度)は様子を変えます(劣化セルは 2 個、3 個、1 個。各セルが 四倍の大きさになるため、公開セルは約三分の一に減り、09-18 では 694 に対して 226):セルサイズは手法上の選択であり、オブザーバトリーは 0.5 度に固定しています。
- 登録済みの劣化セルがない日(10-03)でも、グリッド以外の 22 バリアントのうち 9 で劣化セルが出ます(
low_nic_maxを 6 にすると最大 8 個):「劣化セルなし」という判定も、しきい値に依存する記述です。
独立した外部の参照で確認された代替案となるバリアントはなく、この調査によって登録済みの規則は変わりません。この調査は、どのパラメーター値が正しいかを示すものではありません。