arrayview 0.47.1__tar.gz → 0.47.4__tar.gz

{arrayview-0.47.1 → arrayview-0.47.4}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.5
 Name: arrayview
-Version: 0.47.1
+Version: 0.47.4
 Summary: Fast multi-dimensional array viewer
 Project-URL: Home, https://github.com/oscarvanderheide/arrayview
 Project-URL: Source, https://github.com/oscarvanderheide/arrayview

{arrayview-0.47.1 → arrayview-0.47.4}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "hatchling.build"
 
 [project]
 name = "arrayview"
-version = "0.47.1"
+version = "0.47.4"
 description = "Fast multi-dimensional array viewer"
 readme = { file = "README.md", content-type = "text/markdown" }
 requires-python = ">=3.12"

{arrayview-0.47.1 → arrayview-0.47.4}/src/arrayview/_io.py

@@ -667,7 +667,7 @@ class _FileSeries:
     """Lazy 4D/5D view over a directory of same-shape array files.
 
     Like ``_NiftiSeries`` but works with any supported file format
-    (.npy, .npz, .zarr, .pt/.pth, .h5/.hdf5, .tif/.tiff, .mat).
+    (.npy, .npz, .zarr, .pt/.pth, .h5/.hdf5, .tif/.tiff, .mat, .gipl).
     Uses ``load_data`` for each volume — no nibabel reorientation.
     """
 
@@ -940,6 +940,8 @@ def _get_ext(filepath):
         return ".nii.gz"
     if lower.endswith(".zarr.zip"):
         return ".zarr.zip"
+    if lower.endswith(".gipl.gz"):
+        return ".gipl.gz"
     return os.path.splitext(lower)[1]
 
 
@@ -950,7 +952,7 @@ def _is_supported_collection_path(filepath):
 def _strip_array_ext(filepath):
     base = os.path.basename(filepath)
     lower = base.lower()
-    for suffix in (".nii.gz", ".zarr.zip"):
+    for suffix in (".nii.gz", ".zarr.zip", ".gipl.gz"):
         if lower.endswith(suffix):
             return base[: -len(suffix)]
     return os.path.splitext(base)[0]
@@ -1738,6 +1740,8 @@ def load_data(filepath, key=None, *, progress=None):
         import tifffile
 
         return tifffile.imread(filepath)
+    elif filepath.lower().endswith((".gipl", ".gipl.gz")):
+        return _load_gipl(filepath)
     elif filepath.endswith(".mat"):
         try:
             import scipy.io
@@ -1783,10 +1787,68 @@ def load_data(filepath, key=None, *, progress=None):
             return data
         raise ValueError(
             "Unsupported format. Supported: .npy, .npz, .nii/.nii.gz, .zarr, "
-            ".pt/.pth, .h5/.hdf5, .tif/.tiff, .mat"
+            ".pt/.pth, .h5/.hdf5, .tif/.tiff, .mat, .gipl/.gipl.gz"
         )
 
 
+# GIPL (Guys Image Processing Lab) image types -> numpy dtypes. Data is
+# big-endian; complex types store interleaved real/imaginary components.
+_GIPL_DTYPES = {
+    1: ">u1",  # binary, one value per byte
+    7: ">i1",
+    8: ">u1",
+    15: ">i2",
+    16: ">u2",
+    31: ">u4",
+    32: ">i4",
+    64: ">f4",
+    65: ">f8",
+    144: ">i2",  # complex short
+    160: ">i4",  # complex int
+    192: ">f4",  # complex float
+    193: ">f8",  # complex double
+}
+_GIPL_COMPLEX = frozenset([144, 160, 192, 193])
+_GIPL_MAGIC = (4026526128, 719555000)
+_GIPL_HEADER_SIZE = 256
+
+
+def _load_gipl(filepath):
+    """Read a .gipl/.gipl.gz file into an array indexed (x, y, z, t).
+
+    Trailing singleton dimensions are dropped, matching how a NIfTI volume
+    of the same data is shown.
+    """
+    import gzip
+
+    opener = gzip.open if filepath.lower().endswith(".gz") else open
+    with opener(filepath, "rb") as f:
+        raw = f.read()
+    if len(raw) < _GIPL_HEADER_SIZE:
+        raise ValueError(f"{filepath!r} is too short to be a GIPL file.")
+    header = raw[:_GIPL_HEADER_SIZE]
+    magic = int(np.frombuffer(header, ">u4", count=1, offset=252)[0])
+    if magic not in _GIPL_MAGIC:
+        raise ValueError(f"{filepath!r} is not a GIPL file (bad magic number).")
+    dims = [int(d) for d in np.frombuffer(header, ">u2", count=4, offset=0)]
+    image_type = int(np.frombuffer(header, ">u2", count=1, offset=8)[0])
+    if image_type not in _GIPL_DTYPES:
+        raise ValueError(f"{filepath!r} has unsupported GIPL image type {image_type}.")
+    dims = [max(d, 1) for d in dims]
+    is_complex = image_type in _GIPL_COMPLEX
+    count = int(np.prod(dims)) * (2 if is_complex else 1)
+    data = np.frombuffer(raw, _GIPL_DTYPES[image_type], count=count, offset=_GIPL_HEADER_SIZE)
+    data = data.astype(data.dtype.newbyteorder("="))
+    if is_complex:
+        data = data.astype(np.float64 if data.itemsize == 8 else np.float32)
+        data = data[0::2] + 1j * data[1::2]
+    # x varies fastest on disk.
+    arr = data.reshape(dims[::-1]).transpose(3, 2, 1, 0)
+    while arr.ndim > 2 and arr.shape[-1] == 1:
+        arr = arr[..., 0]
+    return np.ascontiguousarray(arr)
+
+
 def _tensor_to_numpy(obj, source="object"):
     """Convert a tensor-like object (PyTorch, etc.) to a numpy array."""
     if isinstance(obj, np.ndarray):
