matlab大地坐标转换程序

利用matlab编写的大地坐标BLH和空间直角坐标XYZ转换程序

function [ XYZ ] = LBH2XYZ（ afLBH ）
%  此函数用于计算某点的空间直角坐标
%  a为椭球的长半轴长度，f为扁率
%  LB分别为某点的大地经纬度，采用弧度制
%  H为大地高，单位为米
%  XYZ为空间直角坐标，单位为米

L=L/180*pi;
B=B/180*pi;
f=1/f;
e=sqrt（2*f-f^2）;
N=a/（sqrt（1-e*e*sin（B）*sin（B）））;
X=（N+H）*cos（B）*cos（L）;
Y=（N+H）*cos（B）*sin（L）;
Z=（N*（1-e^2）+H）*sin（B）;

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         383  2016-05-04 21:09  matlab大地坐标正反算\LBH2XYZ.m
     文件        1063  2016-04-24 17:22  matlab大地坐标正反算\XYZ2LBH.m