setiastrosuitepro 1.6.4__py3-none-any.whl → 1.7.1.post2__py3-none-any.whl

setiastro/saspro/star_alignment.py

@@ -631,7 +631,10 @@ class StellarAlignmentDialog(QDialog):
         self.setWindowFlag(Qt.WindowType.Window, True)
         self.setWindowModality(Qt.WindowModality.NonModal)
         self.setModal(False)
-        #self.setAttribute(Qt.WidgetAttribute.WA_DeleteOnClose, True)
+        try:
+            self.setAttribute(Qt.WidgetAttribute.WA_DeleteOnClose, True)
+        except Exception:
+            pass  # older PyQt6 versions
 
         self.settings = settings
         self.parent_window = parent
@@ -4583,11 +4586,12 @@ def load_api_key():
 class MosaicMasterDialog(QDialog):
     def __init__(self, settings: QSettings, parent=None, image_manager=None,
                  doc_manager=None, wrench_path=None, spinner_path=None,
-                 list_open_docs_fn=None):                                 # ← add param
+                 list_open_docs_fn=None):
         super().__init__(parent)
         self.settings = settings
         self.image_manager = image_manager
         self._docman = doc_manager or getattr(parent, "doc_manager", None)
+
         # same pattern as StellarAlignmentDialog
         if list_open_docs_fn is None:
             cand = getattr(parent, "_list_open_docs", None)
@@ -4599,19 +4603,22 @@ class MosaicMasterDialog(QDialog):
         self.setWindowFlag(Qt.WindowType.Window, True)
         self.setWindowModality(Qt.WindowModality.NonModal)
         self.setModal(False)
+
         self.wrench_path = wrench_path
         self.spinner_path = spinner_path
 
         self.resize(600, 400)
-        self.loaded_images = []  
+        self.loaded_images = []
         self.final_mosaic = None
         self.weight_mosaic = None
         self.wcs_metadata = None  # To store mosaic WCS header
         self.astap_exe = self.settings.value("paths/astap", "", type=str)
+
         # Variables to store stretching parameters:
         self.stretch_original_mins = []
         self.stretch_original_medians = []
         self.was_single_channel = False
+
         self.initUI()
 
     def initUI(self):
@@ -4632,12 +4639,11 @@ class MosaicMasterDialog(QDialog):
         layout.addWidget(instructions)
 
         btn_layout = QHBoxLayout()
-        # Button to add image from disk
+
         add_btn = QPushButton("Add Image")
         add_btn.clicked.connect(self.add_image)
         btn_layout.addWidget(add_btn)
 
-        # New button to add an image from one of the ImageManager slots
         add_from_view_btn = QPushButton("Add from View")
         add_from_view_btn.setToolTip("Add an image from any open View")
         add_from_view_btn.clicked.connect(self.add_image_from_view)
@@ -4659,25 +4665,80 @@ class MosaicMasterDialog(QDialog):
         save_btn.clicked.connect(self.save_mosaic_to_new_view)
         btn_layout.addWidget(save_btn)
 
-        # Add the wrench button for settings.
         wrench_btn = QPushButton()
-        wrench_btn.setIcon(QIcon(self.wrench_path))
+        if self.wrench_path:
+            wrench_btn.setIcon(QIcon(self.wrench_path))
         wrench_btn.setToolTip("Mosaic Settings")
         wrench_btn.clicked.connect(self.openSettings)
         btn_layout.addWidget(wrench_btn)
 
         layout.addLayout(btn_layout)
 
-        # Horizontal sizer for checkboxes.
+        # ------------------------------------------------------------------
+        # Mode checkboxes (mutually exclusive: Seestar vs WCS-only)
+        # ------------------------------------------------------------------
         checkbox_layout = QHBoxLayout()
+
         self.forceBlindCheckBox = QCheckBox("Force Blind Solve (ignore existing WCS)")
         checkbox_layout.addWidget(self.forceBlindCheckBox)
