UTM坐标与wgs84坐标转换关系

最新推荐文章于 2024-12-12 17:20:27 发布

雪天枫

最新推荐文章于 2024-12-12 17:20:27 发布

阅读量2.9k

点赞数 2

CC 4.0 BY-SA版权

分类专栏： C++ 文章标签：经纬度与UTM左边转换

本文链接：https://blog.youkuaiyun.com/zhaodeming000/article/details/112004638

C++ 专栏收录该内容

3 篇文章

订阅专栏

该博客介绍了如何使用Proj库在C++中实现经纬度到UTM坐标的转换，以及UTM坐标回转为经纬度的过程。通过创建`Projection`类，定义了初始化球体和投影的方法，以及转换函数。博客详细展示了类的定义、构造函数、析构函数以及转换方法的实现，便于理解并应用到地理坐标系统的转换中。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

类定义需要安装Proj库

#ifndef PROJECTION_H
#define PROJECTION_H

// #include <projects.h>
#define ACCEPT_USE_OF_DEPRECATED_PROJ_API_H
#include <proj_api.h>

class Projection
{
private:
    projPJ m_pj_wgs84;  //球体:wgs84、投影:经纬度投影
    projPJ m_pj_utm;    //球体:wgs84、投影UTM
    int m_nUTMZone;     //utm带号
    bool m_bNorth;      //是否为北半球

public:
    Projection();
    ~Projection();
    bool wgs84ToUTM(const double dLongi, const double dLati, double &dX, double &dY);
    bool utmToWgs84(const double dX, const double dY, double &dLongi, double &dLati);

};

#endif // PROJECTION_H

#include <iostream>
#include <string>
using namespace std;

Projection::Projection():
    m_nUTMZone(-10000), m_bNorth(false), m_pj_utm(NULL)
{
    const char *wgs84 = "+proj=longlat +datum=WGS84 +no_defs ";//GPS所用坐标系,EPSG:4326
    m_pj_wgs84 = pj_init_plus(wgs84);
}
Projection::~Projection()
{
    if(m_pj_utm != NULL)
    {
        pj_free(m_pj_utm);
        m_pj_utm = NULL;
    }
    if(m_pj_wgs84 != NULL)
    {
        pj_free(m_pj_wgs84);
        m_pj_wgs84 = NULL;
    }
}
/**
* function：经纬度转为UTM坐标
* input：经度，维度，x轴，y轴
*/
bool Projection::wgs84ToUTM(const double dLongi, const double dLati, double &dX, double &dY)
{
  
    int nUTMZone = (int)((dLongi+186.0)/6.0);
    bool bNorth = dLati > 0 ? true : false;    
    if(m_nUTMZone != nUTMZone || m_bNorth != bNorth)
    {
        //"+proj=utm +zone=50 +datum=WGS84 +units=m +no_defs";
        // "+proj=utm +zone=48 +south +datum=WGS84 +units=m +no_defs"
        std::string qstrUTM = "+proj=utm +zone=" + std::to_string(nUTMZone);
        if(!bNorth)
        {
            qstrUTM += " +south ";
        }
        qstrUTM += " +datum=WGS84 +units=m +no_defs";
        if(m_pj_utm != NULL)
        {
            pj_free(m_pj_utm);
            m_pj_utm = NULL;
        }
        string strUTM = qstrUTM;
        const char *pUTM = strUTM.c_str();

        m_pj_utm = pj_init_plus(pUTM);
        

        m_nUTMZone = nUTMZone;
        m_bNorth = bNorth;
    }
    dX = dLongi * DEG_TO_RAD;
    dY = dLati * DEG_TO_RAD;
    pj_transform(m_pj_wgs84, m_pj_utm, 1, 1, &dX, &dY, NULL);

    return true;
}
/**
* function：UTM坐标转为经纬度
* input：x轴，y轴，经度，维度
*/
bool Projection::utmToWgs84(const double dX, const double dY, double &dLongi, double &dLati)
{
    if(m_pj_utm == NULL)
    {
        return false;
    }
    dLongi = dX;
    dLati = dY;
    pj_transform(m_pj_utm, m_pj_wgs84, 1, 1, &dLongi, &dLati, NULL);
    dLongi = dLongi * RAD_TO_DEG;
    dLati = dLati * RAD_TO_DEG;

    return  true;
}