A novel alternative to classify tissues from <Emphasis Type="Italic">T</Emphasis><Subscript>1</Subscript> and <Emphasis Type="Italic">T</Emphasis><Subscript>2</Subscript> relaxation times for prostate MRI

Objective

To segment and classify the different attenuation regions from MRI at the pelvis level using the T ₁ and T ₂ relaxation times and anatomical knowledge as a first step towards the creation of PET/MR attenuation maps.

Materials and methods

Relaxation times were calculated by fitting the pixel-wise intensities of acquired T ₁- and T ₂-weighted images from eight men with inversion-recovery and multi-echo multi-slice spin-echo sequences. A decision binary tree based on relaxation times was implemented to segment and classify fat, muscle, prostate, and air (within the body). Connected component analysis and an anatomical knowledge-based procedure were implemented to localize the background and bone.

Results

Relaxation times at 3 T are reported for fat (T ₁ = 385 ms, T ₂ = 121 ms), muscle (T ₁ = 1295 ms, T ₂ = 40 ms), and prostate (T ₁ = 1700 ms, T ₂ = 80 ms). The relaxation times allowed the segmentation–classification of fat, prostate, muscle, and air, and combined with anatomical knowledge, they allowed classification of bone. The good segmentation–classification of prostate [mean Dice similarity score (mDSC) = 0.70] suggests a viable implementation in oncology and that of fat (mDSC = 0.99), muscle (mDSC = 0.99), and bone (mDSCs = 0.78) advocates for its implementation in PET/MR attenuation correction.

Conclusion

Our method allows the segmentation and classification of the attenuation-relevant structures required for the generation of the attenuation map of PET/MR systems in prostate imaging: air, background, bone, fat, muscle, and prostate.

     

3D <Emphasis Type="Italic">T</Emphasis><Subscript>2</Subscript>-weighted imaging at 7T using dynamic k<Subscript>T</Subscript>-points on single-transmit MRI systems

Objectives

For turbo spin echo (TSE) sequences to be useful at ultra-high field, they should ideally employ an RF pulse train compensated for the B ₁ ⁺ inhomogeneity. Previously, it was shown that a single k_T-point pulse designed in the small tip-angle regime can replace all the pulses of the sequence (static k_T-points). This work demonstrates that the B ₁ ⁺ dependence of T ₂-weighted imaging can be further mitigated by designing a specific k_T-point pulse for each pulse of a 3D TSE sequence (dynamic k_T-points) even on single-channel transmit systems

Materials and methods

By combining the spatially resolved extended phase graph formalism (which calculates the echo signals throughout the sequence) with a gradient descent algorithm, dynamic k_T-points were optimized such that the difference between the simulated signal and a target was minimized at each echo. Dynamic k_T-points were inserted into the TSE sequence to acquire in vivo images at 7T.

Results

The improvement provided by the dynamic k_T-points over the static k_T-point design and conventional hard pulses was demonstrated via simulations. Images acquired with dynamic k_T-points showed systematic improvement of signal and contrast at 7T over regular TSE—especially in cerebellar and temporal lobe regions without the need of parallel transmission.

Conclusion

Designing dynamic k_T-points for a 3D TSE sequence allows the acquisition of T ₂-weighted brain images on a single-transmit system at ultra-high field with reduced dropout and only mild residual effects due to the B ₁ ⁺ inhomogeneity.

     

Imaging lung function using rapid dynamic acquisition of <Emphasis Type="Italic">T</Emphasis><Subscript>1</Subscript>-maps during oxygen enhancement

This paper describes imaging of lung function with oxygen-enhanced MRI using dynamically acquired T ₁ parameter maps, which allows an accurate, quantitative assessment of time constants of T ₁-enhancement and therefore lung function. Eight healthy volunteers were examined on a 1.5-T whole-body scanner. Lung T ₁-maps based on an IR Snapshot FLASH technique (TE = 1.4 ms, TR = 3.5 ms, FA = 7 ^∘) were dynamically acquired from each subject. Without waiting for full relaxation between subsequent acquisition of T ₁-maps, one T ₁-map was acquired every 6.7 s. For comparison, all subjects underwent a standard pulmonary function test (PFT). Oxygen wash-in and wash-out time course curves of T ₁ relaxation rate (R ₁)-enhancement were obtained and time constants of oxygen wash-in (w _in) and wash-out (w _out) were calculated. Averaged over the whole right lung, the mean w _out was 43.90 ± 10.47 s and the mean (w _in) was 51.20 ± 15.53 s, thus about 17% higher in magnitude. Wash-in time constants correlated strongly with forced expired volume in one second in percentage of the vital capacity (FEV₁ % VC) and with maximum expiratory flow at 25% vital capacity (MEF25), whereas wash-out time constants showed only weak correlation. Using oxygen-enhanced rapid dynamic acquisition of T ₁-maps, time course curves of R ₁-enhancement can be obtained. With w _in and w _out two new parameters for assessing lung function are available. Therefore, the proposed method has the potential to provide regional information of pulmonary function in various lung diseases.     

