Optimizing the Image Sampling#
Introduction#
This example was written to help users better understand the subtleties in improving image sampling for dithered data. One of the powers of the drizzling algorithm is that, given properly dithered images, it can restore much of the information lost due to undersampled images (Fruchter and Hook, 2002).
This work is based on ISR ACS 2015-01, which contains a more detailed discussion than presented here.
In practice, this requires the use of Astrodrizzle
task within the Drizzlepac package. This example will cover the following topics:
1. Description of the Data
2. Run astrodrizzle several times using different settings for the 'final_pixfrac' and 'final_scale' parameters
3. Compare and evaluate results of using different 'final_pixfrac' and 'final_scale' values
import glob
import shutil
import os
from astropy.io import fits
from astroquery.mast import Observations
from astropy.visualization import LogStretch,ImageNormalize,LinearStretch
from astropy import wcs
import drizzlepac
import matplotlib.pyplot as plt
import numpy as np
%matplotlib inline
The following task in the stsci.skypac package can be run with TEAL:
skymatch
The following tasks in the drizzlepac package can be run with TEAL:
astrodrizzle config_testbed imagefindpars mapreg
photeq pixreplace pixtopix pixtosky
refimagefindpars resetbits runastrodriz skytopix
tweakback tweakreg updatenpol
1. Description of the data#
Four images of the spiral galaxy NGC 3370 (Program 11570) were acquired using the WFC3/IR F160W imaging filter. The WFC3-IR-DITHERBOX-MIN
dither pattern designed to provide optimal sampling of the PSF was used.
Data are downloaded using the astroquery
API to access the MAST archive. The astroquery.mast
documentation has more examples for how to find and download data from MAST. For this example, the caibrated FLT exposures from visit 19 will be downloaded.
# Querying MAST and downloading files.
obs_ids = ['ib1f19010']
obs_table = Observations.query_criteria(obs_id=obs_ids, obstype='all')
download_tab = Observations.download_products(obs_table['obsid'], mrp_only=False, productSubGroupDescription=['FLT'])
Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:HST/product/ib1f19l6q_flt.fits to ./mastDownload/HST/ib1f19l6q/ib1f19l6q_flt.fits ...
[Done]
Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:HST/product/ib1f19l7q_flt.fits to ./mastDownload/HST/ib1f19l7q/ib1f19l7q_flt.fits ...
[Done]
Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:HST/product/ib1f19l9q_flt.fits to ./mastDownload/HST/ib1f19l9q/ib1f19l9q_flt.fits ...
[Done]
Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:HST/product/ib1f19laq_flt.fits to ./mastDownload/HST/ib1f19laq/ib1f19laq_flt.fits ...
[Done]
Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:HST/product/hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits to ./mastDownload/HST/hst_11570_19_wfc3_ir_f160w_ib1f19la/hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits ...
[Done]
Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:HST/product/hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits to ./mastDownload/HST/hst_11570_19_wfc3_ir_f160w_ib1f19l6/hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits ...
[Done]
Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:HST/product/hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits to ./mastDownload/HST/hst_11570_19_wfc3_ir_f160w_ib1f19l7/hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits ...
[Done]
Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:HST/product/hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits to ./mastDownload/HST/hst_11570_19_wfc3_ir_f160w_ib1f19l9/hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits ...
[Done]
# Cleaning up directories after dowloading from MAST
for file in download_tab['Local Path']:
os.rename(file, os.path.basename(file))
shutil.rmtree('mastDownload')
2. Dithering#
For the drizzle algorithm to work optimally, observations need to be dithered so that the PSF is optimally sampled. The code below creates plots to show the dithering of each of the images in the association. The plot on the left shows how the images were dithered on the sky (the POSTARG). The plot on the right shows how the offsets translate to pixel phase (e.g. subpixel position).
flt_files = glob.glob('*flt.fits')
plate_scale = fits.getval(flt_files[0], 'idcscale',ext=1)
print('Plate scale: {:>6.4f}'.format(plate_scale))
postarg1 = np.empty(len(flt_files), dtype=np.float)
postarg2 = np.empty(len(flt_files), dtype=np.float)
x_phase = np.empty(len(flt_files), dtype=np.float)
y_phase = np.empty(len(flt_files), dtype=np.float)
for i, im in enumerate(flt_files):
with fits.open(im) as hdu:
postarg1[i] = hdu[0].header['postarg1']
postarg2[i] = hdu[0].header['postarg2']
x_phase[i] = abs(np.modf(postarg1[i] / plate_scale)[0])
y_phase[i] = abs(np.modf(postarg2[i] / plate_scale)[0])
fig, ax = plt.subplots(1, 2, figsize=(12, 6))
ax[0].plot(postarg1, postarg2, 'k+-', ms=15, lw=1)
ax[0].set_xlabel('arcsec', fontsize='large')
ax[0].set_ylabel('arcsec', fontsize='large')
ax[1].plot(x_phase, y_phase, 'k+-', ms=15, lw=1)
ax[1].set_xlim(-0.05, 1)
ax[1].set_ylim(-0.05, 1)
ax[1].set_xlabel('X pixel phase', fontsize='large')
ax[1].set_ylabel('Y pixel phase', fontsize='large')
Plate scale: 0.1283
Text(0, 0.5, 'Y pixel phase')
3. Drizzling#
Before drizzling can be performed, a plate scale for the output image should be chosen. This will usually be dictated by the science needs. In theory, critical sampling of a PSF occurs with 2.355 pixels. The FWHM of the WFC3/IR detector at 1600 nm is ~0.151” (Table 7.5 of WFC Instrument Handbook). Using these values, a plate scale of 0.065”/pix will be selected for the drizzled output frame, approximately half of the native plate scale (0.1283”/pix).
The mechanism by which the drizzle algorithm improves sampling in the output image is by shrinking the input pixels before adding them to the output pixel grid (see Figure 2 in Fruchter and Hook (2002)). In practice this process is controlled in the AstroDrizzle
package by adjusting the final_pixfrac
parameter, which is the fractional size of the pixel to be used.
Below AstroDrizzle
is called twice. The first time, no changes are made to the plate scale or the size of the drop. In the second call, the plate scale and pixfrac are reduced. For a four point dither stack, a final_pixfrac
value of 0.8 is usually optimal.
drizzlepac.astrodrizzle.AstroDrizzle('*flt.fits',
output='f160w_noopt',
runfile='',
context=False,
build=True,
preserve=False,
clean=True,
driz_sep_bits='64,512',
final_bits='64,512',
final_wcs=True,
final_rot=0.)
drizzlepac.astrodrizzle.AstroDrizzle('*flt.fits',
output='f160w_opt',
runfile='',
context=False,
build=True,
preserve=False,
clean=True,
driz_sep_bits='64,512',
final_bits='64,512',
final_pixfrac=0.8,
final_wcs=True,
final_rot=0.,
final_scale=0.065)
No trailer file created...
AstroDrizzle log file: final.tra
AstroDrizzle Version 3.5.1 started at: 16:52:31.634 (19/04/2024)
==== Processing Step Initialization started at 16:52:31.635 (19/04/2024)
##############################################################################
# #
# “minmed” is highly recommended for three images, #
# and is good for four to six images, #
# but should be avoided for ten or more images. #
# #
##############################################################################
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 591.9026595870587 535.72854962001
CD1_1 CD1_2 : -3.562500079472858e-05 1.7530017893904253e-22
CD2_1 CD2_2 : -1.241796017022117e-25 3.562500079472858e-05
NAXIS : 1184 1072
********************************************************************************
*
* Estimated memory usage: up to 89 Mb.
* Output image size: 1184 X 1072 pixels.
* Output image file: ~ 14 Mb.
* Cores available: 4
*
********************************************************************************
==== Processing Step Initialization finished at 16:52:32.664 (19/04/2024)
==== Processing Step Static Mask started at 16:52:32.665 (19/04/2024)
==== Processing Step Static Mask finished at 16:52:32.787 (19/04/2024)
==== Processing Step Subtract Sky started at 16:52:32.787 (19/04/2024)
***** skymatch started on 2024-04-19 16:52:32.986962
Version 1.0.9
'skymatch' task will apply computed sky differences to input image file(s).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
'MDRIZSKY' header keyword will represent sky value *computed* from data.
----- User specified keywords: -----
Sky Value Keyword: 'MDRIZSKY'
Data Units Keyword: 'BUNIT'
----- Input file list: -----
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l6_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l7_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l9_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19la_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l6q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l6q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l7q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l7q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l9q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l9q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19laq_flt.fits'
EXT: 'SCI',1; MASK: ib1f19laq_skymatch_mask_sci1.fits[0]
----- Sky statistics parameters: -----
statistics function: 'median'
lower = None
upper = None
nclip = 5
lsigma = 4.0
usigma = 4.0
binwidth = 0.1
----- Data->Brightness conversion parameters for input files: -----
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l6q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l7q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l9q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19laq_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
----- Computing sky values requested image extensions (detector chips): -----
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
* Image: 'ib1f19l6q_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'ib1f19l7q_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'ib1f19l9q_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'ib1f19laq_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
***** skymatch ended on 2024-04-19 16:52:33.318744
TOTAL RUN TIME: 0:00:00.331782
==== Processing Step Subtract Sky finished at 16:52:33.419 (19/04/2024)
==== Processing Step Separate Drizzle started at 16:52:33.419 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 544.5 481.5
CD1_1 CD1_2 : -3.542984305540133e-05 3.723829036629122e-06
CD2_1 CD2_2 : 3.723829036629122e-06 3.542984305540133e-05
NAXIS : 1089 963
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Writing out image to disk:
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_wht.fits
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_wht.fits
-Generating simple FITS output: ib1f19l6q_single_sci.fits
-Generating simple FITS output: ib1f19l7q_single_sci.fits
-Generating simple FITS output: ib1f19laq_single_sci.fits
-Generating simple FITS output: ib1f19l9q_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l7q_single_sci.fits
Writing out image to disk:
Deleted all instances of WCS with key B in extensions [0]
ib1f19l6q_single_sci.fits
Writing out image to disk:
ib1f19laq_single_sci.fits
Writing out image to disk:
ib1f19l9q_single_sci.fits
Writing out image to disk:
ib1f19l7q_single_wht.fits
Writing out image to disk:
ib1f19l6q_single_wht.fits
Writing out image to disk:
ib1f19laq_single_wht.fits
Writing out image to disk:
ib1f19l9q_single_wht.fits
==== Processing Step Separate Drizzle finished at 16:52:35.480 (19/04/2024)
==== Processing Step Create Median started at 16:52:35.481 (19/04/2024)
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits' is 464.93669456453324
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits' is 437.5836524801294
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits' is 427.0147483274913
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits' is 424.6136825690036
reference sky value for image 'ib1f19l6q_flt.fits' is 464.93669456453324
reference sky value for image 'ib1f19l7q_flt.fits' is 437.5836524801294
reference sky value for image 'ib1f19l9q_flt.fits' is 427.0147483274913
reference sky value for image 'ib1f19laq_flt.fits' is 424.6136825690036
Saving output median image to: 'f160w_noopt_med.fits'
==== Processing Step Create Median finished at 16:52:35.924 (19/04/2024)
==== Processing Step Blot started at 16:52:35.92 (19/04/2024)
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Blot: creating blotted image: ib1f19l6q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l6q_sci1_blt.fits
Writing out image to disk: ib1f19l6q_sci1_blt.fits
Blot: creating blotted image: ib1f19l7q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l7q_sci1_blt.fits
Writing out image to disk: ib1f19l7q_sci1_blt.fits
Blot: creating blotted image: ib1f19l9q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l9q_sci1_blt.fits
Writing out image to disk: ib1f19l9q_sci1_blt.fits
Blot: creating blotted image: ib1f19laq_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19laq_sci1_blt.fits
Writing out image to disk: ib1f19laq_sci1_blt.fits
==== Processing Step Blot finished at 16:52:37.463 (19/04/2024)
==== Processing Step Driz_CR started at 16:52:37.463 (19/04/2024)
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_crmask.fits
Creating output: ib1f19l6q_sci1_crmask.fits
Creating output: ib1f19l7q_sci1_crmask.fits
Creating output: ib1f19l9q_sci1_crmask.fits
Creating output: ib1f19laq_sci1_crmask.fits
==== Processing Step Driz_CR finished at 16:52:39.794 (19/04/2024)
==== Processing Step Final Drizzle started at 16:52:39.803 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 591.9026595870587 535.72854962001
CD1_1 CD1_2 : -3.562500079472858e-05 1.7530017893904253e-22
CD2_1 CD2_2 : -1.241796017022117e-25 3.562500079472858e-05
NAXIS : 1184 1072
-Generating multi-extension output file: f160w_noopt_drz.fits
Deleted all instances of WCS with key A in extensions [1]
Deleted all instances of WCS with key B in extensions [1]
Deleted all instances of WCS with key C in extensions [1]
Writing out to disk: f160w_noopt_drz.fits
==== Processing Step Final Drizzle finished at 16:52:42.669 (19/04/2024)
AstroDrizzle Version 3.5.1 is finished processing at 16:52:42.670 (19/04/2024).
-------------------- --------------------
Step Elapsed time
-------------------- --------------------
Initialization 1.0297 sec.
Static Mask 0.1220 sec.
Subtract Sky 0.6316 sec.
Separate Drizzle 2.0608 sec.
Create Median 0.4430 sec.
Blot 1.5381 sec.
Driz_CR 2.3310 sec.
Final Drizzle 2.8663 sec.
==================== ====================
Total 11.0226 sec.
No trailer file saved...
No trailer file created...
AstroDrizzle log file: final.tra
AstroDrizzle Version 3.5.1 started at: 16:52:42.714 (19/04/2024)
==== Processing Step Initialization started at 16:52:42.714 (19/04/2024)
##############################################################################
# #
# “minmed” is highly recommended for three images, #
# and is good for four to six images, #
# but should be avoided for ten or more images. #
# #
##############################################################################
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
********************************************************************************
*
* Estimated memory usage: up to 257 Mb.
* Output image size: 2336 X 2115 pixels.
* Output image file: ~ 56 Mb.
* Cores available: 4
*
********************************************************************************
==== Processing Step Initialization finished at 16:52:43.775 (19/04/2024)
==== Processing Step Static Mask started at 16:52:43.776 (19/04/2024)
==== Processing Step Static Mask finished at 16:52:43.894 (19/04/2024)
==== Processing Step Subtract Sky started at 16:52:43.89 (19/04/2024)
***** skymatch started on 2024-04-19 16:52:43.995602
Version 1.0.9
'skymatch' task will apply computed sky differences to input image file(s).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
'MDRIZSKY' header keyword will represent sky value *computed* from data.
----- User specified keywords: -----
Sky Value Keyword: 'MDRIZSKY'
Data Units Keyword: 'BUNIT'
----- Input file list: -----
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l6_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l7_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l9_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19la_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l6q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l6q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l7q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l7q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l9q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l9q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19laq_flt.fits'
EXT: 'SCI',1; MASK: ib1f19laq_skymatch_mask_sci1.fits[0]
----- Sky statistics parameters: -----
statistics function: 'median'
lower = None
upper = None
nclip = 5
lsigma = 4.0
usigma = 4.0
binwidth = 0.1
----- Data->Brightness conversion parameters for input files: -----
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l6q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l7q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l9q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19laq_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
----- Computing sky values requested image extensions (detector chips): -----
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
* Image: 'ib1f19l6q_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'ib1f19l7q_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'ib1f19l9q_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'ib1f19laq_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
***** skymatch ended on 2024-04-19 16:52:44.303264
TOTAL RUN TIME: 0:00:00.307662
==== Processing Step Subtract Sky finished at 16:52:44.402 (19/04/2024)
==== Processing Step Separate Drizzle started at 16:52:44.403 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 544.5 481.5
CD1_1 CD1_2 : -3.542984305540133e-05 3.723829036629122e-06
CD2_1 CD2_2 : 3.723829036629122e-06 3.542984305540133e-05
NAXIS : 1089 963
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_wht.fits
Writing out image to disk:
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_wht.fits
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_wht.fits
-Generating simple FITS output: ib1f19l6q_single_sci.fits
-Generating simple FITS output: ib1f19l7q_single_sci.fits
-Generating simple FITS output: ib1f19l9q_single_sci.fits
-Generating simple FITS output: ib1f19laq_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
Deleted all instances of WCS with key B in extensions [0]
ib1f19l6q_single_sci.fits
Writing out image to disk:
ib1f19l7q_single_sci.fits
Writing out image to disk:
Writing out image to disk:
ib1f19l9q_single_sci.fits
ib1f19laq_single_sci.fits
Writing out image to disk:
ib1f19l6q_single_wht.fits
Writing out image to disk:
ib1f19l7q_single_wht.fits
Writing out image to disk:
Writing out image to disk:
ib1f19l9q_single_wht.fits
ib1f19laq_single_wht.fits
==== Processing Step Separate Drizzle finished at 16:52:46.454 (19/04/2024)
==== Processing Step Create Median started at 16:52:46.455 (19/04/2024)
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits' is 464.93669456453324
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits' is 437.5836524801294
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits' is 427.0147483274913
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits' is 424.6136825690036
reference sky value for image 'ib1f19l6q_flt.fits' is 464.93669456453324
reference sky value for image 'ib1f19l7q_flt.fits' is 437.5836524801294
reference sky value for image 'ib1f19l9q_flt.fits' is 427.0147483274913
reference sky value for image 'ib1f19laq_flt.fits' is 424.6136825690036
Saving output median image to: 'f160w_opt_med.fits'
==== Processing Step Create Median finished at 16:52:46.887 (19/04/2024)
==== Processing Step Blot started at 16:52:46.888 (19/04/2024)
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Blot: creating blotted image: ib1f19l6q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l6q_sci1_blt.fits
Writing out image to disk: ib1f19l6q_sci1_blt.fits
Blot: creating blotted image: ib1f19l7q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l7q_sci1_blt.fits
Writing out image to disk: ib1f19l7q_sci1_blt.fits
Blot: creating blotted image: ib1f19l9q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l9q_sci1_blt.fits
Writing out image to disk: ib1f19l9q_sci1_blt.fits
Blot: creating blotted image: ib1f19laq_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19laq_sci1_blt.fits
Writing out image to disk: ib1f19laq_sci1_blt.fits
==== Processing Step Blot finished at 16:52:48.534 (19/04/2024)
==== Processing Step Driz_CR started at 16:52:48.5 (19/04/2024)
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_crmask.fits
Creating output: ib1f19l6q_sci1_crmask.fits
Creating output: ib1f19l7q_sci1_crmask.fits
Creating output: ib1f19laq_sci1_crmask.fits
Creating output: ib1f19l9q_sci1_crmask.fits
==== Processing Step Driz_CR finished at 16:52:50.770 (19/04/2024)
==== Processing Step Final Drizzle started at 16:52:50.779 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
-Generating multi-extension output file: f160w_opt_drz.fits
Deleted all instances of WCS with key A in extensions [1]
Deleted all instances of WCS with key B in extensions [1]
Deleted all instances of WCS with key C in extensions [1]
Writing out to disk: f160w_opt_drz.fits
==== Processing Step Final Drizzle finished at 16:52:54.404 (19/04/2024)
AstroDrizzle Version 3.5.1 is finished processing at 16:52:54.405 (19/04/2024).
-------------------- --------------------
Step Elapsed time
-------------------- --------------------
Initialization 1.0612 sec.
Static Mask 0.1176 sec.
Subtract Sky 0.5077 sec.
Separate Drizzle 2.0515 sec.
Create Median 0.4323 sec.
Blot 1.6461 sec.
Driz_CR 2.2360 sec.
Final Drizzle 3.6254 sec.
==================== ====================
Total 11.6777 sec.
No trailer file saved...
4. Results#
The drizzled science and weight images produced from the first call to AstroDrizzle with no optimization of the plate scale and pixfrac are plotted below.
with fits.open('f160w_noopt_drz.fits') as hdu:
im1wcs = wcs.WCS(hdu[1].header)
sci1 = hdu[1].data
wht1 = hdu[2].data
norm1 = ImageNormalize(sci1, vmin=-0.2, vmax=25, stretch=LogStretch())
fig, ax = plt.subplots(1, 2, figsize=(16, 8), subplot_kw={'projection':im1wcs})
ax[0].imshow(sci1, norm=norm1, cmap='gray', origin='lower')
ax[1].imshow(wht1, cmap='gray', origin='lower')
<matplotlib.image.AxesImage at 0x7f4f1a2dd950>
The drizzled science image is on the left and the associated weight image is on the right, both without optimization of the plate scale and pixfrac.
To compare, the figure plotted below shows close ups of the same part of the sky from the two drizzled products.
radeclims = wcs.utils.pixel_to_skycoord([825, 875], [930, 980], im1wcs)
with fits.open('f160w_opt_drz.fits') as hdu:
im2wcs = wcs.WCS(hdu[1].header)
sci2 = hdu[1].data
wht2 = hdu[2].data
norm2 = ImageNormalize(sci2, vmin=-0.05, vmax=25, stretch=LogStretch())
fig = plt.figure(figsize=(16, 8))
ax = fig.add_subplot(1, 2, 1, projection=im1wcs)
ax.imshow(sci1, norm=norm2, cmap='gray')
ax.set_xlim(radeclims.to_pixel(im1wcs)[0])
ax.set_ylim(radeclims.to_pixel(im1wcs)[1])
ax = fig.add_subplot(1, 2, 2, projection=im2wcs)
ax.imshow(sci2, norm=norm2, cmap='gray')
ax.set_xlim(radeclims.to_pixel(im2wcs)[0])
ax.set_ylim(radeclims.to_pixel(im2wcs)[1])
(1835.448117879787, 1934.1019662353683)
The image on the left is without improved plate scale and pixfrac, and shows that the detector undersamples the PSF. The image on the right is the image with improved parameters where the resolution is greatly improved.
5. Optimizing the final_pixfrac
parameter#
While the optimized final_pixfrac
in the example above was chosen from experience with using a four-point dither, the value that should be used is not known a priori. The value could be different depending on several factors. For example, if the number of images is greater than four, the value used for final_pixfrac
could be smaller since more images are available to fill in holes in the output grid. On the other hand, it is possible that no dithering was used during the observations. In that case, final_pixfrac
should be left at 1.0, since shrinking the size of the drop could be detrimental.
Below, a series of experiments will be run to determine the best final_pixfrac
value for the selected output plate scale, by varying pixfrac in steps of 0.1 over a range of values from 0.1 to 1.0.
pixfracs = np.arange(0.1, 1.1, 0.1)
for pixfrac in pixfracs:
outname = 'f160w_{:.1f}'.format(pixfrac)
drizzlepac.astrodrizzle.AstroDrizzle('*flt.fits',
output=outname,
runfile='',
context=False,
build=True,
preserve=False,
clean=True,
driz_sep_bits='64,512',
final_bits='64,512',
final_pixfrac=pixfrac,
final_wcs=True,
final_rot=0.,
final_scale=0.065)
No trailer file created...
AstroDrizzle log file: final.tra
AstroDrizzle Version 3.5.1 started at: 16:52:56.415 (19/04/2024)
==== Processing Step Initialization started at 16:52:56.415 (19/04/2024)
##############################################################################
# #
# “minmed” is highly recommended for three images, #
# and is good for four to six images, #
# but should be avoided for ten or more images. #
# #
##############################################################################
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
********************************************************************************
*
* Estimated memory usage: up to 257 Mb.
* Output image size: 2336 X 2115 pixels.
* Output image file: ~ 56 Mb.
* Cores available: 4
*
********************************************************************************
==== Processing Step Initialization finished at 16:52:57.553 (19/04/2024)
==== Processing Step Static Mask started at 16:52:57.554 (19/04/2024)
==== Processing Step Static Mask finished at 16:52:57.670 (19/04/2024)
==== Processing Step Subtract Sky started at 16:52:57.671 (19/04/2024)
***** skymatch started on 2024-04-19 16:52:57.771980
Version 1.0.9
'skymatch' task will apply computed sky differences to input image file(s).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
'MDRIZSKY' header keyword will represent sky value *computed* from data.
----- User specified keywords: -----
Sky Value Keyword: 'MDRIZSKY'
Data Units Keyword: 'BUNIT'
----- Input file list: -----
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l6_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l7_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l9_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19la_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l6q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l6q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l7q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l7q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l9q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l9q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19laq_flt.fits'
EXT: 'SCI',1; MASK: ib1f19laq_skymatch_mask_sci1.fits[0]
----- Sky statistics parameters: -----
statistics function: 'median'
lower = None
upper = None
nclip = 5
lsigma = 4.0
usigma = 4.0
binwidth = 0.1
----- Data->Brightness conversion parameters for input files: -----
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l6q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l7q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l9q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19laq_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
----- Computing sky values requested image extensions (detector chips): -----
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
* Image: 'ib1f19l6q_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'ib1f19l7q_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'ib1f19l9q_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'ib1f19laq_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
***** skymatch ended on 2024-04-19 16:52:58.080736
TOTAL RUN TIME: 0:00:00.308756
==== Processing Step Subtract Sky finished at 16:52:58.181 (19/04/2024)
==== Processing Step Separate Drizzle started at 16:52:58.182 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 544.5 481.5
CD1_1 CD1_2 : -3.542984305540133e-05 3.723829036629122e-06
CD2_1 CD2_2 : 3.723829036629122e-06 3.542984305540133e-05
NAXIS : 1089 963
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_wht.fits
-Generating simple FITS output: ib1f19l6q_single_sci.fits
-Generating simple FITS output: ib1f19l7q_single_sci.fits
-Generating simple FITS output: ib1f19l9q_single_sci.fits
-Generating simple FITS output: ib1f19laq_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l6q_single_sci.fits
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
Writing out image to disk:
ib1f19l9q_single_sci.fits
ib1f19l7q_single_sci.fits
Writing out image to disk:
Writing out image to disk:
ib1f19laq_single_sci.fits
ib1f19l6q_single_wht.fits
Writing out image to disk:
ib1f19l9q_single_wht.fits
Writing out image to disk:
ib1f19l7q_single_wht.fits
Writing out image to disk:
ib1f19laq_single_wht.fits
==== Processing Step Separate Drizzle finished at 16:53:00.239 (19/04/2024)
==== Processing Step Create Median started at 16:53:00.240 (19/04/2024)
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits' is 464.93669456453324
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits' is 437.5836524801294
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits' is 427.0147483274913
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits' is 424.6136825690036
reference sky value for image 'ib1f19l6q_flt.fits' is 464.93669456453324
reference sky value for image 'ib1f19l7q_flt.fits' is 437.5836524801294
reference sky value for image 'ib1f19l9q_flt.fits' is 427.0147483274913
reference sky value for image 'ib1f19laq_flt.fits' is 424.6136825690036
Saving output median image to: 'f160w_0.1_med.fits'
==== Processing Step Create Median finished at 16:53:00.668 (19/04/2024)
==== Processing Step Blot started at 16:53:00.668 (19/04/2024)
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Blot: creating blotted image: ib1f19l6q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l6q_sci1_blt.fits
Writing out image to disk: ib1f19l6q_sci1_blt.fits
Blot: creating blotted image: ib1f19l7q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l7q_sci1_blt.fits
Writing out image to disk: ib1f19l7q_sci1_blt.fits
Blot: creating blotted image: ib1f19l9q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l9q_sci1_blt.fits
Writing out image to disk: ib1f19l9q_sci1_blt.fits
Blot: creating blotted image: ib1f19laq_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19laq_sci1_blt.fits
Writing out image to disk: ib1f19laq_sci1_blt.fits
==== Processing Step Blot finished at 16:53:02.203 (19/04/2024)
==== Processing Step Driz_CR started at 16:53:02.203 (19/04/2024)
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_crmask.fits
Creating output: ib1f19l6q_sci1_crmask.fits
Creating output: ib1f19l7q_sci1_crmask.fits
Creating output: ib1f19l9q_sci1_crmask.fits
Creating output: ib1f19laq_sci1_crmask.fits
==== Processing Step Driz_CR finished at 16:53:04.449 (19/04/2024)
==== Processing Step Final Drizzle started at 16:53:04.461 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
-Generating multi-extension output file: f160w_0.1_drz.fits
Deleted all instances of WCS with key A in extensions [1]
Deleted all instances of WCS with key B in extensions [1]
Deleted all instances of WCS with key C in extensions [1]
Writing out to disk: f160w_0.1_drz.fits
==== Processing Step Final Drizzle finished at 16:53:06.935 (19/04/2024)
AstroDrizzle Version 3.5.1 is finished processing at 16:53:06.935 (19/04/2024).
-------------------- --------------------
Step Elapsed time
-------------------- --------------------
Initialization 1.1384 sec.
Static Mask 0.1161 sec.
Subtract Sky 0.5108 sec.
Separate Drizzle 2.0567 sec.
Create Median 0.4282 sec.
Blot 1.5344 sec.
Driz_CR 2.2454 sec.
Final Drizzle 2.4739 sec.
==================== ====================
Total 10.5039 sec.
No trailer file saved...
No trailer file created...
AstroDrizzle log file: final.tra
AstroDrizzle Version 3.5.1 started at: 16:53:06.983 (19/04/2024)
==== Processing Step Initialization started at 16:53:06.983 (19/04/2024)
##############################################################################
# #
# “minmed” is highly recommended for three images, #
# and is good for four to six images, #
# but should be avoided for ten or more images. #
# #
##############################################################################
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
********************************************************************************
*
* Estimated memory usage: up to 257 Mb.
* Output image size: 2336 X 2115 pixels.
* Output image file: ~ 56 Mb.
* Cores available: 4
*
********************************************************************************
==== Processing Step Initialization finished at 16:53:08.139 (19/04/2024)
==== Processing Step Static Mask started at 16:53:08.140 (19/04/2024)
==== Processing Step Static Mask finished at 16:53:08.261 (19/04/2024)
==== Processing Step Subtract Sky started at 16:53:08.262 (19/04/2024)
***** skymatch started on 2024-04-19 16:53:08.365315
Version 1.0.9
'skymatch' task will apply computed sky differences to input image file(s).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
'MDRIZSKY' header keyword will represent sky value *computed* from data.
----- User specified keywords: -----
Sky Value Keyword: 'MDRIZSKY'
Data Units Keyword: 'BUNIT'
----- Input file list: -----
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l6_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l7_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l9_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19la_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l6q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l6q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l7q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l7q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l9q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l9q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19laq_flt.fits'
EXT: 'SCI',1; MASK: ib1f19laq_skymatch_mask_sci1.fits[0]
----- Sky statistics parameters: -----
statistics function: 'median'
lower = None
upper = None
nclip = 5
lsigma = 4.0
usigma = 4.0
binwidth = 0.1
----- Data->Brightness conversion parameters for input files: -----
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l6q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l7q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l9q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19laq_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
----- Computing sky values requested image extensions (detector chips): -----
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
* Image: 'ib1f19l6q_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'ib1f19l7q_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'ib1f19l9q_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'ib1f19laq_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
***** skymatch ended on 2024-04-19 16:53:08.677293
TOTAL RUN TIME: 0:00:00.311978
==== Processing Step Subtract Sky finished at 16:53:08.78 (19/04/2024)
==== Processing Step Separate Drizzle started at 16:53:08.781 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 544.5 481.5
CD1_1 CD1_2 : -3.542984305540133e-05 3.723829036629122e-06
CD2_1 CD2_2 : 3.723829036629122e-06 3.542984305540133e-05
NAXIS : 1089 963
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_wht.fits
-Generating simple FITS output: ib1f19l6q_single_sci.fits
-Generating simple FITS output: ib1f19l7q_single_sci.fits
-Generating simple FITS output: ib1f19l9q_single_sci.fits
-Generating simple FITS output: ib1f19laq_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l6q_single_sci.fits
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l7q_single_sci.fits
Writing out image to disk:
ib1f19l9q_single_sci.fits
Writing out image to disk:
ib1f19laq_single_sci.fits
Writing out image to disk:
ib1f19l6q_single_wht.fits
Writing out image to disk:
ib1f19l7q_single_wht.fits
Writing out image to disk:
ib1f19l9q_single_wht.fits
Writing out image to disk:
ib1f19laq_single_wht.fits
==== Processing Step Separate Drizzle finished at 16:53:10.83 (19/04/2024)
==== Processing Step Create Median started at 16:53:10.840 (19/04/2024)
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits' is 464.93669456453324
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits' is 437.5836524801294
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits' is 427.0147483274913
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits' is 424.6136825690036
reference sky value for image 'ib1f19l6q_flt.fits' is 464.93669456453324
reference sky value for image 'ib1f19l7q_flt.fits' is 437.5836524801294
reference sky value for image 'ib1f19l9q_flt.fits' is 427.0147483274913
reference sky value for image 'ib1f19laq_flt.fits' is 424.6136825690036
Saving output median image to: 'f160w_0.2_med.fits'
==== Processing Step Create Median finished at 16:53:11.274 (19/04/2024)
==== Processing Step Blot started at 16:53:11.27 (19/04/2024)
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Blot: creating blotted image: ib1f19l6q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l6q_sci1_blt.fits
Writing out image to disk: ib1f19l6q_sci1_blt.fits
Blot: creating blotted image: ib1f19l7q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l7q_sci1_blt.fits
Writing out image to disk: ib1f19l7q_sci1_blt.fits
Blot: creating blotted image: ib1f19l9q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l9q_sci1_blt.fits
Writing out image to disk: ib1f19l9q_sci1_blt.fits
Blot: creating blotted image: ib1f19laq_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19laq_sci1_blt.fits
Writing out image to disk: ib1f19laq_sci1_blt.fits
==== Processing Step Blot finished at 16:53:12.81 (19/04/2024)
==== Processing Step Driz_CR started at 16:53:12.814 (19/04/2024)
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_crmask.fits
Creating output: ib1f19l6q_sci1_crmask.fits
Creating output: ib1f19l7q_sci1_crmask.fits
Creating output: ib1f19l9q_sci1_crmask.fits
Creating output: ib1f19laq_sci1_crmask.fits
==== Processing Step Driz_CR finished at 16:53:15.066 (19/04/2024)
==== Processing Step Final Drizzle started at 16:53:15.079 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
-Generating multi-extension output file: f160w_0.2_drz.fits
Deleted all instances of WCS with key A in extensions [1]
Deleted all instances of WCS with key B in extensions [1]
Deleted all instances of WCS with key C in extensions [1]
Writing out to disk: f160w_0.2_drz.fits
==== Processing Step Final Drizzle finished at 16:53:17.726 (19/04/2024)
AstroDrizzle Version 3.5.1 is finished processing at 16:53:17.727 (19/04/2024).
-------------------- --------------------
Step Elapsed time
-------------------- --------------------
Initialization 1.1565 sec.
Static Mask 0.1214 sec.
Subtract Sky 0.5189 sec.
Separate Drizzle 2.0582 sec.
Create Median 0.4339 sec.
Blot 1.5394 sec.
Driz_CR 2.2514 sec.
Final Drizzle 2.6474 sec.
==================== ====================
Total 10.7271 sec.
No trailer file saved...
No trailer file created...
AstroDrizzle log file: final.tra
AstroDrizzle Version 3.5.1 started at: 16:53:17.775 (19/04/2024)
==== Processing Step Initialization started at 16:53:17.775 (19/04/2024)
##############################################################################
# #
# “minmed” is highly recommended for three images, #
# and is good for four to six images, #
# but should be avoided for ten or more images. #
# #
##############################################################################
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
********************************************************************************
*
* Estimated memory usage: up to 257 Mb.
* Output image size: 2336 X 2115 pixels.
* Output image file: ~ 56 Mb.
* Cores available: 4
*
********************************************************************************
==== Processing Step Initialization finished at 16:53:18.831 (19/04/2024)
==== Processing Step Static Mask started at 16:53:18.832 (19/04/2024)
==== Processing Step Static Mask finished at 16:53:18.951 (19/04/2024)
==== Processing Step Subtract Sky started at 16:53:18.951 (19/04/2024)
***** skymatch started on 2024-04-19 16:53:19.163246
Version 1.0.9
'skymatch' task will apply computed sky differences to input image file(s).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
'MDRIZSKY' header keyword will represent sky value *computed* from data.
----- User specified keywords: -----
Sky Value Keyword: 'MDRIZSKY'
Data Units Keyword: 'BUNIT'
----- Input file list: -----
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l6_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l7_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l9_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19la_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l6q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l6q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l7q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l7q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l9q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l9q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19laq_flt.fits'
EXT: 'SCI',1; MASK: ib1f19laq_skymatch_mask_sci1.fits[0]
----- Sky statistics parameters: -----
statistics function: 'median'
lower = None
upper = None
nclip = 5
lsigma = 4.0
usigma = 4.0
binwidth = 0.1
----- Data->Brightness conversion parameters for input files: -----
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l6q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l7q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l9q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19laq_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
----- Computing sky values requested image extensions (detector chips): -----
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
* Image: 'ib1f19l6q_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'ib1f19l7q_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'ib1f19l9q_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'ib1f19laq_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
***** skymatch ended on 2024-04-19 16:53:19.471472
TOTAL RUN TIME: 0:00:00.308226
==== Processing Step Subtract Sky finished at 16:53:19.57 (19/04/2024)
==== Processing Step Separate Drizzle started at 16:53:19.574 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 544.5 481.5
CD1_1 CD1_2 : -3.542984305540133e-05 3.723829036629122e-06
CD2_1 CD2_2 : 3.723829036629122e-06 3.542984305540133e-05
NAXIS : 1089 963
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Writing out image to disk:
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_wht.fits
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_wht.fits
-Generating simple FITS output: ib1f19l6q_single_sci.fits
-Generating simple FITS output: ib1f19l7q_single_sci.fits
-Generating simple FITS output: ib1f19l9q_single_sci.fits
-Generating simple FITS output: ib1f19laq_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l6q_single_sci.fits
Writing out image to disk:
ib1f19l7q_single_sci.fits
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l9q_single_sci.fits
Writing out image to disk:
ib1f19laq_single_sci.fits
Writing out image to disk:
ib1f19l6q_single_wht.fits
Writing out image to disk:
ib1f19l7q_single_wht.fits
Writing out image to disk:
ib1f19l9q_single_wht.fits
Writing out image to disk:
ib1f19laq_single_wht.fits
==== Processing Step Separate Drizzle finished at 16:53:21.631 (19/04/2024)
==== Processing Step Create Median started at 16:53:21.631 (19/04/2024)
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits' is 464.93669456453324
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits' is 437.5836524801294
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits' is 427.0147483274913
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits' is 424.6136825690036
reference sky value for image 'ib1f19l6q_flt.fits' is 464.93669456453324
reference sky value for image 'ib1f19l7q_flt.fits' is 437.5836524801294
reference sky value for image 'ib1f19l9q_flt.fits' is 427.0147483274913
reference sky value for image 'ib1f19laq_flt.fits' is 424.6136825690036
Saving output median image to: 'f160w_0.3_med.fits'
==== Processing Step Create Median finished at 16:53:22.062 (19/04/2024)
==== Processing Step Blot started at 16:53:22.063 (19/04/2024)
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Blot: creating blotted image: ib1f19l6q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l6q_sci1_blt.fits
Writing out image to disk: ib1f19l6q_sci1_blt.fits
Blot: creating blotted image: ib1f19l7q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l7q_sci1_blt.fits
Writing out image to disk: ib1f19l7q_sci1_blt.fits
Blot: creating blotted image: ib1f19l9q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l9q_sci1_blt.fits
Writing out image to disk: ib1f19l9q_sci1_blt.fits
Blot: creating blotted image: ib1f19laq_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19laq_sci1_blt.fits
Writing out image to disk: ib1f19laq_sci1_blt.fits
==== Processing Step Blot finished at 16:53:23.607 (19/04/2024)
==== Processing Step Driz_CR started at 16:53:23.608 (19/04/2024)
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_crmask.fits
Creating output: ib1f19l6q_sci1_crmask.fits
Creating output: ib1f19l7q_sci1_crmask.fits
Creating output: ib1f19l9q_sci1_crmask.fits
Creating output: ib1f19laq_sci1_crmask.fits
==== Processing Step Driz_CR finished at 16:53:25.866 (19/04/2024)
==== Processing Step Final Drizzle started at 16:53:25.878 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
-Generating multi-extension output file: f160w_0.3_drz.fits
Deleted all instances of WCS with key A in extensions [1]
Deleted all instances of WCS with key B in extensions [1]
Deleted all instances of WCS with key C in extensions [1]
Writing out to disk: f160w_0.3_drz.fits
==== Processing Step Final Drizzle finished at 16:53:28.648 (19/04/2024)
AstroDrizzle Version 3.5.1 is finished processing at 16:53:28.649 (19/04/2024).
-------------------- --------------------
Step Elapsed time
-------------------- --------------------
Initialization 1.0562 sec.
Static Mask 0.1194 sec.
Subtract Sky 0.6217 sec.
Separate Drizzle 2.0566 sec.
Create Median 0.4311 sec.
Blot 1.5443 sec.
Driz_CR 2.2578 sec.
Final Drizzle 2.7699 sec.
==================== ====================
Total 10.8571 sec.
No trailer file saved...
No trailer file created...
AstroDrizzle log file: final.tra
AstroDrizzle Version 3.5.1 started at: 16:53:28.695 (19/04/2024)
==== Processing Step Initialization started at 16:53:28.695 (19/04/2024)
##############################################################################
# #
# “minmed” is highly recommended for three images, #
# and is good for four to six images, #
# but should be avoided for ten or more images. #
# #
##############################################################################
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
********************************************************************************
*
* Estimated memory usage: up to 257 Mb.
* Output image size: 2336 X 2115 pixels.
* Output image file: ~ 56 Mb.
* Cores available: 4
*
********************************************************************************
==== Processing Step Initialization finished at 16:53:29.750 (19/04/2024)
==== Processing Step Static Mask started at 16:53:29.750 (19/04/2024)
==== Processing Step Static Mask finished at 16:53:29.870 (19/04/2024)
==== Processing Step Subtract Sky started at 16:53:29.871 (19/04/2024)
***** skymatch started on 2024-04-19 16:53:29.973655
Version 1.0.9
'skymatch' task will apply computed sky differences to input image file(s).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
'MDRIZSKY' header keyword will represent sky value *computed* from data.
----- User specified keywords: -----
Sky Value Keyword: 'MDRIZSKY'
Data Units Keyword: 'BUNIT'
----- Input file list: -----
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l6_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l7_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l9_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19la_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l6q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l6q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l7q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l7q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l9q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l9q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19laq_flt.fits'
EXT: 'SCI',1; MASK: ib1f19laq_skymatch_mask_sci1.fits[0]
----- Sky statistics parameters: -----
statistics function: 'median'
lower = None
upper = None
nclip = 5
lsigma = 4.0
usigma = 4.0
binwidth = 0.1
----- Data->Brightness conversion parameters for input files: -----
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l6q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l7q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l9q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19laq_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
----- Computing sky values requested image extensions (detector chips): -----
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
* Image: 'ib1f19l6q_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'ib1f19l7q_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'ib1f19l9q_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'ib1f19laq_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
***** skymatch ended on 2024-04-19 16:53:30.288516
TOTAL RUN TIME: 0:00:00.314861
==== Processing Step Subtract Sky finished at 16:53:30.391 (19/04/2024)
==== Processing Step Separate Drizzle started at 16:53:30.391 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 544.5 481.5
CD1_1 CD1_2 : -3.542984305540133e-05 3.723829036629122e-06
CD2_1 CD2_2 : 3.723829036629122e-06 3.542984305540133e-05
NAXIS : 1089 963
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
Writing out image to disk:
Deleted all instances of WCS with key C in extensions [0]
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
Writing out image to disk:
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_wht.fits
-Generating simple FITS output: ib1f19l6q_single_sci.fits
-Generating simple FITS output: ib1f19l7q_single_sci.fits
-Generating simple FITS output: ib1f19l9q_single_sci.fits
-Generating simple FITS output: ib1f19laq_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l7q_single_sci.fits
Writing out image to disk:
ib1f19l6q_single_sci.fits
Writing out image to disk:
ib1f19l9q_single_sci.fits
Writing out image to disk:
ib1f19laq_single_sci.fits
Writing out image to disk:
Writing out image to disk:
ib1f19l7q_single_wht.fits
ib1f19l6q_single_wht.fits
Writing out image to disk:
ib1f19l9q_single_wht.fits
Writing out image to disk:
ib1f19laq_single_wht.fits
==== Processing Step Separate Drizzle finished at 16:53:32.449 (19/04/2024)
==== Processing Step Create Median started at 16:53:32.450 (19/04/2024)
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits' is 464.93669456453324
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits' is 437.5836524801294
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits' is 427.0147483274913
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits' is 424.6136825690036
reference sky value for image 'ib1f19l6q_flt.fits' is 464.93669456453324
reference sky value for image 'ib1f19l7q_flt.fits' is 437.5836524801294
reference sky value for image 'ib1f19l9q_flt.fits' is 427.0147483274913
reference sky value for image 'ib1f19laq_flt.fits' is 424.6136825690036
Saving output median image to: 'f160w_0.4_med.fits'
==== Processing Step Create Median finished at 16:53:32.883 (19/04/2024)
==== Processing Step Blot started at 16:53:32.883 (19/04/2024)
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Blot: creating blotted image: ib1f19l6q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l6q_sci1_blt.fits
Writing out image to disk: ib1f19l6q_sci1_blt.fits
Blot: creating blotted image: ib1f19l7q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l7q_sci1_blt.fits
Writing out image to disk: ib1f19l7q_sci1_blt.fits
Blot: creating blotted image: ib1f19l9q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l9q_sci1_blt.fits
Writing out image to disk: ib1f19l9q_sci1_blt.fits
Blot: creating blotted image: ib1f19laq_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19laq_sci1_blt.fits
Writing out image to disk: ib1f19laq_sci1_blt.fits
==== Processing Step Blot finished at 16:53:34.541 (19/04/2024)
==== Processing Step Driz_CR started at 16:53:34.542 (19/04/2024)
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_crmask.fits
Creating output: ib1f19l6q_sci1_crmask.fits
Creating output: ib1f19l7q_sci1_crmask.fits
Creating output: ib1f19l9q_sci1_crmask.fits
Creating output: ib1f19laq_sci1_crmask.fits
==== Processing Step Driz_CR finished at 16:53:36.793 (19/04/2024)
==== Processing Step Final Drizzle started at 16:53:36.80 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
-Generating multi-extension output file: f160w_0.4_drz.fits
Deleted all instances of WCS with key A in extensions [1]
Deleted all instances of WCS with key B in extensions [1]
Deleted all instances of WCS with key C in extensions [1]
Writing out to disk: f160w_0.4_drz.fits
==== Processing Step Final Drizzle finished at 16:53:39.607 (19/04/2024)
AstroDrizzle Version 3.5.1 is finished processing at 16:53:39.608 (19/04/2024).
-------------------- --------------------
Step Elapsed time
-------------------- --------------------
Initialization 1.0549 sec.
Static Mask 0.1199 sec.
Subtract Sky 0.5202 sec.
Separate Drizzle 2.0573 sec.
Create Median 0.4329 sec.
Blot 1.6585 sec.
Driz_CR 2.2513 sec.
Final Drizzle 2.8009 sec.
==================== ====================
Total 10.8960 sec.
No trailer file saved...
No trailer file created...
AstroDrizzle log file: final.tra
AstroDrizzle Version 3.5.1 started at: 16:53:39.654 (19/04/2024)
==== Processing Step Initialization started at 16:53:39.654 (19/04/2024)
##############################################################################
# #
# “minmed” is highly recommended for three images, #
# and is good for four to six images, #
# but should be avoided for ten or more images. #
# #
##############################################################################
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
********************************************************************************
*
* Estimated memory usage: up to 257 Mb.
* Output image size: 2336 X 2115 pixels.
* Output image file: ~ 56 Mb.
* Cores available: 4
*
********************************************************************************
==== Processing Step Initialization finished at 16:53:40.810 (19/04/2024)
==== Processing Step Static Mask started at 16:53:40.810 (19/04/2024)
==== Processing Step Static Mask finished at 16:53:40.930 (19/04/2024)
==== Processing Step Subtract Sky started at 16:53:40.930 (19/04/2024)
***** skymatch started on 2024-04-19 16:53:41.035227
Version 1.0.9
'skymatch' task will apply computed sky differences to input image file(s).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
'MDRIZSKY' header keyword will represent sky value *computed* from data.
----- User specified keywords: -----
Sky Value Keyword: 'MDRIZSKY'
Data Units Keyword: 'BUNIT'
----- Input file list: -----
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l6_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l7_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l9_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19la_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l6q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l6q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l7q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l7q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l9q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l9q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19laq_flt.fits'
EXT: 'SCI',1; MASK: ib1f19laq_skymatch_mask_sci1.fits[0]
----- Sky statistics parameters: -----
statistics function: 'median'
lower = None
upper = None
nclip = 5
lsigma = 4.0
usigma = 4.0
binwidth = 0.1
----- Data->Brightness conversion parameters for input files: -----
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l6q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l7q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l9q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19laq_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
----- Computing sky values requested image extensions (detector chips): -----
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
* Image: 'ib1f19l6q_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'ib1f19l7q_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'ib1f19l9q_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'ib1f19laq_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
***** skymatch ended on 2024-04-19 16:53:41.347400
TOTAL RUN TIME: 0:00:00.312173
==== Processing Step Subtract Sky finished at 16:53:41.450 (19/04/2024)
==== Processing Step Separate Drizzle started at 16:53:41.451 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 544.5 481.5
CD1_1 CD1_2 : -3.542984305540133e-05 3.723829036629122e-06
CD2_1 CD2_2 : 3.723829036629122e-06 3.542984305540133e-05
NAXIS : 1089 963
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_wht.fits
Writing out image to disk:
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_wht.fits
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_wht.fits
-Generating simple FITS output: ib1f19l6q_single_sci.fits
-Generating simple FITS output: ib1f19l7q_single_sci.fits
-Generating simple FITS output: ib1f19l9q_single_sci.fits
-Generating simple FITS output: ib1f19laq_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l7q_single_sci.fits
Writing out image to disk:
Writing out image to disk:
Writing out image to disk:
ib1f19l6q_single_sci.fits
ib1f19l9q_single_sci.fits
ib1f19laq_single_sci.fits
Writing out image to disk:
ib1f19l7q_single_wht.fits
Writing out image to disk:
Writing out image to disk:
ib1f19l9q_single_wht.fits
ib1f19l6q_single_wht.fits
Writing out image to disk:
ib1f19laq_single_wht.fits
==== Processing Step Separate Drizzle finished at 16:53:43.506 (19/04/2024)
==== Processing Step Create Median started at 16:53:43.507 (19/04/2024)
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits' is 464.93669456453324
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits' is 437.5836524801294
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits' is 427.0147483274913
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits' is 424.6136825690036
reference sky value for image 'ib1f19l6q_flt.fits' is 464.93669456453324
reference sky value for image 'ib1f19l7q_flt.fits' is 437.5836524801294
reference sky value for image 'ib1f19l9q_flt.fits' is 427.0147483274913
reference sky value for image 'ib1f19laq_flt.fits' is 424.6136825690036
Saving output median image to: 'f160w_0.5_med.fits'
==== Processing Step Create Median finished at 16:53:43.937 (19/04/2024)
==== Processing Step Blot started at 16:53:43.937 (19/04/2024)
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Blot: creating blotted image: ib1f19l6q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l6q_sci1_blt.fits
Writing out image to disk: ib1f19l6q_sci1_blt.fits
Blot: creating blotted image: ib1f19l7q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l7q_sci1_blt.fits
Writing out image to disk: ib1f19l7q_sci1_blt.fits
Blot: creating blotted image: ib1f19l9q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l9q_sci1_blt.fits
Writing out image to disk: ib1f19l9q_sci1_blt.fits
Blot: creating blotted image: ib1f19laq_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19laq_sci1_blt.fits
Writing out image to disk: ib1f19laq_sci1_blt.fits
==== Processing Step Blot finished at 16:53:45.59 (19/04/2024)
==== Processing Step Driz_CR started at 16:53:45.592 (19/04/2024)
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_crmask.fits
Creating output: ib1f19l6q_sci1_crmask.fits
Creating output: ib1f19l7q_sci1_crmask.fits
Creating output: ib1f19l9q_sci1_crmask.fits
Creating output: ib1f19laq_sci1_crmask.fits
==== Processing Step Driz_CR finished at 16:53:47.864 (19/04/2024)
==== Processing Step Final Drizzle started at 16:53:47.877 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
-Generating multi-extension output file: f160w_0.5_drz.fits
Deleted all instances of WCS with key A in extensions [1]
Deleted all instances of WCS with key B in extensions [1]
Deleted all instances of WCS with key C in extensions [1]
Writing out to disk: f160w_0.5_drz.fits
==== Processing Step Final Drizzle finished at 16:53:50.819 (19/04/2024)
AstroDrizzle Version 3.5.1 is finished processing at 16:53:50.820 (19/04/2024).
-------------------- --------------------
Step Elapsed time
-------------------- --------------------
Initialization 1.1553 sec.
Static Mask 0.1196 sec.
Subtract Sky 0.5195 sec.
Separate Drizzle 2.0553 sec.
Create Median 0.4300 sec.
Blot 1.6543 sec.
Driz_CR 2.2724 sec.
Final Drizzle 2.9429 sec.
==================== ====================
Total 11.1493 sec.
No trailer file saved...
No trailer file created...
AstroDrizzle log file: final.tra
AstroDrizzle Version 3.5.1 started at: 16:53:50.869 (19/04/2024)
==== Processing Step Initialization started at 16:53:50.869 (19/04/2024)
##############################################################################
# #
# “minmed” is highly recommended for three images, #
# and is good for four to six images, #
# but should be avoided for ten or more images. #
# #
##############################################################################
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
********************************************************************************
*
* Estimated memory usage: up to 257 Mb.
* Output image size: 2336 X 2115 pixels.
* Output image file: ~ 56 Mb.
* Cores available: 4
*
********************************************************************************
==== Processing Step Initialization finished at 16:53:52.027 (19/04/2024)
==== Processing Step Static Mask started at 16:53:52.02 (19/04/2024)
==== Processing Step Static Mask finished at 16:53:52.148 (19/04/2024)
==== Processing Step Subtract Sky started at 16:53:52.149 (19/04/2024)
***** skymatch started on 2024-04-19 16:53:52.253140
Version 1.0.9
'skymatch' task will apply computed sky differences to input image file(s).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
'MDRIZSKY' header keyword will represent sky value *computed* from data.
----- User specified keywords: -----
Sky Value Keyword: 'MDRIZSKY'
Data Units Keyword: 'BUNIT'
----- Input file list: -----
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l6_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l7_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l9_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19la_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l6q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l6q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l7q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l7q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l9q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l9q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19laq_flt.fits'
EXT: 'SCI',1; MASK: ib1f19laq_skymatch_mask_sci1.fits[0]
----- Sky statistics parameters: -----
statistics function: 'median'
lower = None
upper = None
nclip = 5
lsigma = 4.0
usigma = 4.0
binwidth = 0.1
----- Data->Brightness conversion parameters for input files: -----
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l6q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l7q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l9q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19laq_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
----- Computing sky values requested image extensions (detector chips): -----
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
* Image: 'ib1f19l6q_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'ib1f19l7q_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'ib1f19l9q_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'ib1f19laq_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
***** skymatch ended on 2024-04-19 16:53:52.571528
TOTAL RUN TIME: 0:00:00.318388
==== Processing Step Subtract Sky finished at 16:53:52.676 (19/04/2024)
==== Processing Step Separate Drizzle started at 16:53:52.677 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 544.5 481.5
CD1_1 CD1_2 : -3.542984305540133e-05 3.723829036629122e-06
CD2_1 CD2_2 : 3.723829036629122e-06 3.542984305540133e-05
NAXIS : 1089 963
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
Writing out image to disk:
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_wht.fits
-Generating simple FITS output: ib1f19l6q_single_sci.fits
-Generating simple FITS output: ib1f19l7q_single_sci.fits
-Generating simple FITS output: ib1f19l9q_single_sci.fits
-Generating simple FITS output: ib1f19laq_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l6q_single_sci.fits
Writing out image to disk:
Writing out image to disk:
ib1f19l9q_single_sci.fits
ib1f19l7q_single_sci.fits
Writing out image to disk:
ib1f19laq_single_sci.fits
Writing out image to disk:
Writing out image to disk:
ib1f19l6q_single_wht.fits
ib1f19l9q_single_wht.fits
Writing out image to disk:
ib1f19l7q_single_wht.fits
Writing out image to disk:
ib1f19laq_single_wht.fits
==== Processing Step Separate Drizzle finished at 16:53:54.733 (19/04/2024)
==== Processing Step Create Median started at 16:53:54.734 (19/04/2024)
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits' is 464.93669456453324
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits' is 437.5836524801294
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits' is 427.0147483274913
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits' is 424.6136825690036
reference sky value for image 'ib1f19l6q_flt.fits' is 464.93669456453324
reference sky value for image 'ib1f19l7q_flt.fits' is 437.5836524801294
reference sky value for image 'ib1f19l9q_flt.fits' is 427.0147483274913
reference sky value for image 'ib1f19laq_flt.fits' is 424.6136825690036
Saving output median image to: 'f160w_0.6_med.fits'
==== Processing Step Create Median finished at 16:53:55.164 (19/04/2024)
==== Processing Step Blot started at 16:53:55.164 (19/04/2024)
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Blot: creating blotted image: ib1f19l6q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l6q_sci1_blt.fits
Writing out image to disk: ib1f19l6q_sci1_blt.fits
Blot: creating blotted image: ib1f19l7q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l7q_sci1_blt.fits
Writing out image to disk: ib1f19l7q_sci1_blt.fits
Blot: creating blotted image: ib1f19l9q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l9q_sci1_blt.fits
Writing out image to disk: ib1f19l9q_sci1_blt.fits
Blot: creating blotted image: ib1f19laq_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19laq_sci1_blt.fits
Writing out image to disk: ib1f19laq_sci1_blt.fits
==== Processing Step Blot finished at 16:53:56.707 (19/04/2024)
==== Processing Step Driz_CR started at 16:53:56.708 (19/04/2024)
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_crmask.fits
Creating output: ib1f19l6q_sci1_crmask.fits
Creating output: ib1f19l7q_sci1_crmask.fits
Creating output: ib1f19l9q_sci1_crmask.fits
Creating output: ib1f19laq_sci1_crmask.fits
==== Processing Step Driz_CR finished at 16:53:58.996 (19/04/2024)
==== Processing Step Final Drizzle started at 16:53:59.008 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
-Generating multi-extension output file: f160w_0.6_drz.fits
Deleted all instances of WCS with key A in extensions [1]
Deleted all instances of WCS with key B in extensions [1]
Deleted all instances of WCS with key C in extensions [1]
Writing out to disk: f160w_0.6_drz.fits
==== Processing Step Final Drizzle finished at 16:54:02.249 (19/04/2024)
AstroDrizzle Version 3.5.1 is finished processing at 16:54:02.249 (19/04/2024).
-------------------- --------------------
Step Elapsed time
-------------------- --------------------
Initialization 1.1580 sec.
Static Mask 0.1202 sec.
Subtract Sky 0.5276 sec.
Separate Drizzle 2.0561 sec.
Create Median 0.4300 sec.
Blot 1.5429 sec.
Driz_CR 2.2883 sec.
Final Drizzle 3.2404 sec.
==================== ====================
Total 11.3635 sec.
No trailer file saved...
No trailer file created...
AstroDrizzle log file: final.tra
AstroDrizzle Version 3.5.1 started at: 16:54:02.29 (19/04/2024)
==== Processing Step Initialization started at 16:54:02.297 (19/04/2024)
##############################################################################
# #
# “minmed” is highly recommended for three images, #
# and is good for four to six images, #
# but should be avoided for ten or more images. #
# #
##############################################################################
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
********************************************************************************
*
* Estimated memory usage: up to 257 Mb.
* Output image size: 2336 X 2115 pixels.
* Output image file: ~ 56 Mb.
* Cores available: 4
*
********************************************************************************
==== Processing Step Initialization finished at 16:54:03.362 (19/04/2024)
==== Processing Step Static Mask started at 16:54:03.363 (19/04/2024)
==== Processing Step Static Mask finished at 16:54:03.483 (19/04/2024)
==== Processing Step Subtract Sky started at 16:54:03.484 (19/04/2024)
***** skymatch started on 2024-04-19 16:54:03.696696
Version 1.0.9
'skymatch' task will apply computed sky differences to input image file(s).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
'MDRIZSKY' header keyword will represent sky value *computed* from data.
----- User specified keywords: -----
Sky Value Keyword: 'MDRIZSKY'
Data Units Keyword: 'BUNIT'
----- Input file list: -----
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l6_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l7_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l9_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19la_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l6q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l6q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l7q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l7q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l9q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l9q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19laq_flt.fits'
EXT: 'SCI',1; MASK: ib1f19laq_skymatch_mask_sci1.fits[0]
----- Sky statistics parameters: -----
statistics function: 'median'
lower = None
upper = None
nclip = 5
lsigma = 4.0
usigma = 4.0
binwidth = 0.1
----- Data->Brightness conversion parameters for input files: -----
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l6q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l7q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l9q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19laq_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
----- Computing sky values requested image extensions (detector chips): -----
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
* Image: 'ib1f19l6q_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'ib1f19l7q_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'ib1f19l9q_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'ib1f19laq_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
***** skymatch ended on 2024-04-19 16:54:04.029743
TOTAL RUN TIME: 0:00:00.333047
==== Processing Step Subtract Sky finished at 16:54:04.135 (19/04/2024)
==== Processing Step Separate Drizzle started at 16:54:04.135 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 544.5 481.5
CD1_1 CD1_2 : -3.542984305540133e-05 3.723829036629122e-06
CD2_1 CD2_2 : 3.723829036629122e-06 3.542984305540133e-05
NAXIS : 1089 963
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
Writing out image to disk:
Deleted all instances of WCS with key C in extensions [0]
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_wht.fits
-Generating simple FITS output: ib1f19l6q_single_sci.fits
-Generating simple FITS output: ib1f19l7q_single_sci.fits
-Generating simple FITS output: ib1f19l9q_single_sci.fits
-Generating simple FITS output: ib1f19laq_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l6q_single_sci.fits
Writing out image to disk:
Deleted all instances of WCS with key B in extensions [0]
ib1f19l9q_single_sci.fits
Writing out image to disk:
ib1f19l7q_single_sci.fits
Writing out image to disk:
Writing out image to disk:
ib1f19laq_single_sci.fits
Writing out image to disk:
ib1f19l6q_single_wht.fits
ib1f19l9q_single_wht.fits
Writing out image to disk:
ib1f19l7q_single_wht.fits
Writing out image to disk:
ib1f19laq_single_wht.fits
==== Processing Step Separate Drizzle finished at 16:54:06.195 (19/04/2024)
==== Processing Step Create Median started at 16:54:06.196 (19/04/2024)
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits' is 464.93669456453324
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits' is 437.5836524801294
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits' is 427.0147483274913
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits' is 424.6136825690036
reference sky value for image 'ib1f19l6q_flt.fits' is 464.93669456453324
reference sky value for image 'ib1f19l7q_flt.fits' is 437.5836524801294
reference sky value for image 'ib1f19l9q_flt.fits' is 427.0147483274913
reference sky value for image 'ib1f19laq_flt.fits' is 424.6136825690036
Saving output median image to: 'f160w_0.7_med.fits'
==== Processing Step Create Median finished at 16:54:06.624 (19/04/2024)
==== Processing Step Blot started at 16:54:06.625 (19/04/2024)
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Blot: creating blotted image: ib1f19l6q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l6q_sci1_blt.fits
Writing out image to disk: ib1f19l6q_sci1_blt.fits
Blot: creating blotted image: ib1f19l7q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l7q_sci1_blt.fits
Writing out image to disk: ib1f19l7q_sci1_blt.fits
Blot: creating blotted image: ib1f19l9q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l9q_sci1_blt.fits
Writing out image to disk: ib1f19l9q_sci1_blt.fits
Blot: creating blotted image: ib1f19laq_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19laq_sci1_blt.fits
Writing out image to disk: ib1f19laq_sci1_blt.fits
==== Processing Step Blot finished at 16:54:08.182 (19/04/2024)
==== Processing Step Driz_CR started at 16:54:08.182 (19/04/2024)
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_crmask.fits
Creating output: ib1f19l6q_sci1_crmask.fits
Creating output: ib1f19l7q_sci1_crmask.fits
Creating output: ib1f19l9q_sci1_crmask.fits
Creating output: ib1f19laq_sci1_crmask.fits
==== Processing Step Driz_CR finished at 16:54:10.42 (19/04/2024)
==== Processing Step Final Drizzle started at 16:54:10.440 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
-Generating multi-extension output file: f160w_0.7_drz.fits
Deleted all instances of WCS with key A in extensions [1]
Deleted all instances of WCS with key B in extensions [1]
Deleted all instances of WCS with key C in extensions [1]
Writing out to disk: f160w_0.7_drz.fits
==== Processing Step Final Drizzle finished at 16:54:13.929 (19/04/2024)
AstroDrizzle Version 3.5.1 is finished processing at 16:54:13.930 (19/04/2024).
-------------------- --------------------
Step Elapsed time
-------------------- --------------------
Initialization 1.0651 sec.
Static Mask 0.1205 sec.
Subtract Sky 0.6509 sec.
Separate Drizzle 2.0600 sec.
Create Median 0.4280 sec.
Blot 1.5571 sec.
Driz_CR 2.2445 sec.
Final Drizzle 3.4898 sec.
==================== ====================
Total 11.6159 sec.
No trailer file saved...
No trailer file created...
AstroDrizzle log file: final.tra
AstroDrizzle Version 3.5.1 started at: 16:54:13.977 (19/04/2024)
==== Processing Step Initialization started at 16:54:13.977 (19/04/2024)
##############################################################################
# #
# “minmed” is highly recommended for three images, #
# and is good for four to six images, #
# but should be avoided for ten or more images. #
# #
##############################################################################
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
********************************************************************************
*
* Estimated memory usage: up to 257 Mb.
* Output image size: 2336 X 2115 pixels.
* Output image file: ~ 56 Mb.
* Cores available: 4
*
********************************************************************************
==== Processing Step Initialization finished at 16:54:15.035 (19/04/2024)
==== Processing Step Static Mask started at 16:54:15.03 (19/04/2024)
==== Processing Step Static Mask finished at 16:54:15.158 (19/04/2024)
==== Processing Step Subtract Sky started at 16:54:15.158 (19/04/2024)
***** skymatch started on 2024-04-19 16:54:15.262437
Version 1.0.9
'skymatch' task will apply computed sky differences to input image file(s).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
'MDRIZSKY' header keyword will represent sky value *computed* from data.
----- User specified keywords: -----
Sky Value Keyword: 'MDRIZSKY'
Data Units Keyword: 'BUNIT'
----- Input file list: -----
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l6_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l7_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l9_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19la_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l6q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l6q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l7q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l7q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l9q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l9q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19laq_flt.fits'
EXT: 'SCI',1; MASK: ib1f19laq_skymatch_mask_sci1.fits[0]
----- Sky statistics parameters: -----
statistics function: 'median'
lower = None
upper = None
nclip = 5
lsigma = 4.0
usigma = 4.0
binwidth = 0.1
----- Data->Brightness conversion parameters for input files: -----
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l6q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l7q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l9q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19laq_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
----- Computing sky values requested image extensions (detector chips): -----
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
* Image: 'ib1f19l6q_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'ib1f19l7q_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'ib1f19l9q_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'ib1f19laq_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
***** skymatch ended on 2024-04-19 16:54:15.577828
TOTAL RUN TIME: 0:00:00.315391
==== Processing Step Subtract Sky finished at 16:54:15.803 (19/04/2024)
==== Processing Step Separate Drizzle started at 16:54:15.809 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 544.5 481.5
CD1_1 CD1_2 : -3.542984305540133e-05 3.723829036629122e-06
CD2_1 CD2_2 : 3.723829036629122e-06 3.542984305540133e-05
NAXIS : 1089 963
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_wht.fits
-Generating simple FITS output: ib1f19l6q_single_sci.fits
-Generating simple FITS output: ib1f19l7q_single_sci.fits
-Generating simple FITS output: ib1f19l9q_single_sci.fits
-Generating simple FITS output: ib1f19laq_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l6q_single_sci.fits
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
Writing out image to disk:
ib1f19l7q_single_sci.fits
ib1f19l9q_single_sci.fits
Writing out image to disk:
ib1f19laq_single_sci.fits
Writing out image to disk:
ib1f19l6q_single_wht.fits
Writing out image to disk:
Writing out image to disk:
ib1f19l7q_single_wht.fits
ib1f19l9q_single_wht.fits
Writing out image to disk:
ib1f19laq_single_wht.fits
==== Processing Step Separate Drizzle finished at 16:54:17.869 (19/04/2024)
==== Processing Step Create Median started at 16:54:17.87 (19/04/2024)
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits' is 464.93669456453324
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits' is 437.5836524801294
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits' is 427.0147483274913
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits' is 424.6136825690036
reference sky value for image 'ib1f19l6q_flt.fits' is 464.93669456453324
reference sky value for image 'ib1f19l7q_flt.fits' is 437.5836524801294
reference sky value for image 'ib1f19l9q_flt.fits' is 427.0147483274913
reference sky value for image 'ib1f19laq_flt.fits' is 424.6136825690036
Saving output median image to: 'f160w_0.8_med.fits'
==== Processing Step Create Median finished at 16:54:18.300 (19/04/2024)
==== Processing Step Blot started at 16:54:18.300 (19/04/2024)
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Blot: creating blotted image: ib1f19l6q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l6q_sci1_blt.fits
Writing out image to disk: ib1f19l6q_sci1_blt.fits
Blot: creating blotted image: ib1f19l7q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l7q_sci1_blt.fits
Writing out image to disk: ib1f19l7q_sci1_blt.fits
Blot: creating blotted image: ib1f19l9q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l9q_sci1_blt.fits
Writing out image to disk: ib1f19l9q_sci1_blt.fits
Blot: creating blotted image: ib1f19laq_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19laq_sci1_blt.fits
Writing out image to disk: ib1f19laq_sci1_blt.fits
==== Processing Step Blot finished at 16:54:19.837 (19/04/2024)
==== Processing Step Driz_CR started at 16:54:19.838 (19/04/2024)
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_crmask.fits
Creating output: ib1f19l6q_sci1_crmask.fits
Creating output: ib1f19l7q_sci1_crmask.fits
Creating output: ib1f19l9q_sci1_crmask.fits
Creating output: ib1f19laq_sci1_crmask.fits
==== Processing Step Driz_CR finished at 16:54:22.08 (19/04/2024)
==== Processing Step Final Drizzle started at 16:54:22.095 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
-Generating multi-extension output file: f160w_0.8_drz.fits
Deleted all instances of WCS with key A in extensions [1]
Deleted all instances of WCS with key B in extensions [1]
Deleted all instances of WCS with key C in extensions [1]
Writing out to disk: f160w_0.8_drz.fits
==== Processing Step Final Drizzle finished at 16:54:25.816 (19/04/2024)
AstroDrizzle Version 3.5.1 is finished processing at 16:54:25.817 (19/04/2024).
-------------------- --------------------
Step Elapsed time
-------------------- --------------------
Initialization 1.0573 sec.
Static Mask 0.1225 sec.
Subtract Sky 0.6448 sec.
Separate Drizzle 2.0597 sec.
Create Median 0.4304 sec.
Blot 1.5368 sec.
Driz_CR 2.2455 sec.
Final Drizzle 3.7206 sec.
==================== ====================
Total 11.8175 sec.
No trailer file saved...
No trailer file created...
AstroDrizzle log file: final.tra
AstroDrizzle Version 3.5.1 started at: 16:54:25.864 (19/04/2024)
==== Processing Step Initialization started at 16:54:25.864 (19/04/2024)
##############################################################################
# #
# “minmed” is highly recommended for three images, #
# and is good for four to six images, #
# but should be avoided for ten or more images. #
# #
##############################################################################
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
********************************************************************************
*
* Estimated memory usage: up to 257 Mb.
* Output image size: 2336 X 2115 pixels.
* Output image file: ~ 56 Mb.
* Cores available: 4
*
********************************************************************************
==== Processing Step Initialization finished at 16:54:26.909 (19/04/2024)
==== Processing Step Static Mask started at 16:54:26.909 (19/04/2024)
==== Processing Step Static Mask finished at 16:54:27.029 (19/04/2024)
==== Processing Step Subtract Sky started at 16:54:27.03 (19/04/2024)
***** skymatch started on 2024-04-19 16:54:27.131784
Version 1.0.9
'skymatch' task will apply computed sky differences to input image file(s).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
'MDRIZSKY' header keyword will represent sky value *computed* from data.
----- User specified keywords: -----
Sky Value Keyword: 'MDRIZSKY'
Data Units Keyword: 'BUNIT'
----- Input file list: -----
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l6_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l7_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l9_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19la_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l6q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l6q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l7q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l7q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l9q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l9q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19laq_flt.fits'
EXT: 'SCI',1; MASK: ib1f19laq_skymatch_mask_sci1.fits[0]
----- Sky statistics parameters: -----
statistics function: 'median'
lower = None
upper = None
nclip = 5
lsigma = 4.0
usigma = 4.0
binwidth = 0.1
----- Data->Brightness conversion parameters for input files: -----
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l6q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l7q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l9q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19laq_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
----- Computing sky values requested image extensions (detector chips): -----
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
* Image: 'ib1f19l6q_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'ib1f19l7q_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'ib1f19l9q_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'ib1f19laq_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
***** skymatch ended on 2024-04-19 16:54:27.449389
TOTAL RUN TIME: 0:00:00.317605
==== Processing Step Subtract Sky finished at 16:54:27.552 (19/04/2024)
==== Processing Step Separate Drizzle started at 16:54:27.553 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 544.5 481.5
CD1_1 CD1_2 : -3.542984305540133e-05 3.723829036629122e-06
CD2_1 CD2_2 : 3.723829036629122e-06 3.542984305540133e-05
NAXIS : 1089 963
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_wht.fits
Writing out image to disk:
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_wht.fits
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_wht.fits
-Generating simple FITS output: ib1f19l6q_single_sci.fits
-Generating simple FITS output: ib1f19l7q_single_sci.fits
-Generating simple FITS output: ib1f19l9q_single_sci.fits
-Generating simple FITS output: ib1f19laq_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l7q_single_sci.fits
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l6q_single_sci.fits
Writing out image to disk:
ib1f19l9q_single_sci.fits
Writing out image to disk:
ib1f19laq_single_sci.fits
Writing out image to disk:
ib1f19l7q_single_wht.fits
Writing out image to disk:
ib1f19l9q_single_wht.fits
Writing out image to disk:
ib1f19l6q_single_wht.fits
Writing out image to disk:
ib1f19laq_single_wht.fits
==== Processing Step Separate Drizzle finished at 16:54:29.612 (19/04/2024)
==== Processing Step Create Median started at 16:54:29.613 (19/04/2024)
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits' is 464.93669456453324
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits' is 437.5836524801294
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits' is 427.0147483274913
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits' is 424.6136825690036
reference sky value for image 'ib1f19l6q_flt.fits' is 464.93669456453324
reference sky value for image 'ib1f19l7q_flt.fits' is 437.5836524801294
reference sky value for image 'ib1f19l9q_flt.fits' is 427.0147483274913
reference sky value for image 'ib1f19laq_flt.fits' is 424.6136825690036
Saving output median image to: 'f160w_0.9_med.fits'
==== Processing Step Create Median finished at 16:54:30.044 (19/04/2024)
==== Processing Step Blot started at 16:54:30.044 (19/04/2024)
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Blot: creating blotted image: ib1f19l6q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l6q_sci1_blt.fits
Writing out image to disk: ib1f19l6q_sci1_blt.fits
Blot: creating blotted image: ib1f19l7q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l7q_sci1_blt.fits
Writing out image to disk: ib1f19l7q_sci1_blt.fits
Blot: creating blotted image: ib1f19l9q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l9q_sci1_blt.fits
Writing out image to disk: ib1f19l9q_sci1_blt.fits
Blot: creating blotted image: ib1f19laq_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19laq_sci1_blt.fits
Writing out image to disk: ib1f19laq_sci1_blt.fits
==== Processing Step Blot finished at 16:54:31.697 (19/04/2024)
==== Processing Step Driz_CR started at 16:54:31.698 (19/04/2024)
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_crmask.fits
Creating output: ib1f19l6q_sci1_crmask.fits
Creating output: ib1f19l7q_sci1_crmask.fits
Creating output: ib1f19laq_sci1_crmask.fits
Creating output: ib1f19l9q_sci1_crmask.fits
==== Processing Step Driz_CR finished at 16:54:33.979 (19/04/2024)
==== Processing Step Final Drizzle started at 16:54:33.99 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
-Generating multi-extension output file: f160w_0.9_drz.fits
Deleted all instances of WCS with key A in extensions [1]
Deleted all instances of WCS with key B in extensions [1]
Deleted all instances of WCS with key C in extensions [1]
Writing out to disk: f160w_0.9_drz.fits
==== Processing Step Final Drizzle finished at 16:54:37.829 (19/04/2024)
AstroDrizzle Version 3.5.1 is finished processing at 16:54:37.829 (19/04/2024).
-------------------- --------------------
Step Elapsed time
-------------------- --------------------
Initialization 1.0450 sec.
Static Mask 0.1199 sec.
Subtract Sky 0.5225 sec.
Separate Drizzle 2.0593 sec.
Create Median 0.4309 sec.
Blot 1.6533 sec.
Driz_CR 2.2814 sec.
Final Drizzle 3.8369 sec.
==================== ====================
Total 11.9492 sec.
No trailer file saved...
No trailer file created...
AstroDrizzle log file: final.tra
AstroDrizzle Version 3.5.1 started at: 16:54:37.878 (19/04/2024)
==== Processing Step Initialization started at 16:54:37.879 (19/04/2024)
##############################################################################
# #
# “minmed” is highly recommended for three images, #
# and is good for four to six images, #
# but should be avoided for ten or more images. #
# #
##############################################################################
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
********************************************************************************
*
* Estimated memory usage: up to 257 Mb.
* Output image size: 2336 X 2115 pixels.
* Output image file: ~ 56 Mb.
* Cores available: 4
*
********************************************************************************
==== Processing Step Initialization finished at 16:54:39.041 (19/04/2024)
==== Processing Step Static Mask started at 16:54:39.042 (19/04/2024)
==== Processing Step Static Mask finished at 16:54:39.161 (19/04/2024)
==== Processing Step Subtract Sky started at 16:54:39.162 (19/04/2024)
***** skymatch started on 2024-04-19 16:54:39.263611
Version 1.0.9
'skymatch' task will apply computed sky differences to input image file(s).
NOTE: Computed sky values WILL NOT be subtracted from image data ('subtractsky'=False).
'MDRIZSKY' header keyword will represent sky value *computed* from data.
----- User specified keywords: -----
Sky Value Keyword: 'MDRIZSKY'
Data Units Keyword: 'BUNIT'
----- Input file list: -----
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l6_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l7_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19l9_skymatch_mask_sci1.fits[0]
** Input image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits'
EXT: 'SCI',1; MASK: hst_11570_19_wfc3_ir_f160w_ib1f19la_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l6q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l6q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l7q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l7q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19l9q_flt.fits'
EXT: 'SCI',1; MASK: ib1f19l9q_skymatch_mask_sci1.fits[0]
** Input image: 'ib1f19laq_flt.fits'
EXT: 'SCI',1; MASK: ib1f19laq_skymatch_mask_sci1.fits[0]
----- Sky statistics parameters: -----
statistics function: 'median'
lower = None
upper = None
nclip = 5
lsigma = 4.0
usigma = 4.0
binwidth = 0.1
----- Data->Brightness conversion parameters for input files: -----
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l6q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l7q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19l9q_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
* Image: ib1f19laq_flt.fits
EXT = 'SCI',1
Data units type: COUNT-RATE
Conversion factor (data->brightness): 60.797431635711504
----- Computing sky values requested image extensions (detector chips): -----
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
* Image: 'ib1f19l6q_flt.fits['SCI',1]' -- SKY = 56.203832965924214 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.924444 NEW MDRIZSKY = 0.924444
* Image: 'ib1f19l7q_flt.fits['SCI',1]' -- SKY = 52.8972627889636 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.870058 NEW MDRIZSKY = 0.870058
* Image: 'ib1f19l9q_flt.fits['SCI',1]' -- SKY = 51.619641705121005 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.849043 NEW MDRIZSKY = 0.849043
* Image: 'ib1f19laq_flt.fits['SCI',1]' -- SKY = 51.32938907415446 (brightness units)
Sky change (data units):
- EXT = 'SCI',1 delta(MDRIZSKY) = 0.844269 NEW MDRIZSKY = 0.844269
***** skymatch ended on 2024-04-19 16:54:39.586315
TOTAL RUN TIME: 0:00:00.322704
==== Processing Step Subtract Sky finished at 16:54:39.692 (19/04/2024)
==== Processing Step Separate Drizzle started at 16:54:39.692 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 544.5 481.5
CD1_1 CD1_2 : -3.542984305540133e-05 3.723829036629122e-06
CD2_1 CD2_2 : 3.723829036629122e-06 3.542984305540133e-05
NAXIS : 1089 963
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key C in extensions [0]
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_sci.fits
Writing out image to disk:
Deleted all instances of WCS with key C in extensions [0]
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_sci.fits
Writing out image to disk:
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l6_single_wht.fits
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_sci.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l7_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19l9_single_wht.fits
Writing out image to disk:
hst_11570_19_wfc3_ir_f160w_ib1f19la_single_wht.fits
-Generating simple FITS output: ib1f19l6q_single_sci.fits
-Generating simple FITS output: ib1f19l7q_single_sci.fits
-Generating simple FITS output: ib1f19l9q_single_sci.fits
-Generating simple FITS output: ib1f19laq_single_sci.fits
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key A in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Deleted all instances of WCS with key B in extensions [0]
Writing out image to disk:
ib1f19l6q_single_sci.fits
Writing out image to disk:
ib1f19l9q_single_sci.fits
Writing out image to disk:
Writing out image to disk:
ib1f19l7q_single_sci.fits
ib1f19laq_single_sci.fits
Writing out image to disk:
Writing out image to disk:
ib1f19l6q_single_wht.fits
ib1f19l9q_single_wht.fits
Writing out image to disk:
ib1f19l7q_single_wht.fits
Writing out image to disk:
ib1f19laq_single_wht.fits
==== Processing Step Separate Drizzle finished at 16:54:41.748 (19/04/2024)
==== Processing Step Create Median started at 16:54:41.749 (19/04/2024)
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits' is 464.93669456453324
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits' is 437.5836524801294
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits' is 427.0147483274913
reference sky value for image 'hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits' is 424.6136825690036
reference sky value for image 'ib1f19l6q_flt.fits' is 464.93669456453324
reference sky value for image 'ib1f19l7q_flt.fits' is 437.5836524801294
reference sky value for image 'ib1f19l9q_flt.fits' is 427.0147483274913
reference sky value for image 'ib1f19laq_flt.fits' is 424.6136825690036
Saving output median image to: 'f160w_1.0_med.fits'
==== Processing Step Create Median finished at 16:54:42.178 (19/04/2024)
==== Processing Step Blot started at 16:54:42.178 (19/04/2024)
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l6_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l7_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19l9_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_blt.fits
Blot: creating blotted image: hst_11570_19_wfc3_ir_f160w_ib1f19la_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Writing out image to disk: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_blt.fits
Blot: creating blotted image: ib1f19l6q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l6q_sci1_blt.fits
Writing out image to disk: ib1f19l6q_sci1_blt.fits
Blot: creating blotted image: ib1f19l7q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l7q_sci1_blt.fits
Writing out image to disk: ib1f19l7q_sci1_blt.fits
Blot: creating blotted image: ib1f19l9q_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19l9q_sci1_blt.fits
Writing out image to disk: ib1f19l9q_sci1_blt.fits
Blot: creating blotted image: ib1f19laq_flt.fits[sci,1]
Using default C-based coordinate transformation...
-Generating simple FITS output: ib1f19laq_sci1_blt.fits
Writing out image to disk: ib1f19laq_sci1_blt.fits
==== Processing Step Blot finished at 16:54:43.716 (19/04/2024)
==== Processing Step Driz_CR started at 16:54:43.716 (19/04/2024)
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l6_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l7_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19l9_sci1_crmask.fits
Creating output: hst_11570_19_wfc3_ir_f160w_ib1f19la_sci1_crmask.fits
Creating output: ib1f19l6q_sci1_crmask.fits
Creating output: ib1f19l9q_sci1_crmask.fits
Creating output: ib1f19l7q_sci1_crmask.fits
Creating output: ib1f19laq_sci1_crmask.fits
==== Processing Step Driz_CR finished at 16:54:45.967 (19/04/2024)
==== Processing Step Final Drizzle started at 16:54:45.980 (19/04/2024)
WCS Keywords
Number of WCS axes: 2
CTYPE : 'RA---TAN' 'DEC--TAN'
CRVAL : 161.7693181775511 17.27337010993609
CRPIX : 1167.3829659085666 1056.5471233856526
CD1_1 CD1_2 : -1.8055555555555558e-05 -1.5854340902423206e-23
CD2_1 CD2_2 : -5.2878726791243137e-23 1.8055555555555558e-05
NAXIS : 2336 2115
-Generating multi-extension output file: f160w_1.0_drz.fits
Deleted all instances of WCS with key A in extensions [1]
Deleted all instances of WCS with key B in extensions [1]
Deleted all instances of WCS with key C in extensions [1]
Writing out to disk: f160w_1.0_drz.fits
==== Processing Step Final Drizzle finished at 16:54:50.161 (19/04/2024)
AstroDrizzle Version 3.5.1 is finished processing at 16:54:50.162 (19/04/2024).
-------------------- --------------------
Step Elapsed time
-------------------- --------------------
Initialization 1.1628 sec.
Static Mask 0.1196 sec.
Subtract Sky 0.5298 sec.
Separate Drizzle 2.0562 sec.
Create Median 0.4286 sec.
Blot 1.5372 sec.
Driz_CR 2.2514 sec.
Final Drizzle 4.1814 sec.
==================== ====================
Total 12.2669 sec.
No trailer file saved...
When evaluating what value to use for final_pixfrac
, THERE IS NO SINGLE METRIC THAT INDICATES WHAT VALUE TO USE, and several factors should be taken into account. The general philosophy is that the chosen value should improve the resolution of the image as much as possible, without causing any adverse effects.
The first thing to look for is an excessive number of holes in the science and weight images. The figure below shows the central region of the science and weight images produced by three different final_pixfrac
values.
with fits.open('f160w_0.1_drz.fits') as hdu1:
sci1 = hdu1[1].data
wht1 = hdu1[2].data
with fits.open('f160w_0.8_drz.fits') as hdu2:
sci2 = hdu2[1].data
wht2 = hdu2[2].data
with fits.open('f160w_1.0_drz.fits') as hdu3:
sci3 = hdu3[1].data
wht3 = hdu3[2].data
fig,ax = plt.subplots(2, 3, figsize=(16, 10), sharex=True, sharey=True)
fig.subplots_adjust(hspace=0, wspace=0)
norm3 = ImageNormalize(wht2, vmin=400, vmax=600, stretch=LinearStretch())
ax[0, 0].imshow(sci1, norm=norm2, cmap='gray')
ax[0, 0].text(1640, 1920, 'final_pixfrac=0.1', fontsize='20', color='gold')
ax[0, 1].imshow(sci2, norm=norm2, cmap='gray')
ax[0, 1].text(1640, 1920, 'final_pixfrac=0.8', fontsize='20', color='gold')
ax[0, 2].imshow(sci3, norm=norm2, cmap='gray')
ax[0, 2].text(1640, 1920, 'final_pixfrac=1.0', fontsize='20', color='gold')
ax[1, 0].imshow(wht1, norm=norm3, cmap='gray')
ax[1, 1].imshow(wht2, norm=norm3, cmap='gray')
ax[1, 2].imshow(wht3, norm=norm3, cmap='gray')
ax[0, 0].set_xlim(radeclims.to_pixel(im2wcs)[0])
ax[0, 0].set_ylim(radeclims.to_pixel(im2wcs)[1])
(1835.448117879787, 1934.1019662353683)
This figure above shows the central region of the science and weight images produced by final_pixfrac
values (left to right) of 0.1, 0.8 and 1.0. The top row is the science frame, the bottom is the weight image. The science image with the smallest final_pixfrac
value of 0.1 shows a noisy background and holes in the image where no input pixels fall into the output grid because pixfrac is too small. Inspection of the weight map corresponding to the smallest final_pixfrac
value shows many places with weights of zero, indicating that a final_pixfrac
value of 0.1 is clearly too small and was only included in this example for illustrative purposes.
Another piece of information that can be useful is the amount of noise in the weight image. As suggested in the HST Dither Handbook Chapter 5 and in the DrizzlePac Handbook Section 6.3.3, statistics performed on the drizzled weight image should yield a RMS/median value less than ~0.2. This threshold controls the trade-off between improving image resolution versus increasing background noise due to pixel resampling.
The figure below shows the RMS/median as a function of final_pixfrac
. One should take care to use the same region in the weight image as the region where the object of interest is located in the science image. If one is using the entire image for scientific analysis, then one should measure the statistics of the weight image where there is more variance.
whtlist = glob.glob('f160w_0.[0-9]*drz.fits')
std_med = np.empty(len(whtlist), dtype=float)
fraclist = np.empty(len(whtlist), dtype=float)
xlims = radeclims.to_pixel(im1wcs)[0].astype(int)
ylims = radeclims.to_pixel(im1wcs)[1].astype(int)
#Loop that measures statistics, also some information gathering
for i, im in enumerate(whtlist):
with fits.open(im) as hdu:
hdr = hdu[0].header
wht = hdu[2].data
if i == 0:
target = hdr['TARGNAME']
scale = str(hdr['D001SCAL'])
nimg = hdr['NDRIZIM']
wht_std = np.std(wht[ylims[0] : ylims[1], xlims[0] : xlims[1]])
wht_med = np.median(wht[ylims[0] : ylims[1], xlims[0] : xlims[1]])
std_med[i] = wht_std / wht_med
fraclist[i] = hdr['D001PIXF']
#Plotting commands
plt.clf()
plt.xlim(1.025, -0.025)
plt.ylim(0., 0.5)
plt.scatter(fraclist, std_med, s=50)
plt.axhline(0.2, ls='--', lw=3, c='r')
plt.xlabel('pixfrac', fontsize=18)
plt.ylabel('rms/median', fontsize=18)
plt.text(0.1, 0.45, target, fontsize=16, horizontalalignment='right')
plt.text(0.1, 0.40, 'Scale=' + str(scale) + '"', fontsize=16, horizontalalignment='right')
plt.text(0.1, 0.35, str(nimg) + ' images', fontsize=16, horizontalalignment='right')
Text(0.1, 0.35, '8 images')
The RMS/median increases steadily up to pixfrac of 0.6, where there is a change in the slope of the function, indicating that something odd might be happening beyond that point. Even though the RMS remains below 20% of the median all the way to final_pixfrac
of 0.1, the higher slope in the curve combined with the information from inspection of the science images themselvse (the evidence of holes) indicates that the pixfrac should remain >0.6.
6. Final thoughts#
The final_pixfrac
value has to be small enough to avoid degrading the combined output image, but large enough that when all images are ”dropped” onto the final frame, coverage of the output frame is fairly uniform. In general, final_pixfrac
should be slightly larger than the final output scale to allow some ”spillover” to adjacent pixels. This will help avoid ”holes” in the final product when a given pixel has been flagged as ”bad” in several frames. As a rule of thumb, statistics performed on the drizzled weight image in the region of interest should yield an RMS value (standard deviation) that is less than 20% of the median value. This threshold is a balance between the benefits of improving the image resolution at the expense of increasing noise in the background.
About this Notebook#
Author: R. Avila, STScI ACS Team
Updated: December 12, 2018