Generate time line graphs in processData-ascii

Port timeline graph generation from processData-graph.

analyze/processData-ascii.py

@@ -1,19 +1,26 @@
 import sys
-
+import time
+import graph
+import Image,ImageDraw
 
 def run(filename):
 
 	f = open(filename)
 	counts = 0;
 	minutes = 0;
+	now = 0;
 	largest = 0;
 	smallest = 0xFFFFFF
-
 	bins = {}
+	start = 0;
+	cpm = []
+	maxcpm = 0;
+	lastcount = 0
+	minuteData = []
 
 	print("scanning for largest ADC value")	
 	while True:
-		line = f.readline()
+		line = f.readline().strip()
 		if (len(line) < 2):
 			break
 		if not line.startswith("time"):
@@ -32,9 +39,16 @@ def run(filename):
 			break
 		if line.startswith("time"):
 			minutes = minutes + 1
+			diff = counts - lastcount
+			cpm.append(minuteData)
+			if len(minuteData) > maxcpm:
+				maxcpm = len(minuteData)
+			minuteData = []
+			lastcount = counts
 			print str(minutes) + " min"
 			continue
 		data = int(line)
+		minuteData.append(data)
 		binNum = int(4096 * (data - smallest) / float(datarange))
 		if bins.has_key(binNum):
 			bins[binNum] = bins[binNum] + 1
@@ -47,15 +61,29 @@ def run(filename):
 
 	outputFilename = filename.split(".")[0] + "_bins.txt";
 	fbin = open(outputFilename,'w')
+	print("Duration: %d minutes\n" % minutes);
+	if minutes > 0:
+		print("Counts: %d (%f CPM)\n" % (counts, counts/minutes));
+	else:
+		print("Counts: %d (%f CPM)\n" % (counts, counts));
+
 	keys = bins.keys()
 	keys.sort()
 	for key in keys:
 		fbin.write(str(key) + " " + str(bins[key]) + "\n")
 	fbin.close()
 
+	g = graph.Graph(1000,1000,0,minutes,0,maxcpm)
+	for t in range(minutes):
+		#g.graphBar(t,cpm[t])
+		#print len(cpm[t])
+		g.graphColorBar(t,len(cpm[t]),[float(f)/float(0xFFFFF) for f in cpm[t]])
+	g.drawLabels("time (min) since %s [%d,%d]" % (time.strftime("%a, %d %b %Y %H:%M:%S +0000", time.gmtime(start)),0,minutes), "pulses recorded [%d,%d]" % (0,maxcpm))
+	g.save(filename.split(".")[0]+"_timeline.bmp")
 
 	if (minutes > 0):
-		print "avg cpm: " + str(counts/minutes)
+		print "avg cpm: %f" % (counts/minutes)
+		print "max cpm: %d" % (maxcpm)
 
 	return outputFilename