WGS-84坐标系转大地坐标系

%将雷达坐标转到大地坐标系下
%

function [BBBLLLHHH]=zt_xyz_BLH（xi_Axi_Bxi_Cwidth）


% clear;
% clc;
% 
% A=geotiffread（‘cut1.tif‘）;
% info=geotiffinfo（‘cut1.tif‘）;
% [widthlongth]=size（A）;
% [filenameObspathnameObs]=uigetfile（‘*.*par‘‘当天的par文件‘）;
% [xi_Axi_Bxi_CT_ST_E]=zt_nihe（filenameObspathnameObs）; 



% A=geotiffread（‘cut1.tif‘）;
% info=geotiffinfo（‘cut1.tif‘）;
% [widthlongth]=size（A）;
% BLH=zt_BLH（Ainfowidthlongth）;    %脚本命令
% BBB=BLH（::1）;
% LLL=BLH（::2）;
% HHH=BLH（::3）;
% [zt_XYZN]=zt_BLH_xyz（BLHwidthlongth）;   %脚本命令



%function [zt_XYZN]=zt_xyz_BLH（BLHwidthlongth）


ee=0.00669437999013; 
a=6378137;  


% zt_XYZ_x=zeros（widthlongth）;
% zt_XYZ_y=zeros（widthlongth）;
% zt_XYZ_z=zeros（widthlongth）;
% B=zeros（widthlongth）;
% L=zeros（widthlongth）;
% H=zeros（widthlongth）;
% for i=1:width
%     for j=1:longth        
%  %将经度和纬度转换为弧度
%  BB=BLH（ij1）*pi/180;
%  LL=BLH（ij2）*pi/180;
% 
% %将直角坐标系转换为大地坐标系 
% N1=a/sqrt（1-ee*sin（BB）*sin（BB））;  %椭球的卯酉圈的半径计算
% 
% zt_XYZ_x（ij）=（N1+BLH（ij3））*cos（BB）*cos（LL）; %坐标的x
% zt_XYZ_y（ij）=（N1+BLH（ij3））*cos（BB）*sin（LL）; %坐标的