C++实现百度坐标(BD09)及GCJ02与WGS84之间的转换

更新时间：2023年03月03日 09:55:28 作者：毓树麟风

这篇文章主要为大家详细介绍了C++实现百度坐标(BD09)及GCJ02与WGS84之间的转换的方法，文中的示例代码讲解详细，希望对大家有所帮助

实现代码

mainwindow.h

#ifndef MAINWINDOW_H
#define MAINWINDOW_H

#include <QMainWindow>
#include <math.h>
typedef struct _POSITION
{
    double longitude;
    double latitude;
}POSITION;

namespace Ui {
class MainWindow;
}

class MainWindow : public QMainWindow
{
    Q_OBJECT

public:
    explicit MainWindow(QWidget *parent = 0);
    ~MainWindow();
private slots:
    double translate_lon(double lon,double lat);
    double translate_lat(double lon,double lat);
    bool outof_China(double lon,double lat);
    POSITION bd09togcj02(double bd_lon, double bd_lat);
    POSITION gcj02tobd09(double gcj_lon,double gcj_lat);
    POSITION gcj02towgs84(double gcj_lon,double gcj_lat);
    POSITION wgs84togcj02(double wgs_lon,double wgs_lat);
    void pushbutton1();
    void pushbutton2();
    void pushbutton3();
private:
    Ui::MainWindow *ui;
    POSITION bd_pos;
    POSITION gcj_pos;
    POSITION wgs_pos;
    double x_PI = 3.14159265358979323846 * 3000.0 / 180.0;
    double PI = 3.1415926535897932384626;
    double a = 6378245.0;
    double ee = 0.00669342162296594323;
};

#endif // MAINWINDOW_H

mainwindow.cpp

#include "mainwindow.h"
#include "ui_mainwindow.h"

