Matlab代码批处理中国地面气象日值数据集（2400站点）及数据共享

地理遥感信息服务平台

于 2022-10-18 11:24:53 发布

阅读量1k

点赞数 1

分类专栏：气象站点数据太阳辐射文章标签：数据库

本文链接：https://blog.csdn.net/m0_66875622/article/details/127384860

版权

气象站点数据同时被 2 个专栏收录

2 篇文章 2 订阅

订阅专栏

太阳辐射

2 篇文章 0 订阅

订阅专栏

       气象数据一直是一个价值较高的数据，它被广泛用于各个领域的研究当中。地理遥感生态网平台提供的气象数据包括有气温、气压、相对湿度、降水、蒸发、风向风速、日照等多种指标，但是包含了这些全部指标的气象数据却较难获取，即使获取到了也不能随意分享。

1级目录

文件名

PRS

SURF_CLI_CHN_MUL_DAY-PRS-10004-YYYYMM.TXT（本站气压）

TEM

SURF_CLI_CHN_MUL_DAY-TEM-12001-YYYYMM.TXT（气温）

RHU

SURF_CLI_CHN_MUL_DAY-RHU-13003-YYYYMM.TXT（相对湿度）

PRE

SURF_CLI_CHN_MUL_DAY-PRE-13011-YYYYMM.TXT（降水）

EVP

SURF_CLI_CHN_MUL_DAY-EVP-13240-YYYYMM.TXT（蒸发）

WIN

SURF_CLI_CHN_MUL_DAY-WIN-11002-YYYYMM.TXT（风向风速）

SSD

SURF_CLI_CHN_MUL_DAY-SSD-14032-YYYYMM.TXT（日照）

GST

SURF_CLI_CHN_MUL_DAY-GST-12030-0cm-YYYYMM.TXT(0cm地温)

     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