wfc3tools 1.6.0__py3-none-any.whl → 1.6.1__py3-none-any.whl

wfc3tools/calwf3.py

@@ -29,11 +29,13 @@ import subprocess
 
 # STSCI
 from stsci.tools import parseinput
+
 from .util import error_code
 
 
-def calwf3(input=None, printtime=False, save_tmp=False,
-           verbose=False, debug=False, parallel=True, version=False, log_func=print):
+def calwf3(
+    input=None, printtime=False, save_tmp=False, verbose=False, debug=False, parallel=True, version=False, log_func=print
+):
     """
     Run the calwf3.e executable as from the shell.
 
@@ -137,33 +139,32 @@ def calwf3(input=None, printtime=False, save_tmp=False,
 
     """
 
-    call_list = ['calwf3.e']
+    call_list = ["calwf3.e"]
     return_code = None
 
     if printtime:
-        call_list.append('-t')
+        call_list.append("-t")
 
     if save_tmp:
-        call_list.append('-s')
+        call_list.append("-s")
 
     if verbose:
-        call_list.append('-v')
+        call_list.append("-v")
 
     if version:
-        call_list.append('--version')
+        call_list.append("--version")
 
     if debug:
-        call_list.append('-d')
+        call_list.append("-d")
 
     if not parallel:
-        call_list.append('-1')
+        call_list.append("-1")
 
     infiles, dummy = parseinput.parseinput(input)
     if (len(parseinput.irafglob(input)) == 0) and not version:
         raise IOError("No valid image specified")
     if len(parseinput.irafglob(input)) > 1:
-        raise IOError("calwf3 can only accept 1 file for"
-                      "input at a time: {0}".format(infiles))
+        raise IOError("calwf3 can only accept 1 file forinput at a time: {0}".format(infiles))
 
     for image in infiles:
         if not os.path.exists(image):
