diff --git a/__pycache__/live_wire.cpython-310.pyc b/__pycache__/live_wire.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..d4bbfbe2e7917ef01fe041e1ec231ac761b3e6e3
Binary files /dev/null and b/__pycache__/live_wire.cpython-310.pyc differ
diff --git a/__pycache__/live_wire.cpython-312.pyc b/__pycache__/live_wire.cpython-312.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..8404b8c9798df68288011af419140c3a717972bb
Binary files /dev/null and b/__pycache__/live_wire.cpython-312.pyc differ
diff --git a/modules/__pycache__/advancedSettingsWidget.cpython-310.pyc b/modules/__pycache__/advancedSettingsWidget.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..7236b63a0e28e1ce6261fd142586150be4764839
Binary files /dev/null and b/modules/__pycache__/advancedSettingsWidget.cpython-310.pyc differ
diff --git a/modules/__pycache__/circleEditorGraphicsView.cpython-310.pyc b/modules/__pycache__/circleEditorGraphicsView.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..bf701598b85bf4473791e050917cea81199453b3
Binary files /dev/null and b/modules/__pycache__/circleEditorGraphicsView.cpython-310.pyc differ
diff --git a/modules/__pycache__/circleEditorWidget.cpython-310.pyc b/modules/__pycache__/circleEditorWidget.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..15985fdff54132c56f69d348837b1226bb88fa7d
Binary files /dev/null and b/modules/__pycache__/circleEditorWidget.cpython-310.pyc differ
diff --git a/modules/__pycache__/circle_edge_kernel.cpython-310.pyc b/modules/__pycache__/circle_edge_kernel.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..ee86e4ef66d3ce22af22a0a8c402d0374cd819a3
Binary files /dev/null and b/modules/__pycache__/circle_edge_kernel.cpython-310.pyc differ
diff --git a/modules/__pycache__/compute_cost_image.cpython-310.pyc b/modules/__pycache__/compute_cost_image.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..e755561615f0c50661393a5a4e40992dd71e5376
Binary files /dev/null and b/modules/__pycache__/compute_cost_image.cpython-310.pyc differ
diff --git a/modules/__pycache__/compute_disk_size.cpython-310.pyc b/modules/__pycache__/compute_disk_size.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..3383772746e94d58575ed8f6aec432708e812f0e
Binary files /dev/null and b/modules/__pycache__/compute_disk_size.cpython-310.pyc differ
diff --git a/modules/__pycache__/draggableCircleItem.cpython-310.pyc b/modules/__pycache__/draggableCircleItem.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..cf4a7fe166a02e860231ece0675a8968b64f8c61
Binary files /dev/null and b/modules/__pycache__/draggableCircleItem.cpython-310.pyc differ
diff --git a/modules/__pycache__/find_path.cpython-310.pyc b/modules/__pycache__/find_path.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..1b19c0f67ced5e24602bf64bf5ab8a63a1e30ea1
Binary files /dev/null and b/modules/__pycache__/find_path.cpython-310.pyc differ
diff --git a/modules/__pycache__/imageGraphicsView.cpython-310.pyc b/modules/__pycache__/imageGraphicsView.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..c7704ad73e53c223bbd845ff2e4e55ab3800bcc6
Binary files /dev/null and b/modules/__pycache__/imageGraphicsView.cpython-310.pyc differ
diff --git a/modules/__pycache__/labeledPointItem.cpython-310.pyc b/modules/__pycache__/labeledPointItem.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..6e7ffced62a9d9085d9c896338dd5bb829dd993b
Binary files /dev/null and b/modules/__pycache__/labeledPointItem.cpython-310.pyc differ
diff --git a/modules/__pycache__/load_image.cpython-310.pyc b/modules/__pycache__/load_image.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..4e867587aa08a678807752452ad7df52ecea128e
Binary files /dev/null and b/modules/__pycache__/load_image.cpython-310.pyc differ
diff --git a/modules/__pycache__/mainWindow.cpython-310.pyc b/modules/__pycache__/mainWindow.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..f202a4a0c3b8081acdb334f6de92fcb4b3009718
Binary files /dev/null and b/modules/__pycache__/mainWindow.cpython-310.pyc differ
diff --git a/modules/__pycache__/panZoomGraphicsView.cpython-310.pyc b/modules/__pycache__/panZoomGraphicsView.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..7c0393a74aa13a847ec4d0a2153cb6cfe456c9f3
Binary files /dev/null and b/modules/__pycache__/panZoomGraphicsView.cpython-310.pyc differ
diff --git a/modules/__pycache__/preprocess_image.cpython-310.pyc b/modules/__pycache__/preprocess_image.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..32055ff7fcd436221cc6f7f3293f35ea7f856c72
Binary files /dev/null and b/modules/__pycache__/preprocess_image.cpython-310.pyc differ
diff --git a/modules/advancedSettingsWidget.py b/modules/advancedSettingsWidget.py
new file mode 100644
index 0000000000000000000000000000000000000000..c11186352d905f3ca41d6f3abacdae3b0415ddba
--- /dev/null
+++ b/modules/advancedSettingsWidget.py
@@ -0,0 +1,158 @@
+from PyQt5.QtWidgets import (
+ QPushButton, QVBoxLayout, QWidget,
+ QSlider, QLabel, QGridLayout, QSizePolicy
+)
+from PyQt5.QtGui import QPixmap, QImage
+from PyQt5.QtCore import Qt
+import numpy as np
+
+class AdvancedSettingsWidget(QWidget):
+ """
+ Shows toggle rainbow, circle editor, line smoothing slider, contrast slider,
+ plus two image previews (contrasted-blurred and cost).
+ The images should maintain aspect ratio upon resize.
+ """
+ def __init__(self, main_window, parent=None):
+ super().__init__(parent)
+ self._main_window = main_window
+
+ self._last_cb_pix = None # store QPixmap for contrasted-blurred
+ self._last_cost_pix = None # store QPixmap for cost
+
+ main_layout = QVBoxLayout()
+ self.setLayout(main_layout)
+
+ # A small grid for controls
+ controls_layout = QGridLayout()
+
+ # 1) Rainbow toggle
+ self.btn_toggle_rainbow = QPushButton("Toggle Rainbow")
+ self.btn_toggle_rainbow.clicked.connect(self._on_toggle_rainbow)
+ controls_layout.addWidget(self.btn_toggle_rainbow, 0, 0)
+
+ # 2) Circle editor
+ self.btn_circle_editor = QPushButton("Calibrate Kernel Size")
+ self.btn_circle_editor.clicked.connect(self._main_window.open_circle_editor)
+ controls_layout.addWidget(self.btn_circle_editor, 0, 1)
+
+ # 3) Line smoothing slider + label
+ self._lab_smoothing = QLabel("Line smoothing (3)")
+ controls_layout.addWidget(self._lab_smoothing, 1, 0)
+ self.line_smoothing_slider = QSlider(Qt.Horizontal)
+ self.line_smoothing_slider.setRange(3, 51)
+ self.line_smoothing_slider.setValue(3)
+ self.line_smoothing_slider.valueChanged.connect(self._on_line_smoothing_slider)
+ controls_layout.addWidget(self.line_smoothing_slider, 1, 1)
+
+ # 4) Contrast slider + label
+ self._lab_contrast = QLabel("Contrast (0.01)")
+ controls_layout.addWidget(self._lab_contrast, 2, 0)
+ self.contrast_slider = QSlider(Qt.Horizontal)
+ self.contrast_slider.setRange(1, 20)
+ self.contrast_slider.setValue(1) # i.e. 0.01
+ self.contrast_slider.setSingleStep(1)
+ self.contrast_slider.valueChanged.connect(self._on_contrast_slider)
+ controls_layout.addWidget(self.contrast_slider, 2, 1)
+
+ main_layout.addLayout(controls_layout)
+
+ # We'll set a minimum width so that the main window expands
+ # rather than overlapping the image
+ self.setMinimumWidth(350)
+
+ # Now a vertical layout for the two images, each with a label above it
+ images_layout = QVBoxLayout()
+
+ # 1) Contrasted-blurred label + image
+ self.label_cb_title = QLabel("Contrasted Blurred Image")
+ self.label_cb_title.setAlignment(Qt.AlignCenter)
+ images_layout.addWidget(self.label_cb_title)
+
+ self.label_contrasted_blurred = QLabel()
+ self.label_contrasted_blurred.setAlignment(Qt.AlignCenter)
+ self.label_contrasted_blurred.setSizePolicy(QSizePolicy.Expanding, QSizePolicy.Expanding)
+ images_layout.addWidget(self.label_contrasted_blurred)
+
+ # 2) Cost image label + image
+ self.label_cost_title = QLabel("Current COST IMAGE")
+ self.label_cost_title.setAlignment(Qt.AlignCenter)
+ images_layout.addWidget(self.label_cost_title)
+
+ self.label_cost_image = QLabel()
+ self.label_cost_image.setAlignment(Qt.AlignCenter)
+ self.label_cost_image.setSizePolicy(QSizePolicy.Expanding, QSizePolicy.Expanding)
+ images_layout.addWidget(self.label_cost_image)
+
+ main_layout.addLayout(images_layout)
+
+ def showEvent(self, event):
+ """ When shown, ask parent to resize to accommodate. """
+ super().showEvent(event)
+ if self.parentWidget():
+ self.parentWidget().adjustSize()
+
+ def resizeEvent(self, event):
+ """
+ Keep the images at correct aspect ratio by re-scaling
+ our stored pixmaps to the new label sizes.
+ """
+ super().resizeEvent(event)
+ self._update_labels()
+
+ def _update_labels(self):
+ if self._last_cb_pix is not None:
+ scaled_cb = self._last_cb_pix.scaled(
+ self.label_contrasted_blurred.size(),
+ Qt.KeepAspectRatio,
+ Qt.SmoothTransformation
+ )
+ self.label_contrasted_blurred.setPixmap(scaled_cb)
+
+ if self._last_cost_pix is not None:
+ scaled_cost = self._last_cost_pix.scaled(
+ self.label_cost_image.size(),
+ Qt.KeepAspectRatio,
+ Qt.SmoothTransformation
+ )
+ self.label_cost_image.setPixmap(scaled_cost)
+
+ def _on_toggle_rainbow(self):
+ self._main_window.toggle_rainbow()
+
+ def _on_line_smoothing_slider(self, value):
+ self._lab_smoothing.setText(f"Line smoothing ({value})")
+ self._main_window.image_view.set_savgol_window_length(value)
+
+ def _on_contrast_slider(self, value):
+ clip_limit = value / 100.0
+ self._lab_contrast.setText(f"Contrast ({clip_limit:.2f})")
+ self._main_window.update_contrast(clip_limit)
+
+ def update_displays(self, contrasted_img_np, cost_img_np):
+ """
+ Called by main_window to refresh the two images in the advanced panel.
+ We'll store them as QPixmaps, then do the re-scale in _update_labels().
+ """
+ cb_pix = self._np_array_to_qpixmap(contrasted_img_np)
+ cost_pix = self._np_array_to_qpixmap(cost_img_np, normalize=True)
+
+ self._last_cb_pix = cb_pix
+ self._last_cost_pix = cost_pix
+ self._update_labels()
+
+ def _np_array_to_qpixmap(self, arr, normalize=False):
+ if arr is None:
+ return None
+ arr_ = arr.copy()
+ if normalize:
+ mn, mx = arr_.min(), arr_.max()
+ if abs(mx - mn) < 1e-12:
+ arr_[:] = 0
+ else:
+ arr_ = (arr_ - mn) / (mx - mn)
+ arr_ = np.clip(arr_, 0, 1)
+ arr_255 = (arr_ * 255).astype(np.uint8)
+
+ h, w = arr_255.shape
+ qimage = QImage(arr_255.data, w, h, w, QImage.Format_Grayscale8)
+ return QPixmap.fromImage(qimage)
diff --git a/modules/circleEditorGraphicsView.py b/modules/circleEditorGraphicsView.py
new file mode 100644
index 0000000000000000000000000000000000000000..d3a35173d58d3a4b0b09b558e0c187c12750e057
--- /dev/null
+++ b/modules/circleEditorGraphicsView.py
@@ -0,0 +1,49 @@
+from PyQt5.QtWidgets import QGraphicsView
+from panZoomGraphicsView import PanZoomGraphicsView
+from PyQt5.QtCore import Qt
+from draggableCircleItem import DraggableCircleItem
+
+# A specialized PanZoomGraphicsView for the circle editor
+class CircleEditorGraphicsView(PanZoomGraphicsView):
+ def __init__(self, circle_editor_widget, parent=None):
+ super().__init__(parent)
+ self._circle_editor_widget = circle_editor_widget
+
+ def mousePressEvent(self, event):
+ if event.button() == Qt.LeftButton:
+ # Check if user clicked on the circle item
+ clicked_item = self.itemAt(event.pos())
+ if clicked_item is not None:
+ # climb up parent chain
+ it = clicked_item
+ while it is not None and not hasattr(it, "boundingRect"):
+ it = it.parentItem()
+
+ if isinstance(it, DraggableCircleItem):
+ # Let normal item-dragging occur, no pan
+ return QGraphicsView.mousePressEvent(self, event)
+ super().mousePressEvent(event)
+
+ def wheelEvent(self, event):
+ """
+ If the mouse is hovering over the circle, we adjust the circle's radius
+ instead of zooming the image.
+ """
+ pos_in_widget = event.pos()
+ item_under = self.itemAt(pos_in_widget)
+ if item_under is not None:
+ it = item_under
+ while it is not None and not hasattr(it, "boundingRect"):
+ it = it.parentItem()
+
+ if isinstance(it, DraggableCircleItem):
+ delta = event.angleDelta().y()
+ step = 1 if delta > 0 else -1
+ old_r = it.radius()
+ new_r = max(1, old_r + step)
+ it.set_radius(new_r)
+ self._circle_editor_widget.update_slider_value(new_r)
+ event.accept()
+ return
+
+ super().wheelEvent(event)
diff --git a/modules/circleEditorWidget.py b/modules/circleEditorWidget.py
new file mode 100644
index 0000000000000000000000000000000000000000..aef3101258b59c15ecb1a9475204800be90ad02a
--- /dev/null
+++ b/modules/circleEditorWidget.py
@@ -0,0 +1,91 @@
+from PyQt5.QtWidgets import (
+ QGraphicsScene, QGraphicsPixmapItem, QPushButton,
+ QHBoxLayout, QVBoxLayout, QWidget, QSlider, QLabel
+)
+from PyQt5.QtGui import QFont
+from PyQt5.QtCore import Qt, QRectF, QSize
+from circleEditorGraphicsView import CircleEditorGraphicsView
+from draggableCircleItem import DraggableCircleItem
+
+class CircleEditorWidget(QWidget):
+ def __init__(self, pixmap, init_radius=20, done_callback=None, parent=None):
+ super().__init__(parent)
+ self._pixmap = pixmap
+ self._done_callback = done_callback
+ self._init_radius = init_radius
+
+ layout = QVBoxLayout(self)
+ self.setLayout(layout)
+
+ #
+ # 1) ADD A CENTERED LABEL ABOVE THE IMAGE, WITH BIGGER FONT
+ #
+ label_instructions = QLabel("Scale the dot to be of the size of your ridge")
+ label_instructions.setAlignment(Qt.AlignCenter)
+ big_font = QFont("Arial", 20)
+ big_font.setBold(True)
+ label_instructions.setFont(big_font)
+ layout.addWidget(label_instructions)
+
+ #
+ # 2) THE SPECIALIZED GRAPHICS VIEW THAT SHOWS THE IMAGE
+ #
+ self._graphics_view = CircleEditorGraphicsView(circle_editor_widget=self)
+ self._scene = QGraphicsScene(self)
+ self._graphics_view.setScene(self._scene)
+ layout.addWidget(self._graphics_view)
+
+ # Show the image
+ self._image_item = QGraphicsPixmapItem(self._pixmap)
+ self._scene.addItem(self._image_item)
+
+ # Put circle in center
+ cx = self._pixmap.width() / 2
+ cy = self._pixmap.height() / 2
+ self._circle_item = DraggableCircleItem(cx, cy, radius=self._init_radius, color=Qt.red)
+ self._scene.addItem(self._circle_item)
+
+ # Fit in view
+ self._graphics_view.setSceneRect(QRectF(self._pixmap.rect()))
+ self._graphics_view.fitInView(self._image_item, Qt.KeepAspectRatio)
+
+ #
+ # 3) CONTROLS BELOW
+ #
+ bottom_layout = QHBoxLayout()
+ layout.addLayout(bottom_layout)
+
+ # label + slider
+ self._lbl_size = QLabel(f"size ({self._init_radius})")
+ bottom_layout.addWidget(self._lbl_size)
+
+ self._slider = QSlider(Qt.Horizontal)
+ self._slider.setRange(1, 200)
+ self._slider.setValue(self._init_radius)
+ bottom_layout.addWidget(self._slider)
+
+ # done button
+ self._btn_done = QPushButton("Done")
+ bottom_layout.addWidget(self._btn_done)
+
+ # Connect signals
+ self._slider.valueChanged.connect(self._on_slider_changed)
+ self._btn_done.clicked.connect(self._on_done_clicked)
+
+ def _on_slider_changed(self, value):
+ self._circle_item.set_radius(value)
+ self._lbl_size.setText(f"size ({value})")
+
+ def _on_done_clicked(self):
+ final_radius = self._circle_item.radius()
+ if self._done_callback is not None:
+ self._done_callback(final_radius)
+
+ def update_slider_value(self, new_radius):
+ self._slider.blockSignals(True)
+ self._slider.setValue(new_radius)
+ self._slider.blockSignals(False)
+ self._lbl_size.setText(f"size ({new_radius})")
+
+ def sizeHint(self):
+ return QSize(800, 600)
diff --git a/modules/circle_edge_kernel.py b/modules/circle_edge_kernel.py
new file mode 100644
index 0000000000000000000000000000000000000000..4270fba0203817b67eaf879f011078a981e94ca6
--- /dev/null
+++ b/modules/circle_edge_kernel.py
@@ -0,0 +1,38 @@
+import numpy as np
+
+def circle_edge_kernel(k_size=5, radius=None):
+ """
+ Create a k_size x k_size array whose values increase linearly
+ from 0 at the center to 1 at the circle boundary (radius).
+
+ Parameters
+ ----------
+ k_size : int
+ The size (width and height) of the kernel array.
+ radius : float, optional
+ The circle's radius. By default, set to (k_size-1)/2.
+
+ Returns
+ -------
+ kernel : 2D numpy array of shape (k_size, k_size)
+ The circle-edge-weighted kernel.
+ """
+ if radius is None:
+ # By default, let the radius be half the kernel size
+ radius = (k_size - 1) / 2
+
+ # Create an empty kernel
+ kernel = np.zeros((k_size, k_size), dtype=float)
+
+ # Coordinates of the center
+ center = radius # same as (k_size-1)/2 if radius is default
+
+ # Fill the kernel
+ for y in range(k_size):
+ for x in range(k_size):
+ dist = np.sqrt((x - center)**2 + (y - center)**2)
+ if dist <= radius:
+ # Weight = distance / radius => 0 at center, 1 at boundary
+ kernel[y, x] = dist / radius
+
+ return kernel
\ No newline at end of file
diff --git a/modules/compute_cost_image.py b/modules/compute_cost_image.py
new file mode 100644
index 0000000000000000000000000000000000000000..3c88768bd224759efb98c0b5f369f7a878ceca04
--- /dev/null
+++ b/modules/compute_cost_image.py
@@ -0,0 +1,41 @@
+from skimage.feature import canny
+from scipy.signal import convolve2d
+from compute_disk_size import compute_disk_size
+from load_image import load_image
+from preprocess_image import preprocess_image
+from circle_edge_kernel import circle_edge_kernel
+import numpy as np
+
+def compute_cost_image(path: str, user_radius: int, sigma: int = 3, clip_limit: float = 0.01) -> np.ndarray:
+ """
+ Compute the cost image for a given image path, user radius, and optional parameters.
+
+ Args:
+ path: The path to the image file.
+ user_radius: The radius of the disk.
+ sigma: The standard deviation for Gaussian smoothing.
+ clip_limit: The limit for contrasting the image.
+
+ Returns:
+ The cost image as a NumPy array.
+ """
+ disk_size = compute_disk_size(user_radius)
+
+ # Load image
+ image = load_image(path)
+
+ # Apply smoothing
+ smoothed_img = preprocess_image(image, sigma=sigma, clip_limit=clip_limit)
+
+ # Apply Canny edge detection
+ canny_img = canny(smoothed_img)
+
+ # Perform disk convolution
+ binary_img = canny_img
+ kernel = circle_edge_kernel(k_size=disk_size)
+ convolved = convolve2d(binary_img, kernel, mode='same', boundary='fill')
+
+ # Create cost image
+ cost_img = (convolved.max() - convolved)**4 # Invert edges: higher cost where edges are stronger
+
+ return cost_img
\ No newline at end of file
diff --git a/modules/compute_disk_size.py b/modules/compute_disk_size.py
new file mode 100644
index 0000000000000000000000000000000000000000..c4cb24d66a0d4c84587d84eff3bc342606d5e74f
--- /dev/null
+++ b/modules/compute_disk_size.py
@@ -0,0 +1,4 @@
+import numpy as np
+
+def compute_disk_size(user_radius, upscale_factor=1.2):
+ return int(np.ceil(upscale_factor * 2 * user_radius + 1) // 2 * 2 + 1)
\ No newline at end of file
diff --git a/modules/downscale.py b/modules/downscale.py
new file mode 100644
index 0000000000000000000000000000000000000000..2387e26bfd74f79861cfbcc3c58336f22554be65
--- /dev/null
+++ b/modules/downscale.py
@@ -0,0 +1,30 @@
+import cv2
+
+# Currently not implemented
+def downscale(img, points, scale_percent):
+ """
+ Downsample `img` to `scale_percent` size and scale the given points accordingly.
+ Returns (downsampled_img, (scaled_seed, scaled_target)).
+ """
+ if scale_percent == 100:
+ return img, (tuple(points[0]), tuple(points[1]))
+ else:
+ # Compute new dimensions
+ width = int(img.shape[1] * scale_percent / 100)
+ height = int(img.shape[0] * scale_percent / 100)
+ new_dimensions = (width, height)
+
+ # Downsample
+ downsampled_img = cv2.resize(img, new_dimensions, interpolation=cv2.INTER_AREA)
+
+ # Scaling factors
+ scale_x = width / img.shape[1]
+ scale_y = height / img.shape[0]
+
+ # Scale the points (x, y)
+ seed_xy = tuple(points[0])
+ target_xy = tuple(points[1])
+ scaled_seed_xy = (int(seed_xy[0] * scale_x), int(seed_xy[1] * scale_y))
+ scaled_target_xy = (int(target_xy[0] * scale_x), int(target_xy[1] * scale_y))
+
+ return downsampled_img, (scaled_seed_xy, scaled_target_xy)
\ No newline at end of file
diff --git a/modules/draggableCircleItem.py b/modules/draggableCircleItem.py
new file mode 100644
index 0000000000000000000000000000000000000000..8d0e6fa15fe823b36a54bd89a96b1ca4fe0b49b2
--- /dev/null
+++ b/modules/draggableCircleItem.py
@@ -0,0 +1,33 @@
+from PyQt5.QtWidgets import QGraphicsEllipseItem
+from PyQt5.QtGui import QPen, QBrush
+from PyQt5.QtCore import Qt
+
+class DraggableCircleItem(QGraphicsEllipseItem):
+ def __init__(self, x, y, radius=20, color=Qt.red, parent=None):
+ super().__init__(0, 0, 2*radius, 2*radius, parent)
+ self._r = radius
+
+ pen = QPen(color)
+ brush = QBrush(color)
+ self.setPen(pen)
+ self.setBrush(brush)
+
+ # Enable item-based dragging
+ self.setFlags(QGraphicsEllipseItem.ItemIsMovable |
+ QGraphicsEllipseItem.ItemIsSelectable |
+ QGraphicsEllipseItem.ItemSendsScenePositionChanges)
+
+ # Position so that (x, y) is the center
+ self.setPos(x - radius, y - radius)
+
+ def set_radius(self, r):
+ old_center = self.sceneBoundingRect().center()
+ self._r = r
+ self.setRect(0, 0, 2*r, 2*r)
+ new_center = self.sceneBoundingRect().center()
+ diff_x = old_center.x() - new_center.x()
+ diff_y = old_center.y() - new_center.y()
+ self.moveBy(diff_x, diff_y)
+
+ def radius(self):
+ return self._r
\ No newline at end of file
diff --git a/modules/export_path.py b/modules/export_path.py
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/modules/find_path.py b/modules/find_path.py
new file mode 100644
index 0000000000000000000000000000000000000000..4ae9c7fa5afc77a9a29cf214956f55fbfd225fc7
--- /dev/null
+++ b/modules/find_path.py
@@ -0,0 +1,30 @@
+from skimage.graph import route_through_array
+
+def find_path(cost_image: "numpy.ndarray", points: list) -> list:
+ """
+ Find the optimal path through a cost image between two points.
+
+ Parameters:
+ cost_image (numpy.ndarray): A 2D array representing the cost of traversing each pixel.
+ points (list): A list containing two tuples, each representing the
+ (row, column) coordinates of the seed and target points.
+
+ Returns:
+ list: A list of (row, column) tuples representing the path from the seed to the target point.
+
+ Raises:
+ ValueError: If the points list does not contain exactly two points.
+ """
+ if len(points) != 2:
+ raise ValueError("Points should be a list of 2 points: seed and target.")
+
+ seed_rc, target_rc = points
+
+ path_rc, cost = route_through_array(
+ cost_image,
+ start=seed_rc,
+ end=target_rc,
+ fully_connected=True
+ )
+
+ return path_rc
diff --git a/modules/imageGraphicsView.py b/modules/imageGraphicsView.py
new file mode 100644
index 0000000000000000000000000000000000000000..387f74191613a4fff6aed65374157dab17d6599c
--- /dev/null
+++ b/modules/imageGraphicsView.py
@@ -0,0 +1,512 @@
+import math
+from scipy.signal import savgol_filter
+from PyQt5.QtWidgets import QGraphicsScene, QGraphicsPixmapItem
+from PyQt5.QtGui import QPixmap, QColor
+from PyQt5.QtCore import Qt, QRectF
+import numpy as np
+from panZoomGraphicsView import PanZoomGraphicsView
+from labeledPointItem import LabeledPointItem
+from find_path import find_path
+
+
+class ImageGraphicsView(PanZoomGraphicsView):
+ """
+ A custom QGraphicsView for displaying and interacting with an image.
+
+ This class extends PanZoomGraphicsView to provide additional functionality
+ for loading images, adding labeled anchor points, and computing paths
+ between points based on a cost image.
+
+ Attributes:
+ scene (QGraphicsScene): The graphics scene for displaying items.
+ image_item (QGraphicsPixmapItem): The item for displaying the loaded image.
+ anchor_points (list): List of tuples representing anchor points (x, y).
+ point_items (list): List of LabeledPointItem objects for anchor points.
+ full_path_points (list): List of QGraphicsEllipseItems representing the path.
+ _full_path_xy (list): List of coordinates for the entire path.
+ dot_radius (int): Radius of the anchor points.
+ path_radius (int): Radius of the path points.
+ radius_cost_image (int): Radius for lowering cost in the cost image.
+ _img_w (int): Width of the loaded image.
+ _img_h (int): Height of the loaded image.
+ _mouse_pressed (bool): Indicates if the mouse is pressed.
+ _press_view_pos (QPoint): Position of the mouse press event.
+ _drag_threshold (int): Threshold for detecting drag events.
+ _was_dragging (bool): Indicates if a drag event occurred.
+ _dragging_idx (int): Index of the point being dragged.
+ _drag_offset (tuple): Offset for dragging points.
+ _drag_counter (int): Counter for drag events.
+ cost_image_original (np.ndarray): Original cost image.
+ cost_image (np.ndarray): Current cost image.
+ _rainbow_enabled (bool): Indicates if rainbow coloring is enabled.
+ _savgol_window_length (int): Window length for Savitzky-Golay smoothing.
+ """
+
+ def __init__(self, parent=None):
+ super().__init__(parent)
+ self.scene = QGraphicsScene(self)
+ self.setScene(self.scene)
+
+ # Image display
+ self.image_item = QGraphicsPixmapItem()
+ self.scene.addItem(self.image_item)
+
+ self.anchor_points = []
+ self.point_items = []
+ self.full_path_points = []
+ self._full_path_xy = []
+
+ self.dot_radius = 4
+ self.path_radius = 1
+ self.radius_cost_image = 2
+ self._img_w = 0
+ self._img_h = 0
+
+ self._mouse_pressed = False
+ self._press_view_pos = None
+ self._drag_threshold = 5
+ self._was_dragging = False
+ self._dragging_idx = None
+ self._drag_offset = (0, 0)
+ self._drag_counter = 0
+
+ # Cost images
+ self.cost_image_original = None
+ self.cost_image = None
+
+ # Rainbow toggle => start with OFF
+ self._rainbow_enabled = False
+
+ # Smoothing parameters
+ self._savgol_window_length = 7
+
+ def set_rainbow_enabled(self, enabled: bool):
+ """Enable rainbow coloring of the path."""
+ self._rainbow_enabled = enabled
+ self._rebuild_full_path()
+
+ def toggle_rainbow(self):
+ """Toggle rainbow coloring of the path."""
+ self._rainbow_enabled = not self._rainbow_enabled
+ self._rebuild_full_path()
+
+ def set_savgol_window_length(self, wlen: int):
+ """Set the window length for Savitzky-Golay smoothing."""
+ wlen = max(3, wlen)
+ if wlen % 2 == 0:
+ wlen += 1
+ self._savgol_window_length = wlen
+
+ self._rebuild_full_path()
+
+ # --------------------------------------------------------------------
+ # LOADING
+ # --------------------------------------------------------------------
+ def load_image(self, path: str):
+ """Load an image from a file path."""
+ pixmap = QPixmap(path)
+ if not pixmap.isNull():
+ self.image_item.setPixmap(pixmap)
+ self.setSceneRect(QRectF(pixmap.rect()))
+
+ self._img_w = pixmap.width()
+ self._img_h = pixmap.height()
+
+ self._clear_all_points()
+ self.resetTransform()
+ self.fitInView(self.image_item, Qt.KeepAspectRatio)
+
+ # By default, add S/E
+ s_x, s_y = 0.15 * self._img_w, 0.5 * self._img_h
+ e_x, e_y = 0.85 * self._img_w, 0.5 * self._img_h
+ self._insert_anchor_point(-1, s_x, s_y, label="S", removable=False, z_val=100, radius=6)
+ self._insert_anchor_point(-1, e_x, e_y, label="E", removable=False, z_val=100, radius=6)
+
+ # --------------------------------------------------------------------
+ # ANCHOR POINTS
+ # --------------------------------------------------------------------
+ def _insert_anchor_point(self, idx, x: float, y: float, label="", removable=True,
+ z_val=0, radius=4):
+ """Insert an anchor point at a specific index."""
+ x_clamped = self._clamp(x, radius, self._img_w - radius)
+ y_clamped = self._clamp(y, radius, self._img_h - radius)
+
+ if idx < 0:
+ # Insert before E if there's at least 2 anchors
+ if len(self.anchor_points) >= 2:
+ idx = len(self.anchor_points) - 1
+ else:
+ idx = len(self.anchor_points)
+
+ self.anchor_points.insert(idx, (x_clamped, y_clamped))
+ color = Qt.green if label in ("S", "E") else Qt.red
+ item = LabeledPointItem(x_clamped, y_clamped,
+ label=label, radius=radius, color=color,
+ removable=removable, z_value=z_val)
+ self.point_items.insert(idx, item)
+ self.scene.addItem(item)
+
+ def _add_guide_point(self, x, y):
+ """Add a guide point to the path."""
+ x_clamped = self._clamp(x, self.dot_radius, self._img_w - self.dot_radius)
+ y_clamped = self._clamp(y, self.dot_radius, self._img_h - self.dot_radius)
+
+ self._revert_cost_to_original()
+
+ if not self._full_path_xy:
+ self._insert_anchor_point(-1, x_clamped, y_clamped,
+ label="", removable=True, z_val=1, radius=self.dot_radius)
+ else:
+ self._insert_anchor_between_subpath(x_clamped, y_clamped)
+
+ self._apply_all_guide_points_to_cost()
+ self._rebuild_full_path()
+
+ def _insert_anchor_between_subpath(self, x_new: float, y_new: float ):
+ """Insert an anchor point between existing anchor points."""
+ # If somehow we have no path yet
+ if not self._full_path_xy:
+ self._insert_anchor_point(-1, x_new, y_new)
+ return
+
+ # Find nearest point in the current full path
+ best_idx = None
+ best_d2 = float('inf')
+ for i, (px, py) in enumerate(self._full_path_xy):
+ dx = px - x_new
+ dy = py - y_new
+ d2 = dx*dx + dy*dy
+ if d2 < best_d2:
+ best_d2 = d2
+ best_idx = i
+
+ if best_idx is None:
+ self._insert_anchor_point(-1, x_new, y_new)
+ return
+
+ def approx_equal(xa, ya, xb, yb, tol=1e-3):
+ """Check if two points are approximately equal."""
+ return (abs(xa - xb) < tol) and (abs(ya - yb) < tol)
+
+ def is_anchor(coord):
+ """Check if a point is an anchor point."""
+ cx, cy = coord
+ for (ax, ay) in self.anchor_points:
+ if approx_equal(ax, ay, cx, cy):
+ return True
+ return False
+
+ # Walk left
+ left_anchor_pt = None
+ i_l = best_idx
+ while i_l >= 0:
+ px, py = self._full_path_xy[i_l]
+ if is_anchor((px, py)):
+ left_anchor_pt = (px, py)
+ break
+ i_l -= 1
+
+ # Walk right
+ right_anchor_pt = None
+ i_r = best_idx
+ while i_r < len(self._full_path_xy):
+ px, py = self._full_path_xy[i_r]
+ if is_anchor((px, py)):
+ right_anchor_pt = (px, py)
+ break
+ i_r += 1
+
+ # If we can't find distinct anchors on left & right,
+ # just insert before E.
+ if not left_anchor_pt or not right_anchor_pt:
+ self._insert_anchor_point(-1, x_new, y_new)
+ return
+ if left_anchor_pt == right_anchor_pt:
+ self._insert_anchor_point(-1, x_new, y_new)
+ return
+
+ # Convert anchor coords -> anchor_points indices
+ left_idx = None
+ right_idx = None
+ for i, (ax, ay) in enumerate(self.anchor_points):
+ if approx_equal(ax, ay, left_anchor_pt[0], left_anchor_pt[1]):
+ left_idx = i
+ if approx_equal(ax, ay, right_anchor_pt[0], right_anchor_pt[1]):
+ right_idx = i
+
+ if left_idx is None or right_idx is None:
+ self._insert_anchor_point(-1, x_new, y_new)
+ return
+
+ # Insert between them
+ if left_idx < right_idx:
+ insert_idx = left_idx + 1
+ else:
+ insert_idx = right_idx + 1
+
+ self._insert_anchor_point(insert_idx, x_new, y_new, label="", removable=True,
+ z_val=1, radius=self.dot_radius)
+
+ # --------------------------------------------------------------------
+ # COST IMAGE
+ # --------------------------------------------------------------------
+ def _revert_cost_to_original(self):
+ if self.cost_image_original is not None:
+ self.cost_image = self.cost_image_original.copy()
+
+ def _apply_all_guide_points_to_cost(self):
+ if self.cost_image is None:
+ return
+ for i, (ax, ay) in enumerate(self.anchor_points):
+ if self.point_items[i].is_removable():
+ self._lower_cost_in_circle(ax, ay, self.radius_cost_image)
+
+ def _lower_cost_in_circle(self, x_f: float, y_f: float, radius: int):
+ """Lower the cost in a circle centered at (x_f, y_f)."""
+ if self.cost_image is None:
+ return
+ h, w = self.cost_image.shape
+ row_c = int(round(y_f))
+ col_c = int(round(x_f))
+ if not (0 <= row_c < h and 0 <= col_c < w):
+ return
+ global_min = self.cost_image.min()
+ r_s = max(0, row_c - radius)
+ r_e = min(h, row_c + radius + 1)
+ c_s = max(0, col_c - radius)
+ c_e = min(w, col_c + radius + 1)
+ for rr in range(r_s, r_e):
+ for cc in range(c_s, c_e):
+ dist = math.sqrt((rr - row_c)**2 + (cc - col_c)**2)
+ if dist <= radius:
+ self.cost_image[rr, cc] = global_min
+
+ # --------------------------------------------------------------------
+ # PATH BUILDING
+ # --------------------------------------------------------------------
+ def _rebuild_full_path(self):
+ """Rebuild the full path based on the anchor points."""
+ for item in self.full_path_points:
+ self.scene.removeItem(item)
+ self.full_path_points.clear()
+ self._full_path_xy.clear()
+
+ if len(self.anchor_points) < 2 or self.cost_image is None:
+ return
+
+ big_xy = []
+ for i in range(len(self.anchor_points) - 1):
+ xA, yA = self.anchor_points[i]
+ xB, yB = self.anchor_points[i + 1]
+ sub_xy = self._compute_subpath_xy(xA, yA, xB, yB)
+ if i == 0:
+ big_xy.extend(sub_xy)
+ else:
+ if len(sub_xy) > 1:
+ big_xy.extend(sub_xy[1:])
+
+ if len(big_xy) >= self._savgol_window_length:
+ arr_xy = np.array(big_xy)
+ smoothed = savgol_filter(
+ arr_xy,
+ window_length=self._savgol_window_length,
+ polyorder=2,
+ axis=0
+ )
+ big_xy = smoothed.tolist()
+
+ self._full_path_xy = big_xy[:]
+
+ n_points = len(big_xy)
+ for i, (px, py) in enumerate(big_xy):
+ fraction = i / (n_points - 1) if n_points > 1 else 0
+ color = Qt.red
+ if self._rainbow_enabled:
+ color = self._rainbow_color(fraction)
+
+ path_item = LabeledPointItem(px, py, label="",
+ radius=self.path_radius,
+ color=color,
+ removable=False,
+ z_value=0)
+ self.full_path_points.append(path_item)
+ self.scene.addItem(path_item)
+
+ # Keep anchor labels on top
+ for p_item in self.point_items:
+ if p_item._text_item:
+ p_item.setZValue(100)
+
+ def _compute_subpath_xy(self, xA: float, yA: float, xB: float, yB: float):
+ """Compute a subpath between two points."""
+ if self.cost_image is None:
+ return []
+ h, w = self.cost_image.shape
+ rA, cA = int(round(yA)), int(round(xA))
+ rB, cB = int(round(yB)), int(round(xB))
+ rA = max(0, min(rA, h - 1))
+ cA = max(0, min(cA, w - 1))
+ rB = max(0, min(rB, h - 1))
+ cB = max(0, min(cB, w - 1))
+ try:
+ path_rc = find_path(self.cost_image, [(rA, cA), (rB, cB)])
+ except ValueError as e:
+ print("Error in find_path:", e)
+ return []
+ # Convert from (row, col) to (x, y)
+ return [(c, r) for (r, c) in path_rc]
+
+ def _rainbow_color(self, fraction: float):
+ """Get a rainbow color."""
+ hue = int(300 * fraction)
+ saturation = 255
+ value = 255
+ return QColor.fromHsv(hue, saturation, value)
+
+ # --------------------------------------------------------------------
+ # MOUSE EVENTS
+ # --------------------------------------------------------------------
+ def mouse_press_event(self, event):
+ if event.button() == Qt.LeftButton:
+ self._mouse_pressed = True
+ self._was_dragging = False
+ self._press_view_pos = event.pos()
+
+ idx = self._find_item_near(event.pos(), threshold=10)
+ if idx is not None:
+ self._dragging_idx = idx
+ self._drag_counter = 0
+ scene_pos = self.mapToScene(event.pos())
+ px, py = self.point_items[idx].get_pos()
+ self._drag_offset = (scene_pos.x() - px, scene_pos.y() - py)
+ self.setCursor(Qt.ClosedHandCursor)
+ return
+
+ elif event.button() == Qt.RightButton:
+ self._remove_point_by_click(event.pos())
+
+ super().mouse_press_event(event)
+
+ def mouse_move_event(self, event):
+ """Handle mouse move events for dragging a point or dragging the view"""
+ if self._dragging_idx is not None:
+ scene_pos = self.mapToScene(event.pos())
+ x_new = scene_pos.x() - self._drag_offset[0]
+ y_new = scene_pos.y() - self._drag_offset[1]
+
+ r = self.point_items[self._dragging_idx]._r
+ x_clamped = self._clamp(x_new, r, self._img_w - r)
+ y_clamped = self._clamp(y_new, r, self._img_h - r)
+ self.point_items[self._dragging_idx].set_pos(x_clamped, y_clamped)
+
+ self._drag_counter += 1
+ # Update path every 4 moves
+ if self._drag_counter >= 4:
+ self._drag_counter = 0
+ self._revert_cost_to_original()
+ self._apply_all_guide_points_to_cost()
+ self.anchor_points[self._dragging_idx] = (x_clamped, y_clamped)
+ self._rebuild_full_path()
+ else:
+ if self._mouse_pressed and (event.buttons() & Qt.LeftButton):
+ dist = (event.pos() - self._press_view_pos).manhattanLength()
+ if dist > self._drag_threshold:
+ self._was_dragging = True
+
+ super().mouse_move_event(event)
+
+ def mouse_release_event(self, event):
+ super().mouse_release_event(event)
+ if event.button() == Qt.LeftButton and self._mouse_pressed:
+ self._mouse_pressed = False
+ self.setCursor(Qt.ArrowCursor)
+
+ if self._dragging_idx is not None:
+ idx = self._dragging_idx
+ self._dragging_idx = None
+ self._drag_offset = (0, 0)
+ newX, newY = self.point_items[idx].get_pos()
+ self.anchor_points[idx] = (newX, newY)
+ self._revert_cost_to_original()
+ self._apply_all_guide_points_to_cost()
+ self._rebuild_full_path()
+ else:
+ # No drag => add point
+ if not self._was_dragging:
+ scene_pos = self.mapToScene(event.pos())
+ x, y = scene_pos.x(), scene_pos.y()
+ self._add_guide_point(x, y)
+
+ self._was_dragging = False
+
+ def _remove_point_by_click(self, view_pos: "QPoint"):
+ idx = self._find_item_near(view_pos, threshold=10)
+ if idx is None:
+ return
+ if not self.point_items[idx].is_removable():
+ return
+
+ self.scene.removeItem(self.point_items[idx])
+ self.point_items.pop(idx)
+ self.anchor_points.pop(idx)
+
+ self._revert_cost_to_original()
+ self._apply_all_guide_points_to_cost()
+ self._rebuild_full_path()
+
+ def _find_item_near(self, view_pos: "QPoint", threshold=10):
+ scene_pos = self.mapToScene(view_pos)
+ x_click, y_click = scene_pos.x(), scene_pos.y()
+
+ closest_idx = None
+ min_dist = float('inf')
+ for i, itm in enumerate(self.point_items):
+ d = itm.distance_to(x_click, y_click)
+ if d < min_dist:
+ min_dist = d
+ closest_idx = i
+ if closest_idx is not None and min_dist <= threshold:
+ return closest_idx
+ return None
+
+ # --------------------------------------------------------------------
+ # UTILS
+ # --------------------------------------------------------------------
+ def _clamp(self, val, mn, mx):
+ return max(mn, min(val, mx))
+
+ def _clear_all_points(self):
+
+ for it in self.point_items:
+ self.scene.removeItem(it)
+ self.point_items.clear()
+ self.anchor_points.clear()
+
+ for p in self.full_path_points:
+ self.scene.removeItem(p)
+ self.full_path_points.clear()
+ self._full_path_xy.clear()
+
+ def clear_guide_points(self):
+ """Clear all guide points."""
+ i = 0
+ while i < len(self.anchor_points):
+ if self.point_items[i].is_removable():
+ self.scene.removeItem(self.point_items[i])
+ del self.point_items[i]
+ del self.anchor_points[i]
+ else:
+ i += 1
+
+ for it in self.full_path_points:
+ self.scene.removeItem(it)
+ self.full_path_points.clear()
+ self._full_path_xy.clear()
+
+ self._revert_cost_to_original()
+ self._apply_all_guide_points_to_cost()
+ self._rebuild_full_path()
+
+ def get_full_path_xy(self):
+ """Returns the entire path as a list of (x, y) coordinates."""
+ return self._full_path_xy
diff --git a/modules/labeledPointItem.py b/modules/labeledPointItem.py
new file mode 100644
index 0000000000000000000000000000000000000000..861a95fb402306c1696c85be33903595c1235b9b
--- /dev/null
+++ b/modules/labeledPointItem.py
@@ -0,0 +1,97 @@
+import math
+from PyQt5.QtWidgets import QGraphicsEllipseItem, QGraphicsTextItem
+from PyQt5.QtGui import QPen, QBrush, QColor, QFont
+from PyQt5.QtCore import Qt
+
+
+class LabeledPointItem(QGraphicsEllipseItem):
+ """
+ A QGraphicsEllipseItem subclass that represents a labeled point in a 2D space.
+
+ This class creates a circular point.
+ The point can be customized with different colors, sizes, and labels, and can
+ be marked as removable.
+
+ Attributes:
+ x (float): The x-coordinate of the point.
+ y (float): The y-coordinate of the point.
+ label (str): The label text for the point.
+ radius (int): The radius of the point.
+ color (QColor): The color of the point.
+ removable (bool): Indicates if the point can be removed.
+ z_value (float): The z-value of the point for stacking order.
+ parent (QGraphicsItem): The parent QGraphicsItem, if any.
+ """
+
+ def __init__(self, x: float, y: float, label: str ="", radius:int =4,
+ color=Qt.red, removable=True, z_value=0, parent=None):
+ super().__init__(0, 0, 2*radius, 2*radius, parent)
+ self._x = x
+ self._y = y
+ self._r = radius
+ self._removable = removable
+
+ pen = QPen(color)
+ brush = QBrush(color)
+ self.setPen(pen)
+ self.setBrush(brush)
+ self.setZValue(z_value)
+
+ self._text_item = None
+ if label:
+ self._text_item = QGraphicsTextItem(self)
+ self._text_item.setPlainText(label)
+ self._text_item.setDefaultTextColor(QColor("black"))
+ font = QFont("Arial", 14)
+ font.setBold(True)
+ self._text_item.setFont(font)
+ self._scale_text_to_fit()
+
+ self.set_pos(x, y)
+
+ def _scale_text_to_fit(self):
+ """Scales the text to fit inside the circle."""
+ if not self._text_item:
+ return
+ self._text_item.setScale(1.0)
+ circle_diam = 2 * self._r
+ raw_rect = self._text_item.boundingRect()
+ text_w = raw_rect.width()
+ text_h = raw_rect.height()
+ if text_w > circle_diam or text_h > circle_diam:
+ scale_factor = min(circle_diam / text_w, circle_diam / text_h)
+ self._text_item.setScale(scale_factor)
+ self._center_label()
+
+ def _center_label(self):
+ """Centers the text inside the circle."""
+ if not self._text_item:
+ return
+ ellipse_w = 2 * self._r
+ ellipse_h = 2 * self._r
+ raw_rect = self._text_item.boundingRect()
+ scale_factor = self._text_item.scale()
+ scaled_w = raw_rect.width() * scale_factor
+ scaled_h = raw_rect.height() * scale_factor
+ tx = (ellipse_w - scaled_w) * 0.5
+ ty = (ellipse_h - scaled_h) * 0.5
+ self._text_item.setPos(tx, ty)
+
+ def set_pos(self, x, y):
+ """Positions the circle so its center is at (x, y)."""
+ self._x = x
+ self._y = y
+ self.setPos(x - self._r, y - self._r)
+
+ def get_pos(self):
+ """Returns the (x, y) coordinates of the center of the circle."""
+ return (self._x, self._y)
+
+ def distance_to(self, x_other, y_other):
+ """Returns the Euclidean distance from the center
+ of the circle to another circle."""
+ return math.sqrt((self._x - x_other)**2 + (self._y - y_other)**2)
+
+ def is_removable(self):
+ """Returns True if the point is removable, False otherwise."""
+ return self._removable
diff --git a/modules/load_image.py b/modules/load_image.py
new file mode 100644
index 0000000000000000000000000000000000000000..d2fa7eadf42d6098f0cc14bee50a4f015d9d0240
--- /dev/null
+++ b/modules/load_image.py
@@ -0,0 +1,13 @@
+import cv2
+
+def load_image(path: str) -> "numpy.ndarray":
+ """
+ Loads an image from the specified file path in grayscale mode.
+
+ Args:
+ path (str): The file path to the image.
+
+ Returns:
+ numpy.ndarray: The loaded grayscale image.
+ """
+ return cv2.imread(path, cv2.IMREAD_GRAYSCALE)
\ No newline at end of file
diff --git a/modules/main.py b/modules/main.py
new file mode 100644
index 0000000000000000000000000000000000000000..e1f76dbb1deeac59c4ef8488c7ab23d0788dd750
--- /dev/null
+++ b/modules/main.py
@@ -0,0 +1,13 @@
+import sys
+from PyQt5.QtWidgets import QApplication
+from mainWindow import MainWindow
+
+def main():
+ app = QApplication(sys.argv)
+ window = MainWindow()
+ window.show()
+ sys.exit(app.exec_())
+
+
+if __name__ == "__main__":
+ main()
\ No newline at end of file
diff --git a/modules/mainWindow.py b/modules/mainWindow.py
new file mode 100644
index 0000000000000000000000000000000000000000..45ca7e8716b4b891b32614580df2e6b3d2240e32
--- /dev/null
+++ b/modules/mainWindow.py
@@ -0,0 +1,251 @@
+import math
+import numpy as np
+from scipy.signal import savgol_filter
+from PyQt5.QtWidgets import (
+ QMainWindow, QPushButton, QHBoxLayout,
+ QVBoxLayout, QWidget, QFileDialog
+)
+from PyQt5.QtGui import QPixmap, QImage
+from compute_cost_image import compute_cost_image
+from preprocess_image import preprocess_image
+from advancedSettingsWidget import AdvancedSettingsWidget
+from imageGraphicsView import ImageGraphicsView
+from circleEditorWidget import CircleEditorWidget
+
+class MainWindow(QMainWindow):
+ def __init__(self):
+ super().__init__()
+ self.setWindowTitle("Test GUI")
+
+ self._last_loaded_pixmap = None
+ self._circle_calibrated_radius = 6
+ self._last_loaded_file_path = None
+
+ # For the contrast slider
+ self._current_clip_limit = 0.01
+
+ # Outer widget + layout
+ self._main_widget = QWidget()
+ self._main_layout = QHBoxLayout(self._main_widget)
+
+ # The "left" part: container for the image area + its controls
+ self._left_panel = QVBoxLayout()
+
+ # We'll make a container widget for the left panel, so we can set stretches:
+ self._left_container = QWidget()
+ self._left_container.setLayout(self._left_panel)
+
+ # Now we add them to the main layout with 70%:30% ratio
+ self._main_layout.addWidget(self._left_container, 7) # 70%
+
+ # We haven't added the advanced widget yet, but we'll do so with ratio=3 => 30%
+ self._advanced_widget = AdvancedSettingsWidget(self)
+ # Hide it initially
+ self._advanced_widget.hide()
+ self._main_layout.addWidget(self._advanced_widget, 3)
+
+ self.setCentralWidget(self._main_widget)
+
+ # The image view
+ self.image_view = ImageGraphicsView()
+ self._left_panel.addWidget(self.image_view)
+
+ # Button row
+ btn_layout = QHBoxLayout()
+ self.btn_load_image = QPushButton("Load Image")
+ self.btn_load_image.clicked.connect(self.load_image)
+ btn_layout.addWidget(self.btn_load_image)
+
+ self.btn_export_path = QPushButton("Export Path")
+ self.btn_export_path.clicked.connect(self.export_path)
+ btn_layout.addWidget(self.btn_export_path)
+
+ self.btn_clear_points = QPushButton("Clear Points")
+ self.btn_clear_points.clicked.connect(self.clear_points)
+ btn_layout.addWidget(self.btn_clear_points)
+
+ # "Advanced Settings" toggle
+ self.btn_advanced = QPushButton("Advanced Settings")
+ self.btn_advanced.setCheckable(True)
+ self.btn_advanced.clicked.connect(self._toggle_advanced_settings)
+ btn_layout.addWidget(self.btn_advanced)
+
+ self._left_panel.addLayout(btn_layout)
+
+ self.resize(1000, 600)
+ self._old_central_widget = None
+ self._editor = None
+
+ def _toggle_advanced_settings(self, checked):
+ if checked:
+ self._advanced_widget.show()
+ else:
+ self._advanced_widget.hide()
+ # Force re-layout
+ self.adjustSize()
+
+ def open_circle_editor(self):
+ """ Replace central widget with circle editor. """
+ if not self._last_loaded_pixmap:
+ print("No image loaded yet! Cannot open circle editor.")
+ return
+
+ old_widget = self.takeCentralWidget()
+ self._old_central_widget = old_widget
+
+ init_radius = self._circle_calibrated_radius
+ editor = CircleEditorWidget(
+ pixmap=self._last_loaded_pixmap,
+ init_radius=init_radius,
+ done_callback=self._on_circle_editor_done
+ )
+ self._editor = editor
+ self.setCentralWidget(editor)
+
+ def _on_circle_editor_done(self, final_radius):
+ self._circle_calibrated_radius = final_radius
+ print(f"Circle Editor done. Radius = {final_radius}")
+
+ if self._last_loaded_file_path:
+ cost_img = compute_cost_image(
+ self._last_loaded_file_path,
+ self._circle_calibrated_radius,
+ clip_limit=self._current_clip_limit
+ )
+ self.image_view.cost_image_original = cost_img
+ self.image_view.cost_image = cost_img.copy()
+ self.image_view._apply_all_guide_points_to_cost()
+ self.image_view._rebuild_full_path()
+ self._update_advanced_images()
+
+ editor_widget = self.takeCentralWidget()
+ if editor_widget is not None:
+ editor_widget.setParent(None)
+
+ if self._old_central_widget is not None:
+ self.setCentralWidget(self._old_central_widget)
+ self._old_central_widget = None
+
+ if self._editor is not None:
+ self._editor.deleteLater()
+ self._editor = None
+
+ def toggle_rainbow(self):
+ self.image_view.toggle_rainbow()
+
+ def load_image(self):
+ options = QFileDialog.Options()
+ file_path, _ = QFileDialog.getOpenFileName(
+ self, "Open Image", "",
+ "Images (*.png *.jpg *.jpeg *.bmp *.tif)",
+ options=options
+ )
+ if file_path:
+ self.image_view.load_image(file_path)
+
+ cost_img = compute_cost_image(
+ file_path,
+ self._circle_calibrated_radius,
+ clip_limit=self._current_clip_limit
+ )
+ self.image_view.cost_image_original = cost_img
+ self.image_view.cost_image = cost_img.copy()
+
+ pm = QPixmap(file_path)
+ if not pm.isNull():
+ self._last_loaded_pixmap = pm
+
+ self._last_loaded_file_path = file_path
+ self._update_advanced_images()
+
+ def update_contrast(self, clip_limit):
+ self._current_clip_limit = clip_limit
+ if self._last_loaded_file_path:
+ cost_img = compute_cost_image(
+ self._last_loaded_file_path,
+ self._circle_calibrated_radius,
+ clip_limit=clip_limit
+ )
+ self.image_view.cost_image_original = cost_img
+ self.image_view.cost_image = cost_img.copy()
+ self.image_view._apply_all_guide_points_to_cost()
+ self.image_view._rebuild_full_path()
+
+ self._update_advanced_images()
+
+ def _update_advanced_images(self):
+ if not self._last_loaded_pixmap:
+ return
+ pm_np = self._qpixmap_to_gray_float(self._last_loaded_pixmap)
+ contrasted_blurred = preprocess_image(
+ pm_np,
+ sigma=3,
+ clip_limit=self._current_clip_limit
+ )
+ cost_img_np = self.image_view.cost_image
+ self._advanced_widget.update_displays(contrasted_blurred, cost_img_np)
+
+ def _qpixmap_to_gray_float(self, qpix):
+ img = qpix.toImage()
+ img = img.convertToFormat(QImage.Format_ARGB32)
+ ptr = img.bits()
+ ptr.setsize(img.byteCount())
+ arr = np.frombuffer(ptr, np.uint8).reshape((img.height(), img.width(), 4))
+ rgb = arr[..., :3].astype(np.float32)
+ gray = rgb.mean(axis=2) / 255.0
+ return gray
+
+ def export_path(self):
+ """
+ Exports the path as a CSV in the format: x, y, TYPE,
+ ensuring that each anchor influences exactly one path point.
+ """
+ full_xy = self.image_view.get_full_path_xy()
+ if not full_xy:
+ print("No path to export.")
+ return
+
+ # We'll consider each anchor point as "USER-PLACED".
+ # But unlike a distance-threshold approach, we assign each anchor
+ # to exactly one closest path point.
+ anchor_points = self.image_view.anchor_points
+
+ # For each anchor, find the index of the closest path point
+ user_placed_indices = set()
+ for ax, ay in anchor_points:
+ min_dist = float('inf')
+ closest_idx = None
+ for i, (px, py) in enumerate(full_xy):
+ dist = math.hypot(px - ax, py - ay)
+ if dist < min_dist:
+ min_dist = dist
+ closest_idx = i
+ if closest_idx is not None:
+ user_placed_indices.add(closest_idx)
+
+ # Ask user for the CSV filename
+ options = QFileDialog.Options()
+ file_path, _ = QFileDialog.getSaveFileName(
+ self, "Export Path", "",
+ "CSV Files (*.csv);;All Files (*)",
+ options=options
+ )
+ if not file_path:
+ return
+
+ import csv
+ with open(file_path, 'w', newline='') as csvfile:
+ writer = csv.writer(csvfile)
+ writer.writerow(["x", "y", "TYPE"])
+
+ for i, (x, y) in enumerate(full_xy):
+ ptype = "USER-PLACED" if i in user_placed_indices else "PATH"
+ writer.writerow([x, y, ptype])
+
+ print(f"Exported path with {len(full_xy)} points to {file_path}")
+
+ def clear_points(self):
+ self.image_view.clear_guide_points()
+
+ def closeEvent(self, event):
+ super().closeEvent(event)
\ No newline at end of file
diff --git a/modules/panZoomGraphicsView.py b/modules/panZoomGraphicsView.py
new file mode 100644
index 0000000000000000000000000000000000000000..ed9947fb58a554164c297602241df7824049fe68
--- /dev/null
+++ b/modules/panZoomGraphicsView.py
@@ -0,0 +1,59 @@
+from PyQt5.QtWidgets import QGraphicsView, QSizePolicy
+from PyQt5.QtCore import Qt
+
+class PanZoomGraphicsView(QGraphicsView):
+ """
+ A QGraphicsView subclass that supports panning and zooming with the mouse.
+
+ Attributes:
+ _panning (bool): Indicates whether panning is currently active.
+ _pan_start (QPoint): The starting point of the panning action.
+
+ Methods:
+ wheel_event(event): Zooms in or out based on the mouse wheel movement.
+ mouse_press_event(event): Starts panning if the left mouse button is pressed.
+ mouse_move_event(event): Translates the scene if panning is active.
+ mouse_release_event(event): Ends panning when the left mouse button is released.
+ """
+ def __init__(self, parent=None):
+ super().__init__(parent)
+ self.setDragMode(QGraphicsView.NoDrag) # We'll handle panning manually
+ self.setTransformationAnchor(QGraphicsView.AnchorUnderMouse)
+ self._panning = False
+ self._pan_start = None
+
+ # Expands layout
+ self.setSizePolicy(QSizePolicy.Expanding, QSizePolicy.Expanding)
+
+ def wheel_event(self, event):
+ """ Zoom in/out with mouse wheel. """
+ zoom_in_factor = 1.25
+ zoom_out_factor = 1 / zoom_in_factor
+ if event.angleDelta().y() > 0:
+ self.scale(zoom_in_factor, zoom_in_factor)
+ else:
+ self.scale(zoom_out_factor, zoom_out_factor)
+ event.accept()
+
+ def mouse_press_event(self, event):
+ """ If left button: Start panning (unless overridden). """
+ if event.button() == Qt.LeftButton:
+ self._panning = True
+ self._pan_start = event.pos()
+ self.setCursor(Qt.ClosedHandCursor)
+ super().mouse_press_event(event)
+
+ def mouse_move_event(self, event):
+ """ If panning, translate the scene. """
+ if self._panning and self._pan_start is not None:
+ delta = event.pos() - self._pan_start
+ self._pan_start = event.pos()
+ self.translate(delta.x(), delta.y())
+ super().mouse_move_event(event)
+
+ def mouse_release_event(self, event):
+ """ End panning. """
+ if event.button() == Qt.LeftButton:
+ self._panning = False
+ self.setCursor(Qt.ArrowCursor)
+ super().mouse_release_event(event)
diff --git a/modules/preprocess_image.py b/modules/preprocess_image.py
new file mode 100644
index 0000000000000000000000000000000000000000..351988f1b00389c592d35e1913a0ca1221cb7224
--- /dev/null
+++ b/modules/preprocess_image.py
@@ -0,0 +1,23 @@
+from skimage.filters import gaussian
+from skimage import exposure
+
+
+def preprocess_image(image: "np.ndarray", sigma: int = 3, clip_limit: float = 0.01) -> "np.ndarray":
+ """
+ Preprocess the input image by applying histogram equalization and Gaussian smoothing.
+
+ Args:
+ image: (ndarray): Input image to be processed.
+ sigma: (float, optional): Standard deviation for Gaussian kernel. Default is 3.
+ clip_limit: (float, optional): Clipping limit for contrast enhancement. Default is 0.01.
+ Returns:
+ ndarray: The preprocessed image.
+ """
+ # Applies histogram equalization to enhance contrast
+ image_contrasted = exposure.equalize_adapthist(
+ image, clip_limit=clip_limit)
+
+ # Applies smoothing
+ smoothed_img = gaussian(image_contrasted, sigma=sigma)
+
+ return smoothed_img