from psfsubtraction import data
psfbase, image, image2 = data.gaussian_PSFs()

from psfsubtraction.fitpsf import fitters
my_fitter = fitters.SimpleSubtraction(psfbase, image)
residual = my_fitter.remove_psf()

import matplotlib.pylab as plt
fig = plt.figure(figsize=(8, 4))
ax1 = fig.add_subplot(121)
im1 = ax1.imshow(image, interpolation='nearest', cmap=plt.cm.hot)
ax1.set_title('image')
plt.colorbar(im1, ax=ax1)
ax2 = fig.add_subplot(122)
im2 = ax2.imshow(residual, interpolation='nearest', cmap=plt.cm.hot)
plt.colorbar(im2, ax=ax2)
ax2.set_title('Subtracted image')