熱線電話: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是個非常有用的工具,值得總結常用的地學工具。

數據分析咨詢請掃描二維碼

若不方便掃碼,搜微信號:CDAshujufenxi

數據分析師資訊
更多

OK
客服在線
立即咨詢
日韩人妻系列无码专区视频,先锋高清无码,无码免费视欧非,国精产品一区一区三区无码
客服在線
立即咨詢