Automatic organ-dose-estimation for patient-specific computed tomography scans

Patent metadata

Field	Value
Publication number	US-10098606-B2
Application number	US-201615055647-A
Country	US
Kind code	B2
Filing date	Feb 29, 2016
Priority date	Feb 29, 2016
Publication date	Oct 16, 2018
Grant date	Oct 16, 2018

Abstract

In accordance with at least some embodiments of the present disclosure, a process for calculating patient-specific organ dose is presented. The process may include constructing a computed tomography (CT) volume based on CT images generated by a CT scanner. The process may include segmenting the CT volume into a plurality of organ regions, generating a material density map for the CT volume based on Hounsfield Unit (HU) values, and generating a dose distribution map for the CT volume based on the material density map by simulating particles emitting from the CT scanner and flowing through the CT volume. The process may further generate a dose value delivered to a specific organ region of the plurality of organ regions based on the dose distribution map.

First claim

We claim:

1. A method for calculating a patient-specific organ dose, comprising:
constructing a computed tomography (CT) volume based on CT images generated by a CT scanner;
segmenting the CT volume into a plurality of organ regions, wherein the CT volume contains a plurality of voxels, and each of the plurality of organ regions contains a subset of voxels selected from the plurality of voxels;
generating a material density map for the CT volume based on the plurality of voxels' Hounsfield Unit (HU) values, wherein the material density map defines a corresponding material type and density for each of the plurality of voxels;
generating a dose distribution map for the CT volume based on the material density map by simulating the CT scanner based on the CT scanner's configuration, simulating macroscopic behavior of particles emitting from the CT scanner and flowing through the plurality of voxels, and constructing a Boltzmann Transport Equation (BTE) using the macroscopic behavior of the particles and the material density map, wherein the dose distribution map stores a plurality of three-dimensional (3D) particle fluence distribution values among the plurality of voxels; and
generating a dose value delivered to a specific organ region of the plurality of organ regions based on the dose distribution map.

2. The method as recited in claim 1 , wherein the segmenting of the CT volume into the plurality of organ regions comprising:
identifying the plurality of organ regions in the CT volume based on the plurality of voxels' HU values.

3. The method as recited in claim 1 , wherein the segmenting of the CT volume into the plurality of organ regions comprising:
identifying the plurality of organ regions in the CT volume based on an atlas which contains previously-segmented organ regions associated with a previous CT volume.

4. The method as recited in claim 3 , wherein the segmenting of the CT volume into the plurality of organ regions further comprising:
performing a deformable registration operation to map the previously-segmented organ regions in the previous CT volume to the plurality of organ regions in the CT volume.

5. The method as recited in claim 1 , wherein the simulating of the CT scanner comprising:
using a fluence grid associated with the CT scanner to provide a spatially-dependent particle fluence and energy spectra for the particles.

6. The method as recited in claim 1 , wherein the CT scanner's configuration accounts for beam energy (kVp), beam current (mA), exposure time, collimation, and/or beam filtering as a function of gantry angle and patient position.

7. The method as recited in claim 1 , wherein the simulating of the macroscopic behavior of particles emitting from the CT scanner comprising:
simulating the particles being absorbed or scattered in the plurality of voxels.

8. The method as recited in claim 7 , wherein the simulating of the macroscopic behavior of particles emitting from the CT scanner further comprising:
transporting the particles emitting from the CT scanner through the plurality of voxels to calculate a set of scattering sources.

9. The method as recited in claim 8 , wherein the simulating of the macroscopic behavior of particles emitting from the CT scanner further comprising:
transporting the particles from the set of scattering sources across the plurality of voxels to calculate angular flux in the plurality of voxels.

10. The method as recited in claim 9 , wherein the simulating of the macroscopic behavior of particles emitting from the CT scanner further comprising:
iterating through the transporting of the particles until a converged solution is obtained in calculating the angular flux in the plurality of voxels.

11. The method as recited in claim 1 , wherein the generating of the dose distribution map further comprising:
generating the dose distribution map by using a deterministic method to solve the BTE and calculate the 3D particle fluence distribution values among the plurality of voxels.

12. The method as recited in claim 11 , wherein the generating of the dose distribution map further comprising:
obtaining a 3D particle fluence distribution value associated with a specific voxel in the plurality of voxels from the deterministic method's solution; and
storing the 3D particle fluence distribution value associated with the specific voxel in the dose distribution map.

13. The method as recited in claim 1 , wherein the generating of the dose value delivered to the specific organ region comprising:
generating the organ dose value by accumulating corresponding dose values for the subset of voxels from the dose distribution map and based on the 3D particle fluence distribution values for the subset of voxels stored in the dose distribution map.

14. The method as recited in claim 1 , wherein the CT scanner is a cone-beam CT scanner and the CT volume is a cone-beam CT volume constructed based on a plurality of x-ray projections generated by the CT scanner.

15. The method as recited in claim 1 , wherein the particles emitting from the CT scanner are x-ray photons originated from a plurality of x-ray projection sources.

16. A method for calculating patient-specific organ dose, comprising:
segmenting a computed tomography (CT) volume originated from a CT scanner into a plurality of organ regions, wherein the CT volume contains a plurality of voxels, and each of the plurality of organ regions contains a subset of voxels selected from the plurality of voxels;
determining material types and densities for the plurality of voxels in the CT volume based on the plurality of voxels' Hounsfield Unit (HU) value;
generating a dose distribution map for the CT volume by simulating macroscopic behavior of particles emitting from the CT scanner and flowing through the plurality of voxels and performing a Boltzmann Transport Equation (BTE) calculation with a deterministic solver based on the macroscopic behavior of the particles, the CT scanner's configuration, and the plurality of voxels' material types and densities, wherein the dose distribution map stores a plurality of dose values associated with the plurality of voxels; and
generating a dose value for a specific organ region selected from the plurality of organ regions, based on the subset of voxels associated with the specific organ region and the subset of voxels' dose values stored in the dose distribution map.

17. The method as recited in claim 16 , wherein the segmenting of the CT volume into the plurality of organ regions comprising:
identifying the plurality of organ regions in the CT volume based on an atlas which contains previously-segmented organ regions associated with a previous CT volume.

18. The method as recited in claim 17 , wherein the segmenting of the CT volume into the plurality of organ regions further comprising:
performing a deformable registration operation to map the previously-segmented organ regions in the previous CT volume to the plurality of organ regions in the CT volume.

19. The method as recited in claim 16 , wherein the performing the BTE calculation with the deterministic solver comprising:
obtaining a 3D particle fluence distribution value associated with a specific voxel in the plurality of voxels from the deterministic solver's solution; and
storing the 3D particle fluence distribution value associated with the specific voxel in the dose distribution map.

20. A system configured to calculate patient-specific organ dose, comprising:
a memory system containing a set of ins

Assignees

• Varian Med Sys Inc

Inventors

• Star-Lack Josh
• Wang Adam
• Maslowski Alexander

Classifications

• A61B6/542 Primary

  involving control of exposure · CPC title

• G06T2207/30004

  Biomedical image processing · CPC title

• G06T7/11

  Region-based segmentation · CPC title

• A61B6/032

  Transmission computed tomography [CT] · CPC title

• G06T7/0012

  Biomedical image inspection · CPC title