Michal J. Wesolowski

A laboratory based edge-illumination x-ray phase-contrast imaging setup with two-directional sensitivity

We report on a preliminary laboratory based x-ray phase-contrast imaging system capable of achieving two-directional phase sensitivity, thanks to the use of L-shaped apertures. We show that in addition to apparent absorption, two-directional differential phase images of an object can be quantitatively retrieved by using only three input images. We also verify that knowledge of the phase derivatives along both directions allows for straightforward phase integration with no streak artefacts, a known problem common to all differential phase techniques. In addition, an analytical method for 2-directional dark field retrieval is proposed and experimentally demonstrated.

A simple method for determining split renal function from dynamic (99m)Tc-MAG3 scintigraphic data.

PurposeCommonly used methods for determining split renal function (SRF) from dynamic scintigraphic data require extrarenal background subtraction and additional correction for intrarenal vascular activity. The use of these additional regions of interest (ROIs) can produce inaccurate results and be challenging, e.g. if the heart is out of the camera field of view. The purpose of this study was to evaluate a new method for determining SRF called the blood pool compensation (BPC) technique, which is simple to implement, does not require extrarenal background correction and intrinsically corrects for intrarenal vascular activity.MethodsIn the BPC method SRF is derived from a parametric plot of the curves generated by one blood-pool and two renal ROIs. Data from 107 patients who underwent 99mTc-MAG3 scintigraphy were used to determine SRF values. Values calculated using the BPC method were compared to those obtained with the integral (IN) and Patlak-Rutland (PR) techniques using Bland-Altman plotting and Passing-Bablok regression. The interobserver variability of the BPC technique was also assessed for two observers.ResultsThe SRF values obtained with the BPC method did not differ significantly from those obtained with the PR method and showed no consistent bias, while SRF values obtained with the IN method showed significant differences with some bias in comparison to those obtained with either the PR or BPC method. No significant interobserver variability was found between two observers calculating SRF using the BPC method.ConclusionThe BPC method requires only three ROIs to produce reliable estimates of SRF, was simple to implement, and in this study yielded statistically equivalent results to the PR method with appreciable interobserver agreement. As such, it adds a new reliable method for quality control of monitoring relative kidney function.

Assessment of freeware programs for the reconstruction of tomography datasets obtained with a monochromatic synchrotron-based X-ray source

Synchrotron-based in-line phase-contrast computed tomography (PC-CT) allows soft tissue to be imaged with sub-gross resolution and has potential to be used as a diagnostic tool. The reconstruction and processing of in-line PC-CT datasets is a computationally demanding task; thus, an efficient and user-friendly software program is desirable. Four freeware programs (NRecon, PITRE, H-PITRE and Athabasca Recon) were compared for the availability of features such as dark- and flat-field calibration, beam power normalization, ring artifact removal, and alignment tools for optimizing image quality. An in-line PC-CT projection dataset (3751 projections, 180° rotation, 10.13 mm × 0.54 mm) was collected from a formalin-fixed canine prostate at the Biomedical Imaging and Therapy Bending Magnet (BMIT-BM) beamline of the Canadian Light Source. This dataset was processed with each of the four software programs and usability of the program was evaluated. Efficiency was assessed by how each program maximized computer processing power during computation. Athabasca Recon had the least-efficient memory usage, least user-friendly interface, and lacked a ring artifact removal feature. NRecon, PITRE and H-PITRE produced similar quality images, but the Athabasca Recon reconstruction suffered from the lack of a native ring remover algorithm. The 64-bit version of NRecon uses GPU (graphics processor unit) memory for accelerated processing and is user-friendly, but does not provide necessary parameters for in-line PC-CT data, such as dark-field and flat-field correction and beam power normalization. PITRE has many helpful features and tools, but lacks a comprehensive user manual and help section. H-PITRE is a condensed version of PITRE and maximizes computer memory for efficiency. To conclude, NRecon has fewer imaging processing tools than PITRE and H-PITRE, but is ideal for less experienced users due to a simple user interface. Based on the quality of reconstructed images, efficient use of computer memory and parameter availability, H-PITRE was the preferred of the four programs compared.

The early plasma concentration of 51Cr-EDTA in patients with cirrhosis and ascites: a comparison of three models.

