全国2400站点地面气候资料日值数据集（V3.0）数据说明与数据处理

1）中国地面气候资料日值数据集（V3.0）数据说明

"中国地面气候资料日值数据集（V3.0）"数据是按年-月存放的, 每年每月的所有站点数据每一项气象要素都放在一个文本文件中,全国站点空间位置以及每个文本的详细内容如下：

(1)SURF_CLI_CHN_MUL_DAY-TEM-12001-201501.TXT

气温数据TEM, 包括站点编号(sid), 纬度(lat), 经度(long), 高程(elev), 年(year), 月(month),日(day), 平均气温(tmean), 最高气温(tmax), 最低气温(tmin)

(2)SURF_CLI_CHN_MUL_DAY-RHU-13003-201501.TXT

湿度数据RHU, 包括站点编号(sid), 纬度(lat), 经度(long), 高程(elev), 年(year), 月(month),日(day), 平均相对湿度(rhmean), 最小相对湿度(rhmin)

(3)SURF_CLI_CHN_MUL_DAY-WIN-11002-201501.TXT

风速数据WIN, 包括站点编号(sid), 纬度(lat), 经度(long), 高程(elev), 年(year), 月(month),日(day),日平均风速、日最大风速、最大风速的风向、日极大风速、极大风速的风向

(4)SURF_CLI_CHN_MUL_DAY-SSD-14032-201501.TXT

日照数据SSD, 包括站点编号(sid), 纬度(lat), 经度(long), 高程(elev),年(year), 月(month),日(day), 日照时数(dh)

(5)SURF_CLI_CHN_MUL_DAY-PRE-13011-201501.TXT

降水数据PRE, 包括站点编号(sid), 纬度(lat), 经度(long), 高程(elev), 年(year), 月(month),日(day), 20-8 时降水量、8-20 时降水量、20-20 时降水量

(6)SURF_CLI_CHN_MUL_DAY-EVP-13240-201501.TXT

蒸发数据EVP, 包括站点编号(sid), 纬度(lat), 经度(long), 高程(elev), 年(year), 月(month),日(day), 小型蒸发量(evp1), 大型蒸发量(evp2)

(7)SURF_CLI_CHN_MUL_DAY-PRS-10004-201501.TXT

气压数据PRS, 包括站点编号(sid), 纬度(lat), 经度(long), 高程(elev), 年(year), 月(month),日(day), 平均气压(pmean), 最高气压(pmax), 最低气压(pmin)

(8)SURF_CLI_CHN_MUL_DAY-GST-12030-0cm-201501.TXT

地温数据GST, 包括站点编号(sid), 纬度(lat), 经度(long), 海拔(elev), 年(year), 月(month),日(day), 平均地表气温(stmean), 最高地表气温(stmax), 最低地表气温(stmin)

#特征值说明#：

①台站海拔高度 +100000 当台站海拔高度为估测值时，在估测数据基础上加100000

②各要素项 32766 数据缺测或无观测任务

③气压日极值 +20000 气压极值取自定时值，在原值上加20000

④日最小相对湿度 +300 最小相对湿度取自定时值，在原值上加300

⑤风速 +1000 当风速超过仪器测量上限时，在上限数据基础上加1000

⑥风向 1-17 用数字表示风向方位，17表示静风；+100 当表示风向为八风向时，在原值上加100；90X 风向出现X个时，风向数据用个数X表示；95X 风向至少出现X个时，风向数据用个数X表示

