Effects of different CT value assignment methods on dose calculations in radiotherapy planning for brain metastases.

BACKGROUND: To study the effects of different CT value assignment methods on the dose calculations in radiotherapy plans for brain metastases, this study will provide a reference for radiotherapy planning design based on MR images.
METHODS: All fifty recruited patients underwent CT and MR simulated localization the same day as, but prior to, three dimensional conformal radiation therapy (3D-CRT) or intensity modulated radiation therapy (IMRT) for brain metastases. After rigid registration of both the CT and MR images, the main tissues and organs were delineated on the CT and MR images. The average CT value of each tissue or organ was calculated. Three groups of pseudo-CT were generated by three CT value assignment methods: (I) the whole tissue was assigned 140 HU; (II) cavity, bone and other tissues were assigned -700, 700 and 20 HU, respectively; (III) tissue- and organ-specific CT values were given. The dose distribution was recalculated based on the three groups of pseudo-CT to obtain Plan2, Plan3 and Plan4, accordingly. The resultant radiotherapy plans were considered the original plan (Plan1). Then, the dosimetric differences between these three plans and Plan1 were compared.
RESULTS: The average pseudo-CT values of bone and cavity were 731.7±69.3 and -725.5±26.1 HU, respectively. The range of average soft-tissue CT values was from -70 to 70 HU. The dose distribution between Plan1 and Plan2, Plan3 or Plan4 showed some differences, and the differences decreased in turn. The differences in the maximum dose of the lenses can reach 5.0%, 1.5% and 1.2%, respectively, while the differences in other dose parameters (maximum dose, mean dose and D 98% to the PTV, D 5% of the brainstem, and maximum dose of the brainstem, corpus callosum, left eye, right eye) were basically less than 2.0%, 1.2% and 0.8%, respectively. This shows that in the CT value assignment method, the dose calculation error can be greatly reduced by assigning the value to the bone and cavity separately, and if the different soft tissues are distinguished, the error of the dose calculation can be further reduced by more than 30%. In the pixel-by-pixel dosimetric comparison, the areas of more than 1% dose difference between Plan1 and Plan3 as well as Plan4 were mainly distributed near skin while those between Plan1 and Plan2 were mainly distributed at the bone, cavity, bone and soft tissues junction, and the skin near the field.
CONCLUSIONS: In summary, a scheme for assigning specific CT values to MRI-based radiotherapy is established. The scheme will provide patients with a dose-free radiotherapy plan. Through the calculation of the differences between the new plans and the old plan, it is found that our scheme can basically control the dose error below 0.8% to meet the clinical requirements. 2020 Translational Cancer Research. All rights reserved.

Introduction

Brain metastasis is a secondary malignant tumor, mainly from lung cancer and breast cancer, with a usually rapid onset, short course of disease, and poor therapeutic response (1-3). Three-dimensional conformal radiotherapy (3D-CRT) and intensity-modulated radiotherapy (IMRT) are commonly used in the treatment of multiple or whole-brain brain metastases. As the general standard of radiotherapy simulation localization and treatment planning design, CT provides spatial location and electron density information for calculating radiotherapy doses (4). CT imaging is sensitive to hard tissue, while for target delineation, it is difficult to distinguish the tumor boundary in CT images, which is the key factor restricting the efficacy of radiotherapy (5,6).

In the radiotherapy of head and neck tumors, MR images provide good soft tissue contrast and multiplanar imaging (6,7). The fusion of MR and CT images greatly improves the accuracy of target mapping (8,9). This may lead to other registration errors at the same time, however. If only the MR images are used to design the radiotherapy plan, this error can be avoided and redundant CT simulation positioning can be reduced to simplify the radiotherapy process. However, the lack of electronic density information is a major obstacle to treatment planning based on MR images alone (10,11). To overcome this difficulty, there have been some reports in recent years wherein the whole area is given uniform density, or at least the bones, cavities and soft tissues are segmented and given uniform density in MR images for treatment planning (12,13). However, there are some differences in CT values (or electron density values) between different soft tissues. The dose calculation error caused by giving uniform CT values to all soft tissues, including target volumes and organs at risk (OARs), which play a key role in the radiotherapy plan evaluation, may have a serious adverse effect and a certain degree of blindness. For this reason, we compared the effects of an overall uniform CT value assignment method, a bone, air and soft tissue CT value assignment method, and a method for assigning different CT values to different tissues on the calculation of the planned dose in 3D-CRT or IMRT for brain metastases.

Methods

Patient data

