koordtrf.m

这是国外关于卫星导航方面一书的源代码

function []=koordtrf(pkt)
% -----------------------------------------------------------------------------
% GPSLab: KOORDTRF.M          coordinate-transformation from/to
%                             - geocentric-cartesian
%                             - geodetic-geographic
%                             - Gauss-Krueger
%                             bzgl. der Ellipsoide WGS84, GRS80, Bessel
%
% calls: BL2GK.M, GK2BL.M, XYZ2BLH.M, BLH2XYZ.M
%
% loads the files S1HEAD und S2HEAD for default values
% -----------------------------------------------------------------------------
% (c) iapg 1998-2002 zeb

if nargin==0
   pkt=1;
end

% Punktbezug f黵 default-Wert

if pkt==1
   load s1head;
elseif pkt==2
   load s2head;
   s1head=s2head;
   clear s2head;
else
         errordlg(['wrong input parameter.' ...
            ' value has to be 1 (Reference) or 2 (Rover).'], ...   
   						 'GPSLab: Break');
         return;
end


% Ellipsoidauswahl

what_ell=questdlg('select reference ellipsoid:', ...
                         'GPS Lab: coordinate trafo', ...
                         'WGS84','GRS80','Bessel','cancel');
   
     switch what_ell,
       case 'WGS84', 
         ell=('wgs84');
       case 'GRS80',
         ell=('grs80');
       case 'Bessel',
         ell=('bessel');
      case 'cancel',
         errordlg(['break at ellipsoid selection in coordinate trafo.' ...
            ' Due to this fact the coordinate trafo was stopped.'], ...   
   						 'GPSLab: break');
         return;
     end % switch

% Typ der Koordinaten-Trafo
     
     what_trafo=questdlg('select a coordinate trafo:', ...
                         'GPS Lab: coordinate trafo', ...
                         'XYZ > BLh + GK','BLh > XYZ + GK','GK > BLh + XYZ','cancel');
   
     switch what_trafo,
        
     case 'XYZ > BLh + GK', 
          
prompt={'X [m]','Y [m]','Z [m]','reference meridian [deg]'};
def={num2str(s1head(1)),num2str(s1head(2)),num2str(s1head(3)),'12'};
titel='GPS Lab: coordinate trafo';
lineNo=[1;1;1;1];
answer=inputdlg(prompt,titel,lineNo,def);

if size(answer)==[0 0] % Der Fall, da