-        # New Seestar Mode checkbox:
+
         self.seestarCheckBox = QCheckBox("Seestar Mode")
-        self.seestarCheckBox.setToolTip("Wwen enabled, images are aligned iteratively using astroalign without plate solving.")
+        self.seestarCheckBox.setToolTip(
+            "When enabled, images are aligned iteratively using astroalign without plate solving."
+        )
         checkbox_layout.addWidget(self.seestarCheckBox)
+
+        self.wcsOnlyCheckBox = QCheckBox("Disable Star Alignment (WCS placement only)")
+        self.wcsOnlyCheckBox.setToolTip(
+            "Skips astroalign + refined alignment.\n"
+            "Panels are only reprojected into the mosaic celestial-sphere frame using WCS, then blended.\n"
+            "Useful when panels have little/no overlap but have valid WCS."
+        )
+        checkbox_layout.addWidget(self.wcsOnlyCheckBox)
+
         layout.addLayout(checkbox_layout)
 
+        # Persisted WCS-only
+        _settings = QSettings("SetiAstro", "SASpro")
+        self.wcsOnlyCheckBox.setChecked(_settings.value("mosaic/wcs_only", False, type=bool))
+
+        # Helpers ----------------------------------------------------------
+        def _set_checked(cb: QCheckBox, checked: bool):
+            cb.blockSignals(True)
+            cb.setChecked(checked)
+            cb.blockSignals(False)
+
+        def _sync_wcs_only_ui():
+            # WCS-only disables transform selection (because refinement is unused)
+            wcs_only = self.wcsOnlyCheckBox.isChecked()
+            if hasattr(self, "transform_combo"):
+                self.transform_combo.setEnabled(not wcs_only)
+                self.transform_combo.setToolTip(
+                    "" if not wcs_only else "Disabled because WCS-only placement is enabled."
+                )
+
+        def _on_seestar_changed(state: int):
+            # If Seestar is turned ON, force WCS-only OFF
+            if self.seestarCheckBox.isChecked():
+                _set_checked(self.wcsOnlyCheckBox, False)
+            _sync_wcs_only_ui()
+
+        def _on_wcs_only_changed(state: int):
+            # Persist setting first (this one is user-visible preference)
+            QSettings("SetiAstro", "SASpro").setValue("mosaic/wcs_only", self.wcsOnlyCheckBox.isChecked())
+
+            # If WCS-only is turned ON, force Seestar OFF
+            if self.wcsOnlyCheckBox.isChecked():
+                _set_checked(self.seestarCheckBox, False)
+            _sync_wcs_only_ui()
+
+        # Wire handlers (NO lambdas, no duplicates)
+        self.seestarCheckBox.stateChanged.connect(_on_seestar_changed)
+        self.wcsOnlyCheckBox.stateChanged.connect(_on_wcs_only_changed)
+
+        # ------------------------------------------------------------------
+        # Other controls
+        # ------------------------------------------------------------------
         self.normalizeCheckBox = QCheckBox("Normalize images (median match)")
         self.normalizeCheckBox.setChecked(True)
         layout.addWidget(self.normalizeCheckBox)
@@ -4701,24 +4762,21 @@ class MosaicMasterDialog(QDialog):
             "Precise — Full WCS: astropy.reproject per channel; slowest, most exact."
         )
 
-        # Persist user choice
-        _settings = QSettings("SetiAstro", "SASpro")
-        _default_mode = _settings.value("mosaic/reproject_mode",
-                                        "Fast — SIP-aware (Exact Remap)")
-        if _default_mode not in [self.reprojectModeCombo.itemText(i) for i in range(self.reprojectModeCombo.count())]:
+        _default_mode = _settings.value("mosaic/reproject_mode", "Fast — SIP-aware (Exact Remap)")
+        valid_modes = [self.reprojectModeCombo.itemText(i) for i in range(self.reprojectModeCombo.count())]
+        if _default_mode not in valid_modes:
             _default_mode = "Fast — SIP-aware (Exact Remap)"
         self.reprojectModeCombo.setCurrentText(_default_mode)
         self.reprojectModeCombo.currentTextChanged.connect(
             lambda t: QSettings("SetiAstro", "SASpro").setValue("mosaic/reproject_mode", t)
         )
 
