example script convert to XISF median modify FITS header