Juan E. Ortuño

Geometrical and Monte Carlo projectors in 3D PET reconstruction

PURPOSE In the present work, the authors compare geometrical and Monte Carlo projectors in detail. The geometrical projectors considered were the conventional geometrical Siddon ray-tracer (S-RT) and the orthogonal distance-based ray-tracer (OD-RT), based on computing the orthogonal distance from the center of image voxel to the line-of-response. A comparison of these geometrical projectors was performed using different point spread function (PSF) models. The Monte Carlo-based method under consideration involves an extensive model of the system response matrix based on Monte Carlo simulations and is computed off-line and stored on disk. METHODS Comparisons were performed using simulated and experimental data of the commercial small animal PET scanner rPET. RESULTS The results demonstrate that the orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions yield better images in terms of contrast and spatial resolution than those obtained after using the conventional method and the multiray-based S-RT. Furthermore, the Monte Carlo-based method yields slight improvements in terms of contrast and spatial resolution with respect to these geometrical projectors. CONCLUSIONS The orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions represent satisfactory alternatives to factorizing the system matrix or to the conventional on-the-fly ray-tracing methods for list-mode reconstruction, where an extensive modeling based on Monte Carlo simulations is unfeasible.

Efficient methodologies for system matrix modelling in iterative image reconstruction for rotating high-resolution PET

A fully 3D iterative image reconstruction algorithm has been developed for high-resolution PET cameras composed of pixelated scintillator crystal arrays and rotating planar detectors, based on the ordered subsets approach. The associated system matrix is precalculated with Monte Carlo methods that incorporate physical effects not included in analytical models, such as positron range effects and interaction of the incident gammas with the scintillator material. Custom Monte Carlo methodologies have been developed and optimized for modelling of system matrices for fast iterative image reconstruction adapted to specific scanner geometries, without redundant calculations. According to the methodology proposed here, only one-eighth of the voxels within two central transaxial slices need to be modelled in detail. The rest of the system matrix elements can be obtained with the aid of axial symmetries and redundancies, as well as in-plane symmetries within transaxial slices. Sparse matrix techniques for the non-zero system matrix elements are employed, allowing for fast execution of the image reconstruction process. This 3D image reconstruction scheme has been compared in terms of image quality to a 2D fast implementation of the OSEM algorithm combined with Fourier rebinning approaches. This work confirms the superiority of fully 3D OSEM in terms of spatial resolution, contrast recovery and noise reduction as compared to conventional 2D approaches based on rebinning schemes. At the same time it demonstrates that fully 3D methodologies can be efficiently applied to the image reconstruction problem for high-resolution rotational PET cameras by applying accurate pre-calculated system models and taking advantage of the systems symmetries.

A multichannel watershed-based algorithm for supervised texture segmentation

Segmentation of image regions based on their texture is a standard problem in image analysis. Once a set of texture features is selected, several algorithms can be applied to segment the image into regions. This paper presents an extension of the watershed algorithm using a vector gradient and multivariate region merging methods. The algorithm uses a set of texture images, and it only depends on an adjustable parameter. Results are presented on a standard set of synthetic images and on textured medical ones, using different texture parameters and merging criteria.

Real-Time Digital Timing in Positron Emission Tomography

Positron emission tomography (PET) requires accurate timing of scintillation events to properly discriminate between coincident and noncoincident pairs. The traditional solution to timing is based on custom application specific integrated circuits (ASIC) designs, whose cost may not be justified in the design of experimental small animal PET scanners. The new generation of PET scanners introduces the idea of continuous sampling of the detected scintillation pulse, replacing event-triggered acquisition front-ends. This approach enables new options to the timing procedure based on digital processing of the sampled pulse signal. This work proposes a time stamping algorithm based on the optically matched filter and compares the potential performance benefits of this approach versus other FIR-based timing algorithms, some of which have been already implemented by different authors. Results show that the coincidence timing resolution may be as low as 1.5 ns without the need of expensive high-speed converters when the proper signal processing is applied.

DCE@urLAB: a dynamic contrast-enhanced MRI pharmacokinetic analysis tool for preclinical data

BackgroundDCE@urLAB is a software application for analysis of dynamic contrast-enhanced magnetic resonance imaging data (DCE-MRI). The tool incorporates a friendly graphical user interface (GUI) to interactively select and analyze a region of interest (ROI) within the image set, taking into account the tissue concentration of the contrast agent (CA) and its effect on pixel intensity.ResultsPixel-wise model-based quantitative parameters are estimated by fitting DCE-MRI data to several pharmacokinetic models using the Levenberg-Marquardt algorithm (LMA). DCE@urLAB also includes the semi-quantitative parametric and heuristic analysis approaches commonly used in practice. This software application has been programmed in the Interactive Data Language (IDL) and tested both with publicly available simulated data and preclinical studies from tumor-bearing mouse brains.ConclusionsA user-friendly solution for applying pharmacokinetic and non-quantitative analysis DCE-MRI in preclinical studies has been implemented and tested. The proposed tool has been specially designed for easy selection of multi-pixel ROIs. A public release of DCE@urLAB, together with the open source code and sample datasets, is available at http://www.die.upm.es/im/archives/DCEurLAB/.

