pulserver 0.3.0__tar.gz → 0.4.0__tar.gz

{pulserver-0.3.0 → pulserver-0.4.0}/.github/workflows/wheels.yml

@@ -30,7 +30,7 @@ jobs:
         with:
           fetch-depth: 0 # setuptools_scm reads the tags
 
-      - uses: pypa/cibuildwheel@v4.2.1
+      - uses: pypa/cibuildwheel@v4.3.1
         env:
           CIBW_BUILD: ${{ inputs.cibw-build }}
           CIBW_SKIP: "*-musllinux_*"

{pulserver-0.3.0 → pulserver-0.4.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pulserver
-Version: 0.3.0
+Version: 0.4.0
 Summary: Orchestrator for MR acquisitions: sequence design, scanner preparation and reconstruction.
 Keywords: mri,pulseq,mrd,ismrmrd,reconstruction,scanner
 Author-Email: Matteo Cencini <matteo.cencini@gmail.com>
@@ -21,7 +21,7 @@ Project-URL: Documentation, https://pulserver.github.io/pulserver/
 Project-URL: Issues, https://github.com/pulserver/pulserver/issues
 Project-URL: Source, https://github.com/pulserver/pulserver
 Requires-Python: >=3.10
-Requires-Dist: pypulseqpp>=0.0.40
+Requires-Dist: pypulseqpp>=0.0.42
 Requires-Dist: numpy>=2.2.6
 Requires-Dist: ismrmrd>=1.15.0
 Requires-Dist: pydicom>=3.0.2
@@ -36,7 +36,7 @@ Requires-Dist: pytest>=9.1.1; extra == "dev"
 Requires-Dist: pytest-cov>=7.1.0; extra == "dev"
 Requires-Dist: pytest-xdist>=3; extra == "dev"
 Requires-Dist: coverage>=7.16.2; extra == "dev"
-Requires-Dist: ruff==0.16.10; extra == "dev"
+Requires-Dist: ruff==0.17.0; extra == "dev"
 Requires-Dist: pre-commit>=4.6.2; extra == "dev"
 Requires-Dist: build; extra == "dev"
 Provides-Extra: doc

{pulserver-0.3.0 → pulserver-0.4.0}/docs/api/recon.md

@@ -70,3 +70,4 @@ sensitivity maps, noise prewhitening and coil compression are the subject of
 | {obj}`~pulserver.recon.ReconData` | The readouts of one reconstruction unit: imaging and calibration k-space, counters and waveforms. |
 | {obj}`~pulserver.recon.ReconBuffer` | K-space of one encoding space, filled one acquisition at a time, with the acquisition nearest the k-space centre. |
 | {obj}`~pulserver.recon.ReconResult` | Image array the runtime packages as an MRD image or DICOM dataset, with the acquisition its geometry is taken from and the rescale of its DICOM pixels. |
+| {obj}`~pulserver.recon.figure_rgb` | Return a matplotlib figure drawn to the pixels of an RGB result, sent as a DICOM secondary capture. |

{pulserver-0.3.0 → pulserver-0.4.0}/docs/developer-guide/internals/virtual-scanner.md

@@ -227,9 +227,14 @@ the models show none of the dielectric effects of a wavelength comparable to
 the head.
 
 A pulse plays on every transmit channel, weighted by the coil's default shim,
-the unit weights that bring the channels into phase at the isocentre, where the
-pulse then has its nominal amplitude: each cube sees the transmit field
-$b^+(\mathbf{r}) = \sum_c w_c\,s^+_c(\mathbf{r})$, for the weights $w_c$.
+the weights of magnitude $1/\sqrt{N}$ for $N$ channels that bring the channels
+into phase at the isocentre, where the pulse then has its nominal amplitude:
+each cube sees the transmit field
+$b^+(\mathbf{r}) = \sum_c w_c\,s^+_c(\mathbf{r})$, for the weights $w_c$. The
+squares of the weights sum to 1, so the channels together carry the power of
+one channel at the pulse's amplitude, and an RF shim is written in the same
+normalisation: a circularly polarised shim of $N$ channels has weights of
+magnitude $1/\sqrt{N}$.
 
 A console given field maps takes the coils, and the subject's field, from
 electromagnetic and susceptibility models of BrainWeb's head instead
@@ -506,8 +511,8 @@ The maps' transmit coils come with the VOPs mariepy compresses from the same
 fields, and every design is made under those of the exam's transmit coil: its
 limits name the VOP file, the default shim, the drive of every channel per
 hertz of a pulse's amplitude with which the channels' fields reach that
-amplitude at the isocentre, $2 / (\gamma \sum_c |B^-_c(\mathbf{0})|)$ in the
-maps' unit of drive, and the head and local SAR limits of IEC 60601-2-33's
+amplitude at the isocentre through the default shim,
+$2\sqrt{N} / (\gamma \sum_c |B^-_c(\mathbf{0})|)$ in the maps' unit of drive, and the head and local SAR limits of IEC 60601-2-33's
 normal operating mode. The IR cache then reports each subsequence's SAR against
 the reference pulse in that coil ({doc}`../../explanations/safety-checks`).
 

{pulserver-0.3.0 → pulserver-0.4.0}/docs/user-guide/reconstruction-client.md

@@ -49,13 +49,44 @@ connection, so the client can finish sending and read the reason.
 | `ExamID`, a `userParameterString`, or else `studyInformation.studyInstanceUID`, or else `studyInformation.studyID` | Optional | The exam; the series of one exam share its exam cache, and a series without one is an exam of its own |
 | `fov_offset_mm`, a `userParameterString` | Optional | The position of the field-of-view centre, three millimetres along the logical readout, phase and slice axes separated by spaces; a difference from the offset the design was converted at is applied to the samples as a receive phase. Without it, the converted offset |
 | `subjectInformation`, `studyInformation`, `measurementInformation`, `acquisitionSystemInformation` | Optional | Copied into the DICOM datasets; `relativeTablePosition` is not |
-| `measurementInformation.measurementID` | A non-negative integer when `systemVendor` names GE | The DICOM series number on a GE system |
+| `measurementInformation.measurementID` | A non-negative integer when `systemVendor` names GE | The DICOM series number on a GE system: of the first series the reconstruction returns, and `100 × measurementID + k` of the `k`-th after it |
+| `dicom:<attribute>`, a `userParameterString` | Optional | A DICOM attribute of every image ({ref}`dicom-attributes`) |
 | `encoding` | Optional | Replaced by one encoding space per subsequence, and one more for the navigator readouts of a subsequence that has them. The reconstruction space is the sequence's matrix size and field of view, one slice deep for a stack of 2D slices; the encoded space of a Cartesian space has the readout widened to the full echo, oversampling included, and its field of view in proportion. The encoding limits are the extremes of the counters the readouts reach, centred on the k-space centre line and partition the sequence defines. A matrix size or field of view the sequence does not define is kept from the client's encoding at the same index |
 | `sequenceParameters` | Optional | TR, TE, TI and flip angles replaced by those the sequence defines; TR and TE it does not define are measured by pypulseqpp's `test_report_dict`, flip angles it does not define are the distinct values of `rf_flip_angles` |
 
 A header that the ISMRMRD schema does not accept, or that names no design the
 store holds for its scan, is refused before any acquisition is read.
 
+(dicom-attributes)=
+### DICOM attributes
+
+A client states the DICOM attributes the MRD header has no element for as
+`userParameterString` entries named `dicom:<attribute>`. The DICOM builder
+reads them once, with the rest of the header, into the attributes every image
+of the series starts from; an entry whose value names a quantity of the image
+is filled in for each image instead. An image's own meta attributes take
+precedence over both. `<attribute>` is one of:
+
+| Form | Example | Resolved by |
+| --- | --- | --- |
+| A keyword of the DICOM standard | `dicom:ImagedNucleus` = `23Na` | pydicom's dictionary |
+| An explicit tag and value representation | `dicom:(0011,1011)FL` = `{location}` | Itself; a private group's creator is one more entry, `dicom:(0011,0010)LO` |
+| A keyword of a private dictionary | `dicom:PulseSequenceLabel` = `pulserver` | The file `dicom:dictionary` names, in MATLAB's format, one `(gggg,eeee) VR Keyword VM` line per attribute |
+
+A value separates multiple values with `\` and may name, in braces, a quantity
+of the image it is written into; `{{` and `}}` are literal braces:
+
+| Name | Value |
+| --- | --- |
+| `slice_location` | The image centre along the plane normal, in mm: the `SliceLocation` written by default |
+| `location` | The image centre along the patient axis the normal lies closest to, in mm, positive towards right, anterior or superior; for an oblique plane, along the normal turned towards those |
+| `location_letter` | The side `location` lies on: `R` or `L`, `A` or `P`, `S` or `I`, and `+` or `-` for an oblique plane |
+| `instance_number`, `series_number` | The image's `InstanceNumber` and `SeriesNumber` |
+
+An attribute that resolves to no tag, for want of the dictionary, is logged and
+left out. `pulserver replay --dicom-dictionary` names another dictionary for a
+series replayed where the header's is not.
+
 ## Acquisitions
 
 The proxy matches acquisitions to readouts by their position in the stream, so

{pulserver-0.3.0 → pulserver-0.4.0}/docs/user-guide/reconstruction-plugins.md

@@ -284,7 +284,7 @@ receives carries the message. `required=False` returns `None` instead of
 raising.
 
 (dicom-pixel-values)=
-## DICOM pixel values
+## DICOM output
 
 A DICOM dataset stores integer pixels, made when a result with `dicom=True` is
 converted from the image the plugin returned. The dataset's rescale maps them
@@ -331,6 +331,32 @@ with one mapping per DICOM series, a series being the images that share
 
 `RescaleType` is `US`, unspecified: the units are those of the reconstruction.
 
+A DICOM image is written in the radiological display orientation of the plane
+it lies closest to: an axial image with the patient's right on the left and
+anterior at the top, a coronal one with the right on the left and superior at
+the top, a sagittal one with anterior on the left and superior at the top. Its
+pixels are flipped and transposed to that orientation, never resampled; the MRD
+image keeps the orientation the reconstruction returned.
+
+Each `series_index` is a DICOM series of its own. The first to arrive takes the
+series identity the header states; each later one a series UID derived from it
+and, on a GE system, the series number `100 × measurementID + k` for the
+`k`-th after the first.
+
+A plot, a fit or a report reaches the scanner as an RGB result: a `(y, x, 3)`
+array of values from 0 to 255, which {func}`~pulserver.recon.figure_rgb` draws
+from a matplotlib figure, stored as a DICOM secondary capture.
+
+```python
+return recon.ReconResult(
+    recon.figure_rgb(figure),
+    image_type="rgb",
+    series_index=10,
+    dicom=True,
+    attributes={"SeriesDescription": "T1 fit"},
+)
+```
+
 ## Diffusion tables
 
 The proxy states the diffusion table a design names in the header's
@@ -365,12 +391,17 @@ process ID, level and logger name.
 {meth}`~pulserver.recon.ReconPlugin.run` reconstructs an ISMRMRD HDF5 file in
 the calling process, through the same hooks. A file recorded as the scanner
 sends it, such as `pulserver scan --mrd raw.h5` or the proxy's `--save-data`
-writes, is enriched from its design store first:
+writes, is enriched first from the design it was played from, found in a design
+store or given as its sequence file, with its manifest beside it:
 
 ```python
 images = PLUGIN.run("raw.h5", store="designs")
+images = PLUGIN.run("raw.h5", design="designs/scan/9c74a8aa28527cf231/sequence.seq")
 ```
 
+`pulserver replay` does the same from the command line
+({ref}`series-replay`).
+
 Calling the plugin on an {class}`~pulserver.mrd.AcquisitionBucket` runs
 `startup`, passes every acquisition to `receive`, reconstructs the units still
 open, runs `finish`, and returns the last output. A header need only describe

{pulserver-0.3.0 → pulserver-0.4.0}/docs/user-guide/running.md

@@ -167,7 +167,7 @@ gradient heating from it.
 | `vop_file` | `.mat` or `.npz` file of VOPs and the head SAR matrix of each body model, at a path the design calls can read; its `safety_factor` multiplies local SAR |
 | `vop_head_limit`, `vop_local_limit` | The scanner's head and local SAR limits in W/kg in its current operating mode |
 | `vop_drive_per_hz` | Channel drive per Hz of RF amplitude, one value or one per channel separated by spaces; a scale common to every channel cancels. Equal drives when left out |
-| `vop_default_shim` | Magnitude and phase in rad of each channel's weight for a pulse played without an RF shim, separated by spaces; equal weights when left out |
+| `vop_default_shim` | Magnitude and phase in rad of each channel's weight for a pulse played without an RF shim, separated by spaces; equal weights whose squares sum to 1 when left out |
 | `vop_coil` | Transmit configuration the scanner reports; the VOP file's `transmit` metadata must name the same one |
 | `acoustic_file` | HDF5 file of the gradient coil's acoustic transfer function, at a path the design calls can read; the peak and A-weighted average sound pressure levels are held to the limits of IEC 60601-2-33 and written into the cache |
 | `acoustic_interval_us` | Sampling interval in µs that the bins of `acoustic_file` refer to; given with it |
@@ -305,10 +305,10 @@ place of a new reconstruction.
 
 With `--save-data`, each series is also written to
 `<save-data>/mrd_<timestamp>.h5` as the scanner sends it, header first and each
-acquisition as it arrives, before it is enriched. A plugin run on that file
-against the store reconstructs it offline exactly as the proxy did
-({ref}`offline-reconstruction`). The queue directory holds a series only until
-it is reconstructed.
+acquisition as it arrives, before it is enriched. `pulserver replay` runs a
+plugin on that file and reconstructs it offline as the proxy did
+({ref}`series-replay`). The queue directory holds a series only until it is
+reconstructed.
 
 ### Proxy lifetime
 
@@ -387,6 +387,27 @@ the proxy forwarded it, enriched. Its stream is neither authenticated nor
 encrypted, and `--host` names the interface the proxy reaches it on. It stops
 on `SIGINT` or `SIGTERM` and waits for the series it is running.
 
+(series-replay)=
+## Series replay
+
+```bash
+pulserver replay SERIES --plugin FILE (--store DIR | --design SEQ) [--config TEXT] [--dicom] [--dicom-dictionary FILE] [--output DIR]
+```
+
+`pulserver replay` reconstructs one series kept by `--save-data` in the calling
+process, with no proxy running. The series is enriched from the design it was
+played from, at the field-of-view offset its header states, as the proxy
+enriches it, and runs through the plugin's hooks
+({ref}`offline-reconstruction`). The design is the one the header names in the
+store `--store`, or the one whose sequence file, the first of its chain, is
+`--design`. That file must be one the manifest beside it records, and the
+manifest must name the design the header names. A series exported beside a
+copy of its design's files replays with no store. Images go to `images.h5` in `--output`, and DICOM datasets to the files
+they are named by; with `--dicom` every image is written as DICOM, as the
+proxy's `--dicom` writes it, and `--dicom-dictionary` names the private
+dictionary of the header's DICOM attributes ({ref}`dicom-attributes`). An exception raised by the plugin is not caught, so its
+traceback names the line that raised it.
+
 ## Playout validation
 
 ```bash

{pulserver-0.3.0 → pulserver-0.4.0}/pyproject.toml

@@ -24,7 +24,7 @@ classifiers = [
 ]
 
 dependencies = [
-  "pypulseqpp>=0.0.40",
+  "pypulseqpp>=0.0.42",
   "numpy>=2.2.6",
   "ismrmrd>=1.15.0",
   "pydicom>=3.0.2",
@@ -45,7 +45,7 @@ dev = [
   "coverage>=7.16.2",
   # Pinned, not ranged: a formatter that differs between a contributor and
   # CI makes "one entry point" a lie. Dependabot proposes the bumps.
-  "ruff==0.16.10",
+  "ruff==0.17.0",
   "pre-commit>=4.6.2",
   "build",
 ]

{pulserver-0.3.0 → pulserver-0.4.0}/src/c/include/pulseg/pulseg_types.h

@@ -846,7 +846,7 @@ typedef struct pulseg_block_info
     int rf_num_samples;    /**< samples per channel, -1 if absent */
     int rf_duration_us;    /**< RF shape duration (us), -1 if absent */
     int rf_is_complex;     /**< 1 if phase shape exists           */
-    int rf_uniform_raster; /**< 1 if time shape present           */
+    int rf_uniform_raster; /**< 1 if samples sit on the RF raster  */
 
     /* ADC */
     int has_adc;      /**< 1 if ADC acquisition active       */

{pulserver-0.3.0 → pulserver-0.4.0}/src/c/structure/pulseg_getters.c

@@ -1301,13 +1301,19 @@ static int block_has_rf(const pulseg_collection *coll, int seg_idx, int blk_idx)
     return (bdef->rf_id != -1) ? 1 : 0;
 }
 
+static int get_rf_num_samples(const pulseg_collection *coll, int seg_idx, int blk_idx);
+
+/* 1 when the RF samples sit at the centres of consecutive RF raster periods:
+ * no time shape, or one that places them there, as a writer may state the
+ * default explicitly. */
 static int block_rf_has_uniform_raster(const pulseg_collection *coll, int seg_idx, int blk_idx)
 {
     const pulseg_sequence_descriptor *desc;
     const pulseg_virtual_segment *seg;
-    int local_blk;
+    int local_blk, i, n, uniform;
     const pulseg_base_block *bdef;
     const pulseg_rf_definition *rdef;
+    float *t, r;
 
     if (!resolve_block(coll, &desc, &seg, &local_blk, seg_idx, blk_idx))
         return -1;
@@ -1317,7 +1323,19 @@ static int block_rf_has_uniform_raster(const pulseg_collection *coll, int seg_id
         return -1;
 
     rdef = &desc->rf_definitions[bdef->rf_id];
-    return (rdef->time_shape_id == 0) ? 1 : 0;
+    if (rdef->time_shape_id == 0)
+        return 1;
+
+    t = pulseg_get_rf_time_us(coll, seg_idx, blk_idx);
+    if (!t)
+        return 0;
+    n = get_rf_num_samples(coll, seg_idx, blk_idx);
+    r = desc->rf_raster_us;
+    uniform = (n > 0 && r > 0.0f);
+    for (i = 0; uniform && i < n; ++i)
+        uniform = fabs((double)t[i] - (i + 0.5) * r) <= 1e-3 * r;
+    PULSEG_FREE(t);
+    return uniform;
 }
 
 static float *alloc_uniform_time_us(int num_samples, float raster_us)

{pulserver-0.3.0 → pulserver-0.4.0}/src/pulserver/_cli.py

@@ -11,6 +11,7 @@ _USAGE = (
     "       pulserver validate SEQ [--vendor NAME] [--played FILE] ...\n"
     "       pulserver proxy --store DIR --port N (--plugins DIR | --forward HOST:PORT) ...\n"
     "       pulserver recon --plugins DIR --port N ...\n"
+    "       pulserver replay SERIES --plugin FILE (--store DIR | --design SEQ) ...\n"
 )
 
 
@@ -38,6 +39,10 @@ def main(argv: list[str] | None = None) -> int:
 
         recon(argv[1:])
         return 0
+    if argv and argv[0] == "replay":
+        from .recon._replay import main as replay
+
+        return replay(argv[1:])
     if argv and argv[0] == "validate":
         from .validate._command import main as check
 

{pulserver-0.3.0 → pulserver-0.4.0}/src/pulserver/ir/_checks.py

@@ -59,7 +59,7 @@ class CheckLimits:
         the ratios.
     default_shim
         Complex channel weights of a pulse played without an RF shim, and of
-        the reference pulse; equal weights when None.
+        the reference pulse; equal weights whose squares sum to 1 when None.
     vop_coil
         The transmit configuration the scanner reports. When given, the VOP
         file's ``metadata["transmit"]`` must name the same one.

{pulserver-0.3.0 → pulserver-0.4.0}/src/pulserver/recon/__init__.py

@@ -8,6 +8,7 @@ from ._calibration import (
 )
 from ._correct import gradient_unwarped
 from ._exam_image import ExamImage
+from ._figure import figure_rgb
 from ._loader import load_plugin
 from .gadgets import AsymmetricEcho, Prewhiten, RemoveReadoutOversampling
 from .plugin import (
@@ -46,6 +47,7 @@ __all__ = [
     "ReconResult",
     "RemoveReadoutOversampling",
     "coil_maps",
+    "figure_rgb",
     "gradient_unwarped",
     "load_plugin",
 ]

pulserver-0.4.0/src/pulserver/recon/_figure.py

@@ -0,0 +1,31 @@
+"""A matplotlib figure as the pixels of an RGB reconstruction result."""
+
+from __future__ import annotations
+
+__all__ = ["figure_rgb"]
+
+import io
+from typing import Any
+
+import numpy as np
+
+
+def figure_rgb(figure: Any, dpi: float | None = None) -> np.ndarray:
+    """Return a matplotlib figure drawn to a ``(y, x, 3)`` ``uint8`` array.
+
+    Pass it as ``ReconResult(figure_rgb(fig), image_type="rgb")`` to send a
+    plot, a fit or a report to the scanner as a secondary capture.
+
+    Parameters
+    ----------
+    figure : matplotlib.figure.Figure
+        The figure; drawn as ``savefig`` draws it, whatever its canvas.
+    dpi : float, default=None
+        Resolution to draw at; the figure's own when not given.
+    """
+    dpi = figure.dpi if dpi is None else dpi
+    buffer = io.BytesIO()
+    figure.savefig(buffer, format="rgba", dpi=dpi)
+    width = round(figure.get_figwidth() * dpi)
+    rgba = np.frombuffer(buffer.getvalue(), dtype=np.uint8).reshape(-1, width, 4)
+    return rgba[..., :3].copy()

pulserver-0.4.0/src/pulserver/recon/_replay.py

@@ -0,0 +1,75 @@
+"""Reconstruct a kept series in this process: ``pulserver replay SERIES --plugin FILE``."""
+
+from __future__ import annotations
+
+import argparse
+from pathlib import Path
+
+
+def _parser() -> argparse.ArgumentParser:
+    parser = argparse.ArgumentParser(
+        prog="pulserver replay",
+        description="Reconstruct a series as the scanner sent it, in this process, "
+        "enriched from the design it was played from as the proxy enriches it.",
+    )
+    parser.add_argument(
+        "series",
+        type=Path,
+        help="ISMRMRD file of the series, as the proxy's --save-data keeps it",
+    )
+    parser.add_argument(
+        "--plugin", type=Path, required=True, help="reconstruction plugin file"
+    )
+    played = parser.add_mutually_exclusive_group(required=True)
+    played.add_argument(
+        "--store", type=Path, help="design store holding the design the header names"
+    )
+    played.add_argument(
+        "--design",
+        type=Path,
+        help="the design's sequence file, the first of its chain, with its "
+        "manifest beside it",
+    )
+    parser.add_argument(
+        "--config", help="configuration payload, as the scanner's stream carries it"
+    )
+    parser.add_argument(
+        "--dicom",
+        action="store_true",
+        help="write every image as DICOM, as the proxy's --dicom does",
+    )
+    parser.add_argument(
+        "--dicom-dictionary",
+        type=Path,
+        help="private DICOM dictionary for the header's dicom: attributes, in "
+        "place of the one the header names",
+    )
+    parser.add_argument(
+        "--output",
+        type=Path,
+        default=Path(),
+        help="directory the reconstruction's images are written to",
+    )
+    return parser
+
+
+def main(argv: list[str] | None = None) -> int:
+    """Run ``pulserver replay`` with ``argv``; return its exit status.
+
+    Images go to ``images.h5`` in ``--output`` and DICOM datasets to the files
+    they are named by. A plugin's exception is not caught, so its traceback
+    names the line that raised it.
+    """
+    from ..virtual._command import _received
+    from ._loader import load_plugin
+
+    args = _parser().parse_args(argv)
+    emitted = load_plugin(args.plugin).run(
+        str(args.series),
+        config=args.config,
+        store=None if args.store is None else str(args.store),
+        design=None if args.design is None else str(args.design),
+        dicom=args.dicom,
+        dicom_dictionary=args.dicom_dictionary,
+    )
+    return _received(emitted, args.output)

{pulserver-0.3.0 → pulserver-0.4.0}/src/pulserver/recon/_runtime/application.py

@@ -157,6 +157,13 @@ def _make_images(
 ) -> tuple[list[ismrmrd.Image], int]:
     data = as_numpy(result.data)
     unit = _unit_header(result, bucket, context.header)
+    if result.image_type.lower() == "rgb":
+        if data.ndim != 3 or data.shape[-1] != 3:
+            raise ValueError(f"an RGB result is (y, x, 3), not {data.shape}")
+        planes = np.clip(np.moveaxis(data, -1, 0), 0, 255).astype(np.uint16)
+        rgb = dataclasses.replace(result, data=planes[:, None])
+        image, next_image_index = _make_image(rgb, unit, next_image_index, app_name)
+        return [image], next_image_index
     if data.ndim != 3 or data.shape[0] == 1:
         image, next_image_index = _make_image(result, unit, next_image_index, app_name)
         return [image], next_image_index
@@ -297,6 +304,7 @@ def _image_type(name: str) -> int:
         "real": ismrmrd.IMTYPE_REAL,
         "imaginary": ismrmrd.IMTYPE_IMAG,
         "complex": ismrmrd.IMTYPE_COMPLEX,
+        "rgb": ismrmrd.IMTYPE_RGB,
     }
     try:
         return types[name.lower()]

{pulserver-0.3.0 → pulserver-0.4.0}/src/pulserver/recon/_runtime/exam.py

@@ -277,7 +277,7 @@ def exam_directory_name(exam_id: Hashable) -> str:
     """Return the name of an exam's directory under the exam root.
 
     An identity :func:`resolve_exam_id` returns is named ``<source>-<value>``,
-    ``exam-5084`` for the ``ExamID`` 5084, so a design call can name the
+    ``exam-42`` for the ``ExamID`` 42, so a design call can name the
     directory of the exam it is made in from the exam's number. Any other
     identity, and a value that is not a plain file name, is named by the
     SHA-256 of its ``repr``.