MainWindow::MainWindow(QWidget *parent) :
    QMainWindow(parent),
    ui(new Ui::MainWindow)
{
    ui->setupUi(this);
    setWindowTitle("Position Translate");
    connect(ui->pushButton  ,SIGNAL(clicked()),this,SLOT(pushbutton1()));
    connect(ui->pushButton_2,SIGNAL(clicked()),this,SLOT(pushbutton2()));
    connect(ui->pushButton_3,SIGNAL(clicked()),this,SLOT(pushbutton3()));
}
MainWindow::~MainWindow()
{
    delete ui;
}
bool MainWindow::outof_China(double lon, double lat)
{
    return(lon<72.004 || lon>137.8374 || lat<0.8293 || lat >55.8271 || false);
}
POSITION MainWindow::bd09togcj02(double bd_lon, double bd_lat)
{
    double x = bd_lon - 0.0065;
    double y = bd_lat - 0.006;
    double z = sqrt(x*x + y*y) - 0.00002*sin(y*x_PI);
    double theta = atan2(y,x) - 0.000003*cos(x*x_PI);
    gcj_pos.longitude = z*cos(theta);
    gcj_pos.latitude = z*sin(theta);
    return gcj_pos;
}
POSITION MainWindow::gcj02tobd09(double gcj_lon, double gcj_lat)
{
     double z = sqrt(gcj_lon*gcj_lon + gcj_lat*gcj_lat) + 0.00002*sin(gcj_lat * x_PI);
     double theta = atan2(gcj_lat,gcj_lon) + 0.000003 * cos(gcj_lon * x_PI);
     bd_pos.longitude = z*cos(theta) + 0.0065;
     bd_pos.latitude = z*sin(theta) + 0.006;
     return bd_pos;
}
double MainWindow::translate_lon(double lon, double lat)
{
    double ret = 300.0 + lon +2.0*lat + 0.1*lon*lon +0.1*lon*lat + 0.1*sqrt(abs(lon));
    ret += (20.0 * sin(6.0*lon*PI) + 20.0*sin(2.0*lon*PI)) *2.0 / 3.0;
    ret += (20.0 * sin(lon*PI) + 40.0*sin(lon/3.0 *PI)) *2.0 /3.0;
    ret += (150 * sin(lon/12.0 *PI) + 300.0*sin(lon/30.0 * PI)) *2.0 /3.0;
    return ret;
}
double MainWindow::translate_lat(double lon, double lat)
{
    double ret = -100 + 2.0*lon + 3.0*lat + 0.2*lat*lat + 0.1*lon*lat + 0.2*sqrt((abs(lon)));
    ret += (20.0 *sin(6.0*lon*PI) + 20*sin(2.0*lon*PI)) *2.0 /3.0;
    ret += (20.0 *sin(lat*PI) + 40.0*sin(lat/3.0*PI)) *2.0 /3.0;
    ret += (160.0*sin(lat/12.0*PI) + 320.0*sin(lat/30.0 *PI)) *2.0 /3.0;
    return ret;
}
POSITION MainWindow::gcj02towgs84(double gcj_lon, double gcj_lat)
{
    if(outof_China(gcj_lon,gcj_lat))
    {
        wgs_pos.longitude = gcj_lon;
        wgs_pos.latitude = gcj_lat;
        return wgs_pos;
    }
    else
    {
        double dlat = translate_lat(gcj_lon - 105.0,gcj_lat -35.0);
        double dlon = translate_lon(gcj_lon - 105.0,gcj_lat -35.0);
        double radlat = gcj_lat/180.0 *PI;
        double magic = sin(radlat);
        magic = 1 - ee*magic*magic;
        double squrtmagic = sqrt(magic);
        dlon = (dlon *180.0)/(a/squrtmagic*cos(radlat)*PI);
        dlat = (dlat *180.0)/((a*(1-ee))/(magic * squrtmagic)*PI);
        wgs_pos.longitude = gcj_lon - dlon;
        wgs_pos.latitude = gcj_lat - dlat;
        return wgs_pos;
    }
}
POSITION MainWindow::wgs84togcj02(double wgs_lon, double wgs_lat)
{
    if(outof_China(wgs_lon,wgs_lat))
    {
        gcj_pos.longitude = wgs_lon;
        gcj_pos.latitude = wgs_lat;
        return gcj_pos;
    }
    else
    {
        double dlat = translate_lat(wgs_lon - 105.0,wgs_lat - 35.0);
        double dlon = translate_lon(wgs_lon - 105.0,wgs_lat - 35.0);
        double radlat = wgs_lat/180.0 * PI;
        double magic = sin(radlat);
        magic = 1 - ee*magic*magic;
        double squrtmagic = sqrt(magic);
        dlon = (dlon *180.0)/(a/squrtmagic*cos(radlat)*PI);
        dlat = (dlat *180.0)/((a*(1-ee))/(magic * squrtmagic)*PI);
        gcj_pos.longitude = wgs_lon + dlon;
        gcj_pos.latitude = wgs_lat +dlat;
        return gcj_pos;
    }
}
/*************************************************************************
 *
 * 读取 BD09 坐标，转换为 WGS84 和 GCJ02 坐标
 *
 ***********************************************************************/
void MainWindow::pushbutton1()
{
    double bd09_lon = ui->lineEdit_bd09lon->text().toDouble();
    double bd09_lat = ui->lineEdit_bd09lat->text().toDouble();
    double gcj02_lon = bd09togcj02(bd09_lon,bd09_lat).longitude;
    double gcj02_lat = bd09togcj02(bd09_lon,bd09_lat).latitude;
    double wgs84_lon = gcj02towgs84(gcj02_lon,gcj02_lat).longitude;
    double wgs84_lat =gcj02towgs84(gcj02_lon,gcj02_lat).latitude;
    ui->lineEdit_gcj02lon->setText(QString::number(gcj02_lon,'d',9));
    ui->lineEdit_gcj02lat->setText(QString::number(gcj02_lat,'d',9));
    ui->lineEdit_wgs84lon->setText(QString::number(wgs84_lon,'d',9));
    ui->lineEdit_wgs84lat->setText(QString::number(wgs84_lat,'d',9));
}
/*************************************************************************
 *
 * 读取 GCJ02 坐标，转换为 WGS84 和 BD09 坐标
 *
 ***********************************************************************/
