2016-08-22 60 views
0

我想知道如果任何人都可以提出建議,NCEP.interp失敗「錯誤,如果(格式(dt.f, 「%間%d%H:%M:%S」)

我試圖提取降雨率從RNCEP數據庫使用包RNCEP和函數NCEP.interp

從RNCEP數據庫中的特定樣本降雨率首先,我從.xlsx導入工作表,其中一列包含Sample_Date格式爲%Y-%m-% d%H:%M:%S

其次,我用NCEP.interp這似乎當它提供了以下錯誤並停止到18%運行偉大:

Error in if (format(dt.f, "%m-%d %H:%M:%S") > "12-31 17:59:59") { : 
missing value where TRUE/FALSE needed 

也有一系列參照下載長度,例如警告:

50: In download.file(paste("http://www.esrl.noaa.gov/psd/thredds/dodsC/Datasets/ncep.reanalysis", ... : downloaded length 522 != reported length 200 

但是,我不認爲警告是一個問題(包括完整性)。我的完整代碼是:

install.packages("RNCEP", dependencies =TRUE) 
library(RNCEP) 

SampPre6 <- read.xlsx("C:....Rainfall.xlsx", 
       sheetName = "1") 
SampPre6$Sample_Date <- strptime(as.character(SampPre6$Sample_Date), "%Y-%m-%d %H:%M:%S") 
format(SampPre6$Sample_Date, "%Y-%m-%d %H:%M:%S") 

PratePre6 <- NCEP.interp(variable='prate.sfc', level='gaussian', 
        lat=SampPre6$Latitude,lon=SampPre6$Longitude, 
        dt=SampPre6$Sample_Date,reanalysis2=TRUE, 
        keep.unpacking.info = TRUE) 

的數據的一個例子是:

Sample_ID   Sample_Date Latitude Longitude 
1  151 2013-03-15 23:29:59 31.39725 121.5068 
2  177 2013-04-13 00:14:59 31.39695 121.5067 
3  500 2013-04-13 06:59:59 23.26420 113.5147 
4  501 2013-04-13 08:15:00 23.26168 113.5214 
5  502 2013-04-13 09:00:00 23.27464 113.5529 
6  247 2013-04-26 19:44:59 22.45018 114.1511 
7  246 2013-04-26 21:14:59 22.45174 114.1346 
8  245 2013-04-26 22:30:00 22.43341 114.1168 
9  300 2013-05-18 00:59:59 31.39733 121.5067 
10  584 2013-06-04 01:15:00 23.18805 113.4910 

總之,任何人都可以闡明其使處理停止的錯誤消息的任何光?

由於提前, 伊恩

回答

1

的DT參數應該是類型「性格」這在幫助文件中指定?NCEP.interp

運行NCEP.interp(之前的),嘗試包括下面...

SampPre6 $ Sample_Date < - 格式(SampPre6 $ Sample_Date, 「%Y-%間%d%H:%M:%S」)

+0

喜邁克爾,這似乎已經奏效。非常感謝。 –

相關問題