function?ELL=cart2ell（CARTellipsFileOut）

%?CART2ELL?performs?transformation?from?cartesian?coordinates?to?ellipsoidal?coordinates?
%?
%?ELL=cart2ell（CARTellipsFileOut）
%?
%?Also?necessary:???Ellipsoids.mat???（see?beneath）
%?
%?Inputs:??CART??Right-handed?cartesian?coordinates?as?nx3-matrix?（XYZ）?[m]
%????????????????3xn-matrices?are?allowed.?Be?careful?with?3x3-matrices!
%????????????????CART?may?also?be?a?file?name?with?ASCII?data?to?be?processed.?No?point?IDs?only
%????????????????coordinates?as?if?it?was?a?matrix.
%
%????????ellips??The?underlying?ellipsoid?as?string?in?lower?case?letters?default?if?omitted?or?set
%????????????????to?[]?is?‘besseldhdn‘
%????????????????See?Ellipsoids.m?for?details.
%???????????????????
%???????FileOut??File?to?write?the?output?to.?If?omitted?no?output?file?is?generated.
%
%?Outputs:??ELL??nx3-matrix?with?ellipsoidal?coordinates?（longitude?latitude?ell.?height）
%????????????????in?[degree?m]
%????????????????Southern?hemisphere?is?signalled?by?negative?latitude.

%?Author:
%?Peter?Wasmeier?Technical?University?of?Munich
%?p.wasmeier@bv.tum.de
%?Jan?18?2006

%%?Do?some?input?checking

%?Load?input?file?if?specified
if?ischar（CART）
????CART=load（CART）;
end

%?Check?input?sizes
if?????（size（CART1）~=3）&&（size（CART2）~=3）?error（‘Coordinate?list?CART?must?be?a?nx3-matrix!‘）
elseif?（size（CART1）==3）&&（size（CART2）~=3）?CART=CART‘;
end

%?Defaults
if?nargin<3?FileOut=‘‘;end
if?nargin<2?||?isempty（ellips）?ellips=‘besseldhdn‘;end

%?Load?ellipsoids
load?Ellipsoids;
if?~exist（ellips‘var‘）?error（[‘Ellipsoid?‘ellips‘?is?not?defined?in?Ellipsoids.mat?-?check?your?definitions!.‘]）
end
eval（[‘ell=‘ellips‘;‘]）;

%%?Do?calculations

ELL=zeros（size（CART））;
%?Longitude
ELL（:1）=atan2（CART（:2）CART（:1））*180/pi;
ELL（ELL（:1）<0:）=ELL（ELL（:1）<0:）+360;

%?Latitude
B0=atan2（CART（:3）sqrt（CART（:1）.^2+CART（:2）.^2））;
B=100*ones（size（B0））;
e2=（ell.a^2-ell.b^2）/ell.a^2;
while（any（abs（B-B0）>1e-10））
????N=ell.a./sqrt（1-e2*sin（B0）.^2）;
????h=sqrt（CART（:1）.^2+CART（:2）.^2）./cos（B0）-N;
????B=B0;
????B0=atan（（CART（:3）./sqrt（CART（:1）.^2+CART（:2）.^2））.*（1-e2*N./（N+h））.^（-1））;
end
ELL（:2）=B*180/pi;
ELL（:3）=h;

%%?Write?output?to?file?if?specified

if?~isempty（FileOut）
????fid=fopen（FileOut‘w+‘）;
????fprintf（fid‘%12.10f??%12.10f??%12.10f\n‘ELL‘）;
????fclose（fid）;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2420??2018-11-16?08:06??GeodeticToolbox\cart2ell.m
?????文件????????6035??2018-11-16?08:06??GeodeticToolbox\Contents.m
?????文件????????2761??2018-11-16?08:06??GeodeticToolbox\d1trafo.m
?????文件????????2858??2018-11-16?08:06??GeodeticToolbox\d2affinetrafo.m
?????文件????????3360??2018-11-16?08:06??GeodeticToolbox\d2projectivetrafo.m
?????文件????????3016??2018-11-16?08:06??GeodeticToolbox\d2trafo.m
?????文件????????2991??2018-11-16?08:06??GeodeticToolbox\d3affinetrafo.m
?????文件????????2577??2018-11-16?08:06??GeodeticToolbox\d3projectivetrafo.m
?????文件????????4447??2018-11-16?08:06??GeodeticToolbox\d3trafo.m
?????文件????????1412??2018-11-16?08:06??GeodeticToolbox\deg2dms.m
?????文件????????1970??2018-11-16?08:06??GeodeticToolbox\dms2deg.m
?????文件????????2433??2018-11-16?08:06??GeodeticToolbox\ell2cart.m
?????文件????????4519??2018-11-16?08:06??GeodeticToolbox\ell2lambertcc.m
?????文件????????7385??2018-11-16?08:06??GeodeticToolbox\ell2tm.m
?????文件???????11494??2018-11-16?08:06??GeodeticToolbox\ell2utm.m
?????文件?????????550??2018-11-16?08:06??GeodeticToolbox\Ellipsoids.m
?????文件????????1086??2018-11-16?08:06??GeodeticToolbox\Ellipsoids.mat
?????文件??????810614??2018-11-16?08:06??GeodeticToolbox\Geodetic?Transformations?Toolbox?-?Examples?and?Background.pdf
?????文件????????5134??2018-11-16?08:06??GeodeticToolbox\helmert1d.m
?????文件????????7162??2018-11-16?08:06??GeodeticToolbox\helmert2d.m
?????文件???????14052??2018-11-16?08:06??GeodeticToolbox\helmert3d.m
?????文件????????4350??2018-11-16?08:06??GeodeticToolbox\helmertaffine2d.m
?????文件????????5156??2018-11-16?08:06??GeodeticToolbox\helmertaffine3d.m
?????文件????????4500??2018-11-16?08:06??GeodeticToolbox\helmertprojective2d.m
?????文件????????4623??2018-11-16?08:06??GeodeticToolbox\helmertprojective3d.m
?????文件???????10821??2018-11-16?08:06??GeodeticToolbox\itrstrafo.m
?????文件????????4618??2018-11-16?08:06??GeodeticToolbox\lambertcc2ell.m
?????文件????????1314??2018-11-16?08:06??GeodeticToolbox\license.txt
?????文件????????5188??2018-11-16?08:06??GeodeticToolbox\molodenskytrafo.m
?????文件????????6976??2018-11-16?08:06??GeodeticToolbox\ntv2trafo.m
?????文件????????5141??2018-11-16?08:06??GeodeticToolbox\Projections.m
............此處省略10個文件信息