IGRAv2（探空站）数据介绍三：Matlab读取IGRA探空站数据

0. 码字不易，点赞加关注WZZHHH，详细代码在最后获取。

1.数据下载

IGRA下载地址https://www1.ncdc.noaa.gov/pub/data/igra/data/data-por/。IGRA数据的的介绍在我前面的博文里面IGRAv2（探空站）数据介绍二（附txt翻译下载链接和如何下载）-CSDN博客和IGRAv2（探空站）数据介绍一（不断更新，看到哪，更到哪）-CSDN博客

2.数据形式

IGRA为txt文件，按照年年份每日两次或者三次依次排列。我以2021年为例对气压、位势高度、温度、相对湿度、露点温度、风向和风速进行提取。

3.保存形式

提取数据为两列元胞数组存到mat文件中，每个站点存一个mat，txt名称为站点名字-年份。 

元胞数组：第一列是数据（气压、位势高度、温度、相对湿度、露点温度、风向、风速）；第二列为时间(年-月-日-小时)，如下图。

4.详细且免费的代码请关注WZZHHH，下面是部分代码展示。

% 读取IGRA的文件
% 基于matlab2022写的
% IGRA下载地址https://www1.ncdc.noaa.gov/pub/data/igra/data/data-por/
% 里面有很多站点已经弃用，并且该网址包括建站以后的所有数据集


% 提取数据存到"data_igra"，按照时间进行排列，最后存到mat文件中
% 每个站点存一个txt，txt名称为站点名字-年份
% 左边是数据（气压、位势高度、温度、相对湿度、露点温度、风向、风速）
% 右边是时间(年-月-日-小时)
clc,clear;


%% -----------------------------需要修改的参数--------------------------
% IGRA路径
path_igra = 'D:\DATA\1\';

% 数据存放路径
save_igra = 'D:\DATA\2\';

% IGRA数据包含了建站以来数据，挑选所需年份
year = '2021'; % 以2021年为例

%% ------------------------------------------------------------------------


% 路径下所有的IGRA的txt文件
igra_d = dir([path_igra,'*.txt']);

% 对a进行赋值，便于读取对应年份所需的数据
aa=1000;
n=1;
% 对txt文件逐个读取
for i = 1:size(igra_d,1)

    fid = fopen([path_igra igra_d(i).name]);

    % 逐行读取
    
    while ~feof(fid)
        line = fgetl(fid);

        % 依次读取每行，字符找到年
        % 年字符位置（14-17）、月（19-20）、日（22-23）、小时（25-26）
        % IGRA详情请阅读博客https://blog.csdn.net/WZZHHH/article/details/135239582
        if contains(line(14:17),year)
            % 只有所需年份时候aa为0