Fifty patients with brain metastases who underwent radiotherapy at Shandong Cancer Hospital from July 2018 to July 2019 were selected. Each patient underwent CT simulation localization and MR simulation localization before radiotherapy. There were 29 males and 21 females, aged 44 to 78 years, with a median age of 60 years.

Postural fixation and simulated localization

Prior to radiotherapy, all patients were treated in the supine position with a fixed head and neck thermoplastic mask fixed for CT simulated positioning on large bore spiral CT device (Philips Healthcare Company, Netherlands). On the same day, MR simulation positioning was performed using a 3T MR system (Discovery 750w, GE Healthcare, the USA) equipped with 6ch Neuro Flex coil. T2 propeller MR sequence was applied for brain imaging with the corresponding scan parameters of TR =13,312 ms, TE =113.5 ms, matrix = 416×416, FOV = 260×260 mm. The scan time was 2 minutes 40 seconds. The scanning range was above the neck and the scanning thickness was 3 mm.

Delineation of tumor volume and OARs

The images from the CT and MR simulation positioning processes were introduced into the Eclipse 13.5 treatment planning system (Varian Healthcare Company, the USA), where they were registered rigidly. Then, the whole head (the whole scanning range of CT images), the target volume—gross tumor volume (GTV), planning target volume (PTV), and the OARs—the brainstem, spinal cord, eyeball, lens, optic nerve, and parotid glands—were delineated on the CT and MR images. Each patient was delineated by the same clinician.

Plan designs

The physicists developed the 5–7 field coplanar IMRT plans or 3D-CRT plans in the Eclipse 13.5 program system, with both an IMRT plan and a 3D-CRT plan developed for 10 patients, for a total of 60 plans. According to the type of plan, there were 37 IMRT plans and 23 3D-CRT plans, and according to the type of target volume, there were 26 local target volume plans and 34 whole-brain target volume plans. The prescription dose of the PTV was 40 to 60 Gy and the fractionation dose was 1.8 to 2.5 Gy/every fraction, 5 fractions a week.

Tissue delineation and CT value analysis

Bone, air cavities, the corpus callosum, the tissue of the cerebrum and brainstem, and the skin were further delineated on MR images. The remaining tissue was defined as other soft tissues. The average CT values of each organ and tissue were calculated and recorded in Eclipse 13.5.

CT values assignment to tissues and dose recalculation

Based on the CT values obtained in the previous step, the tissues of the original CT were reassigned CT values. CT1 consisted of the original CT; in CT2, the whole tissue was set to 140 HU; in CT3, cavities, bone and other tissues were set to −700, 700 and 20 HU, respectively; in CT4, cavities were set to 700 HU, bone to 700 HU, the brainstem to 30 HU, the corpus callosum to 12 HU, the brain (excluding the brainstem and corpus callosum) to 35 HU, the parotid gland to −15 HU, skin to −50 HU, eyeball to 25 HU, the optic nerve to 10 HU, and lens to 60 HU. Then, the original plan (Plan1) was copied to CT2, CT3 and CT4 for recalculating the doses, and the obtained plans were recorded as Plan2, Plan3 and Plan4, respectively.

Dosimetric comparison between the new plans and Plan1

The dose calculation errors of these new plans were assessed through comparison of the conformity indexes (CIs) and heterogeneity indexes (HIs) of the PTVs comparison, point dose comparison, dose-volume histogram (DVH) analysis and pixel dose comparison

The value of CI is calculated as follows (14):

In Eq. [1], Vt represents the target volume, Vt,ref represents the target volume of the reference isodose, and Vref represents the total volume of all the areas of the reference isodose. HI is calculated as follows (14):

In Eq. [2], D2%, D98%, and Dmean represent the 2% maximum volume dose, the 98% maximum volume dose and the mean dose of the target volume, respectively. The closer to 1 the CI is, the better the fitness is, and the closer to 0 the HI is, the better the uniformity is.

The parameter for the point dose comparison is the dose of the isocentric point (DISO).

The parameters for DVH analysis are the maximum dose (Dmax) of the PTV, the mean dose (Dmean) of the PTV, D98% of the PTV, Dmean of the head, D5% of the brainstem, and Dmax of the brainstem, corpus callosum, left eye, right eye, left lens, right lens and spinal cord.

Pixel dose comparison

The CT and MR images and the four sets of plans were imported into a Raystation 7.0 workstation (RaySearch Healthcare Company, Sweden). Dose difference silhouettes of Plan2, Plan3, Plan4 and Plan1 were made, and the pixel doses were compared.

