SKIN DOSE MAPPING IN INTERVENTIONAL CARDIOLOGY - A PRACTICAL SOLUTION (Radiation Protection Dosimetry) (2020)
SKIN DOSE MAPPING IN INTERVENTIONAL CARDIOLOGY - A PRACTICAL SOLUTION (Radiation Protection Dosimetry) (2020)
SKIN DOSE MAPPING IN INTERVENTIONAL CARDIOLOGY - A PRACTICAL SOLUTION (Radiation Protection Dosimetry) (2020)
1093/rpd/ncaa002
Received 11 July 2019; revised 25 November 2019; editorial decision 2 January 2020; accepted 6 January 2020
Numerous cases of radiation-induced tissue reactions following interventional cardiology (IC) procedures have been reported,
resulting in the need for an optimized and personalized dosimetry. At present, there are many fluoroscopy units without Digital
Imaging and Communications in Medicine (DICOM) Radiation Dose Structured Report globally installed. Many of these have
not been updated yet, and may never be, therefore, the main objectives of this paper are to develop an offline skin dose mapping
application, which uses DICOM headers for the peak skin dose (PSD) assessment and to compare the PSD assessment results
to XR-RV3 Gafchromic film for common IC procedures. The mean deviation between the measured and the calculated PSD
was 8.7 ± 26.3%. Simulated skin dose map showed good matching with XR-RV3 Gafchromic film. The skin dose mapping
application presented in this paper is an elegant solution and a suitable alternative to XR-RV3 Gafchromic film.
© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For Permissions, please email: [email protected]
M. KRAJINOVIĆ ET AL.
2
SKIN DOSE MAPPING IN INTERVENTIONAL CARDIOLOGY
using the position of the X-ray tube focal spot and the Application validation
distance from the source to detector. In order to find
XR-RV3 Gafchromic films calibration was per-
the affected patient area, it is necessary to find four
formed in the Secondary Standard Dosimetry
corners of detector. The affected skin location by X-
Laboratory of Vinča Institute of Nuclear Sciences
ray field is then found by finding intersection points
between the patient plane and the lines going through (VINČA). Calibration setup was done as described
X-ray tube focal spot and four corners of detector. in the paper of Farah et al.(23) , where comprehensive
Once four points in the patient plane are determined, characterization of XR-RV3 Gafchromic films was
a polygon can be defined, and thus the affected skin conducted. Calibration was carried out free-in-
area is found. air, where films were placed on a 0.5 cm-thick
PMMA support providing negligible backscatter and
irradiated at 1 m from the focal spot. Films from the
Skin dose calculation same batch were cut in small pieces of 2 x 2 cm and
After the affected skin location is identified, dose is aligned to the central flat region of the X-ray field.
calculated using Equation 1, where Ka,r is air kerma All film pieces were oriented so that the yellow side
at the IRP corresponding to fluoroscopic and cine of the film was in direction of the focal spot and dose
acquisition of each irradiation, CF is calibration fac- delivery was monitored using an ionization chamber.
tor for Ka,r determined annually from measurements The dose interval used to form the calibration curve
made during the annual quality control tests, dIRP is ranged from 0 to 5 Gy. Scanning was done after 24 h
the distance form the source to the Interventional of exposure. The Hewlett-Packard Scanjet G3110
Reference Point (60 cm in this case), dpatient is the was used for scanning the films. Color images were
distance from the source to each pixel of the affected acquired using 16 bits per channel (RGB images) with
area, BSF is backscatter factor interpolated from a resolution of 600 PPI in TIF format. Films were
the coefficients of Benmakhlouf et al.(21) , MAEC analyzed in the Python programming language using
is the ratio of mass energy absorption coefficients only the red channel of the scanned image, as this
for skin to air interpolated from the NIST Physical ensured maximum sensitivity. In order to calculate the
Reference Data Library (skin is modeled as ICRU PSD, a square region of interest (ROI) was formed
four component)(22) : with sides of 225 px (ROI was ∼1 cm2 ). The ROI
thus formed was shifted along the image, and then
the mean value of the red channel inside ROI was
2 calculated, resulting in the mean pixel value (MPV)
dIRP
Dose = Ka,r × CF × × BSF × MAEC. of the red image component. Two unirradiated film
dpatient pieces were used in order to obtain the unexposed
(1) MPV. The calibration curve was obtained using
polynomial
fitting based on reflectance, defined as
log MPVunexposed /MPVexposed , because it showed
Beam attenuation due to the presence of table
and table pad is included in CF. Since the fluoro- the best performance(23) .
scopic acquisitions are not recorded (only cine acqui- During cardiac procedures XR-RV3 Gafchromic
sitions are recorded), it is assumed that fluoroscopic films were placed below the patient in the region of
kerma is equally spread on each projection used for the upper part of the patient’s torso, with the longer
cine acquisition. Fluoroscopic kerma is calculated by side of the film along the cranial-caudal axis, and the
subtracting the total cine Ka,r from the total Ka,r , yellow side of the film oriented in direction of the
which can be found in a proprietary dose report or X-ray tube. PSD was obtained using the same Python
recorded manually after the procedure. routine used in film calibration. The ROI with the
3
M. KRAJINOVIĆ ET AL.
1 67 167 103
2 59 156 60
3 55 185 110
4 55 180 96
5 74 186 90
6 47 174 85
7 65 175 87
8 78 183 89
9 53 168 85
10 53 186 99
Mean ± Std (minimum–maximum) 60 ± 10 (47–78) 176 ± 10 (156–186) 90 ± 13 (60–110)
lowest MPV value was used in order to find the PSD. • Entering cumulative dose, i.e. cumulative air
It should be pointed out that the film does not disturb kerma at the IRP (mandatory to enter in order
interventional procedure in any way, nor does it create to obtain skin dose map).
additional discomfort to the patient. • Displaying skin dose map.
The application generates a skin dose map fast
Geometrical comparison between XR-RV3 (e.g. it takes a few seconds to generate a skin dose map
Gafchromic film and simulation from 25 DICOM files). The skin dose map is gener-
ated in the right side of GUI every time the ‘SKIN
During the coding process geometrical verification DOSE MAP’ button is pressed. A colorbar is shown
was done in order to validate simulated skin dose right to the patient skin dose map for fast visual
map. XR-RV3 Gafchromic films were scanned as inspection of skin dose distribution. The calculated
8 bits per channel color images and were converted PSD is shown in the left corner of the skin dose map.
to grayscale images. Exposure, contrast and levels
were modified in order to enhance visual difference
between the exposed and the unexposed part of the XR-RV3 Gafchromic films calibration
scanned images. The measurement scale is shown on XR-RV3 Gafchromic calibration curve for determin-
all images for better comparison. ing skin dose is shown in Figure 3. Fitting equation is
as follows:
RESULTS
Patient population and fluoroscopy Ka,r Fit = p1 × Ref 3 + p2 × Ref 2 + p3 × Ref + p4 , (2)
Overall, ten PCI procedures (eight men and two
women) were analyzed in this paper. The mean age, where p1 , p2 , p3 , p4 are coefficients of a third-degree
height and weight of patients are shown in Table 2. polynomial and Ref stands for reflectance. Coeffi-
Average contribution of fluoroscopic Ka,r to total cients have the following values: p1 = 5.234 × 105 , p2
Ka,r was 57.2%, whereas average contribution of cine = -7.841 × 104 , p3 = 1.358 × 104 , p4 = -24.88.
Ka,r to total Ka,r was 42.8%. The mean total Ka,r , mean
fluoroscopic Ka,r , mean cine Ka,r and total PKA of
procedures are shown in Table 3. Geometrical comparison between XR-RV3
Gafchromic film and simulation
Skin dose mapping application Geometrical comparison between XR-RV3 Gafchro-
mic films and simulation for two different procedures
The GUI of developed application is shown in
is shown in Figure 4. Clearly there is resemblance
Figure 2. This application enables:
which proves that geometrical solution is accurate. It
• Choosing DICOM files to be used for generating should be noted that the scanner used for geometrical
skin dose map. comparison was not used for calibration purposes
• Entering patient height and weight for patient- due to the presence of horizontal lines on scanned
specific skin dose mapping. images.
4
SKIN DOSE MAPPING IN INTERVENTIONAL CARDIOLOGY
Patient Total Ka,r (mGy) Fluoroscopic Ka,r Cine Ka,r (mGy) Total PKA
(mGy) (μGym2 )
Figure 4. Geometrical comparison between XR-RV3 Gafchromic films (left) and simulations (right).
6
SKIN DOSE MAPPING IN INTERVENTIONAL CARDIOLOGY
Gafchromic films are used as the gold standard in polymethyl methacrylate (PMMA) phantom using
this study; however, it should be noted that skin Siemens Artis Zee and GE Innova IGS interventional
dose assessment expanded uncertainty using XR-RV3 systems, respectively.
Gafchromic films is 9% (k = 2) for tightly controlled The recognized limitations of the solution pre-
measurement conditions, adequate laboratory cali- sented in this paper are explained in the following
brations and well-defined readout protocol. A more text. Firstly, DICOM header does not contain the
realistic expanded uncertainty is 41% (k = 2), where data regarding the table movement in a 3D space,
a well-defined laboratory calibration is performed, therefore several assumptions are made, including: (1)
while other influencing parameters related to clinical the table does not move in the lateral and longitudinal
application of dosimetry films are less controlled(23) . direction, and (2) the table height is constant, so the
It should also be pointed out that the calculated skin table does not move in vertical direction. Lack of this
dose has non-negligible uncertainty. Only calibration information contributes to the uncertainty of a skin
factor, CF, contributes to expanded uncertainty dose map. DICOM RDSR, however, contains all the
with 17% (k = 2) according to annual quality necessary information regarding the table movement
control tests; therefore, it is safe to assume that skin in 3D space(5) . Secondly, the 3D modeling of the
dose assessment expanded uncertainty using this patient would inevitably represent a more realistic
application is beyond 20% (k = 2). dose assessment, especially for larger angulations and
Since there are no geometric and dosimetric data oblique projections. The size of the patient as well
on each fluoroscopic exposure, it is assumed that flu- as the exact position of the patient on the table can
oroscopic kerma is equally spread on each cine acqui- additionally make the dose assessment more precise.
sition, but in reality, it could be different. DICOM Finally, skin dose map is generated after the proce-
RDSR is a solution to these problems, since it keeps dure; however, it would be more practical to have an
data for each radiation exposure. online skin dose mapping solution.
To the best of our knowledge, only three 2D
applications verified with Gafchromic XR-RV3 films
CONCLUSION
in clinical situations were developed. Greffier et al.(18)
accomplished average difference of 3.4 ± 21.1% This paper presents an originally developed IC
between Gafchromic XR-RV3 films and the appli- skin dose mapping application for many legacy
cation when fluoroscopic kerma was equally spread fluoroscopy units without DICOM RDSR globally
on each cine acquisition. Bordier et al.(9) compared a installed which have not been updated yet, and may
calculation method that produces dose maps using never be. The application uses standardized DICOM
Gafchromic XR-RV3 films on a water phantom headers to generate a skin dose map and represents
and a body phantom and found that dose values a suitable alternative to XR-RV3 Gafchromic film in
agreed within 24.9% accuracy. Habib Geryes et al.(7) monitoring delivered dose during IC procedures. The
accomplished average difference of 10 ± 7% and goal of future research is implementation of DICOM
9 ± 7% between the calculated and the measured RDSR and development of skin dose maps for 3D
PSD values for 34 test conditions performed on human models.
7
M. KRAJINOVIĆ ET AL.