陆地水循环分享 http://blog.sciencenet.cn/u/afantics

博文

中国地面降水/气温0.5°×0.5°格点数据集提取

已有 353 次阅读 2021-10-20 20:15 |个人分类:水文方法|系统分类:科研笔记

中国地面降水/气温0.5°×0.5°格点数据集提取

如有疑问或者获取软件、案例数据,请加QQ群790411341


1.数据检索

在中国气象数据网,找到两个数据集,以逐月的为例。通过检索研究时间,获取相应订单,可能需要相关等级的账号。

图片.png


2.数据下载

在数据筐中下载相应文件,文件为txt文件,如下所示。

图片.png

将网址复制到迅雷,可下载一系列txt文件。

图片.png


3.查找所需站点

随便打开一个txt,复制到Excel中可发现,数据为一个中国的矩阵。表头各项的含义为:

第一行"NCOLS 128"表示实体数据有128列;

第二行"NROWS 72"表示实体数据有72行;

第三行"XLLCORNER 72"表示数据最左下方格点单元的经度范围是72°-72.5°E

第四行"YLLCORNER 18"表示数据最左下方格点单元的纬度范围是18°-18.5°N

第五行"CELLSIZE 0.5"表示网格是0.5°×0.5°的;

第六行"NODATA_value -9999.0"表示中国区域以外的值用-9999.0表示。


以左下角的单元格为例,单元格的左下顶点的经纬度为(72,18),单元格边长为0.5°,因此左下角格点的经纬度为(72.25°,18.25°),或者说(72°15′,18°15′),因此我们可以通过给单元格设定相应行列标签,以快速定位所需的格点所在的行列数

图片.png图片.png

以东北松花江流域的某子流域为例,添加进图层后可见所研究流域内有3个格点,再加上左下角的一个,本次共提取这四个点的降水数据,获取四个格点的经纬度如下:

127.2542.25   127.7542.25

127.2541.75   127.7541.75

图片.png

Excel中快速定位到四个格点所在的行列,分别是2425行,111列和112列。

图片.png


4.提取数据

打开matlab,输入代码如下

close all; clear all;clc

data = dir('C:\Users\ASUS\Desktop\降水格点数据2014-2020'); %这块将文件夹的地址修正

dataname={data.name};

[s,v] = listdlg('PromptString','中国地面降水',...

    'SelectionMode','multiple','ListSize',[550,350],'liststring',dataname);

rain=zeros(1,max(size(s))-2);

for i = 1:max(size(s))

    fid = fopen(char(strcat('C:\Users\ASUS\Desktop\降水格点数据2014-2020/',dataname(i+2))),'rt'); %这块将文件夹的地址修正

    tline1=fgetl(fid);tline2=fgetl(fid);tline3=fgetl(fid);

    tline4=fgetl(fid);tline5=fgetl(fid);tline6=fgetl(fid);

    [rainall,count]=fscanf(fid,'%f',[128,72]);

    rainall=rainall'

    rain(i)=rainall(24,111); %这块改为我们需要的格点的行列数

    fclose(fid)

 end

运行,选择需要提取的文件,确定。


图片.png

运行完成后找到rain,打开,即为所需的数据,数据为一行,复制粘贴到excel中转置即可得到列数据,这样,就能得到全国任意位置的降水数据了。

图片.png

如果是1961-2013年的降水数据(逐日或逐月的),推荐下载一个数据集:

国家青藏高原科学数据中心注册,填写相关信息后可下载本数据集,

由于下载使用FTP协议,推荐使用FlashGet 下载。

数据集的名称为:Dataset of Gridded Daily Precipitation in China (1961-2013)





https://wap.sciencenet.cn/blog-3459054-1308735.html

上一篇:使用EToCalculator 计算潜在蒸散发详细教程
下一篇:《水文模型》--徐宗学

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2021-12-2 02:24

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部