Simulation de radiographies à partir d'images tomodensitométriques pour l'enseignement de l'anatomie radiographique en médecine vétérinaire

L'un des principaux défis de l'interprétation radiographique réside dans la compréhension de l’anatomie radiographique, laquelle est intrinsèquement liée à la disposition tridimensionnelle des structures anatomiques et à l’impact du positionnement du tube radiogène vis-à-vis de ces structures lors de l'acquisition de l'image. Traditionnellement, des radiographies obtenues selon des projections standard sont employées pour enseigner l'anatomie radiographique en médecine vétérinaire. La tomodensitométrie − ou communément appelée CT (Computed Tomography) − partage plusieurs des caractéristiques de la radiographie en ce qui a trait à la génération des images. À l’aide d'un plug-in spécialement développé (ORS Visual ©), la matrice contenant les images CT est déformée pour reproduire les effets géométriques propres au positionnement du tube et du détecteur vis-à-vis du patient radiographié, tout particulièrement les effets de magnification et de distorsion. Afin d'évaluer le rendu des images simulées, différentes régions corporelles ont été imagées au CT chez deux chiens, un chat et un cheval, avant d'être radiographiées suivant des protocoles d'examens standards. Pour valider le potentiel éducatif des simulations, dix radiologistes certifiés ont comparé à l'aveugle neuf séries d'images radiographiques simulées aux séries radiographiques standard. Plusieurs critères ont été évalués, soient le grade de visualisation des marqueurs anatomiques, le réalisme et la qualité radiographique des images, le positionnement du patient et le potentiel éducatif de celles-ci pour différents niveaux de formation vétérinaire. Les résultats généraux indiquent que les images radiographiques simulées à partir de ce modèle sont suffisamment représentatives de la réalité pour être employées dans l’enseignement de l’anatomie radiographique en médecine vétérinaire.Understanding radiographic anatomy, which is intimately linked to the comprehension of tridimensional anatomy and the impact of patient, radiographic tube and x-ray detector positioning, represents a challenge for students. Traditionally, radiographs obtained under specific angles of projection have been used for teaching radiographic anatomy. Computed tomography (CT) shares several features with radiography with regard to image production. A plug-in was developed for a DICOM viewer (ORS visual ©) simulating radiographs using CT datasets. This plug-in distorts the CT image matrix to reproduce the magnification and distortion effects that take place in radiographs due to the variations in radiographic tube, patient and detector positioning and angulation. In order to test this model, specific body parts of two dogs, one cat and one horse were radiographed and CT-scanned. The CT datasets were used to generate a total of nine series of radiographic simulations that could be compared to corresponding standard radiographic projections. Ten board-certified veterinary radiologists blindly scored several parameters in these image series, including the visualization of specific anatomical landmarks, image realism and quality, patient positioning, and the educational potential for students and veterinarians of variable degree of veterinary training Overall results indicate that simulated radiographs are representative enough to be used to teach several concepts of image formation and radiographic anatomy in veterinary radiology

Biomechanical importance of proximal human femur morphology and mechanics in orthopaedic purposes

Bone morphology is essential in orthopedic surgery to perform precise preoperative planning and surgery as well as to appropriately design optimal medical implants. In this study we provided a database of surgically important morphological parameters of proximal human femur for orthopedic and biomedical research purposes (study 1), indicated accuracy of the 3D reconstructed images in comparison with the optical 3D scan of real human femur (study 2), and reported the accuracy and reliability of the developed image-based finite element model in comparison with the experimental results (study 3)

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

Statistical Image Reconstruction and Motion Estimation for Image-Guided Radiotherapy.

Image reconstruction and motion estimation are very important for image-guided radiotherapy (IGRT). Three-dimensional reconstruction of patient anatomy using X-ray computed tomography (CT) allows identification of the location of a tumor prior to treatment. The locations of tumorsmay change during actual treatment due to movement such as respiratory motion. Motion estimation helps optimize the accuracy and precision of radiotherapy so that more of the normal surrounding tissue can be spared. This dissertation addresses several important issues related to these two core components of IGRT. Firstly, we developed two new separable footprint (SF) projector methods for X-ray conebeam CT. The SF projectors approximate the voxel footprint functions as 2D separable functions. The SF-TR projector uses trapezoid functions in the transaxial direction and rectangular functions in the axial direction, whereas the SF-TT projector uses trapezoid functions in both directions. Both SF projector methods are more accurate than the distance-driven (DD) projector, which is a current state-of-the-art method in the field. The SF-TT projector is more accurate than the SF-TR projector for rays associated with large cone angles. In addition, the SF-TR projector has similar computation speed with the DD projector and the SF-TT projector is about two times slower. Secondly, we proposed a statistical penalized weighted least-squares (PWLS) method with edge-preserving regularization to reconstruct two basis materials from a single-energy CT scan acquired with differential filtration, such as a split filter or a bow-tie filter. It requires only the use of suitable filters between the X-ray tube and the patient. For both filtration methods, the proposed PWLS method reconstructed soft tissue and bone images with lower RMS errors, reduced the beam-hardening artifacts much more effectively and produced lower noise, as compared with the traditional non-iterative Joseph and Spital method. Thirdly, we conducted an objective characterization of the influence of rotational arc length on accuracy of motion estimation for projection-to-volume targeting during rotational therapy. Simulations illustrate the potential accuracy of limited-angle projection-to-volume alignment. Registration accuracy can be sensitive to angular center, tends to be lower along direction of the projection set, and tends to decrease away from the rotation center.Ph.D.Electrical Engineering: SystemsUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/86254/1/yonglong_1.pd

Spline projection-based volume-to-image registration

This thesis focuses on the rigid-body registration of a three-dimensional model of an object to a set of its two-dimensional projections. The main contribution is the development of two registration algorithms that use a continuous model of the volume based on splines, either in the space domain or in the frequency domain. This allows for a well-defined gradient of the dissimilarity measure, which is a necessary condition for efficient and accurate registration. The first part of the thesis contains a review of the literature on volume-to- image registration. Then, we discuss data interpolation in the space domain and in the frequency domain. The basic concepts of our registration strategy are given in the second part of the thesis. We present a novel one-step approach for fast ray casting to simulate space-based volume projections. We also discuss the use of the central-slice theorem to simulate frequency-based volume projections. Then, we consider the question of the registration robustness. To improve the robustness of the space-based approach, we apply a multiresolution optimization strategy where spline-based data pyramids are processed in coarse-to-fine fashion, which improves speed as well. To improve the robustness of the frequency-based registration, we apply a coarse-to-fine strategy that involves weights in the frequency domain. In the third part, we apply our space-based algorithm to computer-assisted orthopedic surgery while adapting it to the perspective projection model. We show that the registration accuracy achieved using the orthopedic data is consistent with the current standards. Then, we apply our frequency-based registration to three-dimensional electron-microscopy application. We show that our algorithm can be used to obtain a refined solution with respect to currently available algorithms. The novelty of our approach is in dealing with a continuous space of geometric parameters, contrary to the standard methods which deal with quantized parameters. We conclude that our continuous parameter space leads to better registration accuracy. Last, we compare the performance of the frequency-based algorithm with that of the space-based algorithm in the context of electron microscopy. With these data, we observe that frequency-based registration algorithm outperforms the space-based one, which we attribute to the suitability of interpolation in the frequency domain when dealing with strictly space-limited data