基于GNSS经纬度数据的速度估算不确定性分析方法咨询

核心思路：误差传播法 #

利用位置误差的各向同性假设，通过误差传播推导速度的置信区间，避免直接加减距离误差导致的低速负值问题。

步骤1：坐标预处理（推荐） #

将经纬度转换为平面直角坐标（如UTM投影），短距离下平面距离与球面距离误差可忽略，且便于误差计算。若保留经纬度，直接用球面距离公式计算Δs即可，误差传播逻辑一致。

步骤2：基础速度计算 #

对相邻采样点(t₁, lat₁, lon₁)和(t₂, lat₂, lon₂)：

• 计算时间差：Δt = t₂ - t₁（假设时间无误差）
• 计算两点距离：Δs = √[(x₂-x₁)² + (y₂-y₁)²]（平面坐标）或球面距离
• 平均速度估算：v = Δs / Δt

步骤3：推导速度的标准差 #

已知每个位置点的水平位置标准差σ = 5m，且x、y方向误差独立、各向同性：

• 两点距离的误差来自两个位置点的独立误差，因此距离的标准差为：
  σ_Δs = σ * √2 ≈ 7.07m
• 由于时间无误差，速度的标准差直接由距离误差除以时间差得到：
  σ_v = σ_Δs / Δt = (5*√2)/Δt

步骤4：生成95%置信区间 #

95%置信度对应2σ范围，因此每个速度点的置信区间为：
[v - 2*σ_v, v + 2*σ_v]

步骤5：低速场景处理 #

当估算速度v接近0时，置信区间下限可能出现负值——这是GNSS位置误差的合理体现，说明低速时无法准确判断运动方向或是否静止。若需美化图表，可将下限截断为0，但必须在图注中说明这是近似处理，实际置信区间包含负值。

步骤6：绘制速度-时间图 #

• 以采样时间为横轴，估算速度为纵轴绘制主曲线
• 绘制上下两条置信边界曲线（或填充区间），标注为「95%置信区间（2σ）」

内容的提问来源于stack exchange，提问作者JJS