UTM坐标与wgs84坐标转换关系

类定义 需要安装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;
}