这篇文章主要实现XYZ转BLH及XYZ转NEU两个功能,是为了简化一下有需要的同志的工作,原理方面不过多赘述,直接将代码附上
function [B,L,H] = xyz2blh(X,Y,Z)
%function:
% 将坐标由空间直角坐标系转换为地心坐标系
% input:
% X、Y、Z-->空间直角坐标系坐标
% output:
% B、L、H-->地心坐标系坐标
%所用常数(CGCS2000坐标系)
a=6378137;
% e2=0.00669438002290; %第一偏心率的平方
e2=0.0066943799013; %第一偏心率的平方(WGS84)
%1.计算大地经度L
if (X==0) && (Y>0)
L=90;
else