In-plane Spatial Resolution Research Articles

Laminar CBV and BOLD response-characteristics over time and space in the human primary somatosensory cortex at 7T.

Uncovering the cortical representation of the body has been at the core of human brain mapping for decades, with special attention given to the digits. In the last decade, advances in functional magnetic resonance imaging (fMRI) technologies have opened the possibility of noninvasively unraveling the 3rd dimension of digit representations in humans along cortical layers. In laminar fMRI it is common to combine the use of the highly sensitive blood oxygen level dependent (BOLD) contrast with cerebral blood volume sensitive measurements, like vascular space occupancy (VASO), that are more specific to the underlying neuronal populations. However, the spatial and temporal VASO response characteristics across cortical depth to passive stimulation of the digits are still unknown. Therefore, we characterized haemodynamic responses to vibrotactile stimulation of individual digit-tips across cortical depth at 0.75 mm in-plane spatial resolution using BOLD and VASO fMRI at 7T. We could identify digit-specific regions of interest (ROIs) in putative Brodmann area 3b, following the known anatomical organization. In the ROIs, the BOLD response increased towards the cortical surface due to the draining vein effect, while the VASO response was more shifted towards middle cortical layers, likely reflecting bottom-up input from the thalamus, as expected. Interestingly, we also found slightly negative BOLD and VASO responses for non-preferred digits in the ROIs, potentially indicating neuronal surround inhibition. Finally, we explored the temporal signal dynamics for BOLD and VASO as a function of distance from activation peaks resulting from stimulation of contralateral digits. With this analysis, we showed a triphasic response consisting of an initial peak and a subsequent negative deflection during stimulation, followed by a positive post-stimulus response in BOLD and to some extent in VASO. While similar responses were reported with invasive methods in animal models, here we demonstrate a potential neuronal excitation-inhibition mechanism in a center-surround architecture across layers in the human somatosensory cortex. Given that, unlike in animals, human experiments do not rely on anesthesia and can readily implement extensive behavioral testing, obtaining this effect in humans is an important step towards further uncovering the functional significance of the different aspects of the triphasic response.

Image quality evaluation in deep-learning-based CT noise reduction using virtual imaging trial methods: Contrast-dependent spatial resolution.

Deep-learning-based image reconstruction and noise reduction methods (DLIR) have been increasingly deployed in clinical CT. Accurate image quality assessment of these methods is challenging as the performance measured using physical phantoms may not represent the true performance of DLIR in patients since DLIR is trained mostly on patient images. In this work, we aim to develop a patient-data-based virtual imaging trial framework and, as a first application, use it to measure the spatial resolution properties of a DLIR method. The patient-data-based virtual imaging trial framework consists of five steps: (1) insertion of lesions into projection domain data using the acquisition geometry of the patient exam to simulate different lesion characteristics; (2) insertion of noise into projection domain data using a realistic photon statistical model of the CT system to simulate different dose levels; (3) creation of DLIR-processed images from projection or image data; (4) creation of ensembles of DLIR-processed patient images from a large number of noise and lesion realizations; and (5) evaluation of image quality using ensemble DLIR images. This framework was applied to measure the spatial resolution of a ResNet based deep convolutional neural network (DCNN) trained on patient images. Lesions in a cylindrical shape and different contrast levels (-500, -100, -50, -20, -10 HU) were inserted to the lower right lobe of the liver in a patient case. Multiple dose levels were simulated (50%, 25%, 12.5%). Each lesion and dose condition had 600 noise realizations. Multiple reconstruction and denoising methods were used on all the noise realizations, including the original filtered-backprojection (FBP), iterative reconstruction (IR), and the DCNN method with three different strength setting (DCNN-weak, DCNN-medium, and DCNN-strong). Mean lesion signal was calculated by performing ensemble averaging of all the noise realizations for each lesion and dose condition and then subtracting the lesion-present images from the lesion absent images. Modulation transfer functions (MTFs) both in-plane and along the z-axis were calculated based on the mean lesion signals. The standard deviations of MTFs at each condition were estimated with bootstrapping: randomly sampling (with replacement) all the DLIR/FBP/IR images from the ensemble data (600 samples) at each condition. The impact of varying lesion contrast, dose levels, and denoising strengths were evaluated. Statistical analysis with paired t-test was used to compare the z-axis and in-plane spatial resolution of five algorithms for five different contrasts and three dose levels. The in-plane and z-axis spatial resolution degradation of DCNN becomes more severe as the contrast or radiation dose decreased, or DCNN denoising strength increased. In comparison with FBP, a 59.5% and 4.1% reduction of in-plane and z-axis MTF (in terms of spatial frequencies at 50% MTF), respectively, was observed at low contrast (-10 HU) for DCNN with the highest denoising strength at 25% routine dose level. When the dose level reduces from 50% to 12.5% of routine dose, the in-plane and z-axis MTFs reduces from 92.1% to 76.3%, and from 98.9% to 95.5%, respectively, at contrast of -100 HU, using FBP as the reference. For most conditions of contrasts and dose levels, significant differences were found among the five algorithms, with the following relationship in both in-plane and cross-plane spatial resolution: FBP>DCNN-Weak>IR>DCNN-Medium>DCNN-Strong. The spatial resolution difference among algorithms decreases at higher contrast or dose levels. A patient-data-based virtual imaging trial framework was developed and applied to measuring the spatial resolution properties of a DCNN noise reduction method at different contrast and dose levels using real patient data. As with other non-linear image reconstruction and post-processing techniques, the evaluated DCNN method degraded the in-plane and z-axis spatial resolution at lower contrast levels, lower radiation dose, and higher denoising strength.

