广告
返回顶部
首页 > 资讯 > 后端开发 > 其他教程 >C++怎么实现百度坐标及GCJ02与WGS84之间的转换
  • 661
分享到

C++怎么实现百度坐标及GCJ02与WGS84之间的转换

2023-07-05 08:07:12 661人浏览 独家记忆
摘要

这篇文章主要讲解了“c++怎么实现百度坐标及GCJ02与WGS84之间的转换”,文中的讲解内容简单清晰,易于学习与理解,下面请大家跟着小编的思路慢慢深入,一起来研究和学习“C++怎么实现百度坐标及GCJ02与WGS84之间的转换”吧!实现代

这篇文章主要讲解了“c++怎么实现百度坐标及GCJ02与WGS84之间的转换”,文中的讲解内容简单清晰,易于学习与理解,下面请大家跟着小编的思路慢慢深入,一起来研究和学习“C++怎么实现百度坐标及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_OBJECTpublic:    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;    }}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));}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));}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++怎么实现百度坐标及GCJ02与WGS84之间的转换

感谢各位的阅读,以上就是“C++怎么实现百度坐标及GCJ02与WGS84之间的转换”的内容了,经过本文的学习后,相信大家对C++怎么实现百度坐标及GCJ02与WGS84之间的转换这一问题有了更深刻的体会,具体使用情况还需要大家实践验证。这里是编程网,小编将为大家推送更多相关知识点的文章,欢迎关注!

--结束END--

本文标题: C++怎么实现百度坐标及GCJ02与WGS84之间的转换

本文链接: https://www.lsjlt.com/news/350624.html(转载时请注明来源链接)

有问题或投稿请发送至: 邮箱/279061341@qq.com    QQ/279061341

本篇文章演示代码以及资料文档资料下载

下载Word文档到电脑,方便收藏和打印~

下载Word文档
猜你喜欢
软考高级职称资格查询
编程网,编程工程师的家园,是目前国内优秀的开源技术社区之一,形成了由开源软件库、代码分享、资讯、协作翻译、讨论区和博客等几大频道内容,为IT开发者提供了一个发现、使用、并交流开源技术的平台。
  • 官方手机版

  • 微信公众号

  • 商务合作