PurposeVisualization of medical images on a virtual reality (VR) head-mounted display (HMD) requires binocular fusion of a stereoscopic pair of graphical views. However, current image quality assessment on VR HMDs for medical applications has been primarily limited to time-consuming monocular optical bench measurement on a single eyepiece.ApproachAs an alternative to optical bench measurement to quantify the image quality on VR HMDs, we developed a WebXR test platform to perform contrast perceptual experiments that can be used for binocular image quality assessment. We obtained monocular and binocular contrast sensitivity responses (CSRs) from participants on a Meta Quest 2 VR HMD using varied interpupillary distance (IPD) configurations.ResultsThe perceptual result shows that contrast perception on VR HMDs is primarily affected by optical aberration of the VR HMD. As a result, monocular CSR degrades at a high spatial frequency greater than 4 cycles per degree when gazing at the periphery of the display field of view, especially for mismatched IPD settings consistent with optical bench measurements. On the contrary, binocular contrast perception is dominated by the monocular view with superior image quality measured by the contrast.ConclusionsWe developed a test platform to investigate monocular and binocular contrast perception by performing perceptual experiments. The test method can be used to evaluate monocular and/or binocular image quality on VR HMDs for potential medical applications without extensive optical bench measurements.
Purpose: Unlike conventional displays, medical image perception in medical extended reality (MXR) applications often involves relative motion between the digital content and the subject with additional sources of noise in the spatial and temporal domains that affects the MXR image quality.
Methods: We describe a spatiotemporal image perception model with static and dynamic signal and noise configurations. The 3D spatiotemporal noise is decomposed into 2D spatial noise and time-dependent noise with motion. The noise in the temporal domain is categorized into time-invariant fixed pattern noise (FPN) and temporal noise that varies per display frame. Visual integration of the moving signal and noise emulates the spatiotemporal image perception of dynamic detection targets in a smooth-pursuit event. A target detection model is implemented to compute the detectability of both low-contrast and high-resolution signal-known-exactly/background-known-exactly (SKE/BKE) targets in various static and dynamic imaging configurations using a non-pre-whitening model observer with eye filter (NPWE).
Results: Smooth pursuit of a moving target suppresses the high-frequency dynamic resolution and noise in the orientation tangential to the motion trajectory. For the dynamic signal and noise configuration, the reduction of both resolution and high-frequency noise results in similar target detectability compared to the reference static image perception. On the other hand, the visibility of a moving target with static FPN is enhanced due to noise aliasing. Visual integration for approximately 33 ms of time-variant temporal noise at 90 Hz display refresh rate reduces the effective noise compared to the FPN by temporal fusion of noise in neighboring display frames.
Conclusion: Spatiotemporal integration of dynamic signal and noise can potentially affect image quality. Complete assessment of image quality in MXR devices needs to consider the contributions from 3D spatiotemporal characteristics.
Purpose: We investigated the feasibility of detection and quantification of bone marrow edema (BME) using dual-energy (DE) Cone-Beam CT (CBCT) with a dual-layer flat panel detector (FPD) and three-material decomposition. Methods: A realistic CBCT system simulator was applied to study the impact of detector quantization, scatter, and spectral calibration errors on the accuracy of fat-water-bone decompositions of dual-layer projections. The CBCT system featured 975 mm source-axis distance, 1,362 mm source-detector distance and a 430 × 430 mm2 dual-layer FPD (top layer: 0.20 mm CsI:Tl, bottom layer: 0.55 mm CsI:Tl; a 1 mm Cu filter between the layers to improve spectral separation). Tube settings were 120 kV (+2 mm Al, +0.2 mm Cu) and 10 mAs per exposure. The digital phantom consisted of a 160 mm water cylinder with inserts containing mixtures of water (volume fraction ranging 0.18 to 0.46) - fat (0.5 to 0.7) - Ca (0.04 to 0.12); decreasing fractions of fat indicated increasing degrees of BME. A two-stage three-material DE decomposition was applied to DE CBCT projections: first, projection-domain decomposition (PDD) into fat-aluminum basis, followed by CBCT reconstruction of intermediate base images, followed by image-domain change of basis into fat, water and bone. Sensitivity to scatter was evaluated by i) adjusting source collimation (12 to 400 mm width) and ii) subtracting various fractions of the true scatter from the projections at 400 mm collimation. The impact of spectral calibration was studied by shifting the effective beam energy (± 2 keV) when creating the PDD lookup table. We further simulated a realistic BME imaging framework, where the scatter was estimated using a fast Monte Carlo (MC) simulation from a preliminary decomposition of the object; the object was a realistic wrist phantom with an 0.85 mL BME stimulus in the radius. Results: The decomposition is sensitive to scatter: approx. <20 mm collimation width or <10% error of scatter correction in a full field-of-view setting is needed to resolve BME. A mismatch in PDD decomposition calibration of ± 1 keV results in ~25% error in fat fraction estimates. In the wrist phantom study with MC scatter corrections, we were able to achieve ~0.79 mL true positive and ~0.06 mL false positive BME detection (compared to 0.85 mL true BME volume). Conclusions: Detection of BME using DE CBCT with dual-layer FPD is feasible, but requires scatter mitigation, accurate scatter estimation, and robust spectral calibration.
Purpose: We compare the effects of scatter on the accuracy of areal bone mineral density (BMD) measurements obtained using two flat-panel detector (FPD) dual-energy (DE) imaging configurations: a dual-kV acquisition and a dual-layer detector. Methods: Simulations of DE projection imaging were performed with realistic models of x-ray spectra, scatter, and detector response for dual-kV and dual-layer configurations. A digital body phantom with 4 cm Ca inserts in place of vertebrae (concentrations 50 - 400 mg/mL) was used. The dual-kV configuration involved an 80 kV low-energy (LE) and a 120 kV high-energy (HE) beam and a single-layer, 43x43 cm FPD with a 650 μm cesium iodide (CsI) scintillator. The dual-layer configuration involved a 120 kV beam and an FPD consisting of a 200 μm CsI layer (LE data), followed by a 1 mm Cu filter, and a 550 μm CsI layer (HE data). We investigated the effects of an anti-scatter grid (13:1 ratio) and scatter correction. For the correction, the sensitivity to scatter estimation error (varied ±10% of true scatter distribution) was evaluated. Areal BMD was estimated from projection-domain DE decomposition. Results: In the gridless dual-kV setup, the scatter-to-primary ratio (SPR) was similar for the LE and HE projections, whereas in the gridless dual layer setup, the SPR was ~26% higher in the LE channel (top CsI layer) than in the HE channel (bottom layer). Because of the resulting bias in LE measurements, the conventional projection-domain DE decomposition could not be directly applied to dual-layer data; this challenge persisted even in the presence of a grid. In contrast, DE decomposition of dual-kV data was possible both without and with the grid; the BMD error of the 400 mg/mL insert was -0.4 g/cm2 without the grid and +0.3 g/cm2 with the grid. The dual-layer FPD configuration required accurate scatter correction for DE decomposition: a -5% scatter estimation error resulted in -0.1 g/cm2 BMD error for the 50 mg/mL insert and a -0.5 g/cm2 BMD error for the 400 mg/mL with a grid, compared to <0.1 g/cm2 for all inserts in a dual-kV setup with the same scatter estimation error. Conclusion: This comparative study of quantitative performance of dual-layer and dual-kV FPD-based DE imaging indicates the need for accurate scatter correction in the dual-layer setup due to increased susceptibility to scatter errors in the LE channel.
We investigate an image-based strategy to compensate for cardiac motion-induced artifacts in Digital Chest Tomosynthesis (DCT). We apply the compensation to conventional unidirectional vertical “↕” scan DCT and to a multidirectional circular trajectory "O" providing improved depth resolution. Propagation of heart motion into the lungs was simulated as a dynamic deformation. The studies investigated a range of motion propagation distances and scan times. Projection-domain retrospective gating was used to detect heart phases. Sparsely sampled reconstructions of each phase were deformably aligned to yield a motion compensated image with reduced sampling artifacts. The proposed motion compensation mitigates artifacts and blurring in DCT images both for “↕” and "O" scan trajectories. Overall, the “O” orbit achieved the same or better nodule structural similarity index in than the conventional “↕” orbit. Increasing the scan time improved the sampling of individual phase reconstructions.
Purpose: We investigate the feasibility of slot-scan dual-energy x-ray absorptiometry (DXA) on robotic x-ray platforms capable of synchronized source and detector translation. This novel approach will enhance the capabilities of such platforms to include quantitative assessment of bone quality using areal bone mineral density (aBMD), normally obtained only with a dedicated DXA scanner. Methods: We performed simulation studies of a robotized x-ray platform that enables fast linear translation of the x-ray source and flat-panel detector (FPD) to execute slot-scan dual-energy (DE) imaging of the entire spine. Two consecutive translations are performed to acquire the low-energy (LE, 80 kVp) and high-energy (HE, 120 kVp) data in <15 sec total time. The slot views are corrected with convolution-based scatter estimation and backprojected to yield tiled long-length LE and HE radiographs. Projection-based DE decomposition is applied to the tiled radiographs to yield (i) aBMD measurements in bone, and (ii) adipose content measurement in bone-free regions. The feasibility of achieving accurate aBMD estimates was assessed using a high-fidelity simulation framework with a digital body phantom and a realistic bone model covering a clinically relevant range of mineral densities. Experiments examined the effects of slot size (1 – 20 cm), scatter correction, and patient size/adipose content (waist circumference: 77 – 95 cm) on the accuracy and reproducibility of aBMD. Results: The proposed combination of backprojection-based tiling of the slot views and DE decomposition yielded bone density maps of the spine that were free of any apparent distortions. The x-ray scatter increased with slot width, leading to aBMD errors ranging from 0.2 g/cm2 for a 5 cm slot to 0.7 g/cm2 for a 20 cm slot when no scatter correction was applied. The convolution-based correction reduced the aBMD error to within 0.02 g/cm2 for slot widths <10 cm. Reproducible aBMD measurements across a range of body sizes (aBMD variability <0.1 g/cm2) were achieved by applying a calibration based on DE adipose thickness estimates from peripheral body sites. Conclusion: The feasibility of accurate and reproducible aBMD measurements on an FPD-based x-ray platform was demonstrated using DE slot scan trajectories, backprojection-domain decomposition, scatter correction, and adipose precorrection.
Purpose: We investigate cone-beam CT (CBCT) imaging protocols and scan orbits for 3D cervical spine imaging on a twin-robotic x-ray imaging system (Multitom Rax). Tilted circular scan orbits are studied to assess potential benefits in visualization of lower cervical vertebrae, in particular in low-dose imaging scenarios. Methods: The Multitom Rax system enables flexible scan orbit design by using two robotic arms to independently move the x-ray source and detector. We investigated horizontal and tilted circular scan orbits (up to 45° tilt) for 3D imaging of the cervical spine. The studies were performed using an advanced CBCT simulation framework involving GPU accelerated x-ray scatter estimation and accurate modeling of x-ray source, detector and noise. For each orbit, the x-ray scatter and scatter-to-primary ratio (SPR) were evaluated; cervical spine image quality was characterized by analyzing the contrast-to-noise ratio (CNR) for each vertebrae. Performance evaluation was performed for a range of scan exposures (263 mAs/scan – 2.63 mAs/scan) and standard and dedicated low dose reconstruction protocols. Results: The tilted orbit reduces scatter and increases primary detector signal for lower cervical vertebrae because it avoids ray paths crossing through both shoulders. Orbit tilt angle of 35° was found to achieve a balanced performance in visualization of upper and lower cervical spine. Compared with a flat orbit, using the optimized 35° tilted orbit reduces lateral projection SPR at the C7 vertebra by <40%, and increases CNR by 220% for C6 and 76% for C7. Adequate visualization of the vertebrae with CNR <1 was achieved for scan exposures as low as 13.2 mAs / scan, corresponding to ~3 mGy absorbed spine dose. Conclusion: Optimized tilted scan orbits are advantageous for CBCT imaging of the cervical spine. The simulation studies presented here indicate that CBCT image quality sufficient for evaluation of spine alignment and intervertebral joint spaces might be achievable at spine doses below 5 mGy.
Purpose: We optimize scan orbits and acquisition protocols for 3D imaging of the weight-bearing spine on a twin-robotic x-ray system (Multitom Rax). An advanced Cone-Beam CT (CBCT) simulation framework is used for systematic optimization and evaluation of protocols in terms of scatter, noise, imaging dose, and task-based performance in 3D image reconstructions. Methods: The x-ray system uses two robotic arms to move an x-ray source and a 43×43 cm2 flat-panel detector around an upright patient. We investigate two classes of candidate scan orbits, both with the same source-axis distance of 690 mm: circular scans with variable axis-detector distance (ADD, air gap) ranging from 400 to 800 mm, and elliptical scans, where the ADD smoothly changes between the anterior-posterior view (ADDAP) and the lateral view (ADDLAT). The study involved elliptical orbits with fixed ADDAP of 400 mm and variable ADDLAT, ranging 400 to 800 mm. Scans of human lumbar spine were simulated using a framework that included accelerated Monte Carlo scatter estimation and realistic models of the x-ray source and detector. In the current work, x-ray fluence was held constant across all imaging configurations, corresponding to 0.5 mAs/frame. Performance of circular and elliptical orbits was compared in terms of scatter and scatter-to-primary ratio (SPR) in projections, and contrast, noise, contrast-to-noise ratio (CNR), and truncation (field of view, FOV) in 3D image reconstructions. Results: The highest mean SPR was found in lateral views, ranging from ~5 at ADD of 300 mm to ~1.2 at ADD of 800 mm. Elliptical scans enabled image acquisition with reduced lateral SPR and almost constant SPR across projection angles. The improvement in contrast across the investigated range of air gaps (due to reduction in scatter) was ~2.3x for circular orbits and ~1.9x for elliptical orbits. The increase in noise associated with increased ADD was more pronounced for circular scans (~2x) compared to elliptical scans (~1.5x). The circular orbit with the best CNR performance (ADD=600 mm) yielded ~10% better CNR than the best elliptical orbit (ADDLAT=600 mm); however, the elliptical orbit increased FOV by ~16%. Conclusion: The flexible imaging geometry of the robotic x-ray system enables development of highly optimized scan orbits. Imaging of the weight-bearing spine could benefit from elliptical detector trajectories to achieve improved tradeoffs in scatter reduction, noise, and truncation.
Access to the requested content is limited to institutions that have purchased or subscribe to SPIE eBooks.
You are receiving this notice because your organization may not have SPIE eBooks access.*
*Shibboleth/Open Athens users─please
sign in
to access your institution's subscriptions.
To obtain this item, you may purchase the complete book in print or electronic format on
SPIE.org.
INSTITUTIONAL Select your institution to access the SPIE Digital Library.
PERSONAL Sign in with your SPIE account to access your personal subscriptions or to use specific features such as save to my library, sign up for alerts, save searches, etc.