New embedded digital front-end for high resolution PET scanner

This work describes a new digital front-end for a high-resolution low-cost animal PET scanner which is currently under development. The advances in flexibility and size of modern FPGAs together with the release of new tools enable the integration of most of the front-end electronics in a single FPGA. The implemented system includes a small 32-bit RISC processor, several peripherals attached to the internal buses and a special DSP unit closely attached to the processor which is dedicated to the detection of the gamma events. On top of these, a small footprint real time operating system abstracts the underlying hardware, providing the mechanisms to combine on-chip slow control and data streaming.

A new PET prototype for proton therapy: comparison of data and Monte Carlo simulations

Ion beam therapy is a valuable method for the treatment of deep-seated and radio-resistant tumors thanks to the favorable depth-dose distribution characterized by the Bragg peak. Hadrontherapy facilities take advantage of the specific ion range, resulting in a highly conformal dose in the target volume, while the dose in critical organs is reduced as compared to photon therapy. The necessity to monitor the delivery precision, i.e. the ion range, is unquestionable, thus different approaches have been investigated, such as the detection of prompt photons or annihilation photons of positron emitter nuclei created during the therapeutic treatment. Based on the measurement of the induced β+ activity, our group has developed various in-beam PET prototypes: the one under test is composed by two planar detector heads, each one consisting of four modules with a total active area of 10 × 10 cm2. A single detector module is made of a LYSO crystal matrix coupled to a position sensitive photomultiplier and is read-out by dedicated frontend electronics. A preliminary data taking was performed at the Italian National Centre for Oncological Hadron Therapy (CNAO, Pavia), using proton beams in the energy range of 93–112 MeV impinging on a plastic phantom. The measured activity profiles are presented and compared with the simulated ones based on the Monte Carlo FLUKA package.

The role of elastic restoring forces in right ventricular filling

AIMS The physiological determinants of RV diastolic function remain poorly understood. We aimed to quantify the contribution of elastic recoil to RV filling and determine its sensitivity to interventricular interaction. METHODS AND RESULTS High-fidelity pressure-volume loops and simultaneous 3-dimensional ultrasound sequences were obtained in 13 pigs undergoing inotropic modulation, volume overload, and acute pressure overload induced by endotoxin infusion. Using a validated method, we isolated elastic restoring forces from ongoing relaxation using conventional pressure-volume data. The RV contracted below the equilibrium volume in >75% of the data sets. Consequently, elastic recoil generated strong sub-atmospheric passive pressure at the onset of diastole [-3 (-4 to -2) mmHg at baseline]. Stronger restoring suction pressure was related to a shorter isovolumic relaxation period, a higher rapid filling fraction, and lower atrial pressures (all P < 0.05). Restoring forces were mostly determined by the position of operating volumes around the equilibrium volume. By this mechanism, the negative inotropic effect of beta-blockade reduced and sometimes abolished restoring forces. During acute pressure overload, restoring forces initially decreased, but recovered at advanced stages. This biphasic response was related to alterations of septal curvature induced by changes in the diastolic LV-RV pressure balance. The constant of elastic recoil was closely related to the constant of passive stiffness (R = 0.69). CONCLUSION The RV works as a suction pump, exploiting contraction energy to facilitate filling by means of strong elastic recoil. Restoring forces are influenced by the inotropic state and RV conformational changes mediated by direct ventricular interdependence.

3D-OSEM transition matrix for high resolution PET imaging with modeling of the gamma-event detection

We present an efficient methodology for the calculation of the transition matrix for 3D-OSEM iterative image reconstruction, including a model of the gamma-event detection in crystal with photoelectric and Compton scatter interaction. The method is adapted for high resolution PET cameras composed of pixelated scintillator crystal arrays and with axial symmetry. 2D-OSEM algorithm, in combination with rebinning methods such as SSRB and FORE, can also be performed using a subset of this transition matrix

Projector model for efficient list-mode reconstruction in PET scanners with parallel planar detectors

We have developed a new projector model specifically tailored for fast list-mode tomographic reconstructions in Positron emission tomography (PET) scanners with parallel planar detectors. The model provides an accurate estimation of the probability distribution of coincidence events defined by pairs of scintillating crystals. This distribution is parameterized with 2D elliptical Gaussian functions defined in planes perpendicular to the main axis of the tube of response (TOR). The parameters of these Gaussian functions have been obtained by fitting Monte Carlo simulations that include positron range, acolinearity of gamma rays, as well as detector attenuation and scatter effects. The proposed model has been applied efficiently to list-mode reconstruction algorithms. Evaluation with Monte Carlo simulations over a rotating high resolution PET scanner indicates that this model allows to obtain better recovery to noise ratio in OSEM (ordered-subsets, expectation-maximization) reconstruction, if compared to list-mode reconstruction with symmetric circular Gaussian TOR model, and histogram-based OSEM with precalculated system matrix using Monte Carlo simulated models and symmetries.