Image quality evaluation for a clinical organ-targeted PET camera.

A newly developed clinical organ-targeted Positron Emission Tomography (PET) system (also known as Radialis PET) is tested with a set of standardized and custom tests previously used to evaluate the performance of Positron Emission Mammography (PEM) systems. Imaging characteristics impacting standardized uptake value (SUV) and detectability of small lesions, namely spatial resolution, linearity, uniformity, and recovery coefficients, are evaluated. In-plane spatial resolution was measured as 2.3 mm ± 0.1 mm, spatial accuracy was 0.1 mm, and uniformity measured with flood field and NEMA NU-4 phantom was 11.7% and 8.3% respectively. Selected clinical images are provided as reference to the imaging capabilities under different clinical conditions such as reduced activity of 2-[fluorine-18]-fluoro-2-deoxy-D-glucose (18F-FDG) and time-delayed acquisitions. SUV measurements were performed for selected clinical acquisitions to demonstrate a capability for quantitative image assessment of different types of cancer including for invasive lobular carcinoma with comparatively low metabolic activity. Quantitative imaging performance assessment with phantoms demonstrates improved contrast recovery and spill-over ratio for this PET technology when compared to other commercial organ-dedicated PET systems with similar spatial resolution. Recovery coefficients were measured to be 0.21 for the 1 mm hot rod and up to 0.89 for the 5 mm hot rod of NEMA NU-4 Image Quality phantom. Demonstrated ability to accurately reconstruct activity in tumors as small as 5 mm suggests that the Radialis PET technology may be well suited for emerging clinical applications such as image guided assessment of response to neoadjuvant systemic treatment (NST) in lesions smaller than 2 cm. Also, our results suggest that, while spatial resolution greatly influences the partial volume effect which degrades contrast recovery, optimized count rate performance and image reconstruction workflow may improve recovery coefficients for systems with comparable spatial resolution. We emphasize that recovery coefficient should be considered as a primary performance metric when a PET system is used for accurate lesion size or radiotracer uptake assessments.

Pre-excitation gradients for eddy current nulled convex optimized diffusion encoding (Pre-ENCODE).

To evaluate the use of pre-excitation gradients for eddy current-nulled convex optimized diffusion encoding (Pre-ENCODE) to mitigate eddy current-induced image distortions in diffusion-weighted MRI (DWI). DWI sequences using monopolar (MONO), ENCODE, and Pre-ENCODE were evaluated in terms of the minimum achievable echo time (TE ) and eddy current-induced image distortions using simulations, phantom experiments, and in vivo DWI in volunteers ( ). Pre-ENCODE provided a shorter TE than MONO (71.0 17.7ms vs. 77.6 22.9ms) and ENCODE (71.0 17.7ms vs. 86.2 14.2ms) in 100 of the simulated cases for a commercial 3T MRI system with b-values ranging from 500 to 3000 s/mm and in-plane spatial resolutions ranging from 1.0 to 3.0mm . Image distortion was estimated by intravoxel signal variance between diffusion encoding directions near the phantom edges and was significantly lower with Pre-ENCODE than with MONO (10.1 vs. 22.7 , ) and comparable to ENCODE (10.1 vs. 10.4 , ). In vivo measurements of apparent diffusion coefficients were similar in global brain pixels (0.37 [0.28,1.45] mm /s vs. 0.38 [0.28,1.45] mm /s, ) and increased in edge brain pixels (0.80 [0.17,1.49] mm /s vs. 0.70 [0.18,1.48] mm /s, ) for MONO compared to Pre-ENCODE. Pre-ENCODE mitigated eddy current-induced image distortions for diffusion imaging with a shorter TE than MONO and ENCODE.

