An external field prior for the hidden Potts model with application to cone-beam computed tomography

In images with low contrast-to-noise ratio (CNR), the information gain from the observed pixel values can be insufficient to distinguish foreground objects. A Bayesian approach to this problem is to incorporate prior information about the objects into a statistical model. A method for representing spatial prior information as an external field in a hidden Potts model is introduced. This prior distribution over the latent pixel labels is a mixture of Gaussian fields, centred on the positions of the objects at a previous point in time. It is particularly applicable in longitudinal imaging studies, where the manual segmentation of one image can be used as a prior for automatic segmentation of subsequent images. The method is demonstrated by application to cone-beam computed tomography (CT), an imaging modality that exhibits distortions in pixel values due to X-ray scatter. The external field prior results in a substantial improvement in segmentation accuracy, reducing the mean pixel misclassification rate for an electron density phantom from 87% to 6%. The method is also applied to radiotherapy patient data, demonstrating how to derive the external field prior in a clinical context.

Development of a single photon counting computed tomography system using MPGDs

The development of computed tomography systems with energy resolving detectors is a current challenge in medical physics and biomedical engineering. A computed tomography system of this kind allows getting complementary informations relatively to conventional systems, that can help the medical diagnosis, being of great interest in medicine. The work described in this thesis is related to the development of a computed tomography system using micropattern gaseous detectors, which allow storing, simultaneously, information about the interaction position and the energy of each single photon that interacts with the detector. This kind of detectors has other advantages concerning the cost and characteristics of operation when compared with solid state detectors. Tomographic acquisitions were performed using a MicroHole & Strip Plate based detector, which allowed reconstructing cross-sectional images using energy windows, applying the energy weighting technique and performing multi-slice and tri-dimensional reconstructions. The contrast-to-noise ratio was improved by 31% by applying the energy weighting technique, comparing with the corresponding image obtained with the current medical systems. A prototype of a computed tomography with flexibility to change the detector was developed, making it possible to apply different detectors based on Thick-COBRA. Several images acquired with these detectors are presented and demonstrate their applicability in X-ray imaging. When operating in NeCH4, the detector allowed a charge gain of 8 104, an energy resolution of 20% (full width at half maximum at 8 keV), a count rate of 1 106 Hz/mm2, a very stable operation (gain fluctuations below 5%) and a spacial resolution of 1.2 mm for an energy photon of 3.6 keV. Operating the detector in pure Kr allowed increasing the detection efficiency and achieving a charge gain of 2 104, an energy resolution of 32% (full width at half maximum at 22 keV), a count rate of 1 105 Hz/mm2, very stable operation and a spatial resolution of 500 m. The software already existing in the group was improved and tools to correct geometric misalignments of the system were also developed. The reconstructions obtained after geometrical correction are free of artefacts due to the referred misalignments.

Automatic X-ray screening of aircraft hold luggage

R. Benjamin (1995) addressed the application of the “object 3D” X ray reconstruction technique for electronically “unpacking” suspect items, when screening aircraft luggage. However, there is no satisfactory solution to the mass screening of hold luggage. Computed Tomography, CT, entails excessive radiation dosages, and its rate of throughput is quite inadequate. A novel variant of “object 3D” is therefore put forward, adapting some of the technology of existing cabin luggage screening systems-but on a substantially larger scale-which does achieve the required throughput at an acceptable radiation dosage and cost.

Evaluating craniofacial asymmetry with digital cephalometric images and cone-beam computed tomography

Introduction: The aim of this study was to evaluate craniofacial asymmetry by using 2-dimensional (2D) poster-oanterior cephalometric images, 3-dimensional cone-beam computed tomography (CBCT), and physical measurements (gold standard). Methods: Ten dry human skulls were assessed, and radiopaque markers were placed on 17 skeletal landmarks. Twenty linear measurements were taken on each side to compare the right and left sides and to compare these measurements with the physical measurements made with a digital caliper. To acquire the 2D posteroanterior radiographs, an Extraoral Phosphor Storage Plate (Air Techniques, Chicago, Ill) was used as the image receptor with a Eureka x-ray-Duocon Machlett unit (Machlett Laboratores, Chicago, Ill). Three-dimensional imaging data were acquired from a CB MercuRay (Hitachi Medical, Tokyo, Japan). Results: on average, the right side was larger than the left for most of the 20 distances evaluated in the digital 2D and the CBCT images, and there was poor agreement between the digital 2D images and the physical measurements (kappa = 0.0609) and almost perfect agreement (kappa = 0.92) between the CBCT and physical measurements when individual measurements were considered. Conclusions: Human skulls, with no apparent asymmetry, had some differences between the right and left sides, with dominance for the right side but with no clinical significance. CBCT can better evaluate craniofacial morphology when compared with digital 2D images. (Am J Orthod Dentofacial Orthop 2011; 139: e523-e531)