Objectives The aim of the study was to determine which of three two-parameter fitting functions (exponential, linear-log, and negative-power function of time) most accurately models early chromium-51-EDTA (51Cr-EDTA) plasma concentration data prior to 120 min in patients with cirrhosis and ascites and understand how these fitting functions affect the calculation of the area under the plasma concentration curve (AUC). Methods A bolus, antecubital intravenous injection of 2.6 MBq of 51Cr-EDTA was given to 13 patients with cirrhosis and ascites. Up to 16 blood samples were drawn at time points ranging from 5 to 1440 min following injection. The concentration data prior to 120 min were used as reference data. Early time concentration values, estimated by fitting exponential, linear-log, and negative-power functions of time to the time samples at 120, 180, and 240 min, were then compared with reference data. The AUC was calculated for each patient using the exponential, Bröchner-Mortensen-corrected exponential, and linear-log functions, and these values were compared. Results The withheld, observed plasma concentrations were (a) most accurately estimated by linear-log functions (Wilcoxon P=0.4548), (b) significantly underestimated by exponential functions (Wilcoxon P=0.0002), and (c) significantly overestimated by negative-power functions (Wilcoxon P=0.0034). The relative errors when ranked from best to worst were those for the linear-log (12.0%, 9.0%), exponential (22.9%, 14.2%), and negative-power (31.9%, 48.4%) functions of time, respectively (median, interquartile range). For each patient, the values for AUC calculated by the exponential function differed significantly (range=3.4–15.3%, median=8.3%) from those calculated by the corrected Bröchner-Mortensen exponential, as to a lesser extent did those values calculated using linear-log functions (range=0.4–8.0%, median=3.0%). Conclusion In patients with cirrhosis, linear-log functions were significantly more accurate than exponential or power functions in estimating early time plasma concentrations (<120 min). However, the improved linear-log early time plasma concentration model does not provide as much correction to the total AUC as does the corrected Bröchner-Mortensen exponential method. This is likely because of the large contribution of late time data to the AUC, and future work is suggested to explore the late time fit problem.

Accurate and precise plasma clearance measurement using four 99mTc-DTPA plasma samples over 4 h.

ObjectivesGlomerular filtration rate can be measured as the plasma clearance (CL) of a glomerular filtration rate marker despite body fluid disturbances using numerous, prolonged time samples. We desire a simplified technique without compromised accuracy and precision. Materials and methodsWe compared CL values derived from two plasma concentration curve area methods – (a) biexponential fitting [CL (E2)] and (b) Tikhonov adaptively regularized gamma variate fitting [CL (Tk-GV)] – for 4 versus 8 h time samplings from 412 99mTc-DTPA studies in 142 patients, mostly paediatric patients, with suspected fluid disturbances. ResultsCL (Tk-GV) from four samples/4 h and from nine samples/8 h, both accurately and precisely agreed with the standard, which was taken to be nine samples/8 h CL from (noncompartmental) numerical integration [CL (NI)]. The E2 method, four samples/4 h, and nine samples/8 h median CL values significantly overestimated the CL (NI) values by 4.9 and 3.8%, respectively. ConclusionCompared with the standard, CL (E2) from four samples/4 h and from nine samples/8 h proved to be the most inaccurate and imprecise method examined, and can be replaced by better methods for calculating CL. The CL (Tk-GV) can be used to reduce sampling time in half from 8 to 4 h and from nine to four samples for a precise and accurate, yet more easily tolerated and simplified test.

Time Varying Apparent Volume of Distribution and Drug Half-Lives Following Intravenous Bolus Injections

We present a model that generalizes the apparent volume of distribution and half-life as functions of time following intravenous bolus injection. This generalized model defines a time varying apparent volume of drug distribution. The half-lives of drug remaining in the body vary in time and become longer as time elapses, eventually converging to the terminal half-life. Two example fit models were substituted into the general model: biexponential models from the least relative concentration error, and gamma variate models using adaptive regularization for least relative error of clearance. Using adult population parameters from 41 studies of the renal glomerular filtration marker 169Yb-DTPA, simulations of extracellular fluid volumes of 5, 10, 15 and 20 litres and plasma clearances of 40 and 100 ml/min were obtained. Of these models, the adaptively obtained gamma variate models had longer times to 95% of terminal volume and longer half-lives.