使用自己的数据在openair中绘制轨迹

我想使用我自己制作的数据（使用Hysplit和R/openair）在openair软件包（R）中制作轨迹图。使用自己的数据在openair中绘制轨迹

我收到一个错误： “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)`

备用码我曾尝试 - 由“ProcTraj”创建.Rdata文件中读取作为“的OpenAir”手册附录中详细。当在Rstudio运行结束，将所得 'MYDATA' 看起来相同的（如下所示）：

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 $日期[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"))

我回到原来的错误。

MYDATA：

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

来源

2016-03-17 lost in R

似乎与your.Rdata文件内容没有问题。 mydata的日期格式不应该是一个问题。然而，你可以尝试使用：

mydata$date=as.character(mydata$date)  
mydata$date=as.POSIXct(mydata$date,tz="GMT")

来源

2016-08-26 06:29:01 Nick

感谢尼克，我刚才还看见您的回复......我结束了我实现数据分开越过国际日期变更线，并绘出了从该行的每一面的数据，这似乎工作。 –

使用自己的数据在openair中绘制轨迹

回答

相关问题