SWIMRT: a graphical user interface using sliding window algorithm to construct a fluence map machine file

Quantified difference of the collapsed cone convolution (CCC) and Monte Carlo (MC) algorithms based on DVH and gamma analysis for cervical cancer radiation therapy

OBJECTIVE

To quantify the difference between the (collapsed cone convolution) CCC algorithm and the (Monte Carlo) MC algorithm and remind that the planners should pay attention to some possible uncertainties of the two algorithms when employing the two algorithms.

METHODS

Thirty patients' cervical cancer VMAT plans were designed with a Pinnacle TPS (Philips) and divided equally into two groups: the simple group (SG, target volume was only the PTV) and the complex group (CG, target volume included the PTV and PGTV). The plans from the Pinnacle TPS were transferred to the Monaco TPS (Elekta). The plans' parameters all remained unchanged, and the dose was recalculated. Gamma passing rates (GPRs) obtained from dose distribution from Pinnacle TPS compared with that from Monaco TPS with SNC software based on three triaxial planes (transverse, sagittal and coronal). GPRs and DVH were used to quantify the difference between the CCC algorithm in pinnacle TPS and the MC algorithm in Monaco TPS.

RESULTS

Among the statistical dose indexes in DVHs from the Pinnacle and Monaco TPSs, there were 7(7/15) dose indexes difference with statistically significant differences in the SG, and 10(10/18) dose indexes difference with statistically significant differences in the CG. With 3%/3 mm criterion, the most (5/6) GPRs were greater than 95% from the SG and CG. But with 2%/2 mm criterion, the most (5/6) GPRs were less than 90% from the two groups. In addition, we found that GPRs were also related to the selected triaxial planes and the complexity of the plan (GPRs varied with the SG and CG).

CONCLUSIONS

Obvious difference between the CCC and MC algorithms from Pinnacle and Monaco TPS. DVH maybe better than 2D gamma analysis on quantifying difference of the CCC and MC algorithms. Some attention should be paid to the uncertainty of the TPS algorithm, especially when the indicator on the DVH is at the critical point of the threshold value, because the algorithm used may overestimate or underestimate the DVH indicator.

Dosimetric Impact on the Flattening Filter and Addition of Gold Nanoparticles in Radiotherapy: A Monte Carlo Study on Depth Dose Using the 6 and 10 MV FFF Photon Beams

PURPOSE

This phantom study investigated through Monte Carlo simulation how the dose enhancement varied with depth, when gold nanoparticles (NPs) were added using the flattening filter-free (FFF) photon beams in gold NP-enhanced radiotherapy.

METHOD

A phantom with materials varying from pure water to a mixture of water and gold NPs at different concentrations (3-40 mg/mL) were irradiated by the 6 and 10 MV flattening filter (FF) and FFF photon beams. Monte Carlo simulations were carried out to determine the depth doses along the central beam axis of the phantom up to a depth of 40 cm. The dose enhancement ratio (DER) and FFF enhancement ratio (FFFER) were calculated based on the Monte Carlo results.

RESULTS

The DER values were found decreased with an increase of depth and increase of NP concentration in the phantom. For the maximum NP concentration of 40 mg/mL, the DER values decreased 6.9, 12, 4.6 and 7.2% at a phantom depth from 2 to 40 cm, using the 6 MV FF, 6 MV FFF, 10 MV FF and 10 MV FFF photon beams, respectively. The maximum DER values for the 6 MV beams were 1.08 (FF) and 1.14 (FFF), while those for the 10 MV beams were 1.04 (FF) and 1.07 (FFF). When the FF was removed from the linear accelerator head, the FFFER showed a more significant increase of dose enhancement for the 6 MV beams (1.057) than the 10 MV (1.031).

CONCLUSION

From the DER and FFFER values based on the Monte Carlo results, it is concluded that the dose enhancement with depth was dependent on the NP and beam variables, namely, NP concentration, presence of FF in the beam and beam energy. Dose enhancement was more significant when using the lower photon beam energy (i.e., 6 MV), FFF photon beam and higher NP concentration in the study.

Dosimetric variations in calculation grid size in prostate VMAT: a dose-volume histogram analysis using the Gaussian error function

Background

