napari-tmidas 0.2.1__py3-none-any.whl → 0.2.4__py3-none-any.whl

napari_tmidas/_roi_colocalization.py

@@ -11,7 +11,6 @@ and can optionally calculate sizes of these overlapping regions.
 import concurrent.futures
 
 # contextlib is used to suppress exceptions
-import contextlib
 import csv
 import os
 from collections import defaultdict
@@ -41,6 +40,37 @@ from qtpy.QtWidgets import (
 from skimage import measure
 
 
+def convert_semantic_to_instance_labels(image, connectivity=None):
+    """
+    Convert semantic labels (where all objects have the same value) to instance labels.
+
+    Args:
+        image: Label image that may contain semantic labels
+        connectivity: Connectivity for connected component analysis (1, 2, or None for full)
+
+    Returns:
+        Image with instance labels (each connected component gets unique label)
+    """
+    if image is None or np.all(image == 0):
+        return image
+
+    # Get unique non-zero values
+    unique_labels = np.unique(image[image != 0])
+
+    # Quick check: if there's only one unique non-zero value, it's definitely semantic
+    # Otherwise, apply connected components to the entire mask at once (much faster)
+    if len(unique_labels) == 1:
+        # Single semantic label - just label connected components of the binary mask
+        mask = image > 0
+        return measure.label(mask, connectivity=connectivity)
+    else:
+        # Multiple labels - could be instance or semantic
+        # Apply connected components to entire non-zero region at once
+        # This is MUCH faster than iterating over each label value
+        mask = image > 0
+        return measure.label(mask, connectivity=connectivity)
+
+
 def longest_common_substring(s1, s2):
     """Finds the longest common substring between two strings."""
     matcher = SequenceMatcher(None, s1, s2)
@@ -100,9 +130,9 @@ def group_files_by_common_substring(file_lists, channels):
 
         # If matches were found for all channels, add them to the group
         if len(matched_files) == len(channels):
-            # Strip suffixes from the common substring
-            stripped_common_substring = common_substring.rsplit("_", 1)[0]
-            groups[stripped_common_substring] = {
+            # Use the full common substring as the key (don't strip it yet)
+            # This prevents different file pairs from overwriting each other
+            groups[common_substring] = {
                 channel: file_lists[channel][
                     base_files[channel].index(matched_files[channel])
                 ]
@@ -134,6 +164,15 @@ class ColocalizationWorker(QThread):
         get_sizes=False,
         size_method="median",
         output_folder=None,
+        channel2_is_labels=True,
+        channel3_is_labels=True,
+        count_positive=False,
+        threshold_method="percentile",
+        threshold_value=75.0,
+        save_images=True,
+        convert_to_instances_c2=False,
+        convert_to_instances_c3=False,
+        count_c2_positive_for_c3=False,
     ):
         super().__init__()
         self.file_pairs = file_pairs
@@ -141,6 +180,15 @@ class ColocalizationWorker(QThread):
         self.get_sizes = get_sizes
         self.size_method = size_method
         self.output_folder = output_folder
+        self.channel2_is_labels = channel2_is_labels
+        self.channel3_is_labels = channel3_is_labels
+        self.count_positive = count_positive
+        self.threshold_method = threshold_method
+        self.threshold_value = threshold_value
+        self.save_images = save_images
+        self.convert_to_instances_c2 = convert_to_instances_c2
+        self.convert_to_instances_c3 = convert_to_instances_c3
+        self.count_c2_positive_for_c3 = count_c2_positive_for_c3
         self.stop_requested = False
         self.thread_count = max(1, (os.cpu_count() or 4) - 1)  # Default value
 
@@ -167,30 +215,137 @@ class ColocalizationWorker(QThread):
                 header = [
                     "Filename",
                     f"{self.channel_names[0]}_label_id",
-                    f"{self.channel_names[1]}_in_{self.channel_names[0]}_count",
                 ]
 
-                if self.get_sizes:
+                # Add c2_label_id column if in individual mode
+                if (
+                    self.size_method == "individual"
+                    and self.channel2_is_labels
+                ):
+                    header.append(f"{self.channel_names[1]}_label_id")
+
+                # Add Channel 2 columns based on whether it's labels or intensity
+                if self.channel2_is_labels:
+                    if self.size_method != "individual":
+                        # Aggregate mode: add count of C2 labels in C1
+                        header.append(
+                            f"{self.channel_names[1]}_in_{self.channel_names[0]}_count"
+                        )
+                    # Individual mode: no aggregate C2 statistics needed
+                else:
+                    # Channel 2 is intensity: add intensity statistics
                     header.extend(
                         [
-                            f"{self.channel_names[0]}_size",
-                            f"{self.channel_names[1]}_in_{self.channel_names[0]}_size",
+                            f"{self.channel_names[1]}_in_{self.channel_names[0]}_mean",
+                            f"{self.channel_names[1]}_in_{self.channel_names[0]}_median",
+                            f"{self.channel_names[1]}_in_{self.channel_names[0]}_std",
+                            f"{self.channel_names[1]}_in_{self.channel_names[0]}_max",
                         ]
                     )
 
+                if self.get_sizes:
+                    header.append(f"{self.channel_names[0]}_size")
+                    if self.channel2_is_labels:
+                        if self.size_method == "individual":
+                            # Individual mode: one row per c2 label with its size
+                            header.append(f"{self.channel_names[1]}_size")
+                        else:
+                            header.append(
+                                f"{self.channel_names[1]}_in_{self.channel_names[0]}_size"
+                            )
+
                 if len(self.channel_names) == 3:
-                    header.extend(
-                        [
-                            f"{self.channel_names[2]}_in_{self.channel_names[1]}_in_{self.channel_names[0]}_count",
-                            f"{self.channel_names[2]}_not_in_{self.channel_names[1]}_but_in_{self.channel_names[0]}_count",
-                        ]
-                    )
+                    if self.channel2_is_labels and self.channel3_is_labels:
+                        # Both ch2 and ch3 are labels - original 3-channel label mode
+                        header.extend(
+                            [
+                                f"{self.channel_names[2]}_in_{self.channel_names[1]}_in_{self.channel_names[0]}_count",
+                                f"{self.channel_names[2]}_not_in_{self.channel_names[1]}_but_in_{self.channel_names[0]}_count",
+                            ]
+                        )
 
-                    if self.get_sizes:
+                        if self.get_sizes:
+                            if self.size_method == "individual":
+                                # Individual mode: c3 size within each c2 label
+                                header.append(
+                                    f"{self.channel_names[2]}_in_{self.channel_names[1]}_size"
+                                )
+                            else:
+                                header.extend(
+                                    [
+                                        f"{self.channel_names[2]}_in_{self.channel_names[1]}_in_{self.channel_names[0]}_size",
+                                        f"{self.channel_names[2]}_not_in_{self.channel_names[1]}_but_in_{self.channel_names[0]}_size",
+                                    ]
+                                )
+
+                        # Add positive counting columns if requested
+                        if self.count_c2_positive_for_c3:
+                            header.extend(
+                                [
+                                    f"{self.channel_names[1]}_in_{self.channel_names[0]}_positive_for_{self.channel_names[2]}_count",
+                                    f"{self.channel_names[1]}_in_{self.channel_names[0]}_negative_for_{self.channel_names[2]}_count",
+                                    f"{self.channel_names[1]}_in_{self.channel_names[0]}_percent_positive_for_{self.channel_names[2]}",
+                                ]
+                            )
+                    elif (
+                        self.channel2_is_labels and not self.channel3_is_labels
+                    ):
+                        # Ch2 is labels, Ch3 is intensity
+                        if self.size_method == "individual":
+                            # Individual mode: intensity statistics per c2 label
+                            header.extend(
+                                [
+                                    f"{self.channel_names[2]}_in_{self.channel_names[1]}_mean",
+                                    f"{self.channel_names[2]}_in_{self.channel_names[1]}_median",
+                                    f"{self.channel_names[2]}_in_{self.channel_names[1]}_std",
+                                ]
+                            )
+                            # Note: positive counting is not available in individual mode
+                            # Individual mode provides per-label statistics, not aggregate counts
+                        else:
+                            # Aggregate statistics (original behavior)
+                            header.extend(
+                                [
+                                    f"{self.channel_names[2]}_in_{self.channel_names[1]}_in_{self.channel_names[0]}_mean",
+                                    f"{self.channel_names[2]}_in_{self.channel_names[1]}_in_{self.channel_names[0]}_median",
+                                    f"{self.channel_names[2]}_in_{self.channel_names[1]}_in_{self.channel_names[0]}_std",
+                                    f"{self.channel_names[2]}_in_{self.channel_names[1]}_in_{self.channel_names[0]}_max",
+                                    f"{self.channel_names[2]}_not_in_{self.channel_names[1]}_but_in_{self.channel_names[0]}_mean",
+                                    f"{self.channel_names[2]}_not_in_{self.channel_names[1]}_in_{self.channel_names[0]}_median",
+                                    f"{self.channel_names[2]}_not_in_{self.channel_names[1]}_but_in_{self.channel_names[0]}_std",
+                                    f"{self.channel_names[2]}_not_in_{self.channel_names[1]}_but_in_{self.channel_names[0]}_max",
+                                ]
+                            )
+
+                            # Add positive counting columns if requested (only in aggregate mode)
+                            if self.count_positive:
+                                header.extend(
+                                    [
+                                        f"{self.channel_names[1]}_in_{self.channel_names[0]}_positive_for_{self.channel_names[2]}_count",
+                                        f"{self.channel_names[1]}_in_{self.channel_names[0]}_negative_for_{self.channel_names[2]}_count",
+                                        f"{self.channel_names[1]}_in_{self.channel_names[0]}_percent_positive_for_{self.channel_names[2]}",
+                                        f"{self.channel_names[2]}_threshold_used",
+                                    ]
+                                )
+                    elif (
+                        not self.channel2_is_labels and self.channel3_is_labels
+                    ):
+                        # Ch2 is intensity, Ch3 is labels
+                        header.append(
+                            f"{self.channel_names[2]}_in_{self.channel_names[0]}_count"
+                        )
+                        if self.get_sizes:
+                            header.append(
+                                f"{self.channel_names[2]}_in_{self.channel_names[0]}_size"
+                            )
+                    else:
+                        # Both Ch2 and Ch3 are intensity
                         header.extend(
                             [
-                                f"{self.channel_names[2]}_in_{self.channel_names[1]}_in_{self.channel_names[0]}_size",
-                                f"{self.channel_names[2]}_not_in_{self.channel_names[1]}_but_in_{self.channel_names[0]}_size",
+                                f"{self.channel_names[2]}_in_{self.channel_names[0]}_mean",
+                                f"{self.channel_names[2]}_in_{self.channel_names[0]}_median",
+                                f"{self.channel_names[2]}_in_{self.channel_names[0]}_std",
+                                f"{self.channel_names[2]}_in_{self.channel_names[0]}_max",
                             ]
                         )
 
@@ -284,6 +439,17 @@ class ColocalizationWorker(QThread):
             image_c2 = tifffile.imread(filepath_c2)
             image_c3 = tifffile.imread(filepath_c3) if filepath_c3 else None
 
+            # Debugging: Check if images are identical (possible file selection error)
+            if image_c3 is not None and np.array_equal(image_c2, image_c3):
+                print(
+                    "WARNING: Channel 2 and Channel 3 contain IDENTICAL data!"
+                )
+                print(f"  C2 file: {os.path.basename(filepath_c2)}")
+                print(f"  C3 file: {os.path.basename(filepath_c3)}")
+                print(
+                    "  This likely indicates the same file was selected for both channels."
+                )
+
             # Ensure all images have the same shape
             if image_c1.shape != image_c2.shape:
                 raise ValueError(
@@ -318,6 +484,21 @@ class ColocalizationWorker(QThread):
         self, filename, image_c1, image_c2, image_c3=None
     ):
         """Process colocalization between channels"""
+        # Convert semantic labels to instance labels if requested
+        if (
+            self.convert_to_instances_c2
+            and self.channel2_is_labels
+            and image_c2 is not None
+        ):
+            image_c2 = convert_semantic_to_instance_labels(image_c2)
+
+        if (
+            self.convert_to_instances_c3
+            and self.channel3_is_labels
+            and image_c3 is not None
+        ):
+            image_c3 = convert_semantic_to_instance_labels(image_c3)
+
         # Get unique label IDs in image_c1
         label_ids = self.get_nonzero_labels(image_c1)
 
@@ -326,17 +507,48 @@ class ColocalizationWorker(QThread):
         if self.get_sizes:
             roi_sizes = self.calculate_all_rois_size(image_c1)
 
+        # Handle intensity images for channel 2 and 3
+        image_c2_intensity = None
+        image_c3_intensity = None
+
+        if not self.channel2_is_labels:
+            image_c2_intensity = image_c2
+
+        if image_c3 is not None and not self.channel3_is_labels:
+            image_c3_intensity = image_c3
+
         # Process each label
         csv_rows = []
         results = []
 
         for label_id in label_ids:
-            row = self.process_single_roi(
-                filename, label_id, image_c1, image_c2, image_c3, roi_sizes
+            row_or_rows = self.process_single_roi(
+                filename,
+                label_id,
+                image_c1,
+                image_c2,
+                image_c3,
+                roi_sizes,
+                image_c2_intensity,
+                image_c3_intensity,
             )
-            csv_rows.append(row)
 
-            # Extract results as dictionary
+            # Check if we got multiple rows (individual mode) or single row
+            if len(row_or_rows) == 0:
+                # No rows returned (e.g., no C2 labels in this C1 ROI in individual mode)
+                continue
+            elif isinstance(row_or_rows[0], list):
+                # Multiple rows returned (one per c2 label)
+                for row in row_or_rows:
+                    csv_rows.append(row)
+                    # For individual mode, skip creating result_dict (simplified visualization)
+                continue
+            else:
+                # Single row returned
+                row = row_or_rows
+                csv_rows.append(row)
+
+            # Extract results as dictionary (only for non-individual mode)
             result_dict = {"label_id": label_id, "ch2_in_ch1_count": row[2]}
 
             idx = 3
@@ -346,15 +558,65 @@ class ColocalizationWorker(QThread):
                 idx += 2
 
             if image_c3 is not None:
-                result_dict["ch3_in_ch2_in_ch1_count"] = row[idx]
-                result_dict["ch3_not_in_ch2_but_in_ch1_count"] = row[idx + 1]
-                idx += 2
-
-                if self.get_sizes:
-                    result_dict["ch3_in_ch2_in_ch1_size"] = row[idx]
-                    result_dict["ch3_not_in_ch2_but_in_ch1_size"] = row[
+                # Map CSV row columns to result_dict depending on channel modes
+                if self.channel2_is_labels and self.channel3_is_labels:
+                    # Both ch2 and ch3 are labels: two counts (in c2 & not in c2)
+                    result_dict["ch3_in_ch2_in_ch1_count"] = row[idx]
+                    result_dict["ch3_not_in_ch2_but_in_ch1_count"] = row[
                         idx + 1
                     ]
+                    idx += 2
+
+                    if self.get_sizes:
+                        result_dict["ch3_in_ch2_in_ch1_size"] = row[idx]
+                        result_dict["ch3_not_in_ch2_but_in_ch1_size"] = row[
+                            idx + 1
+                        ]
+                elif self.channel2_is_labels and not self.channel3_is_labels:
+                    # ch2 labels, ch3 intensity: many intensity stats were appended
+                    # Map the first group of intensity stats to ch3_in_ch2_in_ch1_* keys
+                    result_dict["ch3_in_ch2_in_ch1_mean"] = row[idx]
+                    result_dict["ch3_in_ch2_in_ch1_median"] = row[idx + 1]
+                    result_dict["ch3_in_ch2_in_ch1_std"] = row[idx + 2]
+                    result_dict["ch3_in_ch2_in_ch1_max"] = row[idx + 3]
+                    result_dict["ch3_not_in_ch2_but_in_ch1_mean"] = row[
+                        idx + 4
+                    ]
+                    result_dict["ch3_not_in_ch2_but_in_ch1_median"] = row[
+                        idx + 5
+                    ]
+                    result_dict["ch3_not_in_ch2_but_in_ch1_std"] = row[idx + 6]
+                    result_dict["ch3_not_in_ch2_but_in_ch1_max"] = row[idx + 7]
+                    idx += 8
+
+                    # If positive counting (intensity mode) appended extra columns
+                    if self.count_positive:
+                        result_dict["ch2_in_ch1_positive_for_ch3_count"] = row[
+                            idx
+                        ]
+                        result_dict["ch2_in_ch1_negative_for_ch3_count"] = row[
+                            idx + 1
+                        ]
+                        result_dict["ch2_in_ch1_percent_positive_for_ch3"] = (
+                            row[idx + 2]
+                        )
+                        result_dict["ch3_threshold_used"] = row[idx + 3]
+                        idx += 4
+                elif not self.channel2_is_labels and self.channel3_is_labels:
+                    # ch2 intensity, ch3 labels: single count (ch3 in ch1)
+                    result_dict["ch3_in_ch1_count"] = row[idx]
+                    idx += 1
+
+                    if self.get_sizes:
+                        result_dict["ch3_in_ch1_size"] = row[idx]
+                        idx += 1
+                else:
+                    # Both channels are intensity: map intensity stats
+                    result_dict["ch3_in_ch1_mean"] = row[idx]
+                    result_dict["ch3_in_ch1_median"] = row[idx + 1]
+                    result_dict["ch3_in_ch1_std"] = row[idx + 2]
+                    result_dict["ch3_in_ch1_max"] = row[idx + 3]
+                    idx += 4
 
             results.append(result_dict)
 
@@ -368,66 +630,230 @@ class ColocalizationWorker(QThread):
         return output
 
     def process_single_roi(
-        self, filename, label_id, image_c1, image_c2, image_c3, roi_sizes
+        self,
+        filename,
+        label_id,
+        image_c1,
+        image_c2,
+        image_c3,
+        roi_sizes,
+        image_c2_intensity=None,
+        image_c3_intensity=None,
     ):
-        """Process a single ROI for colocalization analysis."""
+        """Process a single ROI for colocalization analysis.
+
+        Returns:
+            list or list of lists: Single row for non-individual mode,
+                                  list of rows (one per c2 label) for individual mode
+        """
         # Create masks once
         mask_roi = image_c1 == label_id
-        mask_c2 = image_c2 != 0
 
-        # Calculate counts
-        c2_in_c1_count = self.count_unique_nonzero(
-            image_c2, mask_roi & mask_c2
-        )
+        # Check if we should create individual rows for each c2 label
+        if (
+            self.size_method == "individual"
+            and self.channel2_is_labels
+            and (
+                image_c3 is not None
+                and not self.channel3_is_labels
+                or self.get_sizes
+            )
+        ):
+            # Individual mode: return one row per c2 label
+            return self._process_individual_c2_labels(
+                filename,
+                label_id,
+                image_c1,
+                image_c2,
+                image_c3,
+                roi_sizes,
+                mask_roi,
+                image_c2_intensity,
+                image_c3_intensity,
+            )
 
         # Build the result row
-        row = [filename, int(label_id), c2_in_c1_count]
+        row = [filename, int(label_id)]
+
+        # Handle Channel 2 based on whether it's labels or intensity
+        if self.channel2_is_labels:
+            mask_c2 = image_c2 != 0
+            # Calculate counts
+            c2_in_c1_count = self.count_unique_nonzero(
+                image_c2, mask_roi & mask_c2
+            )
+            row.append(c2_in_c1_count)
+        else:
+            # Channel 2 is intensity - calculate intensity statistics
+            intensity_img = (
+                image_c2_intensity
+                if image_c2_intensity is not None
+                else image_c2
+            )
+            stats_c2 = self.calculate_intensity_stats(intensity_img, mask_roi)
+            row.extend(
+                [
+                    stats_c2["mean"],
+                    stats_c2["median"],
+                    stats_c2["std"],
+                    stats_c2["max"],
+                ]
+            )
 
         # Add size information if requested
         if self.get_sizes:
             size = roi_sizes.get(int(label_id), 0)
-            c2_in_c1_size = self.calculate_coloc_size(
-                image_c1, image_c2, label_id
-            )
-            row.extend([size, c2_in_c1_size])
+            row.append(size)
+
+            if self.channel2_is_labels:
+                # Calculate aggregate size (non-individual mode handles this)
+                c2_in_c1_size = self.calculate_coloc_size(
+                    image_c1, image_c2, label_id
+                )
+                row.append(c2_in_c1_size)
 
         # Handle third channel if present
         if image_c3 is not None:
-            mask_c3 = image_c3 != 0
+            if self.channel2_is_labels and self.channel3_is_labels:
+                # Both ch2 and ch3 are labels - original 3-channel label mode
+                mask_c2 = image_c2 != 0
+                mask_c3 = image_c3 != 0
+
+                # Calculate third channel statistics
+                c3_in_c2_in_c1_count = self.count_unique_nonzero(
+                    image_c3, mask_roi & mask_c2 & mask_c3
+                )
+                c3_not_in_c2_but_in_c1_count = self.count_unique_nonzero(
+                    image_c3, mask_roi & ~mask_c2 & mask_c3
+                )
 
-            # Calculate third channel statistics
-            c3_in_c2_in_c1_count = self.count_unique_nonzero(
-                image_c3, mask_roi & mask_c2 & mask_c3
-            )
-            c3_not_in_c2_but_in_c1_count = self.count_unique_nonzero(
-                image_c3, mask_roi & ~mask_c2 & mask_c3
-            )
+                row.extend(
+                    [c3_in_c2_in_c1_count, c3_not_in_c2_but_in_c1_count]
+                )
 
-            row.extend([c3_in_c2_in_c1_count, c3_not_in_c2_but_in_c1_count])
+                # Add size information for third channel if requested
+                if self.get_sizes:
+                    # Calculate aggregate sizes (non-individual mode)
+                    c3_in_c2_in_c1_size = self.calculate_coloc_size(
+                        image_c1,
+                        image_c2,
+                        label_id,
+                        mask_c2=True,
+                        image_c3=image_c3,
+                    )
+                    c3_not_in_c2_but_in_c1_size = self.calculate_coloc_size(
+                        image_c1,
+                        image_c2,
+                        label_id,
+                        mask_c2=False,
+                        image_c3=image_c3,
+                    )
+                    row.extend(
+                        [c3_in_c2_in_c1_size, c3_not_in_c2_but_in_c1_size]
+                    )
 
-            # Add size information for third channel if requested
-            if self.get_sizes:
-                c3_in_c2_in_c1_size = self.calculate_coloc_size(
-                    image_c1,
-                    image_c2,
-                    label_id,
-                    mask_c2=True,
-                    image_c3=image_c3,
+                # Count C2 objects positive for C3 if requested
+                if self.count_c2_positive_for_c3:
+                    positive_counts = self.count_c2_positive_for_c3_labels(
+                        image_c2, image_c3, mask_roi
+                    )
+                    row.extend(
+                        [
+                            positive_counts["c2_positive_for_c3_count"],
+                            positive_counts["c2_negative_for_c3_count"],
+                            positive_counts["c2_percent_positive_for_c3"],
+                        ]
+                    )
+            elif self.channel2_is_labels and not self.channel3_is_labels:
+                # Ch2 is labels, Ch3 is intensity
+                mask_c2 = image_c2 != 0
+                intensity_img = (
+                    image_c3_intensity
+                    if image_c3_intensity is not None
+                    else image_c3
+                )
+
+                # Calculate aggregate intensity statistics (non-individual mode)
+                # Calculate intensity where c2 is present in c1
+                mask_c2_in_c1 = mask_roi & mask_c2
+                stats_c2_in_c1 = self.calculate_intensity_stats(
+                    intensity_img, mask_c2_in_c1
+                )
+
+                # Calculate intensity where c2 is NOT present in c1
+                mask_not_c2_in_c1 = mask_roi & ~mask_c2
+                stats_not_c2_in_c1 = self.calculate_intensity_stats(
+                    intensity_img, mask_not_c2_in_c1
+                )
+
+                # Add intensity statistics to row
+                row.extend(
+                    [
+                        stats_c2_in_c1["mean"],
+                        stats_c2_in_c1["median"],
+                        stats_c2_in_c1["std"],
+                        stats_c2_in_c1["max"],
+                        stats_not_c2_in_c1["mean"],
+                        stats_not_c2_in_c1["median"],
+                        stats_not_c2_in_c1["std"],
+                        stats_not_c2_in_c1["max"],
+                    ]
                 )
-                c3_not_in_c2_but_in_c1_size = self.calculate_coloc_size(
-                    image_c1,
-                    image_c2,
-                    label_id,
-                    mask_c2=False,
-                    image_c3=image_c3,
+
+                # Count positive Channel 2 objects if requested
+                if self.count_positive:
+                    positive_counts = self.count_positive_objects(
+                        image_c2,
+                        intensity_img,
+                        mask_roi,
+                        self.threshold_method,
+                        self.threshold_value,
+                    )
+                    row.extend(
+                        [
+                            positive_counts["positive_c2_objects"],
+                            positive_counts["negative_c2_objects"],
+                            positive_counts["percent_positive"],
+                            positive_counts["threshold_used"],
+                        ]
+                    )
+            elif not self.channel2_is_labels and self.channel3_is_labels:
+                # Ch2 is intensity, Ch3 is labels - count Ch3 objects in Ch1
+                mask_c3 = image_c3 != 0
+                c3_in_c1_count = self.count_unique_nonzero(
+                    image_c3, mask_roi & mask_c3
+                )
+                row.append(c3_in_c1_count)
+
+                if self.get_sizes:
+                    c3_in_c1_size = self.calculate_coloc_size(
+                        image_c1, image_c3, label_id
+                    )
+                    row.append(c3_in_c1_size)
+            else:
+                # Both Ch2 and Ch3 are intensity - just add Ch3 stats to Ch1 ROIs
+                intensity_img = (
+                    image_c3_intensity
+                    if image_c3_intensity is not None
+                    else image_c3
+                )
+                stats_c3 = self.calculate_intensity_stats(
+                    intensity_img, mask_roi
+                )
+                row.extend(
+                    [
+                        stats_c3["mean"],
+                        stats_c3["median"],
+                        stats_c3["std"],
+                        stats_c3["max"],
+                    ]
                 )
-                row.extend([c3_in_c2_in_c1_size, c3_not_in_c2_but_in_c1_size])
 
         return row
 
     def save_output_image(self, results, file_pair):
         """Generate and save visualization of colocalization results"""
-        if not self.output_folder:
+        if not self.output_folder or not self.save_images:
             return
 
         try:
@@ -435,18 +861,9 @@ class ColocalizationWorker(QThread):
             filepath_c1 = file_pair[0]  # Channel 1
             image_c1 = tifffile.imread(filepath_c1)
 
-            # Try to load channel 2 as well
-            try:
-                # filepath_c2 = file_pair[1]  # Channel 2
-                # image_c2 = tifffile.imread(filepath_c2)
-                has_c2 = True
-            except (FileNotFoundError, IndexError):
-                has_c2 = False
-
-            # Try to load channel 3 if available
-            has_c3 = False
-            if len(file_pair) > 2:
-                contextlib.suppress(FileNotFoundError, IndexError)
+            # Check if we have channel 2 and 3
+            has_c2 = len(file_pair) > 1
+            has_c3 = len(file_pair) > 2
 
             # Create output filename
             channels_str = "_".join(self.channel_names)
@@ -561,6 +978,323 @@ class ColocalizationWorker(QThread):
 
         return int(size)
 
+    def calculate_intensity_stats(self, intensity_image, mask):
+        """
+        Calculate intensity statistics for a masked region.
+
+        Args:
+            intensity_image: Raw intensity image
+            mask: Boolean mask defining the region
+
+        Returns:
+            dict: Dictionary with mean, median, std, max, min intensity
+        """
+        # Get intensity values within the mask
+        intensity_values = intensity_image[mask]
+
+        if len(intensity_values) == 0:
+            return {
+                "mean": 0.0,
+                "median": 0.0,
+                "std": 0.0,
+                "max": 0.0,
+                "min": 0.0,
+            }
+
+        stats = {
+            "mean": float(np.mean(intensity_values)),
+            "median": float(np.median(intensity_values)),
+            "std": float(np.std(intensity_values)),
+            "max": float(np.max(intensity_values)),
+            "min": float(np.min(intensity_values)),
+        }
+
+        return stats
+
+    def _process_individual_c2_labels(
+        self,
+        filename,
+        c1_label_id,
+        image_c1,
+        image_c2,
+        image_c3,
+        roi_sizes,
+        mask_roi,
+        image_c2_intensity,
+        image_c3_intensity,
+    ):
+        """
+        Process each C2 label individually, returning one row per C2 label.
+
+        Returns:
+            list of lists: One row per C2 label within the C1 ROI
+        """
+        # Get all unique Channel 2 objects in the ROI
+        c2_in_roi = image_c2 * mask_roi
+        c2_labels = np.unique(c2_in_roi)
+        c2_labels = c2_labels[c2_labels != 0]  # Remove background
+
+        if len(c2_labels) == 0:
+            # No c2 labels in this ROI, return empty list (no rows)
+            return []
+
+        rows = []
+        c1_size = (
+            roi_sizes.get(int(c1_label_id), 0) if self.get_sizes else None
+        )
+
+        for c2_label in sorted(c2_labels):
+            # Start row with filename, c1_label_id, c2_label_id
+            row = [filename, int(c1_label_id), int(c2_label)]
+
+            # Get mask for this specific Channel 2 object within the ROI
+            mask_c2_obj = (image_c2 == c2_label) & mask_roi
+
+            # Add size information if requested
+            if self.get_sizes:
+                row.append(c1_size)  # C1 ROI size (same for all c2 labels)
+
+                # C2 label size
+                c2_size = int(np.count_nonzero(mask_c2_obj))
+                row.append(c2_size)
+
+            # Handle C3 channel based on its type
+            if image_c3 is not None:
+                if self.channel3_is_labels:
+                    # C3 is labels: count unique C3 labels in this C2 label
+                    mask_c3 = image_c3 != 0
+                    mask_c3_in_c2 = mask_c2_obj & mask_c3
+                    c3_count = self.count_unique_nonzero(
+                        image_c3, mask_c3_in_c2
+                    )
+                    row.append(c3_count)
+
+                    # Add C3 size if requested
+                    if self.get_sizes:
+                        c3_size = int(np.count_nonzero(mask_c3_in_c2))
+                        row.append(c3_size)
+                else:
+                    # C3 is intensity: calculate intensity statistics in this C2 label
+                    intensity_img = (
+                        image_c3_intensity
+                        if image_c3_intensity is not None
+                        else image_c3
+                    )
+                    intensity_in_obj = intensity_img[mask_c2_obj]
+
+                    if len(intensity_in_obj) > 0:
+                        c3_mean = float(np.mean(intensity_in_obj))
+                        c3_median = float(np.median(intensity_in_obj))
+                        c3_std = float(np.std(intensity_in_obj))
+                    else:
+                        c3_mean = 0.0
+                        c3_median = 0.0
+                        c3_std = 0.0
+
+                    row.extend([c3_mean, c3_median, c3_std])
+
+            rows.append(row)
+
+        return rows
+
+    def calculate_individual_c2_intensities(
+        self, image_c2, intensity_c3, mask_roi
+    ):
+        """
+        Calculate individual Channel 3 intensity values for each Channel 2 label.
+
+        Args:
+            image_c2: Label image of Channel 2 (e.g., nuclei)
+            intensity_c3: Intensity image of Channel 3
+            mask_roi: Boolean mask for the ROI from Channel 1
+
+        Returns:
+            dict: Dictionary mapping c2_label_id -> intensity value (mean of c3 in that c2 label)
+        """
+        # Get all unique Channel 2 objects in the ROI
+        c2_in_roi = image_c2 * mask_roi
+        c2_labels = np.unique(c2_in_roi)
+        c2_labels = c2_labels[c2_labels != 0]  # Remove background
+
+        individual_values = {}
+        for c2_label in c2_labels:
+            # Get mask for this specific Channel 2 object within the ROI
+            mask_c2_obj = (image_c2 == c2_label) & mask_roi
+
+            # Get intensity values of Channel 3 in this Channel 2 object
+            intensity_in_obj = intensity_c3[mask_c2_obj]
+
+            if len(intensity_in_obj) > 0:
+                # Use mean intensity as the representative value for this c2 label
+                individual_values[int(c2_label)] = float(
+                    np.mean(intensity_in_obj)
+                )
+            else:
+                individual_values[int(c2_label)] = 0.0
+
+        return individual_values
+
+    def calculate_individual_c2_sizes(self, image_c2, mask_roi, image_c3=None):
+        """
+        Calculate individual sizes for each Channel 2 label within the ROI.
+
+        Args:
+            image_c2: Label image of Channel 2
+            mask_roi: Boolean mask for the ROI from Channel 1
+            image_c3: Optional Channel 3 image (if provided, calculate c3 size within each c2 label)
+
+        Returns:
+            dict: Dictionary mapping c2_label_id -> size (pixel count)
+        """
+        # Get all unique Channel 2 objects in the ROI
+        c2_in_roi = image_c2 * mask_roi
+        c2_labels = np.unique(c2_in_roi)
+        c2_labels = c2_labels[c2_labels != 0]  # Remove background
+
+        individual_sizes = {}
+        for c2_label in c2_labels:
+            # Get mask for this specific Channel 2 object within the ROI
+            mask_c2_obj = (image_c2 == c2_label) & mask_roi
+
+            if image_c3 is not None:
+                # Calculate c3 size within this c2 label
+                mask_with_c3 = mask_c2_obj & (image_c3 != 0)
+                size = int(np.count_nonzero(mask_with_c3))
+            else:
+                # Calculate c2 label size
+                size = int(np.count_nonzero(mask_c2_obj))
+
+            individual_sizes[int(c2_label)] = size
+
+        return individual_sizes
+
+    def count_positive_objects(
+        self,
+        image_c2,
+        intensity_c3,
+        mask_roi,
+        threshold_method="percentile",
+        threshold_value=75.0,
+    ):
+        """
+        Count Channel 2 objects that are positive for Channel 3 signal.
+
+        Args:
+            image_c2: Label image of Channel 2 (e.g., nuclei)
+            intensity_c3: Intensity image of Channel 3 (e.g., KI67)
+            mask_roi: Boolean mask for the ROI from Channel 1
+            threshold_method: 'percentile' or 'absolute'
+            threshold_value: Threshold value (0-100 for percentile, or absolute intensity)
+
+        Returns:
+            dict: Dictionary with counts and threshold info
+        """
+        # Get all unique Channel 2 objects in the ROI
+        c2_in_roi = image_c2 * mask_roi
+        c2_labels = np.unique(c2_in_roi)
+        c2_labels = c2_labels[c2_labels != 0]  # Remove background
+
+        if len(c2_labels) == 0:
+            return {
+                "total_c2_objects": 0,
+                "positive_c2_objects": 0,
+                "negative_c2_objects": 0,
+                "percent_positive": 0.0,
+                "threshold_used": 0.0,
+            }
+
+        # Calculate threshold
+        if threshold_method == "percentile":
+            # Calculate threshold from all Channel 3 intensity values within ROI where Channel 2 exists
+            mask_c2_in_roi = c2_in_roi > 0
+            intensity_in_c2 = intensity_c3[mask_c2_in_roi]
+            if len(intensity_in_c2) > 0:
+                threshold = float(
+                    np.percentile(intensity_in_c2, threshold_value)
+                )
+            else:
+                threshold = 0.0
+        else:  # absolute
+            threshold = threshold_value
+
+        # Count positive objects
+        positive_count = 0
+        for label_id in c2_labels:
+            # Get mask for this specific Channel 2 object
+            mask_c2_obj = (image_c2 == label_id) & mask_roi
+
+            # Get mean intensity of Channel 3 in this Channel 2 object
+            intensity_in_obj = intensity_c3[mask_c2_obj]
+            if len(intensity_in_obj) > 0:
+                mean_intensity = float(np.mean(intensity_in_obj))
+                if mean_intensity >= threshold:
+                    positive_count += 1
+
+        total_count = int(len(c2_labels))
+        negative_count = total_count - positive_count
+        percent_positive = (
+            (positive_count / total_count * 100) if total_count > 0 else 0.0
+        )
+
+        return {
+            "total_c2_objects": total_count,
+            "positive_c2_objects": positive_count,
+            "negative_c2_objects": negative_count,
+            "percent_positive": percent_positive,
+            "threshold_used": threshold,
+        }
+
+    def count_c2_positive_for_c3_labels(self, image_c2, image_c3, mask_roi):
+        """
+        Count Channel 2 objects that contain at least one Channel 3 object (label-based).
+
+        Args:
+            image_c2: Label image of Channel 2 (e.g., nuclei)
+            image_c3: Label image of Channel 3 (e.g., Ki67 spots)
+            mask_roi: Boolean mask for the ROI from Channel 1
+
+        Returns:
+            dict: Dictionary with positive/negative counts and percentage
+        """
+        # Get all unique Channel 2 objects in the ROI
+        c2_in_roi = image_c2 * mask_roi
+        c2_labels = np.unique(c2_in_roi)
+        c2_labels = c2_labels[c2_labels != 0]  # Remove background
+
+        if len(c2_labels) == 0:
+            return {
+                "total_c2_objects": 0,
+                "c2_positive_for_c3_count": 0,
+                "c2_negative_for_c3_count": 0,
+                "c2_percent_positive_for_c3": 0.0,
+            }
+
+        # Count how many C2 objects contain at least one C3 object
+        positive_count = 0
+        for c2_label in c2_labels:
+            # Get mask for this specific Channel 2 object
+            mask_c2_obj = (image_c2 == c2_label) & mask_roi
+
+            # Check if any C3 objects overlap with this C2 object
+            c3_in_c2 = image_c3[mask_c2_obj]
+            c3_labels_in_c2 = np.unique(c3_in_c2[c3_in_c2 != 0])
+
+            if len(c3_labels_in_c2) > 0:
+                positive_count += 1
+
+        total_count = int(len(c2_labels))
+        negative_count = total_count - positive_count
+        percent_positive = (
+            (positive_count / total_count * 100) if total_count > 0 else 0.0
+        )
+
+        return {
+            "total_c2_objects": total_count,
+            "c2_positive_for_c3_count": positive_count,
+            "c2_negative_for_c3_count": negative_count,
+            "c2_percent_positive_for_c3": percent_positive,
+        }
+
     def stop(self):
         """Request worker to stop processing"""
         self.stop_requested = True
@@ -726,6 +1460,12 @@ class ColocalizationAnalysisWidget(QWidget):
         self.ch1_folder = QLineEdit()
         self.ch1_pattern = QLineEdit()
         self.ch1_pattern.setPlaceholderText("*_labels.tif")
+        self.ch1_pattern.setToolTip(
+            "Glob pattern for matching files. Wildcards:\n"
+            "* = any characters (e.g., *_labels.tif)\n"
+            "? = single character (e.g., *_labels?.tif)\n"
+            "[seq] = character in sequence (e.g., *_labels[0-9]*.tif for _labels1, _labels23, etc.)"
+        )
         self.ch1_browse = QPushButton("Browse...")
         self.ch1_browse.clicked.connect(lambda: self.browse_folder(0))
 
@@ -741,6 +1481,12 @@ class ColocalizationAnalysisWidget(QWidget):
         self.ch2_folder = QLineEdit()
         self.ch2_pattern = QLineEdit()
         self.ch2_pattern.setPlaceholderText("*_labels.tif")
+        self.ch2_pattern.setToolTip(
+            "Glob pattern for matching files. Wildcards:\n"
+            "* = any characters (e.g., *_labels.tif)\n"
+            "? = single character (e.g., *_labels?.tif)\n"
+            "[seq] = character in sequence (e.g., *_labels[0-9]*.tif for _labels1, _labels23, etc.)"
+        )
         self.ch2_browse = QPushButton("Browse...")
         self.ch2_browse.clicked.connect(lambda: self.browse_folder(1))
 
@@ -754,8 +1500,15 @@ class ColocalizationAnalysisWidget(QWidget):
         # Channel 3 (optional)
         self.ch3_label = QLabel("Channel 3 (Optional):")
         self.ch3_folder = QLineEdit()
+        self.ch3_folder.textChanged.connect(lambda: self.update_ch3_controls())
         self.ch3_pattern = QLineEdit()
         self.ch3_pattern.setPlaceholderText("*_labels.tif")
+        self.ch3_pattern.setToolTip(
+            "Glob pattern for matching files. Wildcards:\n"
+            "* = any characters (e.g., *_labels.tif or *.tif for intensity)\n"
+            "? = single character (e.g., *_labels?.tif)\n"
+            "[seq] = character in sequence (e.g., *_labels[0-9]*.tif for _labels1, _labels23, etc.)"
+        )
         self.ch3_browse = QPushButton("Browse...")
         self.ch3_browse.clicked.connect(lambda: self.browse_folder(2))
 
@@ -771,26 +1524,61 @@ class ColocalizationAnalysisWidget(QWidget):
 
         # Get sizes option
         self.get_sizes_checkbox = QCheckBox("Calculate Region Sizes")
+        self.get_sizes_checkbox.toggled.connect(self._on_get_sizes_changed)
         options_layout.addRow(self.get_sizes_checkbox)
 
+        # Save images option
+        self.save_images_checkbox = QCheckBox("Save Output Images")
+        self.save_images_checkbox.setChecked(
+            False
+        )  # Default to not saving images
+        self.save_images_checkbox.setToolTip(
+            "Save visualization images showing colocalization results.\n"
+            "Uncheck to only generate CSV output (faster)."
+        )
+        options_layout.addRow(self.save_images_checkbox)
+
         # Size calculation method
         self.size_method_layout = QHBoxLayout()
         self.size_method_label = QLabel("Size Calculation Method:")
         self.size_method_median = QCheckBox("Median")
         self.size_method_median.setChecked(True)
+        self.size_method_median.setToolTip(
+            "Aggregate mode: One row per C1 ROI\n"
+            "Size = median size of C2 objects in this C1 ROI"
+        )
         self.size_method_sum = QCheckBox("Sum")
+        self.size_method_sum.setToolTip(
+            "Aggregate mode: One row per C1 ROI\n"
+            "Size = total size of all C2 objects in this C1 ROI"
+        )
+        self.size_method_individual = QCheckBox("Individual")
+        self.size_method_individual.setToolTip(
+            "Individual mode: One row per C2 object (not per C1 ROI)\n"
+            "Each C2 label gets its own row with individual stats.\n"
+            "\n"
+            "Note: Positive counting is disabled in Individual mode.\n"
+            "You get per-object C3 mean/median/std values instead."
+        )
 
         # Connect to make them mutually exclusive
         self.size_method_median.toggled.connect(
             lambda checked: (
-                self.size_method_sum.setChecked(not checked)
+                self._update_size_method_checkboxes("median", checked)
                 if checked
                 else None
             )
         )
         self.size_method_sum.toggled.connect(
             lambda checked: (
-                self.size_method_median.setChecked(not checked)
+                self._update_size_method_checkboxes("sum", checked)
+                if checked
+                else None
+            )
+        )
+        self.size_method_individual.toggled.connect(
+            lambda checked: (
+                self._update_size_method_checkboxes("individual", checked)
                 if checked
                 else None
             )
@@ -799,8 +1587,158 @@ class ColocalizationAnalysisWidget(QWidget):
         self.size_method_layout.addWidget(self.size_method_label)
         self.size_method_layout.addWidget(self.size_method_median)
         self.size_method_layout.addWidget(self.size_method_sum)
+        self.size_method_layout.addWidget(self.size_method_individual)
         options_layout.addRow(self.size_method_layout)
 
+        # Initially disable size method controls
+        self._update_size_method_controls_state()
+
+        # Channel 2 mode selection
+        self.ch2_is_labels_checkbox = QCheckBox(
+            "Channel 2 is Labels (uncheck for intensity)"
+        )
+        self.ch2_is_labels_checkbox.setChecked(True)
+        self.ch2_is_labels_checkbox.setToolTip(
+            "Check: C2 contains labeled objects (e.g., nuclei segmentation)\n"
+            "  → Counts C2 objects in C1 ROIs\n"
+            "\n"
+            "Uncheck: C2 contains intensity values (e.g., fluorescence)\n"
+            "  → Calculates C2 intensity statistics in C1 ROIs"
+        )
+        self.ch2_is_labels_checkbox.toggled.connect(self.on_ch2_mode_changed)
+        options_layout.addRow(self.ch2_is_labels_checkbox)
+
+        # Channel 3 mode selection
+        self.ch3_is_labels_checkbox = QCheckBox(
+            "Channel 3 is Labels (uncheck for intensity)"
+        )
+        self.ch3_is_labels_checkbox.setChecked(True)
+        self.ch3_is_labels_checkbox.setEnabled(
+            False
+        )  # Disabled until ch3 folder is set
+        self.ch3_is_labels_checkbox.setToolTip(
+            "COMMON MODES:\n"
+            "\n"
+            "C2=Labels + C3=Intensity (UNCHECKED):\n"
+            "  → Measures C3 intensity within C2 objects\n"
+            "  → Also measures C3 where C2 doesn't exist\n"
+            "  → Use Individual mode for per-C2-object stats\n"
+            "  → Use Aggregate + Count Positive to threshold C2 objects\n"
+            "\n"
+            "C2=Labels + C3=Labels (CHECKED):\n"
+            "  → Counts C3 objects within C2 objects\n"
+            "  → Counts C3 objects outside C2 but in C1\n"
+            "\n"
+            "C2=Intensity + C3=Intensity (UNCHECKED):\n"
+            "  → Measures both C2 and C3 intensity in C1 ROIs"
+        )
+        self.ch3_is_labels_checkbox.toggled.connect(self.on_ch3_mode_changed)
+        options_layout.addRow(self.ch3_is_labels_checkbox)
+
+        # Semantic to instance label conversion options
+        self.convert_c2_checkbox = QCheckBox(
+            "Convert C2 Semantic to Instance Labels"
+        )
+        self.convert_c2_checkbox.setChecked(False)
+        self.convert_c2_checkbox.setToolTip(
+            "Enable if C2 contains semantic labels (all objects have same value).\n"
+            "This will find connected components and assign unique labels to each object."
+        )
+        options_layout.addRow(self.convert_c2_checkbox)
+
+        self.convert_c3_checkbox = QCheckBox(
+            "Convert C3 Semantic to Instance Labels"
+        )
+        self.convert_c3_checkbox.setChecked(False)
+        self.convert_c3_checkbox.setToolTip(
+            "Enable if C3 contains semantic labels (all objects have same value).\n"
+            "This will find connected components and assign unique labels to each object."
+        )
+        self.convert_c3_checkbox.setEnabled(
+            False
+        )  # Disabled until ch3 folder is set
+        options_layout.addRow(self.convert_c3_checkbox)
+
+        # Count C2 positive for C3 (both labels)
+        self.count_c2_positive_checkbox = QCheckBox(
+            "Count C2 Objects Positive for C3 (both labels)"
+        )
+        self.count_c2_positive_checkbox.setChecked(False)
+        self.count_c2_positive_checkbox.setEnabled(False)
+        self.count_c2_positive_checkbox.setToolTip(
+            "When both C2 and C3 are labels, count how many C2 objects contain\n"
+            "at least one C3 object (binary: positive/negative per C2 object)."
+        )
+        self.ch3_is_labels_checkbox.toggled.connect(
+            self.update_c2_positive_state
+        )
+        self.ch2_is_labels_checkbox.toggled.connect(
+            self.update_c2_positive_state
+        )
+        options_layout.addRow(self.count_c2_positive_checkbox)
+
+        # Positive counting option (only for intensity mode)
+        self.count_positive_checkbox = QCheckBox(
+            "Count Positive C2 Objects (Aggregate mode only)"
+        )
+        self.count_positive_checkbox.setEnabled(False)  # Disabled initially
+        self.count_positive_checkbox.setToolTip(
+            "Available when: C2 is labels AND C3 is intensity AND aggregate mode\n"
+            "\n"
+            "Counts how many C2 objects are 'positive' for C3 signal based on\n"
+            "mean C3 intensity within each C2 object vs. threshold.\n"
+            "\n"
+            "Outputs: positive_count, negative_count, percent_positive, threshold\n"
+            "\n"
+            "Note: Individual mode provides per-object C3 stats instead of counts."
+        )
+        self.count_positive_checkbox.toggled.connect(
+            self.on_count_positive_changed
+        )
+        options_layout.addRow(self.count_positive_checkbox)
+
+        # Threshold method selection
+        self.threshold_layout = QHBoxLayout()
+        self.threshold_label = QLabel("Threshold Method:")
+        self.threshold_percentile = QCheckBox("Percentile")
+        self.threshold_percentile.setChecked(True)
+        self.threshold_absolute = QCheckBox("Absolute")
+        self.threshold_percentile.setEnabled(False)
+        self.threshold_absolute.setEnabled(False)
+
+        # Connect to make them mutually exclusive
+        self.threshold_percentile.toggled.connect(
+            lambda checked: (
+                self.threshold_absolute.setChecked(not checked)
+                if checked
+                else None
+            )
+        )
+        self.threshold_absolute.toggled.connect(
+            lambda checked: (
+                self.threshold_percentile.setChecked(not checked)
+                if checked
+                else None
+            )
+        )
+
+        self.threshold_layout.addWidget(self.threshold_label)
+        self.threshold_layout.addWidget(self.threshold_percentile)
+        self.threshold_layout.addWidget(self.threshold_absolute)
+        options_layout.addRow(self.threshold_layout)
+
+        # Threshold value input
+        self.threshold_value_layout = QHBoxLayout()
+        self.threshold_value_label = QLabel("Threshold Value:")
+        self.threshold_value_input = QLineEdit("75.0")
+        self.threshold_value_input.setPlaceholderText(
+            "e.g., 75 for 75th percentile"
+        )
+        self.threshold_value_input.setEnabled(False)
+        self.threshold_value_layout.addWidget(self.threshold_value_label)
+        self.threshold_value_layout.addWidget(self.threshold_value_input)
+        options_layout.addRow(self.threshold_value_layout)
+
         layout.addLayout(options_layout)
 
         # Output folder selection
@@ -898,6 +1836,27 @@ class ColocalizationAnalysisWidget(QWidget):
                 self.ch2_folder.setText(folder)
             elif channel_index == 2:
                 self.ch3_folder.setText(folder)
+                # Enable ch3 checkbox when folder is set
+                self.update_ch3_controls()
+
+    def update_ch3_controls(self):
+        """Enable/disable channel 3 controls based on whether folder is set"""
+        ch3_folder = self.ch3_folder.text().strip()
+        # Expand user path (~/...) and normalize the path
+        if ch3_folder:
+            ch3_folder = os.path.expanduser(ch3_folder)
+            ch3_folder = os.path.abspath(ch3_folder)
+        has_ch3 = bool(ch3_folder and os.path.isdir(ch3_folder))
+
+        # Enable/disable ch3 checkbox
+        self.ch3_is_labels_checkbox.setEnabled(has_ch3)
+
+        # Enable/disable ch3 semantic conversion checkbox
+        self.convert_c3_checkbox.setEnabled(has_ch3)
+
+        # Update positive counting state based on ch3 availability
+        self.update_positive_counting_state()
+        self.update_c2_positive_state()
 
     def browse_output(self):
         """Browse for output folder"""
@@ -910,6 +1869,82 @@ class ColocalizationAnalysisWidget(QWidget):
         if folder:
             self.output_folder.setText(folder)
 
+    def _update_size_method_checkboxes(self, selected, checked):
+        """Make size method checkboxes mutually exclusive"""
+        if not checked:
+            return
+        if selected == "median":
+            self.size_method_sum.setChecked(False)
+            self.size_method_individual.setChecked(False)
+        elif selected == "sum":
+            self.size_method_median.setChecked(False)
+            self.size_method_individual.setChecked(False)
+        elif selected == "individual":
+            self.size_method_median.setChecked(False)
+            self.size_method_sum.setChecked(False)
+
+    def _on_get_sizes_changed(self, checked):
+        """Enable/disable size method controls based on get_sizes checkbox"""
+        self._update_size_method_controls_state()
+
+    def _update_size_method_controls_state(self):
+        """Update the enabled state of size method controls"""
+        enabled = self.get_sizes_checkbox.isChecked()
+        self.size_method_label.setEnabled(enabled)
+        self.size_method_median.setEnabled(enabled)
+        self.size_method_sum.setEnabled(enabled)
+        self.size_method_individual.setEnabled(enabled)
+
+    def on_ch2_mode_changed(self, checked):
+        """Handle channel 2 mode change (labels vs intensity)"""
+        # When ch2 is intensity, positive counting doesn't apply
+        # Positive counting only works when ch2 is labels and ch3 is intensity
+        self.update_positive_counting_state()
+
+    def on_ch3_mode_changed(self, checked):
+        """Handle channel 3 mode change (labels vs intensity)"""
+        # Enable/disable positive counting based on mode
+        # Positive counting only available when ch2 is labels and ch3 is intensity
+        self.update_positive_counting_state()
+
+    def update_positive_counting_state(self):
+        """Update the state of positive counting controls based on ch2 and ch3 modes"""
+        # Positive counting only works when ch2 is labels and ch3 is intensity
+        ch2_is_labels = self.ch2_is_labels_checkbox.isChecked()
+        ch3_is_labels = self.ch3_is_labels_checkbox.isChecked()
+
+        # Enable positive counting only when ch2 is labels AND ch3 is intensity
+        can_count_positive = ch2_is_labels and not ch3_is_labels
+        self.count_positive_checkbox.setEnabled(can_count_positive)
+
+        if not can_count_positive:
+            # Disable positive counting if conditions aren't met
+            self.count_positive_checkbox.setChecked(False)
+
+    def on_count_positive_changed(self, checked):
+        """Handle positive counting checkbox change"""
+        # Enable/disable threshold controls
+        self.threshold_percentile.setEnabled(checked)
+        self.threshold_absolute.setEnabled(checked)
+        self.threshold_value_input.setEnabled(checked)
+
+    def update_c2_positive_state(self):
+        """Update the state of C2 positive for C3 counting based on ch2 and ch3 modes"""
+        # Get folder and mode states
+        ch3_folder = self.ch3_folder.text().strip()
+        has_ch3 = bool(ch3_folder and os.path.isdir(ch3_folder))
+
+        # C2 positive counting only works when BOTH ch2 and ch3 are labels
+        ch2_is_labels = self.ch2_is_labels_checkbox.isChecked()
+        ch3_is_labels = self.ch3_is_labels_checkbox.isChecked()
+
+        # Enable only when ch3 exists AND both are labels
+        can_count_c2_positive = has_ch3 and ch2_is_labels and ch3_is_labels
+        self.count_c2_positive_checkbox.setEnabled(can_count_c2_positive)
+
+        if not can_count_c2_positive:
+            self.count_c2_positive_checkbox.setChecked(False)
+
     def find_matching_files(self):
         """Find matching files across channels using the updated grouping function."""
         # Get channel folders and patterns
@@ -922,6 +1957,17 @@ class ColocalizationAnalysisWidget(QWidget):
         ch3_folder = self.ch3_folder.text().strip()
         ch3_pattern = self.ch3_pattern.text().strip() or "*_labels.tif"
 
+        # Expand user paths (~/...) and normalize paths
+        if ch1_folder:
+            ch1_folder = os.path.expanduser(ch1_folder)
+            ch1_folder = os.path.abspath(ch1_folder)
+        if ch2_folder:
+            ch2_folder = os.path.expanduser(ch2_folder)
+            ch2_folder = os.path.abspath(ch2_folder)
+        if ch3_folder:
+            ch3_folder = os.path.expanduser(ch3_folder)
+            ch3_folder = os.path.abspath(ch3_folder)
+
         # Validate required folders
         if not ch1_folder or not os.path.isdir(ch1_folder):
             self.status_label.setText(
@@ -941,12 +1987,36 @@ class ColocalizationAnalysisWidget(QWidget):
         ch1_files = sorted(glob.glob(os.path.join(ch1_folder, ch1_pattern)))
         ch2_files = sorted(glob.glob(os.path.join(ch2_folder, ch2_pattern)))
 
+        # Check if files were found
+        if not ch1_files:
+            self.status_label.setText(
+                f"No files matching pattern '{ch1_pattern}' found in Channel 1 folder"
+            )
+            self.match_label.setText("No matching files found")
+            self.analyze_button.setEnabled(False)
+            return
+
+        if not ch2_files:
+            self.status_label.setText(
+                f"No files matching pattern '{ch2_pattern}' found in Channel 2 folder"
+            )
+            self.match_label.setText("No matching files found")
+            self.analyze_button.setEnabled(False)
+            return
+
         # Check if third channel is provided
         use_ch3 = bool(ch3_folder and os.path.isdir(ch3_folder))
         if use_ch3:
             ch3_files = sorted(
                 glob.glob(os.path.join(ch3_folder, ch3_pattern))
             )
+            if not ch3_files:
+                self.status_label.setText(
+                    f"No files matching pattern '{ch3_pattern}' found in Channel 3 folder"
+                )
+                self.match_label.setText("No matching files found")
+                self.analyze_button.setEnabled(False)
+                return
         else:
             ch3_files = []
 
@@ -997,11 +2067,37 @@ class ColocalizationAnalysisWidget(QWidget):
 
         # Get settings
         get_sizes = self.get_sizes_checkbox.isChecked()
-        size_method = (
-            "median" if self.size_method_median.isChecked() else "sum"
-        )
+        save_images = self.save_images_checkbox.isChecked()
+
+        # Determine size method
+        if self.size_method_median.isChecked():
+            size_method = "median"
+        elif self.size_method_sum.isChecked():
+            size_method = "sum"
+        elif self.size_method_individual.isChecked():
+            size_method = "individual"
+        else:
+            size_method = "median"  # Default fallback
+
         output_folder = self.output_folder.text().strip()
 
+        # Get new settings for channel mode and positive counting
+        channel2_is_labels = self.ch2_is_labels_checkbox.isChecked()
+        channel3_is_labels = self.ch3_is_labels_checkbox.isChecked()
+        count_positive = self.count_positive_checkbox.isChecked()
+        threshold_method = (
+            "percentile"
+            if self.threshold_percentile.isChecked()
+            else "absolute"
+        )
+
+        # Get threshold value
+        try:
+            threshold_value = float(self.threshold_value_input.text())
+        except ValueError:
+            threshold_value = 75.0  # Default value
+            self.threshold_value_input.setText("75.0")
+
         # Create output folder if it doesn't exist and is specified
         if output_folder:
             try:
@@ -1032,13 +2128,54 @@ class ColocalizationAnalysisWidget(QWidget):
         self.analyze_button.setEnabled(False)
         self.cancel_button.setEnabled(True)
 
+        # Determine actual channel names based on folder names
+        # file_pairs contains tuples of (ch1_file, ch2_file) or (ch1_file, ch2_file, ch3_file)
+        ch1_folder = self.ch1_folder.text().strip()
+        ch2_folder = self.ch2_folder.text().strip()
+        ch3_folder = self.ch3_folder.text().strip()
+
+        # Expand paths for consistent handling
+        if ch1_folder:
+            ch1_folder = os.path.expanduser(ch1_folder)
+            ch1_folder = os.path.abspath(ch1_folder)
+        if ch2_folder:
+            ch2_folder = os.path.expanduser(ch2_folder)
+            ch2_folder = os.path.abspath(ch2_folder)
+        if ch3_folder:
+            ch3_folder = os.path.expanduser(ch3_folder)
+            ch3_folder = os.path.abspath(ch3_folder)
+
+        active_channel_names = [
+            os.path.basename(ch1_folder) if ch1_folder else "CH1",
+            os.path.basename(ch2_folder) if ch2_folder else "CH2",
+        ]
+
+        # Only add third channel if it exists
+        num_channels = len(self.file_pairs[0]) if self.file_pairs else 2
+        if num_channels == 3 and ch3_folder:
+            active_channel_names.append(os.path.basename(ch3_folder))
+
+        # Get conversion settings
+        convert_to_instances_c2 = self.convert_c2_checkbox.isChecked()
+        convert_to_instances_c3 = self.convert_c3_checkbox.isChecked()
+        count_c2_positive_for_c3 = self.count_c2_positive_checkbox.isChecked()
+
         # Create worker thread
         self.worker = ColocalizationWorker(
             self.file_pairs,
-            self.channel_names,
+            active_channel_names,
             get_sizes,
             size_method,
             output_folder,
+            channel2_is_labels,
+            channel3_is_labels,
+            count_positive,
+            threshold_method,
+            threshold_value,
+            save_images,
+            convert_to_instances_c2,
+            convert_to_instances_c3,
+            count_c2_positive_for_c3,
         )
 
         # Set thread count
@@ -1061,7 +2198,7 @@ class ColocalizationAnalysisWidget(QWidget):
         # Create results widget if needed
         if not self.results_widget:
             self.results_widget = ColocalizationResultsWidget(
-                self.viewer, self.channel_names
+                self.viewer, active_channel_names
             )
             self.viewer.window.add_dock_widget(
                 self.results_widget,