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.

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 system's symmetries.

Temas Avanzados en Teoría Electromagnética

Colección de monografías sobre temas avanzados en el área de los sistemas de guiado de ondas. El contenido es adecuado para estudiantes que, habiendo seguido un curso básico, sientan especial predilección por el tema, para aquellos que inician su formación de postgrado en el área de electromagnetismo, y para titulados con experiencia que precisen, en un momento dado, una información más completa acerca de ciertos problemas ajenos a su actividad previa.

New diagnostic and therapeutic approaches to treat ventricular tachycardias originating at the summit of the left ventricle role of merged hemodynamic-MRI and alternative ablation sources

The left ventricular (LV) summit is the most common site of idiopathic epicardial LV arrhythmias and frequently represents a diagnostic and a therapeutic challenge.1 We present a case of sustained monomorphic ventricular tachycardia (SMVT) originating at the LV summit that underwent failed cryosurgical epicardial ablation and was successfully treated with the aid of merged hemodynamic and contrast-enhanced MRI (CE-MRI).

Massively parallelizable list-mode reconstruction using a Monte Carlo-based elliptical Gaussian model

Purpose: A fully three-dimensional (3D) massively parallelizable list-mode ordered-subsets expectation-maximization (LM-OSEM) reconstruction algorithm has been developed for high-resolution PET cameras. System response probabilities are calculated online from a set of parameters derived from Monte Carlo simulations. The shape of a system response for a given line of response (LOR) has been shown to be asymmetrical around the LOR. This work has been focused on the development of efficient region-search techniques to sample the system response probabilities, which are suitable for asymmetric kernel models, including elliptical Gaussian models that allow for high accuracy and high parallelization efficiency. The novel region-search scheme using variable kernel models is applied in the proposed PET reconstruction algorithm. Methods: A novel region-search technique has been used to sample the probability density function in correspondence with a small dynamic subset of the field of view that constitutes the region of response (ROR). The ROR is identified around the LOR by searching for any voxel within a dynamically calculated contour. The contour condition is currently defined as a fixed threshold over the posterior probability, and arbitrary kernel models can be applied using a numerical approach. The processing of the LORs is distributed in batches among the available computing devices, then, individual LORs are processed within different processing units. In this way, both multicore and multiple many-core processing units can be efficiently exploited. Tests have been conducted with probability models that take into account the noncolinearity, positron range, and crystal penetration effects, that produced tubes of response with varying elliptical sections whose axes were a function of the crystal's thickness and angle of incidence of the given LOR. The algorithm treats the probability model as a 3D scalar field defined within a reference system aligned with the ideal LOR. Results: This new technique provides superior image quality in terms of signal-to-noise ratio as compared with the histogram-mode method based on precomputed system matrices available for a commercial small animal scanner. Reconstruction times can be kept low with the use of multicore, many-core architectures, including multiple graphic processing units. Conclusions: A highly parallelizable LM reconstruction method has been proposed based on Monte Carlo simulations and new parallelization techniques aimed at improving the reconstruction speed and the image signal-to-noise of a given OSEM algorithm. The method has been validated using simulated and real phantoms. A special advantage of the new method is the possibility of defining dynamically the cut-off threshold over the calculated probabilities thus allowing for a direct control on the trade-off between speed and quality during the reconstruction.

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.

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

Background DCE@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. Results Pixel-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. Conclusions A 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/ webcite.

Utilidad de los mapas de intensidad de señal basados en secuencias de realce tardía 3D para caracterizar el sustrato arritmogénico en pacientes con taquicardía ventricular isquémica

La escara miocárdica constituye el sustrato arritmogénico de los pacientes con TV isquémica. Los mapas de intensidad de señal (IS) obtenidos a partir de mágenes de realce tardío (RT) pueden ser una alternativa no invasiva a los mapas electroanatómicos para estudiar la arquitectura de la escara. El objetivo de este trabajo es describir las características de las escaras observadas en los pacientes (p) con TV isquémicas.

Registration and fusion of contrast-enhanced MRI myocardial substrate maps and X-ray angiograms

The aim of this work is to provide the necessary methods to register and fuse the endo-epicardial signal intensity (SI) maps extracted from contrast-enhanced magnetic resonance imaging (ceMRI) with X-ray coronary ngiograms using an intrinsic registrationbased algorithm to help pre-planning and guidance of catheterization procedures. Fusion of angiograms with SI maps was treated as a 2D-3D pose estimation, where each image point is projected to a Plücker line, and the screw representation for rigid motions is minimized using a gradient descent method. The resultant transformation is applied to the SI map that is then projected and fused on each angiogram. The proposed method was tested in clinical datasets from 6 patients with prior myocardial infarction. The registration procedure is optionally combined with an iterative closest point algorithm (ICP) that aligns the ventricular contours segmented from two ventriculograms.