The impact of fibre orientation on <Emphasis Type="Italic">T</Emphasis><Subscript>1</Subscript>-relaxation and apparent tissue water content in white matter

Objective

Recent MRI studies have shown that the orientation of nerve fibres relative to the main magnetic field affects the R₂*(= 1/T₂*) relaxation rate in white matter (WM) structures. The underlying physical causes have been discussed in several studies but are still not completely understood. However, understanding these effects in detail is of great importance since this might serve as a basis for the development of new diagnostic tools and/or improve quantitative susceptibility mapping techniques. Therefore, in addition to the known angular dependence of R₂*, the current study investigates the relationship between fibre orientation and the longitudinal relaxation rate, R₁ (= 1/T₁), as well as the apparent water content.

Materials and methods

For a group of 16 healthy subjects, a series of gradient echo, echo-planar and diffusion weighted images were acquired at 3T from which the decay rates, the apparent water content and the diffusion direction were reconstructed. The diffusion weighted data were used to determine the angle between the principle fibre direction and the main magnetic field to examine the angular dependence of R₁ and apparent water content.

Results

The obtained results demonstrate that both parameters depend on the fibre orientation and exhibit a positive correlation with the angle between fibre direction and main magnetic field.

Conclusion

These observations could be helpful to improve and/or constrain existing biophysical models of brain microstructure by imposing additional constraints resulting from the observed angular dependence R₁ and apparent water content in white matter.

     

Torque ripple minimization with optimized <Emphasis Type="Italic">i</Emphasis><Subscript><Emphasis Type="Italic">sx</Emphasis></Subscript> estimation in vector control of current source inverter induction motor

In variable speed applications of high-power induction motors, current source inverters are usually used instead of voltage source inverters. Power switches delay and switching losses in high-power ranges, prevents inverters to operate properly in high frequencies. So the voltage and current waveforms of the motor supply consist of high harmonic content causing a high-torque ripple. One of the parameters that influences the torque ripple most is the stator flux reference value. Motors do not usually operate in their nominal load. Therefore, the optimum value of the flux is not equal to the nominal value. In this paper, an optimal i _sx reference for vector control system of current source inverter fed induction motor is produced by applying a simple estimator. Simulation results show that torque ripple has been reduced, especially in the loads less than their nominal power, as significantly as power factor has been improved. Also flux reduction causes to noticeable reduction in core and switching losses.     

Estimation of metabolite <Emphasis Type="Italic">T</Emphasis><Subscript>1</Subscript> relaxation times using tissue specific analysis,signal averaging and bootstrapping from magnetic resonance spectroscopic imaging data

Object A novel method of estimating metabolite T ₁ relaxation times using MR spectroscopic imaging (MRSI) is proposed. As opposed to conventional single-voxel metabolite T ₁ estimation methods, this method investigates regional and gray matter (GM)/white matter (WM) differences in metabolite T ₁ by taking advantage of the spatial distribution information provided by MRSI. Material and methods The method, validated by Monte Carlo studies, involves a voxel averaging to preserve the GM/WM distribution, a non-linear least squares fit of the metabolite T ₁ and an estimation of its standard error by bootstrapping. It was applied in vivo to estimate the T ₁ of N-acetyl compounds (NAA), choline, creatine and myo-inositol in eight normal volunteers, at 1.5 T, using a short echo time 2D-MRSI slice located above the ventricles. Results WM-T _1,NAA was significantly (P < 0.05) longer in anterior regions compared to posterior regions of the brain. The anterior region showed a trend of a longer WM T ₁ compared to GM for NAA, creatine and myo-Inositol. Lastly, accounting for the bootstrapped standard error estimate in a group mean T ₁ calculation yielded a more accurate T ₁ estimation. Conclusion The method successfully measured in vivo metabolite T ₁ using MRSI and can now be applied to diseased brain.     

