#!/usr/bin/Rscript
# Billboard Battery diagram

#### Settings ####
main.hour.interval = 3
main.anz.intervals = 16
battery.file.name = "/home/user/.billboardscripts/batterydata.csv"
#battery.file.name = "batterydata.csv"
plot.file.name = "/home/user/.billboardscripts/battery-diagram.bmp"
#plot.file.name = "battery-diagram.bmp"
plot.width = 400
plot.height = 150
plot.background = "black"
plot.foreground = "white"
plot.battery.linecolor  = "green"
plot.charge.color = "yellow"
plot.tempcolor = "red"
plot.voltagecolor = "cyan"
plot.day.mark.color = "red"

# Get time
now = Sys.time()
this.hour = strftime(now, format="%H")
this.date = strftime(now, format="%F")

#### Get next main hour ####
next.main.hour = as.numeric(this.hour) + 1 # Start with this hour

# Increase until next main hour is reached
while(next.main.hour %% main.hour.interval != 0) {
	next.main.hour = (next.main.hour + 1) %% 24 # Increase
	}

# Get time objects
next.main.time = strptime(x=paste(this.date,next.main.hour),format="%F %H")
if(next.main.hour==0) { next.main.time = next.main.time+86400 }


#### Generate axis sequence ####
axis.sequence = 
	seq(from=next.main.time,
			by=paste("-",main.hour.interval," hours",sep=""),
			length=main.anz.intervals)


#### Read battery data ####
# Read actual state from system
aktstate = as.numeric(system("qdbus --system org.freedesktop.Hal /org/freedesktop/Hal/devices/bme org.freedesktop.Hal.Device.GetProperty battery.charge_level.percentage",intern=TRUE))

# Read the current charging state, TRUE if currently charging, FALSE unless
aktcharge = system("bmestat|awk '/charging state/ { print $3}'",intern=TRUE) == "STARTED"

# Read the current battery temperature
akttemp = as.numeric(system("bmestat|awk '/temperature/ {print $3}'",intern=T))

# Read current voltage
aktvoltage = as.numeric(system("bmestat|awk '/battery cur. voltage/ {print $4}'",intern=T)) / 100


# Create dataframe
batterystate = data.frame(
	TIME=now,
	BATTERY=aktstate,
	CHARGING=aktcharge,
	TEMP=akttemp,
	VOLTAGE=aktvoltage
	)

# Read older entries
if(file.exists(battery.file.name)){
	# Read the data
	batterystateold = read.table(file=battery.file.name,sep="\t",header=TRUE)
	# Create columns that do not exist
	for (col in colnames(batterystate)) if(!col %in% colnames(batterystateold)) batterystateold[col] = NA
	for (col in colnames(batterystateold)) if(!col %in% colnames(batterystate)) batterystate[col] = NA
	# Concatenate data
	batterystate = rbind(batterystate,batterystateold)
	rm(batterystateold)
	}

# Sort the data
batterystate = batterystate[order(batterystate$TIME),]

# Writeout the new data
# This only writes the last period
#write.table(file=battery.file.name,x = batterystate[which(batterystate$TIME>(min(axis.sequence))-main.hour.interval*60*60),],row.names = F,quote = F,sep="\t")
write.table(file=battery.file.name,x = batterystate,row.names = F,quote = F,sep="\t") 

# Convert date
batterystate$TIME = strptime(x = batterystate$TIME, format="%F %H:%M:%S")

# Selection to plot
selection = batterystate$TIME>=min(axis.sequence)-60*60*main.hour.interval
TIME     = batterystate$TIME[selection]
BATTERY  = batterystate$BATTERY[selection]
CHARGING = batterystate$CHARGING[selection]
TEMP     = batterystate$TEMP[selection]
VOLTAGE  = batterystate$VOLTAGE[selection]

rm(selection) # Free memory
rm(batterystate) # Free memory

#### Plot ####
# Open a graphics device
bmp(file=plot.file.name,width=plot.width,height=plot.height)

# Set graphical parameters
par(bg=plot.background,col=plot.battery.linecolor,mar=c(2,2.4,0.3,0.1)+0.01)

# Plot basic
plot(x=axis.sequence,y=rep(1,length(axis.sequence)), type="n",
		 axes=F, ylim=c(0,100))

# Draw axes
axis(side=2, # Y-Axis
		 las=2, # Write text 90 degrees
		 col=plot.foreground,col.axis=plot.foreground, # Color
		 at=seq(from=0,to=100,by=25) # Labels
		 )

# Draw box
box(col=plot.foreground)

# Draw time axis
axis.POSIXct(side=1, # X-Axis
					at=axis.sequence, # Sequence
					format="%H", # Format
					col=plot.foreground, col.axis=plot.foreground, # Color
					padj=-0.5
					)

# Draw a grid
# Horizontal lines
abline(h=seq(from=0,to=100,by=25),col="darkgrey",lty=2)
# Vertical lines
abline(v=axis.sequence,col="darkgrey",lty=2)

# Mark days
abline(v=axis.sequence[which(as.numeric(strftime(axis.sequence,format="%H"))==0)],
			 col=plot.day.mark.color,
			 lty=2)

# Draw the voltage curve
Y = VOLTAGE
lines(
	x=TIME,
	Y,
	col=plot.voltagecolor,
	lwd=4,
	pch=20,
	cex=0.1,
	)


# Draw the temperature curve
Y = TEMP
lines(
	x=TIME,
	Y,
	col=plot.tempcolor,
	lwd=4,
	pch=20,
	cex=0.1,
	)


# Draw the normal curve
Y = BATTERY
lines(
	x=TIME,
	Y,
	col=plot.battery.linecolor,
	lwd=4,
	pch=20,
	cex=0.1,
	)

#	Draw the charging curve on top
Y = BATTERY
Y[CHARGING!=TRUE|is.na(CHARGING)] = NA
lines(
	x=TIME,
	Y,
	col=plot.charge.color,
	lwd=4,
	pch=20,
	cex=0.1,
	)

# Close graphics device
dev.off()

#!/bin/sh
while true; do
/home/user/.billboardscripts/Billboard-battery.R
sleep 300
done

R> source('NokiaN9-battery-period.R')
Period: 24.02245 hours

# How long is battery life on the Nokia N9?
# Find the period of charge of battery

# Read CSV file
BAT <- read.csv("NokiaN9-batterydata-5months.csv", sep = "\t")

# Interpolate on regular grid
APPROX <- approx(BAT$TIME, BAT$BATTERY,xout = seq(f = min(BAT$TIME), t = max(BAT$TIME), len = length(BAT$TIME)))

# Get spectrum
SPECTRUM <- spectrum(APPROX$y,plot=F)

# Find maximum period
# First find frequency of greatest amplitude
# Invert this frequency to get the period
# Multiply this period that is in units of timesteps with the timestep to get the period in seconds
PERIOD <- 1/SPECTRUM$freq[which.max(SPECTRUM$spec)] * unique(zapsmall(diff(APPROX$x)))

# Output period in hours
cat("Period:", PERIOD / 60 / 60,"hours")