From 96708b2693eb6199085d40fa194c016436feb3c4 Mon Sep 17 00:00:00 2001
From: Bryan Steinbach <bsteinbach@berkeley.edu>
Date: Sun, 3 Oct 2010 14:50:49 -0700
Subject: [PATCH] starcam-plot clean up

---
 week4/problem1/starcam-plot.py | 40 ++++++++++++++++++++++++++--------------
 1 file changed, 26 insertions(+), 14 deletions(-)

diff --git a/week4/problem1/starcam-plot.py b/week4/problem1/starcam-plot.py
index cdb9d43..9daf2b9 100755
--- a/week4/problem1/starcam-plot.py
+++ b/week4/problem1/starcam-plot.py
@@ -1,50 +1,62 @@
 #!/usr/bin/env python
-
+'''
+AY250 HW4 Problem 1
+Make a nice plot
+Author:  Bryan Steinbach
+'''
 import numpy as np
 import pylab
 
 azs,dazs,model_dazs,els,dels,model_dels = np.load('starcamdata-20100828.npy')
 
-axfs = 16
+nsource = 15
+total_time = 6.1
+allowed_params = 'npae ca ia ie tf ta te'
+
+axfs = 10
 def azelplot():
 	pylab.figure()
 
+	# Matplotlib doesn't have intelligent layout, so we have to manually pull in the subplots to make room for title and legend
+	pylab.subplots_adjust(bottom=0.16,top=0.89,wspace=0.25,hspace=0.25)
+
 	pylab.title('Star camera pointing offsets')
 
 	pylab.subplot(221)
-	pylab.scatter(azs,dazs)
+	pylab.scatter(azs,dazs,facecolor='white')
 	pylab.scatter(azs,model_dazs,marker='+')
-	pylab.xlabel('Az (degrees)',fontsize=axfs)
+#	pylab.xlabel('Az (degrees)',fontsize=axfs)
 	pylab.ylabel('$\Delta$Az (arcminutes)',fontsize=axfs)
 	pylab.grid()
 
 	pylab.subplot(222)
-	pylab.scatter(els,dazs)
+	pylab.scatter(els,dazs,facecolor='white')
 	pylab.scatter(els,model_dazs,marker='+')
-	pylab.xlabel('El (degrees)',fontsize=axfs)
-	pylab.ylabel('$\Delta$Az (arcminutes)',fontsize=axfs)
+#	pylab.xlabel('El (degrees)',fontsize=axfs)
+#	pylab.ylabel('$\Delta$Az (arcminutes)',fontsize=axfs)
 	pylab.grid()
 
 	pylab.subplot(223)
-	pylab.scatter(azs,dels)
+	pylab.scatter(azs,dels,facecolor='white')
 	pylab.scatter(azs,model_dels,marker='+')
 	pylab.xlabel('Az (degrees)',fontsize=axfs)
 	pylab.ylabel('$\Delta$El (arcminutes)',fontsize=axfs)
 	pylab.grid()
 
 	pylab.subplot(224)
-	l1 = pylab.scatter(els,dels)
+	l1 = pylab.scatter(els,dels,facecolor='white')
 	l2 = pylab.scatter(els,model_dels,marker='+')
 	pylab.xlabel('El (degrees)',fontsize=axfs)
-	pylab.ylabel('$\Delta$El (arcminutes)',fontsize=axfs)
+#	pylab.ylabel('$\Delta$El (arcminutes)',fontsize=axfs)
 	pylab.grid()
 
-	pylab.suptitle('Star Camera Pointing Model  8/28/10  %d stars  %.1f hours\nModel free parameters: %s'%(nsource,total_time,allowed_params),fontsize=20)
+	pylab.suptitle('Star Camera Pointing Model  8/28/10  %d stars  %.1f hours\nModel free parameters: %s'%(nsource,total_time,allowed_params),fontsize=14)
+
+	pylab.figlegend((l1,l2),('Measured star position - boresight','Model'),'lower center',prop={'size':10})
 
-	pylab.figlegend((l1,l2),('Measured star position - boresight','Model'),'lower center',prop={'size':16})
 
-	pylab.savefig('starcam-20100828.png')
+	pylab.savefig('starcam-20100828.png',figsize=((10,8)),dpi=250)
 
 
 azelplot()
-pylab.show()
+#pylab.show()
-- 
2.11.4.GIT