gt_0.6.3
所属分类：GPS编程
开发工具：matlab
文件大小：27371KB
下载次数：28
上传日期：2011-05-04 18:03:45
上 传 者：
qwer_ws
说明： 进行GPS数据处理的软件,shi yong matlab编程的一个具有很强功能的GPS处理软件
( for the GPS it is a GPS solutiong software,it come form the project of GPS with some university)
文件列表:
gt_0.6.3\bpfilt.m (1053, 2006-07-08)
gt_0.6.3\caltomjd.dll (31232, 2006-07-08)
gt_0.6.3\caltomjd.m (640, 2006-07-08)
gt_0.6.3\clktornx.m (4855, 2006-07-21)
gt_0.6.3\contents.m (6449, 2006-07-08)
gt_0.6.3\data\ants_alias.txt (559, 2006-07-08)
gt_0.6.3\data\contents.m (1001, 2006-07-08)
gt_0.6.3\data\dcbs_p1c1.m (5807, 2006-07-08)
gt_0.6.3\data\ephem2000.405 (1881264, 2002-03-07)
gt_0.6.3\data\geoid_egm96.mat (4157296, 2005-06-19)
gt_0.6.3\data\geoid_gsi2000.mat (8653464, 2005-06-17)
gt_0.6.3\data\gshhs_c.mat (250736, 2005-03-30)
gt_0.6.3\data\gshhs_i.mat (9508248, 2005-03-30)
gt_0.6.3\data\gshhs_l.mat (1434784, 2005-03-30)
gt_0.6.3\data\gshhs_w.mat (142248, 2005-03-30)
gt_0.6.3\data\igs_01.pcv (36781, 2004-09-17)
gt_0.6.3\data\igs_05.atx (1357362, 2006-05-08)
gt_0.6.3\data\oload_gsi.blq (687159, 2005-04-24)
gt_0.6.3\data\oload_igs.blq (140896, 2006-04-17)
gt_0.6.3\data\prm_gpssrcs.m (12023, 2006-07-01)
gt_0.6.3\data\rcvslist_igs_hr_58.txt (348, 2006-07-20)
gt_0.6.3\data\rcvs_noncc.txt (660, 2006-07-08)
gt_0.6.3\data\rcvs_params.txt (80624, 2006-07-07)
gt_0.6.3\data\rcvs_poss.txt (70174, 2006-07-08)
gt_0.6.3\data\sats_params.txt (4245, 2006-07-07)
gt_0.6.3\data\utc_tai.txt (1273, 2006-07-06)
gt_0.6.3\data\vssver.scc (368, 2006-07-22)
gt_0.6.3\dbldiff.dll (33792, 2006-07-08)
gt_0.6.3\dbldiff.m (979, 2006-07-08)
gt_0.6.3\eceftogeod.dll (36352, 2006-07-08)
gt_0.6.3\eceftogeod.m (626, 2006-07-08)
gt_0.6.3\eceftosatf.m (1071, 2006-07-08)
gt_0.6.3\ecsftoecef.dll (442368, 2006-07-08)
gt_0.6.3\ecsftoecef.m (1295, 2006-07-08)
gt_0.6.3\ecsftosatf.dll (30720, 2006-07-08)
gt_0.6.3\ecsftosatf.m (745, 2006-07-08)
gt_0.6.3\editlist.m (11713, 2006-07-20)
gt_0.6.3\editmap.m (7839, 2006-07-08)
gt_0.6.3\eletostate.dll (33792, 2006-07-08)
gt_0.6.3\eletostate.m (820, 2006-07-08)
... ...
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: 2 $ $Date: 06/07/08 14:21 $
% 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
近期下载者：
相关文件：
收藏者：