1

我正在尝试使用臭氧数据进行光谱密度分析,但遇到了错误:

ts(tser, frequency = x.fsamp, start = x.start) 中的错误:“ts”对象必须有一个或多个观察值

我不确定问题是什么,因为听起来文件中的条目数有问题,但文件肯定不是空的。什么可能导致这个问题?代码如下:

library(openair)
library(psd)
filedir <-"C:/Users/dfmcg/Documents/Thesis files/ALL_GPMP_O3_Met"
myfiles <-c(list.files(path = filedir))
paste(filedir,myfiles,sep = '/')
npsfiles<-c(paste(filedir,myfiles,sep = '/'))

  for (i in npsfiles[1:3]){
    x <- substr(i,55,61)
    y<-paste(paste('C:/Users/dfmcg/Documents/Thesis files/nps images',x,sep='/'), '.png', sep='')
    png(filename = y)
    timeozone<-import(i,date="DATE",date.format = "%m/%d/%Y %H",header=TRUE,na.strings="-999")

    psdcore(timeozone, ntaper = , 
            na.action = "-999", plot = TRUE,
            refresh = TRUE)
    dev()
  }
4

0 回答 0