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
總之,任何人都可以闡明其使處理停止的錯誤消息的任何光?
由於提前, 伊恩
喜邁克爾,這似乎已經奏效。非常感謝。 –