地固坐标系各坐标系的转换

最新推荐文章于 2023-12-26 16:41:26 发布
最新推荐文章于 2023-12-26 16:41:26 发布
阅读量250 收藏 1
点赞数 2
分类专栏： 大地测量学 文章标签： 算法 数学建模
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/yitian_00/article/details/132133288
版权

c++qt窗口程序下的空间直角坐标系、大地坐标系和子午面直角坐标系的相互转换

公式算法前往：

地球坐标系的各坐标系的转换算法_进击中的小龙的博客-CSDN博客

项目头文件convert.h（窗口类定义）

#pragma once

#include <QtWidgets/QWidget>
#include "ui_convert.h"

namespace Ui {
	class Form;
	class convert;
}

class convert : public QWidget
{
    Q_OBJECT

public:
    convert(QWidget *parent = nullptr);

    ~convert();
	
private  slots:
	void compute();

private:
    Ui::convert ui;
};

主函数main.cpp

#include "convert.h"
#include <QtWidgets/QApplication>

int main(int argc, char *argv[])
{
    QApplication a(argc, argv);
	convert w;
	w.show();
	return a.exec();
}

Ui设计

 窗口类构造函数，及自定义计算函数compute()的实现均在convert.cpp中

//---------convert.cpp----------
#include "convert.h"
#include"ui_convert.h"
#include<qpushbutton.h>
#include<qmessagebox.h>
#include<qdebug.h>
#include<cmath>


#define PI acos(-1)
convert::convert(QWidget *parent)
    : QWidget(parent)
{
    ui.setupUi(this);

	//计算函数调用compute函数,此处使用qt4版本的信号和槽函数，qt5可以兼容
	connect(ui.pushButton,SIGNAL(clicked()), this, SLOT(compute()));
	//清空功能
	connect(ui.pushButton_2, &QPushButton::clicked, this, [=]() {
		ui.lineEdit->clear();
		ui.lineEdit_2->clear();
		ui.lineEdit_3->clear();
		ui.lineEdit_4->clear();
		ui.lineEdit_5->clear();
		ui.lineEdit_6->clear();
		ui.lineEdit_7->clear();
		ui.lineEdit_8->clear();
		ui.lineEdit_9->clear();
		ui.lineEdit_10->clear();
		ui.lineEdit_11->clear();
		ui.lineEdit_12->clear();
		ui.lineEdit_13->clear();
		ui.lineEdit_14->clear();
		ui.lineEdit_15->clear();

	});
	//退出功能
	connect(ui.pushButton_3, &QPushButton::clicked, this, &QWidget::close);


}

convert::~convert()
{}
 