Influence of anatomical location on CT numbers in cone beam computed tomography

Objective: To assess the influence of anatomical location on computed tomography (CT) numbers in mid- and full field of view (FOV) cone beam computed tomography (CBCT) scans. Study Design: Polypropylene tubes with varying concentrations of dipotassium hydrogen phosphate (K2HPO4) solutions (50-1200 mg/mL) were imaged within the incisor, premolar, and molar dental sockets of a human skull phantom. CBCT scans were acquired using the NewTom 3G and NewTom 5G units. The CT numbers of the K2HPO 4 phantoms were measured, and the relationship between CT numbers and K2HPO4 concentration was examined. The measured CT numbers of the K2HPO4 phantoms were compared between anatomical sites. Results: At all six anatomical locations, there was a strong linear relationship between CT numbers and K2HPO4 concentration (R 2 > 0.93). However, the absolute CT numbers varied considerably with the anatomical location. Conclusion: The relationship between CT numbers and object density is not uniform through the dental arch on CBCT scans. © 2013 Elsevier Inc.

Assessment of CT numbers in limited and medium field-of-view scans taken using Accuitomo 170 and Veraviewepocs 3De cone-beam computed tomography scanners

To assess the influence of anatomic location on the relationship between computed tomography (CT) number and X-ray attenuation in limited and medium field-of-view (FOV) scans. Materials and Methods Tubes containing solutions with different concentrations of K2HPO4 were placed in the tooth sockets of a human head phantom. Cone-beam computed tomography (CBCT) scans were acquired, and CT numbers of the K2HPO4 solutions were measured. The relationship between CT number and K2HPO4 concentration was examined by linear regression analyses. Then, the variation in CT number according to anatomic location was examined. Results The relationship between K2HPO4 concentration and CT number was strongly linear. The slopes of the linear regressions for the limited FOVs were almost 2-fold lower than those for the medium FOVs. The absolute CT number differed between imaging protocols and anatomic locations. Conclusion There is a strong linear relationship between X-ray attenuation and CT number. The specific imaging protocol and anatomic location of the object strongly influence this relationship.

Mandibular cortical bone evaluation on cone beam computed tomography images of patients with bisphosphonate related osteonecrosis of the jaws