void MainWindow::pushbutton2()
{
    double gcj02_lon = ui->lineEdit_gcj02lon->text().toDouble();
    double gcj02_lat = ui->lineEdit_gcj02lat->text().toDouble();
    double bd09_lon = gcj02tobd09(gcj02_lon,gcj02_lat).longitude;
    double bd09_lat = gcj02tobd09(gcj02_lon,gcj02_lat).latitude;
    double wgs84_lon = gcj02towgs84(gcj02_lon,gcj02_lat).longitude;
    double wgs84_lat = gcj02towgs84(gcj02_lon,gcj02_lat).latitude;
    ui->lineEdit_bd09lon->setText(QString::number(bd09_lon,'d',9));
    ui->lineEdit_bd09lat->setText(QString::number(bd09_lat,'d',9));
    ui->lineEdit_wgs84lon->setText(QString::number(wgs84_lon,'d',9));
    ui->lineEdit_wgs84lat->setText(QString::number(wgs84_lat,'d',9));
}
/*************************************************************************
 *
 * 读取 WGS84 坐标，转换为 GCJ02 和 BD09 坐标
 *
 ***********************************************************************/
void MainWindow::pushbutton3()
{
    double wgs84_lon = ui->lineEdit_wgs84lon->text().toDouble();
    double wgs84_lat = ui->lineEdit_wgs84lat->text().toDouble();
    double gcj02_lon = wgs84togcj02(wgs84_lon,wgs84_lat).longitude;
    double gcj02_lat = wgs84togcj02(wgs84_lon,wgs84_lat).latitude;
    double bd09_lon = gcj02tobd09(gcj02_lon,gcj02_lat).longitude;
    double bd09_lat = gcj02tobd09(gcj02_lon,gcj02_lat).latitude;
    ui->lineEdit_bd09lon->setText(QString::number(bd09_lon,'d',9));
    ui->lineEdit_bd09lat->setText(QString::number(bd09_lat,'d',9));
    ui->lineEdit_gcj02lon->setText(QString::number(gcj02_lon,'d',9));
    ui->lineEdit_gcj02lat->setText(QString::number(gcj02_lat,'d',9));
}

算法来源于网上搜集，运行结果如下：

到此这篇关于C++实现百度坐标(BD09)及GCJ02与WGS84之间的转换的文章就介绍到这了,更多相关C++坐标转换内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！

您可能感兴趣的文章:

枚举和宏的区别详细解析
枚举常量具有类型，但宏没有类型，枚举变量具有与普通变量相同的诸如作用域、值等性质，但宏没有，宏不是语言的一部分，它是一种预处理替换符
2013-09-09
浅谈C#中List<T>对象的深度拷贝问题
下面小编就为大家带来一篇浅谈C#中List<T>对象的深度拷贝问题。小编觉得挺不错的，现在就分享给大家，也给大家做个参考。一起跟随小编过来看看吧
2017-01-01
Qt TCP网络通信学习
用于数据传输的低层网络协议，多个物联网协议都是基于TCP协议的，这篇文章为大家介绍了Qt TCP网络通信，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
2022-08-08
C语言实现贪吃蛇游戏代码
大家好，本篇文章主要讲的是C语言实现贪吃蛇游戏代码，感兴趣的同学赶快来看一看吧，对你有帮助的话记得收藏一下
2022-02-02
C++ LeetCode543题解二叉树直径
这篇文章主要为大家介绍了C++ LeetCode543题解二叉树直径，有需要的朋友可以借鉴参考下，希望能够有所帮助，祝大家多多进步，早日升职加薪
2022-12-12
C++如何获取鼠标点击位置
这篇文章主要介绍了C++如何获取鼠标点击位置问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教
2023-08-08
C++使用read()和write()读写二进制文件
以文本形式读写文件和以二进制形式读写文件的区别,并掌握了用重载的 >> 和 << 运算符实现以文本形式读写文件,在此基础上,本节将讲解如何以二进制形式读写文件
2023-10-10
C语言中的字符型数据与ASCII码表
这篇文章主要介绍了C语言中的字符型数据与ASCII码表，具有很好的参考价值，希望对大家有所帮助。如有错误或未考虑完全的地方，望不吝赐教
2023-01-01
cin.get()和cin.getline()之间的区别
以下是对cin.get()和cin.getline()的区别进行了详细的分析介绍，需要的朋友可以过来参考下，希望对大家有所帮助
2013-09-09
浅析C语言字中的符串格式化显示
以下是对C语言字中的符串格式化显示进行了详细的分析介绍，需要的朋友可以过来参考下
2013-08-08

C++实现百度坐标(BD09)及GCJ02与WGS84之间的转换

相关文章

最新评论

大家感兴趣的内容

最近更新的内容

常用在线小工具