Metal artifacts from titanium and steel screws in CT, 1.5T and 3T MR images of the tibial Pilon: A quantitative assessment in 3D

Radiographs are commonly used to assess articular reduction of the distal tibia (pilon) fractures postoperatively, but may reveal malreductions inaccurately. While Magnetic Resonance Imaging (MRI) and Computed Tomography (CT) are potential 3D alternatives they generate metal-related artifacts. This study aims to quantify the artifact size from orthopaedic screws using CT, 1.5T and 3T MRI data. Three screws were inserted into one intact human cadaver ankle specimen proximal to and along the distal articular surface, then CT, 1.5T and 3T MRI scanned. Four types of screws were investigated: titanium alloy (TA), stainless steel (SS) (Ø = 3.5 mm), cannulated TA (CTA) and cannulated SS (CSS)(Ø = 4.0 mm, Ø empty core = 2.6 mm). 3D artifact models were reconstructed using adaptive thresholding. The artifact size was measured by calculating the perpendicular distance from the central screw axis to the boundary of the artifact in four anatomical directions with respect to the distal tibia. The artifact sizes (in the order of TA, SS, CTA and CSS) from CT were 2.0 mm, 2.6 mm, 1.6 mm and 2.0 mm; from 1.5T MRI they were 3.7 mm, 10.9 mm, 2.9 mm, and 9 mm; and 3T MRI they were 4.4 mm, 15.3 mm, 3.8 mm, and 11.6 mm respectively. Therefore, CT can be used as long as the screws are at a safe distance of about 2 mm from the articular surface. MRI can be used if the screws are at least 3 mm away from the articular surface except SS and CSS. Artifacts from steel screws were too large thus obstructed the pilon from being visualised in MRI. Significant differences (P < 0.05) were found in the size of artifacts between all imaging modalities, screw types and material types, except 1.5T versus 3T MRI for the SS screws (P = 0.063). CTA screws near the joint surface can improve postoperative assessment in CT and MRI. MRI presents a favourable non-ionising alternative when using titanium hardware. Since these factors may influence the quality of postoperative assessment, potential improvements in operative techniques should be considered.

Assessment of the effects of CT dose in averaged x-ray CT images of a dose-sensitive polymer gel

The signal-to-noise ratio achievable in x-ray computed tomography (CT) images of polymer gels can be increased by averaging over multiple scans of each sample. However, repeated scanning delivers a small additional dose to the gel which may compromise the accuracy of the dose measurement. In this study, a NIPAM-based polymer gel was irradiated and then CT scanned 25 times, with the resulting data used to derive an averaged image and a "zero-scan" image of the gel. Comparison between these two results and the first scan of the gel showed that the averaged and zero-scan images provided better contrast, higher contrast-to- noise and higher signal-to-noise than the initial scan. The pixel values (Hounsfield units, HU) in the averaged image were not noticeably elevated, compared to the zero-scan result and the gradients used in the linear extrapolation of the zero-scan images were small and symmetrically distributed around zero. These results indicate that the averaged image was not artificially lightened by the small, additional dose delivered during CT scanning. This work demonstrates the broader usefulness of the zero-scan method as a means to verify the dosimetric accuracy of gel images derived from averaged x-ray CT data.

Bayesian computational methods for spatial analysis of images

This thesis introduces a new way of using prior information in a spatial model and develops scalable algorithms for fitting this model to large imaging datasets. These methods are employed for image-guided radiation therapy and satellite based classification of land use and water quality. This study has utilized a pre-computation step to achieve a hundredfold improvement in the elapsed runtime for model fitting. This makes it much more feasible to apply these models to real-world problems, and enables full Bayesian inference for images with a million or more pixels.

A MatLAB Based Virtual Phantom for 2D Electrical Impedance Tomography (MatVP2DEIT): Studying the Medical Electrical Impedance Tomography Reconstruction in Computer

