Surfing the optimization space of a multiple-GPU parallel implementation of a X-ray tomography reconstruction algorithm

The increasing popularity of massively parallel architectures based on accelerators have opened up the possibility of significantly improving the performance of X-ray computed tomography (CT) applications towards achieving real-time imaging. However, achieving this goal is a challenging process, as most CT applications have not been designed for exploiting the amount of parallelism existing in these architectures. In this paper we present the massively parallel implementation and optimization of Mangoose(++), a CT application for reconstructing 3D volumes from 20 images collected by scanners based on cone-beam geometry. The main contribution of this paper are the following. First, we develop a modular application design that allows to exploit the functional parallelism inside the application and to facilitate the parallelization of individual application phases. Second, we identify a set of optimizations that can be applied individually and in combination for optimally deploying the application on a massively parallel multi-GPU system. Third, we present a study of surfing the optimization space of the modularized application and demonstrate that a significant benefit can be obtained from employing the adequate combination of application optimizations. (C) 2014 Elsevier Inc. All rights reserved.This work was partially funded by the Spanish Ministry of Science and Technology under the grant TIN2010-16497, the AMIT project (CEN-20101014) from the CDTI-CENIT program, RECAVA-RETIC Network (RD07/0014/2009), projects TEC2010-21619-C04-01, TEC2011-28972-C02-01, and PI11/00616 from the Spanish Ministerio de Ciencia e Innovacion, ARTEMIS program (S2009/DPI-1802), from the Comunidad de Madrid

A Modular Approach to Lung Nodule Detection from Computed Tomography Images Using Artificial Neural Networks and Content Based Image Representation

Lung cancer is one of the most lethal cancer types. Research in computer aided detection (CAD) and diagnosis for lung cancer aims at providing effective tools to assist physicians in cancer diagnosis and treatment to save lives. In this dissertation, we focus on developing a CAD framework for automated lung cancer nodule detection from 3D lung computed tomography (CT) images. Nodule detection is a challenging task that no machine intelligence can surpass human capability to date. In contrast, human recognition power is limited by vision capacity and may suffer from work overload and fatigue, whereas automated nodule detection systems can complement expert’s efforts to achieve better detection performance. The proposed CAD framework encompasses several desirable properties such as mimicking physicians by means of geometric multi-perspective analysis, computational efficiency, and the most importantly producing high performance in detection accuracy. As the central part of the framework, we develop a novel hierarchical modular decision engine implemented by Artificial Neural Networks. One advantage of this decision engine is that it supports the combination of spatial-level and feature-level information analysis in an efficient way. Our methodology overcomes some of the limitations of current lung nodule detection techniques by combining geometric multi-perspective analysis with global and local feature analysis. The proposed modular decision engine design is flexible to modifications in the decision modules; the engine structure can adopt the modifications without having to re-design the entire system. The engine can easily accommodate multi-learning scheme and parallel implementation so that each information type can be processed (in parallel) by the most adequate learning technique of its own. We have also developed a novel shape representation technique that is invariant under rigid-body transformation and we derived new features based on this shape representation for nodule detection. We implemented a prototype nodule detection system as a demonstration of the proposed framework. Experiments have been conducted to assess the performance of the proposed methodologies using real-world lung CT data. Several performance measures for detection accuracy are used in the assessment. The results show that the decision engine is able to classify patterns efficiently with very good classification performance

The NASA SBIR product catalog

The purpose of this catalog is to assist small business firms in making the community aware of products emerging from their efforts in the Small Business Innovation Research (SBIR) program. It contains descriptions of some products that have advanced into Phase 3 and others that are identified as prospective products. Both lists of products in this catalog are based on information supplied by NASA SBIR contractors in responding to an invitation to be represented in this document. Generally, all products suggested by the small firms were included in order to meet the goals of information exchange for SBIR results. Of the 444 SBIR contractors NASA queried, 137 provided information on 219 products. The catalog presents the product information in the technology areas listed in the table of contents. Within each area, the products are listed in alphabetical order by product name and are given identifying numbers. Also included is an alphabetical listing of the companies that have products described. This listing cross-references the product list and provides information on the business activity of each firm. In addition, there are three indexes: one a list of firms by states, one that lists the products according to NASA Centers that managed the SBIR projects, and one that lists the products by the relevant Technical Topics utilized in NASA's annual program solicitation under which each SBIR project was selected

