Home Medicine Experimental computer tomograph
Article Open Access

Experimental computer tomograph

  • D. Heinemann EMAIL logo , A. Keller and D. Jannek
Published/Copyright: September 12, 2015

Abstract

The computed tomography is one of the most important medical instruments, allowing the non-invasive visualization of cross sections which are free from superpositions. Since 2000 an experimental computer tomo-graph of the third generation for the purpose of education and research was set up and further developed. Besides the mechanical construction design reconstruction algorithms, including certain corrections of the measured data were developed and implemented. In 2013 iterative reconstruction methods were investigated and implemented for advanced reconstructions and dose reduction using various ray tracing algorithms. The new reconstruction technique leads to improvements in image quality and low dose reconstructions.

1 Experimental CT

The computer tomograph of the Institute for Biomedical Engineering and Informatics is an inverse third generation. While the original third generation CTs have a rotary gantry, our system is designed with a turntable, where the examined object is placed between a fixed source-detector-arrangement (see Figure 1). This design reduces the mechanical construction effort, while preserving the basic principle of the third CT generation. Using a linear stage the turntable can be moved vertically through the acquisition geometry. This allows for multiple slice or spiral acquisition for 3D reconstructions. As the detector an industrial X-ray sensible CCD line with optimal quantum effi-ciency, high spatial resolution, stability of measured data over large acquisition periods and calibration possibilities is used. Due to the system geometry objects with a maximum diameter of 200 mm can be measured and reconstructed. The detector has a spatial resolution of 5 lp/mm, a sensitivity range of 40-160 keV and a sensitive detector length of 214 mm. The measured raw data has a range of values from 0 to 4095 which corresponds to a dynamic range of 12 bit. The resulting sinogram has to be normalized since the raw data represents just an intensity distribution. The correction is based on the attenuation law:

(1)Icorr=lnII0.
Figure 1 Components and setup of the experimental CT
Figure 1

Components and setup of the experimental CT

Each intensity value of the detector is corrected by this formula, where I0 is the intensity of the unattenuated radiation. In consequence of the low-pass characteristic the simple backprojection has to be extended by filtering the projection data, thus the filtered backprojection is obtained for image reconstruction. More detailed information about the system and its behaviour are given in [16].

2 Iterative reconstruction and hardware acceleration

Iterative reconstruction techniques are as old as the computed tomography itself. Due to restrictions in computer technology they were not practical used over the last decades. Recent developments made these techniques more and more applicable. Since 2009 they can be found in the clinical environment for X-ray computed tomography.

2.1 Gereral advantages

One of the most important reasons for using iterative reconstruction techniques is the optimization of diagnostic relevant image information and patient dose. Using the filtered backprojection image noise and dose are interrelated via a so called “root-factor”. While doubling the dose, the noise is only reduced by the square root of 2. If the half of the noise should remain, then the dose has to be 4 times greater. Iterative methods are able to undock the noise from the spatial resolution and can reduce the dose while maintaining or even enhance the image quality und thus the diagnostic relevance.

2.2 Iterative reconstruction with the experimental CT

The basic principle of iterative reconstruction methods can be seen in Figure 2 and consists of a simulation of the measuring process by calculating synthetic projection data. Comparing the synthetic sinogram with original projection data delivers an error between those datasets. Using this error the initial reconstruction is corrected. This step is weighted by a relaxation parameter which controls the convergence of the method. The obtained error in radon space between synthetic and original data can be used for image correction. This can easily be done by using the simple backprojection with the error data. To calculate this error data the synthetic sinogram has to be created. The synthetic projections are the result of specific ray tracing algorithms. A simple and a more accurate version are described in the following sections. The calculation of synthetic data and image correction is performed in loops, thus the iterative methods are more computationally intensive as the simple or filtered backprojection.

Figure 2 Basic principle of iterative reconstruction
Figure 2

Basic principle of iterative reconstruction

2.2.1 Simple ray tracing

This version of ray tracing does not consider a ray between X-ray source and detector but rays between source and image pixels. Each ray defines a detector element where the current pixel value gives a contribution to (see Figure 3). Via geometric considerations the position of the detector element can be calculated using the formula:

(2)t=th+RFODX·FDD

FOD defines the focus-object-distance, FDD the focus-detector-distance,R describes the distance between pixel midpoint and normal beam, while X determines the distance from the center of rotation to the perpendicular of the current pixel to the normal ray and th is the half of the detector length. These variables get different values for each projection angle φ:

(3)R(xa,ya,ϕ)=R1+R2=cos(φ)·ya+sin(φ)·xa
(4)X(xa,ya,φ)=X1X2=cos(φ)·xasin(φ)·ya

Inserting these expressions in the formula above one obtains:

(5)t=tn+(cos(φ)·ya+sin(φ)·xa)·FDDFOD(cos(φ)·xasin(φ)·ya)

Subsequent the current pixel value is added to the detector element value at position t.

Figure 3 Simple ray tracing
Figure 3

Simple ray tracing

2.2.2 Exact ray tracing

This more accurate version simulates rays from the source to the detector and calculates the pixels which influence the detector value (see Figure 4). At the beginning the entrance point of the projection ray into the image is calculated. xa holds the value of the left image border if the ray hits the image at this edge and Rg describes whether the ray intersects the upper or the lower half of the image. Due to the consideration of one certain detector element, t1 is known. The projection rays deflection with respect to the normal ray and to the x-axis is described as:

(6)β1=arcsin(t1FDD)andα1=φβ1
Figure 4 Exact ray tracing
Figure 4

Exact ray tracing

Through the knowledge of the spatial source position (X,Y ) the exact intersection of the ray can be determined:

(7)X=cos(φ)·FODandY=sin(φ)·FOD
(8)Rg=Yyi1
(9)yi1=tan(α1)·xi1andxi1=Xxa

The ray intersects the image at the left border if yi1 is between the upper and lower image border. If not the next boarder is checked and Xg is examined. The calculations of α2 and β2 are analogous.

(10)Xg=Xxi2
(11)xi2=tan(α2)·yi2andyi2=Yya

Through additional knowledge of pixel size and resolution the first intesected pixel can be determined. Starting from this pixel the projection ray is followed and the values are added to the corresponding detector element, weighted by the length of the projection ray in each pixel. Using this weighting causes a slight error in the reprojected data. As Figure 4 shows only pixels on the projection ray influence the detector value and pixels between rays are missed. This problem can be solved by using area weights like the large dark area in Figure 4.

2.3 Reconstruction using hardware acceleration

The simple sequential processing gives a proper opportunity to compare old and modern computer technology. In the previous version of the reconstruction computer an AMD XP 2200 processor (1 x 1.8 GHz) was installed, while the current system uses an Intel Core i7 950 (4 x 3.07 GHz). The systems deliver a processing time of 50 seconds for the old and 3 seconds for the modern system. As a first real acceleration method OpenMP is used. OpenMP (Open Multi-Processing) gives the opportunity to run code on all available CPU cores in parallel. By parallelizing all for-loops the processing was accelerated by 40%. NVIDIAs CUDA is used to create high performance reconstructions. Modern graphics cards have hundreds of processing units which operate in parallel. The code was designed in such a way that each processing unit of the graphics card runs calculations for only one pixel. By using this technology the reconstruction was done in 42 milliseconds. This means an acceleration by the factor 1190.

3 Results

