rahix/openflexure-microscope-server
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Do blur in PIL

Browse source
This commit is contained in:
Julian Stirling 2025-10-26 14:59:13 +00:00
parent 94c9a0e7dd
commit b8bdb430b7
1 changed files with 16 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

23
src/openflexure_microscope_server/things/camera/simulation.py
View file

@ -15,8 +15,7 @@ import time
import cv2 import cv2
import numpy as np import numpy as np
from PIL import Image from PIL import Image, ImageFilter
from scipy.ndimage import gaussian_filter
import labthings_fastapi as lt import labthings_fastapi as lt
@ -209,11 +208,9 @@ class SimulatedCamera(BaseCamera):
canvas = self.canvas if self._show_sample else self.blank_canvas canvas = self.canvas if self._show_sample else self.blank_canvas
# Use npx to make each 1d index list 3D # Use npx to make each 1d index list 3D
focused_image = canvas[np.ix_(x_indices, y_indices, z_indices)] focused_image = canvas[np.ix_(x_indices, y_indices, z_indices)]
image = gaussian_filter(
focused_image, image = fast_pil_blur(focused_image, sigma=np.abs(pos[2]) / 5)
sigma=np.abs(pos[2]) / 5,
axes=(0, 1),
)
if image.shape != self.shape: if image.shape != self.shape:
raise ValueError( raise ValueError(
f"Image shape {image.shape} does not match intended shape {self.shape}" f"Image shape {image.shape} does not match intended shape {self.shape}"
@ -393,3 +390,15 @@ class SimulatedCamera(BaseCamera):
def manual_camera_settings(self) -> list[PropertyControl]: def manual_camera_settings(self) -> list[PropertyControl]:
"""The camera settings to expose as property controls in the settings panel.""" """The camera settings to expose as property controls in the settings panel."""
return [property_control_for(self, "noise_level", label="Noise Level")] return [property_control_for(self, "noise_level", label="Noise Level")]
def fast_pil_blur(array: np.ndarray, sigma: float) -> np.ndarray:
"""Apply Gaussian blur using PIL (faster than scipy)."""
if sigma < 0.5:
return array # no visible blur needed
img_pil = Image.fromarray(array.astype(np.uint8))
img_pil = img_pil.filter(ImageFilter.GaussianBlur(radius=sigma))
# Convert back to NumPy array
return np.array(img_pil, dtype=array.dtype)