%將雷達坐標轉到大地坐標系下
%

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）;?%坐標的