Study of OSEM with different subsets in grating-based X-ray differential phase-contrast imaging.

Impressive developments in X-ray imaging are associated with X-ray phase contrast computed tomography based on grating interferometry, a technique that provides increased contrast compared with conventional absorption-based imaging. A new "single-step" method capable of separating phase information from other contributions has been recently proposed. This approach not only simplifies data-acquisition procedures, but, compared with the existing phase step approach, significantly reduces the dose delivered to a sample. However, the image reconstruction procedure is more demanding than for traditional methods and new algorithms have to be developed to take advantage of the "single-step" method. In the work discussed in this paper, a fast iterative image reconstruction method named OSEM (ordered subsets expectation maximization) was applied to experimental data to evaluate its performance and range of applicability. The OSEM algorithm with different subsets was also characterized by comparison of reconstruction image quality and convergence speed. Computer simulations and experimental results confirm the reliability of this new algorithm for phase-contrast computed tomography applications. Compared with the traditional filtered back projection algorithm, in particular in the presence of a noisy acquisition, it furnishes better images at a higher spatial resolution and with lower noise. We emphasize that the method is highly compatible with future X-ray phase contrast imaging clinical applications.

Investigation of discrete imaging models and iterative image reconstruction in differential X-ray phase-contrast tomography.

Differential X-ray phase-contrast tomography (DPCT) refers to a class of promising methods for reconstructing the X-ray refractive index distribution of materials that present weak X-ray absorption contrast. The tomographic projection data in DPCT, from which an estimate of the refractive index distribution is reconstructed, correspond to one-dimensional (1D) derivatives of the two-dimensional (2D) Radon transform of the refractive index distribution. There is an important need for the development of iterative image reconstruction methods for DPCT that can yield useful images from few-view projection data, thereby mitigating the long data-acquisition times and large radiation doses associated with use of analytic reconstruction methods. In this work, we analyze the numerical and statistical properties of two classes of discrete imaging models that form the basis for iterative image reconstruction in DPCT. We also investigate the use of one of the models with a modern image reconstruction algorithm for performing few-view image reconstruction of a tissue specimen.

Anthropomorphic model observer performance in three-dimensional detection task for low-contrast computed tomography.

X-ray medical imaging is increasingly becoming three-dimensional (3-D). The dose to the population and its management are of special concern in computed tomography (CT). Task-based methods with model observers to assess the dose-image quality trade-off are promising tools, but they still need to be validated for real volumetric images. The purpose of the present work is to evaluate anthropomorphic model observers in 3-D detection tasks for low-contrast CT images. We scanned a low-contrast phantom containing four types of signals at three dose levels and used two reconstruction algorithms. We implemented a multislice model observer based on the channelized Hotelling observer (msCHO) with anthropomorphic channels and investigated different internal noise methods. We found a good correlation for all tested model observers. These results suggest that the msCHO can be used as a relevant task-based method to evaluate low-contrast detection for CT and optimize scan protocols to lower dose in an efficient way.

Statistical Inversion Theory in X-ray Tomography

In this thesis the X-ray tomography is discussed from the Bayesian statistical viewpoint. The unknown parameters are assumed random variables and as opposite to traditional methods the solution is obtained as a large sample of the distribution of all possible solutions. As an introduction to tomography an inversion formula for Radon transform is presented on a plane. The vastly used filtered backprojection algorithm is derived. The traditional regularization methods are presented sufficiently to ground the Bayesian approach. The measurements are foton counts at the detector pixels. Thus the assumption of a Poisson distributed measurement error is justified. Often the error is assumed Gaussian, altough the electronic noise caused by the measurement device can change the error structure. The assumption of Gaussian measurement error is discussed. In the thesis the use of different prior distributions in X-ray tomography is discussed. Especially in severely ill-posed problems the use of a suitable prior is the main part of the whole solution process. In the empirical part the presented prior distributions are tested using simulated measurements. The effect of different prior distributions produce are shown in the empirical part of the thesis. The use of prior is shown obligatory in case of severely ill-posed problem.

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.