Objectives The objective of this study was to develop a technique for detecting cortical bone dimensional changes in patients with bisphosphonate-related osteonecrosis of the jaw (BRONJ). Study Design Subjects with BRONJ who had cone-beam computed tomography imaging were selected, with age- and gender-matched controls. Mandibular cortical bone measurements to detect bisphosphonate-related cortical bone changes were made inferior to mental foramen, in 3 different ways: within a fixed sized rectangle, in a rectangle varying with the cortical height, and a ratio between area and height. Results Twelve BRONJ cases and 66 controls were evaluated. The cortical bone measurements were significantly higher in cases than controls for all 3 techniques. The bone measurements were strongly associated with BRONJ case status (odds ratio 3.36-7.84). The inter-rater reliability coefficients were high for all techniques (0.71-0.90). Conclusions Mandibular cortical bone measurement is a potentially useful tool in the detection of bone dimensional changes caused by bisphosphonates. Long-term administration of bisphosphonates (BPs) affects bone quality and metabolism following accumulation in bone.1 Since the first cases of bisphosphonate-related osteonecrosis of the jaw (BRONJ) were published in 2003,2 there has been a search for factors that can predict the onset of the condition. Oral and intravenous BPs reduce bone resorption, increase mineral content of bone, and alter bony architecture.3, 4, 5 and 6 Previous studies have demonstrated these changes both radiographically and following histologic analysis.1, 3, 7, 8, 9 and 10 The BP-related jaw changes may present radiological features, such as thickening of lamina dura and cortical borders, diffuse sclerosis, and narrowing of the mandibular canal3 and 11; however, oral radiographs of patients taking BPs do not consistently show radiographic changes to the jaws.11 and 12 The challenge is to find imaging tools that could improve the detection of changes in the bone associated with BP use. Various skeletal radiographic features associated with BRONJ in conventional periapical and panoramic radiographs, computed tomography, magnetic resonance imaging, and nuclear bone scanning have been described.3, 8, 9, 10 and 11 There has also been a search for BP-related quantitative methods for the evaluation of radiographic images, to avoid observer subjectivity in interpretation. Factors thought to be important include trabecular and cortical structure, and bone mineralization.4 Consequently, measurable bone data have been reported in subjects taking BPs through many techniques, including bone density, architecture, and cortical bone thickness.1, 4, 7 and 13 Trabecular microarchitecture of postmenopausal women has been evaluated with noninvasive techniques, such as high-resolution magnetic resonance images showing less deterioration of the bone 1 year after initiation of oral BP therapy.4 A decrease in bone turnover and a trend for an increase in the bone wall thickness has been detected by histomorphometry in subjects taking BPs.1 Alterations in the cortical structure of the second metacarpal have been detected in digital x-ray radiogrammetry of postmenopausal women treated with BPs.7 Mandibular cortical width may be measured on dental panoramic radiographs, and it has been suggested as a screening tool for referring patients for bone densitometry for osteoporosis investigation.14 and 15 Inhibition of the intracortical bone remodeling in the mandible of mice taking BPs has been reported.16 Thus, imaging evaluation of the mandibular cortical bone could be a biologically plausible way to detect BP bone alterations. Computed tomography can assess both cortical and trabecular bone characteristics. Cone-beam computed tomography (CBCT) can provide 3-dimensional information, while using lower doses and costing less than conventional CT. The CBCT images have been studied as a tool for the measurement of trabecular bone in patients with BRONJ.13 Therefore, cortical bone measurements on CBCT of the jaws might also help to understand bone changes in patients with BRONJ. There is no standard in quantifying dimensional changes of mandibular cortical bone. We explored several different approaches to take into consideration possible changes in length, area, and volume. These led to the 3 techniques developed in this study. This article reports a matched case-control study in which mandibular cortical bone was measured on CBCT images of subjects with BRONJ and controls. The aim of the study was to explore the usefulness of 3 techniques for detecting mandibular cortical bone dimensional changes caused by BP.

Assessing apical transportation in curved canals: comparison between cross-sections and micro-computed tomography

The aim of this study was to compare two methods of assessing apical transportation in curved canals after rotary instrumentation, namely, cross-sections and micro-computed tomography (mu CT). Thirty mandibular molars were divided into two groups and prepared according to the requirements of each method. In G1 (cross-sections), teeth were embedded in resin blocks and sectioned at 2.0, 3.5, and 5.0 mm from the anatomic apex. Pre- and postoperative sections were photographed and analyzed. In G2 (mu CT), teeth were embedded in a rubber-base impression material and scanned before and after instrumentation. Mesiobuccal canals were instrumented with the Twisted File (TF) system (SybronEndo, Orange, USA), and mesiolingual canals, with the Endo Sequence (ES) system (Brasseler, Savannah, USA). Images were reconstructed, and sections corresponding to distances 2.0, 3.5, and 5.0 mm from the anatomic apex were selected for comparison. Data were analyzed using Mann-Whitney's test at a 5% significance level. The TF and ES instruments produced little deviation from the root canal center, with no statistical difference between them (P > 0.05). The canal transportation results were significantly lower (0.056 mm) in G2 than in G1 (0.089 mm) (p = 0.0012). The mu CT method was superior to the cross-section method, especially in view of its ability to preserve specimens and provide results that are more closely related to clinical situations.

"Drug mules" as a radiological challenge: Sensitivity and specificity in identifying internal cocaine in body packers, body pushers and body stuffers by computed tomography, plain radiography and Lodox