Advanced capabilities for planar X-ray systems

Mención Internacional en el título de doctorThe past decades have seen a rapid evolution towards the use of digital detectors in radiology and a more flexible robotized movement of the system components, X-ray tube and detector. This evolution opened the possibility for incorporating advanced capabilities in these planar X-ray systems, and for providing new valuable diagnostic information compared to the previous technology. Some of the current challenges for radiography are to obtain more quantitative images and to reduce the inherent superposition of tissues because of the 2D nature of the technique. Dual energy radiography, based on the acquisition of two images at different source voltages, enables a separate characterization of soft tissue and bone structures. Its benefits over conventional radiography have been proven in different applications, since it improves information content without adding significant extra acquisition time or radiation dose. In a different direction, a really disruptive advance would be to obtain 3D imaging with systems designed just for planar images. The incorporation of tomographic capabilities into these systems would have to deal with the acquisition of a limited number of projections, with non-standard geometrical configurations. This thesis presents original contributions in these two directions: dual energy radiography and 3D imaging with X-ray systems designed for planar imaging. The work is framed in a line of research of the Biomedical Imaging and Instrumentation Group from the Bioengineering and Aerospace Department of University Carlos III de Madrid working jointly with the University Hospital Gregorio Marañón, focused on the advance of radiology systems. This research line is carried out in collaboration with the group of Computer Architecture, Communications and Systems (ARCOS), from the same university, the Imaging Research Laboratory (IRL) of the University of Washington and the research center CREATIS, France. The research has a clear focus on technology transfer to the industry through the company Sedecal, a Spanish multinational among the 10 best world companies in the medical imaging field. The first contribution of this thesis is a complete novel protocol to incorporate dual energy capabilities that enable quantitative planar studies. The proposal is based on the use of a preliminary calibration with a very simple and low-cost phantom formed by two parts that represent soft tissue and bone equivalent materials. This calibration is performed automatically with no strict placement requirements. Compared to current Dual-energy X-ray Absorptiometry (DXA) systems, 1) it provides real mass-thickness values directly, enabling quantitative planar studies instead of relative comparisons, and 2) it is based on an automatic preliminary calibration without the need of interaction of an experienced technician. The second contribution is a novel protocol for the incorporation of tomographic capabilities into X-ray systems originally intended for planar imaging. For this purpose, we faced three main challenges. First, the geometrical trajectory of equipment follows non-standard circular orbits, thus posing severe difficulties for reconstruction. To handle this, the proposed protocol comprises a new geometrical calibration procedure that estimates all the system parameters per-projection. Second, the reconstruction of a limited number of projections from a reduced angular span leads to severe artifacts when using conventional reconstruction methods. To deal with these limited-view data, the protocol includes a novel advanced reconstruction method that incorporates the surface information of the sample, which can be extracted with a 3D light surface scanner. These data are introduced as an imposed constraint following the Split Bregman formulation. The restriction of the search space by exploiting the surface-based support becomes crucial for a complete recovery of the external contour of the sample and surroundings when the angular span is extremely reduced. The modular, efficient and flexible design followed for its implementation allows for the reconstruction of limited-view data with non-standard trajectories. Third, the optimization of the acquisition protocols has not yet explored with these systems. This thesis includes a study of the optimum acquisition protocols that allowed us to identify the possibilities and limitations of these planar systems. Using the surface-constrained method, it is possible to reduce the total number of projections up to 33% and the angular span down to 60 degrees. The contributions of this thesis open the way to provide depth and quantitative information very valuable for the improvement of radiological diagnosis. This could impact considerably the clinical practice, where conventional radiology is still the imaging modality most used, accounting for 80-90% of the total medical imaging exams. These advances open the possibility of new clinical applications in scenarios where 1) the reduction of the radiation dose is key, such as lung cancer screening or Pediatrics, according to the ALARA criteria (As Low As Reasonably Achievable), 2) a CT system is not usable due to movement limitations, such as during surgery or in an ICU and 3) where costs issues complicate the availability of CT systems, such as rural areas or underdeveloped countries. The results of this thesis has a clear application in the industry, since it is part of a proof of concept of the new generation of planar X-ray systems that will be commercialized worldwide by the company SEDECAL (Madrid, Spain).Los últimos años están viendo un rápido avance de los sistemas de radiología hacia el uso de detectores digitales y a una mayor flexibilidad de movimientos de los principales componentes del sistema, el tubo de rayos X y el detector. Esta evolución abre la posibilidad de incorporar capacidades avanzadas en sistemas de imagen plana por rayos X proporcionando nueva información valiosa para el diagnóstico. Dos retos en radiografía son obtener imágenes cuantitativas y reducir la superposición de tejidos debida a la naturaleza proyectiva de la técnica. La radiografía de energía dual, basada en la adquisición de dos imágenes a diferente kilovoltaje, permite obtener imágenes de tejido blando y hueso por separado. Los beneficios de esta técnica que aumenta la cantidad de información sin añadir un tiempo de adquisición o de dosis de radiación extra significativos frente al uso de radiografía convencional, han sido demostrados en diferentes aplicaciones. En otra dirección, un avance realmente disruptivo sería la obtención de imagen 3D con sistemas diseñados únicamente para imagen plana. La incorporación de capacidades tomográficas en estos sistemas tendría que lidiar con la adquisición de un número limitado de proyecciones siguiendo trayectorias no estándar. Esta tesis presenta contribuciones originales en esas dos direcciones: radiografía de energía dual e imagen 3D con sistemas de rayos X diseñados para imagen plana. El trabajo se encuadra en una línea de investigación del grupo de Imagen Biomédica e Instrumentación del Departamento de Bioingeniería e Ingeniería Aerospacial de la Universidad Carlos III de Madrid junto con el Hospital Universitario Gregorio Marañon, centrada en el avance de sistemas de radiología. Esta línea de investigación se desarollada en colaboración con el grupo Computer Architecture, Communications and Systems (ARCOS), de la misma universidad, el grupo Imaging Research Laboratory (IRL) de la Universidad de Washington y el centro de investigación CREATIS, de Francia. Se trata de una línea de investigación con un claro enfoque de transferencia tecnológica a la industria a través de la compañía SEDECAL, una multinacional española de entre las 10 líderes del mundo en el campo de la radiología. La primera contribución de esta tesis es un protocolo completo para incorporar capacidades de energía dual que permitan estudios cuantitativos de imagen plana. La propuesta se basa en una calibración previa con un maniquí simple y de bajo coste formado por dos materiales equivalentes de tejido blando y hueso respectivamente. Comparado con los sistemas actuales DXA (Dual-energy X-ray Absorptiometry), 1) proporciona valores reales de tejido atravesado, 2) se basa en una calibración automática que no requiere la interacción de un técnico con gran experiencia. La segunda contribución es un protocolo nuevo para la incorporación de capacidades tomográficas en sistemas de rayos X originariamente diseñados para imagen plana. Para ello, nos enfrentamos a tres principales dificultades. En primer lugar, las trayectorias que pueden seguir la fuente y el detector en estos sistemas no constituyen órbitas circulares estándares, lo que plantea retos importantes en la caracterización geométrica. Para solventarlo, el protocolo propuesto incluye una calibración geométrica que estima todos los parámetros geométricos del sistema para cada proyección. En segundo lugar, la reconstrucción de un número limitado de proyecciones adquiridas en un rango angular reducido da lugar a artefactos graves cuando se reconstruye con algoritmos convencionales. Para lidiar con estos datos de ángulo limitado, el protocolo incluye un nuevo método avanzado de reconstrucción que incorpora la información de superficie de la muestra, que se puede se obtener con un escáner 3D. Esta información se impone como una restricción siguiendo la formulación de Split Bregman, para compensar la falta de datos. La restricción del espacio de búsqueda a través de la explotación del soporte basado en superficie, es crucial para una recuperación completa del contorno externo de la muestra cuando el rango angular es extremadamente pequeño. El diseño modular, eficiente y flexible de la implementación propuesta permite reconstruir datos de ángulo limitado obtenidos con posiciones de fuente y detector no estándar. En tercer lugar, hasta la fecha, no se ha explorado la optimización del protocolo de adquisición con estos sistemas. Esta tesis incluye un estudio de los protocolos óptimos de adquisición que permitió identificar las posibilidades y limitaciones de estos sistemas de imagen plana. Gracias al método de reconstrucción basado en superficie, es posible reducir el número total de proyecciones hasta el 33% y el rango angular hasta 60 grados. Las contribuciones de esta tesis abren la posibilidad de proporcionar información de profundidad y cuantitativa muy valiosa para la mejora del diagnóstico radiológico. Esto podría impactar considerablemente en la práctica clínica, donde la radiología convencional es todavía la modalidad de imagen más utilizada, abarcando el 80- 90% del total de los exámenes de imagen médica. Estos avances abren la posibilidad de nuevas aplicaciones clínicas en escenarios donde 1) la reducción de la dosis de radiación es clave, como en screening de cáncer de pulmón, de acuerdo con el criterio ALARA (As Low As Reasonably Achievable), 2) no se puede usar un sistema TAC por limitaciones de movimiento como en cirugía o UCI, o 3) el coste limita la disponibilidad de sistemas TAC, como en zonas rurales o en países subdesarrollados. Los resultados de esta tesis presentan una clara aplicación industrial, ya que son parte de un prototipo de la nueva generación de sistemas planos de rayos X que serán distribuidos mundialmente por la compañía SEDECAL.This thesis has been developed as part of several research projects with public funding: - DPI2016-79075-R. ”Nuevos escenarios de tomografía por rayos X”, IP: Mónica Abella García, Ministerio de Economía y Competitividad, 01/01/2017-31/12/2019, 147.620 e. - ”Nuevos escenarios de tomografía por rayos X (NEXT) DPI2016-79075-R. Ministerio de Economía”, Industria y Competitividad. (Universidad Carlos III de Madrid). 30/12/2016-29/12/2019. 147.620 e. (…) - FP7-IMI-2012 (GA-115337), ”PreDict-TB: Model-based preclinical development of anti-tuberculosis drug combinations”. FP7-IMI - Seventh Framework Programme (EC-EFPIA). Unión Europea. (Universidad Carlos III de Madrid). 01/05/2012-31/10/2017. (…) - TEC2013-47270-R, ”Avances en Imagen Radiológica (AIR)”, Ministerio de Economía y Competitividad”, 01/01/2014-31/12/2016. IP: Mónica Abella Garcia and Manuel Desco Menéndez. 160.204 e (…) - RTC-2014-3028-1, ”Nuevos Escenarios Clínicos con Radiología Avanzada (NECRA)”, Ministerio de Economía y Competitividad, 01/06/2014-31/12/2016 IP: Mónica Abella García. 2014-2016. 219.458,96 e - IDI-20130301, ”Nuevo sistema integral de radiografía (INNPROVE: INNovative image PROcessing in medicine and VEterinary)”, IP: Mónica Abella García and Manuel Desco Menéndez. Ministerio de Economía y Competitividad. Subcontratación CDTI, 14/01/2013-31/03/2015. Total: 1.860.629e (UC3M: 325.000e). (Art. 83) - IPT-2012-0401-300000 INNPACTO 2012, ”Tecnologías para Procedimientos Intraoperatorios Seguros y Precisos. XIORT. MINECO. (Universidad Carlos III de Madrid). 01/01/2013-31/12/2015.Programa Oficial de Doctorado en Ingeniería MatemáticaPresidente: Doménec Ros Puig.- Secretario: Cyril Riddell.- Vocal: Yannick Boursie

