Cleanup
Change to astropy instead of pyfits and pywcs
Refactor to astropy and add pesta