arrayview 0.27.0__py3-none-any.whl → 0.27.2__py3-none-any.whl

arrayview/_routes_rendering.py

@@ -25,6 +25,7 @@ from arrayview._render import (
 )
 from arrayview._session import HEAVY_OP_LIMIT_BYTES, SESSIONS
 from arrayview._synthetic_mri import (
+    _qmri_adjust_vmin_vmax,
     qmri_display_slice,
     render_qmri_mosaic_rgba,
     synthetic_qmri_slice,
@@ -256,6 +257,8 @@ def register_rendering_routes(app, *, get_session_or_404) -> None:
                     vmin_override=vmin_override,
                     vmax_override=vmax_override,
                 )
+                if qmri_role:
+                    vmin, vmax = _qmri_adjust_vmin_vmax(vmin, vmax, qmri_role)
         render_ms = (time.perf_counter() - render_t0) * 1000.0
         encode_t0 = time.perf_counter()
         img = _pil_image().fromarray(rgba[:, :, :3], mode="RGB")

arrayview/_routes_websocket.py

@@ -31,6 +31,7 @@ from arrayview._session import (
 )
 from arrayview._vectorfield import _compute_vfield_arrows
 from arrayview._synthetic_mri import (
+    _qmri_adjust_vmin_vmax,
     qmri_display_slice,
     render_qmri_mosaic_rgba,
     synthetic_qmri_slice,
@@ -319,6 +320,8 @@ def register_websocket_routes(app) -> None:
                         vmin_override=vmin_override,
                         vmax_override=vmax_override,
                     )
+                    if qmri_role:
+                        vmin, vmax = _qmri_adjust_vmin_vmax(vmin, vmax, qmri_role)
 
                     overlay_sid = msg.get("overlay_sid")
                     overlay_colors = msg.get("overlay_colors")

arrayview/_stdio_server.py

@@ -52,6 +52,7 @@ from arrayview._render import (
 )
 from arrayview._session import SESSIONS, Session
 from arrayview._synthetic_mri import (
+    _qmri_adjust_vmin_vmax,
     qmri_display_slice,
     render_qmri_mosaic_rgba,
     synthetic_qmri_slice,
@@ -638,6 +639,8 @@ def _build_slice_payload(msg: dict) -> bytes:
             session, raw, complex_mode, dr, log_scale,
             vmin_override=vmin_override, vmax_override=vmax_override,
         )
+        if qmri_role:
+            vmin, vmax = _qmri_adjust_vmin_vmax(vmin, vmax, qmri_role)
 
         # Overlay compositing
         overlay_sid = msg.get("overlay_sid")

arrayview/_synthetic_mri.py

@@ -48,7 +48,8 @@ def normalize_relaxation_ms(data: np.ndarray, role: str | None = None) -> np.nda
     if finite.size == 0:
         return arr
     central = float(np.nanmedian(finite))
-    if 0.0 < central <= 5.0:
+    seconds_threshold = 10.0 if role == "t1" else 5.0
+    if 0.0 < central <= seconds_threshold:
         return (arr * np.float32(1000.0)).astype(np.float32)
     return arr
 
@@ -145,6 +146,21 @@ def synthetic_window(data: np.ndarray, contrast: str) -> tuple[float, float]:
     return vmin, vmax
 
 
+def _qmri_adjust_vmin_vmax(vmin: float, vmax: float, role: str) -> tuple[float, float]:
+    """Adjust vmin/vmax for display according to qMRI role conventions."""
+    if not role:
+        return vmin, vmax
+    if role in ("t1", "t2"):
+        return 0.0, float(int(np.ceil(vmax)))
+    if role == "pd":
+        return 0.0, vmax
+    if role == "b1":
+        return round(float(vmin), 2), round(float(vmax), 2)
+    if role == "phase":
+        return -3.14, 3.14
+    return vmin, vmax
+
+
 def render_qmri_mosaic_rgba(
     session,
     dim_x: int,
@@ -174,6 +190,7 @@ def render_qmri_mosaic_rgba(
     else:
         vmin = float(np.percentile(all_data, 1))
         vmax = float(np.percentile(all_data, 99))
+    vmin, vmax = _qmri_adjust_vmin_vmax(vmin, vmax, role)
     grid, _, _ = _build_mosaic_grid(
         all_data, session.shape[dim_z], cols=mosaic_cols
     ) if mosaic_cols is not None else _build_mosaic_grid(all_data, session.shape[dim_z])