Quantitative Analysis of Three-Dimensional Cone-Beam Computed Tomography Using Image Quality Phantoms

In the clinical setting, weight-bearing static 2D radiographic imaging and supine 3D radiographic imaging modalities are used to evaluate radiographic changes such as, joint space narrowing, subchondral sclerosis, and osteophyte formation. These respective imaging modalities cannot distinguish between tissues with similar densities (2D imaging), and do not accurately represent functional joint loading (supine 3D imaging). Recent advances in cone-beam CT (CBCT) have allowed for scanner designs that can obtain weight-bearing 3D volumetric scans. The purpose of this thesis was to analyze, design, and implement advanced imaging techniques to quantify image quality parameters of reconstructed image volumes generated by a commercially-available CBCT scanner, and a novel ceiling-mounted CBCT scanner. In addition, imperfections during rotation of the novel ceiling-mounted CBCT scanner were characterized using a 3D printed calibration object with a modification to the single marker bead method, and prospective geometric calibration matrices

Towards breast tomography with synchrotron radiation at Elettra: First images

The aim of the SYRMA-CT collaboration is to set-up the first clinical trial of phase-contrast breast CT with synchrotron radiation (SR). In order to combine high image quality and low delivered dose a number of innovative elements are merged: a CdTe single photon counting detector, state-of-the-art CT reconstruction and phase retrieval algorithms. To facilitate an accurate exam optimization, a Monte Carlo model was developed for dose calculation using GEANT4. In this study, high isotropic spatial resolution (120 μm)3 CT scans of objects with dimensions and attenuation similar to a human breast were acquired, delivering mean glandular doses in the range of those delivered in clinical breast CT (5-25 mGy). Due to the spatial coherence of the SR beam and the long distance between sample and detector, the images contain, not only absorption, but also phase information from the samples. The application of a phase-retrieval procedure increases the contrast-to-noise ratio of the tomographic images, while the contrast remains almost constant. After applying the simultaneous algebraic reconstruction technique to low-dose phase-retrieved data sets (about 5 mGy) with a reduced number of projections, the spatial resolution was found to be equal to filtered back projection utilizing a four fold higher dose, while the contrast-to-noise ratio was reduced by 30%. These first results indicate the feasibility of clinical breast CT with SR