Electrical Impedance Tomography (EIT) is a computerized medical imaging technique which reconstructs the electrical impedance images of a domain under test from the boundary voltage-current data measured by an EIT electronic instrumentation using an image reconstruction algorithm. Being a computed tomography technique, EIT injects a constant current to the patient's body through the surface electrodes surrounding the domain to be imaged (Omega) and tries to calculate the spatial distribution of electrical conductivity or resistivity of the closed conducting domain using the potentials developed at the domain boundary (partial derivative Omega). Practical phantoms are essentially required to study, test and calibrate a medical EIT system for certifying the system before applying it on patients for diagnostic imaging. Therefore, the EIT phantoms are essentially required to generate boundary data for studying and assessing the instrumentation and inverse solvers a in EIT. For proper assessment of an inverse solver of a 2D EIT system, a perfect 2D practical phantom is required. As the practical phantoms are the assemblies of the objects with 3D geometries, the developing of a practical 2D-phantom is a great challenge and therefore, the boundary data generated from the practical phantoms with 3D geometry are found inappropriate for assessing a 2D inverse solver. Furthermore, the boundary data errors contributed by the instrumentation are also difficult to separate from the errors developed by the 3D phantoms. Hence, the errorless boundary data are found essential to assess the inverse solver in 2D EIT. In this direction, a MatLAB-based Virtual Phantom for 2D EIT (MatVP2DEIT) is developed to generate accurate boundary data for assessing the 2D-EIT inverse solvers and the image reconstruction accuracy. MatVP2DEIT is a MatLAB-based computer program which simulates a phantom in computer and generates the boundary potential data as the outputs by using the combinations of different phantom parameters as the inputs to the program. Phantom diameter, inhomogeneity geometry (shape, size and position), number of inhomogeneities, applied current magnitude, background resistivity, inhomogeneity resistivity all are set as the phantom variables which are provided as the input parameters to the MatVP2DEIT for simulating different phantom configurations. A constant current injection is simulated at the phantom boundary with different current injection protocols and boundary potential data are calculated. Boundary data sets are generated with different phantom configurations obtained with the different combinations of the phantom variables and the resistivity images are reconstructed using EIDORS. Boundary data of the virtual phantoms, containing inhomogeneities with complex geometries, are also generated for different current injection patterns using MatVP2DEIT and the resistivity imaging is studied. The effect of regularization method on the image reconstruction is also studied with the data generated by MatVP2DEIT. Resistivity images are evaluated by studying the resistivity parameters and contrast parameters estimated from the elemental resistivity profiles of the reconstructed phantom domain. Results show that the MatVP2DEIT generates accurate boundary data for different types of single or multiple objects which are efficient and accurate enough to reconstruct the resistivity images in EIDORS. The spatial resolution studies show that, the resistivity imaging conducted with the boundary data generated by MatVP2DEIT with 2048 elements, can reconstruct two circular inhomogeneities placed with a minimum distance (boundary to boundary) of 2 mm. It is also observed that, in MatVP2DEIT with 2048 elements, the boundary data generated for a phantom with a circular inhomogeneity of a diameter less than 7% of that of the phantom domain can produce resistivity images in EIDORS with a 1968 element mesh. Results also show that the MatVP2DEIT accurately generates the boundary data for neighbouring, opposite reference and trigonometric current patterns which are very suitable for resistivity reconstruction studies. MatVP2DEIT generated data are also found suitable for studying the effect of the different regularization methods on reconstruction process. Comparing the reconstructed image with an original geometry made in MatVP2DEIT, it would be easier to study the resistivity imaging procedures as well as the inverse solver performance. Using the proposed MatVP2DEIT software with modified domains, the cross sectional anatomy of a number of body parts can be simulated in PC and the impedance image reconstruction of human anatomy can be studied.

Extraction of fibre network architecture by X-ray tomography and prediction of elastic properties using an affine analytical model

This paper proposes a method for extracting reliable architectural characteristics from complex porous structures using micro-computed tomography (μCT) images. The work focuses on a highly porous material composed of a network of fibres bonded together. The segmentation process, allowing separation of the fibres from the remainder of the image, is the most critical step in constructing an accurate representation of the network architecture. Segmentation methods, based on local and global thresholding, were investigated and evaluated by a quantitative comparison of the architectural parameters they yielded, such as the fibre orientation and segment length (sections between joints) distributions and the number of inter-fibre crossings. To improve segmentation accuracy, a deconvolution algorithm was proposed to restore the original images. The efficacy of the proposed method was verified by comparing μCT network architectural characteristics with those obtained using high resolution CT scans (nanoCT). The results indicate that this approach resolves the architecture of these complex networks and produces results approaching the quality of nanoCT scans. The extracted architectural parameters were used in conjunction with an affine analytical model to predict the axial and transverse stiffnesses of the fibre network. Transverse stiffness predictions were compared with experimentally measured values obtained by vibration testing. © 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

The geometry and the number of contacts of monodisperse sphere packs using X-ray tomography

The three-dimensional structure of very large samples of monodisperse bead packs is studied by means of X-Ray Computed Tomography. We retrieve the coordinatesofeach bead inthe pack and wecalculate the average coordination number by using the tomographic images to single out the neighbors in contact. The results are compared with the average coordination number obtained in Aste et al. (2005) by using a deconvolution technique. We show that the coordination number increases with the packing fraction, varying between 6.9 and 8.2 for packing fractions between 0.59 and 0.64. © 2005 Taylor & Francis Group.

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.

Computed tomographic-dacryocystography (CT-DCG) of the normal canine nasolacrimal drainage system with three-dimensional reconstruction

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Development and applications of three-dimensional gamma ray tomography system using ray casting volume rendering techniques

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Use of cone-beam tomography and digital subtraction radiography for diagnosis and evaluation of traumatized teeth treated with endodontic surgery and MTA. A case report

The aim of this paper was to report the diagnosis and evaluation of periradicular bone repair by using computed tomography and digital subtraction radiography in an endodontic surgery case treated with mineral trioxide aggregate (MTA). The patient had local swelling and fistula 7years after trauma in the upper jaw. Periapical lesion on tooth #12 and root resorption of tooth #11 were detected radiographically. Endodontic therapy and placement of intracanal medication were carried out. After 3months, with no improvement in the clinical signs and symptoms, computed tomography was performed, showing extensive apical bone resorption on tooth #12 and dental resorption promoting communication of the root canal with the periodontium of tooth #11. The patient was referred to endodontic surgery. After surgery and postoperative periods of 15, 30, 60, 90, and 120days, standardized radiographs were taken, digitized, and subjected to digital subtraction of the images using Adobe Photoshop CS software. Four years later, a cone-beam computed tomography was performed, showing bone repair and absence of root resorption at tooth #11. In this clinical case, digital subtraction radiography was effective for early detection of new bone formation and evolution of repair. © 2011 John Wiley & Sons A/S.

