Я пытаюсь написать свой первый очень базовый c скрипт на языке r. Просмотр различных сообщений и руководств не уведет меня далеко. Проблема, которую я пытаюсь решить, заключается в следующем: у меня есть следующие почасовые данные в файле:
year mo dy time ucomp vcomp temp Rhum rain
2020 07 09 00:00 -1.1 0.5 11.8 54.0 0
2020 07 09 01:00 -1.6 0.3 10.1 62.4 0
2020 07 09 02:00 -1.6 0.4 9.2 67.2 0
2020 07 09 03:00 -2.2 0.4 8.2 70.5 0
2020 07 09 04:00 -2.2 0.7 7.5 75.2 0
2020 07 09 05:00 -1.7 0.6 7.1 80.2 0
2020 07 09 06:00 -0.7 1.0 10.0 74.9 0
2020 07 09 07:00 -2.7 -1.2 15.8 61.4 0
2020 07 09 08:00 -5.4 -3.2 19.8 51.4 0
2020 07 09 09:00 -5.6 -3.7 23.3 43.4 0
2020 07 09 10:00 -5.0 -4.1 26.1 36.0 0
2020 07 09 11:00 -4.5 -4.2 28.2 28.1 0
2020 07 09 12:00 -2.9 -3.7 29.1 25.2 0
Я пытаюсь прочитать эти данные и выбрать строку, соответствующую часу утра, например, в 10h Я хочу скопировать строку
2020 07 09 10:00 -5.0 -4.1 26.1 36.0 0
Скопировав строку, я хочу построить скорость и направление ветра по формуле
windDir <- function(ucomp, vcomp) {
(180 / pi) * atan(u/v) + ifelse(v>0,180,ifelse(u>0,360,0))
}
, затем запишите направление и скорость ветра как E, NE, SE, ... (в зависимости от направления) и другие значения в другом файле.
Ниже показано, как далеко я продвинулся в разработке моего скрипта, а это очень мало.
setwd("/home/zmumba/04_R_Scripts")
library (lubridate)
fileName <- "myhourly.txt"
con <- file(fileName,open="r")
my2day<-sys.date()
while (length(oneLine <- readLines(con, n = 1, warn = FALSE, header="T")) > 0) {
I believe I have to put my code here, but how?
}
close(con)
Любая помощь будет оценена по достоинству.