PURPOSE: The purpose of our study was to retrospectively evaluate the specificity, sensitivity and accuracy of computed tomography (CT), digital radiography (DR) and low-dose linear slit digital radiography (LSDR, Lodox(®)) in the detection of internal cocaine containers. METHODS: Institutional review board approval was obtained. The study collectively consisted of 83 patients (76 males, 7 females, 16-45 years) suspected of having incorporated cocaine drug containers. All underwent radiological imaging; a total of 135 exams were performed: nCT=35, nDR=70, nLSDR=30. An overall calculation of all "drug mules" and a specific evaluation of body packers, pushers and stuffers were performed. The gold standard was stool examination in a dedicated holding cell equipped with a drug toilet. RESULTS: There were 54 drug mules identified in this study. CT of all drug carriers showed the highest diagnostic accuracy 97.1%, sensitivity 100% and specificity 94.1%. DR in all cases was 71.4% accurate, 58.3% sensitive and 85.3% specific. LSDR of all patients with internal cocaine was 60% accurate, 57.9% sensitive and 63.4% specific. CONCLUSIONS: CT was the most accurate test studied. Therefore, the detection of internal cocaine drug packs should be performed by CT, rather than by conventional X-ray, in order to apply the most sensitive exam in the medico-legal investigation of suspected drug carriers. Nevertheless, the higher radiation applied by CT than by DR or LSDR needs to be considered. Future studies should include evaluation of low dose CT protocols in order to address germane issues and to reduce dosage.

The presence of iodinated contrast agents amplifies DNA radiation damage in computed tomography

The purpose of this study was to determine the influence of iodinated contrast agents on the formation of DNA double-strand breaks in vitro in lymphocytes and to verify these results in patients undergoing diagnostic computed tomography examinations. Blood samples were irradiated in vitro in the presence of iodinated X-ray contrast agent. Controls were irradiated without contrast agent. Fourteen patients were investigated using contrast-enhanced computed tomography (CT), and 14 other patients with unenhanced CT. Blood samples were taken prior to and 5 min and 1, 2 and 24 h after the CT examination. In these blood samples the average number of γH2Ax-foci per lymphocyte was enumerated by fluorescence microscopy. Statistical differences between foci numbers developed in the presence and absence of contrast agent were tested using an independent sample t-test. In vitro foci numbers after irradiation were significantly higher when contrast agent was present during irradiation. In vivo, γH2Ax-foci levels were 58% higher in patients undergoing contrast-enhanced CT compared with those undergoing unenhanced CT. In the presence of iodinated contrast agents DNA, damage is increased and the radiation dose is not the only factor affecting the amount of DNA damage. Individual patient characteristics and biological dosimetry applications, e.g. the analysis of γH2Ax-foci, have to be considered.

Maximum intensity projection of cranial computed tomography data for dental identification

Dental radiographs play the major role in the identification of victims in mass casualties besides DNA. Under circumstances such as those caused by the recent tsunami in Asia, it is nearly impossible to document the entire dentition using conventional x-rays as it would be too time consuming. Multislice computed tomography can be used to scan the dentition of a deceased within minutes, and the postprocessing software allows visualization of the data adapted to every possible antemortem x-ray for identification. We introduce the maximum intensity projection of cranial computed tomography data for the purpose of dental identification exemplarily in a case of a burned corpse. As transportable CT scanners already exist, these could be used to support the disaster victim identification teams in the field.

Post-mortem radiological CT identification based on classical ante-mortem X-ray examinations

Radiological identification is important in forensic medicine. Identification using comparison of individualising structures with ante- and post-mortem conventional radiographs has been known for a long time. New radiological procedures such as computed tomography (CT) and magnetic resonance imaging (MRI) are being increasingly used for identification. In this paper, a new comparative approach using various radiological methods is described and its application demonstrated. This new approach is the comparison of ante-mortem conventional radiographs with projected images calculated from post-mortem CT data. The identification procedure will be illustrated with reference to the frontal sinus and the pelvis.

Finite element 3D reconstruction of the pulmonary acinus imaged by synchrotron X-ray tomography

The alveolated structure of the pulmonary acinus plays a vital role in gas exchange function. Three-dimensional (3D) analysis of the parenchymal region is fundamental to understanding this structure-function relationship, but only a limited number of attempts have been conducted in the past because of technical limitations. In this study, we developed a new image processing methodology based on finite element (FE) analysis for accurate 3D structural reconstruction of the gas exchange regions of the lung. Stereologically well characterized rat lung samples (Pediatr Res 53: 72-80, 2003) were imaged using high-resolution synchrotron radiation-based X-ray tomographic microscopy. A stack of 1,024 images (each slice: 1024 x 1024 pixels) with resolution of 1.4 mum(3) per voxel were generated. For the development of FE algorithm, regions of interest (ROI), containing approximately 7.5 million voxels, were further extracted as a working subunit. 3D FEs were created overlaying the voxel map using a grid-based hexahedral algorithm. A proper threshold value for appropriate segmentation was iteratively determined to match the calculated volume density of tissue to the stereologically determined value (Pediatr Res 53: 72-80, 2003). The resulting 3D FEs are ready to be used for 3D structural analysis as well as for subsequent FE computational analyses like fluid dynamics and skeletonization.