-        # Add to layout where the old checkbox lived
         row = QHBoxLayout()
         row.addWidget(self.reprojectModeLabel)
         row.addWidget(self.reprojectModeCombo, 1)
         layout.addLayout(row)
 
-
+        # Transform selection
         self.transform_combo = QComboBox()
         self.transform_combo.addItems([
             "Partial Affine Transform",
@@ -4726,11 +4784,14 @@ class MosaicMasterDialog(QDialog):
             "Homography Transform",
             "Polynomial Warp Based Transform"
         ])
-        # Set the default selection to "Affine Transform" (index 1)
         self.transform_combo.setCurrentIndex(1)
         layout.addWidget(QLabel("Select Transformation Method:"))
         layout.addWidget(self.transform_combo)
 
+        # Now that transform_combo exists, apply WCS-only UI state
+        _sync_wcs_only_ui()
+
+        # List + status + spinner
         self.images_list = QListWidget()
         self.images_list.setSelectionMode(self.images_list.SelectionMode.SingleSelection)
         layout.addWidget(self.images_list)
@@ -4740,13 +4801,15 @@ class MosaicMasterDialog(QDialog):
 
         self.spinnerLabel = QLabel(self)
         self.spinnerLabel.setAlignment(Qt.AlignmentFlag.AlignCenter)
-        self.spinnerMovie = QMovie(self.spinner_path)
-        self.spinnerLabel.setMovie(self.spinnerMovie)
-        self.spinnerLabel.hide() 
+        self.spinnerMovie = QMovie(self.spinner_path) if self.spinner_path else None
+        if self.spinnerMovie:
+            self.spinnerLabel.setMovie(self.spinnerMovie)
+        self.spinnerLabel.hide()
         layout.addWidget(self.spinnerLabel)
 
         self.setLayout(layout)
 
+
     def _target_median_from_first(self, items):
         # Pick a stable target (median of first image after safe clipping)
         if not items:
@@ -5293,6 +5356,7 @@ class MosaicMasterDialog(QDialog):
 
     # ---------- Align (Entry Point) ----------
     def align_images(self):
+        # Seestar mode is its own pipeline
         if self.seestarCheckBox.isChecked():
             self.align_images_seestar_mode()
             return
@@ -5303,28 +5367,33 @@ class MosaicMasterDialog(QDialog):
 
         # Show spinner and start animation.
         self.spinnerLabel.show()
-        self.spinnerMovie.start()
+        if self.spinnerMovie:
+            self.spinnerMovie.start()
         QApplication.processEvents()
 
-        # Step 1: Force blind solve if requested.
+        # ------------------------------------------------------------
+        # 1) Plate solve (unless already solved and not forcing blind)
+        # ------------------------------------------------------------
         force_blind = self.forceBlindCheckBox.isChecked()
-        images_to_process = (self.loaded_images if force_blind
-                            else [item for item in self.loaded_images if item.get("wcs") is None])
+        images_to_process = (
+            self.loaded_images if force_blind
+            else [item for item in self.loaded_images if item.get("wcs") is None]
+        )
 
-        # Process each image for plate solving.
         for item in images_to_process:
-            # Check if ASTAP is set.
+            # Ensure ASTAP path (or fall back to blind solve)
             if not self.astap_exe or not os.path.exists(self.astap_exe):
                 executable_filter = "Executables (*.exe);;All Files (*)" if sys.platform.startswith("win") else "Executables (*)"
                 new_path, _ = QFileDialog.getOpenFileName(self, "Select ASTAP Executable", "", executable_filter)
                 if new_path:
                     self._save_astap_exe_to_settings(new_path)