Imaging performance of phase-contrast breast computed tomography with synchrotron radiation and a CdTe photon-counting detector

Purpose: Within the SYRMA-CT collaboration based at the ELETTRA synchrotron radiation (SR) facility (Trieste, Italy) the authors investigated the imaging performance of the phase-contrast computed tomography (CT) system dedicated to monochromatic in vivo 3D imaging of the female breast, for breast cancer diagnosis. Methods: Test objects were imaged at 38 keV using monochromatic SR and a high-resolution CdTe photon-counting detector. Signal and noise performance were evaluated using modulation transfer function (MTF) and Noise Power Spectrum (NPS). Phase-contrast CT images as well as images obtained after the application of a phase-retrieval algorithm were evaluated. The contrast to noise ratio (CNR) and the capability of detecting test microcalcification clusters and soft masses were explored. Results: For a voxel size of (60 \u3bcm)3, phase-contrast images showed higher spatial resolution (6.7 mm-1 at 10% MTF) than corresponding phase retrieval images (2.5 mm-1). Phase retrieval produced a reduction of the noise level as well as an increase of the CNR by more than one order of magnitude, compared to raw phase-contrast images. CaCO3 microcalcifications with a diameter down to 130 \u3bcm were detected both in phase-contrast and in phase retrieval images of the test object. Conclusions: The investigation on test objects indicates that breast CT with a monochromatic SR source is technically feasible in terms of spatial resolution, image noise and contrast, for in vivo 3D imaging with a dose comparable to that of two-view mammography. Phase-retrieved images showed the best performance in the trade-off between spatial resolution and image noise