@@ -1823,13 +1885,15 @@ _SUPPORTED_EXTS = frozenset(
         ".tif",
         ".tiff",
         ".mat",
+        ".gipl",
+        ".gipl.gz",
         ".dcm",
     ]
 )
 
 # Formats that load the entire array into RAM (no mmap/lazy access).
 # RAM guard only applies to these.
-FULL_LOAD_EXTS = frozenset([".pt", ".pth", ".tif", ".tiff", ".mat"])
+FULL_LOAD_EXTS = frozenset([".pt", ".pth", ".tif", ".tiff", ".mat", ".gipl", ".gipl.gz"])
 
 
 def _peek_file_shape(fpath: str, ext: str):

{arrayview-0.47.1 → arrayview-0.47.4}/src/arrayview/_launcher.py

@@ -6211,17 +6211,18 @@ def arrayview():
             pattern for _name, pattern in _normalize_dir_overlay_specs(args.overlay or [])
         )
         collection_inputs.extend(args.overlay_dir or [])
-        for pattern in collection_inputs:
-            pattern = lexical_abspath(pattern)
-            probe_prefix = scan_root_before_magic(pattern)
-            mount = direct_network_mount(probe_prefix)
-            nested_mount = network_mount_below(probe_prefix)
-            if mount is not None or nested_mount is not None:
-                _print_failure(
-                    "network directory collections are not accessed directly. "
-                    "Copy the collection to local storage before using --match."
-                )
-                sys.exit(1)
+        if os.environ.get("ARRAYVIEW_SKIP_SOURCE_STAGING", "1") == "0":
+            for pattern in collection_inputs:
+                pattern = lexical_abspath(pattern)
+                probe_prefix = scan_root_before_magic(pattern)
+                mount = direct_network_mount(probe_prefix)
+                nested_mount = network_mount_below(probe_prefix)
+                if mount is not None or nested_mount is not None:
+                    _print_failure(
+                        "network directory collections are not accessed directly. "
+                        "Copy the collection to local storage before using --match."
+                    )
+                    sys.exit(1)
         if len(dir_patterns) == 1 and os.path.isdir(dir_patterns[0]):
             dir_patterns = [os.path.join(dir_patterns[0], "**", "*")]
         dir_overlay_specs = [

{arrayview-0.47.1 → arrayview-0.47.4}/src/arrayview/_routes_loading.py

@@ -167,17 +167,17 @@ def register_loading_routes(app, *, notify_shells, setup_rgb) -> None:
                             f"{nested_sibling_mount.mountpoint!r}; copy the series locally first."
                         )
                     }
-        for pattern in dir_patterns or []:
-            pattern_path = lexical_abspath(str(pattern))
-            scan_root = scan_root_before_magic(pattern_path)
-            unsafe_mount = direct_network_mount(scan_root) or network_mount_below(scan_root)
-            if unsafe_mount is not None:
-                return {
-                    "error": (
-                        f"Directory scan would enter network mount "
-                        f"{unsafe_mount.mountpoint!r}; copy the collection locally first."
-                    )
-                }
+            for pattern in dir_patterns or []:
+                pattern_path = lexical_abspath(str(pattern))
+                scan_root = scan_root_before_magic(pattern_path)
+                unsafe_mount = direct_network_mount(scan_root) or network_mount_below(scan_root)
+                if unsafe_mount is not None:
+                    return {
+                        "error": (
+                            f"Directory scan would enter network mount "
+                            f"{unsafe_mount.mountpoint!r}; copy the collection locally first."
+                        )
+                    }
         if requested_sid in CANCELLED_PENDING_SESSIONS:
             CANCELLED_PENDING_SESSIONS.discard(requested_sid)
             raise HTTPException(status_code=409, detail="load request was cancelled")

