## To compute the root mean squared contrast of a fits image:

import pyfits
import numpy as np

def RMS(image):
   data = pyfits.getdata(image,ignore_missing_end=True)
   rms = data.std()/data.mean()
   return rms