Impacto de la geometría del septo en el fenómeno de succión del ventrículo derecho

El impacto del fenómeno de interdependencia ventricular en la función diastólica del VD no es bien conocido. El VD es capaz de facilitar el llenado generando presión (P) negativa al inicio de la diástole pero los mecanismos implicados en este fenómeno están aún por dilucidar. Nos propusimos analizar la contribución de las fuerzas de retroceso elástico al llenado del VD y su relación con la geometría del septo.

Assessment of diastolic chamber properties of the right ventricle by global fitting of pressure-volume data and conformational analysis of 3D + T echocardiographic sequences

Assessment of diastolic chamber properties of the right ventricle by global fitting of pressure-volume data and conformational analysis of 3D + T echocardiographic sequences

Uso de circuitos equivalentes y teoría de filtros para el diseño de superficies selectivas en frecuencia

Frequency selective surfaces (FSS) and reflect/trasmitarrays are mostly designed on the basis ot optimization using an electromagnetic simulator. That is a time consuming method and some decisions have to be taken using simply intuition. The use of equivalent circuits of the scatterers selected for the design allows the substitution of the intuition and most part of the optimization process by the application of the classic rules of filter design. This communication presents all the steps necessary to obtain the equivalent circuit of different square scatterers in a periodic lattice and to implement the desired FSS frequency behaviour calculating the number of layers and the dimensions of the periodic cells. Several examples are included to evaluate the results.

Estudio de hormigón autocompactante con árido reciclado

