观测站结果:30天波罗的海历史数据与敏感性研究
版本kshana-observatory/adsb/v1,建立在kshana-interference-map/adsb/v2之上(参数已在基础方法中预先登记,未作更改)。MODELLED:没有外部参考基准(oracle),没有地面真值,也不作成因归属。“降级”单元格表示该日在该处采样的飞机中有较高比例报告了较低的导航精度,但并不说明原因。未观测或缺失的单元格表示未作评估,这与“无异常”并不相同。
运行内容#
- 数据源:adsb.lol readsb历史数据(ODbL 1.0),
adsblol/globe_history_2026的GitHub发布版v<YYYY>.<MM>.<DD>-planes-readsb-prod-0,由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(依次为南、西、北、东:波罗的海及其沿岸),30个UTC日,自2026-09-15至2026-10-14,0.5度单元格,由
scripts/observatory_backcatalogue.sh处理(逐日进行,每天处理完毕后删除原始文件)。 - 发布规则:最新一天(2026-10-14)在生成时已过其结束时刻24小时以上;否则命令拒绝运行。区域运行仅根据范围框内的报告计算设备基线和当日背景(
method.scope)。 - 文件:120个产品文件(GeoJSON、PNG、事件JSON,以及每日的清单)不在git中。每个文件的哈希值都在
back-catalogue/MANIFEST.json中,可以根据某一天的归档重建该日并进行比对。
某一天的读数示例:2026-09-15当天,归档中有110,619,797条位置记录,分布在79,826个轨迹文件中;其中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为226个,对比694个,因为每个单元格的面积大四倍):单元格大小是一项方法选择,观测台将其固定为0.5度。
- 没有已登记降级单元格的一天(10-03),在22个非网格变体中有9个出现了降级单元格(使用
low_nic_max为6时最多达8个):“没有降级单元格”同样是一项取决于阈值的陈述。
没有任何变体是经独立外部参考基准核对过的替代方案,已登记的规则也不因本研究而改变。本研究不说明哪些参数取值是正确的。