{arrayview-0.47.1 → arrayview-0.47.4}/src/arrayview/_routes_preferences.py

@@ -37,7 +37,7 @@ def _response_payload() -> dict:
     if not native_available:
         defaults["window"]["default"] = "browser"
         defaults["window"]["terminal"] = "browser"
-        for key in ("default", "terminal", "vscode", "jupyter", "julia"):
+        for key in ("terminal", "vscode", "jupyter", "julia"):
             schema["window"][key] = [
                 value for value in schema["window"][key] if value != "native"
             ]

{arrayview-0.47.1 → arrayview-0.47.4}/src/arrayview/_source_safety.py

@@ -53,6 +53,8 @@ _KNOWN_LOCAL_STAGE_FILESYSTEMS = frozenset(
     }
 )
 _SELF_CONTAINED_FILE_EXTENSIONS = (
+    ".gipl",
+    ".gipl.gz",
     ".h5",
     ".hdf5",
     ".mat",

{arrayview-0.47.1 → arrayview-0.47.4}/src/arrayview/_viewer.html

@@ -7086,6 +7086,7 @@
         let _mvPrevZoom = 1.0;
         let _cmvPrevSquareStretch = false;
         let _qmriPrevZoom = 1.0;  // saved on enterQmri, restored on exitQmri (mirrors _mvPrevZoom)
+        let _qmriPrevActiveDim = null;  // scroll dim before qMRI moved it off the map dim; null while rebuilding keeps the original
         let _comparePrevZoom = 1.0;  // saved on enterCompareMode, restored on exitCompareMode
         let _comparePrevNormalView = null; // Phase 4: saved primary View, restored on exitCompareMode
         let _eggsVisible = true;
@@ -35472,6 +35473,7 @@
             if (multiViewActive) exitMultiView();
             if (stackMosaicActive) exitStackMosaic();
             _qmriPrevZoom = userZoom;  // save prior zoom so exitQmri can restore (mirrors enterMultiView)
+            if (_qmriPrevActiveDim === null) _qmriPrevActiveDim = activeDim;
             userZoom = 1.0;
             qvPan.reset();
             // Opening straight into qMRI never draws the main canvas, so its
@@ -35994,6 +35996,10 @@
             wrap.innerHTML = '';
             wrap.className = '';
             userZoom = _qmriPrevZoom;  // restore prior zoom (mirrors exitMultiView line 10521)
+            if (!_qmriSyntheticRebuilding && _qmriPrevActiveDim !== null) {
+                if (_qmriPrevActiveDim < shape.length) activeDim = _qmriPrevActiveDim;
+                _qmriPrevActiveDim = null;
+            }
             if (!_qmriSyntheticRebuilding) updateView();
             _reconcileUI();
             renderInfo();

{arrayview-0.47.1 → arrayview-0.47.4}/src/arrayview/_vscode_extension.py

@@ -58,7 +58,7 @@ _VSCODE_EXT_FRESH_INSTALL = False  # True if we just installed it this session
 _VSCODE_EXT_RELOAD_REQUIRED = False  # installed files are newer than the live host
 _VSCODE_EXT_INSTALL_FAILED = False  # automatic install could not complete safely
 _VSCODE_EXT_NO_LIVE_WINDOW = False  # no live host claims this terminal's window id
-_VSCODE_EXT_VERSION = "0.15.62"  # current bundled extension version
+_VSCODE_EXT_VERSION = "0.15.64"  # current bundled extension version
 _VSCODE_CONFIGURED_PORTS: set[int] = set()
 # Version-skew notices already printed, so a long-lived process (Jupyter, a
 # script opening several arrays) states the mismatch once instead of per call.

{arrayview-0.47.1 → arrayview-0.47.4}/tests/conftest.py

@@ -34,6 +34,21 @@ def isolated_runtime_directory(tmp_path_factory, monkeypatch):
     )
 
 
+@pytest.fixture(autouse=True)
+def no_leaked_fake_sessions():
+    """Drop anything a test put in the shared session table that isn't a session.
+
+    The test server is shared by the whole run, so one stray placeholder makes
+    every later load in every later file fail.
+    """
+    from arrayview._session import SESSIONS, Session
+
+    yield
+    for sid, value in list(SESSIONS.items()):
+        if not isinstance(value, Session):
+            SESSIONS.pop(sid, None)
+
+
 # ---------------------------------------------------------------------------
 # Server
 # ---------------------------------------------------------------------------