5 Author: Bryan Steinbach
10 azs
,dazs
,model_dazs
,els
,dels
,model_dels
= np
.load('starcamdata-20100828.npy')
14 allowed_params
= 'npae ca ia ie tf ta te'
20 # Matplotlib doesn't have intelligent layout, so we have to manually pull in the subplots to make room for title and legend
21 pylab
.subplots_adjust(bottom
=0.16,top
=0.89,wspace
=0.25,hspace
=0.25)
23 pylab
.title('Star camera pointing offsets')
26 pylab
.scatter(azs
,dazs
,facecolor
='white')
27 pylab
.scatter(azs
,model_dazs
,marker
='+')
28 # pylab.xlabel('Az (degrees)',fontsize=axfs)
29 pylab
.ylabel('$\Delta$Az (arcminutes)',fontsize
=axfs
)
33 pylab
.scatter(els
,dazs
,facecolor
='white')
34 pylab
.scatter(els
,model_dazs
,marker
='+')
35 # pylab.xlabel('El (degrees)',fontsize=axfs)
36 # pylab.ylabel('$\Delta$Az (arcminutes)',fontsize=axfs)
40 pylab
.scatter(azs
,dels
,facecolor
='white')
41 pylab
.scatter(azs
,model_dels
,marker
='+')
42 pylab
.xlabel('Az (degrees)',fontsize
=axfs
)
43 pylab
.ylabel('$\Delta$El (arcminutes)',fontsize
=axfs
)
47 l1
= pylab
.scatter(els
,dels
,facecolor
='white')
48 l2
= pylab
.scatter(els
,model_dels
,marker
='+')
49 pylab
.xlabel('El (degrees)',fontsize
=axfs
)
50 # pylab.ylabel('$\Delta$El (arcminutes)',fontsize=axfs)
53 pylab
.suptitle('Star Camera Pointing Model 8/28/10 %d stars %.1f hours\nModel free parameters: %s'%(nsource
,total_time
,allowed_params
),fontsize
=14)
55 pylab
.figlegend((l1
,l2
),('Measured star position - boresight','Model'),'lower center',prop
={'size':10})
58 pylab
.savefig('starcam-20100828.png',figsize
=((10,8)),dpi
=250)