当我运行以下R代码块时:
require(openair)
require(png)
topDir <- "C:/Users/djh/Desktop/WindRoses"
subdirs <- c("Abbotsford_Observations") #, "Vancouver_Observations", "Abbotsford_Modelled", "Vancouver_Modelled")
years <- c(1985) #, 1995, 2001, 2006)
for(i in 1:length(subdirs)){
for(j in 1:length(years)){
wd <- paste(topDir, subdirs[i], years[j], sep="/")
files <- list.files(wd, pattern = "\.out$")
for(k in 1:length(files)){
theData <- data.frame(read.table(paste(wd, files[k], sep="/"), header = TRUE, sep=""))
u <- theData$U10
v <- theData$V10
theData["windSpd"] <- sqrt(u^2 + v^2)
theData["windDir"] <- (270 - (atan2(u/theData$windSpd, v/theData$windSpd)*(180/pi)))
nameSplit <- strsplit(files[k], ".")
png(file=paste(wd, "/", nameSplit[[1]], ".png", sep = ""))
windRose <- windRose(theData, theData$windSpd, theData$windDir, angle = 22.5)
dev.off()
}
}
}
我收到错误:
在整个代码运行后,".subset2(x,i,exact=exact(中的错误:在级别1没有这样的索引">
看看堆栈交换中这个错误的其他实例,它似乎可能与files[k]
的字符串拆分有关,但没有一个答案能解决我的问题。
注意:我曾尝试在nameSplit
上使用unlist
,但它没有解决问题。
我试图拆分的字符串的一个例子是:
wrfout_d04_1985-07-16.ts.abbotsford.out
这肯定是错误发生的地方吗?还是脚本中的其他地方有明显的原因?
我认为,在用于拆分文件名的regex定义中,点应该用双反斜杠屏蔽,就像在分配files
变量的行中一样。此外,如果以后要使用nameSplit[1]
选择拆分字符串的第一部分,则unlist()
命令可能会很有用。
因此,我建议您使用
nameSplit <- unlist(strsplit(files[k], "\."))
看看问题是否仍然存在。
希望这能有所帮助。
问题实际上来自另一行。openair程序包中的windRose命令要求输入风速和风向,因为它们在下一行中,而不是它们在原始问题中的显示方式。
wind_rose <- windRose(theData, ws="windSpd", wd="windDir", angle = 22.5)