WCS FITS conversion to Virtual GDAL header
Temporary location for script to go from a GDAL supported image to GeoFITS
From your python prompt:
>>> import fits2vrt
>>> a=fits2vrt.fitskeys('myimage.fits')
>>> a.fits2vrt()
will create myimage.vrt
To have fits2vrt
working you need to install