+                    self.astap_exe = new_path
                     QMessageBox.information(self, "Mosaic Master", "ASTAP path updated successfully.")
                 else:
-                    QMessageBox.warning(self, "Mosaic Master", "ASTAP path not provided. Falling back to blind solve.")
+                    self.status_label.setText(f"No ASTAP path; blind solving {item['path']}...")
+                    QApplication.processEvents()
                     solved_header = self.perform_blind_solve(item)
                     if solved_header:
-                        # normalize + build WCS with relax=True so SIP is retained
                         solved_header = self._normalize_wcs_header(solved_header, item["image"].shape)
                         item["wcs"] = self._build_wcs(solved_header, item["image"].shape)
                     continue
@@ -5338,49 +5407,53 @@ class MosaicMasterDialog(QDialog):
                 self.status_label.setText(f"ASTAP failed for {item['path']}. Falling back to blind solve...")
                 QApplication.processEvents()
                 solved_header = self.perform_blind_solve(item)
-            else:
-                self.status_label.setText(f"Plate solve successful using ASTAP for {item['path']}.")
 
             if solved_header:
-                # Single, centralized sanitize → keeps SIP and fixes types
                 solved_header = self._normalize_wcs_header(solved_header, item["image"].shape)
                 item["wcs"] = self._build_wcs(solved_header, item["image"].shape)
             else:
-                print(f"Plate solving failed for {item['path']}.")
+                print(f"[Mosaic] Plate solving failed for {item['path']}.")
 
-        # After processing, get all images with valid WCS.
+        # ------------------------------------------------------------
+        # 2) Gather WCS-valid panels
+        # ------------------------------------------------------------
         wcs_items = [x for x in self.loaded_images if x.get("wcs") is not None]
         if not wcs_items:
-            print("No images have WCS, skipping WCS alignment.")
-            self.spinnerMovie.stop()
+            print("[Mosaic] No images have WCS; cannot build WCS mosaic.")
+            if self.spinnerMovie:
+                self.spinnerMovie.stop()
             self.spinnerLabel.hide()
             return
 
+        # ------------------------------------------------------------
+        # 3) Establish mosaic WCS + output bounding box
+        # ------------------------------------------------------------
+        reference_wcs = self._build_wcs(
+            wcs_items[0]["wcs"].to_header(relax=True),
+            wcs_items[0]["image"].shape
+        ).deepcopy()
 
-        # Use the first image's WCS as reference and compute the mosaic bounding box.
-        # (Rebuild with relax=True just in case, then deepcopy.)
-        reference_wcs = self._build_wcs(wcs_items[0]["wcs"].to_header(relax=True), wcs_items[0]["image"].shape).deepcopy()
         min_x, min_y, max_x, max_y = self.compute_mosaic_bounding_box(wcs_items, reference_wcs)
-        mosaic_width = int(max_x - min_x)
+        mosaic_width  = int(max_x - min_x)
         mosaic_height = int(max_y - min_y)
 
         if mosaic_width < 1 or mosaic_height < 1:
-            print("ERROR: Computed mosaic size is invalid. Check WCS or inputs.")
-            self.spinnerMovie.stop()
+            print("[Mosaic] ERROR: Computed mosaic size invalid. Check WCS/inputs.")
+            if self.spinnerMovie:
+                self.spinnerMovie.stop()
             self.spinnerLabel.hide()
             return
 
-        # Adjust the reference WCS so that (min_x, min_y) becomes (0,0).
         mosaic_wcs = reference_wcs.deepcopy()
         mosaic_wcs.wcs.crpix[0] -= min_x
         mosaic_wcs.wcs.crpix[1] -= min_y
-        # keep SIP in the stored header
         self.wcs_metadata = mosaic_wcs.to_header(relax=True)
 
