[SciPy-user] FITS images with header-supplied axes?

Is there any plotting routine in scipy / matplotlib that can plot a fits image with correct WCS coordinates on the axes? I know pyfits can load fits files, astLib has routines to interpret header coordinates, and I think you can make the axes different using matplotlib transforms, but is there anything that puts all three together currently available? Thanks, Adam