Radiotherapy plans were made on MR images

Proper CT values were assigned to different tissues and organs on MR images to make IMRT plans in 10 cases on the Raystation 7.0 workstation. Then the average time of the plan was calculated.

Statistical analysis

The data were expressed as the means ± SDs, using SPSS 20.0 statistical software. Matched-pairs t tests were used to compare the data, and P<0.05 indicated a statistically significant difference.

Results

Average CT values

The average CT values of the head, bone and air cavity were 141.0±25.7, 731.7±69.3 and −725.5±26.1 HU, respectively. The average CT values of the soft tissues were mostly distributed between −20 and 40 HU (the values for the brainstem, corpus callosum, brain tissue, parotid glands, eyeball and optic nerve were 29.1±1.4, 12.0±1.9, 35.4±2.9, −17.2±24.3, 25.8±2.7 and 10.1±9.9 HU, respectively), with the exception of the lens and skin, which were 60.1±6.8 and −51.3±16.8 HU, respectively. The difference in CT values among the different tissues was statistically significant (P<0.05).

CI and HI of the PTV

The HI of Plan2 was 6.4% higher than that of Plan1 in the local target area, and the difference was statistically significant (P<0.05). There was no significant difference in other CIs, HIs between Plan1 and the other plans. See for details.

Table 1

CI and HI of the PTV in the initial plan and in the three CT value assignment plans (mean ± SD)

Projects	Local target cases			Whole-brain target cases
	CI	HI		CI	HI
Plan1	0.792±0.088	0.088±0.025		0.835±0.097	0.180±0.176
Plan2	0.788±0.094	0.093±0.030		0.828±0.109	0.187±0.181
Plan3	0.793±0.090	0.088±0.029		0.833±0.109	0.185±0.177
Plan4	0.792±0.089	0.088±0.028		0.833±0.106	0.183±0.174
P
Plan2 vs. Plan1	0.197	0.008		0.373	0.082
Plan3 vs. Plan1	0.466	0.671		0.329	0.335
Plan4 vs. Plan1	0.487	0.277		0.317	0.476

CI, conformity index; HI, heterogeneity index; PTV, planning target volume.

General dosimetric difference

The dosimetric difference between Plan2, Plan3, Plan4 and Plan1 decreased in turn. For the 3D-CRT plan, the average differences for each dose parameters were −4.8% to 0.6%, −1.2% to 0.1%, and −0.9% to 0.1%, respectively, which were higher than the −3.2% to 0.4%, −0.9% to 0.1%, and −0.8% to 0.1% in the IMRT plan, respectively. The difference in the whole-brain target area cases was larger than that in the local target area cases, see for details.

Table 2

Relative differences in PTV dose parameters between the three CT value assignment plans and the initial plan (%, mean ± SD)

Dose parameters of PTV	Local target cases				Whole-brain target cases
	Plan2 vs. Plan1	Plan3 vs. Plan1	Plan4 vs. Plan1		Plan2 vs. Plan1	Plan3 vs. Plan1	Plan4 vs. Plan1
Dmax	0.37±0.63a	−0.08±0.24	−0.02±0.26		0.52±1.07a	0.06±0.45	0.14±0.46
Dmean	0.21±0.66a	0.14±0.24a	0.06±0.21		−0.32±0.55a	−0.01±0.39	−0.04±0.16
D98%	−0.19±0.68	0.08±0.30	0.04±0.21		−0.50±0.77a	−0.24±0.40a	−0.20±0.34a

a, indicates that the difference is statistically significant (P<0.05). PTV, planning target volume.

Table 3

Relative differences in the regions of interest dose parameters between the three CT value assignment plans and the initial plan (%, mean ± SD)

Regions of interest	Dose parameters	Local target cases				Whole-brain target cases
		Plan2 vs. Plan1	Plan3 vs. Plan1	Plan4 vs. Plan1		Plan2 vs. Plan1	Plan3 vs. Plan1	Plan4 vs. Plan1
Isocentric	DISO	−0.52±0.88	0.22±0.37	−0.06±0.18		−0.92±0.71a	0.12±0.32	−0.08±0.20
Head	Dmean	−0.46±0.61a	−0.31±0.41a	−0.15±0.16a		−0.51±0.53a	−0.38±0.18a	−0.21±0.16a
Brainstem	D5%	−0.81±1.13a	−0.06±0.21	−0.12±0.14		−0.66±0.71a	−0.04±0.18	−0.06±0.05
	Dmax	−0.72±1.06a	−0.21±0.28a	−0.14±0.22a		−0.62±1.13a	−0.22±0.17a	−0.15±0.21a