Intravascular contrast agent <Emphasis Type="Italic">T1</Emphasis> shortening: fast <Emphasis Type="Italic">T1</Emphasis> relaxometry in a carotid volunteer study

OBJECTIVE: To measure the T1 times in blood after the administration of the intravascular contrast agent gadofosveset trisodium in humans. MATERIALS AND METHODS: In a pilot study for parameter optimization, the T1-shortening induced by the injection of a single dose (0.03 mmol/kg body weight) of the MR contrast agent Vasovist (Bayer Schering Pharma AG) was measured at B (0) = 1.5 T as a function of time. In four sessions, T1 measurements were performed in the carotid vein of 9 volunteers up to 30 min after injection. T1 times were measured using a segmented saturation recovery turboFLASH (SSRTFL) pulse sequence with 7 different saturation recovery delay times in a total acquisition time of 20 s. RESULTS: The SSRTFL measurements showed T1 times of about 100 ms immediately after injection, which gradually increased to 175 ms at 30 min. The time curve of the R1 = 1/T1 averaged over all volunteers could be described with an exponential decay with a time constant T = 330 +/- 65 s and an amplitude DeltaR1 = 4.1 +/- 0.3 s(-1), and a constant offset of R1(0) = 5.7 +/- 0.2 s(-1). Mean relaxation values are in excellent agreement with theoretical predictions. CONCLUSION: An analytical expression for the initial T1-shortening of Vasovist was derived which can now be used for optimization of the pulse sequence parameters in clinical studies.     

Dielectric properties and defect mechanisms of (1-<Emphasis Type="Italic">x</Emphasis>)Ba(Fe<Subscript>0.5</Subscript>Nb<Subscript>0.5</Subscript>)O<Subscript>3</Subscript> -<Emphasis Type="Italic">x</Emphasis>BiYbO<Subscript>3</Subscript> ceramics

(1-x)Ba(Fe_0.5Nb_0.5)O₃ -xBiYbO₃ (BFN-xBY) ceramics were prepared by a conventional solid-state reaction method. The dielectric properties and relaxation behavior of BFN-xBY ceramics were analyzed according to dielectric and impedance spectroscopy. Dielectric permittivity of the ceramics increases with increasing temperature below 500 K then remains unchanged up to 700 K, while corresponding loss factor decreases with the increase of temperature below 500 K then increase slowly. Defect compensation mechanism of this system was analyzed in detail. The giant dielectric behavior of the ceramics arises from the internal barrier layer capacitor (IBLC) effect. Polarization effect at insulating grain boundaries between semiconducting grains accompanied by a strong Maxwell-Wagner (MW) relaxation mode. The characteristic of grain boundaries was revealed using impedance spectroscope and the universal dielectric response law.     

Towards accurate and precise <Emphasis Type="Italic">T</Emphasis><Subscript>1</Subscript> and extracellular volume mapping in the myocardium: a guide to current pitfalls and their solutions

Mapping of the longitudinal relaxation time (T ₁) and extracellular volume (ECV) offers a means of identifying pathological changes in myocardial tissue, including diffuse changes that may be invisible to existing T ₁-weighted methods. This technique has recently shown strong clinical utility for pathologies such as Anderson-Fabry disease and amyloidosis and has generated clinical interest as a possible means of detecting small changes in diffuse fibrosis; however, scatter in T ₁ and ECV estimates offers challenges for detecting these changes, and bias limits comparisons between sites and vendors. There are several technical and physiological pitfalls that influence the accuracy (bias) and precision (repeatability) of T ₁ and ECV mapping methods. The goal of this review is to describe the most significant of these, and detail current solutions, in order to aid scientists and clinicians to maximise the utility of T ₁ mapping in their clinical or research setting. A detailed summary of technical and physiological factors, issues relating to contrast agents, and specific disease-related issues is provided, along with some considerations on the future directions of the field.     

Identification of myocardial diffuse fibrosis by 11 heartbeat MOLLI <Emphasis Type="Italic">T</Emphasis><Subscript>1</Subscript> mapping: averaging to improve precision and correlation with collagen volume fraction

