如何使用Wget批量下载IGS站的o文件、精密星历和广播星历（附代码和软件的下载链接）

1.下载IGS站的o文件

1.1.使用wget_o.m文件下载所需的txt文件，更改里面对应的年份等信息。以2021年年积日1-3为例。

clc;clear;

%--------------------------------------------------------------------------
% wget_o.m文件
% 往txt文件中写入利用weget下载数据的地址
% 我只下载2个站点，比较少。
% 下载的比较多，可以弄个Excel文件，用for循环一次下载
%--------------------------------------------------------------------------

%% =======================需要修改的部分================================

%  下面地址，年份，以及时间跨度自己按需改
fidin = fopen('D:\paper_down\rinex\wget_o.txt','w');  % 建立需要下载站点的txt文件
year = '2021';       % 下载开始的年份
firstday =1 ;        % 下载的开始年积日
endday =3 ;          % 下载的结束年积日
file_name_3 = 'o.gz'; % 20年及之前后缀是Z，21年及之后后缀是gz
                 

% 需要下载的两个站点名字
file_name_1 = 'irkj';
file_name_2 = 'cusv';

%% =====================================================================

for J = firstday:endday
    doy = num2str(J,'%03d');
    prefix = ['wget ftp://igs.gnsswhu.cn/pub/gps/data/daily/',year,'/',doy,'/',year(3:4),'o/'];
    file_name = strcat(file_name_1,doy,'0.',year(3:4),file_name_3);%将文件名前缀，后缀和年份，年积日拼接在一起，得到完整的文件名
    IP =  strcat(prefix,file_name);     %将地址与文件名拼接在一起
    fprintf(fidin,'%s\n',IP);
    
end

for J = firstday:endday
    doy = num2str(J,'%03d');
    prefix = ['wget ftp://igs.gnsswhu.cn/pub/gps/data/daily/',year,'/',doy,'/',year(3:4),'o/'];
    file_name = strcat(file_name_2,doy,'0.',year(3:4),file_name_3);%将文件名前缀，后缀和年份，年积日拼接在一起，得到完整的文件名
    IP =  strcat(prefix,file_name);     %将地址与文件名拼接在一起
    fprintf(fidin,'%s\n',IP);
    
end
fclose(fidin);

1.2.生成的“wget_o.txt”文件如下。

1.3.右键把后缀名“txt”改为为“bat”（如果文件没显示后缀，自行百度“电脑文件怎么显示后缀”）。

1.4.双击“wget_o.bat”文件，即可自动下载。

1.5.查看文件下载文件，解压即可。

2.下载精密星历

2.1.精密星历使用的是周，跟o文件的年积日不同，因此下载数据首先要确定下载的时间，进入武汉大学IGS数据中心，点击数据中心，右边选择下载时间，比如我是2021-10-01到2021-10-10。下面会有“周2177 5到2179 0”。注意：我这个代码不智能，会把当周全部下载完，即2177-2179周全部都会下载。

2.2.更改“wget_igs.m”文件的周等信息，点击运行生成txt文件。剩下步骤跟1.3、1.4和1.5一致。

clc;clear;
%--------------------------------------------------------------------------
% 往txt文件中写入利用weget下载数据的地址
% wget_igs.m
%--------------------------------------------------------------------------

%% =======================需要修改的部分================================

%  下面地址，年份，以及时间跨度自己按需改

fidin = fopen('D:\paper_down\igs\wget_igs.txt','w');  % 建立需要下载站点的txt文件
year = '2021' ;     % 下载开始的年份
firstday =2138;  % 下载的开始周
endday =2139;    % 下载的结束周

%% =====================================================================
file_name_3 = '.sp3.Z';
file_name_1 = 'igs';
for J = firstday:endday
    for i = 0:6
        doy = num2str(J);
        j = num2str(i);
        prefix = ['wget ftp://igs.gnsswhu.cn/pub/gps/products/',doy,'/'];
        file_name = strcat(file_name_1,doy,j,file_name_3);%将文件名前缀，后缀和年份，年积日拼接在一起，得到完整的文件名
        IP =  strcat(prefix,file_name);     %将地址与文件名拼接在一起
        fprintf(fidin,'%s\n',IP);
    end
end
fclose(fidin);

3.下载广播星历

3.1.使用“wget_brdc.m”文件下载所需的txt文件，更改里面对应的年份等信息。以2022年年积日1-365为例。后面步骤和1.2、1.3、1.4和1.5一模一样。

clc;clear;
%--------------------------------------------------------------------------
% wget_brdc.m
% 往txt文件中写入利用weget下载数据的地址
%--------------------------------------------------------------------------

%% =======================需要修改的部分================================

%  下面地址，年份，以及时间跨度自己按需改
fidin = fopen('D:\paper_down\brdc\wget_brdc.txt','w');                   %打开文件
year = '2022'; % 下载开始的年份
firstday =1 ;  % 下载的开始年积日
endday =365 ;  % 下载的结束年积日

%% =====================================================================

file_name_3 = 'n.gz';
file_name_1 = 'brdc';
for J = firstday:endday
    doy = num2str(J,'%03d');
    prefix = ['wget http://www.igs.gnsswhu.cn/index.php/Home/Index/download.html?link=pub/gps/data/daily/',year,'/',doy,'/',year(3:4),'n/'];
    file_name = strcat(file_name_1,doy,'0.',year(3:4),file_name_3);%将文件名前缀，后缀和年份，年积日拼接在一起，得到完整的文件名
    IP =  strcat(prefix,file_name);     %将地址与文件名拼接在一起
    fprintf(fidin,'%s\n',IP);
end
fclose(fidin);

0.代码和软件下载链接https://cloud.189.cn/t/im2u6jV7fqmq （访问码：jn0d）