Corpus callosum	Dmax	−0.72±0.82a	0.14±0.18a	−0.03±0.14		−0.62±0.53a	0.06±0.16	−0.03±0.13
Left eyeball	Dmax	−1.12±2.40a	0.11±1.09	−0.16±0.58		0.47±0.83a	−0.26±0.25a	−0.12±0.26
Right eyeball	Dmax	−0.88±1.82a	−0.32±1.11	−0.35±0.42		0.36±0.79a	−0.21±0.32	−0.16±0.21a
Left lens	Dmax	−5.42±7.21a	−0.25±1.24	−0.36±0.92		−3.46±1.48a	−0.91±0.53a	−0.41±0.52a
Right lens	Dmax	−3.26±4.93a	−0.26±0.95	−0.40±0.59a		−3.62±1.82a	−0.80±0.67a	−0.46±0.50a
Spinal cord	Dmax	−0.19±1.17	−0.88±1.42a	−0.64±1.09a		−0.25±1.36	−1.06±0.54a	−0.80±0.51a

a, indicates that the difference is statistically significant (P<0.05).

Dosimetric difference in the PTV

For the local target area, Plan2 and Plan3 were slightly different from Plan1, while there was no significant difference between Plan4 and Plan1. For the whole-brain target area, Plan2 was significantly different from Plan1, in which the Dmax of the PTV increased by (0.52±1.07)%, and Dmean and D98% of the PTV decreased by (0.32±0.55)% and (0.50±0.77)%, respectively. Plan3 and Plan4 were slightly different from Plan1. See for details.

The dosimetric difference in the OARs

Most of the dose parameters of the Plan2 were significantly different from that those of the Plan1, and the difference of Dmax of the lenses of some patients was more than 5.0%. The differences of dose parameters between Plan3, Plan4 and Plan1 were greatly reduced. Except that the Dmax of the left and right lenses and spinal cord can reach 1.5% and 1.2%, respectively. The upper limit of the 95% confidence interval of the dosimetric difference in the other tissues and organs was less than 1.2% and 0.8%, respectively. In particular, the differences in the isocentric point dose, D5% and Dmax of the brainstem and Dmax of the corpus callosum was smaller. See , for details.

Figure 1

Local target case. (A) Is the dose distribution map of Plan1; (B) is the dose difference distribution map between Plan2 and Plan1; (C) is the dose difference distribution map between Plan3 and Plan1; and (D) is the dose difference distribution map between Plan4 and Plan1. As shown in , the areas where the difference was greater than 1% between Plan2, Plan3, Plan4 and Plan1, which decreased in turn, were mainly distributed at the skin near the field, and the pixel doses of Plan4 and Plan1 were basically the same.

Figure 2

Whole-brain target case. (A) Is the dose distribution map of Plan1; (B) is the dose difference distribution map between Plan2 and Plan1; (C) is the dose difference distribution map between Plan3 and Plan1; and (D) is the dose difference distribution map between Plan4 and Plan1. As shown in , the areas where the difference was greater than 1% between Plan2, Plan3, Plan4 and Plan1, which decreased in turn, were mainly distributed at the boundary between bone, air cavity and soft tissue, and the skin near the field. The difference between Plan4 and Plan1 was the smallest but still higher than that in local target cases.

Calculation of the IMRT plan on MR images

The group CT value assignment method was used to develop the IMRT plan on the MR images. Compared with the CT images, the time required to create the plan on the MR images was increased by 10 to 30 minutes. The main reason for the increased workload was the delineation of bone, air, skin and the corpus callosum and the assignment of CT values to tissues and organs. The dose distribution and DVH display are shown in .

Figure 3

Dose distribution of the IMRT plan based on MRI. IMRT, intensity modulated radiation therapy; MRI, magnetic resonance imaging.

Discussion

In recent years, there have been some reports on the use of the CT value assignment method for MRI radiotherapy plan design (12,13), but these previous reports often did not subdivide the soft tissue for assigning CT values, which may lead to blindness in soft tissue assignment. Previous studies have not deeply explored the effect of CT valuation on the uniformity and conformability of target dose nor the effect of some non-OARs (such as skin, bone, etc.) on radiation dose. In addition, there has been no study of CT valuation in relation to brain metastases. The most important feature of brain metastases is that not only local tumor target area irradiation but also whole-brain target area irradiation is very common. Moreover, the distribution of bone, air and soft tissue in the head is relatively less complex, and individual differences are small, which can yield more instructive results. In this study, we further studied the above three aspects.

