Muon tomography imaging improvement using optimized limited angle data

• US 9,639,973 B2
• Filed: 04/03/2015
• Issued: 05/02/2017
• Est. Priority Date: 04/04/2014
• Status: Active Grant

• Current Assignee: Decision Sciences International Corporation
• Original Assignee: Decision Sciences International Corporation
• Inventors: Bai, Chuanyong, Kindem, Joel, Luo, Weidong, Steiger, Matthew, Simon, Sean, Sossong, Michael James
• Primary Examiner(s): Niu, Feng
• Application Number:
  US14/678,921
• Publication Number:
  US 20150287237A1
• Time in Patent Office: 760 Days
• Field of Search: None
• CPC Class Codes: G01N 9/00 Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
• CPC Class Codes: G01V 5/00 Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
• CPC Class Codes: G01V 5/005 using Tomography, e.g. CT or SPECT detector details in CT applications G01T1/2985
• CPC Class Codes: G01V 5/22 Active interrogation, i.e. by irradiating objects or goods using external radiation sources, e.g. using gamma rays or cosmic rays
• CPC Class Codes: G01V 5/226 using tomography
• CPC Class Codes: G06T 11/006 Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
• CPC Class Codes: G06T 11/60 Editing figures and text; Combining figures or text
• CPC Class Codes: G06T 15/08 Volume rendering

Abstract

Disclosed technology can provide a process for generating reconstructed muon image resolution to optimize the use of the limited angular range muon track data collected by a muon tomography system. In one aspect, a process for improving reconstructed muon image resolution for a volume of interest (VOI) imaged by a muon tomography system includes: collecting raw muon track data of cosmic ray muon tracks passing through the VOI; grouping the raw muon track data into two or more subsets of tracks based on at least one angular distribution of the muon tracks in the raw muon track data; generating a set of images of the VOI based on the two or more subsets of tracks; and combining information from the set of reconstructed images and a reconstructed image based on the full set of the raw muon track data to obtain a resulting reconstructed image of the VOI.

Claims

1. A method for generating charged particle tomography images for a volume of interest (VOI) by a charged particle tomography system, the method comprising:

collecting raw track data of tracks of cosmic-ray based charged particles passing through the VOI;

grouping the collected raw track data into subsets of the tracks based on angular distribution of the tracks in a zenith direction in the raw track data by partitioning the collected raw track data into a first subset corresponding to large zenith angles and a second subset corresponding to small zenith angles;

generating a set of images of the VOI based on the grouped subsets of the tracks; and

combining information from at least the generated set of images to obtain an resulting image of the VOI.

2. The method of claim 1, comprising:

generating an image based on full set of the raw track data; and

wherein the combining includes combining information from the generated image based on the full set of the raw track data along with the information from the generated set of images to obtain the resulting image of the VOI.

3. The method of claim 1, wherein the first subset of tracks has a zenith angle range from 0 to X and the second subset of tracks has a zenith angle range from X to 90°

, wherein X is between 20°

to 45°

to achieve a distribution of muon counts in the first and second subsets that limits a noise level of the set of images.

4. The method of claim 3, comprising:

selecting X so that percentage of the tracks in the first subset and percentage of the tracks in the second subset have a predetermined relationship.

5. The method of claim 4, wherein the predetermined relationship includes the two percentages being substantially equal to each other.

6. The method of claim 3, wherein a first reconstructed image of the set of images based on the first subset of the tracks is associated with an object resolution in a horizontal direction parallel to a detector plane in the charged particle tomography system;

and

wherein a second reconstructed image based on the second subset of the tracks is associated with an object resolution in a vertical direction perpendicular to a detector plane in the charged particle tomography system.

7. The method of claim 6, wherein a third reconstructed image based on the full set of the raw track data is associated with an object resolution in the horizontal direction which is at a level between the object resolutions in the horizontal direction for the first reconstructed image and the second reconstructed image;

and

wherein the third reconstructed image based on the full set of the raw track data is associated with an object resolution in the vertical direction at a level between the resolutions in the vertical direction for the first reconstructed image and the second reconstructed image.

8. The method of claim 1, wherein grouping the raw track data into the subsets of the tracks based on the angular distribution includes:

identifying a subset of the tracks which are substantially perpendicular to a direction û

in a detector plane of the charged particle tomography system.

9. The method of claim 8, wherein an image based on the identified subset of the tracks substantially perpendicular to the direction û

is associated with an object resolution in the û

-direction.

10. The method of claim 1, wherein grouping the raw track data into the subsets of the tracks based on the angular distribution includes:

partitioning the collected raw track data into a third subset corresponding to azimuthal angles, wherein the large zenith angles, the small zenith angles, and the azimuthal angles for the first, second, and third subsets respectively are within 45°

of an x-direction in a detector plane of the charged particle tomography system; and

partitioning the raw track data into a fourth subset corresponding to azimuthal angles which are within 45°

of a y-direction in a detector plane of the charged particle tomography system.

11. The method of claim 1, wherein at least two of the subsets of the tracks share a common cosmic-ray based charged particle track.

12. The method of claim 1, wherein the cosmic-ray based charge particles include muons.

13. A method for generating charged particle tomography image for a volume of interest (VOI) by a charged particle tomography system, the method comprising:

collecting raw track data of tracks of cosmic-ray charged particles passing through the VOI;

generating an image of the VOI based on the raw track data;

identifying an orientation of objects in the generated image of the VOI based on the raw track data;

grouping the raw track data into subsets of the tracks based on the identified orientation, wherein grouping the raw track data into the subsets of the tracks includes partitioning the raw track data into a first subset corresponding to large zenith angles and a second subset corresponding to small zenith angles;

generating a set of images of the VOI based on the subsets of the tracks; and

combining information from the generated set of images of the VOI based on the subsets of the tracks and the generated image based on the full set of the track data to obtain a resulting image of the VOI associated with a resolution along the identified orientation.

14. The method of claim 13, wherein grouping the raw track data into the subsets of the tracks based on the identified orientation includes identifying a subset of the tracks substantially perpendicular to the identified orientation.

15. The method of claim 14, wherein an image of the set of images based on the identified subset of tracks substantially perpendicular to the identified orientation is associated with an object resolution in the identified orientation.

16. The method of claim 13, wherein the first subset of the tracks has a zenith angle range from 0 to X and the second subset of the tracks has a zenith angle range from X to 90°

, wherein X is between 20°

or 45°

.

17. The method of claim 13, wherein a first image based on the first subset of tracks provides an improvement of an object resolution in a horizontal direction parallel to a detector plane in the charged particle tomography system;

and wherein a second image based on the second subset of tracks provides an improvement of an object resolution in a vertical direction perpendicular to a detector plane in the charged particle tomography system.

18. The method of claim 13, wherein grouping the raw track data into the subsets of the tracks includes partitioning the raw track data into a first subset corresponding to azimuthal angles in a first range which are within 45°

of an x-direction in a detector plane of the charged particle tomography system; and

a second subset corresponding to azimuthal angles in a second range which are within 45°

of a y-direction in the detector plane of the charged particle tomography system.

19. The method of claim 18, wherein a first reconstructed image based on the first subset of tracks is associated with an object resolution in the y-direction;

and

wherein a second reconstructed image based on the second subset of tracks is associated with an object resolution in the x-direction.

20. The method of claim 13, wherein the cosmic-ray charged particles include muons.