contents.m

Geodetic tools to be used in navigation analysis

% MatLab Geodetic Toolbox
% Version 2.9 (05.03.31)
%
% Copyright (c) 1990-2005 by Michael R. Craymer
% 26 Banting Crescent
% Kanata, Ontario, Canada  K2K 1P4
% Tel. [phone]
% Fax. [phone]
% Email: [e-mail]
%
% Angle Conversions
% deg2rad   - Degrees to radians
% dms2deg   - Degrees,minutes,seconds to degrees
% dms2rad   - Degrees,minutes,seconds to radians
% rad2deg   - Radians to degrees
% rad2dms   - Radians to degrees,minutes,seconds
% rad2sec   - Radians to seconds
% sec2rad   - Seconds to radians
%
% Coordinate Conversions
% ell2utm   - Ellipsoidal (lat,long) to UTM (N,E) coordinates
% ell2xyz   - Ellipsoidal (lat,long) to Cartesian (x,y,z) coodinates
% sph2xyz   - Shperical (az,va,dist) to Cartesian (x,y,z) coordinates
% xyz2sph   - Cartesian (x,y,z) to spherical (az,va,dist) coordinates
% xyz2ell   - Cartesian (x,y,z) to ellipsoidal (lat,long,ht) coordinates
% xyz2ell2  - xyz2ell with Bowring height formula
% xyz2ell3  - xyz2ell using complete Bowring version
%
% Coordinate Transformations
% refell    - Reference ellipsoid definition
% ellradii  - Various radii of curvature
% cct2clg   - Conventional terrestrial to local geodetic cov. matrix
% clg2cct   - Local geodetic to conventional terrestrial cov. matrix
% ct2lg     - Conventional terrestrial (ECEF) to local geodetic (NEU)
% dg2lg     - Differences in Geodetic (lat,lon) to local geodetic (NEU)
% lg2ct     - Local geodetic (NEU) to conventional terrestrial (ECEF)
% lg2dg     - Local geodetic (NEU) to differences in geodetic (lat,lon)
% direct    - Direct geodetic problem (X1,Y1,Z1 + Az,VA,Dist to X2,Y2,Z2)
% inverse   - Inverse geodetic problem (X1,Y1,Z1 + X2,Y2,Z2 to Az,VA,Dist)
% simil     - Similarity transformation (translation,rotation,scale change)
%
% Date Conversions
% cal2jd    - Calendar date to Julian date
% dates     - Converts between different date formats
% doy2jd    - Year and day of year to Julian date
% gps2jd    - GPS week & seconds of week to Julian date
% jd2cal    - Julian date to calenar date
% jd2dow    - Julian date to day of week
% jd2doy    - Julian date to year & day of year
% jd2gps    - Julian date to GPS week & seconds of week
% jd2yr     - Julian date to year & decimal year
% yr2jd     - Year & decimal year to Julian date
%
% Error Ellipses
% errell2   - Computes error ellipse semi-axes and azimuth
% errell3   - Computes error ellipsoid semi-axes, azimuths, inclinations
% plterrel  - Plots error ellipse for covariance matrix
%
% Miscellaneous
% findfixed - Finds fixed station based on 3D covariance matrix
% pltnet    - Plots network of points with labels
%
% Example Scripts
%
% DirInv    - Simple partial GUI script for direct and inverse problems
% DirProb   - Example of direct problem
% Dist3D    - Example to compute incremental 3D distances between points.
% InvProb   - Example of inverse problem
% PltNetEl  - Example plot of network error ellipses
% ToUTM     - Example of conversion from latitude,longitude to UTM

% Version History
%
% 1.0  Created
%
% 1.1  Corrected dms2rad.m & dms2deg.m when converting vectors with both
%      positive and negative angles;
%      Replaced incomplete pl2utm.m with complete version;
%      Added ConvertToUTM for example of conversion from lat,long to UTM;
%      Added caldate & juldate to convert between Gregorian calendar and
%      Julian dates.
%
% 1.2  Added Dist3D for example of computing incremental 3D distances
%      between points.
%
% 1.3  Corrected rad2deg.m function name in file;
%      Modified r2p.m to use atan instead of atan2.
%
% 1.4  Added covct2lg.m & covlg2ct.m.
%
% 1.5  Added Krassovsky & International ellipsoids to refell.m
%
% 1.6  Added xyz2plh2.m to use Bowring height formula;
%      Added xyz2plh3.m to use entire Bowring algorithm.
%
% 1.7  Added errell.m to compute error ellipse parameters
%      Added ploterrell.m to plot error ellipses
%      Added PlotNetEll example to plot network error ellipses
%
% 1.8  Added dg2lg.m & lg2dg.m to convert between geodetic & local geodetic
%      Added rad2sec.m & sec2rad to convert between radians & arc seconds
%      Renamed ct2local.m to ct2lg.m and local2ct to lg2ct
%      Modified refell.m input parameter sequence (e2 before finv)
%
% 1.9  Modified ploterrell.m to add user-defined color
% 1.9a Corrected errell.m to handle equal major/minor axes
% 1.9b Modified ploterrell.m to plot solid line
% 1.9c Modified ploterrell.m for any line type/color
% 1.9d Changed errell to errell2 for 2D error ellipses
%      Added errell3d for computing 3D error ellipsoids
%
% 2.0  Modified for MatLab 4.2c.1
% 2.0a Modified cct2clg.m & clg2cct.m to check sizes of matrices
% 2.0b Modified rad2dms to output negative angles.
% 2.0c Changed names of scripts & ploterrell (plterrel) for PC compatibility
% 2.0d Updated direct & indirect for new function names
% 2.0e Corrected xyz2ell to take abs value of dlat,dh for neg. latitudes
% 2.0f Corrected sec2rad function help description
% 2.0g Noted in comments which functions are vectorized
%      Vectorized dms2deg, dms2rad, rad2deg, rad2dms
% 2.0h Corrected juldate to refer to beginning of UT day not noon
%      Added pltnet for plotting a network of points with labels
% 2.0i Corrected rad2dms for negative angles less than 1 min of arc
%
% 2.1  Moved juldate & caldate to utils toolbox and replaced with new date
%      routines cal2jd, doy2jd, gps2jd, jd2cal, jd2dow, jd2doy, jd2gps and
%      dates
%
% 2.2  Changed rad2deg output +/- radians for compatibility with other funtions
%      Corrected help comments in dms2deg, dms2rad and rad2dms
%
% 2.3  Added rotlg2ct and rotct2lg to give rotation matrix between CT and LG
%      Corrected comments in refell, errell3
%      Correctly sorted functions in this file
% 2.4  Corrected ell2utm for long>180 & switched order of dummy args
%      Added comments in sph2xyz & xyz2sph to specify handedness of xyz system
%
% 2.5  Added Modified Airy reference ellipsoid in refell.
% 2.6  Added findfixed
% 2.7  Corrected gps2jd and jd2gps to begin GPS week numbering at 0
% 2.8  Added ellradii
%      Added TOPEX/POSEIDEN ellipsoid to refell
% 2.9  Modified ell2utm for any arbitrary central meridian (default standard
%      UTM zones)