A short script to add RADESYS keyword to Pan-STARRS1 image cutouts obtained from the MAST archive. This is to help some packages/softwares (like ds9 v8.0) read the WCS info properly. For more information about the issue, see here.
The keyword is added using Astropy's FITS file handling capabilities.
To run it:
$ python ps1_wcs_corr.py ps1_cutout.fits