ORART Optimization example - cerr/CERR GitHub Wiki
ORART Optimization example
This example will work you through to show how to use CERR and IMRTP systems to generate 3D dose distribution based on the prescription of the case.
Here we will use an ORART test case: orart_test_case1_9sept05.mat (CERR planC structure) to demonstrate all steps and corresponding results.
Prescription for orart_test_case1_9sept05.mat
•72 Gy to the CTV 1 3mm structure.
•64 Gy to the CTV 2 3mm structure,
•60 Gy to the CTV 3 3mm structure.
•The mean dose to the parotid glands should be held as low as possible, but not at the expense of an adequate target dose
distribution. Preferably, one Parotid gland at least should be held below 26 Gy.
•The mean dose to the oral cavity should be held as low as possible, but not at the expense of an adequate target dose
distribution.
•The mandible should receive no more than 70 Gy max dose.
•The max to the cord should be 45 Gy (hard constraint),
•The max to the cord_3mm should be 50 Gy (hard constraint),
•The max to the brainstem (brainstem) should be 54 Gy (hard constraint).
•The max to the brainstem expansion (brainstem_3mm) should be 58 Gy (hard constraint).
•An adequate target dose distribution will have:
Min 93% of prescribed dose
Max <115%
1. Load orart_test_case1_9sept05.mat using CERR
300px|Three target volumes: CTV 1 3mm, CTV 2 3mm and CTV 3 3mm
2. Beam Calculation
1> Launch IMRTP module from CERR to create beamlets, click button “Equispaced” to create a 9 equispaced beams with starting angle at 0 along AP direction. Set the beam Energy as 6MV and isodistance as 100 cm
2> Select structures for computation, set up margin as 0.5cm , SampRate as 2 and mark the CTV 1 3mm, CTV 2 3mm and CTV 3 3mm as the target volumes.
3>Under panel “IM Parameter” Set up IM parameters: Choose “QIB” algorithm and select “threshold” as ScatterMethod.
4> under file, select “Recompute & add dosimetry”, click button “go” to compute. It will take 9 minutes to finish computation. The computed dosimetry will be stored in planC, under planC{12}(2).
3. Dose calculation
1> In matlab script: runOptimExample.m, define the prescription doses and relevant anatomical structures, and relative weights
2> run script: runOptimExample(1,1), the computed dose will add into planC structure.
4. Evaluation of resulted dose
1> the resulted dose is shown in the following figure. The max dose is 83.6 Gy. 300px|Resulted Dose Distribution 2> use DVH menu to display the DVHs of some strucures 300px|DVH Menu 300px|DVHs of selected structures
