Implementation seems correct: at least, objectest.py yields
the correct centroids when I create an image from random shapelets

BUT - reconstructed image flux out by a factor of 1.0/(xpixscale*ypixscale)

I think problem lies in how I'm implementing add_to_image in shapelet.py

May just be a case of normalising flux by pixel size
Code now seems to successfully do decomposition in 'real' scale
rather than pixel scale

Gallery of shapelet method not quite perfect

Merge with master once this is perfected?
Modified coefficients.py to allow I/O to files

Wrote simple script which analyses multiple ASCII files

Merge this into master
Conversion to real units successful - has fixed centroiding issues

Rewritten shapelet.py and coefficients.py to reuse code more often

Written some simple scripts to plot images and coefficients, as well as
a basic script to decompose multiple images

Now need scripts to analyse individual coefficients over many images
plot_image.py prints out centroid when plotting
Slightly modified the output nomenclature from decomposition runs
This allows the new plotting scripts to more easily produce graphs
of how individual shapelet coefficients evolve with time

	new file:   plot_allcoefficients_vs_time.py
	new file:   plot_coefficient_vs_time.py
Small mods to objecttest.py to test different additions of shapelets
to an image for decomposition