文件的读取
GIM文件有编辑好的函数读取dcb、sp3、rnx与o文件
read_rinex()函数
函数的输入:文件的输出位置与文件的输入位置
函数的输出:Sites_Info 其中包含了测站的name、doy、coor
-
其中的name与doy可以从文件名字直接获取
-
coor文件需要利用read_rnx3()函数进行获取
根据函数的输出分别介绍程序是如何运行的:
Sites_Info的name与doy的获取
可以看到程序直接利用文件名建立对应的name与doy
for i=1:len
%读取文件名,并在对应的位置写入文件名.mat文件
obsn=list_obs(i).name;
if length(obsn)>25
doy=obsn(15:19);
Sites_Info.name{1,i}=upper(obsn(1:4));
name=[upper(obsn(1:4)) doy '.mat'];
if exist([r_opath,'/',doy],'dir')==0
mkdir([r_opath,'/',doy]);
end
filename1=[r_opath,'/',doy,'/',name];
else
doy=[obsn(10:11),obsn(5:7)];
Sites_Info.name{1,i}=upper(obsn(1:4));
name=[upper(obsn(1:4)) doy '.mat'];
Sites_Info的coor的获取
利用read_rinx()函数进行读取,其构建函数的具体思路为:将观测值文件一行一行的读取
其他在函数中注释的比较清晰,就不再多阐述
想说以下面的代码块(部分)对不同的卫星数进行了分类讨论还是针对rimx文件的部分内容进行
G 18 C1C L1C D1C S1C C1W S1W C2W L2W D2W S2W C2L L2L D2L SYS / # / OBS TYPES
S2L C5Q L5Q D5Q S5Q SYS / # / OBS TYPES
可以看出有多少课卫星
> 2019 10 02 00 00 0.0000000 0 38 读入时间
G20 21605555.997 7 113538031.34007 -2945.184 7 46.886 21605554.710 5 34.952 21605552.267 5 88471231.36605 -2294.952 5 34.952
具体含义还没查到