C++ proj库使用说明与字符串说明

proj使用说明

PROJ github:
https://github.com/OSGeo/PROJ

1. 引用proj头文件

#include "proj_api.h"

2. 构建坐标系定义指针

构建所需的proj_string字符串可参见： https://epsg.io/
字符串搜索用例：
首先在首页搜索栏输入所需坐标系名称，点击search搜索。

之后选择所需的坐标系统。

在页面下方export中选择PROJ4可以在右边获取到proj_string字符串。

字符串各个字段含义可参见：https://blog.csdn.net/u013240519/article/details/113243000

projPJ pWGS84 = nullptr;  // WGS84坐标系指针
projPJ pGKP = nullptr;    //高斯克鲁格坐标系指针
const char* wgs84_string = "+proj=longlat +ellps=WGS84 +datum=WGS84";
pWGS84 = pj_init_plus(wgs84_string);
// 以40度带中央经度120为例
const char* gk_string = "+proj=temrc +ellps=WGS84 +datum=WGS84 +x_0=500000 +lon_0=120 +units=m";
pGKP = pj_init_plus(gk_string);

3. 坐标转换

// GKP to WGS84
bool WGS842GKP(const double dLon/*rad*/, const double dLat/*rad*/, double& dX, double &dY)
{
    dX = dLon;
    dY = dLat;

    int iRev = pj_transform(pWGS84, pGKP, 1, 1, &dX, &dY, nullptr);
    if(iRev != 0)
    {
        return false;
    }
    return true;
}

4. 转换结束或更换proj_string前释放内存

if(pWGS84 != nullptr)
{
    pj_free(pWGS84);
    pWGS84 = nullptr;
}
if(pGKP != nullptr)
{
    pj_free(pGKP);
    pGKP = nullptr;
}