热线电话：13121318867

登录
首页精彩阅读R语言批量处理TRMM降雨数据
R语言批量处理TRMM降雨数据
2018-06-10
收藏

R语言批量处理TRMM降雨数据

最近在使用TRMM遥感降雨数据，下载的源数据是bin格式的，要转成tif格式的才能供后续研究使用，下载的数据文件中提供了MATLAB与Fortran两种读取方式，一开始尝试使用MATLAB读取，发现效率很低；后来经过搜索，使用R可以非常方便快速的实现bin到tif的批量转换

1、准备工作

安装R后需要安装raster、gdal、Rsenal三个包（可能安装带有）


2、批处理脚本

setwd("你的工作空间（文件存放目录）"); library(raster) trmm <- raster(xmn=0, xmx=360, ymn=-50, ymx=50, ncol=1440, nrow=400) files <- list.files(pattern='.bin')for(finfiles) {# for all filestrmm[] <-readBin(f,'double', n=576000, size =4, endian='big') x <- flip(rotate(trmm),'y') x <- writeRaster(x, extension(f,'tif')) }

非常简单，只需要几行代码，即可实现bin->tif的批量转换

3、处理结果

TRMM数据范围为180W-180E，50S-50N；在ArcGIS中打开，可以看到转换后的tif是完全正确可用的。



4、后续计划

4.1 TRMM降雨精度评估

需要对站点降雨插值到TRMM统一分辨率进行日、月、年尺度评价。这里就带来一个问题，插值方法的选取(IDW与Co/kriging)可能IDW简单，现有采用的大都是IDW插值，但是研究表明，站点的降水量与高程显著相关。因此co/kriging（协同克里金）方法精度更高。

4.2 站点降雨插值

无论采用哪种插值方法，都需要编写批处理工具，还好通过搜索，R提供了强大的空间分析包gstat，其中就包括协同克里金方法，后面采用gstat对站点降雨进行批处理

5、总结

R是个非常有用的工具，值得总结常用的地学工具。


CDA数据分析师考试相关入口一览（建议收藏）：

▷ 想报名CDA认证考试，点击>>>   “CDA报名” 了解CDA考试详情；

▷ 想学习CDA考试教材，点击>>> “CDA教材” 了解CDA考试详情；

▷ 想加入CDA考试题库，点击>>> “CDA题库” 了解CDA考试详情；

▷ 想了解CDA考试含金量，点击>>> “CDA含金量” 了解CDA考试详情；

R语言

数据分析咨询请扫描二维码

若不方便扫码，搜微信号：CDAshujufenxi

上一篇区块链+数字货币开启“不一样的”支付方式
下一篇在Python中使用glob模块查找文件路径的方法
CDA考试动态
CDA报考指南
热门栏目
最新资讯
更多
Copyright © 2015-present 北京国富如荷网络科技有限公司 All Rights Reserved.
京ICP备11001960号-9   京公网安备 11010802034615号   经营许可证编号：京B2-20210330
客服在线
立即咨询