-        # Set up accumulators.
+        # ------------------------------------------------------------
+        # 4) Allocate accumulators
+        # ------------------------------------------------------------
         is_color = any(not item["is_mono"] for item in wcs_items)
 
-        # stats for optional "unstretch"
         self.stretch_original_mins = []
         self.stretch_original_medians = []
         self.was_single_channel = (not is_color)
@@ -5391,49 +5464,77 @@ class MosaicMasterDialog(QDialog):
             self.final_mosaic, _ = smart_zeros((mosaic_height, mosaic_width), dtype=np.float32)
         self.weight_mosaic, _ = smart_zeros((mosaic_height, mosaic_width), dtype=np.float32)
 
+        # ------------------------------------------------------------
+        # 5) Normalization target — compute ONCE, apply ALWAYS
+        # ------------------------------------------------------------
+        did_normalize = bool(self.normalizeCheckBox.isChecked())
+
+        # IMPORTANT: compute from RAW first panel image, not from dict list in a way that can drift
+        if did_normalize:
+            try:
+                first_raw = wcs_items[0]["image"]
+                self._mosaic_target_median = self._target_median_from_first([{"image": first_raw}])
+            except Exception:
+                # fallback: compute directly
+                a0 = wcs_items[0]["image"].astype(np.float32, copy=False)
+                if a0.ndim == 3:
+                    a0m = np.mean(a0, axis=2)
+                else:
+                    a0m = a0
+                lo = np.percentile(a0m, 1)
+                hi = np.percentile(a0m, 99)
+                self._mosaic_target_median = float(max(np.median(np.clip(a0m, lo, hi)), 1e-6))
+
+            print(f"[Mosaic] normalization target median = {self._mosaic_target_median:.6g}")
+
+        # Reprojection helpers cache
+        if not hasattr(self, "_H_cache"):
+            self._H_cache = {}
+
+        # WCS-only toggle (no star alignment/refinement)
+        wcs_only = bool(getattr(self, "wcsOnlyCheckBox", None) and self.wcsOnlyCheckBox.isChecked())
+
+        # ------------------------------------------------------------
+        # 6) Main loop: normalize -> reproject -> (optional) star align -> accumulate
+        # ------------------------------------------------------------
         first_image = True
+
         for idx, itm in enumerate(wcs_items):
             arr = itm["image"]
+
             self.status_label.setText(f"Mapping {itm['path']} into mosaic frame...")
             QApplication.processEvents()
 
             img_lin = arr.astype(np.float32, copy=False)
 
-            # --- record original stats for optional "unstretch" ---
+            # --- record original stats (pre-normalization) for optional unstretch ---
             mono_for_stats = img_lin if img_lin.ndim == 2 else np.mean(img_lin, axis=2)
             self.stretch_original_mins.append(float(np.min(mono_for_stats)))
             self.stretch_original_medians.append(float(np.median(mono_for_stats)))
 
-            # 1) optional median normalization only
-            if self.normalizeCheckBox.isChecked():
-                target_med = getattr(self, "_mosaic_target_median", None)
-                if target_med is None:
-                    self._mosaic_target_median = self._target_median_from_first(wcs_items)
-                    target_med = self._mosaic_target_median
-                img_lin = self._normalize_linear(img_lin, target_med)
-
-            # 2) Reprojection (3 modes)
-            if not hasattr(self, "_H_cache"):
-                self._H_cache = {}
+            # --- ALWAYS normalize here if enabled (even in WCS-only) ---
+            if did_normalize:
+                img_lin = self._normalize_linear(img_lin, float(self._mosaic_target_median))
 
+            # --- Reprojection mode ---
             mode = self.reprojectModeCombo.currentText()
 
             if mode.startswith("Fast — SIP"):
