Monte Carlo based treatment planning for neutron capture therapy

• US 5,341,292 A
• Filed: 06/04/1992
• Issued: 08/23/1994
• Est. Priority Date: 06/04/1992
• Status: Expired due to Fees

• Current Assignee: New England Medical Center Hospitals, Inc.
• Original Assignee: New England Medical Center Hospitals, Inc.
• Inventors: Zamenhof, Robert G.
• Primary Examiner(s): Weinhardt, Robert A.
• Application Number:
  US07/893,650
• Time in Patent Office: 810 Days
• Field of Search: 364/413.26, 364/413.22, 364/413.19, 364/413.14, 364/413.13, 128/653.1, 600/4
• US Class Current: 600/425
• CPC Class Codes: A61N 2005/1034 Monte Carlo type methods; particle tracking
• CPC Class Codes: A61N 2005/109 Neutrons
• CPC Class Codes: A61N 5/1031 using a specific method of dose optimization
• CPC Class Codes: G01T 1/2985 In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis); using external radiation sources A61B6/02
• CPC Class Codes: G16H 20/40 relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
• CPC Class Codes: G16H 30/20 for handling medical images, e.g. DICOM, HL7 or PACS
• CPC Class Codes: G16H 30/40 for processing medical images, e.g. editing
• CPC Class Codes: G16H 50/50 for simulation or modelling of medical disorders

Abstract

A method for developing a treatment plan for neutron capture therapy of a target organ including the steps of using a diagnostic cross-sectional scanner to scan the target organ to generate cross-sectional images of the organ; generating a three-dimensional model of the organ from the cross-sectional images, the model including information about internal features of the organ; and using the three-dimensional model to compute expected radiation doses throughout the organ assuming a given arrangement of one or more neutron beams.

Claims

1. A method of developing a treatment plan for neutron capture therapy of a target organ containing a tumor, said method comprising:

using a diagnostic cross-sectional scanner to scan the target organ to generate cross-sectional images of the organ;

generating a three-dimensional model of the organ from the cross-sectional images, said model including information about variations in internal composition of the organ in addition to said tumor; and

computing a three-dimensional array of expected radiation doses throughout the organ assuming a given arrangement of one or more neutron beams, said computing step using the three-dimensional model of the organ to compute said three-dimensional array of expected radiation doses.

2. The method of claim 1 wherein the information about variations in internal composition of the organ is represented in said model as variations in radiation transport properties in the organ.

3. The method of claim 2 wherein said diagnostic cross-sectional scanner is a CT scanner.

4. The method of claim 3 wherein the computing step uses Monte Carlo simulation techniques to compute said three-dimensional array of expected radiation doses.

5. The method of claim 4 wherein said neutron capture therapy uses boron to enhance neutron capture in selected regions of the organ, said method further comprising identifying on said image of the organ a region within the organ that is believed to contain said tumor and wherein said generating step models said region as having a higher boron concentration than the part of the organ outside said region.

6. The method of claim 4 wherein the three-dimensional model is a three-dimensional array of three-dimensional cells and wherein said model generating step further comprises assigning each of said cells a corresponding one of a plurality of predetermined materials based upon a CT image value for that cell.

7. The method of claim 6 wherein the cells are cubes of uniform dimension,

8. The method of claim 6 further comprising:

selecting a cross-section of the organ for which dose information is desired;

computing isocontours for the selected cross-section from the computed three dimensional array of expected radiation doses; and

displaying the computed isocontours over a CT image of the selected cross-section.

9. The method of claim 4 further comprising:

selecting a cross-section of the organ for which dose information is desired;

setting the computed dose values outside of the organ to zero;

after setting the computed dose values outside the organ to zero, computing isocontours for the selected cross-section from the three-dimensional array of expected radiation doses; and

displaying the computing isocontours over a CT image of the selected cross-section.

10. The method of claim 4 wherein each of the cross-sectional images is a matrix of data elements, and said model generating step comprises:

assigning one of a plurality of predetermined materials to each of said data elements of said matrices, thereby generating an intermediate model; and

constructing the three-dimensional model from the intermediate model.

11. The method of claim 10 wherein the intermediate model is a three-dimensional array of subcells, each subcell corresponding to a different data element of the matrices of data elements, and wherein said three-dimensional model is an array of cells and wherein each cell is made up of a corresponding plurality of subcells.

12. The method of claim 11 wherein the model generating step further comprises computing for each given cell a percentage of each of said plurality of the predetermined materials that is assigned to the subcells that make up that given cell.

13. The method of claim 12 wherein the model generating step further comprises assigning to the given cell a composition that corresponds to the computed percentages of said plurality of predetermined materials that are assigned to the subcells that correspond to the given cell.

14. The method of claim 13 wherein the composition contains quantized amounts of said plurality of predetermined materials.

15. The method of claim 14 wherein the plurality of predetermined materials includes bone, soft tissue and air.

16. The method of claim 15 wherein the plurality of predetermined materials further includes tumor tissue.

17. The method of claim 16 wherein said neutron capture therapy uses boron to enhance neutron capture in selected regions of the organ, and wherein the soft tissue and the tumor are specified to contain different amounts of boron.

18. The method of claim 4 wherein the target organ is a brain.

19. A method of developing a treatment plan for boron neutron capture therapy of a patient'"'"'s brain, said method comprising:

using a diagnostic cross-sectional scanner to scan the patient'"'"'s head to generate cross-sectional images of the head;

generating a three-dimensional model of the head from the cross-sectional images, said model including information about variations in radiation transport properties inside the organ;

computing a three-dimensional array of expected radiation doses throughout the head assuming a given arrangement of one or more neutron beams, said computing step using the three-dimensional model of the head to compute said three-dimensional array of expected radiation doses and said computing step using Monte Carlo simulation techniques to compute said three-dimensional array of expected radiation doses;

selecting a cross-section of the patient'"'"'s head for which dose information is desired;

computing isocontours for the selected cross-section from the computed three dimensional array of expected radiation doses; and

displaying the computed isocontours over a CT image of the selected cross-section.