Novel bulk triaxial residual stress mapping in an additive manufactured bridge sample by coupling energy dispersive X-ray diffraction and contour method measurements

A novel approach for determining triaxial residual stress states by coupling energy dispersive X-ray diffraction and contour method measurements is provided and validated in a Ti‐5Al‐5V‐5Mo‐3Cr additive manufactured (AM) bridge sample. Synchrotron X-ray diffraction (SXRD) can provide relatively fine spatial resolution (on the order of 10–100 µm) for mapping 3D elastic strain fields within a sample. However, for samples with dimensions larger than one or two centimeters the path length can get prohibitive as both constant wavelength and energy dispersive SXRD are based upon transmission measurements. As an example, for a plate-like sample geometry where the thickness is limited (less than a centimeter) it is trivial to measure the longitudinal and height direction elastic strain components with excellent in-plane spatial resolution (∼100 µm) and a somewhat lower through-thickness resolution (10–15 times larger) but obtaining the through thickness component is often not possible as the beam path must be parallel to one of the large dimensions of the plate. While small samples of low Z-number alloys (e.g., Ti or Al) allow determination of the three elastic strain components, this is often not the case for higher Z alloys (e.g., Fe or Ni) or large samples where some strain components can be indeterminable. To overcome these limitations, this work applies a combination of synchrotron X-ray diffraction and the contour method (a mechanical relaxation technique), to determine the triaxial stress state. This novel combination is demonstrated and validated in a relatively small additive manufactured sample where all three orthogonal strain components are accessible via X-ray diffraction for stress determination. The paper also explores methods for determining the strain-free lattice parameter, typically obtained from a small stress-free reference sample. This work shows that a small size AM sample is not stress free, producing unreliable magnitudes of strain and stress. Instead, a strain-free lattice parameter is determined using residual stress equilibrium conditions, which gives consistent strain and stress trends for both X-ray diffraction (alone) and the new diffraction-contour coupling technique. This demonstrates that the coupling technique can be confidently applied to samples to determine stress when path length (size or Z-number) prohibit determination of three orthogonal strain components via diffraction.

Simultaneous two-plane flame front detection using PIV based on defocusing.

This Letter presents a simultaneous two-plane flame front detection method. It is based on a standard single-camera single-plane particle image velocimetry (PIV) system in combination with an inexpensive and compact image splitting device. The image splitting optics places images from two depth-offset planes onto the two halves of a camera sensor. A shallow depth of field ensures only one plane is in focus on each half of the sensor. By using a high-pass filter and a novel two-step filter we have devised, the out-of-focus particle images are effectively removed, while the in-focus particle images remain, allowing the turbulent flame fronts on two planes to be detected simultaneously. Our approach could be combined with conventional polarization/wavelength discrimination methods to achieve simultaneous multi-plane flame front reconstruction with similarly high in-plane spatial resolution to single-plane measurement and is suitable for practical combustion devices with limited optical access.

Investigating the coupling between transport and reaction within a catalyst pellet using operando magnetic resonance spectroscopic imaging

In heterogeneous catalysis, the distribution of reactant and product species within a catalyst pellet are of central importance to the catalyst performance, and ultimately to the overall conversion and selectivity achieved in commercial-scale reactors. Despite this importance, fundamental understanding of pellet scale behaviour at operando conditions is limited due to the coupled, multi-scale nature of the underlying transport processes, and the limited techniques capable of mapping the intra-pellet composition at industrially relevant length scales. Here, we demonstrate the coupling between transport and reaction at the pellet scale by using operando magnetic resonance spectroscopic imaging (MRSI) to map the liquid composition within Pd/γ-Al2O3 catalyst pellets at an isotropic in-plane spatial resolution of 86 µm during a styrene hydrogenation reaction. The results reveal that asymmetric liquid wetting at the pellet surface, typical of that occurring in trickle-bed reactors, causes significant composition heterogeneity across a catalyst pellet. Due to the wetting inhomogeneity, commonly used 1D reaction-diffusion models are inadequate for describing the observed composition heterogeneity across the pellet. The results reported here provide unique experimental insight into the coupling between transport and reaction at the pellet scale within an operating reactor.

