Cone-Beam Computed Tomography for Oral and Maxillofacial Imaging

The invention of computed tomography (CT) technique revolutionized diagnostic imaging. Compared to conventional X-ray imaging procedures, CT involves higher radiation doses. Recently, cone-beam CT (CBCT) specifically designed for maxillofacial imaging was introduced. CBCT technique is based on a cone-shaped X-ray beam centered on a two-dimensional (2D) detector. The detector system performs one rotation around the patient, producing a series of 2D images which are then reconstructed in a 3D data set. The contemporary knowledge regarding CBCT and its proper application guides the practitioner for improvement in diagnostic purposes and treatment planning. The aim of this chapter is to focus on the details, advantages, drawbacks, and clinical applications of CBCT as a headmost CT imaging technique in the oral and maxillofacial (OMF) region. The main clinical applications of CBCT in the OMF region are dentistry including dentoalveolar and maxillofacial surgery, orthodontics, endodontics, and periodontics; and otolaryngology. The aforementioned clinical use of CBCT was described in detail with illustrated sample cases. In most of the cases in OMF region, CBCT takes the place of multi-slice CT. Thus, clinicians should know the clinical applications and capabilities of CBCT technique with its drawbacks

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

A 3D computer assisted Orthopedic Surgery Planning approach based on planar radiography

Dissertação de mestrado integrado em Engenharia Biomédica (área de especialização em Informática Médica)The main goal of this work consisted in develop a system to perform the 3D reconstruction of bone models from radiographic images. This system can be then integrated with a commercial software that performs pre-operative planning of orthopedic surgeries. The benefit of performing this 3D reconstruction from planar radiography is that this modality has some advantages over other modalities that perform this reconstruction directly, like CT and MRI. To develop the system it was used radiographic images of the femur obtained from medical image databases online. It was also used a generic model of the femur available in the online repository BEL. This generic model completes the information missing in the radiographic images. It was developed two methods to perform the 3D reconstruction through the deformation of the generic model, one uses triangulation of extracted edge points and the other don't. The first method was not successful, the final model had very low thickness, possibly because the triangulation process was not performed correctly. With the second method it was obtained a 3D bone model of the femur aligned with the radiographic images of the patient and with the same size as the patient's bone. However, the obtained model still needs some adjustment to coincide fully with reality. To perform this is necessary to enhance the deformation step of the model so that it will have the same shape as the patient's bone. The second method is more advantageous because it doesn't need the parameters of the x-ray imaging system. However, it's necessary to enhance the step deformation of this method so that the final model matches patient's anatomy.O principal objetivo deste trabalho consistiu em desenvolver um sistema capaz de realizar a reconstrução 3D de modelos ósseos a partir de imagens radiográficas. Este sistema pode posteriormente ser integrado num produto comercial que realiza o planeamento pré-operativo de cirurgias ortopédicas. O benefício de realizar esta reconstrução 3D a partir de radiografias está relacionado com o facto desta modalidade ter vantagens em relação às outras modalidades que fazem esta reconstrução diretamente, como as modalidades CT e MRI. Para desenvolver este sistema foram usadas imagens radiográficas do fémur obtidas através de bases de dados online de imagens médicas. Também foi usado um modelo genérico do fémur disponível no repositório online BEL. Este modelo genérico completa a informação que está em falta nas imagens radiográficas. Foram desenvolvidos dois métodos, que realizam a reconstrução 3D através da deformação do modelo genérico sendo que num é feita a triangulação de pontos dos contornos e noutro não. O primeiro método não foi bem sucedido, visto que o modelo final tinha uma espessura muito pequena, possivelmente devido ao facto do processo de triangulação não ter sido executado corretamente. Com o segundo método foi obtido um modelo 3D do fémur alinhado com as imagens radiográficas do paciente e com o mesmo tamanho do osso do paciente. No entanto, o modelo obtido carece ainda de alguma afinação de modo a coincidir na íntegra com a realidade. Para fazer isto é necessário melhorar o passo de deformação do modelo, para que este fique com a mesma forma do osso do paciente. O segundo método é mais vantajoso porque não necessita dos parâmetros dos sistema de raios- X. No entanto, é necessário melhorar o passo de deformação deste método para que o modelo final coincida com a anatomia do paciente

Interlandmark measurements from lodox statscan images with application to femoral neck anteversion assessment

