teareduce 0.4.5__py3-none-any.whl → 0.4.7__py3-none-any.whl

teareduce/ctext.py

@@ -1,5 +1,5 @@
 #
-# Copyright 2022-2024 Universidad Complutense de Madrid
+# Copyright 2022-2025 Universidad Complutense de Madrid
 #
 # This file is part of teareduce
 #
@@ -9,8 +9,8 @@
 
 
 def ctext(s=None,
-          fg=None, 
-          bg=None, 
+          fg=None,
+          bg=None,
           under=False,
           rev=False,
           faint=False,
@@ -95,4 +95,3 @@ def ctext(s=None,
     result = f'{final_style}{s}\x1B[0m'
 
     return result
-

teareduce/draw_rectangle.py

@@ -1,5 +1,5 @@
 #
-# Copyright 2022-2024 Universidad Complutense de Madrid
+# Copyright 2022-2025 Universidad Complutense de Madrid
 #
 # This file is part of teareduce
 #
@@ -56,11 +56,11 @@ def draw_rectangle(ax, image_data, x1, x2, y1, y2,
     ax.plot((x1, x2), (y2, y2), color, lw=1)
 
     if text:
-        ax.text((x1+x2)/2, y1+(y2-y1)/8, 
-                '{:.{prec}f}'.format(mean, prec=ndigits), 
-                ha='center', va='center', color=color, fontsize=fontsize)        
+        ax.text((x1+x2)/2, y1+(y2-y1)/8,
+                '{:.{prec}f}'.format(mean, prec=ndigits),
+                ha='center', va='center', color=color, fontsize=fontsize)
         ax.text((x1+x2)/2, y2-(y2-y1)/8,
-                '{:.{prec}f}'.format(std, prec=ndigits), 
+                '{:.{prec}f}'.format(std, prec=ndigits),
                 ha='center', va='top', color=color, fontsize=fontsize)
 
     return mean, std

teareduce/imshow.py

@@ -41,7 +41,7 @@ def imshowme(data, **kwargs):
     return fig, ax, img, cax, cbar
 
 
-def imshow(fig=None, ax=None, data=None,
+def imshow(fig=None, ax=None, data=None, ds9mode=False,
            crpix1=1, crval1=None, cdelt1=None, cunit1=None, cunitx=Unit('Angstrom'),
            xlabel=None, ylabel=None, title=None,
            colorbar=True, cblabel='Number of counts',
@@ -60,6 +60,10 @@ def imshow(fig=None, ax=None, data=None,
         Instance of Axes.
     data : numpy array
         2D array to be displayed.
+    ds9mode : bool
+        If True, the extent parameter is set to
+        [0.5, NAXIS1+0.5, 0.5, NAXIS2+0.5]
+        to mimic the DS9 display.
     crpix1 : astropy.units.Quantity
         Float number providing the CRPIX1 value: the reference pixel
         for which CRVAL1 is given.
@@ -103,13 +107,6 @@ def imshow(fig=None, ax=None, data=None,
     if not isinstance(ax, Axes):
         raise ValueError("Unexpected 'ax' argument")
 
-    # default labels
-    if xlabel is None:
-        xlabel = 'X axis (array index)'
-
-    if ylabel is None:
-        ylabel = 'Y axis (array index)'
-
     wavecalib = False
     if crpix1 is not None and crval1 is not None and cdelt1 is not None and cunit1 is not None:
         if 'extent' in kwargs:
@@ -129,11 +126,25 @@ def imshow(fig=None, ax=None, data=None,
         aspect = 'auto'
         wavecalib = True
     else:
-        if 'extent' in kwargs:
-            extent = kwargs['extent']
-            del kwargs['extent']
+        if ds9mode:
+            if 'extent' in kwargs:
+                raise ValueError('extent parameter can not be used with ds9mode=True')
+            naxis2, naxis1 = data.shape
+            extent = [0.5, naxis1 + 0.5, 0.5, naxis2 + 0.5]
+            if xlabel is None:
+                xlabel = 'X pixel (from 1 to NAXIS1)'
+            if ylabel is None:
+                ylabel = 'Y pixel (from 1 to NAXIS2)'
         else:
-            extent = None
+            if xlabel is None:
+                xlabel = 'X axis (array index)'
+            if ylabel is None:
+                ylabel = 'Y axis (array index)'
+            if 'extent' in kwargs:
+                extent = kwargs['extent']
+                del kwargs['extent']
+            else:
+                extent = None
         if 'aspect' in kwargs:
             aspect = kwargs['aspect']
             del kwargs['aspect']

teareduce/peaks_spectrum.py

@@ -1,6 +1,6 @@
 # -*- coding: utf-8 -*-
 #
-# Copyright 2015-2024 Universidad Complutense de Madrid
+# Copyright 2015-2025 Universidad Complutense de Madrid
 #
 # This file is part of teareduce
 #
@@ -58,12 +58,12 @@ def find_peaks_spectrum(sx, nwinwidth, deltaflux=0, threshold=0, debugplot=False
               'pixels will be ignored')
 
     xpeaks = []  # list to store the peaks
-    
+
     if sx_shape[0] < nwinwidth:
         print('find_peaks_spectrum> sx shape......:', sx_shape)
         print('find_peaks_spectrum> nwinwidth.....:', nwinwidth)
         raise ValueError('sx.shape < nwinwidth')
-    
+
     i = nmed
     while i < sx_shape[0] - nmed:
         if sx[i] > threshold:
@@ -93,7 +93,7 @@ def find_peaks_spectrum(sx, nwinwidth, deltaflux=0, threshold=0, debugplot=False
                 i += 1
         else:
             i += 1
-    
+
     ixpeaks = np.array(xpeaks)
 
     if debugplot:

teareduce/sliceregion.py

@@ -8,7 +8,7 @@
 #
 """Auxiliary classes to handle slicing regions in 1D, 2D, and 3D.
 
-These classes provide a way to define and manipulate slices in a 
+These classes provide a way to define and manipulate slices in a
 consistent manner, following both FITS and Python conventions.
 """
 
@@ -16,9 +16,10 @@ import re
 
 import numpy as np
 
+
 class SliceRegion1D:
     """Store indices for slicing of 1D regions.
-    
+
     The attributes .python and .fits provide the indices following
     the Python and the FITS convention, respectively.
 
@@ -42,7 +43,7 @@ class SliceRegion1D:
         Check if slice 'other' is within the parent slice.
     """
 
-    def __init__(self, region, mode=None):
+    def __init__(self, region, mode=None, naxis1=None):
         """Initialize SliceRegion1D.
 
         Parameters
@@ -53,6 +54,9 @@ class SliceRegion1D:
         mode : str
             Convention mode employed to define the slice.
             The two possible modes are 'fits' and 'python'.
+        naxis1 : int
+            The axis 1 size (length) of the data being sliced.
+            If provided, it is used to validate the slice region.
         """
         if isinstance(region, str):
             pattern = r'^\s*\[\s*\d+\s*:\s*\d+\s*\]\s*$'
@@ -92,6 +96,12 @@ class SliceRegion1D:
         else:
             raise ValueError(errmsg)
 
+        if naxis1 is not None:
+            if not (1 <= self.fits.start <= naxis1):
+                raise ValueError(f'Invalid start={self.fits.start} for naxis1={naxis1}')
+            if not (1 <= self.fits.stop <= naxis1+1):
+                raise ValueError(f'Invalid stop={self.fits.stop} for naxis1={naxis1}')
+
         s = self.fits
         self.fits_section = f'[{s.start}:{s.stop}]'
 
@@ -136,6 +146,10 @@ class SliceRegion1D:
         result = True
         return result
 
+    def length(self):
+        """Return the length of the slice."""
+        return self.python.stop - self.python.start
+
 
 class SliceRegion2D:
     """Store indices for slicing of 2D regions.
@@ -162,7 +176,7 @@ class SliceRegion2D:
     within(other)
         Check if slice 'other' is within the parent slice."""
 
-    def __init__(self, region, mode=None):
+    def __init__(self, region, mode=None, naxis1=None, naxis2=None):
         """Initialize SliceRegion2D.
 
         Parameters
@@ -174,6 +188,14 @@ class SliceRegion2D:
         mode : str
             Convention mode employed to define the slice.
             The two possible modes are 'fits' and 'python'.
+        naxis1 : int
+            The axis 1 size (length) of the data being sliced,
+            assuming the FITS convention.
+            If provided, it is used to validate the slice region.
+        naxis2 : int
+            The axis 2 size (length) of the data being sliced,
+            assuming the FITS convention.
+            If provided, it is used to validate the slice region.
         """
         if isinstance(region, str):
             pattern = r'^\s*\[\s*\d+\s*:\s*\d+\s*,\s*\d+\s*:\s*\d+\s*\]\s*$'
@@ -223,6 +245,17 @@ class SliceRegion2D:
         s1, s2 = self.fits
         self.fits_section = f'[{s1.start}:{s1.stop},{s2.start}:{s2.stop}]'
 
+        if naxis1 is not None:
+            if not (1 <= s1.start <= naxis1):
+                raise ValueError(f'Invalid start={s1.start} for naxis1={naxis1}')
+            if not (1 <= s1.stop <= naxis1):
+                raise ValueError(f'Invalid stop={s1.stop} for naxis1={naxis1}')
+        if naxis2 is not None:
+            if not (1 <= s2.start <= naxis2):
+                raise ValueError(f'Invalid start={s2.start} for naxis2={naxis2}')
+            if not (1 <= s2.stop <= naxis2):
+                raise ValueError(f'Invalid stop={s2.stop} for naxis2={naxis2}')
+
     def __eq__(self, other):
         return self.fits == other.fits and self.python == other.python
 
@@ -268,6 +301,11 @@ class SliceRegion2D:
         result = True
         return result
 
+    def area(self):
+        """Return the area of the slice."""
+        s1, s2 = self.python
+        return (s1.stop - s1.start) * (s2.stop - s2.start)
+
 
 class SliceRegion3D:
     """Store indices for slicing of 3D regions.
@@ -294,7 +332,7 @@ class SliceRegion3D:
     within(other)
         Check if slice 'other' is within the parent slice."""
 
-    def __init__(self, region, mode=None):
+    def __init__(self, region, mode=None, naxis1=None, naxis2=None, naxis3=None):
         """Initialize SliceRegion3D.
 
         Parameters
@@ -306,6 +344,18 @@ class SliceRegion3D:
         mode : str
             Convention mode employed to define the slice.
             The two possible modes are 'fits' and 'python'.
+        naxis1 : int
+            The axis 1 size (length) of the data being sliced,
+            assuming the FITS convention.
+            If provided, it is used to validate the slice region.
+        naxis2 : int
+            The axis 2 size (length) of the data being sliced,
+            assuming the FITS convention.
+            If provided, it is used to validate the slice region.
+        naxis3 : int
+            The axis 3 size (length) of the data being sliced,
+            assuming the FITS convention.
+            If provided, it is used to validate the slice region.
         """
         if isinstance(region, str):
             pattern = r'^\s*\[\s*\d+\s*:\s*\d+\s*,\s*\d+\s*:\s*\d+\s*,\s*\d+\s*:\s*\d+\s*\]\s*$'
@@ -359,6 +409,22 @@ class SliceRegion3D:
         s1, s2, s3 = self.fits
         self.fits_section = f'[{s1.start}:{s1.stop},{s2.start}:{s2.stop},{s3.start}:{s3.stop}]'
 
+        if naxis1 is not None:
+            if not (1 <= s1.start <= naxis1):
+                raise ValueError(f'Invalid start={s1.start} for naxis1={naxis1}')
+            if not (1 <= s1.stop <= naxis1):
+                raise ValueError(f'Invalid stop={s1.stop} for naxis1={naxis1}')
+        if naxis2 is not None:
+            if not (1 <= s2.start <= naxis2):
+                raise ValueError(f'Invalid start={s2.start} for naxis2={naxis2}')
+            if not (1 <= s2.stop <= naxis2):
+                raise ValueError(f'Invalid stop={s2.stop} for naxis2={naxis2}')
+        if naxis3 is not None:
+            if not (1 <= s3.start <= naxis3):
+                raise ValueError(f'Invalid start={s3.start} for naxis3={naxis3}')
+            if not (1 <= s3.stop <= naxis3):
+                raise ValueError(f'Invalid stop={s3.stop} for naxis3={naxis3}')
+
     def __eq__(self, other):
         return self.fits == other.fits and self.python == other.python
 
@@ -407,3 +473,8 @@ class SliceRegion3D:
             return result
         result = True
         return result
+
+    def volume(self):
+        """Return the volume of the slice."""
+        s1, s2, s3 = self.python
+        return (s1.stop - s1.start) * (s2.stop - s2.start) * (s3.stop - s3.start)

teareduce/statsummary.py

@@ -1,5 +1,5 @@
 #
-# Copyright 2022-2024 Universidad Complutense de Madrid
+# Copyright 2022-2025 Universidad Complutense de Madrid
 #
 # This file is part of teareduce
 #
@@ -22,7 +22,7 @@ def statsummary(x=None, rm_nan=False, show=True):
     Parameters
     ----------
     x : numpy array or None
-        Input array with values which statistical properties are 
+        Input array with values which statistical properties are
         requested.
     rm_nan : bool
         If True, filter out NaN values before computing statistics.
@@ -42,9 +42,9 @@ def statsummary(x=None, rm_nan=False, show=True):
 
     # protections
     if x is None:
-        return ['npoints', 'minimum', 'maximum', 
-                'mean', 'median', 'std', 'robust_std', 
-                'percentile16', 'percentile25', 'percentile75', 'percentile84'] 
+        return ['npoints', 'minimum', 'maximum',
+                'mean', 'median', 'std', 'robust_std',
+                'percentile16', 'percentile25', 'percentile75', 'percentile84']
 
     if isinstance(x, np.ndarray):
         xx = np.copy(x.flatten())

teareduce/tests/test_sliceregion.py

@@ -32,6 +32,7 @@ def test_slice_region_creation():
     assert region3d.python == (slice(1, 10, None), slice(2, 20, None), slice(3, 30, None))
     assert region3d.fits_section == '[4:30,3:20,2:10]'
 
+
 def test_slice_values():
     """Test the values of the slices in different modes."""
 
@@ -47,3 +48,69 @@ def test_slice_values():
     array3d = np.arange(24).reshape(3, 4, 2)
     region3d = SliceRegion3D(np.s_[1:3, 2:4, 1:2], mode='python')
     assert np.all(array3d[region3d.python] == np.array([[[13], [15]], [[21], [23]]]))
+
+
+def test_wrong_number_of_dimensions():
+    """Test the creation of a SliceRegion with wrong slices."""
+
+    try:
+        SliceRegion1D(np.s_[1:3, 2:4], mode='python')
+        assert False, "Expected ValueError for 1D slice with 2D input"
+    except ValueError:
+        pass
+
+    try:
+        SliceRegion2D(np.s_[1:3], mode='python')
+        assert False, "Expected ValueError for 2D slice with 1D input"
+    except ValueError:
+        pass
+
+    try:
+        SliceRegion3D(np.s_[1:3, 2:4], mode='python')
+        assert False, "Expected ValueError for 3D slice with 2D input"
+    except ValueError:
+        pass
+
+
+def test_wrong_limits_order():
+    """Test the creation of a SliceRegion with wrong limits."""
+
+    try:
+        SliceRegion1D(np.s_[10:5], mode='python')
+        assert False, "Expected ValueError for 1D slice with start > stop"
+    except ValueError:
+        pass
+
+    try:
+        SliceRegion2D(np.s_[1:3, 5:2], mode='python')
+        assert False, "Expected ValueError for 2D slice with start > stop in second dimension"
+    except ValueError:
+        pass
+
+    try:
+        SliceRegion3D(np.s_[1:3, 2:4, 6:1], mode='python')
+        assert False, "Expected ValueError for 3D slice with start > stop in third dimension"
+    except ValueError:
+        pass
+
+
+def test_limits_out_of_range():
+    """Test the creation of a SliceRegion with limits out of range."""
+
+    try:
+        SliceRegion1D(np.s_[-1:5], mode='python', naxis1=10)
+        assert False, "Expected ValueError for 1D slice with negative start"
+    except ValueError:
+        pass
+
+    try:
+        SliceRegion2D(np.s_[1:3, 2:25], mode='python', naxis1=10, naxis2=20)
+        assert False, "Expected ValueError for 2D slice with stop > naxis in second dimension"
+    except ValueError:
+        pass
+
+    try:
+        SliceRegion3D(np.s_[1:3, 2:4, 3:35], mode='python', naxis1=10, naxis2=20, naxis3=30)
+        assert False, "Expected ValueError for 3D slice with stop > naxis in third dimension"
+    except ValueError:
+        pass

teareduce/version.py

@@ -9,7 +9,7 @@
 #
 """Module to define the version of the teareduce package."""
 
-VERSION = '0.4.5'
+VERSION = '0.4.7'
 
 
 def main():

teareduce/wavecal.py

@@ -1,6 +1,6 @@
 # -*- coding: utf-8 -*-
 #
-# Copyright 2015-2024 Universidad Complutense de Madrid
+# Copyright 2015-2025 Universidad Complutense de Madrid
 #
 # This file is part of teareduce
 #
@@ -30,10 +30,10 @@ from .sliceregion import SliceRegion1D
 
 class TeaWaveCalibration:
     """Auxiliary class to compute and apply wavelength calibration.
-    
+
     It is assumed that the wavelength and spatial directions correspond
     to the X (array columns) and Y (array rows) axes, respectively.
-    
+
     Attributes
     ----------
     ns_window : int
@@ -44,7 +44,7 @@ class TeaWaveCalibration:
         Standard deviation for Gaussian kernel to smooth median spectrum.
         A value of 0 means that no smoothing is performed.
     nx_window : int
-        Number of pixels (spectral direction) of the window where the 
+        Number of pixels (spectral direction) of the window where the
         peaks are sought. It must be odd.
     delta_flux : float
         Minimum difference between the flux at the line center and the
@@ -293,7 +293,7 @@ class TeaWaveCalibration:
             Standard deviation for Gaussian kernel to smooth median spectrum.
             A value of 0 means that no smoothing is performed.
         nx_window : int or None
-            Number of pixels (spectral direction) of the window where the 
+            Number of pixels (spectral direction) of the window where the
             peaks are sought. It must be odd.
         delta_flux : float
             Minimum difference between the flux at the line center and the
@@ -335,7 +335,7 @@ class TeaWaveCalibration:
             self.delta_flux = delta_flux
         if method is not None:
             self.method = method
-        
+
         naxis2, naxis1 = data.shape
 
         if self._naxis1 is None:
@@ -343,7 +343,7 @@ class TeaWaveCalibration:
         else:
             if naxis1 != self._naxis1:
                 raise ValueError(f'Unexpected naxis1: {naxis1}')
-        
+
         if self._naxis2 is None:
             self._naxis2 = naxis2
         else:
@@ -363,7 +363,7 @@ class TeaWaveCalibration:
 
         # initial median spectrum
         xpeaks, ixpeaks, sp_median_smooth = self._find_peaks_scan(
-            data=data, 
+            data=data,
             ns1=ns1,
             ns2=ns2,
             plot_peaks=plot_peaks,
@@ -420,12 +420,12 @@ class TeaWaveCalibration:
                              xpeaks_reference=None,
                              ns_range=None,
                              direction='up',
-                             ns_window = None,
-                             threshold = None,
-                             sigma_smooth = None,
-                             nx_window = None,
-                             delta_flux = None,
-                             method = None,
+                             ns_window=None,
+                             threshold=None,
+                             sigma_smooth=None,
+                             nx_window=None,
+                             delta_flux=None,
+                             method=None,
                              plots=False,
                              title=None,
                              pdf_output=None,
@@ -455,7 +455,7 @@ class TeaWaveCalibration:
             Standard deviation for Gaussian kernel to smooth median spectrum.
             A value of 0 means that no smoothing is performed.
         nx_window : int or None
-            Number of pixels (spectral direction) of the window where the 
+            Number of pixels (spectral direction) of the window where the
             peaks are sought. It must be odd.
         delta_flux : float
             Minimum difference between the flux at the line center and the
@@ -540,7 +540,7 @@ class TeaWaveCalibration:
             color_previous = ['blue', 'cyan']
             fig, ax = plt.subplots(figsize=(15, 15*naxis2/naxis1))
             vmin, vmax = np.percentile(data, [5, 95])
-            img = imshow(fig, ax, data, vmin=vmin, vmax=vmax, cmap='gray', title=title, aspect='auto')
+            imshow(fig, ax, data, vmin=vmin, vmax=vmax, cmap='gray', title=title, aspect='auto')
             # display previously identified lines
             yplot = np.arange(naxis2)[self._valid_scans]
             for i in range(self._nlines_reference):
@@ -549,18 +549,17 @@ class TeaWaveCalibration:
         else:
             fig = None
             ax = None
-            img = None
 
         # search for peaks in the 2D image
         dict_xpeaks = dict()
-        for ns in tqdm(range(ns_min_fits, ns_max_fits + ns_step, ns_step), 
+        for ns in tqdm(range(ns_min_fits, ns_max_fits + ns_step, ns_step),
                        desc='Finding peaks', disable=disable_tqdm):
             ns1 = ns - self.ns_window // 2
             ns1 = max([ns1, min(ns_min_fits, ns_max_fits)])
             ns2 = ns + self.ns_window // 2
             ns2 = min([ns2, max(ns_min_fits, ns_max_fits)])
             xpeaks, ixpeaks, sp_median_smooth = self._find_peaks_scan(
-                data=data, 
+                data=data,
                 ns1=ns1,
                 ns2=ns2,
                 plot_peaks=False,
@@ -585,7 +584,7 @@ class TeaWaveCalibration:
                 xpeaks_predicted = np.median(self._xpeaks_all_lines_array[(ns1-1):ns2, :], axis=0)
             for i in range(self._nlines_reference):
                 value = xpeaks_predicted[i]
-                # if there is no peak near to the expected location 
+                # if there is no peak near to the expected location
                 # (within a distance given by self.nx_window) the line
                 # is probably weak and we need to avoid jumping into
                 # another line
@@ -661,7 +660,7 @@ class TeaWaveCalibration:
         self.peak_wavelengths = wavelengths
 
     @u.quantity_input(xpeaks=u.pixel)
-    def overplot_identified_lines(self, xpeaks, spectrum, 
+    def overplot_identified_lines(self, xpeaks, spectrum,
                                   title=None, fontsize_title=16, fontsize_wave=10,
                                   pdf_output=None, pdf_only=False):
         """Overplot identified lines
@@ -790,7 +789,7 @@ class TeaWaveCalibration:
         xfit = np.arange(self._naxis2)[self._valid_scans]
         if len(xfit) <= self.degree_cdistortion:
             raise ValueError(f'Insufficient number of points to fit a polynomial of degree {self.degree_cdistortion}')
-        for i in tqdm(range(self._nlines_reference), 
+        for i in tqdm(range(self._nlines_reference),
                       desc='Fitting C distortion', disable=disable_tqdm):
             yfit = self._xpeaks_all_lines_array[self._valid_scans, i]
             poly, yres, reject = polfit_residuals_with_sigma_rejection(
@@ -1124,11 +1123,11 @@ class TeaWaveCalibration:
                 raise ValueError('You must set plots=True to make use of pdf_output')
 
         return poly_fits_yx, residual_std_yx, \
-               poly_fits_xy, residual_std_xy, \
-               crval1_linear, cdelt1_linear, crmax1_linear
+            poly_fits_xy, residual_std_xy, \
+            crval1_linear, cdelt1_linear, crmax1_linear
 
     def fit_wavelengths(self, degree_wavecalib=None,
-                        output_filename=None, history_list=None, 
+                        output_filename=None, history_list=None,
                         plots=False, title=None,
                         pdf_output=None, pdf_only=False,
                         silent_mode=False, disable_tqdm=True):
@@ -1320,9 +1319,9 @@ class TeaWaveCalibration:
         # crval1, cdelt1, crmax1, residual_std
         fig, axarr = plt.subplots(nrows=1, ncols=5, figsize=(15, 3))
         axarr = axarr.flatten()
-        for i, item in enumerate(['_array_crval1_linear', 
+        for i, item in enumerate(['_array_crval1_linear',
                                   '_array_cdelt1_linear',
-                                  '_array_crmax1_linear', 
+                                  '_array_crmax1_linear',
                                   '_array_residual_std_wav',
                                   '_array_residual_std_pix']):
             ax = axarr[i]
@@ -1357,8 +1356,7 @@ class TeaWaveCalibration:
         nrows = int(ncoeff / npprow)
         if ncoeff % npprow != 0:
             nrows += 1
-        fig, axarr = plt.subplots(nrows=nrows, ncols=npprow, 
-                                  figsize=(figwidth, 3*nrows))
+        fig, axarr = plt.subplots(nrows=nrows, ncols=npprow, figsize=(figwidth, 3*nrows))
         axarr = axarr.flatten()
         for ax in axarr:
             ax.axis('off')
@@ -1435,7 +1433,7 @@ class TeaWaveCalibration:
 
         old_x_borders_fits = np.arange(naxis1 + 1) + 0.5  # FITS convention
 
-        for k in tqdm(range(naxis2), 
+        for k in tqdm(range(naxis2),
                       desc='Applying wavelength calibration',
                       disable=disable_tqdm):
             poly = Polynomial(self._array_poly_wav[k])
@@ -1678,4 +1676,3 @@ def apply_wavecal_ccddata(infile, wcalibfile, outfile,
             cdelt1=cdelt1,
             title=f'{title_}(UNCERT extension)'
         )
-

teareduce/zscale.py

@@ -1,5 +1,5 @@
 #
-# Copyright 2022-2024 Universidad Complutense de Madrid
+# Copyright 2022-2025 Universidad Complutense de Madrid
 #
 # This file is part of teareduce
 #
@@ -12,8 +12,8 @@ import numpy as np
 
 def zscale(image, factor=0.25):
     """Compute z1 and z2 cuts in a similar way to Iraf.
-    
-    If the total number of pixels is less than 10, the function simply 
+
+    If the total number of pixels is less than 10, the function simply
     returns the minimum and the maximum values.
 
     Parameters
@@ -48,5 +48,5 @@ def zscale(image, factor=0.25):
         z1 = max(z1, q000)
         z2 = q500+(zslope*npixels/2)/factor
         z2 = min(z2, q1000)
-    
+
     return z1, z2

{teareduce-0.4.5.dist-info → teareduce-0.4.7.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: teareduce
-Version: 0.4.5
+Version: 0.4.7
 Summary: Utilities for astronomical data reduction
 Author-email: Nicolás Cardiel <cardiel@ucm.es>
 License: GPL-3.0-or-later