Native myocardial T1 mapping: influence of spatial resolution on quantitative results and reproducibility.

Myocardial mapping techniques can be used to quantitatively assess alterations in myocardial tissue properties. This study aims to evaluate the influence of spatial resolution on quantitative results and reproducibility of native myocardial T1 mapping in cardiac magnetic resonance imaging (MRI). In this cross-sectional study with prospective data collection between October 2019 and February 2020, 50 healthy adults underwent two identical cardiac MRI examinations in the radiology department on the same day. T1 mapping was performed using a MOLLI 5(3)3 sequence with higher (1.4 mm × 1.4 mm) and lower (1.9 mm × 1.9 mm) in-plane spatial resolution. Global quantitative results of T1 mapping were compared between high-resolution and low-resolution acquisitions using paired t-test. Intra-class correlation coefficient (ICC) and Bland-Altman statistics (absolute and percentage differences as means ± SD) were used for assessing test-retest reproducibility. There was no significant difference between global quantitative results acquired with high vs. low-resolution T1 mapping. The reproducibility of global T1 values was good for high-resolution (ICC: 0.88) and excellent for low-resolution T1 mapping (ICC: 0.95, P=0.003). In subgroup analyses, inferior test-retest reproducibility was observed for high spatial resolution in women compared to low spatial resolution (ICC: 0.71 vs. 0.91, P=0.001) and heart rates >77 bpm (ICC: 0.53 vs. 0.88, P=0.004). Apical segments had higher T1 values and variability compared to other segments. Regional T1 values for basal (ICC: 0.81 vs. 0.89, P=0.023) and apical slices (ICC: 0.86 vs. 0.92, P=0.024) showed significantly higher reproducibility in low-resolution compared to high-resolution acquisitions but without differences for midventricular slice (ICC: 0.91 vs. 0.92, P=0.402). Based on our data, we recommend a spatial resolution on the order of 1.9 mm × 1.9 mm for native myocardial T1 mapping using a MOLLI 5(3)3 sequence at 1.5 T particularly in individuals with higher heart rates and women.

A deep unrolled neural network for real-time MRI-guided brain intervention

Accurate navigation and targeting are critical for neurological interventions including biopsy and deep brain stimulation. Real-time image guidance further improves surgical planning and MRI is ideally suited for both pre- and intra-operative imaging. However, balancing spatial and temporal resolution is a major challenge for real-time interventional MRI (i-MRI). Here, we proposed a deep unrolled neural network, dubbed as LSFP-Net, for real-time i-MRI reconstruction. By integrating LSFP-Net and a custom-designed, MR-compatible interventional device into a 3 T MRI scanner, a real-time MRI-guided brain intervention system is proposed. The performance of the system was evaluated using phantom and cadaver studies. 2D/3D real-time i-MRI was achieved with temporal resolutions of 80/732.8 ms, latencies of 0.4/3.66 s including data communication, processing and reconstruction time, and in-plane spatial resolution of 1 × 1 mm2. The results demonstrated that the proposed method enables real-time monitoring of the remote-controlled brain intervention, and showed the potential to be readily integrated into diagnostic scanners for image-guided neurosurgery.

Accuracy of pulmonary venous flow quantification by 2D phase contrast cardiac magnetic resonance

Abstract Funding Acknowledgements Type of funding sources: None. Background Determination of pulmonary venous flow by cardiac magnetic resonance (CMR) is frequently required in patients with complex congenital heart disease (CHD) lesions. Although new advanced imaging techniques like 4D flow are increasingly available, standard 2D phase contrast (2DPC) sequences remain the most frequently used method for quantification of flow. However, little is known on the accuracy of 2DPC flow measurements in the pulmonary veins. Purpose This study therefore aims to investigate the accuracy of pulmonary venous flow quantification by 2DPC. Methods Patients with CHD underwent a CMR on a 1.5T scanner. 2DPC flow sequences (typical temporal resolution 35–50ms, inplane spatial resolution 1.3–1.8 mm2, VENC 80–100) were acquired in the ascending aorta (AA), main, right and left pulmonary arteries (MPA, RPA, LPA), as well as left (LPV) and right (RPV) pulmonary veins. Pulmonary venous flow was compared to pulmonary arterial flow in all patients and to aortic flow in patients without left to right (L-R) shunt using Bland-Altman plots and Pearson’s correlations. Results Twenty CHD patients (age 27±3, range 7; 53 years, 7 women, tetralogy of Fallot N=7, Ross operation N=1; Fontan N=4, transposition of great arteries N=2, atrial septal defect N=2, Turner syndrome 1, Ebstein N=1, valvulopathies N=2) were included. Table 1 shows absolute flow values of the different vessels. Bland-Altman analyses and correlations (Table 2) showed a good agreement of total pulmonary venous (PV) to flow in the ascending aorta (group without shunt) and the sum of flow of both branch pulmonary arteries. However, agreement of total PV with MPA flow as well as of LPV and RPV with the respective branch pulmonary artery flow was reduced. Conclusions Total PV flow determined by 2DPC shows a good agreement with flow in the ascending aorta and the sum of the branch pulmonary arteries. Flow turbulences due to the pathologies of the right ventricular outflow tract probably contribute to the reduced agreement of total PV and MPA flow. However, the reduced agreement of the flow in left/right pulmonary veins and their respective branch pulmonary arteries underlines the need for newer advanced flow quantification techniques to improve accuracy of pulmonary venous flow measurements.