⑦降水量 32700 表示降水"微量"；32XXX XXX为纯雾露霜；31XXX XXX为雨和雪的总量；30XXX XXX为雪量(仅包括雨夹雪，雪暴）

⑧蒸发量 32700 表示蒸发器结冰； +1000 蒸发器中注入的水全部蒸发，在注入的水量数据基础上加1000

⑨ 0cm地温 +10000 实际温度（零上）超仪器上限刻度，在上限数据基础上加10000； -10000 实际温度（零下）超仪器下限刻度，在下限数据基础上减10000

文本文件的原始数据单位:

气温(0.1°C), 湿度(1%), 风速(0.1m/s), 日照(0.1h), 降水量(0.1mm), 蒸发量(0.1mm), 气压(0.1hPa), 纬度(度分), 经度(度分), 高程(0.1m)

标准数据单位：

气温(°C), 湿度(1%), 风速(m/s), 日照(h), 降水(mm), 蒸发(mm), 气压(hPa), 纬度(°), 经度(°), 高程(m)

Matlab：中国地面气象日值数据集提取所需省份全部站点数据
代码结构：
（1）Input：8个变量，每年12个月，每年96个txt文件；
（2）目标：多年循环，提取所需站点；
（3）Output：一个子变量一个表，全年数据放入一个sheet中，不同年份放入不同sheet中
输入数据命名基本为默认：

%%%%

% @LastEditors : 地理遥感生态网平台

% @Date : 2021-12-01 20:00

% @LastEditTime : 2021-12-07 22:00

% @Description : 全系列标准化气象数据的处理txt2excel;注意对每个输出结果进行检查

% @Objective : 按类输出所需省份全部年份n个气象站点数据

%%

clear;

clc;

%% 0. 前期输入输出&定义

% 输入txt根路径

txt_Path='D:\科研数据\metadata\2010_2017\';

% 输出xls根路径

xls_Path='D:\科研数据\metadata\xls2010_2017\';

% 预分配中间变量内存，列数根据元数据类别数目确定(注意按需修改)

A11=zeros(68,11);%3 RHU相对湿度；5 EVP蒸发；

% 按年读取全部txt文件

% year by year

for yy = 2010:2017

year = num2str(yy);

txt_file=fullfile(txt_Path,year,'\');

txt_list=dir(strcat(txt_file,'*.txt'));%列出路径下所有txt文件

% 判别平闰年，确定月天数，计算doy

if (mod(yy,4)==0 && mod(yy,100)~=0 || mod(yy,400)==0)

mon2day=[31 29 31 30 31 30 31 31 30 31 30 31]; % 闰年各月天数

else

mon2day=[31 28 31 30 31 30 31 31 30 31 30 31]; % 平年各月天数

end

doy=sum(mon2day);

% 根据doy预分配输出数据和表头(注意按需修改)

evp_d20=zeros(66,4+doy);

evp_E601=zeros(66,4+doy);

title{1,1}='station';

title{1,2}='lat';

title{1,3}='lon';

title{1,4}='altitude';

for dd=1:doy

title{1,4+dd}=strcat('doy',year,num2str(dd,'%03d'));

end

%% 1. 集成EVP结果

% month by month

i=1;

m=1;%第1个变量，排序根据txt_list来确定

m1=1+12*(m-1);%每个变量的第1个txt

m0=12*(m-1);

m12=12*m;%每个变量的第12个txt

for mm=m1:m12

evp_name=txt_list(mm).name;

evp_mm=load(fullfile(txt_file,evp_name));%原始txt文件读成表

n=1;

ss_1=53*mon2day(mm-m0)+1;%51053第1天

ss_2=(53+68-1)*mon2day(mm-m0)+1;%52313第1天

% day by day

for dm=0:(mon2day(mm-m0)-1)%1个月每天的循环

for ss=(ss_1+dm):mon2day(mm-m0):(ss_2+dm)%68个站点提取同一天

A11(n,:)=evp_mm(ss,:);

n=n+1;

end

A11(61,:)=[];% delete 51886 剔除不需要的站点

A11(66,:)=[];% delete 52267 注意上一个删除后位次67需减去1

evp_d20(:,1:3)=A11(1:66,1:3);

evp_E601(:,1:3)=A11(1:66,1:3);

evp_d20(:,4)=A11(1:66,4)*0.1;%altitude，0.1米

evp_E601(:,4)=A11(1:66,4)*0.1;%altitude，0.1米

evp_d20(:,4+i)=A11(:,8)*0.1;%小型蒸发量，0.1mm

evp_E601(:,4+i)=A11(:,9)*0.1;%大型蒸发量,0.1mm

i=i+1;

clear A11;

n=1;

end

end

% 经纬单位（度、分），转换为度

lat_int=floor(evp_d20(:,2)*0.01);

lat_decimal=(evp_d20(:,2)*0.01-lat_int).*100/60;

lat=lat_int+lat_decimal;

lon_int=floor(evp_d20(:,3)*0.01);

lon_decimal=(evp_d20(:,3)*0.01-lon_int).*100/60;

lon=lon_int+lon_decimal;

evp_d20(:,2)=lat;

evp_E601(:,2)=lat;

evp_d20(:,3)=lon;

evp_E601(:,3)=lon;

% 质量控制，剔除异常值（matlab固有bug，查找替换不完全，采用excel查找替换）

% evp_d20(evp_d20(:,5:end)==evp_d20(1,274))=nan;

% % evp_d20(evp_d20(:,5:end)>32765)=nan;

% evp_E601(evp_E601(:,5:end)==evp_E601(1,5))=nan;

% check1=find(evp_d20==3276.6);

% check2=find(evp_E601==3276.6);

% 输出，sheet名改为年份，所有年份，全年存在一个excel表中

writecell(title,strcat(xls_Path,'EVP_d20.xlsx'),'Sheet',year,'Range','A1');

writematrix(evp_d20,strcat(xls_Path,'EVP_d20.xlsx'),'Sheet',year,'Range','A2');

writecell(title,strcat(xls_Path,'EVP_E601.xlsx'),'Sheet',year,'Range','A1');

writematrix(evp_d20,strcat(xls_Path,'EVP_E601.xlsx'),'Sheet',year,'Range','A2');

clear doy title;

end

通过这样就可以得到逐年的所需站点的数据集了，方便又快捷！如果对您有帮助的话请点赞转发支持下！

数据来源请引用：中国日值气象站点数据集.地理遥感生态网科学数据注册与出版系统

赶紧三连关注下， 获取途径可私信。