Includes abstract.Includes bibliographical references.Clinicians often take measurements between anatomical landmarks on X-ray radiographs for diagnosis and treatment planning, for example in orthopaedics and orthodontics. X-ray images, however, overlap three-dimensional internal structures onto a two-dimensional plane during image formation. Depth information is therefore lost and measurements do not truly reflect spatial relationships. The main aim of this study was to develop an inter-landmark measurement tool for the Lodox Statscan digital radiography system. X-ray stereophotogrammetry was applied to Statscan images to enable three-dimensional point localization for inter-landmark measurement using two-dimensional radiographs. This technique requires images of the anatomical region of interest to be acquired from different perspectives as well as a suitable calibration tool to map image coordinates to real world coordinates. The Statscan is suited to the technique because it is capable of axial rotations for multiview imaging. Three-dimensional coordinate reconstruction and inter-landmark measurements were taken using a planar object and a dry pelvis specimen in order to assess the intra-observer measurement accuracy, reliability and precision. The system yielded average (X, Y, Z) coordinate reconstruction accuracy of (0.08 0.12 0.34) mm and resultant coordinate reconstruction accuracy within 0.4mm (range 0.3mm – 0.6mm). Inter-landmark measurements within 2mm for lengths and 1.80 for angles were obtained, with average accuracies of 0.4mm (range 0.0mm – 2.0 mm) and 0.30 (range 0.0 – 1.8)0 respectively. The results also showed excellent overall precision of (0.5mm, 0.10) and were highly reliable when all landmarks were completely visible in both images. Femoral neck anteversion measurement on Statscan images was also explored using 30 dry right adult femurs. This was done in order to assess the feasibility of the algorithm for a clinical application. For this investigation, four methods were tested to determine the optimal landmarks for measurement and the measurement process involved calculation of virtual landmarks. The method that yielded the best results produced all measurements within 10 of reference values and the measurements were highly reliable with very good precision within 0.10. The average accuracy was within 0.40 (range 0.10 –0.80).In conclusion, X-ray stereophotogrammetry enables accurate, reliable and precise inter-landmark measurements for the Lodox Statscan X-ray imaging system. The machine may therefore be used as an inter-landmark measurement tool for routine clinical applications

Three-Dimensional Biplanar Reconstruction of the Scoliotic Spine for Standard Clinical Setup

Tese de Doutoramento. Engenharia Informática. Faculdade de Engenharia. Universidade do Porto. 201

3D reconstruction of ribcage geometry from biplanar radiographs using a statistical parametric model approach

Rib cage 3D reconstruction is an important prerequisite for thoracic spine modelling, particularly for studies of the deformed thorax in adolescent idiopathic scoliosis. This study proposes a new method for rib cage 3D reconstruction from biplanar radiographs, using a statistical parametric model approach. Simplified parametric models were defined at the hierarchical levels of rib cage surface, rib midline and rib surface, and applied on a database of 86 trunks. The resulting parameter database served to statistical models learning which were used to quickly provide a first estimate of the reconstruction from identifications on both radiographs. This solution was then refined by manual adjustments in order to improve the matching between model and image. Accuracy was assessed by comparison with 29 rib cages from CT scans in terms of geometrical parameter differences and in terms of line-to-line error distance between the rib midlines. Intra and inter-observer reproducibility were determined regarding 20 scoliotic patients. The first estimate (mean reconstruction time of 2’30) was sufficient to extract the main rib cage global parameters with a 95% confidence interval lower than 7%, 8%, 2% and 4° for rib cage volume, antero-posterior and lateral maximal diameters and maximal rib hump, respectively. The mean error distance was 5.4 mm (max 35mm) down to 3.6 mm (max 24 mm) after the manual adjustment step (+3’30). The proposed method will improve developments of rib cage finite element modeling and evaluation of clinical outcomes.This work was funded by Paris Tech BiomecAM chair on subject specific muscular skeletal modeling, and we express our acknowledgments to the chair founders: Cotrel foundation, Société générale, Protéor Company and COVEA consortium. We extend your acknowledgements to Alina Badina for medical imaging data, Alexandre Journé for his advices, and Thomas Joubert for his technical support

Comparative linear accuracy and reliability of cone beam CT derived 2-dimensional and 3-dimensional images constructed using an orthodontic volumetric rendering program.