Studies have shown that the CT value of soft tissue is mainly in the range of −70 to 70 HU, but the CT value of widely distributed bones is often higher than that of 600 HU (15,16). In this study, it was also found that most of the soft tissue CT values of the head ranged from −20 to 40 HU, with the exception of the lens and skin, in which the lens and the skin had the highest and lowest soft tissue CT values, approximately 60 and −50 HU, respectively.

In radiotherapy, the accuracy of dose calculation is crucial, and the wrong dose calculation can mislead clinicians and physicists from choosing the right treatment (17,18). High-density bones absorb more radiation, so the correct assignment of bone electron density will greatly affect the accuracy of dose calculation (19). At the same time, both the tumor target and the OARs are soft tissue, and the fine segmentation of and CT value assignment for each soft tissue may further improve the accuracy of dose calculation. Karotki et al. (12) studied 10 patients with head and neck tumors and found that when all tissues were given a uniform density of 1 g/cm3, the error in dose calculation could reach 4% to 5%. When bone, cavities and other tissues were given uniform densities of 1.5, 0 and 1 g/cm3 respectively, the error in dose calculation was less than 2%. Young et al. (13) studied 10 patients with nasopharyngeal carcinoma and found that when all tissues were given a uniform density of 1 g/cm3, the error in dose calculation could be close to 4%. When bone, cavities and other tissues were assigned separately, the dose calculation error could often be reduced to less than 2%, and only the bone and soft tissue boundary in some patients could reach 2% to 3%. Similarly, in this study, it was found that when all tissues were given uniform CT values (CT2), and when bone, cavities and other tissues were given CT values accordingly (CT3), the differences in the Dmax of the lens can reach 5.0%, and 1.5%, respectively, while the differences in the other dose parameters were basically less than 2.0%, and 1.2%, respectively. We further find that when tissue- and organ-specific CT values were given (CT4), the dose calculation error could be further reduced to less than 0.8%. This shows that in the CT value assignment method, the dose calculation error can be greatly reduced by assigning values to the bone and cavity separately, and if the different soft tissues can be distinguished, the error in the dose calculation can be further reduced by more than 30%. This can meet the clinical dose requirements.

The CI and HI of the target volume are very important for evaluating the treatment plan and cannot be intuitively obtained from other dose parameters. In this study, when different CT value assignment methods were used for dose recalculation, there was no significant difference in the CI and HI of the target volume compared with Plan1 (P>0.05), except for the HI in Plan2 for local target cases. One possible reason for this is that in brain metastases, the CT values of the target volume should be close to 35 HU for the brain tissue, while the average CT value of the head given by the uniform assignment method was 140 HU.

In the pixel and point dose comparison, differences between Plan2, Plan3 and Plan4 and Plan1 of more than 1% decreased in turn in certain areas; for the local target cases, these areas were mainly distributed near the skin near the field, while in the whole-brain target cases, they were mainly distributed at the junction of bone, cavity and soft tissue, as well as near the skin of the field. It was also found that the dosimetric differences in the 3D-CRT plans were larger than those in the IMRT plans, which was similar to the results of the study by Kartutik et al. (20). The possible reason for this is that the IMRT plans have more freedom in performing dose calculations. The differences in the dose calculation of multiple free fields can compensate for each other, which reduces the difference in total dose calculation. At the same time, the dosimetric differences in the whole-brain target cases were larger than those in the local target cases. The reason is that the target area of the former is larger, so the high dose area is closer to bone and air, and the CT value assignment methods have a greater influence on the dose calculation.

In addition, in this study, the tissue- and organ-specific CT value assignment method was used to make 10 IMRT plans based on MR images. A single plan developed on MR images took 10 to 30 minutes more than on CT images. This is mainly because the MR plans require more tissues and organs to be outlined and the CT values to be assigned. However, in view of the increasing improvements in and convenience of automatic profiling software, this difficulty is expected to be overcome in the future.

For the MR plans for brain metastases, we suggest that the tissue- and organ-specific CT value assignment method be used to generate a pseudo-CT for planning design. The CT values of each tissue and organ can be assigned according to the results of this study. Given that there may be physiological differences between different populations and that young patients (less than 45 years old) may also have physiological differences with the case samples in this study, other departments can also study a set of local CT value distributions according to the local actual situation. With the implementation of MR simulation positioning and MR image-guided radiotherapy technology in increasing numbers of departments, radiotherapy planning based on MR images will have increasingly broad prospects.