Application of two segmentation protocols during the processing of virtual images in rapid prototyping: ex vivo study with human dry mandibles

The aim of this study was to evaluate the accuracy of virtual three-dimensional (3D) reconstructions of human dry mandibles, produced from two segmentation protocols (outline only and all-boundary lines).Twenty virtual three-dimensional (3D) images were built from computed tomography exam (CT) of 10 dry mandibles, in which linear measurements between anatomical landmarks were obtained and compared to an error probability of 5 %.The results showed no statistically significant difference among the dry mandibles and the virtual 3D reconstructions produced from segmentation protocols tested (p = 0,24).During the designing of a virtual 3D reconstruction, both outline only and all-boundary lines segmentation protocols can be used.Virtual processing of CT images is the most complex stage during the manufacture of the biomodel. Establishing a better protocol during this phase allows the construction of a biomodel with characteristics that are closer to the original anatomical structures. This is essential to ensure a correct preoperative planning and a suitable treatment.

Wavelet-based algorithm to the evaluation of contrasted hepatocellular carcinoma in CT-images after transarterial chemoembolization

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

SPECT (Single photon emission tomography): Gama câmara, reconstrução tomográfica e características funcional

The SPECT (Single Photon Emission Computed Tomography) systems are part of a medical image acquisition technology which has been outstanding, because the resultant images are functional images complementary to those that give anatomic information, such as X-Ray CT, presenting a high diagnostic value. These equipments acquire, in a non-invasive way, images from the interior of the human body through tomographic mapping of radioactive material administered to the patient. The SPECT systems are based on the Gamma Camera detection system, and one of them being set on a rotational gantry is enough to obtain the necessary data for a tomographic image. The images obtained from the SPECT system consist in a group of flat images that describe the radioactive distribution on the patient. The trans-axial cuts are obtained from the tomographic reconstruction techniques. There are analytic and iterative methods to obtain the tomographic reconstruction. The analytic methods are based on the Fourier Cut Theorem (FCT), while the iterative methods search for numeric solutions to solve the equations from the projections. Within the analytic methods, the filtered backprojection (FBP) method maybe is the simplest of all the tomographic reconstruction techniques. This paper's goal is to present the operation of the SPECT system, the Gamma Camera detection system, some tomographic reconstruction techniques and the requisites for the implementation of this system in a Nuclear Medicine service

Cone Beam Computed Tomographic Evaluation of a Maxillary Alveolar Ridge Reconstruction With Iliac Crest Graft and Implants

The present article discusses an atrophic maxilla reconstruction with iliac crest bone block and particulate grafts and dental implants. Onlay block grafts were used to restore bone volume of the anterior maxilla, whereas bilateral sinus floor augmentation was performed using a particulate graft. Ten months after the grafting surgery, 9 dental implants were placed to rehabilitate the case. Results of a 7-year follow-up were obtained clinically and by cone beam computed tomographic images.

Regional lung perfusion estimated by electrical impedance tomography in a piglet model of lung collapse

Borges JB, Suarez-Sipmann F, Bohm SH, Tusman G, Melo A, Maripuu E, Sandstrom M, Park M, Costa EL, Hedenstierna G, Amato M. Regional lung perfusion estimated by electrical impedance tomography in a piglet model of lung collapse. J Appl Physiol 112: 225-236, 2012. First published September 29, 2011; doi: 10.1152/japplphysiol.01090.2010.-The assessment of the regional match between alveolar ventilation and perfusion in critically ill patients requires simultaneous measurements of both parameters. Ideally, assessment of lung perfusion should be performed in real-time with an imaging technology that provides, through fast acquisition of sequential images, information about the regional dynamics or regional kinetics of an appropriate tracer. We present a novel electrical impedance tomography (EIT)-based method that quantitatively estimates regional lung perfusion based on first-pass kinetics of a bolus of hypertonic saline contrast. Pulmonary blood flow was measured in six piglets during control and unilateral or bilateral lung collapse conditions. The first-pass kinetics method showed good agreement with the estimates obtained by single-photon-emission computerized tomography (SPECT). The mean difference (SPECT minus EIT) between fractional blood flow to lung areas suffering atelectasis was -0.6%, with a SD of 2.9%. This method outperformed the estimates of lung perfusion based on impedance pulsatility. In conclusion, we describe a novel method based on EIT for estimating regional lung perfusion at the bedside. In both healthy and injured lung conditions, the distribution of pulmonary blood flow as assessed by EIT agreed well with the one obtained by SPECT. The method proposed in this study has the potential to contribute to a better understanding of the behavior of regional perfusion under different lung and therapeutic conditions.