Varying the calculation grid size can change the results of dose-volume and radiobiological parameters in a treatment plan, and therefore has an impact on the treatment planning quality assurance.

Purpose

This study investigated the dosimetric influence of the calculation grid size variation in the prostate volumetric modulated arc therapy (VMAT) plan.

Methods and materials

Dose distributions of 10 prostate VMAT plans were acquired using calculation grid sizes of 1–5 mm. Dose-volume histogram (DVH) analysis was carried out to determine the dose-volume variation corresponding to the grid size change using the Gaussian error function (GEF). At the same time, dose-volume points, dose-volume parameters and radiobiological parameters were calculated based on DVHs of targets and organs at risk (OARs) for each grid size.

Results

Comparing percentage variations of GEF parameters between the planning target volume (PTV) and clinical target volume (CTV), GEF parameters of the PTV were found varied more significantly than the CTV. This resulted in larger variations of dose-volume (%ΔCI=40·02 versus 13·55%, %ΔHI=12·45 versus 2·93% and %ΔGI=0·22 versus 0·06%) and radiobiological parameters (%ΔTCP=0·61 versus 0·25% and %ΔEUD=2·11 versus 0·26%) of the PTV compared with CTV. For OARs, the rectal wall showed a larger dose-volume variation than the rectum. However, similar dose-volume variation due to grid size change was not found in the bladder, bladder wall and femur.

Conclusions

Knowing the dosimetric variation in this study is important to the radiotherapy staff in the quality assurance for the prostate VMAT planning.

Empirical Driven Automatic Detection of Lobulation Imaging Signs in Lung CT

Computer-aided detection (CAD) of lobulation can help radiologists to diagnose/detect lung diseases easily and accurately. Compared to CAD of nodule and other lung lesions, CAD of lobulation remained an unexplored problem due to very complex and varying nature of lobulation. Thus, many state-of-the-art methods could not detect successfully. Hence, we revisited classical methods with the capability of extracting undulated characteristics and designed a sliding window based framework for lobulation detection in this paper. Under the designed framework, we investigated three categories of lobulation classification algorithms: template matching, feature based classifier, and bending energy. The resultant detection algorithms were evaluated through experiments on LISS database. The experimental results show that the algorithm based on combination of global context feature and BOF encoding has best overall performance, resulting in F1 score of 0.1009. Furthermore, bending energy method is shown to be appropriate for reducing false positives. We performed bending energy method following the LIOP-LBP mixture feature, the average positive detection per image was reduced from 30 to 22, and F1 score increased to 0.0643 from 0.0599. To the best of our knowledge this is the first kind of work for direct lobulation detection and first application of bending energy to any kind of lobulation work.

Some computer graphical user interfaces in radiation therapy

In this review, five graphical user interfaces (GUIs) used in radiation therapy practices and researches are introduced. They are: (1) the treatment time calculator, superficial X-ray treatment time calculator (SUPCALC) used in the superficial X-ray radiation therapy; (2) the monitor unit calculator, electron monitor unit calculator (EMUC) used in the electron radiation therapy; (3) the multileaf collimator machine file creator, sliding window intensity modulated radiotherapy (SWIMRT) used in generating fluence map for research and quality assurance in intensity modulated radiation therapy; (4) the treatment planning system, DOSCTP used in the calculation of 3D dose distribution using Monte Carlo simulation; and (5) the monitor unit calculator, photon beam monitor unit calculator (PMUC) used in photon beam radiation therapy. One common issue of these GUIs is that all user-friendly interfaces are linked to complex formulas and algorithms based on various theories, which do not have to be understood and noted by the user. In that case, user only needs to input the required information with help from graphical elements in order to produce desired results. SUPCALC is a superficial radiation treatment time calculator using the GUI technique to provide a convenient way for radiation therapist to calculate the treatment time, and keep a record for the skin cancer patient. EMUC is an electron monitor unit calculator for electron radiation therapy. Instead of doing hand calculation according to pre-determined dosimetric tables, clinical user needs only to input the required drawing of electron field in computer graphical file format, prescription dose, and beam parameters to EMUC to calculate the required monitor unit for the electron beam treatment. EMUC is based on a semi-experimental theory of sector-integration algorithm. SWIMRT is a multileaf collimator machine file creator to generate a fluence map produced by a medical linear accelerator. This machine file controls the multileaf collimator to deliver intensity modulated beams for a specific fluence map used in quality assurance or research. DOSCTP is a treatment planning system using the computed tomography images. Radiation beams (photon or electron) with different energies and field sizes produced by a linear accelerator can be placed in different positions to irradiate the tumour in the patient. DOSCTP is linked to a Monte Carlo simulation engine using the EGSnrc-based code, so that 3D dose distribution can be determined accurately for radiation therapy. Moreover, DOSCTP can be used for treatment planning of patient or small animal. PMUC is a GUI for calculation of the monitor unit based on the prescription dose of patient in photon beam radiation therapy. The calculation is based on dose corrections in changes of photon beam energy, treatment depth, field size, jaw position, beam axis, treatment distance and beam modifiers. All GUIs mentioned in this review were written either by the Microsoft Visual Basic.net or a MATLAB GUI development tool called GUIDE. In addition, all GUIs were verified and tested using measurements to ensure their accuracies were up to clinical acceptable levels for implementations.