Inexpensive multi-plane particle image velocimetry based on defocusing: Proof of concept on two-component measurement

This paper presents a method for simultaneous particle image velocimetry (PIV) on parallel planes offset in depth. The method places images from two planes onto a different half of a camera sensor by using image splitting optics with variable optical path lengths. A shallow depth of field is achieved to ensure only one plane is in focus on each half of the sensor. Without needing additional lasers, the method is designed as an inexpensive means to increase the number of measurement plane(s) of single/multi-plane PIV setups and can be combined with existing plane discrimination approaches such as polarization and wavelength. The method is useful for studying instantaneous flow correlations on different planes while retaining high in-plane spatial resolution of typical planar PIV measurement. The measurement uncertainty caused by crosstalk from out-of-focus images is discussed. Experimental results from a laminar flow rig test indicate that the average measurement error of each velocity component is lower than 0.1 pixels per time step, with a 20 mm plane separation in depth and a 35 × 54 mm2 field of view. As an application with varying background scatter and out-of-plane flow motions, in-cylinder flow measurements in an optically accessible internal combustion engine were performed on two swirl planes simultaneously. Characteristics of the proposed method performing stereoscopic PIV measurements will be studied in future.

PO-05-037 EFFICACY AND CLINICAL CHARACTERISTICS OF AVERT-VT: ABLATION AT VIRTUAL-HEART PREDICTED VT

We have previously demonstrated the feasibility of virtual-heart simulations in determining noninvasively the optimal infarct-related ventricular tachycardia (VT) ablation targets and guiding the clinical ablation. We are now proceeding with AVERT-VT, an FDA-approved study for ten prospective patients. To develop the end-to-end personalized computational modeling prediction pipeline enabling AVERT-VT and to assess the corresponding clinical characteristics of the predicted virtual-heart ablation targets. We developed a comprehensive end-to-end computational pipeline enabling AVERT-VT (Fig). 3D LGE-MRI scans with 1.5x1.5mm in-plane spatial resolution and slice thickness <2mm were used to reconstruct the 3D left ventricular (LV) model incorporating deep scar, border (gray) zone, and non-infarcted tissue. In silico rapid pacing from 7 LV regions were performed to induce VTs. The optimal ablation targets were determined by thorough analyses of the wavefront propagation of the induced in silico VTs. A linear ablation connecting a VT ablation target to a transmural deep scar region, the basal region, or the other VT target was performed to eliminate the potential VT pathway formed by the target ablation. Potential emergent VTs in the post-ablation substrate were investigated by repeating the VT induction pacing protocol. This process was iterated until complete VT non-inducibility was achieved. The initial and repeat ablation targets, linear ablations, as well as the LV, right ventricle, infarct area, aortic arch, and left and right coronary cusp geometries, were exported as an ablation strategy plan compatible with the clinical mapping system. The figure presents results from one patient in the clinical study. The ablation plan was successfully imported and co-registered to the patient’s heart. The extent of the model infarct corresponded to the bipolar low-voltage regions. Abnormal electrograms were observed at the proposed targets. 2 VTs were observed during the procedure. Areas of latest activation in sinus rhythm and best pace maps, with a long stimulus to QRS, were found at 2 of the proposed targets. One of the seven proposed targets was not ablated due to the proximity to the left sided His potential. The patient was not inducible of any VT at the end of the procedure. AVERT-VT proposed targets corresponded to abnormal electrograms, and ablating at the targets led to VT non-inducibility.