The realisation of fast X-ray computed tomography using a limited number of projection images for dimensional metrology

Due to the merit of establishing volumetric data, X-ray computed tomography (XCT) is increasingly used as a non-destructive evaluation technique in the quality control of advanced manufactured parts with complex non-line-of-sight features. However, the cost of measurement time and data storage hampers the adoption of the technique in production lines. Commercial fast XCT utilises X-ray detectors with fast detection capability, which can be expensive and result in a large amount of data. This paper discussed a different approach, where fast XCT was realised via the acquisition of a small number of projection images instead of full projection images. An established total variation (TV) algorithm was used to handle the reconstruction. The paper investigates the feasibility of using the TV algorithm in handling a significantly reduced number of projection images for reconstruction. This allows a reduction of measurement time from 52 min to 1 min for a typical industrial XCT system with an exposure time of 1 s per projection. It also enables a reduction of data size proportionally. A test strategy including both quantitative and qualitative test metrics was considered to evaluate the effectiveness of the reconstruction algorithm. The qualitative evaluation includes both the signal to noise ratio and the contrast to noise ratio. The quantitative evaluation was established using reference samples with different internal and external geometries. Simulation data were used in the assessment considering various influence factors, such as X-ray source property and instrument noise. The results demonstrated the possibility of using advanced reconstruction algorithms in handling XCT measurements with a significantly limited number of projection images for dimensional measurements. This work lays down the foundation for conducting fast XCT measurements without the need for instrument alteration or enhancement. Although the reconstruction time required is still considerable, various possibilities to improve this have been discussed.</p