dkist-processing-visp 3.6.3__py3-none-any.whl → 4.0.0__py3-none-any.whl

dkist_processing_visp/tasks/make_movie_frames.py

@@ -79,7 +79,7 @@ class MakeVispMovieFrames(VispTaskBase):
                                     fits_access_class=VispL1FitsAccess,
                                 )
                             )
-                            data = calibrated_frame.data
+                            data = np.nan_to_num(calibrated_frame.data, nan=0)
                             if self.constants.num_raster_steps == 1:
                                 logger.info(
                                     "Only a single raster step found. Making a spectral movie."

dkist_processing_visp/tasks/mixin/corrections.py

@@ -17,6 +17,9 @@ class CorrectionsMixin:
         arrays: Iterable[np.ndarray] | np.ndarray,
         shift: np.ndarray = np.zeros(2),
         angle: float = 0.0,
+        mode: str = "edge",
+        order: int = 5,
+        cval: float = np.nan,
     ) -> Generator[np.ndarray, None, None]:
         """
         Shift and then rotate data.
@@ -35,6 +38,24 @@ class CorrectionsMixin:
         angle : float
             The angle (in radians) between slit hairlines and pixel axes.
 
+        mode : {'constant', 'edge', 'symmetric', 'reflect', 'wrap'}
+            Points outside the boundaries of the input are filled according
+            to the given mode.  Modes match the behaviour of `numpy.pad`.
+
+        order : int
+            The order of interpolation. The order has to be in the range 0-5:
+             - 0: Nearest-neighbor
+             - 1: Bi-linear (default)
+             - 2: Bi-quadratic
+             - 3: Bi-cubic
+             - 4: Bi-quartic
+             - 5: Bi-quintic
+
+        cval : float
+            Used in conjunction with mode 'constant', the value outside
+            the image boundaries.
+
+
         Returns
         -------
         Generator
@@ -46,12 +67,21 @@ class CorrectionsMixin:
             array[np.where(array == np.inf)] = np.max(array[np.isfinite(array)])
             array[np.where(array == -np.inf)] = np.min(array[np.isfinite(array)])
             array[np.isnan(array)] = np.nanmedian(array)
-            translated = affine_transform_arrays(array, translation=-shift, mode="reflect", order=5)
-            yield next(rotate_arrays_about_point(translated, angle=-angle, mode="reflect", order=5))
+            translated = affine_transform_arrays(
+                array, translation=-shift, mode=mode, order=order, cval=cval
+            )
+            yield next(
+                rotate_arrays_about_point(
+                    translated, angle=-angle, mode=mode, order=order, cval=cval
+                )
+            )
 
     @staticmethod
     def corrections_remove_spec_geometry(
-        arrays: Iterable[np.ndarray] | np.ndarray, spec_shift: np.ndarray
+        arrays: Iterable[np.ndarray] | np.ndarray,
+        spec_shift: np.ndarray,
+        cval: float | None = None,
+        order: int = 3,
     ) -> Generator[np.ndarray, None, None]:
         """
         Remove spectral curvature.
@@ -67,6 +97,19 @@ class CorrectionsMixin:
             Array with shape (X), where X is the number of pixels in the spatial dimension.
             This dimension gives the spectral shift.
 
+        order : int
+            The order of interpolation. The order has to be in the range 0-5:
+             - 0: Nearest-neighbor
+             - 1: Bi-linear (default)
+             - 2: Bi-quadratic
+             - 3: Bi-cubic
+             - 4: Bi-quartic
+             - 5: Bi-quintic
+
+        cval : float
+            Used in conjunction with mode 'constant', the value outside
+            the image boundaries.
+
         Returns
         -------
         Generator
@@ -78,8 +121,14 @@ class CorrectionsMixin:
             numy = array.shape[1]
             array_output = np.zeros(array.shape)
             for j in range(numy):
+                if cval is None:
+                    cval = np.nanmedian(array[:, j])
                 array_output[:, j] = spnd.shift(
-                    array[:, j], -spec_shift[j], mode="constant", cval=np.nanmedian(array[:, j])
+                    array[:, j],
+                    -spec_shift[j],
+                    mode="constant",
+                    cval=cval,
+                    order=order,
                 )
             yield array_output
 

dkist_processing_visp/tasks/quality_metrics.py

@@ -162,7 +162,7 @@ class VispL1QualityMetrics(VispTaskBase, QualityMixin):
                             continue
 
                         # compute sensitivity for this Stokes parameter
