diff --git a/example_fit_to_the_cumulus.R b/example_fit_to_the_cumulus.R index 60c10629a56ee931d237519f22fbb5ccdc6662fe..7668ef5dc92260305aba7a6793692292a626ce37 100644 --- a/example_fit_to_the_cumulus.R +++ b/example_fit_to_the_cumulus.R @@ -23,7 +23,7 @@ cumulus$local_hour <- as.numeric(format.POSIXct(cumulus$local, "%H")) ############################################################################### cumulus_lag <- lag(cumulus) -cumulus$distance <- geosphere::distHaversine(cumulus[,2:3],cumulus_lag[,2:3]) +cumulus$distance <- geosphere::distHaversine(cumulus[,2:3],cumulus_lag[,2:3]) / 1000 cumulus$time_lag <- as.numeric(difftime(cumulus$date,cumulus_lag$date,units = "secs")) cumulus$ship.speed <- cumulus$distance / cumulus$time_lag cumulus$heading <- geosphere::bearing(cumulus[,2:3],cumulus_lag[,2:3])