COMPARISON OF 3D T2-WEIGHTED SEQUENCE WITH CONVENTIONAL 2D MRI SEQUENCES IN THE ASSESSMENT OF SHOULDER JOINT PATHOLOGIES

Background: Conventional 2D MR protocol routinely used in the institutions for assessment of shoulder joint pathologies have good tissue contrast and in-plane spatial resolution, but they must be repeated in several planes which takes more examination time, and the thicker slices can miss some pathologies due to partial volume averaging. Instead, procuring thin slices is possible with 3D acquisition and can be post-processed into other planes using multiplanar reformation techniques, thereby reducing overall examination time and provides good patient compliance, fewer motion artefacts. This study's goal is to determine the agreement between conventional 2D sequences and 3D T2 weighted sequence in the assessment of rotator cuff and other shoulder joint pathologies. 50 patients with suspected shoulde Methods: r joint pathologies were examined with 2D and 3D protocol in the same session. In the assessment of shoulder joint pathologies, agreement between T2-weighted 3D sequence and conventional 2D sequences is evaluated with kappa statistics. Diagnostic accuracy of T2-weighted 3D sequence with conventional 2D sequences as standard of reference is evaluated using McNemar and Chi square test. Perfe Results: ct agreement (k=1.00) between 2D and 3D sequences seen in the assessment of full thickness tears of supraspinatus and infraspinatus, partial thickness tears of subscapularis, tendinitis of rotator cuff muscles, biceps tendon, AC joint, bursitis and joint effusion. Strong agreement is seen with assessment of partial thickness tears of supraspinatus (k=0.811), moderate agreement with articular cartilage lesions (k=0.79) and bone marrow edema (k=0.779). Perfect agreement is Conclusion: seen between 2D and 3D sequences in the assessment of most of the rotator cuff and other shoulder pathologies in our series. With 3D sequence being capable of reducing overall examination time, it provides more patient compliance and reduced motion artefacts. We recommend further studies with arthroscopy as gold standard before relying purely on 3D sequence for shoulder evaluation.

Real-time observation of the precursor film for low viscosity liquids spreading on solid substrates by a customized differential laser interference microscopy.

While static wettability is well treated with Young's equation via its static contact angle, theoretical analyses for wetting dynamics are not yet reaching consensus due to a singularity of the spreading forces worked at the vapor/liquid/solid contact line. One plausible explanation to overcome the singularity problem is that there is a so-called precursor film spreading outside the apparent contact line. After its first finding in 1919, many researchers have attempted to visualize its shape. However, because its length and thickness are as small as micrometer and nanometer-order, respectively, its visualization stillremains a challenging issue especially for low-viscosity liquids. In the present study, we developed a differential laser interference microscope, which has a thickness resolution of approximately 2nm at the best, and applied it to the wetting front of 10 cSt of silicone oil spreading on a silicon wafer with an almost constant spreading velocity. As a result, the precursor film of 14µm long and 108nm thick was clearly visualized. While the macro contact line has a finite advancing contact angle of 4.0°, the gradient of the precursor film surface gradually decreased and converged to ~ 0.1° at the micro-contact angle. The shape of the precursor film was independent of the time after the dropping for the range of 600s ± 10%, which is consistent to theoretical estimation. The present study demonstrated that our interferometer simultaneously achieved nanometer thickness resolutions, micrometer in-plane spatial resolution, and at least a millisecond temporal resolution with a simple optical setup.

Single-shot spiral diffusion-weighted imaging at 7T using expanded encoding with compressed sensing.

The expanded encoding model incorporates spatially- and time-varying field perturbations for correction during reconstruction. To date, these reconstructions have used the conjugate gradient method with early stopping used as implicit regularization. However, this approach is likely suboptimal for low-SNR cases like diffusion or high-resolution MRI. Here, we investigate the extent that -wavelet regularization, or equivalently compressed sensing (CS), combined with expanded encoding improves trade-offs between spatial resolution, readout time and SNR for single-shot spiral DWI at 7T. The reconstructions were performed using our open-source graphics processing unit-enabled reconstruction toolbox, "MatMRI," that allows inclusion of the different components of the expanded encoding model, with or without CS. In vivo accelerated single-shot spirals were acquired with five acceleration factors (R) (2×-6×) and three in-plane spatial resolutions (1.5, 1.3, and 1.1 mm). From the in vivo reconstructions, we estimated diffusion tensors and computed fractional anisotropy maps. Then, simulations were used to quantitatively investigate and validate the impact of CS-based regularization on image quality when compared to a known ground truth. In vivo reconstructions revealed improved image quality with retainment of small features when CS was used. Simulations showed that the joint use of the expanded encoding model and CS improves accuracy of image reconstructions (reduced mean-squared error) over the range of R investigated. The expanded encoding model and CS regularization are complementary tools for single-shot spiral diffusion MRI, which enables both higher spatial resolutions and higher R.