@@ -180,7 +181,7 @@ def calwf3(input=None, printtime=False, save_tmp=False,
 
     if log_func is not None:
         for line in proc.stdout:
-            log_func(line.decode('utf8'))
+            log_func(line.decode("utf8"))
 
     return_code = proc.wait()
     ec = error_code(return_code)

wfc3tools/embedsub.py

@@ -1,13 +1,13 @@
-from .sub2full import sub2full
+import numpy
 
 # STDLIB
 from astropy.io import fits
-import os
-import numpy
 
 # STSCI
 from stsci.tools import parseinput
 
+from .sub2full import sub2full
+
 __taskname__ = "embedsub"
 
 
@@ -36,14 +36,13 @@ def embedsub(files):
 
     # process all the input subarrays
     for filename in infiles:
-
         # Make sure the input name conforms to normal style
-        if '_flt' not in filename:
+        if "_flt" not in filename:
             print("Warning: Can't properly parse '%s'; Skipping" % files)
 
         # Extract the root name and build SPT and output file names
-        root = filename[0:filename.find('_flt')]
-        full = root[0:len(root)-1] + 'f_flt.fits'
+        root = filename[0 : filename.find("_flt")]
+        full = root[0 : len(root) - 1] + "f_flt.fits"
 
         try:
             # open input file read-only
@@ -51,14 +50,13 @@ def embedsub(files):
         except EnvironmentError:
             print("Problem opening fits file %s" % (filename))
 
-        detector = flt[0].header['DETECTOR']
-        if 'UVIS' in detector:
+        detector = flt[0].header["DETECTOR"]
+        if "UVIS" in detector:
             uvis = True
 
         # compute subarray corners assuming the raw image location
         x1, x2, y1, y2 = sub2full(filename, fullExtent=True)[0]
-        print("Subarray image section [x1,x2,y1,y2] = [%d:%d,%d:%d]" % (x1, x2,
-                                                                        y1, y2))
+        print("Subarray image section [x1,x2,y1,y2] = [%d:%d,%d:%d]" % (x1, x2, y1, y2))
 
         if uvis:
             xaxis = uvis_full_x
@@ -74,44 +72,44 @@ def embedsub(files):
         err = numpy.zeros([xaxis, yaxis], dtype=numpy.float32)
         dq = numpy.zeros([xaxis, yaxis], dtype=numpy.int16) + 4
 
-        sci[y1-1:y2, x1-1:x2] = flt[1].data
-        err[y1-1:y2, x1-1:x2] = flt[2].data
-        dq[y1-1:y2, x1-1:x2] = flt[3].data
+        sci[y1 - 1 : y2, x1 - 1 : x2] = flt[1].data
+        err[y1 - 1 : y2, x1 - 1 : x2] = flt[2].data
+        dq[y1 - 1 : y2, x1 - 1 : x2] = flt[3].data
 
         if not uvis:
             samp = numpy.zeros([xaxis, yaxis], dtype=numpy.int16)
             time = numpy.zeros([xaxis, yaxis], dtype=numpy.float32)
-            samp[y1-1:y2, x1-1:x2] = flt[4].data
-            time[y1-1:y2, x1-1:x2] = flt[5].data
+            samp[y1 - 1 : y2, x1 - 1 : x2] = flt[4].data
+            time[y1 - 1 : y2, x1 - 1 : x2] = flt[5].data
 
         # Reset a few WCS values to make them appropriate for a
         # full-chip image
-        crpix1 = flt[1].header['CRPIX1']
-        crpix2 = flt[1].header['CRPIX2']
+        crpix1 = flt[1].header["CRPIX1"]
+        crpix2 = flt[1].header["CRPIX2"]
 
-        flt[1].header['sizaxis1'] = yaxis
-        flt[1].header['sizaxis2'] = xaxis
+        flt[1].header["sizaxis1"] = yaxis
+        flt[1].header["sizaxis2"] = xaxis
 
         for i in range(1, 4):
-            if 'CRPIX1' in flt[i].header:
-                flt[i].header['crpix1'] = crpix1 + x1 - 1
-                flt[i].header['crpix2'] = crpix2 + y1 - 1
-            if 'LTV1' in flt[i].header:
-                flt[i].header['ltv1'] = 0.0
-                flt[i].header['ltv2'] = 0.0
+            if "CRPIX1" in flt[i].header:
+                flt[i].header["crpix1"] = crpix1 + x1 - 1
+                flt[i].header["crpix2"] = crpix2 + y1 - 1
+            if "LTV1" in flt[i].header:
+                flt[i].header["ltv1"] = 0.0
+                flt[i].header["ltv2"] = 0.0
 
         # set the header value of SUBARRAY to False since it's now
         # regular size image
-        flt[0].header['SUBARRAY'] = False
-        
+        flt[0].header["SUBARRAY"] = False
+
         # Now write out the SCI, ERR, DQ extensions to the full-chip file
-        flt[1].data = sci 
+        flt[1].data = sci
         flt[2].data = err
         flt[3].data = dq
-        
+
         if not uvis:
             flt[4].data = samp
-            flt[5].data = time    
+            flt[5].data = time
 
         flt.writeto(full, overwrite=False)
 

wfc3tools/pstack.py

@@ -42,16 +42,14 @@ Usage:
 """
 
 # STDLIB
-import os
-from astropy.io import fits
 import numpy as np
+from astropy.io import fits
 from matplotlib import pyplot as plt
 
 plt.ion()
 
 
-def pstack(filename, column=0, row=0, extname="sci", units="counts",
-           title=None, xlabel=None, ylabel=None, plot=True):
+def pstack(filename, column=0, row=0, extname="sci", units="counts", title=None, xlabel=None, ylabel=None, plot=True):
     """
     A function to plot the statistics of one pixels up the IR ramp image.
     Original implementation in the iraf nicmos package. Pixel values here are
@@ -139,28 +137,28 @@ def pstack(filename, column=0, row=0, extname="sci", units="counts",
 
         for i in range(1, nsamp, 1):
             if time:
-                yaxis[i-1] = myfile["SCI", i].header['SAMPTIME']
+                yaxis[i - 1] = myfile["SCI", i].header["SAMPTIME"]
             else:
                 # Numpy is row-major with array indices written row-first
                 # (lexicographical access order)
-                yaxis[i-1] = myfile[extname.upper(), i].data[row, column]
-                xaxis[i-1] = myfile["SCI", i].header['SAMPTIME']
+                yaxis[i - 1] = myfile[extname.upper(), i].data[row, column]
+                xaxis[i - 1] = myfile["SCI", i].header["SAMPTIME"]
 
                 # convert to countrate
                 if "rate" in units.lower() and "/" not in bunit.lower():
-                    exptime = myfile["SCI", i].header['SAMPTIME']
-                    yaxis[i-1] /= exptime
+                    exptime = myfile["SCI", i].header["SAMPTIME"]
+                    yaxis[i - 1] /= exptime
                 # convert to counts
                 if "counts" in units.lower() and "/" in bunit.lower():
-                    exptime = myfile["SCI", i].header['SAMPTIME']
-                    yaxis[i-1] *= exptime
+                    exptime = myfile["SCI", i].header["SAMPTIME"]
+                    yaxis[i - 1] *= exptime
 
     if not ylabel:
         if "rate" in units.lower():
             if "/" in bunit.lower():
                 ylabel = bunit
             else:
-                ylabel = bunit+" per second"
+                ylabel = bunit + " per second"
         else:
             if "/" in bunit:
                 stop_index = bunit.find("/")
@@ -177,8 +175,7 @@ def pstack(filename, column=0, row=0, extname="sci", units="counts",
             plt.xlabel("Sample time")
 
         if not title:
-            title = "%s   Pixel stack for col=%d, row=%d" % (filename, column,
-                                                             row)
+            title = "%s   Pixel stack for col=%d, row=%d" % (filename, column, row)
         plt.title(title)
 
         if time:

wfc3tools/pstat.py

@@ -37,7 +37,7 @@ Usage:
     Note that the arrays are structured in SCI order, so the final exposure is the first element in the array.
 
 .. Warning::
-    The interface to this utility has been updated from previous versions and 
+    The interface to this utility has been updated from previous versions and
     is **not backwards compatible.**  Here is an example to illustrate the "original"
     syntax, the "original syntax corrected for row/column order", and finally the
     "new" syntax which requires column and row sections to be specified as tuples.
@@ -56,17 +56,27 @@ Usage:
 """
 
 # STDLIB
-import os
-from astropy.io import fits
 import numpy as np
+from astropy.io import fits
 from matplotlib import pyplot as plt
 from scipy.stats import mode as mode
 
 plt.ion()
 
 
-def pstat(filename, col_slice=None, row_slice=None, extname="sci", units="counts",
-          stat="midpt", title=None, xlabel=None, ylabel=None, plot=True, overplot=False):
+def pstat(
+    filename,
+    col_slice=None,
+    row_slice=None,
+    extname="sci",
+    units="counts",
+    stat="midpt",
+    title=None,
+    xlabel=None,
+    ylabel=None,
+    plot=True,
+    overplot=False,
+):
     """
     A function to plot the statistics of one or more pixels up an IR ramp.
 
@@ -147,9 +157,9 @@ def pstat(filename, col_slice=None, row_slice=None, extname="sci", units="counts
     bracket_loc = filename.find("[")
 
     # get the base filename and strip off any extra information as necessary
-    if (bracket_loc < 0):
+    if bracket_loc < 0:
         imagename = filename
-    elif (bracket_loc > 0):
+    elif bracket_loc > 0:
         imagename = filename[:bracket_loc]
         print("Any extension name or image section must be specified via parameters.")
         print("Input filename has been stripped of data in brackets, %s" % (imagename))
@@ -169,14 +179,14 @@ def pstat(filename, col_slice=None, row_slice=None, extname="sci", units="counts
     all_cols = False
     if not col_slice:
         all_cols = True
-    elif not(isinstance(col_slice, tuple) and len(col_slice) == 2 and all(isinstance(val, int) for val in col_slice)):
+    elif not (isinstance(col_slice, tuple) and len(col_slice) == 2 and all(isinstance(val, int) for val in col_slice)):
         print("Invalid specification for col_slice which must be a tuple of two integer values.")
         return 0, 0
 
     all_rows = False
     if not row_slice:
         all_rows = True
-    elif not(isinstance(row_slice, tuple) and len(row_slice) == 2 and all(isinstance(val, int) for val in row_slice)):
+    elif not (isinstance(row_slice, tuple) and len(row_slice) == 2 and all(isinstance(val, int) for val in row_slice)):
         print("Invalid specification for row_slice which must be a tuple of two integer values.")
         return 0, 0
 
@@ -207,33 +217,32 @@ def pstat(filename, col_slice=None, row_slice=None, extname="sci", units="counts
 
         for i in range(1, nsamp, 1):
             if "midpt" in stat:
-                yaxis[i-1] = np.median(myfile[extname.upper(), i].data[ystart:yend, xstart:xend])
+                yaxis[i - 1] = np.median(myfile[extname.upper(), i].data[ystart:yend, xstart:xend])
 
             if "mean" in stat:
-                yaxis[i-1] = np.mean(myfile[extname.upper(), i].data[ystart:yend, xstart:xend])
+                yaxis[i - 1] = np.mean(myfile[extname.upper(), i].data[ystart:yend, xstart:xend])
 
             if "mode" in stat:
-                yaxis[i-1] = mode(myfile[extname.upper(), i].data[ystart:yend, xstart:xend],
-                                  axis=None)[0]
+                yaxis[i - 1] = mode(myfile[extname.upper(), i].data[ystart:yend, xstart:xend], axis=None)[0]
 
             if "min" in stat:
-                yaxis[i-1] = np.min(myfile[extname.upper(), i].data[ystart:yend, xstart:xend])
+                yaxis[i - 1] = np.min(myfile[extname.upper(), i].data[ystart:yend, xstart:xend])
 
             if "max" in stat:
-                yaxis[i-1] = np.max(myfile[extname.upper(), i].data[ystart:yend, xstart:xend])
+                yaxis[i - 1] = np.max(myfile[extname.upper(), i].data[ystart:yend, xstart:xend])
 
             if "stddev" in stat:
-                yaxis[i-1] = np.std(myfile[extname.upper(), i].data[ystart:yend, xstart:xend])
+                yaxis[i - 1] = np.std(myfile[extname.upper(), i].data[ystart:yend, xstart:xend])
 
-            exptime = myfile["SCI", i].header['SAMPTIME']
-            xaxis[i-1] = exptime
+            exptime = myfile["SCI", i].header["SAMPTIME"]
+            xaxis[i - 1] = exptime
 
             # convert to countrate
             if "rate" in units.lower() and "/" not in bunit.lower():
-                yaxis[i-1] /= exptime
+                yaxis[i - 1] /= exptime
             # convert to counts
             if "counts" in units.lower() and "/" in bunit.lower():
-                yaxis[i-1] *= exptime
+                yaxis[i - 1] *= exptime
 
     if plot:
         if not overplot:
@@ -251,15 +260,14 @@ def pstat(filename, col_slice=None, row_slice=None, extname="sci", units="counts
                 else:
                     ylabel = bunit
 
-        ylabel += ("   %s" % (stat))
+        ylabel += "   %s" % (stat)
         plt.ylabel(ylabel)
 
         if not xlabel:
             plt.xlabel("Sample time (s)")
 
         if not title:
-            title = "%s   Pixel stats for [%d:%d,%d:%d]" % (imagename, xstart,
-                                                            xend, ystart, yend)
+            title = "%s   Pixel stats for [%d:%d,%d:%d]" % (imagename, xstart, xend, ystart, yend)
         plt.title(title)
         plt.plot(xaxis, yaxis, "+")
         plt.draw()

wfc3tools/sampinfo.py

@@ -74,9 +74,9 @@ Usage:
 """
 
 # STDLIB
-from astropy.io import fits
-import os
+
 import numpy as np
+from astropy.io import fits
 
 # STSCI
 from stsci.tools import parseinput
@@ -132,11 +132,11 @@ def sampinfo(imagelist, add_keys=None, mean=False, median=False):
         ir_list += add_keys
 
     # measure the min and max data
-    if (mean):
-        if (add_keys):
-            if ("DATAMIN" not in add_keys):
+    if mean:
+        if add_keys:
+            if "DATAMIN" not in add_keys:
                 ir_list += ["DATAMIN"]
-            if ("DATAMAX" not in add_keys):
+            if "DATAMAX" not in add_keys:
                 ir_list += ["DATAMAX"]
         else:
             ir_list += ["DATAMIN", "DATAMAX"]
@@ -155,19 +155,18 @@ def sampinfo(imagelist, add_keys=None, mean=False, median=False):
         samp_seq = header0["SAMP_SEQ"]
 
         print("IMAGE\t\t\tNEXTEND\tSAMP_SEQ\tNSAMP\tEXPTIME")
-        print("%s\t%d\t%s\t\t%d\t%f\n" % (image, nextend, samp_seq,
-                                          nsamp, exptime))
+        print("%s\t%d\t%s\t\t%d\t%f\n" % (image, nextend, samp_seq, nsamp, exptime))
         printline = "IMSET\tSAMPNUM"
 
         for key in ir_list:
-            printline += ("\t"+key)
+            printline += "\t" + key
         print(printline)
 
         # loop through all the samples for the image and print stuff as we go
-        for samp in range(1, nsamp+1, 1):
+        for samp in range(1, nsamp + 1, 1):
             printline = ""
             printline += str(samp)
-            printline += ("\t"+str(nsamp-samp))
+            printline += "\t" + str(nsamp - samp)
             for key in ir_list:
                 if "DATAMIN" in key:
                     datamin = True
@@ -176,16 +175,15 @@ def sampinfo(imagelist, add_keys=None, mean=False, median=False):
                     datamax = True
                     datamaxval = np.max(current["SCI", samp].data)
                 try:
-                    printline += ("\t"+str(current["SCI", samp].header[key]))
+                    printline += "\t" + str(current["SCI", samp].header[key])
                 except KeyError:
                     try:
-                        printline += ("\t"+str(current[0].header[key]))
+                        printline += "\t" + str(current[0].header[key])
                     except KeyError as e:
-                        printline += ("\tNA")
-            if (datamin and datamax):
-                printline += ("\tAvgPixel: "+str((dataminval+datamaxval)/2.))
-            if (median):
-                printline += ("\tMedPixel: "+str(np.median(current["SCI",
-                                                           samp].data)))
+                        printline += "\tNA"
+            if datamin and datamax:
+                printline += "\tAvgPixel: " + str((dataminval + datamaxval) / 2.0)
+            if median:
+                printline += "\tMedPixel: " + str(np.median(current["SCI", samp].data))
             print(printline)
         current.close()

wfc3tools/sub2full.py

@@ -36,8 +36,9 @@ Usage:
 """
 
 # STDLIB
-from astropy.io import fits
 import os
+
+from astropy.io import fits
 from stsci.tools import parseinput
 
 
@@ -90,8 +91,7 @@ def sub2full(filename, x=None, y=None, fullExtent=False):
     coords = list()
 
     for f in infiles:
-        spt = os.path.join(os.path.dirname(f), os.path.basename(f)[0:9] +
-                           '_spt.fits')
+        spt = os.path.join(os.path.dirname(f), os.path.basename(f)[0:9] + "_spt.fits")
         uvis_x_size = 2051
         serial_over = 25.0
         ir_overscan = 5.0
@@ -100,16 +100,16 @@ def sub2full(filename, x=None, y=None, fullExtent=False):
         try:
             fd2 = fits.open(spt)
         except (ValueError, IOError) as e:
-            raise ValueError('%s ' % (e))
+            raise ValueError("%s " % (e))
 
         # check for required keywords and close the images
         try:
-            detector = fd2[0].header['SS_DTCTR']
-            subarray = fd2[0].header['SS_SUBAR']
-            xcorner = int(fd2[1].header['XCORNER'])
-            ycorner = int(fd2[1].header['YCORNER'])
-            numrows = int(fd2[1].header['NUMROWS'])
-            numcols = int(fd2[1].header['NUMCOLS'])
+            detector = fd2[0].header["SS_DTCTR"]
+            subarray = fd2[0].header["SS_SUBAR"]
+            xcorner = int(fd2[1].header["XCORNER"])
+            ycorner = int(fd2[1].header["YCORNER"])
+            numrows = int(fd2[1].header["NUMROWS"])
+            numcols = int(fd2[1].header["NUMCOLS"])
             fd2.close()
         except KeyError as e:
             raise KeyError(f"Required header keyword missing: {e}")
@@ -120,7 +120,7 @@ def sub2full(filename, x=None, y=None, fullExtent=False):
         sizaxis1 = numcols
         sizaxis2 = numrows
 
-        if (xcorner == 0 and ycorner == 0):
+        if xcorner == 0 and ycorner == 0:
             cornera1 = 0
             cornera2 = 0
             cornera1a = cornera1 + 1
@@ -128,7 +128,7 @@ def sub2full(filename, x=None, y=None, fullExtent=False):
             cornera2a = cornera2 + 1
             cornera2b = cornera2a + sizaxis2 - 1
         else:
-            if 'UVIS' in detector:
+            if "UVIS" in detector:
                 cornera1 = ycorner
                 cornera2 = uvis_x_size - xcorner - sizaxis2
                 if xcorner >= uvis_x_size:
@@ -152,17 +152,16 @@ def sub2full(filename, x=None, y=None, fullExtent=False):
                 cornera2a = cornera2 + 1
                 cornera2b = cornera2a + sizaxis2 - 11
 
-        if (x or y):
-            if ((not isinstance(x, int) or (not isinstance(y, int)))):
+        if x or y:
+            if not isinstance(x, int) or (not isinstance(y, int)):
                 raise ValueError("Must input integer value for x and y ")
             else:
                 cornera1a = cornera1a + x
                 cornera2a = cornera2a + y
                 fullExtent = False
 
-        if (fullExtent):
-            coords.append((int(cornera1a), int(cornera1b), int(cornera2a),
-                           int(cornera2b)))
+        if fullExtent:
+            coords.append((int(cornera1a), int(cornera1b), int(cornera2a), int(cornera2b)))
         else:
             coords.append((int(cornera1a), int(cornera2a)))
 

wfc3tools/util.py

@@ -7,50 +7,54 @@ __all__ = ["display_help"]
 
 
 def display_help():
-    """ display local html help in a browser window"""
-    url = "http://wfc3tools.readthedocs.io/"
+    """display local html help in a browser window"""
+    url = "https://wfc3tools.readthedocs.io/"
     print(url)
     try:
         import webbrowser
+
         # grab the version that's installed
-        if "dev" not in version('wfc3tools'):
-            url += "en/{0:s}/".format(version('wfc3tools'))
+        if "dev" not in version("wfc3tools"):
+            url += "en/{0:s}/".format(version("wfc3tools"))
         webbrowser.open(url)
     except ImportError:
-        warnings.warn("webbrowser module not installed, see {0:s} help \
-                       pages".format(url))
+        warnings.warn(
+            "webbrowser module not installed, see {0:s} help \
+                       pages".format(url)
+        )
 
 
 def error_code(code=None):
-    """ return the error code text, or all if code is None."""
-
-    codes = {2: "ERROR_RETURN",
-             111: "OUT_OF_MEMORY",
-             114: "OPEN_FAILED",
-             115: "CAL_FILE_MISSING",
-             116: "NOTHIN_TO_DO",
-             117: "KEYWORD_MISSING",
-             118: "ALLOCATION_PROBLEM",
-             119: "HEADER_PROBLEM",
-             120: "SIZE_MISMATCH",
-             130: "CAL_STEP_NOT_DONE",
-             141: "TABLE_ERROR",
-             142: "COLUMN_NOT_FOUND",
-             143: "ELEMENT_NOT_FOUND",
-             144: "ROW_NOT_FOUND",
-             151: "NO_GOOD_DATA",
-             152: "NO_CHIP_FOUND",
-             171: "REF_TOO_SMALL",
-             999: "INTERNAL_ERROR",
-             1001: "INVALID_EXPTIME",
-             1011: "INVALID_FILENAME",
-             1020: "WRITE_FAILED",
-             1021: "INVALID_TEMP_FILE",
-             1023: "FILE_NOT_READABLE",
-             1025: "COPY_NOT_POSSIBLE",
-             1111: "INVALID_VALUE",
-             1030: "UNSUPPORTED_APERTURE",
-             }
+    """return the error code text, or all if code is None."""
+
+    codes = {
+        2: "ERROR_RETURN",
+        111: "OUT_OF_MEMORY",
+        114: "OPEN_FAILED",
+        115: "CAL_FILE_MISSING",
+        116: "NOTHIN_TO_DO",
+        117: "KEYWORD_MISSING",
+        118: "ALLOCATION_PROBLEM",
+        119: "HEADER_PROBLEM",
+        120: "SIZE_MISMATCH",
+        130: "CAL_STEP_NOT_DONE",
+        141: "TABLE_ERROR",
+        142: "COLUMN_NOT_FOUND",
+        143: "ELEMENT_NOT_FOUND",
+        144: "ROW_NOT_FOUND",
+        151: "NO_GOOD_DATA",
+        152: "NO_CHIP_FOUND",
+        171: "REF_TOO_SMALL",
+        999: "INTERNAL_ERROR",
+        1001: "INVALID_EXPTIME",
+        1011: "INVALID_FILENAME",
+        1020: "WRITE_FAILED",
+        1021: "INVALID_TEMP_FILE",
+        1023: "FILE_NOT_READABLE",
+        1025: "COPY_NOT_POSSIBLE",
+        1111: "INVALID_VALUE",
+        1030: "UNSUPPORTED_APERTURE",
+    }
 
     if code is None:
         return codes

wfc3tools/version.py

@@ -28,7 +28,7 @@ version_tuple: VERSION_TUPLE
 commit_id: COMMIT_ID
 __commit_id__: COMMIT_ID
 
-__version__ = version = '1.6.0'
-__version_tuple__ = version_tuple = (1, 6, 0)
+__version__ = version = '1.6.1'
+__version_tuple__ = version_tuple = (1, 6, 1)
 
-__commit_id__ = commit_id = 'g1206755c4'
+__commit_id__ = commit_id = 'gd6931c975'

wfc3tools/wf32d.py

@@ -40,12 +40,23 @@ import subprocess
 
 # STSCI
 from stsci.tools import parseinput
+
 from .util import error_code
 
 
-def wf32d(input, output=None, dqicorr="PERFORM", darkcorr="PERFORM",
-          flatcorr="PERFORM", shadcorr="PERFORM", photcorr="PERFORM",
-          verbose=False, quiet=True, debug=False, log_func=print):
+def wf32d(
+    input,
+    output=None,
+    dqicorr="PERFORM",
+    darkcorr="PERFORM",
+    flatcorr="PERFORM",
+    shadcorr="PERFORM",
+    photcorr="PERFORM",
+    verbose=False,
+    quiet=True,
+    debug=False,
+    log_func=print,
+):
     """
     Call the wf32d.e executable.
 
@@ -115,29 +126,29 @@ def wf32d(input, output=None, dqicorr="PERFORM", darkcorr="PERFORM",
 
     """
 
-    call_list = ['wf32d.e']
+    call_list = ["wf32d.e"]
     return_code = None
 
     if verbose:
-        call_list += ['-v', '-t']
+        call_list += ["-v", "-t"]
 
     if debug:
-        call_list.append('-d')
+        call_list.append("-d")
 
-    if (darkcorr == "PERFORM"):
-        call_list.append('-dark')
+    if darkcorr == "PERFORM":
+        call_list.append("-dark")
 
-    if (dqicorr == "PERFORM"):
-        call_list.append('-dqi')
+    if dqicorr == "PERFORM":
+        call_list.append("-dqi")
 
-    if (flatcorr == "PERFORM"):
-        call_list.append('-flat')
+    if flatcorr == "PERFORM":
+        call_list.append("-flat")
 
-    if (shadcorr == "PERFORM"):
-        call_list.append('-shad')
+    if shadcorr == "PERFORM":
+        call_list.append("-shad")
 
-    if (photcorr == "PERFORM"):
-        call_list.append('-phot')
+    if photcorr == "PERFORM":
+        call_list.append("-phot")
 
     infiles, dummy = parseinput.parseinput(input)
     if "_asn" in input:
@@ -145,8 +156,7 @@ def wf32d(input, output=None, dqicorr="PERFORM", darkcorr="PERFORM",
     if len(parseinput.irafglob(input)) == 0:
         raise IOError("No valid image specified")
     if len(parseinput.irafglob(input)) > 1:
-        raise IOError("wf32d can only accept 1 file for"
-                      "input at a time: {0}".format(infiles))
+        raise IOError("wf32d can only accept 1 file forinput at a time: {0}".format(infiles))
 
     for image in infiles:
         if not os.path.exists(image):
@@ -164,7 +174,7 @@ def wf32d(input, output=None, dqicorr="PERFORM", darkcorr="PERFORM",
     )
     if log_func is not None:
         for line in proc.stdout:
-            log_func(line.decode('utf8'))
+            log_func(line.decode("utf8"))
 
     return_code = proc.wait()
     ec = error_code(return_code)

wfc3tools/wf3ccd.py

@@ -48,12 +48,22 @@ import subprocess
 
 # STSCI
 from stsci.tools import parseinput
+
 from .util import error_code
 
 
-def wf3ccd(input, output=None, dqicorr="PERFORM", atodcorr="PERFORM",
-           blevcorr="PERFORM", biascorr="PERFORM", flashcorr="PERFORM",
-           verbose=False, quiet=True, log_func=print):
+def wf3ccd(
+    input,
+    output=None,
+    dqicorr="PERFORM",
+    atodcorr="PERFORM",
+    blevcorr="PERFORM",
+    biascorr="PERFORM",
+    flashcorr="PERFORM",
+    verbose=False,
+    quiet=True,
+    log_func=print,
+):
     """
     Run the ``wf3ccd.e`` executable as from the shell.
 
@@ -114,26 +124,26 @@ def wf3ccd(input, output=None, dqicorr="PERFORM", atodcorr="PERFORM",
 
     """
 
-    call_list = ['wf3ccd.e']
+    call_list = ["wf3ccd.e"]
     return_code = None
 
     if verbose:
-        call_list += ['-v', '-t']
+        call_list += ["-v", "-t"]
 
-    if (dqicorr == "PERFORM"):
-        call_list.append('-dqi')
+    if dqicorr == "PERFORM":
+        call_list.append("-dqi")
 
-    if (atodcorr == "PERFORM"):
-        call_list.append('-atod')
+    if atodcorr == "PERFORM":
+        call_list.append("-atod")
 
-    if (blevcorr == "PERFORM"):
-        call_list.append('-blev')
+    if blevcorr == "PERFORM":
+        call_list.append("-blev")
 
-    if (biascorr == "PERFORM"):
-        call_list.append('-bias')
+    if biascorr == "PERFORM":
+        call_list.append("-bias")
 
-    if (flashcorr == "PERFORM"):
-        call_list.append('-flash')
+    if flashcorr == "PERFORM":
+        call_list.append("-flash")
 
     infiles, dummy = parseinput.parseinput(input)
     if "_asn" in input:
@@ -141,8 +151,7 @@ def wf3ccd(input, output=None, dqicorr="PERFORM", atodcorr="PERFORM",
     if len(parseinput.irafglob(input)) == 0:
         raise IOError("No valid image specified")
     if len(parseinput.irafglob(input)) > 1:
-        raise IOError("wf3ccd can only accept 1 file for"
-                      "input at a time: {0}".format(infiles))
+        raise IOError("wf3ccd can only accept 1 file forinput at a time: {0}".format(infiles))
 
     for image in infiles:
         if not os.path.exists(image):
@@ -160,7 +169,7 @@ def wf3ccd(input, output=None, dqicorr="PERFORM", atodcorr="PERFORM",
     )
     if log_func is not None:
         for line in proc.stdout:
-            log_func(line.decode('utf8'))
+            log_func(line.decode("utf8"))
 
     return_code = proc.wait()
     ec = error_code(return_code)

wfc3tools/wf3cte.py

@@ -14,7 +14,6 @@ wf3cte:
 """
 
 # STDLIB
-import os.path
 import subprocess
 
 # STSCI
@@ -61,16 +60,16 @@ def wf3cte(input, parallel=True, verbose=False, log_func=print):
 
     """
 
-    call_list = ['wf3cte.e']
+    call_list = ["wf3cte.e"]
 
     if verbose:
-        call_list.append('-v')
+        call_list.append("-v")
 
     if not parallel:
-        call_list.append('-1')
+        call_list.append("-1")
 
     infiles, dummy_out = parseinput.parseinput(input)
-    call_list.append(','.join(infiles))
+    call_list.append(",".join(infiles))
 
     print(call_list)
 
@@ -81,7 +80,7 @@ def wf3cte(input, parallel=True, verbose=False, log_func=print):
     )
     if log_func is not None:
         for line in proc.stdout:
-            log_func(line.decode('utf8'))
+            log_func(line.decode("utf8"))
 
     return_code = proc.wait()
     if return_code != 0:

wfc3tools/wf3ir.py

@@ -42,6 +42,7 @@ import subprocess
 
 # STSCI
 from stsci.tools import parseinput
+
 from .util import error_code
 
 
@@ -96,11 +97,11 @@ def wf3ir(input, output=None, verbose=False, quiet=True, log_func=print):
 
     """
 
-    call_list = ['wf3ir.e']
+    call_list = ["wf3ir.e"]
     return_code = None
 
     if verbose:
-        call_list += ['-v', '-t']
+        call_list += ["-v", "-t"]
 
     infiles, dummy = parseinput.parseinput(input)
     if "_asn" in input:
@@ -108,8 +109,7 @@ def wf3ir(input, output=None, verbose=False, quiet=True, log_func=print):
     if len(parseinput.irafglob(input)) == 0:
         raise IOError("No valid image specified")
     if len(parseinput.irafglob(input)) > 1:
-        raise IOError("wf3ir can only accept 1 file for"
-                      "input at a time: {0}".format(infiles))
+        raise IOError("wf3ir can only accept 1 file forinput at a time: {0}".format(infiles))
 
     for image in infiles:
         if not os.path.exists(image):
@@ -127,7 +127,7 @@ def wf3ir(input, output=None, verbose=False, quiet=True, log_func=print):
     )
     if log_func is not None:
         for line in proc.stdout:
-            log_func(line.decode('utf8'))
+            log_func(line.decode("utf8"))
 
     return_code = proc.wait()
     ec = error_code(return_code)

wfc3tools/wf3rej.py

@@ -39,13 +39,26 @@ import subprocess
 
 # STSCI
 from stsci.tools import parseinput
+
 from .util import error_code
 
 
-def wf3rej(input, output, crrejtab="", scalense=0., initgues="",
-           skysub="", crsigmas="", crradius=0., crthresh=0.,
-           badinpdq=0, crmask=False, shadcorr=False, verbose=False,
-           log_func=print):
+def wf3rej(
+    input,
+    output,
+    crrejtab="",
+    scalense=0.0,
+    initgues="",
+    skysub="",
+    crsigmas="",
+    crradius=0.0,
+    crthresh=0.0,
+    badinpdq=0,
+    crmask=False,
+    shadcorr=False,
+    verbose=False,
+    log_func=print,
+):
     """
     wf3rej, the cosmic-ray rejection and image combination task in calwf3,
     combines CR-SPLIT or REPEAT-OBS exposures into a single image, first
@@ -136,7 +149,7 @@ def wf3rej(input, output, crrejtab="", scalense=0., initgues="",
             raise IOError("Input file not found: {0}".format(image))
 
     # Generate a comma-separated string of the input filenames
-    input = ','.join(infiles)
+    input = ",".join(infiles)
 
     call_list.append(input)
 
@@ -147,46 +160,46 @@ def wf3rej(input, output, crrejtab="", scalense=0., initgues="",
         call_list.append("-v")
         call_list.append("-t")
 
-    if (shadcorr):
+    if shadcorr:
         call_list.append("-shadcorr")
 
-    if (crmask):
+    if crmask:
         call_list.append("-crmask")
 
-    if (crrejtab != ""):
+    if crrejtab != "":
         call_list += ["-table", crrejtab]
 
-    if (scalense != ""):
+    if scalense != "":
         call_list += ["-scale", str(scalense)]
 
-    if (initgues != ""):
+    if initgues != "":
         options = ["min", "med"]
         if initgues not in options:
             raise ValueError("Invalid option for initgues")
         else:
             call_list += ["-init", str(initgues)]
 
-    if (skysub != ""):
+    if skysub != "":
         options = ["none", "mode", "median"]
         if skysub not in options:
             raise ValueError(f"Invalid skysub option {options}: {skysub}")
         else:
             call_list += ["-sky", str(skysub)]
 
-    if (crsigmas != ""):
+    if crsigmas != "":
         call_list += ["-sigmas", str(crsigmas)]
 
-    if (crradius >= 0.):
+    if crradius >= 0.0:
         call_list += ["-radius", str(crradius)]
     else:
         raise ValueError("Invalid crradius specified")
 
-    if (crthresh >= 0.):
+    if crthresh >= 0.0:
         call_list += ["-thresh", str(crthresh)]
     else:
         raise ValueError("Invalid crthresh specified")
 
-    if (badinpdq >= 0):
+    if badinpdq >= 0:
         call_list += ["-pdq", str(badinpdq)]
 
     else:
@@ -199,7 +212,7 @@ def wf3rej(input, output, crrejtab="", scalense=0., initgues="",
     )
     if log_func is not None:
         for line in proc.stdout:
-            log_func(line.decode('utf8'))
+            log_func(line.decode("utf8"))
 
     return_code = proc.wait()
     ec = error_code(return_code)

{wfc3tools-1.6.0.dist-info → wfc3tools-1.6.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: wfc3tools
-Version: 1.6.0
+Version: 1.6.1
 Summary: Python Tools for HST WFC3 Data
 Author: STScI
 Project-URL: Homepage, http://wfc3tools.readthedocs.io/
@@ -23,31 +23,30 @@ Requires-Dist: scipy>=1.8.0
 Provides-Extra: tests
 Requires-Dist: pytest; extra == "tests"
 Requires-Dist: astroquery; extra == "tests"
+Requires-Dist: ci-watson; extra == "tests"
 Provides-Extra: docs
 Requires-Dist: sphinx; extra == "docs"
 Requires-Dist: sphinx-automodapi; extra == "docs"
-Requires-Dist: sphinx-rtd-theme; extra == "docs"
+Requires-Dist: sphinx-rtd-theme>1.2.0; extra == "docs"
 Dynamic: license-file
 
 WFC3TOOLS
 =========
 
+.. image:: https://github.com/spacetelescope/wfc3tools/actions/workflows/tests.yml/badge.svg?branch=main
+    :target: https://github.com/spacetelescope/wfc3tools/actions/workflows/tests.yml
+    :alt: tests
+
 .. image:: https://readthedocs.org/projects/wfc3tools/badge/?version=latest
-    :target: http://wfc3tools.readthedocs.io/en/latest/?badge=latest
+    :target: https://wfc3tools.readthedocs.io/en/latest/?badge=latest
     :alt: Documentation Status
 
-.. image:: http://img.shields.io/badge/powered%20by-AstroPy-orange.svg?style=flat
-    :target: http://www.astropy.org
+.. image:: https://img.shields.io/badge/powered%20by-AstroPy-orange.svg?style=flat
+    :target: https://www.astropy.org
     :alt: Powered by Astropy Badge
 
 
-For more information please see the `online documentation <http://wfc3tools.readthedocs.io/>`_.
-
-You can also display the docs locally after install, import ``wfc3tools`` and then issue the following command to display the help docs in your local browser:
-
-::
-
-    wfc3tools.display_help()
+For more information please see the `online documentation <https://wfc3tools.readthedocs.io/>`_.
 
 To install using pip:
 
@@ -73,19 +72,11 @@ or haven't used github before, please feel free to contact `HST Help Desk <https
 Feedback and feature requests?   Is there something missing you would like
 to see?  Please open an issue or contact `HST Help Desk <https://stsci.service-now.com/hst>`_.
 
-``wfc3tools`` follows the `Astropy Code of Conduct`_ and strives to provide a
+``wfc3tools`` follows the STScI Code of Conduct and strives to provide a
 welcoming community to all of our users and contributors.
 
-Want more information about how to make a contribution?  Take a look at
-the the astropy `contributing`_ and `developer`_ documentation.
-
 
 License
 -------
 
 ``wfc3tools`` is licensed under a 3-clause BSD style license (see the ``LICENSE.txt`` file).
-
-.. _AstroPy: http://www.astropy.org/
-.. _contributing: http://docs.astropy.org/en/stable/index.html#contributing
-.. _developer: http://docs.astropy.org/en/stable/index.html#developer-documentation
-.. _Astropy Code of Conduct:  http://www.astropy.org/about.html#codeofconduct

wfc3tools-1.6.1.dist-info/RECORD

@@ -0,0 +1,19 @@
+wfc3tools/__init__.py,sha256=VrAe32AMZjIB3NgpMYCFJLRvGFno9lmuKI3vXMuEmck,608
+wfc3tools/calwf3.py,sha256=DeevxlnsZuVTYWgOuKf66h7n47BLlozHkOIqEVx3G9E,6765
+wfc3tools/embedsub.py,sha256=P-6mUragRH3G3V5nZMUoUiFZutuORv34R6w2Z4NX-Hk,3826
+wfc3tools/pstack.py,sha256=51vGAOhJwc1ziuNACzl-Be2WkUs_3ZFDXJDMGEKQxD0,6787
+wfc3tools/pstat.py,sha256=m3pNxW4T9RSwTfjnGdKL5ZevtTUa3hDFuU03x9Lhjhw,10097
+wfc3tools/sampinfo.py,sha256=SAJYiwsdFPnAFqb1nUZIXgJ801OJJz1QGN3pGJevLuY,6855
+wfc3tools/sub2full.py,sha256=Ab2wRam7V72_PWHyHeg3tTo3hUWeTqsu94bUMy4VsHw,5561
+wfc3tools/util.py,sha256=3tKKfWDvIVoEkFWXj6RnRr3chhdD-24PSgG-EgJegf0,1723
+wfc3tools/version.py,sha256=uJKvFol3xySBlFsy9vVg4JvWueNldnf8ba4HnPDu0xA,712
+wfc3tools/wf32d.py,sha256=ZxEYY-P-8DHfcYDM59o0td_aqbxEbrsFgXZ_BgbwjWg,5093
+wfc3tools/wf3ccd.py,sha256=2_RuVClcutEVWnl7jUN8-XaJTdV1Tj-CldhAg4e6kKg,5523
+wfc3tools/wf3cte.py,sha256=gNGDK402xyPPiYgmScWZWT4k2Q5MBOJ0EsVJqwVvFu4,2398
+wfc3tools/wf3ir.py,sha256=QLKHXnLqtwl9R_EYJ89ahQFR5SSApmwznTXcDJ4DUG8,4353
+wfc3tools/wf3rej.py,sha256=nBtWNU0H3KIjI9l_ImfoI6Ll4g2wF0Dc4uLynA2oLps,6636
+wfc3tools-1.6.1.dist-info/licenses/LICENSE.txt,sha256=psOKqD89Ruy_VuezmDEHL4M8MFGFePSOSRZ7NnuEyQM,1464
+wfc3tools-1.6.1.dist-info/METADATA,sha256=VJpR5Y6UsC0nGtyys1F1VFV5B-6upFC4vK-cSaaeQvc,2900
+wfc3tools-1.6.1.dist-info/WHEEL,sha256=wUyA8OaulRlbfwMtmQsvNngGrxQHAvkKcvRmdizlJi0,92
+wfc3tools-1.6.1.dist-info/top_level.txt,sha256=v9YyAr2gaseWI3B1PZpcZN1e8yeao5hGnM6-20IbQ50,10
+wfc3tools-1.6.1.dist-info/RECORD,,

{wfc3tools-1.6.0.dist-info → wfc3tools-1.6.1.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (80.9.0)
+Generator: setuptools (80.10.2)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

wfc3tools-1.6.0.dist-info/RECORD

@@ -1,19 +0,0 @@
-wfc3tools/__init__.py,sha256=VrAe32AMZjIB3NgpMYCFJLRvGFno9lmuKI3vXMuEmck,608
-wfc3tools/calwf3.py,sha256=B-kFOtZb5ZHgAOV1OVBE9cfuNPmDfExI8akO8kvJnDw,6794
-wfc3tools/embedsub.py,sha256=0kBeNCEua7VMj2e2cNAMOklolaL60K1id2M3btJxq4w,3883
-wfc3tools/pstack.py,sha256=6ytoFD8pBY2NAd_QuNBjkho17h85QjIndxr4fUN8LeU,6857
-wfc3tools/pstat.py,sha256=VBvFcUgGImd0LgpK7zg5cfrqgyOhL8UchgOoviLIb_g,10151
-wfc3tools/sampinfo.py,sha256=3-3ORCdN2ir26frbtTr2AMwqBO6JcU_SDgjFmBvLkjs,6970
-wfc3tools/sub2full.py,sha256=KrXPaG-XRVvUNxxof11uWYbCeV6bnDXAMfKAcuBItuA,5624
-wfc3tools/util.py,sha256=Q2vha6PKu_5P1C8K0nAj8bjrl1qLCc_91wtJn_ZVHME,1826
-wfc3tools/version.py,sha256=e6EGE-zfqfrGkBHDAkcJH8wzWg-QqsXiHlpX2I5kvxA,712
-wfc3tools/wf32d.py,sha256=ZAGdhxN4pQ8j1TOYuFB1G-ccnullykJpzb4JS80jxJM,5100
-wfc3tools/wf3ccd.py,sha256=xC91xv0cwrRWQx-y0tahOrV6nlMYrh-z3Q0TV6qVwH8,5536
-wfc3tools/wf3cte.py,sha256=CfdKjCoNdR3virbJkdv4CWmvcdEfHDGMVt-WXRs_3lo,2413
-wfc3tools/wf3ir.py,sha256=-U29B8DZTkJSeOCfg61lBdZRkXp7JFGDCCAWojVv_TA,4377
-wfc3tools/wf3rej.py,sha256=2cFzw4MBwy2uNYoWPK-OCuqHpIKTE14aNMDtJ-F8qJg,6624
-wfc3tools-1.6.0.dist-info/licenses/LICENSE.txt,sha256=psOKqD89Ruy_VuezmDEHL4M8MFGFePSOSRZ7NnuEyQM,1464
-wfc3tools-1.6.0.dist-info/METADATA,sha256=-ufUwRfHx2mazpbuAaSbJk3odWQgdZ6QP3agaLVSxr4,3252
-wfc3tools-1.6.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-wfc3tools-1.6.0.dist-info/top_level.txt,sha256=v9YyAr2gaseWI3B1PZpcZN1e8yeao5hGnM6-20IbQ50,10
-wfc3tools-1.6.0.dist-info/RECORD,,