A straightforward procedure to build a non-toxic relaxometry phantom with desired T1 and T2 times at 3T.

OBJECTIVE: To prepare and analyze soy-lecithin-agar gels for non-toxic relaxometry phantoms with tissue-like relaxation times at 3T. METHODS: Phantoms mimicking the relaxation times of various tissues (gray and white matter, kidney cortex and medulla, spleen, muscle, liver) were built and tested with a clinical 3T whole-body MR scanner. Simple equations were derived to calculate the appropriate concentrations of soy lecithin and agar in aqueous solutions to achieve the desired relaxation times. Phantoms were tested for correspondence between measurements and calculated T1 and T2 values, reproducibility, spatial homogeneity, and temporal stability. T1 and T2 mapping techniques and a 3D T1-weighted sequence with high spatial resolution were applied. RESULTS: Except for the liver relaxation phantom, all phantoms were successfully and reproducibly produced. Good agreement was found between the targeted and measured relaxation times. The percentage deviations from the targeted relaxation times were less than 3% for T1 and less than 6.5% for T2. In addition, the phantoms were homogeneous and had little to no air bubbles. However, the phantoms were unstable over time: after a storage period of 4 weeks, mold growth and also changes in relaxation times were detected in almost all phantoms. CONCLUSION: Soy-lecithin-agar gels are a non-toxic material for the construction of relaxometry phantoms with tissue-like relaxation times. They are easy to prepare, inexpensive and allow independent adjustment of T1 and T2. However, there is still work to be done to improve the long-term stability of the phantoms.

Detectability and accuracy of computational measurements of in-silico and physical representations of enlarged perivascular spaces from magnetic resonance images.

BACKGROUND: Magnetic Resonance Imaging (MRI) visible perivascular spaces (PVS) have been associated with age, decline in cognitive abilities, interrupted sleep, and markers of small vessel disease. But the limits of validity of their quantification have not been established. NEW METHOD: We use a purpose-built digital reference object to construct an in-silico phantom for addressing this need, and validate it using a physical phantom. We use cylinders of different sizes as models for PVS. We also evaluate the influence of 'PVS' orientation, and different sets of parameters of the two vesselness filters that have been used for enhancing tubular structures, namely Frangi and RORPO filters, in the measurements' accuracy. RESULTS: PVS measurements in MRI are only a proxy of their true dimensions, as the boundaries of their representation are consistently overestimated. The success in the use of the Frangi filter relies on a careful tuning of several parameters. Alpha= 0.5, beta= 0.5 and c= 500 yielded the best results. RORPO does not have these requirements and allows detecting smaller cylinders in their entirety more consistently in the absence of noise and confounding artefacts. The Frangi filter seems to be best suited for voxel sizes equal or larger than 0.4 mm-isotropic and cylinders larger than 1 mm diameter and 2 mm length. 'PVS' orientation did not affect measurements in data with isotropic voxels. COMPARISON WITH EXISTENT METHODS: Does not apply. CONCLUSIONS: The in-silico and physical phantoms presented are useful for establishing the validity of quantification methods of tubular small structures.

A novel anthropomorphic breathing phantom with a pneumatic MR-safe actuator for tissue deformation studies during MRI and radiotherapy.

A novel MR-safe anthropomorphic torso phantom with an MR-conditional pneumatic respiration system and inflatable lungs for tissue deformation studies is proposed. The phantom consists of a pair of lungs made from sponges encased in flexible polyurethane. The lung phantom also contains a set of silicone tubes of various diameters to mimic the larger vasculature and airways of the lungs. The lungs are surrounded by a plastic ribcage and immersed in a gelatine hydrogel within a flexible polyurethane skin. A plastic pneumatic pump system was constructed to inflate and deflate the lungs. A fibre optic rotary encoder was constructed to determine the volume of displaced air in the lungs. The pneumatic pump and rotary encoder were constructed of plastic materials to allow placement within the bore of the MR scanner with minimal interaction with the magnetic field. Breath-gated scans and rapid imaging scans (2.5 s per image) were taken of the phantom in the stationary state and during inflation/deflation, and with Cartesian and BLADE k-space sampling. It was found that BLADE shows the least motion artifacts during breathing. This phantom and respiration system shows potential for quality assurance of MRI incorporating breathing corrections and for radiotherapy applications in tracking a moving target.

