目录
1、前言
在地球表面进行刚体变换时候,要将具有经纬度、高程和偏北角的坐标信息转换为变换矩阵表达,首先需要了解坐标系之间的转换关系。
通常,我们会将经纬度坐标转换为地心坐标系(ECEF坐标系),然后再根据高程和偏北角进行变换。
2、示例
#include <iostream>
#include <Eigen/Dense>
#include <cmath>
// 定义常量
constexpr double PI = 3.14159265358979323846;
constexpr double EARTH_RADIUS = 6378137.0; // 地球半径,单位:米
// 经纬度转ECEF坐标
Eigen::Vector3d geodeticToECEF(double latitude, double longitude, double altitude) {
double cosLat = cos(latitude * PI / 180);
double sinLat = sin(latitude * PI / 180);
double cosLon = cos(longitude * PI / 180);
double sinLon