PURPOSE: This study investigates the distortion of geometry of catheters and anatomy in acquired U/S images, caused by utilizing various stand-off materials for covering a transrectal bi-planar ultrasound probe in HDR and LDR prostate brachytherapy, biopsy and other interventional procedures. Furthermore, an evaluation of currently established water-bath based quality assurance (QA) procedures is presented. MATERIAL AND METHODS: Image acquisitions of an ultrasound QA setup were carried out at 5 MHz and 7 MHz. The U/S probe was covered by EA 4015 Silicone Standoff kit, or UA0059 Endocavity balloon filled either with water or one of the following: 40 ml of Endosgel(®), Instillagel(®), Ultraschall gel or Space OAR™ gel. The differences between images were recorded. Consequently, the dosimetric impact of the observed image distortion was investigated, using a tissue equivalent ultrasound prostate phantom - Model number 053 (CIRS Inc., Norfolk, VA, USA). RESULTS: By using the EA 4015 Silicone Standoff kit in normal water with sound speed of 1525 m/s, a 3 mm needle shift was observed. The expansion of objects appeared in radial direction. The shift deforms also the PTV (prostate in our case) and other organs at risk (OARs) in the same way leading to overestimation of volume and underestimation of the dose. On the other hand, Instillagel(®) and Space OAR™ "shrinks" objects in an ultrasound image for 0.65 mm and 0.40 mm, respectively. CONCLUSIONS: The use of EA 4015 Silicone Standoff kit for image acquisition, leads to erroneous contouring of PTV and OARs and reconstruction and placement of catheters, which results to incorrect dose calculation during prostate brachytherapy. Moreover, the reliability of QA procedures lies mostly in the right temperature of the water used for accurate simulation of real conditions of transrectal ultrasound imaging.
PURPOSE: This study investigates the distortion of geometry of catheters and anatomy in acquired U/S images, caused by utilizing various stand-off materials for covering a transrectal bi-planar ultrasound probe in HDR and LDR prostate brachytherapy, biopsy and other interventional procedures. Furthermore, an evaluation of currently established water-bath based quality assurance (QA) procedures is presented. MATERIAL AND METHODS: Image acquisitions of an ultrasound QA setup were carried out at 5 MHz and 7 MHz. The U/S probe was covered by EA 4015 Silicone Standoff kit, or UA0059 Endocavity balloon filled either with water or one of the following: 40 ml of Endosgel(®), Instillagel(®), Ultraschall gel or Space OAR™ gel. The differences between images were recorded. Consequently, the dosimetric impact of the observed image distortion was investigated, using a tissue equivalent ultrasound prostate phantom - Model number 053 (CIRS Inc., Norfolk, VA, USA). RESULTS: By using the EA 4015 Silicone Standoff kit in normal water with sound speed of 1525 m/s, a 3 mm needle shift was observed. The expansion of objects appeared in radial direction. The shift deforms also the PTV (prostate in our case) and other organs at risk (OARs) in the same way leading to overestimation of volume and underestimation of the dose. On the other hand, Instillagel(®) and Space OAR™ "shrinks" objects in an ultrasound image for 0.65 mm and 0.40 mm, respectively. CONCLUSIONS: The use of EA 4015 Silicone Standoff kit for image acquisition, leads to erroneous contouring of PTV and OARs and reconstruction and placement of catheters, which results to incorrect dose calculation during prostate brachytherapy. Moreover, the reliability of QA procedures lies mostly in the right temperature of the water used for accurate simulation of real conditions of transrectal ultrasound imaging.
In HDR and LDR prostate brachytherapy, needle implantation is guided by real-time ultrasound (U/S) images acquired by using a transrectal bi-planar probe. This probe is covered by a UA0059 Endocavity balloon (CIVCO, Kalona, Iowa, USA: see Fig. 1A) filled with water. The purpose of this balloon is double. Firstly, to provide smoother transmission (coupling) of ultrasound signals between ultrasound probe and human tissue and thus to improve image quality. Secondly, to move the prostate gland according to template coordinates, by adding or removing the proper amount of water in it. Although being successfully used in clinical routine, this equipment is consumable, and after one use it has to be replaced, increasing the cost of every treatment session. On the other hand, EA 4015 Silicone Standoff kit (CIVCO, Kalona, Iowa, USA: Fig. 1B) seems to be a potential solution to the above mentioned disadvantage, as it can be sterilized or used with condoms and reused. The question arise if EA 4015 Silicone Standoff kit can be utilized instead of UA0059 Endocavity Balloon, without distorting imaging of patient's anatomy and needle positioning. This appears as silicone rubber has different acoustic properties than those of water. Any kind of geometry transformation during image acquisition could possibly lead to wrong contouring, reconstruction and placement of needles, and further, to incorrect dose calculation [1-10].
Fig. 1
Two different ways of covering the U/S probe during the image acquisition: (A) UA0059 Endocavity balloon and (B) EA 4015 Silicone Standoff kit (both by CIVCO, Kalona, Iowa, USA)
Two different ways of covering the U/S probe during the image acquisition: (A) UA0059 Endocavity balloon and (B) EA 4015 Silicone Standoff kit (both by CIVCO, Kalona, Iowa, USA)The second part of this study investigates the impact of other hydro-gel type materials on quality of ultrasound imaging. U/S-imaging algorithms are based on the speed of sound in soft tissue [11, 12]. However, established quality assurance procedures for adjustment of U/S-image geometry and template position are based on a water-bath experimental setup. Therefore, in addition, the purpose of this research is to evaluate the accuracy of water-bath based quality assurance procedures.
Material and methods
Theoretical calculations
Speed of sound in normal water
In this experiment, normal or tapwater (not distilled) was used, which contained some impurities (ion concentration). In contrast to pure water, in which sound velocity can be defined only as a function of temperature, the speed of sound in normal water depends on temperature, salinity, and pressure [13-24]. For the theoretical calculations and the experimental procedure a fixed temperature of 37.7°C was selected, thus simulating rectum conditions. For the calculation of water salinity (concentration of mineral salts dissolved in the water or total dissolved salts: TDS), three measurements of electrical conductivity (EC) were performed. They were made by a GMH 3410 operating manual conductivity measuring instrument (Greisinger electronic GmbH, Regenstauf, Germany®) in a sample of the utilized water, at 25°C (reference temperature). The mean value was used to calculate TDS from equation 1 (see Table 1):
Table 1
Physical properties of the utilized tap water: temperature T, electrical conductivity (EC), total dissolved salts (salinity: TDS), water density r and pressure P. Speed of sound in water as a function of temperature, pressure and salinity was derived from Chen-Millero equation [14] which is valid in a range of temperature between 0 to 40°C, salinity of 0 to 40 parts per thousand and pressure of 0 to 1000 bar
T (°C)
EC (µSm/cm)
S (‰)
Density (g/l)
Pressure (bar)
Sound speed (m/s)
37.7
502
0.321
0.993
0.101
1525.421
Physical properties of the utilized tapwater: temperature T, electrical conductivity (EC), total dissolved salts (salinity: TDS), water density r and pressure P. Speed of sound in water as a function of temperature, pressure and salinity was derived from Chen-Millero equation [14] which is valid in a range of temperature between 0 to 40°C, salinity of 0 to 40 parts per thousand and pressure of 0 to 1000 barHydrostatic pressure was calculated using the Equation 2:where P
is the atmospheric pressure, r is the density of the water, g is the acceleration of gravity and h is the depth in water where the sound velocity needs to be calculated. Density of the water r was calculated using the McCutcheon-Martin equation [13], see Table 1. The speed of sound in water as a function of temperature, pressure and salinity was derived from Chen-Millero equation [14], see Table 1. It is valid in a range that meets our requirements: temperature 0 to 40°C, salinity 0 to 40 parts per thousand, pressure 0 to 1000 bar.
Speed of sound in distilled water
Speed of sound in distilled water at 37.7°C was calculated by the experimental equations found in bibliography [13-24].
Expected image distortion with the use of EA 4015 Silicone Standoff kit
To predict the image distortion, we need additionally the speed of sound in silicone rubber. For the determination of this speed, the results from Hachiya et al. [15] were employed. According to their experimental curve, speed of sound in silicone rubber at 37.7°C corresponds to 981.2 m/s.
Potential calibration errors
Any deviation of speed of sound in water in test tank from the U/S calibration speed of 1540 m/s would lead to inaccurate quality assurance tests. Therefore, an estimation of potential calibration errors was done, assuming that percentage acoustic velocity deviation of each water quality from 1540 m/s have the same percentage radial effect in pulse-echo distance. Extended calculations of theoretical needle displacement are presented in the appendix.
Experimental procedure
Setup
On a High Resolution TPV 061 template (GfM-Medizintechnik mbH, Darmstadt, Germany®), 9 stainless steel trocar needles of 1.5 mm×200 mm (Nucletron B.V., Veenendaal, Netherlands) were placed at the positions that are shown in Fig. 2. The choice of the holes was made in such a way to investigate the “shift effect” in ultrasound image, caused by the materials under test, at the central area and edge positions of the template symmetrically. The distances between the holes were large enough, ensuring no overlapping of the needles caused by shadowing artefacts at the acquired ultrasound images.
Fig. 2
Needle positions (shown with a red circles) on a High Resolution TPV 061 template (GfM–Medizintechnik mbH, Darmstadt, Germany). This setup with 9 stainless steel trocar needles of 1.5 mm×200 mm (Nucletron B.V., Veenendaal, Netherlands) was used to investigate the “shift effect”
Needle positions (shown with a red circles) on a High Resolution TPV 061 template (GfM–Medizintechnik mbH, Darmstadt, Germany). This setup with 9 stainless steel trocar needles of 1.5 mm×200 mm (Nucletron B.V., Veenendaal, Netherlands) was used to investigate the “shift effect”The TPV 061 template was connected to a GfM mechanical stepper (GfM-Medizintechnik mbH, Darmstadt, Germany®) and a CLA 6.5/R10/2x128/BIP-1224 ultrasound probe (Vermon SA, Tours, France®). It was placed vertically in a tank filled with normal water with a temperature of 37.7°C. The temperature was measured by a TFX 422 precision thermometer (ebro Electronic GmbH & Co. KG, Ingolstadt, Germany®) and maintained stable during image acquisitions. This was done in order to simulate human rectum environment, see Fig. 3. Subsequently, the probe and the stepper were connected to LogicScan 128 INT-2Z beamformer (Telemed UAB, Vilnius, Lithuania®) and to BiopSee® 0.9.14 software (MedCom GmbH, Darmstadt, Germany and Pi Medical Ltd. Athens, Greece®).
Fig. 3
Experimental setup (front view): The TPV 061 template was connected to a GfM mechanical stepper and a CLA 6.5/R10/2x128/BIP-1224 ultrasound probe. It was placed vertically in a tank filled with normal water with a temperature of 37.7°C which was continuously observed by a TFX 422 precision thermometer
Experimental setup (front view): The TPV 061 template was connected to a GfM mechanical stepper and a CLA 6.5/R10/2x128/BIP-1224 ultrasound probe. It was placed vertically in a tank filled with normal water with a temperature of 37.7°C which was continuously observed by a TFX 422 precision thermometer
Acquisitions
The tested materials that covered the U/S probe were: 1) EA 4015 Silicone Standoff kit or 2) UA0059 Endocavity balloon filled with one of the following each time: a) water, 40 ml of Endosgel® (Farco-Farma GmbH, Germany®), b) 40 ml of Instillagel® (Farco-Farma GmbH, Germany®), c) 40 ml of Ultraschall gel (Dahlhausen & Co. GmbH, Germany®). Finally, d) Space OAR™ Gel (Augmenix, Waltham, MA, USA®) was injected into a Endocavity balloon and was arranged properly to form a 5.0 mm layer in front of U/S probe curved array (same thickness as Silicone Standoff kit).For each available material, acquisitions were carried out at two frequencies: 5 MHz and 7 MHz. This was done to make sure if the needle displacement is frequency-dependent. Prior to that, a pair of reference acquisitions at both frequencies was retrieved, with the ultrasound probe being “uncovered” in the water. After completing all acquisition cycles, the coordinates of the centre of each needle in the acquired image sets (centre of the needle fingerprint on U/S images) were reconstructed at the base, reference and apex plane. The position of the reconstructed needles from these image sets was compared with the position of reconstructed needles from the reference image set acquired with uncovered U/S probe in water for both frequencies.
Prostate volumetric analysis
Theoretical approach
Based on theoretical calculations of the sound speed in the different materials, it is expected that the use of the Silicone Standoff kit for coverage of the U/S probe during the image set acquisition will lead to volume and anatomical structure deformation and needle displacement. For that reason, the theoretically calculated shift was applied to cylindrical prostate models with various radii and distances from the probe surface. For each prostate model the percentage deviation of the distorted volume from the initial one was plotted in dependence on the distance from probe surface.
Experimental prostate volumetric calculations
In order to verify our theoretical volumetric calculations, the following experiment was designed. Using a tissue equivalent ultrasound prostate phantom-Model number 053 (CIRS Inc., Norfolk, VA, USA®) and previously mentioned equipment, two image acquisitions were performed (see Fig. 4). During the first acquisition, the ultrasound probe was covered by Silicone Standoff kit, and in the second one the Endocavity balloon was used. The amount of water inserted in the balloon was adjusted to form a layer above probe surface equal to the Silicone Standoff physical thickness of 5 mm, hence to include the same geometry. Afterwards, the prostate volume was contoured for each image set. The two PTV volumes were compared. The results acquired from the theoretical prediction were additionally evaluated.
Fig. 4
Tissue equivalent ultrasound prostate phantom-Model number 053 (CIRS Inc., Norfolk, VA, USA) arrangement: The TPV 061 template was connected to a GfM mechanical stepper and a ultrasound probe. 17 stainless steel trocar needles were implanted in the prostate phantom
Tissue equivalent ultrasound prostate phantom-Model number 053 (CIRS Inc., Norfolk, VA, USA) arrangement: The TPV 061 template was connected to a GfM mechanical stepper and a ultrasound probe. 17 stainless steel trocar needles were implanted in the prostate phantom
Dosimetric study
Using the setup described above, 17 stainless steel trocar needles were implanted in the prostate phantom, using the Endocavity balloon live images (Fig. 4). Afterwards, two image acquisitions were carried out in the previously described way: one was made with Silicone Standoff kit and the second with the use of the U/S probe covered by Endocavity balloon filled with water. The position of needles has remained unchanged all the time. The acquired datasets were transferred to Oncentra™ Prostate treatment planning system (Nucletron B.V., Veenendaal, Netherlands®).Next, the treatment plan (P1) was prepared based on the first image set acquired with the U/S probe covered by Silicone Standoff kit. The thickness of the used Silicone Standoff kit was 5 mm. The PTV and OARs were contoured and needles were reconstructed. The plan was calculated utilizing the inverse planning and optimization tool HIPO [25] available in Oncentra Prostate. All dosimetric indices and parameters of the plan satisfied our clinical protocol presented in Table 2. Afterwards, the second image set acquired using the U/S probe covered by Endocavity balloon filled with water was used. The amount of water inserted in the balloon was adjusted to form a layer above probe surface equal to the Silicone Standoff physical thickness of 5 mm. The PTV and OARs were contoured and the needles were reconstructed. The active dwell positions and dwell times calculated in the previously prepared plan (P1) were imported into the reconstructed catheters of present image set (Advance Plan Loading procedure in Oncentra Prostate®). Without any modification of the imported source dwell positions and times, the plan was created and the 3D dose distribution and the DVHs were calculated. This was considered to be a plan (P2) which corresponds to the factual situation. Finally, all dosimetric indices defined in our clinical protocol for both plans P1 and P2 were compared (see Table 2). In addition the Conformal Index [26] (COIN) and the External Index [27] (EI), for both plans were calculated and analysed.
Table 2
Clinical protocol (constraints) for HDR brachytherapy of prostate used at our clinic Dose-Volume parameterValue
Dose-Volume parameter
Value
D90 – Prostate
=100% of prescription dose: 11.5 Gy
V100 – Prostate
=90%
V150 – Prostate
=35%
D10 – Urethra
=115%: 13.2 Gy
D0.1 cm3 – Urethra
=120%: 13.8 Gy
D10 – Bladder and D10 – Rectum
=75%: 8.6 Gy
D0.1 cm3 – Bladder and D0.1 cm3 – Rectum
=80%: 9.2 Gy
D90: The dose that covers 90% of the PTV volume
V100: The percentage of prostate volume (PTV) that has received at least the prescription dose (PD)
V150: The percentage of prostate volume (PTV) that has received more than 150% of the PD
D10: the highest dose covering 10% of the urethra/rectum/bladder volume
D0.1 cm3: the dose for the most exposed 0.1 cm
3 of the urethra/ bladder/rectum
Clinical protocol (constraints) for HDR brachytherapy of prostate used at our clinic Dose-Volume parameterValueD90: The dose that covers 90% of the PTV volumeV100: The percentage of prostate volume (PTV) that has received at least the prescription dose (PD)V150: The percentage of prostate volume (PTV) that has received more than 150% of the PDD10: the highest dose covering 10% of the urethra/rectum/bladder volumeD0.1 cm3: the dose for the most exposed 0.1 cm
3 of the urethra/ bladder/rectumWe have compared COIN values including OARs, for the reference/prescribed dose to PTV. It considers also the conformity of the 3D dose distribution regarding the OARs [26, 28, 29]:The coefficient c
is the fraction of the PTV that is enclosed by the prescription dose. The coefficient c
is the fraction of the volume encompassed by the prescription dose that is covered by PTV. It is a measure of how much tissue outside the PTV is covered by the prescription dose. The c
is given as:where N
OAR is the total number of OARs, is the volume of the i-th OAR, is the dose limit defined for the i-th OAR and (D > ) is the volume of the i-th OAR that receives a dose that exceeds the dose limit .It is important to mention that while the coefficients c
and c
depend on the dose value under consideration, the coefficient c
depends only on the dose limits defined for each of OARs. The EI [27] is the ratio of the normal tissue volume outside of the PTV that receives a dose equal to or greater than the reference dose to the volume of the PTV. In this way the influence of using the Silicone Standoff kit on the 3D-dosimetry can be investigated.
Results
Theoretical calculations of the speed of sound
The speed of sound at 37.7° for the materials we used in our analyses were as follows: for distilled water 1524.9 m/s, normal water 1525.4 m/s, soft tissue 1540.0 m/s, and silicone rubber 981.0 m/s. The speed of sound in soft tissue was considered to be 1540 m/s, as referred in bibliography [11, 12].
Experimental evaluation of effect of materials on U/S image geometry
EA 4015 Silicone Standoff kit
Using a Silicone Standoff kit, we observed a radial needle shift of about 3 mm (Tables 3 and 4) at both frequencies (Figs. 5 and 6). Experimentally evaluated shift was in a good agreement with theoretical estimation, as displayed in Table 3. In addition, we found out that the needle displacement was independent from its distance from the probe and depended only on the wall thickness of the Silicone Standoff kit (Fig. 7). This was precisely what was expected theoretically (see the appendix).
Table 3
Comparison of theoretical calculations to experimental needle displacement for EA 4015 Silicone Standoff kit
Shift (mm)
Freq. (MHz)
Experimental
Theoreticala in vnormal water
Difference (%)a
Theoreticalb in vsoft tissue
Difference (%)b
5
3.03±0.13
2.77
−8.6
2.80
−7.6
7
3.04±0.14
2.77
−8.9
2.80
−7.9
Theoretical needle shift, in case of utilizing the real speed of sound in normal water in U/S software (Appendix)
Theoretical needle shift with actual software (Appendix)
Table 4
Absolute values of the average measured needle shift per standoff material used
Material
5 MHz
7 MHz
Shift (mm)
Shift (mm)
Silicone Standoff kit
3.03±0.11
3.05±0.13
Endocavity balloon filled with:
water
0.26±0.14
0.28±0.21
Endosgel®
0.29±0.16
0.20±0.12
Instillagel®
0.72±0.20
0.62±0.10
Ultrashall gel
0.32±0.25
0.29±0.17
Space OAR™
0.34±0.14
0.37±0.13
Fig. 5
Fused images (BiopSee® 0.9.14 software): reference image (acquisition made using an Endocavity balloon filled with water) and image from acquisition made utilizing a Silicone Standoff kit at reference plane at 5 MHz
Fig. 6
Radial shift of needles in U/S imaging with the use of Silicone Standoff kit at 5 MHz: a radial needle shift of 3 mm in direction from U/S probe is observed. Shift is presented according to the reference image acquisition made using an Endocavity balloon filled with water
Fig. 7
Needle displacement per template hole at 5 MHz and 7 MHz as the consequence of the application of EA 4015 Silicone Standoff kit on a U/S probe during the image acquisition. Influence of frequency on needle displacement is negligible
Fused images (BiopSee® 0.9.14 software): reference image (acquisition made using an Endocavity balloon filled with water) and image from acquisition made utilizing a Silicone Standoff kit at reference plane at 5 MHzRadial shift of needles in U/S imaging with the use of Silicone Standoff kit at 5 MHz: a radial needle shift of 3 mm in direction from U/S probe is observed. Shift is presented according to the reference image acquisition made using an Endocavity balloon filled with waterNeedle displacement per template hole at 5 MHz and 7 MHz as the consequence of the application of EA 4015 Silicone Standoff kit on a U/S probe during the image acquisition. Influence of frequency on needle displacement is negligibleComparison of theoretical calculations to experimental needle displacement for EA 4015 Silicone Standoff kitTheoretical needle shift, in case of utilizing the real speed of sound in normal water in U/S software (Appendix)Theoretical needle shift with actual software (Appendix)Absolute values of the average measured needle shift per standoff material used
UA0059 Endocavity balloon
The use of Endocavity balloon filled with water did not cause any significant shift (0.3 mm) in needle position on images (Table 4). The observed displacement had no specific direction and as such it could be attributed to the uncertainty in reconstruction of the actual position of the needle centre.
Endosgel® in UA0059 Endocavity balloon
Endosgel® along with Endocavity balloon did not affect needle imaging (0.2 mm) (Table 4). No particular direction of the shifts was observed.
Instillagel® in UA0059 Endocavity balloon
The insertion of Instilla gel in an Endocavity balloon had obvious effects in needle U/S imaging. The average disposition of the needles was 0.65 mm in the radial direction and towards the U/S probe (Fig. 8 and Table 4).
Fig. 8
Effects of insertion of Instillagel® in an Endocavity balloon at 5 MHz on the needles disposition: the average disposition of the needles is 0.65 mm in the radial direction and towards the U/S probe
Effects of insertion of Instillagel® in an Endocavity balloon at 5 MHz on the needles disposition: the average disposition of the needles is 0.65 mm in the radial direction and towards the U/S probe
Ultraschall gel in UA0059 Endocavity balloon
Ultraschall gel caused an average image distortion of 0.3 mm. As well as Endosgel, this shift seems to be statistical with no specific direction (Table 4).
Space OAR™ in UA0059 Endocavity balloon
Space OAR seems to compress anatomy, as needles were shifted towards the U/S probe for an average distance of 0.35 mm (Fig. 9 and Table 4).
Fig. 9
Effects of the use of Space OAR™ at 5 MHz on the needles disposition: needles are shifted towards the U/S probe for an average distance of 0.35 mm
Effects of the use of Space OAR™ at 5 MHz on the needles disposition: needles are shifted towards the U/S probe for an average distance of 0.35 mm
Prostate volume effect
Theoretical prostate volume effect with the use of Silicone Standoff kit
We have applied the radial transformation of 3 mm directed away from EA 4015 Silicone Standoff kit to different cylindrical prostate models (an example is shown in Fig. 10). The resulting transformed volume VS was significantly larger from the initial volume VW. Moreover, the size of this effect depended inversely on the prostate dimensions and its distance from the surface of the ultrasound probe. Our results are presented in Fig. 11.
Fig. 10
Prostate shape deformation per slice with the application of EA 4015 Silicone Standoff kit: this example shows radial deformation for an average shape prostate (∼20 mm) and an average distance from probe surface (∼11.8 mm)
Fig. 11
Percent volume deviation of spherical models of prostate of various radii and distances from the probe surface, with the application of EA 4015 Silicone Standoff kit
Prostate shape deformation per slice with the application of EA 4015 Silicone Standoff kit: this example shows radial deformation for an average shape prostate (∼20 mm) and an average distance from probe surface (∼11.8 mm)Percent volume deviation of spherical models of prostate of various radii and distances from the probe surface, with the application of EA 4015 Silicone Standoff kit
Experimental estimation of influence of material on prostate volume
The comparison between the two contoured prostates revealed a deviation in volume of the order of approximately 7.17%. The contouring from the image set from acquisition completed with the help of Endocavity balloon filled with water, resulted with the prostate volume (VW: 61979 mm3) of approximately 7.17% smaller than the one we have obtained from the image set acquired by covering the ultrasound probe by Silicone Standoff kit (VS: 66423 mm3): [(1-VS/VW)×100%]. This result supports our theoretical calculations presented in Fig. 11.Specifically, in our experimental case we used the tissue equivalent ultrasound prostate phantom of elliptical shape with the major dimension of 53 mm and the minor dimension of 37 mm. From here we calculated the equivalent radius of 20 mm [30]. The distance from U/S probe was around 10 mm. The theoretical volume difference for these parameters is 7.23%, see Fig. 11.Experimental analysis of dependence of dosimetric parameters on geometrical deformations caused by usage of different materials for US acquisition.As previously described, we found out that making the acquisition by covering the ultrasound probe by Silicone Standoff kit leads to prostate volume overestimation. We further observed that this geometry alteration of volume and needles possibly leads to dose over/underestimation.
Dose-Volume-Histogram
By calculating the plans P1 and P2 as defined previously in Material and methods, we found out that the large dose underestimation occurs. Evaluated DVH parameters are presented in Table 5. Plan P2 describes the factual situation, and the P1 shows the result in case of using the Silicon Standoff kit.
Table 5
Dose-Volume-Histogram parameters derived for the plans P1 (based on the image set acquisition made with the U/S probe covered by Silicone Standoff kit) and P2 (U/S probe covered with Endocavity balloon filled with water with dwell times calculated and imported from P1). Values in brackets are the relative difference of the parameters between the P1 and P2 in [%]
DVH-Parameter
(P1) EA 4015 Silicone Standoff kit* [%]
(P2) UA0059 Endocavity balloon filled with water with dwell times calculated and imported from (P1)* [%]
Prostate
D90
107.8
111.9 [3.8]
D95
101.4
103.5 [2.1]
D1cm3
287.4
285.5 [−0.7]
V100
95.7
96.5 [0.8]
V120
70.8
77.8 [9.8]
V150
24.7
30.5 [23.6]
V200
7.1
7.9 [10.8]
Urethra
D10
111.7
115.0 [2.9]
D0.1 cm3
114.7
115.5 [0.7]
D1 cm3
110.2
113.2 [2.7]
D1
115.2
118.1 [2.5]
V100
65.4
70.8 [8.3]
Rectum
D10
50.7
51.1 [0.8]
D0.1 cm3
68.7
66.4 [−3.3]
D2 cm3
55.4
52.5 [−5.2]
D1
65.3
64.4 [−1.3]
% of prescribed dose DP
Dose-Volume-Histogram parameters derived for the plans P1 (based on the image set acquisition made with the U/S probe covered by Silicone Standoff kit) and P2 (U/S probe covered with Endocavity balloon filled with water with dwell times calculated and imported from P1). Values in brackets are the relative difference of the parameters between the P1 and P2 in [%]% of prescribed dose DP
COIN
By observing the COIN, we noted the significant fall down of the COIN value in plan P2 (0.860) comparing to plan P1 (0.810), see Table 6. This means that the use of Silicone Standoff kit as the U/S probe coverage causes considerable overestimation of the COIN value.
Table 6
COIN and External Index values for the plans P1 and P2
COIN and COIN factors
(P1) EA 4015 Silicone Standoff kit
(P2) UA0059 Endocavity balloon filled with water with dwell times calculated and imported from (P1)
COIN
0.860
0.810
c1
0.957
0.965
c2
0.899
0.840
c3
1.000
0.999
EI
0.108
0.184
COIN and External Index values for the plans P1 and P2The coefficient c
which is the fraction of the PTV (prostate) that is enclosed by the prescription dose, plays insignificant role in this drop of COIN value (0.957 for P1 to 0.965 for P2). This means that the PTV coverage with prescription dose remains good in both cases. One could expect such a result, as the volume of the prostate used by calculation of plan P2 is lesser than that of P1 based on which the dwell times were calculated. In this case the COIN coefficient c3 has negligible influence on drop of COIN (1.000 for P1 to 0.999 for P2), as the dose of OARs has not crossed over the permitted limit defined in our clinical protocol. Main responsibility for the decrease of COIN has the reduction of a factor c2 (0.899 for P1 to 0.840 for P2). This is the consequence of the significant increase of volume of the healthy tissue outside of PTV covered with prescribed dose.
External index (EI)
Observing the external volume index has just confirmed the results and conclusions we obtained from COIN. The ratio of the normal tissue volume outside of the PTV that receives a dose equal to or greater than the reference dose to the volume of the PTV, is much higher in P2 that was calculated for the original plan P1: 0.184 for P2 compared to 0.108 for P1.
Discussion
The recorded shift of needles and other organ structures is radial, following the form of the sound wave transmission. There was no distance dependence of the shifts, as needles closer to the probe seem to have the same displacement with those in the edge positions. Additionally, no dependence on frequency was observed. Small deviations among frequencies could be attributed to uncertainty in recording the needle centre coordinates. The greatest shift was observed with the application of the Silicone Standoff kit (3 mm). This dislocation was expected from the theoretical calculations, due to the difference in sound speed between water and silicone rubber. Silicone Standoff kit makes the needles to appear in larger distances, causing an expansion of geometry. That indicates that using a high-resolution template type TPV 061, where distance between holes is 2.5 mm, a needle could have deviation of more than one hole from its original position. Moreover, the Silicone Standoff kit has the ability to widen the dimensions of every object included in the ultrasound image, as well as altering its shape. In our study we have shown that this image distortion results in an overestimation of a prostate volume which leads as a consequence to the dose underestimation. Instillagel®, along with the Endocavity balloon, has an adverse effect to the ultrasound image by “shrinking” objects. Geometry has radial shift directed to U/S probe (0.7 mm). This can be attributed to the fact that sound waves travel faster through Instillagel® than through water. Despite the fact that Ultraschall gel effect in needle shift was insignificant (0.3 mm) (Table 4), this gel is not recommended to be used as a balloon-filling material as its insertion in the balloon is quite difficult due to its high viscosity. Space OAR™ also compresses dimensions of appeared objects for an average value of 0.40 mm. This effect has specific direction towards the probe, revealing that speed of sound travels in a higher speed in Space OAR™ than water. The measured shift can be considered as negligible for clinical application when using layers of up to 5 mm thickness considered in our study.Speed of sound at 37.7°C in normal water (1525.4 m/s) is slightly closer to the speed of sound in soft tissue (1540 m/s) than that of pure water (1524.9 m/s). Therefore, both normal and distilled water are suitable for quality assurance tests, with the condition of a stable temperature.To calibrate accurately U/S system for prostate brachytherapy, biopsy and other interventional procedures, medium with the sound speed of 1540 m/s should be available. Goldstein et al.
[31] suggests the use of a 9% ethylene glycol-water mixture instead of tapwater in order to eliminate inaccuracy in QA tests in brachytherapy. Sound speed in the suggested mixture is equal to both soft tissue and U/S calibration sound speed (1540 m/s) at the room temperature. The errors described in the above-mentioned article (3.6%) coincide with the predicted error values (3.8%) of Figs. 12 and 13. According to literature available equations, the speed of sound of soft tissue can be achieved when using the tapwater by raising the water temperature to 47.4°C. An additional way is to increase the total dissolved solids or TDS concentration (add some salts) up to 15‰, while maintaining temperature at 37.7°C. To confirm that the above mentioned TDS value has been reached, electrical conductivity (EC) measurements should be performed at 25°C, until EC becomes equal to 23.5 mS/cm (see Equation 1). If nothing of the above is feasible, the potential calibration errors from Figs. 12 and 13 should be taken into account. This expected calibration error is given in dependence of the speed of sound in tapwater and the corresponding temperature. The values are derived from Chen-Millero equation [14] and further confirmed by Bilaniuk-Wong's [17], Marczak's [19] and Lubbers-Graaff's [22] equations for sound speed calculation.
Fig. 12
Predicted calibration error of U/S system in function with speed of sound in tap water (data derived from Chen-Millero equation [14])
Fig. 13
Predicted calibration error of U/S system in function with temperature in tap water (data derived from Chen-Millero equation [14])
Predicted calibration error of U/S system in function with speed of sound in tapwater (data derived from Chen-Millero equation [14])Predicted calibration error of U/S system in function with temperature in tapwater (data derived from Chen-Millero equation [14])
Conclusions
Our study has shown that the EA 4015 Silicone Standoff kit can not be utilized for image acquisition during HDR or LDR prostate brachytherapy, biopsy and other interventional procedures. It leads to erroneous placement and reconstruction of needles, contouring of PTV and OARs. In case of prostate brachytherapy this results in incorrect dose calculation. Effects of all other hydrogel materials can be consider as negligible in the performance of a safe brachytherapy procedure. Furthermore, the accurate performance of QA checks depends on the quality of the water used. If the acoustic velocity of the test water diverts from 1540 m/s, potential errors may affect accurate calibration of the U/S system. A possible solution to eliminate this deficiency is the choice of the right temperature and salt concentration of the QA water.
Authors: Karl Mikael Kälkner; Josef Nilsson; Marie Lundell; Seymour Levitt; Sten Nilsson Journal: Radiother Oncol Date: 2006-11-27 Impact factor: 6.280
Authors: A A Martinez; I Pataki; G Edmundson; E Sebastian; D Brabbins; G Gustafson Journal: Int J Radiat Oncol Biol Phys Date: 2001-01-01 Impact factor: 7.038
Authors: Christian Kirisits; Mark J Rivard; Dimos Baltas; Facundo Ballester; Marisol De Brabandere; Rob van der Laarse; Yury Niatsetski; Panagiotis Papagiannis; Taran Paulsen Hellebust; Jose Perez-Calatayud; Kari Tanderup; Jack L M Venselaar; Frank-André Siebert Journal: Radiother Oncol Date: 2013-11-30 Impact factor: 6.280