# Clip values above a specified threshold

import numpy
import pyfits
import os

# Changes to the directory where the script is located
abspath = os.path.abspath(__file__)
dname = os.path.dirname(abspath)
os.chdir(dname)

infile = 'M33Clouds.fits'

outfile = 'M33Clouds_clipped.fits'

FitsFile = pyfits.open(infile)

image= FitsFile[0].data
header=FitsFile[0].header


image[image < 0.1] = 0.0

FitsFile = pyfits.writeto(outfile,image,header=header)