DUAL ENERGY COMPUTED TOMOGRAPHY FOR PROTON THERAPY TREATMENT PLANNING

Proton radiation therapy is gaining popularity because of the unique characteristics of its dose distribution, e.g., high dose-gradient at the distal end of the percentage-depth-dose curve (known as the Bragg peak). The high dose-gradient offers the possibility of delivering high dose to the target while still sparing critical organs distal to the target. However, the high dose-gradient is a double-edged sword: a small shift of the highly conformal high-dose area can cause the target to be substantially under-dosed or the critical organs to be substantially over-dosed. Because of that, large margins are required in treatment planning to ensure adequate dose coverage of the target, which prevents us from realizing the full potential of proton beams. Therefore, it is critical to reduce uncertainties in the proton radiation therapy. One major uncertainty in a proton treatment is the range uncertainty related to the estimation of proton stopping power ratio (SPR) distribution inside a patient. The SPR distribution inside a patient is required to account for tissue heterogeneities when calculating dose distribution inside the patient. In current clinical practice, the SPR distribution inside a patient is estimated from the patient’s treatment planning computed tomography (CT) images based on the CT number-to-SPR calibration curve. The SPR derived from a single CT number carries large uncertainties in the presence of human tissue composition variations, which is the major drawback of the current SPR estimation method. We propose to solve this problem by using dual energy CT (DECT) and hypothesize that the range uncertainty can be reduced by a factor of two from currently used value of 3.5%. A MATLAB program was developed to calculate the electron density ratio (EDR) and effective atomic number (EAN) from two CT measurements of the same object. An empirical relationship was discovered between mean excitation energies and EANs existing in human body tissues. With the MATLAB program and the empirical relationship, a DECT-based method was successfully developed to derive SPRs for human body tissues (the DECT method). The DECT method is more robust against the uncertainties in human tissues compositions than the current single-CT-based method, because the DECT method incorporated both density and elemental composition information in the SPR estimation. Furthermore, we studied practical limitations of the DECT method. We found that the accuracy of the DECT method using conventional kV-kV x-ray pair is susceptible to CT number variations, which compromises the theoretical advantage of the DECT method. Our solution to this problem is to use a different x-ray pair for the DECT. The accuracy of the DECT method using different combinations of x-ray energies, i.e., the kV-kV, kV-MV and MV-MV pair, was compared using the measured imaging uncertainties for each case. The kV-MV DECT was found to be the most robust against CT number variations. In addition, we studied how uncertainties propagate through the DECT calculation, and found general principles of selecting x-ray pairs for the DECT method to minimize its sensitivity to CT number variations. The uncertainties in SPRs estimated using the kV-MV DECT were analyzed further and compared to those using the stoichiometric method. The uncertainties in SPR estimation can be divided into five categories according to their origins: the inherent uncertainty, the DECT modeling uncertainty, the CT imaging uncertainty, the uncertainty in the mean excitation energy, and SPR variation with proton energy. Additionally, human body tissues were divided into three tissue groups – low density (lung) tissues, soft tissues and bone tissues. The uncertainties were estimated separately because their uncertainties were different under each condition. An estimate of the composite range uncertainty (2s) was determined for three tumor sites – prostate, lung, and head-and-neck, by combining the uncertainty estimates of all three tissue groups, weighted by their proportions along typical beam path for each treatment site. In conclusion, the DECT method holds theoretical advantages in estimating SPRs for human tissues over the current single-CT-based method. Using existing imaging techniques, the kV-MV DECT approach was capable of reducing the range uncertainty from the currently used value of 3.5% to 1.9%-2.3%, but it is short to reach our original goal of reducing the range uncertainty by a factor of two. The dominant source of uncertainties in the kV-MV DECT was the uncertainties in CT imaging, especially in MV CT imaging. Further reduction in beam hardening effect, the impact of scatter, out-of-field object etc. would reduce the Hounsfeld Unit variations in CT imaging. The kV-MV DECT still has the potential to reduce the range uncertainty further.