void convert::compute()
{

	double L, B, H, X, Y, Z, x, y;
	double a = 6378137;//椭圆长半轴
	double b = 6356752.3141;//椭圆短半轴
	double e_2 = 0.00669438002290;//椭圆第一偏心率e^2
	double e_dot_2 = 0.00673949677548;//椭圆第二偏心率e'^2  
	double c = a * a / b;
	double k = 1 + e_dot_2;

	int pp = 0;
	if (ui.radioButton->isChecked())
	{
		a = 6378245;//椭圆长半轴
		b = 6356863.0187730473;//椭圆短半轴
		e_2 = 0.006693421622966;//椭圆第一偏心率e^2
		e_dot_2 = 0.006738525414683;//椭圆第二偏心率e'^2
		pp = 1;
	}
	else if (ui.radioButton_2->isChecked())//
	{
		a = 6378140;//椭圆长半轴
		b = 6356755.2881575287;//椭圆短半轴
		e_2 = 0.006694384999588;//椭圆第一偏心率e^2
		e_dot_2 = 0.006739501819473;//椭圆第二偏心率e'^2
		pp = 1;
	}
	else if (ui.radioButton_3->isChecked())//
	{
		a = 6378137;//椭圆长半轴
		b = 6356752.3142;//椭圆短半轴
		e_2 = 0.00669437999013;//椭圆第一偏心率e^2
		e_dot_2 = 0.00673949674227;//椭圆第二偏心率e'^2
		pp = 1;
	}
	else if (ui.radioButton_4->isChecked())//
	{
		a = 6378137;//椭圆长半轴
		b = 6356752.3141;//椭圆短半轴
		e_2 = 0.00669438002290;//椭圆第一偏心率e^2
		e_dot_2 = 0.00673949677548;//椭圆第二偏心率e'^2
		pp = 1;
	}
	else if (pp == 0)
	{
		QMessageBox::critical(this, "warning", "Please select an ellipsoid");
	}


	if (!ui.lineEdit->text().isEmpty()&& !ui.lineEdit_2->text().isEmpty()&& !ui.lineEdit_3->text().isEmpty()&&
		!ui.lineEdit_4->text().isEmpty() && !ui.lineEdit_5->text().isEmpty() && !ui.lineEdit_6->text().isEmpty() 
		&& !ui.lineEdit_7->text().isEmpty() )//已知大地坐标系
	{
		L = (ui.lineEdit->text().toDouble() + ui.lineEdit_4->text().toDouble() / 60 + ui.lineEdit_6->text().toDouble() / 3600) / 180 * PI;
		B = (ui.lineEdit_2->text().toDouble() + ui.lineEdit_5->text().toDouble() / 60 + ui.lineEdit_7->text().toDouble() / 3600) / 180 * PI;
		H = ui.lineEdit_3->text().toDouble();
		double W = sqrt(1 - e_2 * sin(B) * sin(B));
		double V = sqrt(1 + e_dot_2 * cos(B ) *cos(B ));
		double N = a / W;
		//计算空间直角坐标系
		X = (N + H)*cos(B ) * cos(L );
		Y = (N + H) * cos(B) * sin(L );
		Z = (N*(1 - e_2) + H)*sin(B);
		ui.lineEdit_13->setText(QString::number(X,'f',4));
		ui.lineEdit_14->setText(QString::number(Y, 'f', 4));
		ui.lineEdit_15->setText(QString::number(Z, 'f', 4));
		//计算子午面直角坐标系
		x = X / cos(L);
		y = Z;
		//打印经度
		ui.lineEdit_8->setText(ui.lineEdit->text());
		ui.lineEdit_9->setText(ui.lineEdit_4->text());
		ui.lineEdit_10->setText(ui.lineEdit_6->text());
		//
		ui.lineEdit_11->setText(QString::number(x, 'f', 4));
		ui.lineEdit_12->setText(QString::number(y, 'f', 4));
	}
	else if (!ui.lineEdit_8->text().isEmpty() && !ui.lineEdit_9->text().isEmpty() && !ui.lineEdit_10->text().isEmpty() && 
		!ui.lineEdit_11->text().isEmpty() && !ui.lineEdit_12->text().isEmpty() )//已知子午面直角坐标系
	{
		L = (ui.lineEdit_8->text().toDouble() + ui.lineEdit_9->text().toDouble() / 60 + ui.lineEdit_10->text().toDouble() / 3600) / 180 * PI;
		x = ui.lineEdit_11->text().toDouble();
		y = ui.lineEdit_12->text().toDouble();
		//计算空间直角坐标系
		X = x * cos(L );
		Y = x * sin(L );
		Z = y;
		ui.lineEdit_13->setText(QString::number(X, 'f', 4));
		ui.lineEdit_14->setText(QString::number(Y, 'f', 4));
		ui.lineEdit_15->setText(QString::number(Z, 'f', 4));
		//计算大地坐标系
		double tan0 = Z /sqrt(X *X + Y * Y);
		double P = c * e_2 / sqrt(X * X + Y * Y);
		double tanB = tan0, cha;//循环迭代
		do
		{
			cha = tanB;
			tanB = tan0 + P * cha / sqrt(k + cha * cha);
		} while (abs(cha - tanB) > pow(10, -10));
		B =atan(tanB);
		double W = sqrt(1 - e_2 * sin(B) * sin(B));
		double V = sqrt(1 + e_dot_2 * cos(B) * cos(B));
		double N = a / W;
		H = sqrt(X * X + Y * Y) / cos(B) - N;
		//打印经度
		L = L * 180 / PI;
		ui.lineEdit->setText(QString::number((int)(L)));
		ui.lineEdit_4->setText(QString::number((int)((L - (int)(L)) * 60)));
		ui.lineEdit_6->setText(QString::number(((L - (int)(L)) * 60 - (int)((L - (int)(L)) * 60)) * 60));
		//打印纬度
		ui.lineEdit_2->setText(QString::number((int)(B)));
		ui.lineEdit_5->setText(QString::number((int)((B - (int)(B)) * 60)));
		ui.lineEdit_7->setText(QString::number(((B- (int)(B)) * 60 - (int)((B - (int)(B)) * 60)) * 60));
		//
		ui.lineEdit_3->setText(QString::number(H,'f',4));
	}
	else if (!ui.lineEdit_13->text().isEmpty() && !ui.lineEdit_14->text().isEmpty() && !ui.lineEdit_15->text().isEmpty())//已知空间直角坐标系
	{
		X = ui.lineEdit_13->text().toDouble();;
		Y = ui.lineEdit_14->text().toDouble();;
		Z = ui.lineEdit_15->text().toDouble();;
		//计算大地坐标系
		L = atan(Y / X);
		double tan0 = Z /sqrt(X * X + Y * Y);
		double P = c * e_2 / sqrt(X * X + Y * Y);
		double tanB = tan0, cha;//循环迭代
		do
		{
			cha = tanB;
			tanB = tan0 + P * cha / sqrt(k + cha * cha);
		} while (abs(cha - tanB) > pow(10, -10));
		B = atan(tanB);
		double W = sqrt(1 - e_2 *sin(B) * sin(B));
		double V = sqrt(1 + e_dot_2 *cos(B) * cos(B));
		double N = a / W;
		H = sqrt(X * X + Y * Y) /cos(B) - N;
		//打印经度
		
		ui.lineEdit->setText(QString::number((int)(L*180/PI)));
		ui.lineEdit_4->setText(QString::number((int)((L * 180 / PI - (int)(L * 180 / PI)) * 60)));
		ui.lineEdit_6->setText(QString::number(((L * 180 / PI - (int)(L * 180 / PI)) * 60 - floor((L * 180 / PI - (int)(L * 180 / PI)) * 60)) * 60));
		//打印纬度
		B = B * 180 / PI;
		ui.lineEdit_2->setText(QString::number((int)(B)));
		ui.lineEdit_5->setText(QString::number((int)((B - (int)(B)) * 60)));
		ui.lineEdit_7->setText(QString::number(((B - (int)(B)) * 60 - (int)((B - (int)(B)) * 60)) * 60));
		//
		ui.lineEdit_3->setText(QString::number(H,'f',4));
		//计算子午面直角坐标系
		x = X / cos(L);
		y = Z;
		L = L * 180 / PI;
		// 打印经度
		ui.lineEdit_8->setText(QString::number((int)(L)));
		ui.lineEdit_9->setText(QString::number((int)((L - (int)(L)) * 60)));
		ui.lineEdit_10->setText(QString::number(((L - (int)(L)) * 60 - (int)((L - (int)(L)) * 60)) * 60));
		//
		ui.lineEdit_11->setText(QString::number(x, 'f', 4));
		ui.lineEdit_12->setText(QString::number(y, 'f', 4));
	}

进击中的小龙
关注
专栏目录
ECI惯性坐标系转ECEF地心地固坐标系
11-06
Convert ECI (CIS, Epoch J2000.0) Coordinates to WGS 84 (CTS, ECEF)
地固系惯性系坐标转换程序
04-03
Fortran写的地固系和惯性系之间的坐标转换程序，学测绘的可能用得着！
地心地固坐标系与大地坐标系下坐标的相互转换
最新发布
weixin_63030308的博客
12-26 1140
在GPS定位计算中，地心地固直角坐标系和大地坐标系之间的坐标经常需要来回相互转换。从大地坐标系（(Φ)，(λ)，）到地心地固坐标系（）的变换公式如下： ————（1.A） ————（1.B） ————（1.C） 其中，是基准椭球体的卯酉圆曲率半径，为椭球偏心率，他们与基准椭球体的长半径和短半径存在如下关系： ————（2） ————（3） 反过来，从地心地固直角坐标（）到大地坐标系（(Φ)，(λ)，）的变换公式为 ————（4.A） ————（4.B） ————（4.C） 其中，和可分
【时间与空间】惯性系与地固系之间的转换
bear_miao的博客
04-13 3249
两个坐标系是指惯性系和地固系。 惯性系是天球参考系；(基本天球坐标系FCCS) 地固系是地球参考系。(基本地球坐标系FTCS) 做此转换的目的是卫星的运动规律在惯性系下描述，但是地球引力、固体潮等相关摄动均是地球参考系下的球谐函数展开形式，需要在地固系下进行求解。 输入为：待分析时刻的UTC时间t，表示为: 输出为：待分析时刻两个坐标系之间的转换矩阵，表示为： 具体地， 转换关系为： 其中，是极移矩阵，是地球自转矩阵，是岁差章动矩阵。 ...
两种地心地固坐标系及其相互转换
m0_68228800的博客
10-25 706
地心地固坐标系及坐标转换
地固系和惯性系之间的转换程序（ECEF2ECIF）
08-14
该程序包计算了IERS Convention 2003中由地固系至惯性系转换转换矩阵，其中岁差和章动采用IAU 2000A模型，所有程序采用Fortran语言编写。另外本程序提供示例程序以表明如何在C下调用采用该程序包。
轨道坐标系到星固坐标系的四元数转换方法-论文
07-08
绘制星体表面资源分布图,需要将深空探测器直接获取的遥感影像从轨道坐标系转换到星固坐标系。为解决通过矢量的几何关系进行转换的两步转换法和一步转换法精度受探测器轨道数据精度影响较大的问题,提出基于四元数的三...
地固直角坐标系转换到大地坐标系
03-06
地固直角坐标系转换到大地坐标系，坐标系转换程序，研究生阶段使用的，真实有效
坐标系转换.rar_matlab 地球_plateza9_地固坐标系_地球坐标转换
07-14
本压缩包提供的MATLAB程序正是为了实现这样的转换，即从地球的地心地固坐标系（Earth-Centered Earth-Fixed，ECEF）转换到经纬度高度（Latitude, Longitude, Height，LLH）坐标系。下面我们将详细探讨相关的知识点。...
Convert.rar_XYZ 转换_xyz _地固坐标系_坐标转换_经纬度转换
07-15
"Convert.rar_XYZ 转换_xyz 地固坐标系_坐标转换_经纬度转换"这个标题揭示了我们讨论的主题是关于地理坐标系统的转换，具体涉及到XYZ坐标系和经纬度之间的相互转换。XYZ坐标系常用于地球固有的直角坐标系统，而...
地心地固坐标系ECEF转地心惯性系ECI matlab代码
11-06
Convert WGS 84 (CTS, ECEF) Coordinates to ECI (CIS, Epoch J2000.0)
地心坐标系转换为地理坐标系
04-15
地心坐标系转换为地理坐标系，MATLAB程序，简单易懂，适合初学者。
GPS从入门到放弃（三） --- GPS坐标系
tyst08的专栏
09-18 1万+
GPS从入门到放弃（三） — GPS坐标系 定位就需要坐标，坐标当然是相对坐标系而言的，我们描述一个物体的位置，首先就需要建立坐标系。 按大类来分，坐标系可以分为惯性坐标系和非惯性坐标系。惯性坐标系是在空间静止或者做匀速直线运动的坐标系，其他都是非惯性坐标系。 GPS涉及到的坐标系大体有五个，在说这五大坐标系之前，我们需要先了解一些基本概念。 基本概念 地极：（Polar）地球自转轴与地球表面的...
导航坐标系：地心惯性坐标系、地心地固坐标系、当地水平坐标系、载体/机体坐标系
主要从事定位、估计、滤波、多传感器融合等研究，主要应用于雷达、红外、声纳、传感器网络等的目标（如无人驾驶、智能终端、移动机器人、传感器节点、飞机等）跟踪、定位和导航
05-19 2万+
https://baijiahao.baidu.com/s?id=1669210631913390392&wfr=spider&for=pc
惯性导航原理（1）：导航坐标系及相互转换
热门推荐
Mua111的博客
06-24 2万+
导航知识整理
7-坐标系及变换
伊伊_f
10-09 5196
坐标系简介
常用坐标系转换实现
qq_35575258的博客
12-04 1万+
坐标系的转换Matlab实现 常用的几种导航坐标系 1. 大地坐标系，WGS84（WorldGeodeticCoordinateSystem1984） 这是为GPS全球定位系统建立的坐标系统。WGS-84坐标系的原点在地球质心，Z轴指向BIH1984.0定义的协定地球极（CTP）方向，X轴指向BIH1984.0的零度子午面和CTP赤道的交点，Y轴和Z、X轴构成右手坐标系。其参数为经度、纬度、海拔高度。 其基本参数如下： 　　　　　　长半径：a=6378137±2（m）； 　　　　　　地球引力和地球质量的乘积
坐标系转换与GPS数据适应
常见的坐标系有地心地固坐标系（如WGS84）、国家坐标系以及投影坐标系等。例如，GPS系统使用的是WGS84坐标系，而许多国家有各自的本地坐标系，这些都需要相互转换。 转换过程中，需要考虑的主要因素有： 1. 椭球体...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值