WGS84经纬度坐标转化成UTM坐标

最新推荐文章于 2023-04-20 15:26:03 发布
最新推荐文章于 2023-04-20 15:26:03 发布
阅读量4.9k 收藏 26
点赞数 6
分类专栏： 无人驾驶 文章标签： 其他
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_29923461/article/details/119949010
版权
本文介绍了无人驾驶和机器人领域中常见的WGS84经纬度坐标与UTM坐标的转换，提供了使用geodesy功能包和GeographicLib库的详细实现。通过设置原点坐标，实现了从经纬度到UTM的转换，并给出了测试案例展示转换过程。
摘要由CSDN通过智能技术生成

WGS84经纬度坐标转化成UTM坐标


  在无人驾驶或者机器人领域，经常会涉及到WGS 84 经纬度坐标与UTM坐标之间的转换，本文总结了两种方式来实现两者的转化。

1 利用geodesy功能包实现转化

  头文件：

#ifndef WGS2UTM_H
#define WGS2UTM_H

#include <iostream>
#include <fstream>
#include <vector>
#include <Eigen/Geometry>
#include <GeographicLib/UTMUPS.hpp>
#include <geographic_msgs/GeoPointStamped.h>
#include <geodesy/utm.h>
#include <ros/ros.h>

using namespace std;

class WGS2UTM
{
    public:

    /**
     * @brief 构造函数
     */
    WGS2UTM();

    /**
     * @brief 获取选取原点的经纬度lla坐标
     * @return  所选取原点的经纬度lla坐标
     */
    inline Eigen::Vector3d get_origin_lla() 
    {
        return origin_LLA_;
    }

    /**
     * @brief 获取选取原点的经纬度utm坐标
     * @return  所选取原点的经纬度utm坐标
     */
    inline Eigen::Vector3d get_origin_tum() 
    {
        return origin_UTM_;
    }

    /***********************************************************/
    /***********************************************************/
    /*****利用geodesy功能包实现WGS84经纬度坐标转化成UTM坐标*****/
    /***********************************************************/
    /***********************************************************/

    /**
     * @brief 把选取的原点的经纬度lla坐标转化为东北天utm坐标
     * @param origin_LLA    输入--原点经纬度lla坐标
     */
    void set_origin_1(Eigen::Vector3d origin_LLA);

    /**
     * @brief 把经纬度lla坐标转化为东北天utm坐标
     * @param LLA    输入--lla坐标
     * @param UTM    输出--utm坐标
     */
    void LLA2UTM_1(const Eigen::Vector3d& LLA, Eigen::Vector3d& UTM);

    /**
     * @brief 把经纬度lla坐标转化为相对于所选取原点的东北天utm坐标
     * @param LLA    输入--lla坐标
     * @param UTM    输出--utm坐标
     */
    void LLA2UTM_origin_1(const Eigen::Vector3d& LLA, Eigen::Vector3d& UTM_origin);

private:

    Eigen::Vector3d origin_LLA_;
    Eigen::Vector3d origin_UTM_;

    double torad = M_PI / 180;
};
#endif //WGS2UTM_H

  源文件：

#include "convert_wgs2utm/wgs2utm.h"

WGS2UTM::WGS2UTM() {}

void WGS2UTM::set_origin_1(Eigen::Vector3d origin_LLA) 
{
    origin_LLA_ = origin_LLA;
    LLA2UTM_1(origin_LLA_, origin_UTM_);
}

void WGS2UTM::LLA2UTM_1(const Eigen::Vector3d& LLA, Eigen::Vector3d& UTM)
{
    geographic_msgs::GeoPointStampedPtr gps_msg(new geographic_msgs::GeoPointStamped());
    gps_msg->position.latitude  = LLA[0];
    gps_msg->position.longitude = LLA[1];
    gps_msg->position.altitude  = LLA[2];
    geodesy::UTMPoint utm_point;
    geodesy::fromMsg(gps_msg->position, utm_point);
    UTM[0] = utm_point.easting;
    UTM[1] = utm_point.northing;
    UTM[2] = utm_point.altitude;
}

