GPS转平面视图坐标的变形问题解决

最新推荐文章于 2023-06-21 11:03:15 发布

原创最新推荐文章于 2023-06-21 11:03:15 发布 · 948 阅读

0 ·

CC 4.0 BY-SA版权

Qt 专栏收录该内容

6 篇文章

订阅专栏

本文讨论了如何将基于经纬度坐标的地图投影到OpenGL视图坐标平面上，并解决了投影过程中产生的变形问题，通过计算正确的比例系数来保持地图显示的准确性。

1. 经纬度的范围

经度0°——180°（东行,标注E）0°——180°（西行,标注W）
纬度0°——90°N、0°——90°S

2. 需求

将基于经纬度坐标的地图绘制到基于opengl视图坐标的平面上

3. 问题

由于地球是个椭球体，每单位的经度和纬度对应的地面距离是不相等的：

纬度变化一度，球面南北方向距离变化：πR/180 ........111.7(km)
经度变化一度，球面东西方向距离变化：πR/180*cos(latitude) ....111.7*cos(latitude) (km)

比如北京的纬度 latitude-B = 40、cos(latitude-B) = 0.766，经度变化1度，则东西方向距离变化 111.7 * cos(40) = 85.567(km)

4. 解决变形问题

如果投影到视图平面的地图长宽比例直接使用经度跨度与纬度跨度之比，投影后的地图会在水平方向拉长，发生变形。

为了保持正确的显示，可以通过以下方法计算平面视图的比例：

viewWidth / viewHeight = longitudeSpan * cos(latidude) / latitudeSpan

如：

double CCoordinateManager::lengthToWidthRatio()
{
    if (CShapeMath::isZero(m_mapRect.width()) || CShapeMath::isZero(m_mapRect.height()))
    {
        return ((double)(m_screenRect.width()) / (double)(m_screenRect.height()));
    }
    else
    {
        return m_mapRect.width() / m_mapRect.height() * qCos(qDegreesToRadians(m_mapRect.y()));
    }
}