如何计算指定位置的GPS（GNSS）定位不确定性？

我正在开发一款用于计算指定位置GNSS（GPS）定位不确定性的软件。多数GNSS接收机会分别上报经度不确定性（Longitude Uncertainty）与纬度不确定性（Latitude Uncertainty）。

我原本假设这可以转化为椭圆不确定性模型，公式如下：

r = LongUncert * cos(bearing) + LatUncert * sin(bearing)

注：此处方位角（bearing）从0°（东）开始逆时针递增，例如90°代表北，180°代表西。

但有人提出，这种上报的不确定性更适合用矩形不确定性模型描述——以纬度不确定性为高度、经度不确定性为宽度，我也认为这种模型具备合理性。

我的核心应用需求是：告知用户当前位置与目标经纬度坐标的距离，并提供以米为单位的±不确定性数值。

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