void WGS2UTM::LLA2UTM_origin_1(const Eigen::Vector3d& LLA, Eigen::Vector3d& UTM_origin)
{
    //经纬度转换为UTM投影
    Eigen::Vector3d UTM;
    LLA2UTM_1(LLA, UTM);

    //UTM2ENU, UTM投影平移至东北天坐标系, 东北天坐标系原点为第一个采样点
    UTM_origin = UTM - origin_UTM_;
}

  测试案例：

#include <iostream>
#include "convert_wgs2utm/wgs2utm.h"

using namespace std;

int main()
{
    WGS2UTM wgs2utm;

    Eigen::Vector3d origin_LLA;
    Eigen::Vector3d origin_UTM;
    origin_LLA[0] = 30.9543508;
    origin_LLA[1] = 121.4795984;
    origin_LLA[2] = 17.0970000;

    wgs2utm.set_origin_1(origin_LLA);
    origin_UTM = wgs2utm.get_origin_tum();
    cout << "origin_UTM: (" <<  origin_UTM[0] << ", " << origin_UTM[1] << ", " << origin_UTM[2] << ")" << endl;
    cout << "******************************" << endl;

    Eigen::Vector3d P_LLA;  // 任意一点的经纬度坐标
    Eigen::Vector3d P_UTM;  // 该点相对于所选取原点的utm坐标

    P_LLA[0] = 30.9543508;
    P_LLA[1] = 121.4755984;
    P_LLA[2] = 17.0970000;
    wgs2utm.LLA2UTM_origin_1(P_LLA, P_UTM);
    cout << "P_UTM: (" <<  P_UTM[0] << ", " << P_UTM[1] << ", " << P_UTM[2] << ")" << endl;
    cout << "******************************" << endl;

}

  输出结果如下：
在这里插入图片描述

2 利用GeographicLib功能包实现转化

  头文件：

#ifndef WGS2UTM_H
#define WGS2UTM_H

#include <iostream>
#include <fstream>
#include <vector>
#include <Eigen/Geometry>
#include <GeographicLib/UTMUPS.hpp>
#include <geographic_msgs/GeoPointStamped.h>
#include <geodesy/utm.h>
#include <ros/ros.h>

using namespace std;

class WGS2UTM
{
    public:

    /**
     * @brief 构造函数
     */
    WGS2UTM();

    /**
     * @brief 获取选取原点的经纬度lla坐标
     * @return  所选取原点的经纬度lla坐标
     */
    inline Eigen::Vector3d get_origin_lla() 
    {
        return origin_LLA_;
    }

    /**
     * @brief 获取选取原点的经纬度utm坐标
     * @return  所选取原点的经纬度utm坐标
     */
    inline Eigen::Vector3d get_origin_tum() 
    {
        return origin_UTM_;
    }

    /***********************************************************/
    /***********************************************************/
    /**利用GeographicLib功能包实现WGS84经纬度坐标转化成UTM坐标**/
    /***********************************************************/
    /***********************************************************/
    
    /**
     * @brief 把选取的原点的经纬度lla坐标转化为东北天utm坐标
     * @param origin_LLA    输入--原点经纬度lla坐标
     */
    void set_origin_2(Eigen::Vector3d origin_LLA);

    /**
     * @brief 把经纬度lla坐标转化为东北天utm坐标
     * @param LLA    输入--lla坐标
     * @param UTM    输出--utm坐标
     */
    void LLA2UTM_2(const Eigen::Vector3d& LLA, Eigen::Vector3d& UTM);

    /**
     * @brief 把经纬度lla坐标转化为相对于所选取原点的东北天utm坐标
     * @param LLA    输入--lla坐标
     * @param UTM    输出--utm坐标
     */
    void LLA2UTM_origin_2(const Eigen::Vector3d& LLA, Eigen::Vector3d& UTM_origin);

private:

    Eigen::Vector3d origin_LLA_;
    Eigen::Vector3d origin_UTM_;

    double torad = M_PI / 180;
};
#endif //WGS2UTM_H

  源文件：

void WGS2UTM::set_origin_2(Eigen::Vector3d origin_LLA) 
{
    origin_LLA_ = origin_LLA;
    LLA2UTM_2(origin_LLA_, origin_UTM_);
}