-                # Exact dense remap (SIP-aware), tiled; keep mono as 2D
                 reprojected = self._warp_via_wcs_remap_exact(
                     img_lin, itm["wcs"], mosaic_wcs, (mosaic_height, mosaic_width), tile=512
-                ).astype(np.float32)
+                ).astype(np.float32, copy=False)
+                reprojected = self._polish_reprojected(reprojected)
                 reproj_red = reprojected[..., 0] if reprojected.ndim == 3 else reprojected
 
             elif mode.startswith("Fast — Homography"):
-                # Single global homography; keep mono as 2D
                 reprojected = self._warp_via_wcs_homography(
                     img_lin, itm["wcs"], mosaic_wcs, (mosaic_height, mosaic_width), H_cache=self._H_cache
-                ).astype(np.float32)
+                ).astype(np.float32, copy=False)
+                reprojected = self._polish_reprojected(reprojected)
                 reproj_red = reprojected[..., 0] if reprojected.ndim == 3 else reprojected
 
             else:
-                # Precise — Full WCS (astropy.reproject); mono stays 2D
+                # Precise — Full WCS
                 if img_lin.ndim == 3:
                     channels = []
                     for c in range(3):
@@ -5443,52 +5544,49 @@ class MosaicMasterDialog(QDialog):
                     reprojected = np.stack(channels, axis=-1)
                     reproj_red = reprojected[..., 0]
                 else:
-                    reproj_red, _ = reproject_interp((img_lin, itm["wcs"]), mosaic_wcs,
-                                                    shape_out=(mosaic_height, mosaic_width))
-                    reprojected = np.nan_to_num(reproj_red, nan=0.0).astype(np.float32)  # 2D mono
-                    # no fake stacking here
-
-            self.status_label.setText(f"WCS map: {itm['path']} processed.")
-            QApplication.processEvents()
-
-            # --- Stellar Alignment ---
-            if not first_image:
-                transform_method = self.transform_combo.currentText()
-                mosaic_gray = (self.final_mosaic if self.final_mosaic.ndim == 2
-                            else np.mean(self.final_mosaic, axis=-1))
-                try:
-                    self.status_label.setText("Computing affine transform with astroalign...")
-                    QApplication.processEvents()
-                    transform_obj, (src_pts, dst_pts) = self._aa_find_transform_with_backoff(reproj_red, mosaic_gray)
-                    transform_matrix = transform_obj.params[0:2, :].astype(np.float32)
-                    self.status_label.setText("Astroalign computed transform successfully.")
-                except Exception as e:
-                    self.status_label.setText(f"Astroalign failed: {e}. Using identity transform.")
-                    transform_matrix = np.eye(2, 3, dtype=np.float32)
+                    rpj, _ = reproject_interp((img_lin, itm["wcs"]), mosaic_wcs,
+                                            shape_out=(mosaic_height, mosaic_width))
+                    reprojected = np.nan_to_num(rpj, nan=0.0).astype(np.float32)
+                    reproj_red = reprojected  # 2D
 
-                A = transform_matrix[:, :2]
-                scale1 = np.linalg.norm(A[:, 0])
-                scale2 = np.linalg.norm(A[:, 1])
-                print(f"Computed affine scales: {scale1:.6f}, {scale2:.6f}")
+            # --- Optional star alignment/refinement (skipped in WCS-only) ---
+            if wcs_only:
+                aligned = reprojected
+                if first_image:
+                    first_image = False
+            else:
+                if first_image:
+                    aligned = reprojected
+                    first_image = False
+                else:
+                    transform_method = self.transform_combo.currentText()
+                    mosaic_gray = (self.final_mosaic if self.final_mosaic.ndim == 2
+                                else np.mean(self.final_mosaic, axis=-1))
 