La tesis doctoral “Estudio de hormigón autocompactante con árido reciclado” realizada dentro del programa de doctorado de la Universidad Politécnica de Madrid “Máster en técnicas experimentales avanzadas en la ingeniería civil”, investiga la sustitución de áridos gruesos naturales por reciclados en hormigones autocompactantes, para demostrar la posibilidad de utilización de este tipo de árido en la fabricación de hormigones autocompactantes. En cuanto a la línea experimental adoptada, la primera fase corresponde a la caracterización de los cementos y de los áridos naturales y reciclados. En ella se han obtenido las principales características físicas y mecánicas. Una vez validadas las características de todos los materiales y adoptada una dosificación de hormigón autocompactante, se han elaborado cuatro dosificaciones con cuatro grados de incorporación de árido reciclado cada una, y una dosificación con seis grados de incorporación de árido reciclado. Fabricándose un total de 22 tipos de hormigón diferentes, sin contar todas las amasadas iniciales hasta la consecución de un hormigón autocompactante. Las cinco dosificaciones se han dividido en dos grupos para poder analizar con mayor grado de definición las características de cada uno. El primer grupo es aquel que contienen los hormigones con diferentes relaciones a/c, que incluye a la muestra A (a/c=0.55), muestra D (a/c=0.50) y muestra E (a/c=0.45). Por el contrario, el segundo grupo dispone de una relación fija de a/c=0.45 pero diferentes relaciones a/c efectivas, ya que algunas de las muestras disponen de un contenido de agua que permite contrarrestar la mayor absorción del árido reciclado. Estando en este grupo la muestra E (sin agua adicional), la muestra H (con presaturación de los áridos) y la muestra I (con un aporte de agua junto con el agua de amasado. Una vez fabricados los hormigones, se pasa a la segunda fase del estudio correspondiente a la caracterización del hormigón en estado fresco. En esta fase se han llevado a cabo los ensayos de escurrimiento, escurrimiento con anillo japonés, ensayo embudo en V y embudo V a los 5 minutos. Todos estos ensayos permiten evaluar la autocompactabilidad del hormigón según el anejo 17 de la EHE-08. La tercera fase del estudio se centra en la caracterización de los hormigones en estado endurecido, evaluando las características resistentes del hormigón. Para ello, se han realizado los ensayos de resistencia a compresión, a tracción, módulo de elasticidad y coeficiente de Poisson. En la cuarta y última fase, se han analizado la durabilidad de los hormigones, debido que a pesar de ofrecer una adecuada autocompactabilidad y resistencia mecánica, se debe de obtener un hormigón con una correcta durabilidad. Para tal fin, se ha determinado la resistencia a la penetración de agua bajo presión y carbonatación de las probetas. Este último ensayo se ha realizado teniendo en cuenta las condiciones del denominado método natural, con una exposición al ambiente de 90 días y 365 días. Con todos estos resultados se elaboraron las conclusiones derivadas de la investigación, demostrándose la posibilidad de fabricación de hormigones autocompactantes con árido reciclado (HACR) con sustituciones de hasta un 40%, e incluso dependiendo de la relación a/c con sustituciones del 60% y el 80%. ABSTRACT The doctoral thesis titled Analysis of self-compacting concrete with recycled aggregates, has been developed in accordance with the doctoral program: Master degree in advanced experimental techniques in civil engineering, at UPM. It investigates the possibility of replacing natural coarse aggregates with recycled coarse aggregates, in the field of self-compacting concrete. The aim of this dissertation is to analyze the possibility of using recycled coarse aggregates in the manufacture process of self-compacting concretes. Regarding the experimental part, the first phase refers to mechanical and physical characterization of some materials such as cement, natural aggregates and recycled aggregates. Once the characteristics of all materials have been validated and the mixing proportions have been adopted, four different mixes are elaborated by using four dosage rates of recycled aggregates in each one of the samples. Moreover, an additional sample consisting of six different dosages of recycled aggregates is considered. A total number of 22 concrete specimens have been manufactured, without including all the initial kneading samples used to obtain this type of self-compacting concrete. The aforementioned mixes have been divided in two different groups to be able to analyze with more definition. The first group is the one in which the concrete contains different values of the water - cement ratio. It includes the next samples: A (w/c=0.55), D (w/c=0.50) and E (w/c=0.45). The second group has a fixed water -cement ratio, w/c=0.45, but a different effective water - cement ratio, since some of the samples have a water content that enables to offset the major absorption of the recycled aggregates, being in this group the mixing E (without additional water), the mixing H (with saturated recycled aggregate) and the mixing I (with an additional water content to the existing kneading water). Once the concrete samples have been manufactured, the following section deals with the characterization of the concrete in fresh conditions. To accomplish this, several characterization tests are carried out such as the slump-flow test, test slump flow with Japanese ring, test V-funnel and V-funnel to 5 minutes. These tests are used to assess the self-compacting conditions according to the annex 17 of the EHE-08 The third phase of the study focuses on the mechanical characterization, the assessment of the strength properties of the concrete such as compressive strength, tensile strength, modulus of elasticity and Poisson´s ratio. Within the fourth and last phase, durability of the concrete is evaluated. This fact is motivated by the need to obtain not only good self-compacting and mechanical strength properties, but also adequate durability conditions. To accomplish the aforementioned durability, resistance of the samples under certain conditions such as water penetration pressure and carbonation, has been obtained. The latter test has been carried out taking into account the natural method, with an exhibition period to the environment of 90 days and 365 days. Through the results coming from this research work, it has been possible to obtain the main conclusions. It has been demonstrated the possibility to manufacture self-compacting concrete by using recycled aggregates with replacement rates up to 40% or, depending on the w/c ratio, rates of 60% and 80% might be reached.

Análisis electromagnético de células TEM mediante técnicas de ajuste de campo

Una célula TEM es un sistema concebido para conseguir en una determinada zona del espacio un campo electromagnético de forma controlada cuyo comportamiento tienda lo más posible al de una onda plana, aprovechando para ello ciertas propiedades asociadas al modo TEM. En este proyecto se analizan, desde el punto de vista electromagnnético, algunos tipos de células TEM. Tras una breve descripción de las diferentes geometrías que éstas pueden tomar y de las aplicaciones asociadas a las mismas, se estudian algunas de las técnicas de ajuste de campo (Resonancia Transversal y Ajuste Modal) con el fin de obtener los conocimientos necesarios para implementar una herramienta computacional que permita la determinación electromagnética completa "Full-Wave"del problema planteado. Una vez desarrollada, se emplearán los resultados derivados de la misma para analizar ciertos fenómenos que se producen en este tipo de sistemas y que nos aportará un conocimiento profundo sobre los mismos con vista a concebir estrategias de diseño que permitan optimizar su funcionamiento.

Design of broadband transmitarray unit cells with comparative study of different numbers of layers

A generalized methodology to design low-profile transmitarray (TA) antennas made of several stacked layers with nonresonant printed phasing elements is presented. A study of the unit cell bandwidth, phase-shift range and tolerances has been conducted considering different numbers of layers. A structure with three metalized layers with capacitive and inductive elements enabling a phase range of nearly 360° and low insertion loss is introduced. A study of the four-layer structure shows improvement in the performance of the unit cells in terms of bandwidth from 2% to more than 20% and a complete phase coverage. Implementations on a flexible substrate of TAs with progressive phase shift operating at 19 GHz are used for validation.