File:Tulilahti Finland illumination night 27 28 07 1959 1.svg

1. 1. illumination at Tulilahti, Heinävesi, Finland
2. famous double murder site and time, night 27-18.1959
   1. times UTC, not local time
3. "R" code
   1. uses "moonlit" package
5. 13.8.2023 0000.0002

1. install and load devtools
2. install.packages("devtools")
3. library(devtools)

1. install moonlit library from github repo
2. install_github("msmielak/moonlit")

1. https://cran.r-project.org/web/packages/suncalc/suncalc.pdf
2. in POSIXct

1. load the moonlit library

library(moonlit)

library(suncalc)

library(ggplot2)

library(lubridate)

1. 1. heinävesi tulilahden murha

paikka1="Tulilahti, Heinävesi"

lon1=28.630556 lat1=62.426389 date1="1959-07-27" date2="1959-07-28" datetime1="1959-07-27 21:59:59" timezone1="UTC"

1. timezone1="Europe/Helsinki"

e1=0.26 t1="1 mins"

1. "1959-7-27 20:00"
2. "1959-7-28 2:00"

datetimes2<-seq(as.POSIXct("1959-7-27 18:00", tz = "UTC"),

as.POSIXct("1959-7-28 2:00", tz = "UTC"), by = "min")

1. print(datetimes2)

1. stop(-1)

1. print("Sun " )

sunpos1<-getSunlightPosition(date = as.POSIXct(datetimes2, tz = "UTC"), lat = lat1, lon = lon1)

1. print("Moon " )

moonpos1<-getMoonPosition(date = as.POSIXct(datetimes2, tz = "UTC"), lat = lat1, lon = lon1)

moonillum1<-getMoonIllumination(date = c(datetimes2), keep = c("fraction", "phase"))

1. 1. phase 0.25 asc half moon crescent 0.5 full moon 0.75 last qarter

sunset01<-getSunlightTimes(date = as.Date(date1, tz = "UTC"), lat = lat1, lon = lon1, keep=c("sunset")) sunset1=as.POSIXct( sunset01['sunset']$sunset, tz = "EET")

sunrise01<-getSunlightTimes(date = as.Date(date2, tz = "UTC"), lat = lat1, lon = lon1, keep=c("sunrise")) sunrise1=as.POSIXct( sunrise01['sunrise']$sunrise, tz = "EET")

print(sunset1) print(sunrise1)

1. stop(-1)

salt1= sunpos1[4]$altitude saz1= sunpos1[4]$azimuth1

1. print(salt1)

salt2=180*salt1/3.1416

saz2=(360*salt1)/3.1416

malt1= moonpos1[4]$altitude maz1= moonpos1[4]$azimuth1

1. print(malt1)

malt2=180*malt1/3.1416

maz2=(360*malt1)/3.1416+360

print("Sun")

1. print(sunpos1)

1. print("Moon")

1. print(moonpos1)

date1=as.POSIXct(datetimes2, tz=timezone1)

intensity1<-calculateMoonlightIntensity(lat1, lon1, datetimes2, e1)

1. calculateMoonlightStatistics(lat, lon, date, e, t, timezone)

1. print(intensity1)

1. print(paikka1)
2. print(c(datetimes2, " UTC") )

1. print("Auringon korkeus, suuntakulma astetta")

1. print(salt2)
2. print(saz2)

1. print("Kuun korkeus, suuntakulma astetta")
2. print(malt2)

1. print(maz2)

1. print("Kuun valaistu osa ")
2. print( moonillum1[2]$fraction )
3. print("Kuun vaihe ")
4. print( moonillum1[3]$phase )

1. print ( "Valaistus lx" )
2. print(intensity1['illumination']$illumination)
3. print (" Josta kuunvaloa lx")
4. print(intensity1['moonlightModel']$moonlightModel)

svg("tulilahti_Finland_illumination_night_27_28_07_1979_1.svg")

1. plot(datetimes2,log10(intensity1['illumination']$illumination),
2. col='orange', type='l', lwd=5,
3. main="Illumination at Tulilahti, Heinävesi, Finland 27.-28.7. 1959",
4. col.main='blue',
5. xlab="Time UTC, 27-28.7 1959",
6. ylab="Log10 (illumination lx)")

print("KKKK")

1. datee <- hour(hm( datetimes2))
2. datee <-hour(datetimes2)

datee0<-datetimes2 datee1<-format(datee0, tz="EET",usetz=TRUE) datee<-as.POSIXct( datee1, tz = "EET")

log10lx<-log10(intensity1['illumination']$illumination)

1. print("pppp")

data2 <- data.frame(datee,log10lx)

1. print("qqq")

names(data2)<-c("Time_EET", "Illum_log_Lx")

print("rrr")

1. stop(-1)

1. Plot

ggplot(data2, aes(x=Time_EET, y=Illum_log_Lx)) +

 geom_line( color="#ff7f00", linewidth=1.5, alpha=0.9, linetype=1) +

ggtitle(label=" Illumination of Sun, night 27.-28.7.1959", subtitle="Tulilahti, Heinävesi, Savo, Finland ")+
 xlab("Time EET") + ylab("Max illumination, log10 of Lx")+

 scale_x_datetime(date_breaks = "1 hours", date_labels = "%H:%M")+
 theme_bw() +
 geom_hline(yintercept=log10(300), linetype="dashed", color = "red")+
 geom_hline(yintercept=log10(3.5), linetype="dashed", color = "red")+
 
  geom_hline(yintercept=log10(0.1), linetype="dashed", color = "blue") +
                                             
  annotate(geom="text", x=sunset1, y=log10(300), label="Sunset", color="black", size=6)+
  annotate(geom="text", x=sunrise1, y=log10(300), label="Sunrise", color="black", size=6)+
  annotate(geom="text", x=as.POSIXct("1959-7-27 22:00", tz = "EET"), y=log10(0.1), label="Full moon if it was", color="black", size=6)+
  annotate(geom="text", x=as.POSIXct("1959-7-27 22:00", tz = "EET"), y=log10(3.5), label="Limit of civil/nautical twilight ", color="black", size=6)+
  annotate(geom="text", x=as.POSIXct("1959-7-28 00:00", tz = "EET"), y=log10(30), label="Civil twilight ", color="blue", size=8)+
  annotate(geom="text", x=as.POSIXct("1959-7-28 00:00", tz = "EET"), y=log10(1), label="Nau. twi.", color="#00007f", size=8)+

  theme(text = element_text(size = 19, family="Sans") )+
theme(axis.text = element_text(size = 16)) +
theme(axis.title = element_text(size = 17)) +
theme(plot.title = element_text(size = 15, face = "bold"))
                                                   

dev.off()

1. statistics1<-calculateMoonlightStatistics(lat1, lon1, datetime1, e1, t1, timezone1)

1. print(statistics1)