一、下载数据

从中国气象数据网获得2014年的中国地面降水日值0.5°×0.5°格点数据集（V2.0），网址：http://data.cma.cn/

数据集实体文件名称：中国降水日值0.5°×0.5°格点数据集文件命名由数据集代码（SURF_CLI_CHN_PRE_DAY_GRID_0.5）、年月日标识（YYYYMMDD）组成。具体形式：SURF_CLI_CHN_PRE_DAY_GRID_0.5-YYYYMMDD.TXT

日降雨格点数据（.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.25°E ，53.75°N），第七行第二列数据网格中心为（72.75°E ，53.75°N），……，数据最后一行最后一列网格中心为（135.75°E ，18.25°N）。降水值保留1位小数。经度单位：度，纬度单位：度，格点降水单位：mm。

二、计算所需站点位置

根据所需样地的经纬度，确定样地所在格点数据网格中的位置；
例如：山东东营(118.75°,37.75°) => (34,93)

三、使用R批量处理

使用R语言将所需格点位置的数据提取出来，组成所需表格

R文件与格点数据文件放在相同目录下，提取结果如下：

所需R代码如下：

#将工作目录设置为数据文件存放的位置

setwd("F:\\qixiang")

#批量读取txt文件的内容
#list.files列出指定目录中的文件
files=list.files(pattern = ".txt")
#获得文件数量
length(files)
#对文件进行批量处理
for(i in 1:length(files))

{
#filename = SURF_CLI_CHN_PRE_DAY_GRID_0.5-20140101.txt
filename=files[i]

#filenamea = 20140101
filename_cut=substr(filename,31,38)

#读取文件名为filename的文件，不检查变量名是否有效，跳过前6行再读取数据
data=read.table(file=filename,check.names=FALSE,skip=6)

#将所需的站点数据存储到data1中
data1<-data.frame(prec=data[34,93])

#names函数：获取数据集中所有变量的名字
names(data1)<-filename_cut

#assign函数：创建变量x，将value赋值给x
assign(x=filename_cut,value=data1)

}

#ls()列出所有对象,将对象名称存储到date中
date<-ls()

#合并对象组成所需表格

qixiang<-get(date[1])

for(i in c(2:(length(date)-6)))

{

qixiang<-cbind(qixiang,get(date[i]))

}

#将所得数据写入csv文件

write.csv(t(qixiang),"qixiang2014.csv")

全年日降雨数据下载与处理教程相关推荐

GOES-16数据下载(保姆级教程)
GOES-16数据下载(保姆级教程) 先上GOES-16数据下载的网址link 1.先进入上述网址,下拉选项可以选择想要下载的数据,这里我选择的是辐亮度数据. 2.选择想要的数据的日期,最多只能选择3 ...
NASA全球GPM免费降雨数据下载（2021年11月更新）
文章目录前言一.申请注册NASA Earth Data账户二.查找数据三.开始下载 0.下载前言 1.站点管理 2.添加Edge浏览器代理总结前言 NASA提供的免费地学数据简直是造福社会 ...
一个ICESat-2数据下载的保姆教程（downthemall）
分享一个icesat-2数据下载的步骤.下载的方法多种多样,由于本人接触py很少,所以就用了downthemall插件下载,过程比较简单,py大佬不用看.筛选数据整个过程不用科学上网,下载过程中科学上 ...
Himawari8-9数据下载和预处理教程
Himawa-8是2015年7月7日凌晨2点(UTC)发射并投入运行,Himawa-9是2016年11月2日发射的,于2017年3月开始后备作业.这两个卫星都计划运行到2029年.早期阶段,Himaw ...
ERA5数据下载和批处理教程
ERA5 再分析数据是最新一代的再分析数据,由欧盟资助的哥白尼气候变化服务(C3S)创建,由 ECMWF 运营.同化了包括全球范围内不同区域和来源的遥感资料.地表与上层大气常规气象资料.它覆盖了 19 ...
ERA5日均数据下载处理
下载ERA5-Land数据及数据处理(每小时转日数据)_是一个橙子呀的博客-CSDN博客_cdsapi怎么安装 ERA5逐日资料下载方法-数据资料-气象家园_气象人自己的家园 Retrieve dai ...
CMIP6数据下载教程大全（含批量）
CMIP6数据下载教程大全(含批量) 注意:CMIP6数据网站属国外网站,比较卡,需耐心等待及其下载. 这里对一般情况下,我们用到的检索条件做个简要的介绍,Variable下选择我们需要的变量,如降雨 ...
NOAA网站气象台站数据下载教程
NOAA网站气象台站数据下载教程注意:本文最初成稿于2020年10月. 功能: 我国的气象数据下载有种种限制,但是我国的气象数据在美国NOAA网站上有共享.在2019年的时候做课题需要逐时气象数据, ...
python下载所有股票日交易数据
python下载所有股票日交易数据一.下载一.下载 #使用的变量a.b.c.d.e def DownStockData(StockCode):#引用 requests文件import reques ...
GEO芯片数据下载和探针ID转换（保姆级教程）
GEO芯片数据下载和探针ID转换(保姆级教程) 一.问题描述探针ID转换数据是否预处理过二.Rstudio的安装(建议阅读,避免后续转换时出错) 安装包的下载安装步骤三.(正文)芯片数据下载 ...

全年日降雨数据下载与处理教程

一、下载数据

二、计算所需站点位置

三、使用R批量处理

全年日降雨数据下载与处理教程相关推荐

最新文章

热门文章