-                self.status_label.setText("Affine alignment computed. Warping image...")
-                QApplication.processEvents()
-                affine_aligned = cv2.warpAffine(reprojected, transform_matrix, (mosaic_width, mosaic_height),
-                                                flags=cv2.INTER_LANCZOS4)
-                aligned = affine_aligned
+                    try:
+                        self.status_label.setText("Computing affine transform with astroalign...")
+                        QApplication.processEvents()
+                        transform_obj, (src_pts, dst_pts) = self._aa_find_transform_with_backoff(reproj_red, mosaic_gray)
+                        transform_matrix = transform_obj.params[0:2, :].astype(np.float32)
+                    except Exception as e:
+                        self.status_label.setText(f"Astroalign failed: {e}. Using identity transform.")
+                        transform_matrix = np.eye(2, 3, dtype=np.float32)
 
-                if transform_method in ["Homography Transform", "Polynomial Warp Based Transform"]:
-                    self.status_label.setText(f"Starting refined alignment using {transform_method}...")
-                    QApplication.processEvents()
-                    refined_result = self.refined_alignment(affine_aligned, mosaic_gray, method=transform_method)
-                    if refined_result is not None:
-                        aligned, best_inliers2 = refined_result
-                        self.status_label.setText(f"Refined alignment succeeded with {best_inliers2} inliers.")
-                    else:
-                        self.status_label.setText("Refined alignment failed; falling back to affine alignment.")
-            else:
-                aligned = reprojected
-                first_image = False
+                    affine_aligned = cv2.warpAffine(
+                        reprojected, transform_matrix, (mosaic_width, mosaic_height),
+                        flags=cv2.INTER_LANCZOS4
+                    )
+                    aligned = affine_aligned
+
+                    if transform_method in ["Homography Transform", "Polynomial Warp Based Transform"]:
+                        self.status_label.setText(f"Starting refined alignment using {transform_method}...")
+                        QApplication.processEvents()
+                        refined_result = self.refined_alignment(affine_aligned, mosaic_gray, method=transform_method)
+                        if refined_result is not None:
+                            aligned, best_inliers2 = refined_result
+                            self.status_label.setText(f"Refined alignment succeeded with {best_inliers2} inliers.")
+                        else:
+                            self.status_label.setText("Refined alignment failed; falling back to affine alignment.")
 
             # If mosaic is color but aligned is mono, expand for accumulation only
             if is_color and aligned.ndim == 2:
@@ -5496,13 +5594,13 @@ class MosaicMasterDialog(QDialog):
 
             gray_aligned = aligned[..., 0] if aligned.ndim == 3 else aligned
 
-            # Compute weight mask
+            # --- Weight mask ---
             binary_mask = (gray_aligned > 0).astype(np.uint8)
             smooth_mask = cv2.distanceTransform(binary_mask, cv2.DIST_L2, 5)
             smooth_mask = (smooth_mask / np.max(smooth_mask)) if np.max(smooth_mask) > 0 else binary_mask.astype(np.float32)
             smooth_mask = cv2.GaussianBlur(smooth_mask, (15, 15), 0)
 
-            # Accumulate
+            # --- Accumulate ---
             if is_color:
                 self.final_mosaic += aligned * smooth_mask[..., np.newaxis]
             else:
@@ -5512,22 +5610,27 @@ class MosaicMasterDialog(QDialog):
             self.status_label.setText(f"Processed: {itm['path']}")
             QApplication.processEvents()
 
-        # Final blending.
-        nonzero_mask = (self.weight_mosaic > 0)
+        # ------------------------------------------------------------
+        # 7) Final blend
+        # ------------------------------------------------------------
         if is_color:
-            self.final_mosaic = np.where(self.weight_mosaic[..., None] > 0,
-                                        self.final_mosaic / self.weight_mosaic[..., None],
-                                        self.final_mosaic)
+            self.final_mosaic = np.where(
+                self.weight_mosaic[..., None] > 0,
+                self.final_mosaic / np.maximum(self.weight_mosaic[..., None], 1e-12),
+                self.final_mosaic
+            )
         else:
-            self.final_mosaic[nonzero_mask] = self.final_mosaic[nonzero_mask] / self.weight_mosaic[nonzero_mask]
+            nz = (self.weight_mosaic > 0)
+            self.final_mosaic[nz] = self.final_mosaic[nz] / np.maximum(self.weight_mosaic[nz], 1e-12)
 
-        print("WCS + Star Alignment Complete.")
-        self.status_label.setText("WCS + Star Alignment Complete. De-Normalizing Mosaic...")
+        self.status_label.setText("Mosaic built. De-normalizing mosaic...")
+        QApplication.processEvents()
 
-        # Call-guard: only unstretch if we normalized AND we actually recorded stats
-        did_normalize = self.normalizeCheckBox.isChecked()
+        # ------------------------------------------------------------
+        # 8) Optional “unstretch” (your existing logic)
+        # ------------------------------------------------------------
         if (did_normalize and
-            hasattr(self, "_mosaic_target_median") and self._mosaic_target_median > 0 and
+            hasattr(self, "_mosaic_target_median") and float(self._mosaic_target_median) > 0 and
             getattr(self, "stretch_original_medians", None) and len(self.stretch_original_medians) > 0 and
             getattr(self, "stretch_original_mins", None) and len(self.stretch_original_mins) > 0):
             self.final_mosaic = self.unstretch_image(self.final_mosaic)
@@ -5535,16 +5638,37 @@ class MosaicMasterDialog(QDialog):
         self.status_label.setText("Final Mosaic Ready.")
         QApplication.processEvents()
 
-        display_image = (np.stack([self.final_mosaic]*3, axis=-1)
+        display_image = (np.stack([self.final_mosaic] * 3, axis=-1)
                         if self.final_mosaic.ndim == 2 else self.final_mosaic)
         display_image = self._autostretch_if_requested(display_image)
-        MosaicPreviewWindow(display_image, title="Final Mosaic", parent=self,
-                            push_cb=self._push_mosaic_to_new_doc).show()
 
-        self.spinnerMovie.stop()
+        MosaicPreviewWindow(
+            display_image,
+            title=("Final Mosaic (WCS-only)" if wcs_only else "Final Mosaic"),
+            parent=self,
+            push_cb=self._push_mosaic_to_new_doc
+        ).show()
+
+        if self.spinnerMovie:
+            self.spinnerMovie.stop()
         self.spinnerLabel.hide()
         QApplication.processEvents()
+
       
+    def _polish_reprojected(self, img: np.ndarray) -> np.ndarray:
+        # 1) kill NaNs/Infs from reprojection
+        out = np.nan_to_num(img, nan=0.0, posinf=0.0, neginf=0.0).astype(np.float32, copy=False)
+
+        # 2) trim tiny negative ringing
+        out[out < 0] = 0.0
+
+        # 3) optional: zero a 1px border to avoid remap edge seams
+        if out.ndim == 2:
+            out[:1, :] = 0; out[-1:, :] = 0; out[:, :1] = 0; out[:, -1:] = 0
+        else:
+            out[:1, :, :] = 0; out[-1:, :, :] = 0; out[:, :1, :] = 0; out[:, -1:, :] = 0
+
+        return out
 
     def debayer_image(self, image, file_path, header):
         from setiastro.saspro.legacy.numba_utils import debayer_raw_fast, debayer_fits_fast  

setiastro/saspro/star_spikes.py

@@ -53,6 +53,10 @@ class StarSpikesDialogPro(QDialog):
         self.setWindowFlag(Qt.WindowType.Window, True)
         self.setWindowModality(Qt.WindowModality.NonModal)
         self.setModal(False)
+        try:
+            self.setAttribute(Qt.WidgetAttribute.WA_DeleteOnClose, True)
+        except Exception:
+            pass  # older PyQt6 versions        
         self.docman = doc_manager
         self.doc = initial_doc or (self.docman.get_active_document() if self.docman else None)
         self.jwstpupil_path = jwstpupil_path