A multi-modality medical imaging head and neck phantom: Part 1. Design and fabrication.

A novel anthropomorphic head and neck phantom which features the emulation of blood flow and perfusion is proposed. The phantom is helpful in both education and research and contains major blood vessels and a porous silicone elastomer brain compartment with microvascular capillary flow. The porous brain compartment is fabricated by use of a novel cast-moulding-dissolution technique. The skull and vertebra are fabricated by a combination of 3D printing and cast-moulding and are tissue equivalent with CT numbers ranging from 1000 HU to 1200 HU. The elastic structure of the phantom allows ultrasound imaging in the neck region. MRI compatible pressure sensors measure the pressure in the carotic arteries and the jugular veins and pulsatile flow is created by use of a peristaltic pump. The pressure-flow dynamics are physiologically realistic and also matches well with computational simulations of porous Darcy flow. The phantom can be used to optimize and validate MRI pulse sequences and protocols for flow imaging, MR angiography, Arterial Spin Labeling (ASL) and dynamic contrast enhanced (DCE) MRI.

A multi-modality medical imaging head and neck phantom: Part 2. Medical imaging.

The head and neck phantom discussed in an accompanying paper (part 1), is imaged with MRI, X-ray CT, PET and ultrasound. MRI scans show a distinct image contrast between the brain compartment and other anatomical regions of the head. The silicone matrix that was used to create a porous brain compartment has a relatively high proton density and a spin-spin relaxation time (T2) that is long enough to provide an MRI signal. While the longitudinal magnetization was found to recover according to a mono-exponential, a bi-exponential decay was observed for the transverse relaxation with a slow T2 relaxation component corresponding to the perfusate and a fast T2 relaxation component corresponding to the silicone. The fraction of the slow T2 relaxation component increases upon perfusion. A dynamic contrast enhanced (DCE) MRI experiment is conducted in which the injection rate of the contrast agent is varied. Parametric DCE maps are created and reveal regional differences in contrast agent kinetics as a result of differences in porosity. The skull, vertebra and the brain compartment are clearly visible on X-ray CT. Dynamic PET scanning has been performed while the carotic arterial input function is monitored by use of a Geiger-Müller counter. Similar regions of perfusion are found in the PET study as in the DCE MRI study. By doping the perfusate with a lipid micelle emulsion, the phantom is applicable for carotic Doppler ultrasound demonstration and validation.

Technical Note: Human tissue-equivalent MRI phantom preparation for 3 and 7 Tesla.

PURPOSE: While test objects (phantoms) in magnetic resonance imaging (MRI) are crucial for sequence development, protocol validation, and quality control, studies on the preparation of phantoms have been scarce, particularly at fields exceeding 3 Tesla. Here, we present a framework for the preparation of phantoms with well-defined T1 and T2 times at 3 and 7 Tesla. METHODS: Phantoms with varying concentrations of agarose and Gd-DTPA were prepared and measured at 3 and 7 Tesla using T1 and T2 mapping techniques. An empirical, polynomial model was constructed that best represents the data at both field strengths, enabling the preparation of new phantoms with specified combinations of both T1 and T2 . Instructions for three different tissue types (brain gray matter, brain white matter, and renal cortex) are presented and validated. RESULTS: T1 times in the samples ranged from 698 to 2820 ms and from 695 to 2906 ms, whereas T2 times ranged from 39 to 227 ms and from 34 to 235 ms for 3 and 7 Tesla scans, respectively. Models for both relaxation times used six parameters to represent the data with an adjusted R² of 0.998 and 0.997 for T1 and T2 , respectively. CONCLUSION: Based on the equations derived from the current study, it is now possible to obtain accurate weight specifications for a test object with desired T1 and T2 relaxation times. This will spare researchers the laborious task of trail-and-error approaches in test object preparation attempts.