-                        data_list.append(np.std(stokes_frame.data) / stokesI_med)
+                        data_list.append(np.nanstd(stokes_frame.data) / stokesI_med)
 
                 all_datetimes.append(Time(np.mean(polarization_data.datetimes), format="mjd").isot)
                 for target, source in zip(

dkist_processing_visp/tasks/science.py

@@ -51,6 +51,9 @@ class CalibrationCollection:
 
         This is done by considering that state offset values computed by the GeometricCalibration task. Any sub-pixel
         overlaps are rounded to the next integer that still guarantees overlap.
+
+        When "start pixels" are mentioned, those are pixels being counted from zero on a given axis in the positive direction.
+        When "end pixels" are mentioned, those are pixels being counted from the end of a given axis in the negative direction.
         """
         logger.info("Computing beam overlap slices")
         # This will be a flat list of (x, y) pairs for all modstates and beams
@@ -63,25 +66,33 @@ class CalibrationCollection:
         logger.info(f"All x shifts: {all_x_shifts}")
         logger.info(f"All y shifts: {all_y_shifts}")
 
-        # The amount we need to "slice in" from the front of the array is simply the maximum positive shift
+        # The amount we need to "slice in" from the start of the array is equivalent to the absolute value of the most negative shift.
         # The call to `np.ceil` ensures that the integer rounding doesn't allow non-overlap regions to leak in.
-        max_x = int(np.ceil(np.max(all_x_shifts)))
-        max_y = int(np.ceil(np.max(all_y_shifts)))
+        start_pixels_to_slice_x = int(np.ceil(abs(np.min(all_x_shifts))))
+        start_pixels_to_slice_y = int(np.ceil(abs(np.min(all_y_shifts))))
 
-        # The amount we need to "chop off" the end of the array is the most negative negative shift.
+        # The amount we need to "chop off" the end of the array is the most positive shift.
         #
         # Here we rely on the fact that the fiducial array's shift is *always* (0, 0)
-        #  (see `geometric.compute_modstate_offset`). Thus, if there are no negative shifts then the following lines
-        #  will result in None. This is required for slicing because array[x:0] is no good. So if the min is 0 then we
-        #  end up with array[x:None] which goes all the way to the end of the array.
+        # (see `geometric.compute_modstate_offset`). Thus, if there are no negative shifts then the following lines
+        # will result in None. This is required for slicing because array[x:0] is no good. So if the max is 0 then we
+        # end up with array[x:None] which goes all the way to the end of the array.
         #
-        # The call to `np.floor` ensures that the integer rounding doesn't allow non-overlap regions to leak in.
-        #  (because more negative slices will cut out more data).
-        min_x = int(np.floor(np.min(all_x_shifts))) or None
-        min_y = int(np.floor(np.min(all_y_shifts))) or None
+        # The call to `np.ceil` ensures that the integer rounding doesn't allow non-overlap regions to leak in.
+        # (because more negative slices will cut out more data).
+        end_pixels_to_slice_x = int(np.ceil(np.max(all_x_shifts))) or None
+        end_pixels_to_slice_y = int(np.ceil(np.max(all_y_shifts))) or None
+
+        # As the pixels to remove from the end of axes is given as a positive number, we need to make it negative for slicing.
+        if end_pixels_to_slice_x is not None:
+            end_pixels_to_slice_x *= -1
 
-        x_slice = slice(max_x, min_x)
-        y_slice = slice(max_y, min_y)
+        if end_pixels_to_slice_y is not None:
+            end_pixels_to_slice_y *= -1
+
+        # Construct the slices
+        x_slice = slice(start_pixels_to_slice_x, end_pixels_to_slice_x)
+        y_slice = slice(start_pixels_to_slice_y, end_pixels_to_slice_y)
 
         return x_slice, y_slice
 
@@ -277,6 +288,7 @@ class ScienceCalibration(
                 for raster_step in range(0, self.constants.num_raster_steps):
                     beam_storage = dict()
                     header_storage = dict()
+                    nan_storage = dict()
                     for beam in range(1, self.constants.num_beams + 1):
                         apm_str = f"{map_scan = }, {raster_step = }, and {beam = }"
                         with self.telemetry_span(f"Basic corrections for {apm_str}"):
@@ -288,6 +300,7 @@ class ScienceCalibration(
                                 (
                                     intermediate_array,
                                     intermediate_header,
+                                    nan_mask,
                                 ) = self.process_polarimetric_modstates(
                                     beam=beam,
                                     raster_step=raster_step,
@@ -299,7 +312,11 @@ class ScienceCalibration(
                                 logger.info(
                                     f"Processing spectrographic observe frames from {apm_str}"
                                 )
-                                intermediate_array, intermediate_header = self.correct_single_frame(
+                                (
+                                    intermediate_array,
+                                    intermediate_header,
+                                    nan_mask,
+                                ) = self.correct_single_frame(
                                     beam=beam,
                                     modstate=1,
                                     raster_step=raster_step,
@@ -310,6 +327,7 @@ class ScienceCalibration(
                                 intermediate_header = self.compute_date_keys(intermediate_header)
                             beam_storage[VispTag.beam(beam)] = intermediate_array
                             header_storage[VispTag.beam(beam)] = intermediate_header
+                            nan_storage[VispTag.beam(beam)] = nan_mask
 
                     with self.telemetry_span("Combining beams"):
                         calibrated = self.combine_beams(beam_storage, header_storage, calibrations)
@@ -318,12 +336,29 @@ class ScienceCalibration(
                         with self.telemetry_span("Correcting telescope polarization"):
                             calibrated = self.telescope_polarization_correction(calibrated)
 
+                    with self.telemetry_span("Combining NaN masks from beams"):
+                        cut_combined_nan_mask = self.combine_and_cut_nan_masks(
+                            list(nan_storage.values()), calibrations
+                        )
+
                     # Save the final output files
                     with self.telemetry_span("Writing calibrated arrays"):
                         self.write_calibrated_array(
-                            calibrated, map_scan=map_scan, calibrations=calibrations
+                            calibrated,
+                            map_scan=map_scan,
+                            calibrations=calibrations,
+                            nan_mask=cut_combined_nan_mask,
                         )
 
+    @staticmethod
+    def combine_and_cut_nan_masks(
+        nan_masks: list[np.ndarray], calibrations: CalibrationCollection
+    ) -> np.ndarray:
+        """Combine two NaN masks into one, cropping the result based on pre-calculated shifts."""
+        combined_nan_mask = np.logical_or.reduce(nan_masks)
+        x_slice, y_slice = calibrations.beams_overlap_slice
+        return combined_nan_mask[x_slice, y_slice]
+
     def process_polarimetric_modstates(
         self,
         beam: int,
@@ -331,7 +366,7 @@ class ScienceCalibration(
         map_scan: int,
         readout_exp_time: float,
         calibrations: CalibrationCollection,
-    ) -> tuple[np.ndarray, fits.Header]:
+    ) -> tuple[np.ndarray, fits.Header, np.ndarray]:
         """
         Process a single polarimetric beam as much as is possible.
 
@@ -343,11 +378,12 @@ class ScienceCalibration(
         ].shape
         array_stack = np.zeros(array_shape + (self.constants.num_modstates,))
         header_stack = []
+        nan_mask_stack = np.zeros(array_shape + (self.constants.num_modstates,))
 
         with self.telemetry_span(f"Correcting {self.constants.num_modstates} modstates"):
             for modstate in range(1, self.constants.num_modstates + 1):
                 # Correct the arrays
-                corrected_array, corrected_header = self.correct_single_frame(
+                corrected_array, corrected_header, nan_mask = self.correct_single_frame(
                     beam=beam,
                     modstate=modstate,
                     raster_step=raster_step,
@@ -358,6 +394,7 @@ class ScienceCalibration(
                 # Add this result to the 3D stack
                 array_stack[:, :, modstate - 1] = corrected_array
                 header_stack.append(corrected_header)
+                nan_mask_stack[:, :, modstate - 1] = nan_mask
 
         with self.telemetry_span("Applying instrument polarization correction"):
             intermediate_array = nd_left_matrix_multiply(
@@ -366,7 +403,8 @@ class ScienceCalibration(
             )
             intermediate_header = self.compute_date_keys(header_stack)
 
-        return intermediate_array, intermediate_header
+        # The modulator state NaN masks are stacked along axis=2 with axis=0 & 1 being the array axes of one modstate
+        return intermediate_array, intermediate_header, np.logical_or.reduce(nan_mask_stack, axis=2)
 
     def combine_beams(
         self,
@@ -440,6 +478,7 @@ class ScienceCalibration(
         calibrated_object: VispL0FitsAccess,
         map_scan: int,
         calibrations: CalibrationCollection,
+        nan_mask: np.ndarray,
     ) -> None:
         """
         Write out calibrated science frames.
@@ -458,6 +497,9 @@ class ScienceCalibration(
         calibrations
             Calibration collection
 
+        nan_mask
+            A mask containing the known areas where data does not exist for both beams
+
         Returns
         -------
         None
@@ -473,7 +515,8 @@ class ScienceCalibration(
             stokes_I_data = calibrated_object.data[:, :, 0]
             for i, stokes_param in enumerate(self.constants.stokes_params):
                 stokes_data = calibrated_object.data[:, :, i]
-                final_data = self.re_dummy_data(stokes_data)
+                nan_masked_data = np.where(nan_mask, np.nan, stokes_data)
+                final_data = self.re_dummy_data(nan_masked_data)
                 pol_header = self.add_L1_pol_headers(final_header, stokes_data, stokes_I_data)
                 self.write_cal_array(
                     data=final_data,
@@ -483,7 +526,8 @@ class ScienceCalibration(
                     map_scan=map_scan,
                 )
         else:  # Only write stokes I
-            final_data = self.re_dummy_data(calibrated_object.data)
+            nan_masked_data = np.where(nan_mask, np.nan, calibrated_object.data)
+            final_data = self.re_dummy_data(nan_masked_data)
             self.write_cal_array(
                 data=final_data,
                 header=final_header,
@@ -500,7 +544,7 @@ class ScienceCalibration(
         map_scan: int,
         readout_exp_time: float,
         calibrations: CalibrationCollection,
-    ) -> tuple[np.ndarray, fits.Header]:
+    ) -> tuple[np.ndarray, fits.Header, np.ndarray]:
         """
         Apply basic corrections to a single frame.
 
@@ -595,7 +639,40 @@ class ScienceCalibration(
             self.corrections_remove_spec_geometry(geo_corrected_array, spec_shift)
         )
 
-        return spectral_corrected_array, observe_object.header
+        nan_mask = self.generate_nan_mask(
+            solar_corrected_array=solar_corrected_array,
+            state_offset=state_offset,
+            angle=angle,
+            spec_shift=spec_shift,
+        )
+
+        return (
+            spectral_corrected_array,
+            observe_object.header,
+            nan_mask,
+        )
+
+    def generate_nan_mask(
+        self,
+        solar_corrected_array: np.ndarray,
+        state_offset: np.ndarray,
+        angle: float,
+        spec_shift: np.ndarray,
+    ) -> np.ndarray:
+        """Calculate the NaN mask through geometric correction to be applied to the final L1 arrays."""
+        # Using a bi-cubic polynomial (order = 3) best converges to the desired result in the underlying fits
+        geo_corrected_with_nan = next(
+            self.corrections_correct_geometry(
+                solar_corrected_array, state_offset, angle, mode="constant", order=3, cval=np.nan
+            )
+        )
+        # Interpolating with nearest neighbor (order = 0) prevents NaN values from "taking over" the whole array
+        spectral_corrected_with_nan = next(
+            self.corrections_remove_spec_geometry(
+                geo_corrected_with_nan, spec_shift, cval=np.nan, order=0
+            )
+        )
+        return np.isnan(spectral_corrected_with_nan)
 
     def telescope_polarization_correction(
         self,

dkist_processing_visp/tests/test_science.py

@@ -364,10 +364,10 @@ def test_science_calibration_task(
         assert header["VSPMAP"] == map_scan
 
         # Check that WCS keys were updated
-        if offsets[1, 0, 0] > 0:
-            assert header["CRPIX2"] == input_header["CRPIX2"] - np.ceil(offsets[1, 0, 0])
-        if offsets[1, 0, 1] > 0:
-            assert header["CRPIX1"] == input_header["CRPIX1"] - np.ceil(offsets[1, 0, 1])
+        if offsets[1, 0, 0] < 0:
+            assert header["CRPIX2"] == input_header["CRPIX2"] - np.ceil(-offsets[1, 0, 0])
+        if offsets[1, 0, 1] < 0:
+            assert header["CRPIX1"] == input_header["CRPIX1"] - np.ceil(-offsets[1, 0, 1])
 
     quality_files = task.read(tags=[Tag.quality("TASK_TYPES")])
     for file in quality_files:
@@ -420,7 +420,7 @@ def test_readout_normalization_correct(
     )
 
     # When:
-    corrected_array, _ = task.correct_single_frame(
+    corrected_array, _, _ = task.correct_single_frame(
         beam=1,
         modstate=1,
         raster_step=1,
@@ -513,7 +513,7 @@ def test_compute_date_keys_compressed_headers(
                     [[1.0, 2.0], [11.0, 10.0], [3.0, 2.0]],  # Beam 2
                 ]
             ),
-            [slice(11, None, None), slice(10, None, None)],
+            [slice(0, -11, None), slice(0, -10, None)],
         ),
         (
             np.array(
@@ -522,7 +522,7 @@ def test_compute_date_keys_compressed_headers(
                     [[-1.0, -2.0], [-11.0, -10.0], [-3.0, -2.0]],  # Beam 2
                 ]
             ),
-            [slice(0, -11, None), slice(0, -10, None)],
+            [slice(11, None, None), slice(10, None, None)],
         ),
         (
             np.array(
@@ -531,7 +531,7 @@ def test_compute_date_keys_compressed_headers(
                     [[1.0, 2.0], [-11.0, 10.0], [-3.0, -2.0]],  # Beam 2
                 ]
             ),
-            [slice(10, -11, None), slice(10, -2, None)],
+            [slice(11, -10, None), slice(2, -10, None)],
         ),
     ],
     ids=["All positive", "All negative", "Positive and negative"],
@@ -581,3 +581,98 @@ def test_combine_beams(
         expected = np.ones((10, 10, 4)) * 2.5
 
     np.testing.assert_array_equal(data, expected)
+
+
+@pytest.mark.parametrize(
+    "shifts",
+    # Shifts have shape (num_beams, num_modstates, 2)
+    # So the inner-most lists below (e.g., [5.0, 6.0]) correspond to [x_shift, y_shit]
+    [
+        np.array(
+            [
+                [[0.0, 0.0], [10.0, 2.0], [5.0, 6.0]],  # Beam 1
+                [[1.0, 2.0], [-11.0, 10.0], [-3.0, -2.0]],  # Beam 2
+            ]
+        ),
+    ],
+    ids=["Positive and negative"],
+)
+def test_combine_and_cut_nan_masks(
+    science_calibration_task, calibration_collection_with_geo_shifts, shifts
+):
+    """
+    Given: A ScienceCalibration task and NaN masks, along with geometric shifts
+    When: Combining the two NaN masks
+    Then: The final mask has NaN values in the correct place and is correctly cropped
+    """
+    nan_1_location = [0, 1]
+    nan_2_location = [50, 50]
+    nan_3_location = [4, 1]
+    nan_4_location = [55, 63]
+    nan_mask_shape = (100, 100)
+    nan_mask_1 = np.zeros(shape=nan_mask_shape)
+    nan_mask_1[nan_1_location[0], nan_1_location[1]] = np.nan
+    nan_mask_1[nan_2_location[0], nan_2_location[1]] = np.nan
+    nan_mask_2 = np.zeros(shape=nan_mask_shape)
+    nan_mask_2[nan_3_location[0], nan_3_location[1]] = np.nan
+    nan_mask_2[nan_4_location[0], nan_4_location[1]] = np.nan
+    task, _, _, _, _, _ = science_calibration_task
+    combined_nan_mask = task.combine_and_cut_nan_masks(
+        nan_masks=[nan_mask_1, nan_mask_2], calibrations=calibration_collection_with_geo_shifts
+    )
+    beam_1_shifts = shifts[0]
+    beam_2_shifts = shifts[1]
+    beam_1_x_shifts = [i[0] for i in beam_1_shifts]
+    beam_2_x_shifts = [i[0] for i in beam_2_shifts]
+    beam_1_y_shifts = [i[1] for i in beam_1_shifts]
+    beam_2_y_shifts = [i[1] for i in beam_2_shifts]
+    x_shifts = beam_1_x_shifts + beam_2_x_shifts
+    y_shifts = beam_1_y_shifts + beam_2_y_shifts
+    assert combined_nan_mask.shape == (
+        nan_mask_shape[0] - (max(x_shifts) - min(x_shifts)),
+        nan_mask_shape[1] - (max(y_shifts) - min(y_shifts)),
+    )
+    # Check that one NaN value from each original mask is present in the combined mask and in the correct place
+    assert (
+        combined_nan_mask[
+            nan_2_location[0] - int(abs(min(x_shifts))), nan_2_location[1] - int(abs(min(y_shifts)))
+        ]
+        == True
+    )
+    assert (
+        combined_nan_mask[
+            nan_4_location[0] - int(abs(min(x_shifts))), nan_4_location[1] - int(abs(min(y_shifts)))
+        ]
+        == True
+    )
+    assert np.sum(combined_nan_mask) == 2  # only two NaN values are in the final mask
+
+
+def test_generate_nan_mask(science_calibration_task, dummy_calibration_collection):
+    """
+    Given: a calibration collection
+    When: calculating the NaN mask to use
+    Then: the mask takes up some, but not all, of the frame size
+    """
+    task, _, _, _, _, _ = science_calibration_task
+    calibration_collection, _, _ = dummy_calibration_collection
+    beam = 1
+    modstate = 1
+    solar_gain_array = calibration_collection.solar_gain[VispTag.beam(beam)][
+        VispTag.modstate(modstate)
+    ]
+    angle = calibration_collection.angle[VispTag.beam(beam)]
+    spec_shift = calibration_collection.spec_shift[VispTag.beam(beam)]
+    state_offset = calibration_collection.state_offset[VispTag.beam(beam)][
+        VispTag.modstate(modstate)
+    ]
+    nan_mask = task.generate_nan_mask(
+        solar_corrected_array=np.random.random(size=solar_gain_array.shape),
+        state_offset=state_offset,
+        angle=angle,
+        spec_shift=spec_shift,
+    )
+    # Some of the mask is marked as NaN but not all
+    assert np.sum(nan_mask) < np.size(nan_mask)
+    # Ensure that only zeroes and ones are in the mask
+    assert set(np.unique(nan_mask)) == {0, 1}

{dkist_processing_visp-3.6.3.dist-info → dkist_processing_visp-4.0.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: dkist-processing-visp
-Version: 3.6.3
+Version: 4.0.0
 Summary: Science processing code for the ViSP instrument on DKIST
 Author-email: NSO / AURA <dkistdc@nso.edu>
 License: BSD-3-Clause
@@ -116,19 +116,19 @@ Requires-Dist: asdf==3.5.0; extra == "frozen"
 Requires-Dist: asdf_standard==1.4.0; extra == "frozen"
 Requires-Dist: asdf_transform_schemas==0.6.0; extra == "frozen"
 Requires-Dist: asgiref==3.10.0; extra == "frozen"
-Requires-Dist: asteval==1.0.6; extra == "frozen"
+Requires-Dist: asteval==1.0.7; extra == "frozen"
 Requires-Dist: astropy==7.0.2; extra == "frozen"
-Requires-Dist: astropy-iers-data==0.2025.11.3.0.38.37; extra == "frozen"
+Requires-Dist: astropy-iers-data==0.2025.11.10.0.38.31; extra == "frozen"
 Requires-Dist: asyncpg==0.30.0; extra == "frozen"
 Requires-Dist: attrs==25.4.0; extra == "frozen"
 Requires-Dist: babel==2.17.0; extra == "frozen"
 Requires-Dist: billiard==4.2.2; extra == "frozen"
 Requires-Dist: blinker==1.9.0; extra == "frozen"
-Requires-Dist: boto3==1.40.65; extra == "frozen"
-Requires-Dist: botocore==1.40.65; extra == "frozen"
+Requires-Dist: boto3==1.40.71; extra == "frozen"
+Requires-Dist: botocore==1.40.71; extra == "frozen"
 Requires-Dist: cachelib==0.13.0; extra == "frozen"
 Requires-Dist: celery==5.3.1; extra == "frozen"
-Requires-Dist: certifi==2025.10.5; extra == "frozen"
+Requires-Dist: certifi==2025.11.12; extra == "frozen"
 Requires-Dist: cffi==2.0.0; extra == "frozen"
 Requires-Dist: charset-normalizer==3.4.4; extra == "frozen"
 Requires-Dist: click==8.3.0; extra == "frozen"
@@ -152,7 +152,7 @@ Requires-Dist: dkist-processing-common==11.8.0; extra == "frozen"
 Requires-Dist: dkist-processing-core==6.0.0; extra == "frozen"
 Requires-Dist: dkist-processing-math==2.2.1; extra == "frozen"
 Requires-Dist: dkist-processing-pac==3.1.1; extra == "frozen"
-Requires-Dist: dkist-processing-visp==3.6.3; extra == "frozen"
+Requires-Dist: dkist-processing-visp==4.0.0; extra == "frozen"
 Requires-Dist: dkist-service-configuration==4.1.7; extra == "frozen"
 Requires-Dist: dkist-spectral-lines==3.0.0; extra == "frozen"
 Requires-Dist: dkist_fits_specifications==4.17.0; extra == "frozen"
@@ -164,10 +164,9 @@ Requires-Dist: fonttools==4.60.1; extra == "frozen"
 Requires-Dist: frozenlist==1.8.0; extra == "frozen"
 Requires-Dist: fsspec==2025.10.0; extra == "frozen"
 Requires-Dist: globus-sdk==3.65.0; extra == "frozen"
-Requires-Dist: google-re2==1.1.20250805; extra == "frozen"
-Requires-Dist: googleapis-common-protos==1.71.0; extra == "frozen"
+Requires-Dist: google-re2==1.1.20251105; extra == "frozen"
+Requires-Dist: googleapis-common-protos==1.72.0; extra == "frozen"
 Requires-Dist: gqlclient==1.2.3; extra == "frozen"
-Requires-Dist: greenlet==3.2.4; extra == "frozen"
 Requires-Dist: grpcio==1.76.0; extra == "frozen"
 Requires-Dist: gunicorn==23.0.0; extra == "frozen"
 Requires-Dist: h11==0.16.0; extra == "frozen"
@@ -254,9 +253,9 @@ Requires-Dist: protobuf==6.33.0; extra == "frozen"
 Requires-Dist: psutil==7.1.3; extra == "frozen"
 Requires-Dist: psycopg2-binary==2.9.11; extra == "frozen"
 Requires-Dist: pycparser==2.23; extra == "frozen"
-Requires-Dist: pydantic==2.12.3; extra == "frozen"
-Requires-Dist: pydantic-settings==2.11.0; extra == "frozen"
-Requires-Dist: pydantic_core==2.41.4; extra == "frozen"
+Requires-Dist: pydantic==2.12.4; extra == "frozen"
+Requires-Dist: pydantic-settings==2.12.0; extra == "frozen"
+Requires-Dist: pydantic_core==2.41.5; extra == "frozen"
 Requires-Dist: pyerfa==2.0.1.5; extra == "frozen"
 Requires-Dist: pyparsing==3.2.5; extra == "frozen"
 Requires-Dist: python-daemon==3.1.2; extra == "frozen"
@@ -300,7 +299,7 @@ Requires-Dist: typing_extensions==4.15.0; extra == "frozen"
 Requires-Dist: tzdata==2025.2; extra == "frozen"
 Requires-Dist: uc-micro-py==1.0.3; extra == "frozen"
 Requires-Dist: uncertainties==3.2.3; extra == "frozen"
-Requires-Dist: universal_pathlib==0.3.4; extra == "frozen"
+Requires-Dist: universal_pathlib==0.3.5; extra == "frozen"
 Requires-Dist: urllib3==2.5.0; extra == "frozen"
 Requires-Dist: vine==5.1.0; extra == "frozen"
 Requires-Dist: voluptuous==0.15.2; extra == "frozen"

{dkist_processing_visp-3.6.3.dist-info → dkist_processing_visp-4.0.0.dist-info}/RECORD

@@ -24,16 +24,16 @@ dkist_processing_visp/tasks/geometric.py,sha256=rvWa13T2_cPQtiMjlWrqMYS5UvL8Xgmi
 dkist_processing_visp/tasks/instrument_polarization.py,sha256=uj7iyzM3CiJcbQeF4eKpk_KCoheXaM4FpDI83GYDld4,25854
 dkist_processing_visp/tasks/l1_output_data.py,sha256=lon6bIUBvURV_7gn1HFGAGVuEiLoOiMaJwfbdhSRW3w,459
 dkist_processing_visp/tasks/lamp.py,sha256=RciNB8zW5fDx5yrsaNZgFP7LsajmSvkTNvbe3Sm6tjc,6060
-dkist_processing_visp/tasks/make_movie_frames.py,sha256=v-i-PqDe2ZCYS4oe8LeEcM19g77tQAOQUcOs8X1Jz8k,7522
+dkist_processing_visp/tasks/make_movie_frames.py,sha256=fw25ksKiJJNS57XV5a7rHpYGcSkYxS2Qf13Fb1UGNpE,7544
 dkist_processing_visp/tasks/parse.py,sha256=3f8LWSKQtuY7n-PdHiw1j55i5l4juBn3pIHqTtp2-rk,6752
-dkist_processing_visp/tasks/quality_metrics.py,sha256=Xz0ft1Y6ocJiVIDKe2fp8YxWsy3ctuGvYI-4hWFL-P8,8130
-dkist_processing_visp/tasks/science.py,sha256=HO1Oeyy_31z50dNMgaV0gbaHEptJy9k0TiGbkgexhe8,30961
+dkist_processing_visp/tasks/quality_metrics.py,sha256=Pw55-PXW0cl39FuNkEQCGGhvI_zMDimwmh-swVPVBD4,8133
+dkist_processing_visp/tasks/science.py,sha256=BiEtgDowExHAXXzfO4BSr34BT0PZ28CdtOJMP4TDA20,34680
 dkist_processing_visp/tasks/solar.py,sha256=QPkrJqUz6gYxFFvNlATMKa4NIEdn4iassYn86qPb60k,27567
 dkist_processing_visp/tasks/visp_base.py,sha256=flUM-dCWkjcGvfaiFWH89DRZXz7NADgShxzpFjFqVEw,1370
 dkist_processing_visp/tasks/write_l1.py,sha256=bsDZ0BwoqpTtS_f_rAzUn7Ra8UvYb-kINQhX6BwwFQw,8796
 dkist_processing_visp/tasks/mixin/__init__.py,sha256=z2nFVvvIzirxklQ9i5-F1nR-WOgcDttYtog_jx4yN5I,12
 dkist_processing_visp/tasks/mixin/beam_access.py,sha256=1VSJkH6yMxCiZWdWOp_RJ37fX5ULMYmB_0_ulT7YJpI,870
-dkist_processing_visp/tasks/mixin/corrections.py,sha256=tsXexnqsYgd14CjGYrsBxGjX7ZTUVrTj6_-2VdKeWP0,5666
+dkist_processing_visp/tasks/mixin/corrections.py,sha256=FhLFgD9ZYLZd3SaC3PFF-szrcs-zmdrUYNDUEK-h7JA,7145
 dkist_processing_visp/tasks/mixin/downsample.py,sha256=SvKzY6HJRn-FeyG7O6HPvyOS5dmMu6uPoWkfnpPXpVw,1344
 dkist_processing_visp/tasks/mixin/line_zones.py,sha256=5jfea9V5RJAi-834z_Y9v4fhlRFJdK1McAqO9X92bZo,4065
 dkist_processing_visp/tests/README.rst,sha256=rnedwwg25c0lB9Me7cT7QNZA17FYlqCu9ZnjQxR5hi0,12502
@@ -54,7 +54,7 @@ dkist_processing_visp/tests/test_map_repeats.py,sha256=9g3NnvSfn1OqxxYYxTFoOIi1U
 dkist_processing_visp/tests/test_parameters.py,sha256=h9EemGJf0b4ma0jLGd321untkhLkhwgo88mIRnSmxXs,4369
 dkist_processing_visp/tests/test_parse.py,sha256=aFBbLBXdz2KJZy62gC1KC049FFfRKYqYVo-Y32mx_E4,20853
 dkist_processing_visp/tests/test_quality.py,sha256=YW24VjEHoILseFIXZBp4-o7egT26mfT1lafzajVjXu8,6905
-dkist_processing_visp/tests/test_science.py,sha256=o9bnm8fZmzRYk5u0c9GgeWPyGAZu2_sbN52pPpvrJUo,20272
+dkist_processing_visp/tests/test_science.py,sha256=rIoHu1hifNYj1niQaDkmyID5BQVzTEgZ5BhsdZOSaWk,23985
 dkist_processing_visp/tests/test_solar.py,sha256=Bpkm59HI0GqrwM8LqIK7GR0pUUW0anboJkePQnofgmA,9471
 dkist_processing_visp/tests/test_trial_create_quality_report.py,sha256=t3de9LMJOqrRaFXAvKV_5sotAfzDR8fZVIrFNRB2I_A,2779
 dkist_processing_visp/tests/test_visp_constants.py,sha256=v8meAXdR6jJntaFtWR-N65Z-3ovr9Vf9_-BwAZv2RBc,1974
@@ -87,7 +87,7 @@ docs/requirements_table.rst,sha256=_HIbwFpDooM5n0JjiDAbFozGfJuX13smtcoujLFN4Gk,2
 docs/science_calibration.rst,sha256=VN_g7xSjN-nbXhlBaFnPCbNcsc_Qu0207jEUfRAjnBE,2939
 docs/scientific_changelog.rst,sha256=01AWBSHg8zElnodCgAq-hMxhk9CkX5rtEENx4iz0sjI,300
 licenses/LICENSE.rst,sha256=piZaQplkzOMmH1NXg6QIdo9wwo9pPCoHkvm2-DmH76E,1462
-dkist_processing_visp-3.6.3.dist-info/METADATA,sha256=OuZTRuLp8xx5ZWIiVWD6TgdIuRYeBd175z0JtZP8jdw,29121
-dkist_processing_visp-3.6.3.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-dkist_processing_visp-3.6.3.dist-info/top_level.txt,sha256=9GHSn-ZMGQxaRNGrPP3HNc5ZkE7ftzluO74Jz5vUSTE,46
-dkist_processing_visp-3.6.3.dist-info/RECORD,,
+dkist_processing_visp-4.0.0.dist-info/METADATA,sha256=4S8qtlI-Cp12PpvTDgD0ZgSA2ZaaOp0QnHtC94RZ4Ak,29073
+dkist_processing_visp-4.0.0.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+dkist_processing_visp-4.0.0.dist-info/top_level.txt,sha256=9GHSn-ZMGQxaRNGrPP3HNc5ZkE7ftzluO74Jz5vUSTE,46
+dkist_processing_visp-4.0.0.dist-info/RECORD,,