Code covered by the BSD License

# Geodetic Toolbox

### Mike Craymer (view profile)

13 Jun 2007 (Updated )

Toolbox for angle, coordinate and date conversions and transformations. Version 2.97.

Clg=cct2clg(Cct,lat,lon)
```function Clg=cct2clg(Cct,lat,lon)
% CCT2CLG  Convert CT covariance matrix to local geodetic.
% Version: 2011-04-05
% Useage:  Clg=cct2clg(Cct,lat,lon)
% Input:   Cct - CT covariance matrix
%          lat - vector of station latitudes (rad)
%          lon - vector of station longitudes (rad)
% Output:  Clg - LG covariance matrix

% Copyright (c) 2011, Michael R. Craymer
% Email: mike@craymer.com

n=length(lat);
if (n*3 ~= max(size(Cct)) )
error('Size of lat,lon does not match size of Cct');
end

for i=1:n
sinlat=sin(lat(i));
coslat=cos(lat(i));
sinlon=sin(lon(i));
coslon=cos(lon(i));
Ji=rotct2lg(lat(i),lon(i));
indi=(i-1)*3+[1:3];
for j=1:n
sinlat=sin(lat(j));
coslat=cos(lat(j));
sinlon=sin(lon(j));
coslon=cos(lon(j));
Jj=rotct2lg(lat(j),lon(j));
indj=(j-1)*3+[1:3];
Clg(indi,indj)=Ji*Cct(indi,indj)*Jj';
end
end
```