A physical phantom for amine chemical exchange saturation transfer (CEST) MRI.

OBJECTIVE: To develop a robust amine chemical exchange saturation transfer (CEST) physical phantom, validate the temporal stability, and create a supporting software for automatic image processing and quality assurance. MATERIALS AND METHODS: The phantom was designed as an assembled laser-cut acrylic rack and 18 vials of phantom solutions, prepared with different pHs, glycine concentrations, and gadolinium concentrations. We evaluated glycine concentrations using ultraviolet absorbance for 70 days and measured the pH, relaxation rates, and CEST contrast for 94 days after preparation. We used Spearman's correlation to determine if glycine degraded over time. Linear regression and Bland-Altman analysis were performed between baseline and follow-up measurements of pH and MRI properties. RESULTS: No degradation of glycine was observed (p > 0.05). The pH and MRI measurements stayed stable for 3 months and showed high consistency across time points (R2 = 1.00 for pH, R1, R2, and CEST contrast), which was further validated by the Bland-Altman plots. Examples of automatically generated reports are provided. DISCUSSION: We designed a physical phantom for amine CEST-MRI, which is easy to assemble and transfer, holds 18 different solutions, and has excellent short-term chemical and MRI stability. We believe this robust phantom will facilitate the development of novel sequences and cross-scanners validations.

Longitudinal acquisition repeatability of MRI radiomics features: An ACR MRI phantom study on two MRI scanners using a 3D T1W TSE sequence.

PURPOSE: The purpose of this study was to quantitatively assess the longitudinal acquisition repeatability of MRI radiomics features in a three-dimensional (3D) T1-weighted (T1W) TSE sequence via a well-controlled prospective phantom study. METHODS: Thirty consecutive daily datasets of an ACR-MRI phantom were acquired on two 1.5T MRI simulators using a 3D T1W TSE sequence. Images were blindly segmented by two observers. Post-acquisition processing was minimized but an intensity discretization (fixed bin size of 25). One hundred and one radiomics features (shape n = 12; first order n = 16; texture n = 73) were extracted. Longitudinal repeatability of each feature was evaluated by Pearson correlation and coefficient of variance (CV68% ). Interobserver feature value agreement was also quantified using intraclass correlation coefficient (ICC) and Bland-Altman analysis. A most repeatable radiomics feature set on both scanners was determined by feature coefficient of variance (CV68% <5%), ICC (>0.75), and the ratio of the interobserver difference to the interobserver mean δ<5%. RESULTS: No trend of radiomics feature value changed with time. Longitudinal feature repeatability CV68% ranged 0.01-38.60% (mean/median: 12.5%/9.9%), and 0.01-40.47%, (8.49%/7.34%) on the scanners A and B. Shape features exhibited significantly better repeatability than first-order and texture features (all P < 0.01). Significant longitudinal repeatability difference was observed in texture features (P < 0.001) between the two scanners, but not in shape and first-order features (P > 0.30). First-order and texture features had smaller interobserver-dependent variation than acquisition-dependent variation. They also showed good interobserver agreement on both scanners (A:ICC = 0.80 ± 0.23; B:ICC = 0.80 ± 0.22), independent of acquisition repeatability. The repeatable radiomics features in common on both scanners, including 12 shape features, 0 first-order features, and 3 texture features, were determined as the most repeatable MRI radiomics feature set. CONCLUSIONS: Radiomics features exhibited heterogeneous longitudinal repeatability, while the shape features were the most repeatable, in this phantom study with a 3D T1W TSE acquisition. The most repeatable radiomics feature set derived in this study should be helpful for the selection of reliable radiomics features in the future clinical use.

Automatic Analysis of ACR Phantom Images in MRI.

