gt_0.6.4
所属分类:GPS编程
开发工具:matlab
文件大小:35655KB
下载次数:842
上传日期:2009-09-30 11:48:05
上 传 者:
ykoky2k
说明: 日本开发的一个gps精密单点定位程序,精度不错
(Japan' s development of a precise point positioning program gps, precision, good)
文件列表:
gt_0.6.4 (0, 2009-05-01)
gt_0.6.4\bpfilt.m (1083, 2008-11-25)
gt_0.6.4\caltomjd.m (670, 2008-11-25)
gt_0.6.4\caltomjd.mexw32 (60928, 2009-05-01)
gt_0.6.4\caltomjd.mexw64 (59392, 2009-05-01)
gt_0.6.4\clktornx.m (4885, 2008-11-25)
gt_0.6.4\contents.m (6781, 2008-12-12)
gt_0.6.4\data (0, 2009-05-01)
gt_0.6.4\data\ants_alias.txt (589, 2008-11-25)
gt_0.6.4\data\contents.m (1001, 2008-11-25)
gt_0.6.4\data\dcbs_p1c1.m (6377, 2008-12-12)
gt_0.6.4\data\ephem2000.405 (1881264, 2008-11-25)
gt_0.6.4\data\geoid_egm96.mat (4157296, 2008-11-25)
gt_0.6.4\data\geoid_gsi2000.mat (8653464, 2008-11-25)
gt_0.6.4\data\gshhs_c.mat (250736, 2008-11-25)
gt_0.6.4\data\gshhs_i.mat (9508248, 2008-11-25)
gt_0.6.4\data\gshhs_l.mat (1434784, 2008-11-25)
gt_0.6.4\data\gshhs_w.mat (142248, 2008-11-25)
gt_0.6.4\data\IGS03P33_RS99.snx (4747241, 2008-12-12)
gt_0.6.4\data\igs05.atx (2067628, 2009-05-01)
gt_0.6.4\data\IGS05.snx (8357000, 2008-12-12)
gt_0.6.4\data\igs05_1400.atx (1452030, 2008-12-12)
gt_0.6.4\data\igs05_1480.atx (1824382, 2008-11-25)
gt_0.6.4\data\igs05_1499.atx (1858299, 2008-12-12)
gt_0.6.4\data\igs05_1502.atx (1860339, 2008-12-12)
gt_0.6.4\data\igs05_1525.atx (2067628, 2009-05-01)
gt_0.6.4\data\igs05_ref.txt (780, 2008-12-12)
gt_0.6.4\data\igs_01.pcv (36781, 2008-11-25)
gt_0.6.4\data\ITRF2000_GPS.SSC.txt (70854, 2008-12-12)
gt_0.6.4\data\ITRF2005_GPS.SSC.txt (104020, 2008-12-12)
gt_0.6.4\data\oload_gsi.blq (687159, 2008-11-25)
gt_0.6.4\data\oload_igs.blq (140896, 2008-11-25)
gt_0.6.4\data\prm_gpssrcs.m (14460, 2009-05-01)
gt_0.6.4\data\prm_userapc.m (335, 2008-12-12)
gt_0.6.4\data\rcvs_noncc.txt (756, 2008-12-12)
gt_0.6.4\data\rcvs_params.txt (81621, 2008-12-12)
gt_0.6.4\data\rcvs_poss.txt (70339, 2008-12-12)
gt_0.6.4\data\sats_params.txt (6375, 2009-05-01)
gt_0.6.4\data\utc_tai.txt (1406, 2008-12-12)
gt_0.6.4\dbldiff.m (1009, 2008-11-25)
... ...
function metprm=readmet(td,time,rcvs,metdir,metsrc)
%-------------------------------------------------------------------------------
% [system] : GpsTools
% [module] : read meterological parameters
% [func] : read meterological parameters
% [argin] : td = day (mjd-gpst)
% time = time vector (sec)
% rcvs = station names
% (metdir) = meteorological parameters directory (default:current)
% (metsrc) = meteorologibal parameters source (default:'mso')
% 'mso'= JMA MSM Online
% [argout] : metprm = meterological parameters
% metprm(n,1,m) = rcvs{n} time(m) pressure (hPa)
% metprm(n,2,m) = rcvs{n} time(m) temperture (C)
% metprm(n,3,m) = rcvs{n} time(m) relative humidity (%)
% [note] :
% [version]: $Revision: 12 $ $Date: 2008-11-25 10:02:15 +0900 (, 25 11 2008) $
% Copyright(c) 2004-2006 by T.Takasu, all rights reserved
% [history]: 05/06/15 0.1 new
%-------------------------------------------------------------------------------
if nargin<4, metdir=''; end
if nargin<5, metsrc='mso'; end
switch metsrc
case {'rms','msm','gsm','rso','mso','gso'}
metprm=ReadGpvData(td,time,rcvs,metdir,metsrc);
otherwise
disp(['warning : met data source error : ',metsrc])
end
% read gpv --------------------------------------------------------------------
function metprm=ReadGpvData(td,time,rcvs,metdir,metsrc)
metprm=repmat(nan,[length(rcvs),3,length(time)]);
for n=1:length(rcvs)
gpos(n,:)=eceftogeod(readpos(td,time(1),rcvs{n},'','approx')');
gpos(n,3)=geodh(gpos(n,:));
end
switch metsrc
case {'mso','gso'}
if strcmp(metsrc,'mso'), tu=6*3600; ftt=0:5; else tu=12*3600; ftt=0:6:6; end
m=0; t=(floor(time(1)/tu):floor(time(end)/tu)+1)*tu;
for n=1:length(t)
for ft=ftt
m=m+1; ts(m)=t(n)+ft*3600;
[pmsl(:,:,m),gprm]=readgpv(td,t(n),'pmsl',metdir,metsrc,0,ft);
[temp(:,:,m)]=readgpv(td,t(n),'temp',metdir,metsrc,0,ft);
[humi(:,:,m)]=readgpv(td,t(n),'humi',metdir,metsrc,0,ft);
if ts(m)>time(end), break, end
end
end
if ~isempty(gprm)
[x,y,z]=meshgrid(1:gprm.nx,1:gprm.ny,ts);
pmsl=double(pmsl);
temp=double(temp);
humi=double(humi);
for n=1:length(rcvs)
[xi,yi]=gmt('lltogrid',gpos(n,2),gpos(n,1),gprm);
pm=shiftdim(interp3(x,y,z,pmsl,xi,yi,time),2);
te=shiftdim(interp3(x,y,z,temp,xi,yi,time),2);
hu=shiftdim(interp3(x,y,z,humi,xi,yi,time),2);
pr=pm.*(1-0.0065*gpos(n,3)./(te+273.15+0.0065*gpos(n,3))).^5.257;
metprm(n,:,:)=[pr,te,hu]';
end
end
otherwise
disp(['warning : met data source error : ',metsrc])
end
近期下载者:
相关文件:
收藏者: