Home Radiopacity assessment of neurovascular implants
Article Open Access

Radiopacity assessment of neurovascular implants

  • Thomas Hoffmann EMAIL logo , Sebastian Gugel , Oliver Beuing and Georg Rose
Published/Copyright: September 30, 2016

Abstract

An experimental method for radiopacity evaluation of neurovascular implants is presented. State of the art methods are described for cardiovasular implants. To extend the methods to neurovascular devices, imaging parameters were determined from clinical radiography data sets and protocols. To reduce background noise, a thresholding method is introduced. This method is evaluated on a setup with four different laser cut and braided stents. The results are compared to current evaluation methods without filtering. It is shown that unfiltered methods underestimate the radiopacity by 0.89 mm ± 0.21 mm (aluminum equivalent). The introduced method is highly repeatable and uses modern clinical imaging technology. Thus an image evaluation under realisitc conditions is possible. Slightly visible vascular devices and implants can be compared with high accuracy.

1 Introduction

Implants for cerebral arteries, like laser cut stents or braided nitinol and flow diverter stents often show poor radiopacity in radiography. The tradeoff between delivery systems with a small diameter and material thickness as a main parameter influencing the radiopacity is the major reason for poor visibility of neurovascular implants. Furthermore, skull bones which highly absorb X-ray, may impair visibility in the target areas of the implants.

But implant visibility is a crucial factor for therapeutic success. Therefore, implant manufacturers face the problem to develop implants with adequate visibility. For this, a reproducible testing setup with realistic properties is needed to evaluate the progress. Currently no common written standard for the evaluation of X-ray visibility of vascular stents is available. Few details were published by Schmidt et al. [1] referring to DIN 13273-7:2003-08 [2] and ASTM F640-12 [3] as well as by Dobránszky et al. in [4]. Different methods for the evaluation of radiopacity of neurovascular implants were based on qualitative phantom studies or contrast measurements with individual setups [5]. These setups do not comply with the imaging standards of a neurovascular intervention and should be adapted. In the following, a test setup and an image processing method for a quantitative and repeatable evaluation will be introduced and tested on different implants.

1.1 State of the art

State of the art methods recommend the use of a flat detector and X-ray system with a tube voltage of 70 kV and a detector-to-source distance of 1150 mm [1]. A scattering body made of aluminum or polymethyl methacrylate (PMMA) is recommended to be used. As described by Schmidt et al. the scattering phantom mimics the attenuation of the human body and decreases the contrast between stent and background. The resulting image is evaluated by contrast calculations, aluminum equivalent or an image histogram. Michelson or Weber contrast is calculated by means of an area including the stent and a reference area right beside the stent. The stent surface area has a strong influence on this method. With this, contrast calculations of stents with small struts (low stent surface area) and high radiopacity will be underestimated by the arithmetic mean, because the stent structure is not separated from the surrounding.

The aluminum equivalent is determined by measuring the mean gray value of the aluminum steps. A linear regression of those values is done and the calculated mean gray value of the implant area is fitted in the regression line. This results in an equivalent aluminum thickness value, which is independent of the composition of the implant and thus enables the comparison of different stents.

Another method is based on a reference radiographic image without stent and a second image with a stent placed in a defined area. The image without stent is subtracted from that with stent. Consequently it is possible to merely calculate the gray value of the stent surface area by using a simple thresholding. The disadvantage of this method is that the distribution of noise cannot be predicted and varies substantially between different images. Since different material or geometric properties of implants only have little impact on the radiopacity, even minor changes and uncertainties in the measurement setup can influence the result.

The aim of our method is to avoid the inaccuracies and uncertainties of current methods.

2 Material and methods

2.1 Setup

Imaging parameters were determined by the observation of stent interventions. The parameters were adapted to a setup based on a clinical angiographic system (Siemens Artis zeego, Siemens Healthcare GmbH, Erlangen, Germany). The setup comprises a scattering body made of a 200 × 200 × 150 × mm3 (l × w × h) polymethyl methacrylate (PMMA) block (see Figure 1). An aluminum step wedge (see Figure 2 for details) is used as a reference object to translate stent attenuation into a universal and comparable value of equivalent aluminum thickness. The step wedge is placed in the center of scattering body and detector.

Figure 1 The mechanical setup consists of a tissue mimicking scattering body, an aluminum reference step wedge and a flat panel C-arm. System parameters were compared to realistic clinical scenarios.
Figure 1

The mechanical setup consists of a tissue mimicking scattering body, an aluminum reference step wedge and a flat panel C-arm. System parameters were compared to realistic clinical scenarios.

Figure 2 Experimental setup with defined image center (black circle), step wedge, PMMA block (APMMA), area for implant placement (AMA) and different reference areas for contrast calculation. The reference area (ASref) for stent contrast calculation has the same geometry as the stent region (AS) in length (lS) and width (dS).
Figure 2

Experimental setup with defined image center (black circle), step wedge, PMMA block (APMMA), area for implant placement (AMA) and different reference areas for contrast calculation. The reference area (ASref) for stent contrast calculation has the same geometry as the stent region (AS) in length (lS) and width (dS).

Several factors attributed to the X-ray system can strongly influence the comparability of different implants and distort the result. For instance, dead pixels of the detector system can overestimate the radiopacity, in particular of small structures. Likewise, an inhomogeneous X-ray exposure, e.g. due to the heel effect, can influence the measurement (see Figure 3). Those site depending parameters were determined and considered in our calculations.

Figure 3 Experimental setup with inhomogenous dose distribution on the detector caused by the Heel effect of the X-ray tube. 1 - Measuring area with nearly homogenous distribution. 2 - Investigated vertical line with listed gray scale values in diagram below. 3 - Aluminum step wedge.
Figure 3

Experimental setup with inhomogenous dose distribution on the detector caused by the Heel effect of the X-ray tube. 1 - Measuring area with nearly homogenous distribution. 2 - Investigated vertical line with listed gray scale values in diagram below. 3 - Aluminum step wedge.

According to the clinical scenario the detector is placed above the scattering body parallel to the surface of the PMMA block. Source to detector distance is adjusted to 1000 mm. The PMMA block is placed on an extra mounting which is connected to the patient table. The mounting prevents that the patient table is in the field-of-view.

2.2 Image system parameters

State of the art angiographic systems operate with complex control systems, which adapt dose, tube voltage and exposure time during the radiographic acquisition to improve the image quality and reduce the patient dose. As the X-ray parameters depend on the imaged objects, it is impossible to perform identical experiments. The eventual consequent X-ray parameter highly depend on the object imaged. Therefore it is necessary to ensure a reproducible placement of the scattering body and the aluminum reference.

Image post-processing filters (gain correction, noise reduction, edge enhancement) were disabled and the copper filter was switched off to detect the unfiltered stent characteristics without manufacturer depending image properties. The tube voltage was set to 70 kV as used in clinical practice. Radiography is done in single shot mode with a quadratic measuring field and a field of view of 200 × 200 mm2.

2.3 Image evaluation

The Michelson contrast measure is used to assess the visibility of the implant. For this purpose, a reference region close to the stent (with the same length and width as the stent region) is defined. The mean gray value GrefM of this region is calculated by the arithmetic mean method of all pixels. As distinguished from state of the art methods, the background is separated by a threshold filtering from the stent to merely calculate the mean gray value of the implant surface structure. The threshold between background and stent structure is defined as the value for the reduction of the pixel number in the reference area down to 1% for 70 kV. The remaining pixels ensure, that dead pixels have no influence to the result. The mean stent gray value GstentM is used for the calculation of the Michelson contrast C to:

(1)C=GrefMGstentMGrefM+GstentM

As reference for the measurement of the aluminum equivalent, the contrast of the single steps of the aluminum wedge is determined by calculating the mean gray values. The area right beside the 0.5 mm aluminum step is used as reference area (see Figure 2). It is identical in height and width as a single step. The respective contrast of each step is divided by the maximum Michelson contrast of 10 mm aluminum and a linear regression is calculated out of this. For calibration the exact height of the used aluminum step is measured and taken into account. The contrast value of the segmented stent is divided by the maximum contrast of the aluminum step wedge and fit into the regression line.

The majority of neurovascular stents contain marker structures made of a divergent material than the stent struts to enhance the attenuation at well-defined points. The markers were considered separately with the described methods to evaluate the stent radiopacity. For this, the area of interest is reduced to the marker structure.

The amount of pixels of the filtered stent, represented in the area of interest, is counted to obtain an estimation of the projected stent surface area. The coverage of the stent is calculated by dividing the pixel number of corpus and marker by the number of pixel in the unfiltered stent region.

2.4 Quality assurance

For the approval processes of medical devices all relevant parameters of a testing should be recorded. For this the following is necessary:

  • calibration protocol of the aluminium step wedge

  • tested detector (without dead pixels)

  • predefined area for stent positioning without disruptive tube/detector effects

  • logfile of the X-ray system

  • original DICOM-data

  • filtered DICOM files

  • diagram of contrast and aluminum equivalent with fitted stent/marker

  • number of structure representing pixels

A software demonstrator was created to perform the calibration, filtering, measurement, calculation, evaluation and data export as described above.

3 Results

Body and marker structures of four different nitinol stents were investigated with the introduced method with and without the thresholding method. One was a laser cut stent (3.5 × 25 mm2, Figure 4, 1A–C), the other three were braided nitinol stents (3.5 × 15 mm2, see Figure 4, 2A–C; 6.0 × 35 mm2, see Figure 4, 3A–C; 4.0 × 33 mm2, see Figure 4, 4A–C). Stent three and four had the largest stent areas, which find expression in the number of pixels. In relation to the size of the corresponding stent area, the relative coverage of stent four is the highest with 73.1%. The aluminum equivalent was calculated for every corpus and one marker structure of each stent (see Table 1).

Figure 4 Investigated Stents. Radiography (A) with the previous described parameters. Filtered stent corpus (B) and marker structures (C) from radiography images. Implant size is not comparable, due to a different scale.
Figure 4

Investigated Stents. Radiography (A) with the previous described parameters. Filtered stent corpus (B) and marker structures (C) from radiography images. Implant size is not comparable, due to a different scale.

Table 1

Evaluation of the investigated stents. Aluminum equivalent was measured for stent one to stent four for corpus (B) and marker (M) structures. The new threshold based calculation was compared to the method without filter. Pixel number and relative coverage illustrate the surface areas and geometries of each stent corpus and marker.

ImplantAl equivalent [mm]Pixel [−]Rel. coverage [%]
w/o filterw filter
B10.000.981911.54
B20.211.04306544.00
B30.301.04923532.89
B40.871.49924173.10
M10.001.11450.36
M20.621.372834.06
M30.021.271880.66
M40.601.45930.74

It was shown that unfiltered images underestimate the equivalent aluminum thickness up to 0.89 mm ± 0.21 mm.

4 Conclusion

In this paper we presented a method for radiopacity assessment of neurovascular implants, which is less prone to inaccuracies than currently used methods. The results are more precise with background filtering. By counting the pixels, representing the stent corpus and marker, information about the object size can be derived. Furthermore, the introduced method employs imaging parameters which are often used in neurointerventional practice. Consequently the method generates realistic data sets with a highly reproducible evaluation and documentation process.

Author’s Statement

Research funding: This work was partly funded by the Federal Ministry of Education and Research in Germany within the Research Campus STIMULATE under grant number 03FO16102A. 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 animal use.

References

[1] Schmidt W, Behrens P, Kamke F, Schmitz KP. Measurement of radiopacity of vascular implants. Biomed Tech (Berl). 2014;59:S1216–9.Search in Google Scholar

[2] DIN 13273-7:2003-08 – Catheters for medical use – Part 7: Determination of the x-ray attenuation of catheters; Requirements and testing.Search in Google Scholar

[3] ASTM F640 – 12 – Standard test methods for determining radiopacity for medical use.Search in Google Scholar

[4] Dobránszky J, Ring G, Bognár E, Kovács R, Bitay E. New method for evaluating the visibility of coronary stents. Acta Polytechnica Hungarica. 2014;11;81–94.10.12700/APH.11.05.2014.05.5Search in Google Scholar

[5] Boese A, Rose G, Friebe M, Hoffmann T, Serowy S, Skalej M, et al. Increasing the visibility of thin NITINOL vascular implants. Current Directions in Biomedical Engineering. 2015;1:503–6.10.1515/cdbme-2015-0120Search in Google Scholar

Published Online: 2016-9-30
Published in Print: 2016-9-1

©2016 Thomas Hoffmann et al., licensee De Gruyter.

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.

Articles in the same Issue

  1. Synthesis and characterization of PIL/pNIPAAm hybrid hydrogels
  2. Novel blood protein based scaffolds for cardiovascular tissue engineering
  3. Cell adhesion and viability of human endothelial cells on electrospun polymer scaffolds
  4. Effects of heat treatment and welding process on superelastic behaviour and microstructure of micro electron beam welded NiTi
  5. Long-term stable modifications of silicone elastomer for improved hemocompatibility
  6. The effect of thermal treatment on the mechanical properties of PLLA tubular specimens
  7. Biocompatible wear-resistant thick ceramic coating
  8. Protection of active implant electronics with organosilicon open air plasma coating for plastic overmolding
  9. Examination of dielectric strength of thin Parylene C films under various conditions
  10. Open air plasma deposited antimicrobial SiOx/TiOx composite films for biomedical applications
  11. Systemic analysis about residual chloroform in PLLA films
  12. A macrophage model of osseointegration
  13. Towards in silico prognosis using big data
  14. Technical concept and evaluation of a novel shoulder simulator with adaptive muscle force generation and free motion
  15. Usability evaluation of a locomotor therapy device considering different strategies
  16. Hypoxia-on-a-chip
  17. Integration of a semi-automatic in-vitro RFA procedure into an experimental setup
  18. Fabrication of MEMS-based 3D-μECoG-MEAs
  19. High speed digital interfacing for a neural data acquisition system
  20. Bionic forceps for the handling of sensitive tissue
  21. Experimental studies on 3D printing of barium titanate ceramics for medical applications
  22. Patient specific root-analogue dental implants – additive manufacturing and finite element analysis
  23. 3D printing – a key technology for tailored biomedical cell culture lab ware
  24. 3D printing of hydrogels in a temperature controlled environment with high spatial resolution
  25. Biocompatibility of photopolymers for additive manufacturing
  26. Biochemical piezoresistive sensors based on pH- and glucose-sensitive hydrogels for medical applications
  27. Novel wireless measurement system of pressure dedicated to in vivo studies
  28. Portable auricular device for real-time swallow and chew detection
  29. Detection of miRNA using a surface plasmon resonance biosensor and antibody amplification
  30. Simulation and evaluation of stimulation scenarios for targeted vestibular nerve excitation
  31. Deep brain stimulation: increasing efficiency by alternative waveforms
  32. Prediction of immediately occurring microsleep events from brain electric signals
  33. Determining cardiac vagal threshold from short term heart rate complexity
  34. Classification of cardiac excitation patterns during atrial fibrillation
  35. An algorithm to automatically determine the cycle length coverage to identify rotational activity during atrial fibrillation – a simulation study
  36. Deriving respiration from high resolution 12-channel-ECG during cycling exercise
  37. Reducing of gradient induced artifacts on the ECG signal during MRI examinations using Wilcoxon filter
  38. Automatic detection and mapping of double potentials in intracardiac electrograms
  39. Modeling the pelvic region for non-invasive pelvic intraoperative neuromonitoring
  40. Postprocessing algorithm for automated analysis of pelvic intraoperative neuromonitoring signals
  41. Best practice: surgeon driven application in pelvic operations
  42. Vasomotor assessment by camera-based photoplethysmography
  43. Classification of morphologic changes in photoplethysmographic waveforms
  44. Novel computation of pulse transit time from multi-channel PPG signals by wavelet transform
  45. Efficient design of FIR filter based low-pass differentiators for biomedical signal processing
  46. Nonlinear causal influences assessed by mutual compression entropy
  47. Comparative study of methods for solving the correspondence problem in EMD applications
  48. fNIRS for future use in auditory diagnostics
  49. Semi-automated detection of fractional shortening in zebrafish embryo heart videos
  50. Blood pressure measurement on the cheek
  51. Derivation of the respiratory rate from directly and indirectly measured respiratory signals using autocorrelation
  52. Left cardiac atrioventricular delay and inter-ventricular delay in cardiac resynchronization therapy responder and non-responder
  53. An automatic systolic peak detector of blood pressure waveforms using 4th order cumulants
  54. Real-time QRS detection using integrated variance for ECG gated cardiac MRI
  55. Preprocessing of unipolar signals acquired by a novel intracardiac mapping system
  56. In-vitro experiments to characterize ventricular electromechanics
  57. Continuous non-invasive monitoring of blood pressure in the operating room: a cuffless optical technology at the fingertip
  58. Application of microwave sensor technology in cardiovascular disease for plaque detection
  59. Artificial blood circulatory and special Ultrasound Doppler probes for detecting and sizing gaseous embolism
  60. Detection of microsleep events in a car driving simulation study using electrocardiographic features
  61. A method to determine the kink resistance of stents and stent delivery systems according to international standards
  62. Comparison of stented bifurcation and straight vessel 3D-simulation with a prior simulated velocity profile inlet
  63. Transient Euler-Lagrange/DEM simulation of stent thrombosis
  64. Automated control of the laser welding process of heart valve scaffolds
  65. Automation of a test bench for accessing the bendability of electrospun vascular grafts
  66. Influence of storage conditions on the release of growth factors in platelet-rich blood derivatives
  67. Cryopreservation of cells using defined serum-free cryoprotective agents
  68. New bioreactor vessel for tissue engineering of human nasal septal chondrocytes
  69. Determination of the membrane hydraulic permeability of MSCs
  70. Climate retainment in carbon dioxide incubators
  71. Multiple factors influencing OR ventilation system effectiveness
  72. Evaluation of an app-based stress protocol
  73. Medication process in Styrian hospitals
  74. Control tower to surgical theater
  75. Development of a skull phantom for the assessment of implant X-ray visibility
  76. Surgical navigation with QR codes
  77. Investigation of the pressure gradient of embolic protection devices
  78. Computer assistance in femoral derotation osteotomy: a bottom-up approach
  79. Automatic depth scanning system for 3D infrared thermography
  80. A service for monitoring the quality of intraoperative cone beam CT images
  81. Resectoscope with an easy to use twist mechanism for improved handling
  82. In vitro simulation of distribution processes following intramuscular injection
  83. Adjusting inkjet printhead parameters to deposit drugs into micro-sized reservoirs
  84. A flexible standalone system with integrated sensor feedback for multi-pad electrode FES of the hand
  85. Smart control for functional electrical stimulation with optimal pulse intensity
  86. Tactile display on the remaining hand for unilateral hand amputees
  87. Effects of sustained electrical stimulation on spasticity assessed by the pendulum test
  88. An improved tracking framework for ultrasound probe localization in image-guided radiosurgery
  89. Improvement of a subviral particle tracker by the use of a LAP-Kalman-algorithm
  90. Learning discriminative classification models for grading anal intraepithelial neoplasia
  91. Regularization of EIT reconstruction based on multi-scales wavelet transforms
  92. Assessing MRI susceptibility artefact through an indicator of image distortion
  93. EyeGuidance – a computer controlled system to guide eye movements
  94. A framework for feedback-based segmentation of 3D image stacks
  95. Doppler optical coherence tomography as a promising tool for detecting fluid in the human middle ear
  96. 3D Local in vivo Environment (LivE) imaging for single cell protein analysis of bone tissue
  97. Inside-Out access strategy using new trans-vascular catheter approach
  98. US/MRI fusion with new optical tracking and marker approach for interventional procedures inside the MRI suite
  99. Impact of different registration methods in MEG source analysis
  100. 3D segmentation of thyroid ultrasound images using active contours
  101. Designing a compact MRI motion phantom
  102. Cerebral cortex classification by conditional random fields applied to intraoperative thermal imaging
  103. Classification of indirect immunofluorescence images using thresholded local binary count features
  104. Analysis of muscle fatigue conditions using time-frequency images and GLCM features
  105. Numerical evaluation of image parameters of ETR-1
  106. Fabrication of a compliant phantom of the human aortic arch for use in Particle Image Velocimetry (PIV) experimentation
  107. Effect of the number of electrodes on the reconstructed lung shape in electrical impedance tomography
  108. Hardware dependencies of GPU-accelerated beamformer performances for microwave breast cancer detection
  109. Computer assisted assessment of progressing osteoradionecrosis of the jaw for clinical diagnosis and treatment
  110. Evaluation of reconstruction parameters of electrical impedance tomography on aorta detection during saline bolus injection
  111. Evaluation of open-source software for the lung segmentation
  112. Automatic determination of lung features of CF patients in CT scans
  113. Image analysis of self-organized multicellular patterns
  114. Effect of key parameters on synthesis of superparamagnetic nanoparticles (SPIONs)
  115. Radiopacity assessment of neurovascular implants
  116. Development of a desiccant based dielectric for monitoring humidity conditions in miniaturized hermetic implantable packages
  117. Development of an artifact-free aneurysm clip
  118. Enhancing the regeneration of bone defects by alkalizing the peri-implant zone – an in vitro approach
  119. Rapid prototyping of replica knee implants for in vitro testing
  120. Protecting ultra- and hyperhydrophilic implant surfaces in dry state from loss of wettability
  121. Advanced wettability analysis of implant surfaces
  122. Patient-specific hip prostheses designed by surgeons
  123. Plasma treatment on novel carbon fiber reinforced PEEK cages to enhance bioactivity
  124. Wear of a total intervertebral disc prosthesis
  125. Digital health and digital biomarkers – enabling value chains on health data
  126. Usability in the lifecycle of medical software development
  127. Influence of different test gases in a non-destructive 100% quality control system for medical devices
  128. Device development guided by user satisfaction survey on auricular vagus nerve stimulation
  129. Empirical assessment of the time course of innovation in biomedical engineering: first results of a comparative approach
  130. Effect of left atrial hypertrophy on P-wave morphology in a computational model
  131. Simulation of intracardiac electrograms around acute ablation lesions
  132. Parametrization of activation based cardiac electrophysiology models using bidomain model simulations
  133. Assessment of nasal resistance using computational fluid dynamics
  134. Resistance in a non-linear autoregressive model of pulmonary mechanics
  135. Inspiratory and expiratory elastance in a non-linear autoregressive model of pulmonary mechanics
  136. Determination of regional lung function in cystic fibrosis using electrical impedance tomography
  137. Development of parietal bone surrogates for parietal graft lift training
  138. Numerical simulation of mechanically stimulated bone remodelling
  139. Conversion of engineering stresses to Cauchy stresses in tensile and compression tests of thermoplastic polymers
  140. Numerical examinations of simplified spondylodesis models concerning energy absorption in magnetic resonance imaging
  141. Principle study on the signal connection at transabdominal fetal pulse oximetry
  142. Influence of Siluron® insertion on model drug distribution in the simulated vitreous body
  143. Evaluating different approaches to identify a three parameter gas exchange model
  144. Effects of fibrosis on the extracellular potential based on 3D reconstructions from histological sections of heart tissue
  145. From imaging to hemodynamics – how reconstruction kernels influence the blood flow predictions in intracranial aneurysms
  146. Flow optimised design of a novel point-of-care diagnostic device for the detection of disease specific biomarkers
  147. Improved FPGA controlled artificial vascular system for plethysmographic measurements
  148. Minimally spaced electrode positions for multi-functional chest sensors: ECG and respiratory signal estimation
  149. Automated detection of alveolar arches for nasoalveolar molding in cleft lip and palate treatment
  150. Control scheme selection in human-machine- interfaces by analysis of activity signals
  151. Event-based sampling for reducing communication load in realtime human motion analysis by wireless inertial sensor networks
  152. Automatic pairing of inertial sensors to lower limb segments – a plug-and-play approach
  153. Contactless respiratory monitoring system for magnetic resonance imaging applications using a laser range sensor
  154. Interactive monitoring system for visual respiratory biofeedback
  155. Development of a low-cost senor based aid for visually impaired people
  156. Patient assistive system for the shoulder joint
  157. A passive beating heart setup for interventional cardiology training
Downloaded on 11.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/cdbme-2016-0118/html
Scroll to top button