function?[XYZ]?=?BLh2XYZ（LATLONh）?

%?Example:
%（BLh）WGS84-（XYZ）WGS84
%?LAT?=?40.9987167395335;
%?LON?=?39.7652393428761;
%?h?=?51.403;

refell?=?1;

switch?refell
????case?1
????????%?IERS?2003?numerical?standards
????????%?ellipsoid?parameters?for?xyz2ellip.m
????????a_tidefree?=?6378136.6;?%m??????Equatorial?radius?of?the?Earth
????????f_tidefree?=?1/298.25642;?????%?Flattening?factor?of?the?Earth
????????a?=?a_tidefree;??%m??????Equatorial?radius?of?the?Earth
????????f?=?f_tidefree;???????%?Flattening?factor?of?the?Earth
????case?2
????????%?GRS?80?（http://www.bkg.bund.de/nn_164850/geodIS/EVRS/EN/References/...
????????%??Definitions/Def__GRS80-pdftemplateId=rawproperty=publication...
????????%??File.pdf/Def_GRS80-pdf.pdf）
????????a_grs80????=?6378137;
????????f_grs80????=?0.00335281068118;
????????a?=?a_grs80;???%m??????Equatorial?radius?of?the?Earth
????????f?=?f_grs80;????????%?Flattening?factor?of?the?Earth
????case?3
????????%?WGS84
????????a=6378137;
????????f=1/298.25722356;
????case?4
????????%?Hayford
????????a=6378388;
????????f=1/297;
end

b?=?a-f*a;

lat=LAT*pi/180;
lon=LON*pi/180;
e2=（a^2-b^2）/a^2;
N=a/sqrt（1-e2*（sin（lat）^2））;
X=（N+h）*cos（lat）*cos（lon）;
Y=（N+h）*cos（lat）*sin（lon）;
Z=（N*（1-e2）+h）*sin（lat）;

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????目錄???????????0??2017-07-10?16:47??MATLABcode\BLh2XYZ\
?????文件????????1290??2012-05-28?18:09??MATLABcode\BLh2XYZ\BLh2XYZ.m
?????目錄???????????0??2017-07-10?16:47??MATLABcode\XYZ2BLh\
?????文件????????1831??2012-05-22?15:41??MATLABcode\XYZ2BLh\xyz2ell.m
?????目錄???????????0??2017-07-10?16:47??MATLABcode\