imcombine and pixel masks

Paul Kerry wrote on Dec 04, 2012

Hi

Our of our users has asked me to report this problem...
--
I'm attempting to median stack HAWK-I images in MEF files. I'm using imcombine and have pixel masks which I wish to use during the stacking.

A master pixel mask has been produced using optmasks, and the mask for each image is produced by shifting this mask spatially using python/pyfits.
The resulting masks seem to be OK and can, for example, be used with fixpix without errors or problems.

When using them in imcombine, however, I get the following error:

Dec 4 12:36: IMCOMBINE
combine = median, scale = none, zero = none, weight = none
reject = pclip, nkeep = 1
pclip = -1., lsigma = 3., hsigma = 3.
blank = 0.
masktype = goodval, maskval = 0
Images Maskfile
products/hawki_step_basic_calib_obj008.fits[extname=CHIP4.INT1] products/hawki_step_detect_obj008_mask.fits[extname=CHIP4.INT1]
products/hawki_step_basic_calib_obj007.fits[extname=CHIP4.INT1] products/hawki_step_detect_obj007_mask.fits[extname=CHIP4.INT1]
products/hawki_step_basic_calib_obj001.fits[extname=CHIP4.INT1] products/hawki_step_detect_obj001_mask.fits[extname=CHIP4.INT1]
products/hawki_step_basic_calib_obj003.fits[extname=CHIP4.INT1] products/hawki_step_detect_obj003_mask.fits[extname=CHIP4.INT1]
products/hawki_step_basic_calib_obj002.fits[extname=CHIP4.INT1] products/hawki_step_detect_obj002_mask.fits[extname=CHIP4.INT1]
products/hawki_step_basic_calib_obj006.fits[extname=CHIP4.INT1] products/hawki_step_detect_obj006_mask.fits[extname=CHIP4.INT1]

Output image = products/hawki_step_compute_bkg_005.fits[append], ncombine = 6
TRAP: 925 PLIO: reference out of bounds on mask
Warning: PLIO: reference out of bounds on mask

This is with IRAF v2.16 on both OS X and linux systems.

Does anyone have any idea why I get this error?

I've uploaded the images I'm using, along with the results of 'lpar imcombine' and the python script I am using to the following location

ftp://hydra.shef.ac.uk/pub/sl/imcombine

I am invoking the python script with

python bkgSubtract.py '*calib*' 3
--

Best regards
Paul.

---

Paul Kerry wrote on Dec 04, 2012

Please ignore the post - the dimensions of the image that were used to make the mask didn't match the images he was trying to mask.

Cheers
Paul.

---

Last post on Dec 04, 2012