【Matlab】 气候资料数据集预处理

图片详情及更加详细教程，请关注微信公众号：夫也的笔记
不啦不啦，今天介绍一种快速处理整理气象数据的方法，首次学起来很费劲，但是掌握以后对后面的各种任务、论文和项目应该会有帮助吧！先搞起来再说吧！声明一下啦！我不是大神，我只是从网上找到资料和教程以后，自己消化学习一下，然后发出来，当做自己的学习笔记方便回忆。文章收费是因为我也是花钱学来的，所以如果你也需要请为知识付费。

本次主要是对中国气象数据网的V3.0地面气象观测数据集进行快速处理，将原先分散的txt各个要素整合到一个列表中，一行代表一个站点。直接从全国的数据集中来提取出来自己想要的站点数据集，该过程也包含了一年缺测数据低于30天的站点数据线性内插，同时参考数据说明对异常值进行了处理。本次以新疆区域内的气象站点为例子，从全国V3.0数据集中提取新疆站点。气候数据预处理过程和代码如下：

% author email: 1154318421@qq.com
station=xlsread('D:\xinjiang\新疆.xlsx');% 一列新疆的站点数据
%
for year=2000:2017
    datasum=[];
   for mon=1:12
       if mon<10
            filename_pre=['F:\Day1951-2012yuanshi\DAY_PRE195101201612\SURF_CLI_CHN_MUL_DAY-PRE-13011-',int2str(year),'0',int2str(mon),'.TXT'];
            filename_temp=['F:\Day1951-2012yuanshi\DAY_TEM195101201612\SURF_CLI_CHN_MUL_DAY-TEM-12001-',int2str(year),'0',int2str(mon),'.TXT'];
            filename_rhu=['F:\Day1951-2012yuanshi\DAY_RHU195101201612\SURF_CLI_CHN_MUL_DAY-RHU-13003-',int2str(year),'0',int2str(mon),'.TXT'];
            filename_SSD=['F:\Day1951-2012yuanshi\DAY_SSD195101201612\SURF_CLI_CHN_MUL_DAY-SSD-14032-',int2str(year),'0',int2str(mon),'.TXT'];
            filename_WIN=['F:\Day1951-2012yuanshi\DAY_WIN195101201612\SURF_CLI_CHN_MUL_DAY-WIN-11002-',int2str(year),'0',int2str(mon),'.TXT'];
            filename_prs=['F:\Day1951-2012yuanshi\DAY_PRS195101201612\SURF_CLI_CHN_MUL_DAY-PRS-10004-',int2str(year),'0',int2str(mon),'.TXT'];
       else
            filename_pre=['F:\Day1951-2012yuanshi\DAY_PRE195101201612\SURF_CLI_CHN_MUL_DAY-PRE-13011-',int2str(year),int2str(mon),'.TXT'];
            filename_temp=['