Objectives

Our objectives involved identifying whether repeated averaging in basal and mid left ventricular myocardial levels improves precision and correlation with collagen volume fraction for 11 heartbeat MOLLI T ₁ mapping versus assessment at a single ventricular level.

Materials and methods

For assessment of T ₁ mapping precision, a cohort of 15 healthy volunteers underwent two CMR scans on separate days using an 11 heartbeat MOLLI with a 5(3)3 beat scheme to measure native T ₁ and a 4(1)3(1)2 beat post-contrast scheme to measure post-contrast T ₁, allowing calculation of partition coefficient and ECV. To assess correlation of T ₁ mapping with collagen volume fraction, a separate cohort of ten aortic stenosis patients scheduled to undergo surgery underwent one CMR scan with this 11 heartbeat MOLLI scheme, followed by intraoperative tru-cut myocardial biopsy. Six models of myocardial diffuse fibrosis assessment were established with incremental inclusion of imaging by averaging of the basal and mid-myocardial left ventricular levels, and each model was assessed for precision and correlation with collagen volume fraction.

Results

A model using 11 heart beat MOLLI imaging of two basal and two mid ventricular level averaged T ₁ maps provided improved precision (Intraclass correlation 0.93 vs 0.84) and correlation with histology (R ² = 0.83 vs 0.36) for diffuse fibrosis compared to a single mid-ventricular level alone. ECV was more precise and correlated better than native T ₁ mapping.

Conclusion

T ₁ mapping sequences with repeated averaging could be considered for applications of 11 heartbeat MOLLI, especially when small changes in native T ₁/ECV might affect clinical management.

     

Stage substitution <Emphasis Type="Italic">E</Emphasis>-<Emphasis Type="Italic">H</Emphasis> circuits for calculating magnetoelectric induction heating plants

An approximate analytical method is proposed to calculate the electromagnetic field of a magnetoelectric induction heating plant using a stage substitution circuit formed based on a comparison of the general laws and methods of the computation of electromagnetic field in a conducting medium and methods of the circuit theory for the sinusoidal current.     

Generalized 12<Emphasis Type="Italic">n</Emphasis>-phase converter

Topological features and calculation relationships are considered for a 12n-phase converter with a multistage configuration of six-phase asymmetric voltage systems of secondary transformer windings. An example of a synthesized circuit of the phase number and shift converter is given for a 24-pulse rectifier recommended for use in networks with distorting factors.     

Hyperpolarized <Superscript>3</Superscript>He lung ventilation imaging with <Emphasis Type="Italic">B</Emphasis><Subscript>1</Subscript>-inhomogeneity correction in a single breath-hold scan

The signal-to-noise ratio (SNR) of hyperpolarized noble gas MR images is sensitive to the flip angle used. Variations in flip angle due to B ₁-inhomogeneity of the RF coil cause intensity variation artifacts in lung ventilation images which may mask or mimic disease. We show these artifacts can be minimized by using the optimal flip angle and corrected if the local flip angle is known. Hyperpolarized ³He lung images were obtained in ten healthy subjects using both a conventional gradient-echo sequence and a new hybrid pulse sequence designed to simultaneously acquire lung ventilation images and corresponding flip-angle maps in comparable imaging time. Flip-angle maps and corrected images were calculated from the hybrid scan and compared with conventional images. The qualitative theoretical dependence of flip angle on SNR was verified. Ventilation images and flip-angle maps were successfully obtained with the hybrid sequence. Corrections to image intensity calculated from the flip-angle maps appeared reasonable for images acquired using an average flip angle near optimal. Use of the optimal flip angle is crucial to the quality of lung ventilation images. Artifactual intensity variations due to RF-coil inhomogeneity may be identified and potentially corrected using our hybrid sequence. Acknowledgement The authors thank John M. Christopher, RT(R)(MR) for valuable assistance performing the studies, and Jaime F. Mata, MS, Jing Cai, MS, and Andrew G. Reish, BS, for excellent operation of the ³He polarization system. We gratefully acknowledge support for this research from the Commonwealth of Virginia Technology Research Fund (Grant No. IN2002-01), Siemens Medical Solutions and Amersham Health.     

The Applicability of Sr-deficient <Emphasis Type="Italic">n</Emphasis>-type SrTiO<Subscript>3</Subscript> for SOFC Anodes