Intensity modulated radiation therapy with irregular multileaf collimated field: a dosimetric study on the penumbra region with different leaf stepping patterns

Using a Varian 21 EX linear accelerator with a multileaf collimator (MLC) of 120 leaves, the penumbra regions of beam profiles within an irregular multileaf collimated fields were studied. MLC fields with different leaf stepping angles from 21.8 degrees to 68.2 degrees were used. Beam profiles in different directions: (1) along the cross-line and in-line axis, (2) along the leaf stepping edges of the field, and (3) parallel to the stepping edges but in the middle of the field, were measured and calculated using Kodak XV radiographic film and Pinnacle3 treatment planning system version 7.4f. These beam profiles were measured and calculated at source to axis distance= 100 cm with 5 cm of solid water slab on top. On the one hand, for both cross-line and in-line beam profiles, the penumbra widths of 20%-80% did not vary with the leaf stepping angles and were about 0.4 cm. On the other hand, the penumbra widths of 10%-90% of the above two profiles varied with the stepping angles and had maximum widths of about 1.9 cm (cross-line) and 1.65 cm (in-line) for stepping angles of 38.7 degrees and 51.3 degrees , respectively. For profiles crossing the "rippled" stepping edges of the field, the penumbra widths (10%-90%) at the regions between two opposite leaves (i.e., profile end at the Y1/Y2 jaw position) decreased with the stepping angles. At the penumbra regions between two leaf edges with the tongue-and-groove structure of the same bank (i.e., profile end at the X1/X2 jaw position), the penumbra widths increased with the stepping angles. When the penumbra widths were measured between two opposite leaf edges and at corners between two leaves, the widths first decreased with the stepping angles and then increased beyond the minimum width point at stepping angle of 45 degrees. The penumbra width (10%-90%) measured at the leaf edge was larger than that at the corner. For the beam profiles calculated using Pinnacle3, it is found that the results agreed well with the measurements along the cross-line and in-line axis, while there was a deviation for the profiles along the leaf stepping edge of the field compared to the film measurements. The measured results in this study can help us to understand the dosimetric effect of the leaf stepping (due to finite leaf width), tongue-and-groove and rounded leaf end structure in the penumbra region of an irregular MLC field. A more dedicated penumbra model can be developed for the treatment planning system.

Measurement for the MLC leaf velocity profile by considering the leaf leakage using a radiographic film

A method to measure the velocity profile of a multi-leaf collimator (MLC) leaf along its travel range using a radiographic film is reported by considering the intra-leaf leakage. A specific dynamic MLC field with leaves travelling from the field edge to the isocentre line was designed. The field was used to expose a radiographic film, which was then scanned, and the dose profile along the horizontal leaf axis was measured. The velocity at a sampling point on the film can be calculated by considering the horizontal distance between the sampling point and the isocentre line, dose at the sampling point, dose rate of the linear accelerator, the total leaf travel time from the field edge to isocentre line and the pre-measured dose rate of leaf leakage. With the leaf velocities and velocity profiles for all MLC leaves measured routinely, a comprehensive and simple QA for the MLC can be set up to test the consistency of the leaf velocity performance which is essential to the IMRT delivery using a sliding window technique.