void WGS2UTM::LLA2UTM_2(const Eigen::Vector3d& LLA, Eigen::Vector3d& UTM)
{
    GeographicLib::Math::real lat = LLA[0];
    GeographicLib::Math::real lon = LLA[1];
    GeographicLib::Math::real alt = LLA[2];
    int UTMZone;
    bool isNorth;   
    GeographicLib::UTMUPS::Forward(lat, lon, UTMZone, isNorth, UTM[0], UTM[1]);
    UTM[2] = LLA[2];
}

void WGS2UTM::LLA2UTM_origin_2(const Eigen::Vector3d& LLA, Eigen::Vector3d& UTM_origin)
{
    //经纬度转换为UTM投影
    Eigen::Vector3d UTM;
    LLA2UTM_2(LLA, UTM);

    //UTM2ENU, UTM投影平移至东北天坐标系, 东北天坐标系原点为第一个采样点
    UTM_origin = UTM - origin_UTM_;
}

  测试案例：

#include <iostream>
#include "convert_wgs2utm/wgs2utm.h"

using namespace std;

int main()
{
    WGS2UTM wgs2utm;

    Eigen::Vector3d origin_LLA;
    Eigen::Vector3d origin_UTM;
    origin_LLA[0] = 30.9543508;
    origin_LLA[1] = 121.4795984;
    origin_LLA[2] = 17.0970000;

    wgs2utm.set_origin_2(origin_LLA);
    origin_UTM = wgs2utm.get_origin_tum();
    cout << "origin_UTM: (" <<  origin_UTM[0] << ", " << origin_UTM[1] << ", " << origin_UTM[2] << ")" << endl;
    cout << "******************************" << endl;

    Eigen::Vector3d P_LLA;  // 任意一点的经纬度坐标
    Eigen::Vector3d P_UTM;  // 该点相对于所选取原点的utm坐标

    P_LLA[0] = 30.9543508;
    P_LLA[1] = 121.4755984;
    P_LLA[2] = 17.0970000;
    wgs2utm.LLA2UTM_origin_2(P_LLA, P_UTM);
    cout << "P_UTM: (" <<  P_UTM[0] << ", " << P_UTM[1] << ", " << P_UTM[2] << ")" << endl;
    cout << "******************************" << endl;

}
Roar冷颜
关注
  • 6
    点赞
  • 踩
  • 26
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
【GlobalMapper精品教程】080：WGS84转UTM投影
「 刘一哥与GIS的故事」
05-12 325
查看源坐标系：双击图层的，图层投影选项卡，数据的已有坐标系为WGS84地理坐标系。选择输出投影为UTM，基准为WGS84，中央子午线为111，点击确定。选择Shapefile格式，点击确定。打开配套案例数据包中的。查看投影后的坐标系统。
Python库 | utm-0.3.0.tar.gz
05-24
资源分类：Python库 所属语言：Python 资源全名：utm-0.3.0.tar.gz 资源来源：官方 安装方法：https://lanzao.blog.csdn.net/article/details/101784059
c# WGS84坐标转换成UTM坐标
buaazx1992的博客
07-03 2472
通过查阅相关资料，有一篇博客里面给出了WGS84坐标转换成UTM坐标的功能，该博主使用matlab程序实现坐标的转换，基于已有的公式，用c#语言实现坐标的转换。 参考的博客连接如下： https://blog.csdn.net/tap880507/article/details/88172120. c#实现坐标转换的代码如下： public static double a = 637...
GPS经纬度坐标转UTM坐标(c++)
Bobsweetie的专栏
06-18 7110
UTM的简单介绍 我们在地理课上学过用经纬度表示在地球上的坐标，这种表示方法适合在一个三维的地球曲面上 描述一个坐标，但是如果我们想在二维的平面上描述一个地理坐标应该怎么做呢？ 首先是将地球曲面投影到平面上 如下图在地球中心为投影中心，向圆柱形进行射线的投影，然后将圆柱形的侧面进行展开，就可以得到一个平面投影的坐标。 如下图所示就是一幅世界的平面地图，正如我们平时看到的地图 2.将平面分成许多个cell，叫做zone UTM的东北坐标只有两个数字，还需要zone的区间号才能转换为经纬度坐标 因此准确
经纬度转换为UTM坐标
weixin_45626706的博客
04-20 4685
UTM（通用横向墨卡托投影）： 是一种以米为单位的坐标系统，用于地图和GPS导航。将地球划分为60个纵向的区域。每个区域宽6度，从赤道开始往南北两级方向划分。每个区域都有一个特定的字母。
WGS84转UTM
01-12
WGS84Z转投影到UTM，内含示例数据，代码中有详细的注释说明，输入数据之后自动计算带号。非常实用。
PROJ将WGS84转UTM坐标
泠山的博客
03-18 2521
PROJ是一个通用的坐标转换软件，它将地理空间坐标从一个坐标参考系统(coordinate reference system, CRS)转换到另一个坐标参考系统。这包括地图投影以及大地测量转换。平面直角坐标系是一种直角坐标系，用于在地图上将地球表面的位置表示为平面上的二维坐标，它是一种局部坐标系，适用于较小区域的地图制图和定位。
WGS84经纬度转UTM坐标系.zip
06-27
WGS84转UTM坐标系，GPS坐标转大地坐标，Delphi实例源码。
简单的WGS84转UTM程序[C++]
HRex39的博客
07-04 2684
目录前言公式传入参数传出参数程序验证 前言 找了不少WGS84转UTM，不少都有些错误，兜兜转转找到一个即插即用的分享一下。 公式 传入参数 longitute(deg) 如:120.75819277 latitude(deg) 如:31.34534646 注意：如果传入的lon和lat是度分形式的记得要进行转换 传出参数 UTME(m) 如： 286733.498379 UTEN(m) 如： 3470047.079936 程序 #include "math.h" constexpr doubl
WGS84 转 UTM（Python代码版）
guoyc439的博客
09-09 1426
基于Python语言，使用pyproj库进行WGS 和UTM的转换。
WGS84与UTM转换
12-27
基于MFC的WGS84与UTM转换,原工程可直接运行,方便GPS相关开发者
WGS84坐标转换成UTM坐标
03-14
将wgs84格式的坐标写入到txt当中，读取之后自动计算代号，并将转换后的UTM坐标写入到txt中
WGS84与UTM转换器
12-27
基于MFC的WGS84与UTM转换器,源代码工程上传,方便GPS开发者
经纬度转UTM
08-02
经纬度转UTM投影换算，使用简单，一看就明白
UTM 坐标转换
06-10
UTM坐标转换大地坐标系转换，详细介绍了UTM到大地坐标系的转换方法，地理信息工作者适用
C# UTM坐标和WGS84坐标转换工具及源码
11-25
C# UTM坐标和WGS84坐标转换工具及源码 工具路径：UTMAndWGS84TransTool\UTMAndWGS84TransTool\bin\Debug\UTMAndWGS84TransTool.exe
WGS84坐标转化为UTM坐标系
06-17
该代码为经纬度在WGS84坐标下转化为UTM坐标的算法，真实可靠
UTM与WGS84坐标系互转C/C++实现
09-28
- `UTM_to_WGS84(int zone, double easting, double northing, double& lat, double& lon)`: 输入UTM坐标，返回WGS84的经纬度。 使用这样的库，开发者可以在C/C++程序中方便地进行坐标转换，满足GIS应用的需求。在...
笛卡尔二维坐标系转换为WGS84经纬度坐标系
最新发布
05-16
要将笛卡尔二维坐标系转换为WGS84经纬度坐标系，需要知道以下信息： 1. 笛卡尔坐标系原点的经纬度； 2. 笛卡尔坐标系中每个点的x、y坐标值； 3. 笛卡尔坐标系的投影方式。 一般情况下，我们可以使用以下步骤来进行转换： 1. 将笛卡尔坐标系中的每个点的x、y坐标值转换为WGS84坐标系中的经度和纬度值。这个转换需要根据笛卡尔坐标系的投影方式进行，常见的有UTM投影和高斯投影等。具体的转换方法可以参考相关的地图学和测量学资料。 2. 将每个点的经度和纬度值转换为度分秒表示法或者十进制表示法，方便人类理解和使用。 需要注意的是，不同的笛卡尔坐标系和投影方式可能会导致转换结果的误差。因此，在实际应用中，需要根据具体情况选择合适的转换方法，并进行误差分析和校正。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Roar冷颜

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值