The purpose of this project was to compare the accuracy and reliability of linear measurements made on 2D projections and 3D reconstructions using Dolphin 3D software (Chatsworth, CA) as compared to direct measurements made on human skulls. The linear dimensions between 6 bilateral and 8 mid-sagittal anatomical landmarks on 23 dentate dry human skulls were measured three times by multiple observers using a digital caliper to provide twenty orthodontic linear measurements. The skulls were stabilized and imaged via PSP digital cephalometry as well as CBCT. The PSP cephalograms were imported into Dolphin (Chatsworth, CA, USA) and the 3D volumetric data set was imported into Dolphin 3D (Version 2.3, Chatsworth, CA, USA). Using Dolphin 3D, planar cephalograms as well as 3D volumetric surface reconstructions were (3D CBCT) generated. The linear measurements between landmarks of each three modalities were then computed by a single observer three times. For 2D measurements, a one way ANOVA for each measurement dimension was calculated as well as a post hoc Scheffe multiple comparison test with the anatomic distance as the control group. 3D measurements were compared to anatomic truth using Student\u27s t test (PiÜ50.05). The intraclass correlation coefficient (ICC) and absolute linear and percentage error were determined as indices of intraobserver reliability. Our results show that for 2D mid sagittal measurements that Simulated LC images are accurate and similar to those from PSP images (except for Ba-Na), and for bilateral measurements simulated LC measurements were similar to PSP but less accurate, underestimating dimensions by between 4.7% to 17%.For 3D volumetric renderings, 2/3 rd of CBCT measurements are statistically different from actual measurements, however this possibly is not clinically relevant

3D approximation of scapula bone shape from 2D X-ray images using landmark-constrained statistical shape model fitting

Two-dimensional X-ray imaging is the dominant imaging modality in low-resource countries despite the existence of three-dimensional (3D) imaging modalities. This is because fewer hospitals in low-resource countries can afford the 3D imaging systems as their acquisition and operation costs are higher. However, 3D images are desirable in a range of clinical applications, for example surgical planning. The aim of this research was to develop a tool for 3D approximation of scapula bone from 2D X-ray images using landmark-constrained statistical shape model fitting. First, X-ray stereophotogrammetry was used to reconstruct the 3D coordinates of points located on 2D X-ray images of the scapula, acquired from two perspectives. A suitable calibration frame was used to map the image coordinates to their corresponding 3D realworld coordinates. The 3D point localization yielded average errors of (0.14, 0.07, 0.04) mm in the X, Y and Z coordinates respectively, and an absolute reconstruction error of 0.19 mm. The second phase assessed the reproducibility of the scapula landmarks reported by Ohl et al. (2010) and Borotikar et al. (2015). Only three (the inferior angle, acromion and the coracoid process) of the eight reproducible landmarks considered were selected as these were identifiable from the two different perspectives required for X-ray stereophotogrammetry in this project. For the last phase, an approximation of a scapula was produced with the aid of a statistical shape model (SSM) built from a training dataset of 84 CT scapulae. This involved constraining an SSM to the 3D reconstructed coordinates of the selected reproducible landmarks from 2D X-ray images. Comparison of the approximate model with a CT-derived ground truth 3D segmented volume resulted in surface-to-surface average distances of 4.28 mm and 3.20 mm, using three and sixteen landmarks respectively. Hence, increasing the number of landmarks produces a posterior model that makes better predictions of patientspecific reconstructions. An average Euclidean distance of 1.35 mm was obtained between the three selected landmarks on the approximation and the corresponding landmarks on the CT image. Conversely, a Euclidean distance of 5.99 mm was obtained between the three selected landmarks on the original SSM and corresponding landmarks on the CT image. The Euclidean distances confirm that a posterior model moves closer to the CT image, hence it reduces the search space for a more exact patient-specific 3D reconstruction by other fitting algorithms

Comparative linear accuracy of cone beam CT derived 3D images in orthodontic analysis.

Objective . To compare the in vitro reliability and accuracy of linear measurements between cephalometric landmarks on CBCT 3D images with varying basis projection images to direct measurements on human skulls. Methods . Sixteen linear dimensions between anatomical sites marked on 19 human skulls were directly measured. Skulls were imaged with CBCT at three settings: 153, 306, and 612 basis projections. The mean absolute error and modality mean of linear measurements between landmarks on 3D images were compared to the anatomic truth. Results . No difference in mean absolute error between the scan settings was found. The average skull absolute error between marked reference points were less than the distances between unmarked reference sites. Conclusion . CBCT measurements were consistent between scan sequences and for direct measurements between marked reference points. Reducing the number of projections for 3D reconstruction did not lead to reduced dimensional accuracy and potentially provides reduced patient radiation exposure