BACKGROUND: Quality Assurance (QA) of Magnetic Resonance Imaging (MRI) system is an essential step to avoid problems in diagnosis when image quality is low. It is considered a patient safety issue. The accreditation program of the American College of Radiology (ACR) includes a standardized image quality measurement protocol. However, it has been shown that human testing by visual inspection is not objective and not reproducible. METHODS: The overall goal of the present paper was to develop and implement a fully automated method for accurate image analysis to increase its objectivity. It can positively impact the QA process by decreasing the reaction time, improving repeatability, and by reducing operator dependency. The proposed QA procedures were applied to ten clinical MRI scanners. The performance of the automated procedure was assessed by comparing the test results with the decisions made by trained MRI technologists according to ACR guidelines. The p-value, correlation coefficient of the manual and automatic measurements were also computed using the Pearson test. RESULTS AND CONCLUSION: Compared to the manual process, the use of the proposed approach can significantly reduce the time requirements while maintaining consistency with manual measurements and furthermore, decrease the subjectivity of the results. Accordingly, a strong correlation was found and the corresponding p-value was much lower than the significance level of 0.05 indicating a good agreement between the two measurements.

Cortical bone vessel identification and quantification on contrast-enhanced MR images.

BACKGROUND: Cortical bone porosity is a major determinant of bone strength. Despite the biomechanical importance of cortical bone porosity, the biological drivers of cortical porosity are unknown. The content of cortical pore space can indicate pore expansion mechanisms; both of the primary components of pore space, vessels and adipocytes, have been implicated in pore expansion. Dynamic contrast-enhanced MRI (DCE-MRI) is widely used in vessel detection in cardiovascular studies, but has not been applied to visualize vessels within cortical bone. In this study, we have developed a multimodal DCE-MRI and high resolution peripheral QCT (HR-pQCT) acquisition and image processing pipeline to detect vessel-filled cortical bone pores. METHODS: For this in vivo human study, 19 volunteers (10 males and 9 females; mean age =63±5) were recruited. Both distal and ultra-distal regions of the non-dominant tibia were imaged by HR-pQCT (82 µm nominal resolution) for bone structure segmentation and by 3T DCE-MRI (Gadavist; 9 min scan time; temporal resolution =30 sec; voxel size 230×230×500 µm3) for vessel visualization. The DCE-MRI was registered to the HR-pQCT volume and the voxels within the MRI cortical bone region were extracted. Features of the DCE data were calculated and voxels were categorized by a 2-stage hierarchical kmeans clustering algorithm to determine which voxels represent vessels. Vessel volume fraction (volume ratio of vessels to cortical bone), vessel density (average vessel count per cortical bone volume), and average vessel volume (mean volume of vessels) were calculated to quantify the status of vessel-filled pores in cortical bone. To examine spatial resolution and perform validation, a virtual phantom with 5 channel sizes and an applied pseudo enhancement curve was processed through the proposed image processing pipeline. Overlap volume ratio and Dice coefficient was calculated to measure the similarity between the detected vessel map and ground truth. RESULTS: In the human study, mean vessel volume fraction was 2.2%±1.0%, mean vessel density was 0.68±0.27 vessel/mm3, and mean average vessel volume was 0.032±0.012 mm3/vessel. Signal intensity for detected vessel voxels increased during the scan, while signal for non-vessel voxels within pores did not enhance. In the validation phantom, channels with diameter 250 µm or greater were detected successfully, with volume ratio equal to 1 and Dice coefficient above 0.6. Both statistics decreased dramatically for channel sizes less than 250 µm. CONCLUSIONS: We have a developed a multi-modal image acquisition and processing pipeline that successfully detects vessels within cortical bone pores. The performance of this technique degrades for vessel diameters below the in-plane spatial resolution of the DCE-MRI acquisition. This approach can be applied to investigate the biological systems associated with cortical pore expansion.

Gellan gum-based gels with tunable relaxation properties for MRI phantoms.