Clinical commissioning of the first point-of-care spectral photon-counting CT for the upper extremities.

Acceptance testing and quality assurance (QA) of computed tomography (CT) scans are of great importance to ensure the appropriate performance of the systems. However, current standards and guidelines do not include a dedicated QA program for spectral photon-counting CT (SPCCT), nor adapted tolerancelevels. To evaluate the technical performance, in terms of image quality and radiation dose, of the first point-of-care SPCCT for the upper extremities (MARS Extremity 5X120, MARS Bioimaging Ltd., Christchurch, New Zealand) and to establish a comprehensive QAprogram. The specific dimensions of the scanner with a 125mm diameter gantry and a small voxel size of 0.1×0.1×0.1mm3 require the use of suitable phantoms and evaluation techniques. Indicators such as CT number accuracy, image noise, uniformity, and slice thickness were assessed to characterize the image quality. The in-plane and longitudinal spatial resolutions were evaluated by means of the modulation transfer function (MTF). Noise power spectra (NPS) were calculated to further evaluate the image noise. Material identification capabilities were assessed using clinically relevant high-Z materials (iodine, gold, gadolinium, and calcium). A 100-mm diameter CTDI-like phantom was used to measure the dose indices. A complete radiation survey was carried out to measure the radiation exposure at different points around thescanner. The proposed QA program is based on international and local recommendations as well as practical experience. It includes standardised CT tests and SPCCT-specific methods. Additional methodologies to further assess the system performance are also presented. Tolerance levels are discussed and revised when appropriate. Both in-plane and longitudinal high spatial resolutions were evidenced by the MTF measurements with 1.8lp· mm-1 and 5.0lp· mm-1 at 10%, respectively. The calculated effective slice thickness ranged between 0.15and 0.16mm for the five energy bins and for a reconstructed voxel size of 0.1×0.1×0.1mm3 . Reference values of the linear attenuation coefficient of water have been calculated and used to assess the CT number uniformity of water. Evaluation of the CT number accuracy and stability of various clinically relevant materials showed excellent spectral correlation and linearity between HU values and concentrations (r2 > 0.99). The NPS showed less noise correlation between slices than within transverse slice, as well as a systematic increase at low spatial frequencies. The volume CT dose index (CTDI ) for a custom-made 100mm diameter phantom was 9.32mGy. Radiation measurements around the scanner showed that it is completely shielded except for the access port, and that no additional protective measures are necessary for thepatient. A routine QA framework for SPCCT systems has been proposed. Image quality and radiation dose were assessed using newly designed phantoms, relevant metrics, and automated algorithms. Baseline values were established and tolerance levels discussed for the MARS SPCCT scanner based on collected data and internationalrecommendations.

Development of an ultrafast brain MR neuronavigation protocol for ventricular shunt placement.

Advancements in MRI technology have provided improved ways to acquire imaging data and to more seamlessly incorporate MRI into modern pediatric surgical practice. One such situation is image-guided navigation for pediatric neurosurgical procedures, including intracranial catheter placement. Image-guided surgery (IGS) requires acquisition of CT or MR images, but the former carries the risk of ionizing radiation and the latter is associated with long scan times and often requires pediatric patients to be sedated. The objective of this project was to circumvent the use of CT and standard-sequence MRI in ventricular neuronavigation by investigating the use of fast MR sequences on the basis of 3 criteria: scan duration comparable to that of CT acquisition, visualization of ventricular morphology, and image registration with surface renderings comparable to standard of care. The aim of this work was to report image development, implementation, and results of registration accuracy testing in healthy subjects. The authors formulated 11 candidate MR sequences on the basis of the standard IGS protocol, and various scan parameters were modified, such as k-space readout direction, partial k-space acquisition, sparse sampling of k-space (i.e., compressed sensing), in-plane spatial resolution, and slice thickness. To evaluate registration accuracy, the authors calculated target registration error (TRE). A candidate sequence was selected for further evaluation in 10 healthy subjects. The authors identified a candidate imaging protocol, termed presurgical imaging with compressed sensing for time optimization (PICO). Acquisition of the PICO protocol takes 25 seconds. The authors demonstrated noninferior TRE for PICO (3.00 ± 0.19 mm) in comparison with the default MRI neuronavigation protocol (3.35 ± 0.20 mm, p = 0.20). The developed and tested sequence of this work allowed accurate intraoperative image registration and provided sufficient parenchymal contrast for visualization of ventricular anatomy. Further investigations will evaluate use of the PICO protocol as a substitute for CT and conventional MRI protocols in ventricular neuronavigation.

