0

自分で作成したデータ (Hysplit と R/openair を使用) を使用して、openair パッケージ (R) で軌跡プロットを作成しようとしています。

「abs(mydata$len) のエラー: 数学関数の数値以外の引数」というエラーが表示されます。

私のコード（csvファイルから直接軌跡データを読み込む）：

#  read in trajectory
mydata <- read.csv("c:/hysplit4/working/2014/2014rml.csv", header = TRUE,
                   stringsAsFactors = FALSE  # I have tried with /without this line
                   )
# make sure dates are read correctly
mydata$date2 <- as.POSIXct(strptime(mydata$date2, format = "%Y-%m-%d %H:%M:%S"))
mydata$date <- as.POSIXct(strptime(mydata$date, format = "%Y-%m-%d %H:%M:%S"))
  as.vector(as.matrix(mydata))    # I have tried with/ without this line
#  plot trajectory for specific dates
library(openair)
library(mapdata)
trajPlot(mydata)`

私が試した代替コード - 「openair」マニュアルの付録に詳述されているように、「ProcTraj」によって作成された .Rdata ファイルを読み取ります。Rstudio で開くと、結果の 'my​​data' は同じように見えます (以下に示すように)。

mydata <- importTraj(site = "SuvaCity", year = 2014, 
local = "C:/HYSPLIT/TrajProc/")  
##  at this point it prints the trajectory data on the screen and it looks fine
trajPlot(mydata)

これにより、次のエラーが発生します。

「as.POSIXlt.character(mydata$date[1]) のエラー: 文字列が標準の明確な形式ではありません」

私が追加した場合：

mydata$date2 <- as.POSIXct(strptime(mydata$date2, format = "%Y-%m-%d %H:%M:%S"))
mydata$date <- as.POSIXct(strptime(mydata$date, format = "%Y-%m-%d %H:%M:%S"))

元のエラーに戻りました。

マイデータ:

receptor,year,month,day,hour,hour.inc,lat,lon,height,pressure,date2,date
1,2014,1,1,0,0,-18.134,178.424,18,1000.2,2014-01-01 00:00:00,2014-01-01 00:00:00
1,2013,12,31,23,-1,-18.089,178.532,30,998.6,2013-12-31 23:00:00,2014-01-01 00:00:00
1,2013,12,31,22,-2,-18.045,178.659,44.5,996.9,2013-12-31 22:00:00,2014-01-01 00:00:00
1,2013,12,31,21,-3,-18.002,178.806,61.5,995.3,2013-12-31 21:00:00,2014-01-01 00:00:00
1,2013,12,31,20,-4,-17.96,178.97,80.6,993.3,2013-12-31 20:00:00,2014-01-01 00:00:00
1,2013,12,31,19,-5,-17.92,179.151,100.5,990.4,2013-12-31 19:00:00,2014-01-01 00:00:00
1,2013,12,31,18,-6,-17.885,179.352,120.3,988,2013-12-31 18:00:00,2014-01-01 00:00:00
4

1 に答える 1

0

.Rdata ファイルの内容に問題はないようです。mydata の日付形式は問題になりません。ただし、次を使用して試すことができます。

mydata$date=as.character(mydata$date)    
mydata$date=as.POSIXct(mydata$date,tz="GMT")
于 2016-08-26T06:29:01.127 に答える