热线电话: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是个非常有用的工具,值得总结常用的地学工具。

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

最新资讯
更多
客服在线
立即咨询