Here we discuss the effect of preparation conditions on structural stability and electrical properties of Sr-deficient n-type SrTiO₃. In particular, an explanation of a wide scatter of conductivity values in Y- and Nb-doped SrTiO₃. reported in the literature is proposed, based on the existing defect chemistry model of n-doped SrTiO₃. It was confirmed that when sintered in air, Sr-deficient SrTiO₃ doped with Nb and/or Y, remains single phase until the solubility limit (e.g., 30% for Nb or 4% for Y). However, when sintered at low po₂, the material transforms from a vacancy compensated to an electronically compensated compound with a strontium deficient second phase. Measured at 800°C in low po₂, the maximum conductivity of these multi-phase compounds was 340 S/cm and 100 S/cm for the Nb-doped and Y-doped sample, respectively. However, the conductivity dropped dramatically to less than 10 S/cm when samples of the same compositions were sintered in air, again measured in reducing atmosphere.     

3D single point imaging with compressed sensing provides high temporal resolution <Emphasis Type="Italic">R</Emphasis><Subscript>2</Subscript>* mapping for in vivo preclinical applications

Objective

Purely phase-encoded techniques such as single point imaging (SPI) are generally unsuitable for in vivo imaging due to lengthy acquisition times. Reconstruction of highly undersampled data using compressed sensing allows SPI data to be quickly obtained from animal models, enabling applications in preclinical cellular and molecular imaging.

Materials and methods

TurboSPI is a multi-echo single point technique that acquires hundreds of images with microsecond spacing, enabling high temporal resolution relaxometry of large-R ₂* systems such as iron-loaded cells. TurboSPI acquisitions can be pseudo-randomly undersampled in all three dimensions to increase artifact incoherence, and can provide prior information to improve reconstruction. We evaluated the performance of CS-TurboSPI in phantoms, a rat ex vivo, and a mouse in vivo.

Results

An algorithm for iterative reconstruction of TurboSPI relaxometry time courses does not affect image quality or R ₂* mapping in vitro at acceleration factors up to 10. Imaging ex vivo is possible at similar acceleration factors, and in vivo imaging is demonstrated at an acceleration factor of 8, such that acquisition time is under 1 h.

Conclusions

Accelerated TurboSPI enables preclinical R ₂* mapping without loss of data quality, and may show increased specificity to iron oxide compared to other sequences.

     

Low Temperature Sintering of BaTi<Subscript>4</Subscript>O<Subscript>9</Subscript>-Based Middle-<Emphasis Type="Italic">k</Emphasis> Dielectric Composition for LTCC Applications

Effect of glass addition on the low-temperature sintering and microwave dielectric properties of BaTi₄O₉-based ceramics were studied to develop the middle-k dielectric composition for the functional substrate of low-temperature co-fired ceramics. When 10 wt% of glass was added, sufficient densification was obtained and the relative density more than 98% was reached at the sintering temperature of 875C. The microwave dielectric properties were k = 32, Q × f = 9000 GHz, and t_cf = 10 ppm/C. As the added amount of glass frit with base dielectric composition, phase changes from BaTi₄O₉ to BaTi₅O₁₁ and Ba₄Ti₁₃O₃₀ was observed, which result in the modification of microwave dielectric properties.     

On high suppression of NO<Subscript><Emphasis Type="Italic">x</Emphasis></Subscript> and CO emissions in gas-turbine plants with combined gas-and-steam cycles

In this work are given results of analyzing processes of production of nitrogen oxides (NO _x ) and afterburning of CO when firing natural gas at combined-cycle gas-turbine plants. It is shown that for suppressing emissions of the said microcomponents it is necessary to lower temperature in hot local zones of the flame in which NO_x is formed, and, in so doing, to avoid chilling of cold flame zones that prevents afterburning of CO. The required lowering of the combustion temperature can be provided by combustion of mixtures of methane with steam, with high mixing uniformity that ensures the same and optimum fraction of the steam “ballast” in each microvolume of the flame. In addition to chilling, the steam ballast makes it possible to maintain a fairly high concentration of hydroxil radicals in the flame zone as well, and this provides high burning out of fuel and reduction in carbon monoxide emissions (active steam ballast). Due to this fact the fraction of steam when firing its mixtures with methane in a gas-turbine plant can be increased up to the weight ratio 4: 1. In this case, the concentrations of NO _x and CO in emissions can be reduced to ultra-low values (less than 3 ppm).     

