Blind deconvolution decreases requirements on temporal resolution of DCE-MRI: Application to 2nd generation pharmacokinetic modeling.

PURPOSE: Dynamic Contrast-Enhanced (DCE) MRI with 2nd generation pharmacokinetic models provides estimates of plasma flow and permeability surface-area product in contrast to the broadly used 1st generation models (e.g. the Tofts models). However, the use of 2nd generation models requires higher frequency with which the dynamic images are acquired (around 1.5 s per image). Blind deconvolution can decrease the demands on temporal resolution as shown previously for one of the 1st generation models. Here, the temporal-resolution requirements achievable for blind deconvolution with a 2nd generation model are studied. METHODS: The 2nd generation model is formulated as the distributed-capillary adiabatic-tissue-homogeneity (DCATH) model. Blind deconvolution is based on Parker's model of the arterial input function. The accuracy and precision of the estimated arterial input functions and the perfusion parameters is evaluated on synthetic and real clinical datasets with different levels of the temporal resolution. RESULTS: The estimated arterial input functions remained unchanged from their reference high-temporal-resolution estimates (obtained with the sampling interval around 1 s) when increasing the sampling interval up to about 5 s for synthetic data and up to 3.6-4.8 s for real data. Further increasing of the sampling intervals led to systematic distortions, such as lowering and broadening of the 1st pass peak. The resulting perfusion-parameter estimation error was below 10% for the sampling intervals up to 3 s (synthetic data), in line with the real data perfusion-parameter boxplots which remained unchanged up to the sampling interval 3.6 s. CONCLUSION: We show that use of blind deconvolution decreases the demands on temporal resolution in DCE-MRI from about 1.5 s (in case of measured arterial input functions) to 3-4 s. This can be exploited in increased spatial resolution or larger organ coverage.

In silico study of pseudoprogression in glioblastoma: collaboration of radiologists and radiation oncologists in the estimation of extent of high dose RT region.

BACKGROUND AND AIM: Oncologists play a vital role in the interpretation of radiographic results in glioblastoma patients. Molecular pathology and information on radiation treatment protocols among others are all important for accurate interpretation of radiology images. One important issue that may arise in interpreting such images is the phenomenon of tumor "pseudoprogression"; oncologists need to be able to distinguish this effect from true disease progression.Exact knowledge about the location of high-dose radiotherapy region is needed for valid determination of pseudoprogression according to RANO (Response Assessment in Neuro-Oncology) criteria in neurooncology. The aim of the present study was to evaluate the radiologists' understanding of a radiotherapy high-dose region in routine clinical practice since radiation oncologists do not always report 3-dimensional isodoses when ordering follow up imaging. METHODS: Eight glioblastoma patients who underwent postresection radiotherapy were included in this study. Four radiologists worked with their pre-radiotherapy planning MR, however, they were blinded to RT target volumes which were defined by radiation oncologists according to current guidelines. The aim was to draw target volume for high dose RT fields (that is the region, where they would consider that there may be a pseudoprogression in future MRI scans). Many different indices describing structure differences were analyzed in comparison with original per-protocol RT target volumes. RESULTS: The median volume for RT high dose field was 277 ccm (range 218 to 401 ccm) as defined per protocol by radiation oncologist and 87 ccm (range 32-338) as defined by radiologists (median difference of paired difference 31%, range 15-112%). The Median Dice index of similarity was 0.46 (range 0.14 - 0.78), the median Hausdorff distance 25 mm. CONCLUSION: Continuing effort to improve education on specific procedures in RT and in radiology as well as automatic tools for exporting RT targets is needed in order to increase specificity and sensitivity in response evaluation.

Distributed capillary adiabatic tissue homogeneity model in parametric multi-channel blind AIF estimation using DCE-MRI.

PURPOSE: One of the main challenges in quantitative dynamic contrast-enhanced (DCE) MRI is estimation of the arterial input function (AIF). Usually, the signal from a single artery (ignoring contrast dispersion, partial volume effects and flow artifacts) or a population average of such signals (also ignoring variability between patients) is used. METHODS: Multi-channel blind deconvolution is an alternative approach avoiding most of these problems. The AIF is estimated directly from the measured tracer concentration curves in several tissues. This contribution extends the published methods of multi-channel blind deconvolution by applying a more realistic model of the impulse residue function, the distributed capillary adiabatic tissue homogeneity model (DCATH). In addition, an alternative AIF model is used and several AIF-scaling methods are tested. RESULTS: The proposed method is evaluated on synthetic data with respect to the number of tissue regions and to the signal-to-noise ratio. Evaluation on clinical data (renal cell carcinoma patients before and after the beginning of the treatment) gave consistent results. An initial evaluation on clinical data indicates more reliable and less noise sensitive perfusion parameter estimates. CONCLUSION: Blind multi-channel deconvolution using the DCATH model might be a method of choice for AIF estimation in a clinical setup.

Absolute ultrasound perfusion parameter quantification of a tissue-mimicking phantom using bolus tracking [Correspondence].

This study presents three methods for absolute quantification in ultrasound perfusion analysis based on bolus tracking. The first two methods deconvolve the perfusion time sequence with a measured AIF, using a nonparametric or a parametric model of the tissue residue function, respectively. The third method is a simplified approach avoiding deconvolution by assuming a narrow AIF. A phantom with a dialyzer filter as a tissue-mimicking model was used for evaluation. Estimated mean transit times and blood volumes were compared with the theoretical values. A match with a maximum error of 12% was achieved.

The precision of DCE-MRI using the tissue homogeneity model with continuous formulation of the perfusion parameters.

The present trend in dynamic contrast-enhanced MRI is to increase the number of estimated perfusion parameters using complex pharmacokinetic models. However, less attention is given to the precision analysis of the parameter estimates. In this paper, the distributed capillary adiabatic tissue homogeneity pharmacokinetic model is extended by the bolus arrival time formulated as a free continuous parameter. With the continuous formulation of all perfusion parameters, it is possible to use standard gradient-based optimization algorithms in the approximation of the tissue concentration time sequences. This new six-parameter model is investigated by comparing Monte-Carlo simulations with theoretically derived covariance matrices. The covariance-matrix approach is extended from the usual analysis of the primary perfusion parameters of the pharmacokinetic model to the analysis of the perfusion parameters derived from the primary ones. The results indicate that the precision of the estimated perfusion parameters can be described by the covariance matrix for signal-to-noise ratio higher than~20dB. The application of the new analysis model on a real DCE-MRI data set is also presented.

Ultrasound perfusion analysis combining bolus-tracking and burst-replenishment.

A new signal model and processing method for quantitative ultrasound perfusion analysis is presented, called bolus-and-burst. The method has the potential to provide absolute values of blood flow, blood volume, and mean transit time. Furthermore, it provides an estimate of the local arterial input function which characterizes the arterial tree, allowing accurate estimation of the bolus arrival time. The method combines two approaches to ultrasound perfusion analysis: bolus-tracking and burst-replenishment. A pharmacokinetic model based on the concept of arterial input functions and tissue residue functions is used to model both the bolus and replenishment parts of the recording. The pharmacokinetic model is fitted to the data using blind deconvolution. A preliminary assessment of the new perfusion-analysis method is presented on clinical recordings.