We acquired a phantom (see Figure 5) with four radiation doses. The acquisition parameters where 100 kV, 1.5 mA. Varying the integration time (0.05 ms; 0.1 ms; 0.15 ms; 0.2 ms results in dose reductions of 25 %, 50 % and 75 %. Figure 5 shows reconstructions with 25 % (first row) and 100 % (second row) dose using the filtered backprojection (left) and the developed iterative reconstruction (right). As mentioned in the previous section the noise directly corresponds with the applied dose what can clearly be seen in Figure 5. The first row shows a significant worse image quality due to increased noise. Using the iterative method the image noise is reduced and spatial resolution is maintained. Even a dose reduction of 75 % gives a better reconstruction result as the filtered backprojection with 100 % dose. Variations of the acquired data and reconstructions obtained from certain parameter settings of the iterative algorithms are demonstrated in [7].

Figure 5 Reconstructions with different doses
Figure 5

Reconstructions with different doses

4 Conclusion

The investigated algorithms and the implementation extend the experimental CT with respect to modern reconstruction techniques and technology. This gives forward looking opportunities in education and research in the field of computed tomography. Even with much lower radiation dose the reconstruction with consistent or better image quality is possible. Significant acceleration of the processing was reached by using modern computer technology. The software and algorithms are still under further development and will give additional improvements in future work.

Author's Statement

  1. Conflict of interest: Authors state no conflict of interest. Material and Methods: Informed consent: Informed consent is not applicable. Ethical approval: The conducted research is not related to either human or animals use.

References

[1] Keller A.: Experimenteller Computertomograph für Ausbildung und Forschung – Teil 1: Modellierung der Querschnittsrekonstruktion. In: medizintechnik 2002 (3), 122: 103–105.Search in Google Scholar

[2] Gross S., Keller A., Krüger U.: Experimenteller Computertomo-graph für Ausbildung und Forschung – Teil 2: Problemanalyse, Konzept, Bildrekonstruktion. In: medizintechnik 2002 (4), 122: 135–140.Search in Google Scholar

[3] Kyriakou Y., Keller A.: Experimenteller Computertomograph für Ausbildung und Forschung – Teil 3: Aufhärtungskorrektur. In: medizintechnik 2004 (1), 124: 15–19.Search in Google Scholar

[4] Keller A., Scholz N., Schmidt F.: Experimenteller Computertomo-graph für Ausbildung und Forschung – Teil 4: Örtliche Dynamik des Zeilendetektors. In: medizintechnik 2005 (4), 125: 131–137.Search in Google Scholar

[5] Keller A., Heinz Chr.: Experimenteller Computertomograph für Ausbildung und Forschung – Teil 5: Schwächungsmessung, CT-Zahlen. In: medizintechnik 2007 (2), 127: 71–75.Search in Google Scholar

[6] Keller A., Herzog S.: Experimenteller Computertomograph für Ausbildung und Forschung – Teil 6: Erfahrungen und Ergebnisse. In: medizintechnik 2011 (6), 131: 224-228.Search in Google Scholar

[7] Heinemann D.: Methodischer Vergleich iterativer Rekonstruktionsverfahren bei der Computertomographie, Master Thesis, Technische Universität Ilmenau 2013: 160.Search in Google Scholar

Published Online: 2015-9-12
Published in Print: 2015-9-1

© 2015 by Walter de Gruyter GmbH, Berlin/Boston

This article is distributed under the terms of the Creative Commons Attribution Non-Commercial License, which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Articles in the same Issue

  1. Research Article
  2. Development and characterization of superparamagnetic coatings
  3. Research Article
  4. The development of an experimental setup to measure acousto-electric interaction signal
  5. Research Article
  6. Stability analysis of ferrofluids
  7. Research Article
  8. Investigation of endothelial growth using a sensors-integrated microfluidic system to simulate physiological barriers
  9. Research Article
  10. Energy harvesting for active implants: powering a ruminal pH-monitoring system
  11. Research Article
  12. New type of fluxgate magnetometer for the heart’s magnetic fields detection
  13. Research Article
  14. Field mapping of ballistic pressure pulse sources
  15. Research Article
  16. Development of a new homecare sleep monitor using body sounds and motion tracking
  17. Research Article
  18. Noise properties of textile, capacitive EEG electrodes
  19. Research Article
  20. Detecting phase singularities and rotor center trajectories based on the Hilbert transform of intraatrial electrograms in an atrial voxel model
  21. Research Article
  22. Spike sorting: the overlapping spikes challenge
  23. Research Article
  24. Separating the effect of respiration from the heart rate variability for cases of constant harmonic breathing
  25. Research Article
  26. Locating regions of arrhythmogenic substrate by analyzing the duration of triggered atrial activities
  27. Research Article
  28. Combining different ECG derived respiration tracking methods to create an optimal reconstruction of the breathing pattern
  29. Research Article
  30. Atrial and ventricular signal averaging electrocardiography in pacemaker and cardiac resynchronization therapy
  31. Research Article
  32. Estimation of a respiratory signal from a single-lead ECG using the 4th order central moments
  33. Research Article
  34. Compressed sensing of multi-lead ECG signals by compressive multiplexing
  35. Research Article
  36. Heart rate monitoring in ultra-high-field MRI using frequency information obtained from video signals of the human skin compared to electrocardiography and pulse oximetry
  37. Research Article
  38. Synchronization in wireless biomedical-sensor networks with Bluetooth Low Energy
  39. Research Article
  40. Automated classification of stages of anaesthesia by populations of evolutionary optimized fuzzy rules
  41. Research Article
  42. Effects of sampling rate on automated fatigue recognition in surface EMG signals
  43. Research Article
  44. Closed-loop transcranial alternating current stimulation of slow oscillations
  45. Research Article
  46. Cardiac index in atrio- and interventricular delay optimized cardiac resynchronization therapy and cardiac contractility modulation
  47. Research Article
  48. The role of expert evaluation for microsleep detection
  49. Research Article
  50. The impact of baseline wander removal techniques on the ST segment in simulated ischemic 12-lead ECGs
  51. Research Article
  52. Metal artifact reduction by projection replacements and non-local prior image integration
  53. Research Article
  54. A novel coaxial nozzle for in-process adjustment of electrospun scaffolds’ fiber diameter
  55. Research Article
  56. Processing of membranes for oxygenation using the Bellhouse-effect
  57. Research Article
  58. Inkjet printing of viable human dental follicle stem cells
  59. Research Article
  60. The use of an icebindingprotein out of the snowflea Hypogastrura harveyi as a cryoprotectant in the cryopreservation of mesenchymal stem cells
  61. Research Article
  62. New NIR spectroscopy based method to determine ischemia in vivo in liver – a first study on rats
  63. Research Article
  64. QRS and QT ventricular conduction times and permanent pacemaker therapy after transcatheter aortic valve implantation
  65. Research Article
  66. Adopting oculopressure tonometry as a transient in vivo rabbit glaucoma model
  67. Research Article
  68. Next-generation vision testing: the quick CSF
  69. Research Article
  70. Improving tactile sensation in laparoscopic surgery by overcoming size restrictions
  71. Research Article
  72. Design and control of a 3-DOF hydraulic driven surgical instrument
  73. Research Article
  74. Evaluation of endourological tools to improve the diagnosis and therapy of ureteral tumors – from model development to clinical application
  75. Research Article
  76. Frequency based assessment of surgical activities
  77. Research Article
  78. “Hands free for intervention”, a new approach for transoral endoscopic surgery
  79. Research Article
  80. Pseudo-haptic feedback in medical teleoperation
  81. Research Article
  82. Feasibility of interactive gesture control of a robotic microscope
  83. Research Article
  84. Towards structuring contextual information for workflow-driven surgical assistance functionalities
  85. Research Article
  86. Towards a framework for standardized semantic workflow modeling and management in the surgical domain
  87. Research Article
  88. Closed-loop approach for situation awareness of medical devices and operating room infrastructure
  89. Research Article
  90. Kinect based physiotherapy system for home use
  91. Research Article
  92. Evaluating the microsoft kinect skeleton joint tracking as a tool for home-based physiotherapy
  93. Research Article
  94. Integrating multimodal information for intraoperative assistance in neurosurgery
  95. Research Article
  96. Respiratory motion tracking using Microsoft’s Kinect v2 camera
  97. Research Article
  98. Using smart glasses for ultrasound diagnostics
  99. Research Article
  100. Measurement of needle susceptibility artifacts in magnetic resonance images
  101. Research Article
  102. Dimensionality reduction of medical image descriptors for multimodal image registration
  103. Research Article
  104. Experimental evaluation of different weighting schemes in magnetic particle imaging reconstruction
  105. Research Article
  106. Evaluation of CT capability for the detection of thin bone structures
  107. Research Article
  108. Towards contactless optical coherence elastography with acoustic tissue excitation
  109. Research Article
  110. Development and implementation of algorithms for automatic and robust measurement of the 2D:4D digit ratio using image data
  111. Research Article
  112. Automated high-throughput analysis of B cell spreading on immobilized antibodies with whole slide imaging
  113. Research Article
  114. Tissue segmentation from head MRI: a ground truth validation for feature-enhanced tracking
  115. Research Article
  116. Video tracking of swimming rodents on a reflective water surface
  117. Research Article
  118. MR imaging of model drug distribution in simulated vitreous
  119. Research Article
  120. Studying the extracellular contribution to the double wave vector diffusion-weighted signal
  121. Research Article
  122. Artifacts in field free line magnetic particle imaging in the presence of inhomogeneous and nonlinear magnetic fields
  123. Research Article
  124. Introducing a frequency-tunable magnetic particle spectrometer
  125. Research Article
  126. Imaging of aortic valve dynamics in 4D OCT
  127. Research Article
  128. Intravascular optical coherence tomography (OCT) as an additional tool for the assessment of stent structures
  129. Research Article
  130. Simple concept for a wide-field lensless digital holographic microscope using a laser diode
  131. Research Article
  132. Intraoperative identification of somato-sensory brain areas using optical imaging and standard RGB camera equipment – a feasibility study
  133. Research Article
  134. Respiratory surface motion measurement by Microsoft Kinect
  135. Research Article
  136. Improving image quality in EIT imaging by measurement of thorax excursion
  137. Research Article
  138. A clustering based dual model framework for EIT imaging: first experimental results
  139. Research Article
  140. Three-dimensional anisotropic regularization for limited angle tomography
  141. Research Article
  142. GPU-based real-time generation of large ultrasound volumes from freehand 3D sweeps
  143. Research Article
  144. Experimental computer tomograph
  145. Research Article
  146. US-tracked steered FUS in a respiratory ex vivo ovine liver phantom
  147. Research Article
  148. Contribution of brownian rotation and particle assembly polarisation to the particle response in magnetic particle spectrometry
  149. Research Article
  150. Preliminary investigations of magnetic modulated nanoparticles for microwave breast cancer detection
  151. Research Article
  152. Construction of a device for magnetic separation of superparamagnetic iron oxide nanoparticles
  153. Research Article
  154. An IHE-conform telecooperation platform supporting the treatment of dementia patients
  155. Research Article
  156. Automated respiratory therapy system based on the ARDSNet protocol with systemic perfusion control
  157. Research Article
  158. Identification of surgical instruments using UHF-RFID technology
  159. Research Article
  160. A generic concept for the development of model-guided clinical decision support systems
  161. Research Article
  162. Evaluation of local alterations in femoral bone mineral density measured via quantitative CT
  163. Research Article
  164. Creating 3D gelatin phantoms for experimental evaluation in biomedicine
  165. Research Article
  166. Influence of short-term fixation with mixed formalin or ethanol solution on the mechanical properties of human cortical bone
  167. Research Article
  168. Analysis of the release kinetics of surface-bound proteins via laser-induced fluorescence
  169. Research Article
  170. Tomographic particle image velocimetry of a water-jet for low volume harvesting of fat tissue for regenerative medicine
  171. Research Article
  172. Wireless medical sensors – context, robustness and safety
  173. Research Article
  174. Sequences for real-time magnetic particle imaging
  175. Research Article
  176. Speckle-based off-axis holographic detection for non-contact photoacoustic tomography
  177. Research Article
  178. A machine learning approach for planning valve-sparing aortic root reconstruction
  179. Research Article
  180. An in-ear pulse wave velocity measurement system using heart sounds as time reference
  181. Research Article
  182. Measuring different oxygenation levels in a blood perfusion model simulating the human head using NIRS
  183. Research Article
  184. Multisegmental fusion of the lumbar spine a curse or a blessing?
  185. Research Article
  186. Numerical analysis of the biomechanical complications accompanying the total hip replacement with NANOS-Prosthetic: bone remodelling and prosthesis migration
  187. Research Article
  188. A muscle model for hybrid muscle activation
  189. Research Article
  190. Mathematical, numerical and in-vitro investigation of cooling performance of an intra-carotid catheter for selective brain hypothermia
  191. Research Article
  192. An ideally parameterized unscented Kalman filter for the inverse problem of electrocardiography
  193. Research Article
  194. Interactive visualization of cardiac anatomy and atrial excitation for medical diagnosis and research
  195. Research Article
  196. Virtualizing clinical cases of atrial flutter in a fast marching simulation including conduction velocity and ablation scars
  197. Research Article
  198. Mesh structure-independent modeling of patient-specific atrial fiber orientation
  199. Research Article
  200. Accelerating mono-domain cardiac electrophysiology simulations using OpenCL
  201. Research Article
  202. Understanding the cellular mode of action of vernakalant using a computational model: answers and new questions
  203. Research Article
  204. A java based simulator with user interface to simulate ventilated patients
  205. Research Article
  206. Evaluation of an algorithm to choose between competing models of respiratory mechanics
  207. Research Article
  208. Numerical simulation of low-pulsation gerotor pumps for use in the pharmaceutical industry and in biomedicine
  209. Research Article
  210. Numerical and experimental flow analysis in centifluidic systems for rapid allergy screening tests
  211. Research Article
  212. Biomechanical parameter determination of scaffold-free cartilage constructs (SFCCs) with the hyperelastic material models Yeoh, Ogden and Demiray
  213. Research Article
  214. FPGA controlled artificial vascular system
  215. Research Article
  216. Simulation based investigation of source-detector configurations for non-invasive fetal pulse oximetry
  217. Research Article
  218. Test setup for characterizing the efficacy of embolic protection devices
  219. Research Article
  220. Impact of electrode geometry on force generation during functional electrical stimulation
  221. Research Article
  222. 3D-based visual physical activity assessment of children
  223. Research Article
  224. Realtime assessment of foot orientation by Accelerometers and Gyroscopes
  225. Research Article
  226. Image based reconstruction for cystoscopy
  227. Research Article
  228. Image guided surgery innovation with graduate students - a new lecture format
  229. Research Article
  230. Multichannel FES parameterization for controlling foot motion in paretic gait
  231. Research Article
  232. Smartphone supported upper limb prosthesis
  233. Research Article
  234. Use of quantitative tremor evaluation to enhance target selection during deep brain stimulation surgery for essential tremor
  235. Research Article
  236. Evaluation of adhesion promoters for Parylene C on gold metallization
  237. Research Article
  238. The influence of metallic ions from CoCr28Mo6 on the osteogenic differentiation and cytokine release of human osteoblasts
  239. Research Article
  240. Increasing the visibility of thin NITINOL vascular implants
  241. Research Article
  242. Possible reasons for early artificial bone failure in biomechanical tests of ankle arthrodesis systems
  243. Research Article
  244. Development of a bending test procedure for the characterization of flexible ECoG electrode arrays
  245. Research Article
  246. Tubular manipulators: a new concept for intracochlear positioning of an auditory prosthesis
  247. Research Article
  248. Investigation of the dynamic diameter deformation of vascular stents during fatigue testing with radial loading
  249. Research Article
  250. Electrospun vascular grafts with anti-kinking properties
  251. Research Article
  252. Integration of temperature sensors in polyimide-based thin-film electrode arrays
  253. Research Article
  254. Use cases and usability challenges for head-mounted displays in healthcare
  255. Research Article
  256. Device- and service profiles for integrated or systems based on open standards
  257. Research Article
  258. Risk management for medical devices in research projects
  259. Research Article
  260. Simulation of varying femoral attachment sites of medial patellofemoral ligament using a musculoskeletal multi-body model
  261. Research Article
  262. Does enhancing consciousness for strategic planning processes support the effectiveness of problem-based learning concepts in biomedical education?
  263. Research Article
  264. SPIO processing in macrophages for MPI: The breast cancer MPI-SNLB-concept
  265. Research Article
  266. Numerical simulations of airflow in the human pharynx of OSAHS patients
Downloaded on 25.12.2025 from https://www.degruyterbrill.com/document/doi/10.1515/cdbme-2015-0072/html
Scroll to top button