为什么ADS-B不直接发经纬度pyModeS CPR压缩位置解码原理全解【免费下载链接】pyModeSPython decoder for Mode S and ADS-B signals项目地址: https://gitcode.com/gh_mirrors/py/pyModeSpyModeS 是一款用 Python 编写的 ADS-B / Mode-S 航空信号解码器。如果你拆过一条 ADS-B 位置报文会发现一个反直觉的事实飞机并不直接广播经纬度而是发两个 17 位的CPR 压缩位置数字。这篇文章用最少代码讲透 ADS-B CPR 压缩位置的原理以及 pyModeS 如何把这两个神秘数字还原成地图上的坐标。 先看一条真实的 ADS-B 位置报文ADS-B 的空中位置消息Type Code 9-18 气压高度 / 20-22 GNSS 高度载荷只有56 位。拆开来看字段位宽作用TCType Code5 位消息类型SS NIC3 位监视状态、完整性等级高度码 AC12 位气压/卫星高度T F2 位时间同步位 CPR 格式位0偶帧 / 1奇帧CPR 纬度17 位压缩后的纬度CPR 经度17 位压缩后的经度一条来自法国上空的真实偶帧消息8D40058B58C901375147EFD09357它的 CPR 纬度是17161CPR 经度是97623。这不是 49.8°N、6.0°E——那是解码后的结果。为什么不直接发经纬度算一笔账就明白了两个 32 位浮点经纬度需要64 位已经超过整个 56 位载荷高度、呼号等字段就没地方放了即便压缩成 24 位定点经纬度全球精度也只有约 5 公里远不够识别跑道和航路ADS-B 工作在 1090 MHz 共享频段每架飞机约 0.5~1 秒广播一次每个字节都要用在刀刃上。所以 DO-260B 标准设计了 CPRCompact Position Reporting紧凑位置报告把位置拆成34 位量化数字 1 位格式位单帧精度达到米级但代价是——单条消息在全球范围内有几十上百个候选位置。️ CPR 的精髓把地球画成一张网格CPR 的核心思想是给地球表面画经纬网格偶帧even纬度切成 60 条带每带高 360°/60 6°奇帧odd纬度切成 59 条带每带高 360°/59 ≈6.10°经度更讲究在赤道上切成 59 份每份约 6.1°越靠近极地经度线越收敛份数也越少——高纬度地区经度圈本身变短格子要变宽才能保持近似等距。这个份数函数就叫NL(纬度)赤道是 59±87° 纬度降为 2。于是 17 位原始值就是格子编号 × 格宽 格内小数位置单帧的 17 位精度6° ÷ 131072 ≈ 5 米。精度足够但不知道落在全球第几个格子——偶帧有 60 个纬度候选 × 若干经度候选。这就是为什么 ADS-B 接收机必须记住飞机飞机持续广播偶帧 奇帧交替到达。偶帧落在 60 分格、奇帧落在 59 分格两套网格的错位恰好把候选位置唯一确定下来——这是 CPR 最巧妙的一笔。 pyModeS 如何还原经纬度三条解码路径pyModeS 把 CPR 数学全部收敛在 src/pyModeS/position/_cpr.py 中提供三种解法按手头信息从少到多排列1️⃣ 偶/奇帧配对解码最常用无需任何参考飞机每次广播都切换 even/odd 格式。收到一对相邻的帧后用公式求纬度j floor(59 × cpr_lat_even − 60 × cpr_lat_odd 0.5)再结合 NL(纬度) 解出经度。实现见 airborne_position_pair它还内置了两道保险两帧解出的纬度若落在不同 NL 纬度带说明配对失败或算出 |纬度| 90° 的物理不可能结果直接判为失败返回 None——真实案例中曾有一架从阿姆斯特丹起飞的飞机差点被解到白令海上空见 tests/test_cpr.py 中的回归测试。⏱️ 配对窗口默认10 秒间隔太久的两帧可能已跨越纬度带边界不再可信。2️⃣ 参考点解码已知道大致位置如果你已经有一个粗略参考点上次位置、雷达航迹、接收机位置单条消息就能解码偶帧参考点在 180 海里内即可。例如给上面那条法国消息一个粗参考 (49.0, 6.0)立刻得到49.8241°N, 6.0679°E验证向量见 tests/test_cpr.py。pyModeS 的PipeDecoder会在配对成功后自动把该位置升级为本地参考之后每一帧都能立即解出坐标不必等奇偶配对。3️⃣ 地面滑行 CPR需要机场参考地面移动报文BDS 0,6TC 5-8用的是更密的90° 纬度带bds06.py 负责提取原始字段经度还有四个 90° 象限候选因此必须给一个 45 海里内的参考点。pyModeS 贴心地内置了机场 ICAO 码数据库src/pyModeS/position/_airports.py一个代码就能换参考点from pyModeS import decode # 真实 DF18 地面滑行报文参考机场 LFBO图卢兹 r decode(903a23ff426a4e65f7487a775d17, surface_refLFBO) print(r[latitude], r[longitude]) # 43.6264, 1.3747 新手上手三步解码出第一条 ADS-B 位置pip install pyModeS3from pyModeS import decode # 偶帧 奇帧一起喂进去CPR 对自动配对 results decode([ 8D40058B58C901375147EFD09357, # even 8D40058B58C904A87F402D3B8C59, # odd ], timestamps[1.0, 2.0]) print(results[1][latitude], results[1][longitude]) # 49.8175 6.0844不想写代码pyModeS 附带modes命令行工具逗号分隔的多条消息共享同一个状态机CPR 对自动解析modes decode 8D40058B58C901375147EFD09357,8D40058B58C904A87F402D3B8C59完整状态机TTL 淘汰、幻影帧过滤、DF20/21 ICAO 校验的深入说明见 docs/pipe.md入门教程见 docs/quickstart.md。 一句话总结ADS-B 不发经纬度是因为56 位载荷装不下全球米级精度的完整坐标。CPR 用60/59 分格错位的数学技巧把位置压缩到 34 位再用飞机自己的连续广播偶帧奇帧或一个参考点还原真相。而 pyModeS 的PipeDecoder把这套繁琐的配对、状态管理、防错逻辑全部自动化——你只需要一行decode()拿到干净的latitude和longitude。【免费下载链接】pyModeSPython decoder for Mode S and ADS-B signals项目地址: https://gitcode.com/gh_mirrors/py/pyModeS创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考 为什么ADS-B不直接发经纬度？pyModeS CPR压缩位置解码原理全解
为什么ADS-B不直接发经纬度pyModeS CPR压缩位置解码原理全解【免费下载链接】pyModeSPython decoder for Mode S and ADS-B signals项目地址: https://gitcode.com/gh_mirrors/py/pyModeSpyModeS 是一款用 Python 编写的 ADS-B / Mode-S 航空信号解码器。如果你拆过一条 ADS-B 位置报文会发现一个反直觉的事实飞机并不直接广播经纬度而是发两个 17 位的CPR 压缩位置数字。这篇文章用最少代码讲透 ADS-B CPR 压缩位置的原理以及 pyModeS 如何把这两个神秘数字还原成地图上的坐标。 先看一条真实的 ADS-B 位置报文ADS-B 的空中位置消息Type Code 9-18 气压高度 / 20-22 GNSS 高度载荷只有56 位。拆开来看字段位宽作用TCType Code5 位消息类型SS NIC3 位监视状态、完整性等级高度码 AC12 位气压/卫星高度T F2 位时间同步位 CPR 格式位0偶帧 / 1奇帧CPR 纬度17 位压缩后的纬度CPR 经度17 位压缩后的经度一条来自法国上空的真实偶帧消息8D40058B58C901375147EFD09357它的 CPR 纬度是17161CPR 经度是97623。这不是 49.8°N、6.0°E——那是解码后的结果。为什么不直接发经纬度算一笔账就明白了两个 32 位浮点经纬度需要64 位已经超过整个 56 位载荷高度、呼号等字段就没地方放了即便压缩成 24 位定点经纬度全球精度也只有约 5 公里远不够识别跑道和航路ADS-B 工作在 1090 MHz 共享频段每架飞机约 0.5~1 秒广播一次每个字节都要用在刀刃上。所以 DO-260B 标准设计了 CPRCompact Position Reporting紧凑位置报告把位置拆成34 位量化数字 1 位格式位单帧精度达到米级但代价是——单条消息在全球范围内有几十上百个候选位置。️ CPR 的精髓把地球画成一张网格CPR 的核心思想是给地球表面画经纬网格偶帧even纬度切成 60 条带每带高 360°/60 6°奇帧odd纬度切成 59 条带每带高 360°/59 ≈6.10°经度更讲究在赤道上切成 59 份每份约 6.1°越靠近极地经度线越收敛份数也越少——高纬度地区经度圈本身变短格子要变宽才能保持近似等距。这个份数函数就叫NL(纬度)赤道是 59±87° 纬度降为 2。于是 17 位原始值就是格子编号 × 格宽 格内小数位置单帧的 17 位精度6° ÷ 131072 ≈ 5 米。精度足够但不知道落在全球第几个格子——偶帧有 60 个纬度候选 × 若干经度候选。这就是为什么 ADS-B 接收机必须记住飞机飞机持续广播偶帧 奇帧交替到达。偶帧落在 60 分格、奇帧落在 59 分格两套网格的错位恰好把候选位置唯一确定下来——这是 CPR 最巧妙的一笔。 pyModeS 如何还原经纬度三条解码路径pyModeS 把 CPR 数学全部收敛在 src/pyModeS/position/_cpr.py 中提供三种解法按手头信息从少到多排列1️⃣ 偶/奇帧配对解码最常用无需任何参考飞机每次广播都切换 even/odd 格式。收到一对相邻的帧后用公式求纬度j floor(59 × cpr_lat_even − 60 × cpr_lat_odd 0.5)再结合 NL(纬度) 解出经度。实现见 airborne_position_pair它还内置了两道保险两帧解出的纬度若落在不同 NL 纬度带说明配对失败或算出 |纬度| 90° 的物理不可能结果直接判为失败返回 None——真实案例中曾有一架从阿姆斯特丹起飞的飞机差点被解到白令海上空见 tests/test_cpr.py 中的回归测试。⏱️ 配对窗口默认10 秒间隔太久的两帧可能已跨越纬度带边界不再可信。2️⃣ 参考点解码已知道大致位置如果你已经有一个粗略参考点上次位置、雷达航迹、接收机位置单条消息就能解码偶帧参考点在 180 海里内即可。例如给上面那条法国消息一个粗参考 (49.0, 6.0)立刻得到49.8241°N, 6.0679°E验证向量见 tests/test_cpr.py。pyModeS 的PipeDecoder会在配对成功后自动把该位置升级为本地参考之后每一帧都能立即解出坐标不必等奇偶配对。3️⃣ 地面滑行 CPR需要机场参考地面移动报文BDS 0,6TC 5-8用的是更密的90° 纬度带bds06.py 负责提取原始字段经度还有四个 90° 象限候选因此必须给一个 45 海里内的参考点。pyModeS 贴心地内置了机场 ICAO 码数据库src/pyModeS/position/_airports.py一个代码就能换参考点from pyModeS import decode # 真实 DF18 地面滑行报文参考机场 LFBO图卢兹 r decode(903a23ff426a4e65f7487a775d17, surface_refLFBO) print(r[latitude], r[longitude]) # 43.6264, 1.3747 新手上手三步解码出第一条 ADS-B 位置pip install pyModeS3from pyModeS import decode # 偶帧 奇帧一起喂进去CPR 对自动配对 results decode([ 8D40058B58C901375147EFD09357, # even 8D40058B58C904A87F402D3B8C59, # odd ], timestamps[1.0, 2.0]) print(results[1][latitude], results[1][longitude]) # 49.8175 6.0844不想写代码pyModeS 附带modes命令行工具逗号分隔的多条消息共享同一个状态机CPR 对自动解析modes decode 8D40058B58C901375147EFD09357,8D40058B58C904A87F402D3B8C59完整状态机TTL 淘汰、幻影帧过滤、DF20/21 ICAO 校验的深入说明见 docs/pipe.md入门教程见 docs/quickstart.md。 一句话总结ADS-B 不发经纬度是因为56 位载荷装不下全球米级精度的完整坐标。CPR 用60/59 分格错位的数学技巧把位置压缩到 34 位再用飞机自己的连续广播偶帧奇帧或一个参考点还原真相。而 pyModeS 的PipeDecoder把这套繁琐的配对、状态管理、防错逻辑全部自动化——你只需要一行decode()拿到干净的latitude和longitude。【免费下载链接】pyModeSPython decoder for Mode S and ADS-B signals项目地址: https://gitcode.com/gh_mirrors/py/pyModeS创作声明：本文部分内容由AI辅助生成（AIGC），仅供参考