Matlab 等间隔读取IGS对流层天顶延迟（ZPD）文件

IGS(International GNSS Service)提供对流层天顶延迟（ZPD）文件，可从cddis分析中心下载（ftp://cddis.gsfc.nasa.gov/gps/products/trop_zpd），其数据间隔为300s，使用MATLAB提取2小时间隔的对流层天顶延迟数据。

以DARW站2016年第1天至第366天数据为例    （注：ZPD文件自2017年第29天起记录发生变化，使用该程序需做调整）

⑴数据提取

clear all;

clc;

day='缺失天数';         %定义变量名day，记录缺失的文件年积日

for n=1:366;                                %zpd文件天数
    ns=num2str(n,'%03d');         % 输值转字符串，固定三位数字（001,010,100）
name=['darw',ns,'0.16zpd'];     % 字符串拼接得到文件名称

ffid=fopen(name,'rt');              % 打开zpd文件，从第一行开始读取，数据存储在cell里

if ffid>0                  % 文件存在，读取数据

i = 1;
while feof(ffid) == 0    
    tline{i,1} = fgetl(ffid);
    i = i+1;
end    
% 提取间隔2h的数据，从0时至