Matlab代码批处理中国地面气象日值数据集(2400站点)及数据共享
气象数据一直是一个价值较高的数据,它被广泛用于各个领域的研究当中。地理遥感生态网平台提供的气象数据包括有气温、气压、相对湿度、降水、蒸发、风向风速、日照等多种指标,但是包含了这些全部指标的气象数据却较难获取,即使获取到了也不能随意分享。
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
通过这样就可以得到逐年的所需站点的数据集了,方便又快捷!如果对您有帮助的话请点赞转发支持下!
数据来源请引用:中国日值气象站点数据集.地理遥感生态网科学数据注册与出版系统
Matlab代码批处理中国地面气象日值数据集(2400站点)及数据共享相关推荐
- Matlab代码批处理全国气象站点日值数据集(2400站点)
1级目录 文件名 PRS SURF_CLI_CHN_MUL_DAY-PRS-10004-YYYYMM.TXT(本站气压) TEM SURF_CLI_CHN_MUL_DAY-TEM-12001-YYYY ...
- 全国2400站点地面气候资料日值数据集(V3.0)数据说明与数据处理
(1)中国地面气候资料日值数据集(V3.0)数据说明 "中国地面气候资料日值数据集(V3.0)"数据是按年-月存放的, 每年每月的所有站点数据每一项气象要素都放在一个文本文件中,全 ...
- 【数据集处理】中国地面气候资料日值数据集(V3.0)
1 数据集介绍-中国地面气候资料日值数据集(V3.0) 中国地面气候资料日值数据集(V3.0)包含以下气象数据资料(包括降水.气温.风速等). 2 数据处理 以湿度数据(RHU) 处理为例,数据文件如 ...
- 利用python处理中国地面气候资料日值数据集(V3.0)
原始文件 数据链接 00年到18年气象数据,但是未进行处理-数据集文档类资源-CSDN下载00年到18年气象数据,但是未进行处理,如果需要处理好的数据(execl和shp),请查看本博主其他更多下载资 ...
- 中国地面气候资料日值数据集(V3.0)
1.数据集信息 数据集中文名称:中国国家级地面气象站基本气象要素日值数据集(V3.0) 数据集代码:SURF_CLI_CHN_MUL_DAY 数据集版本:V3.0 2.数据文件命名由数据集代码(SUR ...
- 中国地面气候资料日值数据集(V3.0)——数据说明与处理
一.数据说明 "中国地面气候资料日值数据集(V3.0)"包含了中国699个基准.基本气象站1951年1月以来本站气压.气温.降水量.蒸发量.相对湿度.风向风速.日照时数和0cm地温 ...
- 中国地面气候资料日值数据集(V3.0)数据说明以及数据处理
(1)<中国地面气候资料/日值数据集(V3.0)>数据说明 "<中国地面气候资料/日值数据集(V3.0)>"数据是按年-月存放的, 每年每月的所有站点数据每 ...
- 利用中国区域地面气象要素驱动数据集(CMFD)和CN05.1气候数据格点化气象数据驱动Biome-BGC在区域上进行模拟
Biome-BGC是利用站点描述数据.气象数据和植被生理生态参数,模拟日尺度碳.水和氮通量的有效模型,其研究的空间尺度可以从点尺度扩展到陆地生态系统. 在Biome-BGC模型中,对于碳的生物量积累, ...
- 独家塞罕坝1951-2012中国国家级地面气象站基本气象要素日值数据集
独家塞罕坝1951-2012中国国家级地面气象站基本气象要素日值数据集 加Q2657801430
最新文章
- LeetCode简单题之相交链表
- 身份证号码验证-支持新的带x身份证
- C# 系统应用之清除Cookies、IE临时文件、历史记录 转载
- iPhone Three20软件引擎之构建开发环境与HelloWorld
- js遮罩层以及移动端的上拉框
- DBMS-数据库设计与E-R模型:E-R模型、约束、E-R图、E-R扩展特性、E-R图转换为关系模式、UML建模...
- Jenkins常用插件
- git 移动分支指针_理解git 中的HEAD指针branch指针
- adb shell操作文件
- Hive学习笔记(二)——数据模型
- Java代理的几种方式
- 识别合格ToB产品经理
- ByteBuffer详解(大概2333)
- Clojure学习笔记(一)——介绍、安装和语法
- vivo X系列为什么能不断拉高手机上限?
- (正则)校验 8-16位,必须含有特殊字符、而大写字母、小写字母、数字至少包含其中两项
- css3 logo 自上而下动画 渐渐出现
- TA游戏推荐:Android益智游戏《戳青蛙》
- 工厂模式与超级英雄的那些事儿
- html中的float含义,float是什么意思?