OBJECT: The research follows the analysis of gellan gum-based gels as novel MRI phantom material with the implementation of a design of experiments model to obtain tunable relaxation properties. MATERIALS AND METHODS: Gellan gum gels doped with newly synthesized superparamagnetic iron oxide nanoparticles (SPIONs) and either MnCl2 or GdCl3 were prepared and scanned from 230 µT to 3 T. Nineteen gel samples were formulated with varying concentrations of contrast agents to determine the linear, quadratic, and interactive effects of the contrast agents by a central composite design of experiment. To inhibit microbial growth in the gels and to enable long-term use, methyl 4­hydroxybenzoate (methylparaben) was utilized. RESULTS: The model containing SPIONs and metal salts relaxivity was analyzed with ANOVA, and the resulting significant coefficients were tabulated. The mathematical model was able to accurately predict the intended relaxation property from the concentration of the contrast agent with adjusted R2 values > 0.97 for longitudinal (R1) relaxation rates and 0.87 for transverse (R2) relaxation rates. CONCLUSION: The gel material maintained physical, chemical, and biological stability for at least four months and contained controllable relaxation properties while maintaining optical clarity.

Image quality assessment of a 1.5T dedicated magnetic resonance-simulator for radiotherapy with a flexible radio frequency coil setting using the standard American College of Radiology magnetic resonance imaging phantom test.

BACKGROUND: A flexible RF coil setting has to be used on an MR-simulator (MR-sim) in the head and neck simulation scan for radiotherapy (RT) purpose, while the image quality might be compromised due to the sub-optimized flexible coil compared to the normal diagnostic radiological (DR) head coil. In this study, we assessed the image quality of an MR-sim by conducting the standard American College of Radiology (ACR) MRI phantom test on a 1.5T MR-sim under RT-setting and comparing it to DR-setting. METHODS: A large ACR MRI phantom was carefully positioned, aligned and scanned 9 times for each under RT- and DR-setting on a 1.5T MR-sim, following the ACR scanning instruction. Images were analyzed following the ACR guidance. Measurement results under two coil settings were quantitatively compared. Inter-observer disagreements under RT-setting between two physicists were compared using Bland-Altman (BA) analysis and intra-class correlation coefficient (ICC). RESULTS: The MR-sim with RT-setting obtained sufficiently good image quality to pass all ACR recommended criteria. No significant difference was found in phantom length accuracy, high-contrast spatial resolution, slice thickness accuracy, slice position accuracy, and percent-signal ghosting. RT-setting significantly under-performed in low-contrast object detectability, while better performed in image intensity uniformity. BA analysis showed that 95% limit of agreement and biases of phantom test measurement under RT-setting between two observers were very small. Excellent inter-observer agreement (ICC >0.75) was achieved in all measurements except for slice thickness accuracy (ICC =0.42, moderate agreement) under RT-setting. CONCLUSIONS: Very good and highly reproducible image quality could be achieved on a 1.5T MR-sim with a flexible coil setting as revealed by the standard ACR MRI phantom test. The flexible RT-setting compromised in image signal-to-noise ratio (SNR) compared to the normal DR-setting, and resulted in reduced low-contrast object detectability.

Isotropic reconstruction of a 4-D MRI thoracic sequence using super-resolution.

PURPOSE: Four-dimensional (4D) thoracic magnetic resonance imaging (MRI) sequences have been shown to successfully monitor both tumor and lungs anatomy. However, a high temporal resolution is required to avoid motion artifacts, which leads to volumes with poor spatial resolution. This article proposes to reconstruct an isotropic 4D MRI thoracic sequence with minimum modifications to the acquisition protocols. This could be an important step toward the use of 4D MRI for thoracic radiotherapy applications. METHODS: In a postacquisition step, three orthogonal 4D anisotropic acquisitions are combined using super-resolution to reconstruct a series of isotropic volumes. A new phantom that simulates lung tumor motion is developed to evaluate the performance of the algorithm. The proposed framework is also applied to real data of a lung cancer patient. RESULTS: Subjective and objective evaluations show clear resolution enhancement and partial volume effect diminution. The isotropic reconstruction of patient data significantly improves both the visualization and segmentation of thoracic structures. CONCLUSIONS: The results presented here are encouraging and suggest that super-resolution can be regarded as an efficient method to improve the resolution of 4D MRI sequences. It produces an isotropic 4D sequence that would be impossible to acquire in practice. Further investigations will be required to evaluate its reproducibility in various clinical applications.