MRI on a new polymeric multichannel membrane for ultrafiltration

Membrane ultrafiltration in new polymeric multi-channel membranes designed for in-out filtration was investigated to get insights into structure, flow and filtration properties. The apparent novelty of the membrane concerns the geometry and configuration of the feed channels.In-situmagnetic resonance imaging (MRI) allows non-invasive and non-destructive investigations with adequate spatial and time resolution. The structure of the new polymeric membrane was measured with an in-plane spatial resolution of 35 µm/pixel revealing first the polymer density distribution over the 19-channel membrane and second the wettability of the fiber and its cavities of different dimensions. MRI was also used to answer questions about flow and consequently feed distribution in the channels. Finally,in-situfiltration of an aqueous solution of sodium alginate was observed which led to deposit formation at the channel’s inner surfaces. The kinetics of this deposit formation was quantified. Backwashing and flushing gave insight into the cleanability of the channels.

Effects of spatial and temporal resolution on cardiovascular magnetic resonance feature tracking measurements using a simple realistic numerical phantom.

To develop a single-slice numerical phantom with known myocardial motion, at several temporal and in-plane spatial resolutions, for testing and comparison of Cardiovascular Magnetic Resonance (CMR) feature tracking (FT) software. The phantom was developed based on CMR acquisitions of one volunteer (acquired cine, tagging cine, T1 map, T2 map, proton density weighted image). The numerical MRI simulator JEMRIS was used, and the phantom was generated at several in-plane spatial resolutions (1.4 × 1.4 mm2 to 3.0 × 3.0 mm2) and temporal resolutions (20 to 40 cardiac phases). Two feature tracking software packages were tested: Medical Image Tracking Toolbox (MITT) and two versions of cvi42 (v5.3.8 and v5.13.7). The effect of resolution on strain results was investigated with reference to ground-truth radial and circumferential strain. Peak radial strain was consistently undermeasured more for cvi42 v5.13.7 than for v5.3.8. Increased pixel size produced a trend of increased difference from ground-truth peak strain, with the largest changes for cvi42 obtained using v5.13.7 between 1.4 × 1.4 mm2 and 3.0 × 3.0 mm2, at 9.17 percentage points (radial) and 8.42 percentage points (circumferential). The results corroborate the presence of intervendor differences in feature tracking results and show the magnitude of strain differences between software versions. This study shows how temporal and in-plane spatial resolution can affect feature tracking with reference to the ground-truth strain of a numerical phantom. Results reaffirm the need for numerical phantom development for the validation and testing of FT software.

MR microscopy of the developing upper extremity of the chicken in ovo using 7 Tesla MRI.

MR microscopy (MRM) is known as ultra-high-field (UHF) magnetic resonance imaging with an in-plane spatial resolution of <100 μm, yields highly resolved non-invasive anatomical imaging and allows longitudinal assessment of embryonic avian development. The aim of the present study was to evaluate the feasibility of in vivo anatomical MRI assessment of the developing upper extremity of the chicken. Thirty-eight fertilized chicken eggs were examined at 7 Tesla acquiring high-resolution T2-weighted images with an in-plane resolution of 74 × 74 μm. To reduce motion artefacts, the eggs were moderately cooled before and during MRI. Development of the upper extremity was anatomically and quantitatively assessed. Chondrification and ossification on MRI were correlated with histological examination. MRM allowed the identification of the embryo from stage D5 onwards. First chondrification of the upper extremity was visible at stage D7, and the differentiation of the forearm was possible from stage D9 throughout the developmental period with excellent correlation to histology. MRM also allowed the differentiation between cortical and medullary bone as well as the detection of chondrified areas. UHF MRM allows the in vivo and in ovo evaluation of the upper limb during embryonic development and provides non-invasive longitudinal anatomical information. This technique allows longitudinal studies of the same embryo during the developmental period and may therefore provide further insights into the development of the upper extremity. With improved coil technique and increasing availability of UHF MR systems, there is great potential regarding several research topics in experimental musculoskeletal radiology.