Correction of <Emphasis Type="Italic">B</Emphasis><Subscript>0</Subscript>-induced geometric distortion variations in prospective motion correction for 7T MRI

Objective

Prospective motion correction can effectively fix the imaging volume of interest. For large motion, this can lead to relative motion of coil sensitivities, distortions associated with imaging gradients and B ₀ field variations. This work accounts for the B ₀ field change due to subject movement, and proposes a method for correcting tissue magnetic susceptibility-related distortion in prospective motion correction.

Materials and methods

The B ₀ field shifts at the different head orientations were characterized. A volunteer performed large motion with prospective motion correction enabled. The acquired data were divided into multiple groups according to the object positions. The correction of B ₀-related distortion was applied to each group of data individually via augmented sensitivity encoding with additionally integrated gradient nonlinearity correction.

Results

The relative motion of the gradients, B ₀ field and coil sensitivities in prospective motion correction results in residual spatial distortion, blurring, and coil artifacts. These errors can be mitigated by the proposed method. Moreover, iterative conjugate gradient optimization with regularization provided superior results with smaller RMSE in comparison to standard conjugate gradient.

Conclusion

The combined correction of B ₀-related distortion and gradient nonlinearity leads to a reduction of residual motion artifacts in prospective motion correction data.

     

Pulsed arterial spin labelling at ultra-high field with a <Emphasis Type="Italic">B</Emphasis><Stack><Subscript>1</Subscript><Superscript>+</Superscript></Stack>-optimised adiabatic labelling pulse

Objective

Arterial spin labelling (ASL) techniques benefit from the increased signal-to-noise ratio and the longer T ₁ relaxation times available at ultra-high field. Previous pulsed ASL studies at 7 T concentrated on the superior regions of the brain because of the larger transmit radiofrequency inhomogeneity experienced at ultra-high field that hinders an adequate inversion of the blood bolus when labelling in the neck. Recently, researchers have proposed to overcome this problem with either the use of dielectric pads, through dedicated transmit labelling coils, or special adiabatic inversion pulses.

Materials and methods

We investigate the performance of an optimised time-resampled frequency-offset corrected inversion (TR-FOCI) pulse designed to cause inversion at much lower peak B ₁ ⁺ . In combination with a PICORE labelling, the perfusion signal obtained with this pulse is compared against that obtained with a FOCI pulse, with and without dielectric pads.

Results

Mean grey matter perfusion with the TR-FOCI was 52.5 ± 10.3 mL/100 g/min, being significantly higher than the 34.6 ± 2.6 mL/100 g/min obtained with the FOCI pulse. No significant effect of the dielectric pads was observed.

Conclusion

The usage of the B ₁ ⁺ -optimised TR-FOCI pulse results in a significantly higher perfusion signal. PICORE–ASL is feasible at ultra-high field with no changes to operating conditions.

     

Comparison of methods for estimation of the intravoxel incoherent motion (IVIM) diffusion coefficient (<Emphasis Type="Italic">D</Emphasis>) and perfusion fraction (<Emphasis Type="Italic">f</Emphasis>)

Objective

Intravoxel incoherent motion (IVIM) shows great potential in many applications, e.g., tumor tissue characterization. To reduce image-quality demands, various IVIM analysis approaches restricted to the diffusion coefficient (D) and the perfusion fraction (f) are increasingly being employed. In this work, the impact of estimation approach for D and f is studied.

Materials and methods

Four approaches for estimating D and f were studied: segmented IVIM fitting, least-squares fitting of a simplified IVIM model (sIVIM), and Bayesian fitting of the sIVIM model using marginal posterior modes or posterior means. The estimation approaches were evaluated in terms of bias and variability as well as ability for differentiation between tumor and healthy liver tissue using simulated and in vivo data.

Results

All estimation approaches had similar variability and ability for differentiation and negligible bias, except for the Bayesian posterior mean of f, which was substantially biased. Combined use of D and f improved tumor-to-liver tissue differentiation compared with using D or f separately.

Discussion

The similar performance between estimation approaches renders the segmented one preferable due to lower numerical complexity and shorter computational time. Superior tissue differentiation when combining D and f suggests complementary biologically relevant information.