SYSTEM AND METHOD FOR DIGITAL VOLUME PROCESSING WITH GPU ACCELERATIONS
A method for processing digital volume image data, executed at least in part on a computer that has an associated graphics processing unit, orders the digital volume image data as a stack of two-dimensional image slices of voxels. A 1:1 mapping of each of the slices is formed, in sequential order, to a corresponding tile in a digital one-row flat volume in a global memory. For each voxel in the digital flat volume, a neighborhood is defined that has, relative to each voxel, four or more adjacent voxels that are within the corresponding tile of each voxel, and one or more adjacent voxels that are within the preceding tile, and one or more adjacent voxels that are within the next tile. Each voxel is updated according to information stored for the adjacent voxels in the defined neighborhood for said each voxel to form a rendered volume image, which is displayed.
1 . A method for processing digital volume image data, the method executed at least in part on a computer that has an associated graphics processing unit, comprising:
ordering the digital volume image data as a stack of two-dimensional image slices, each image slice containing a plurality of voxels;
forming a 1:1 mapping of each of the slices, in sequential order, to a corresponding tile in a digital one-row flat volume in a global memory of the graphics processing unit;
defining, for each voxel of a plurality of voxels in the digital flat volume, a neighborhood that comprises, relative to said each voxel, four or more adjacent voxels that are within the corresponding tile of said each voxel, and one or more adjacent voxels that are within the preceding tile, and one or more adjacent voxels that are within the next tile;
updating said each voxel according to information stored for the adjacent voxels in the defined neighborhood for said each voxel to form a rendered volume image; and
displaying the updated rendered volume image.
2 . The method of claim 1 wherein defining the neighborhood comprises addressing each adjacent voxel using an x-y offset that gives the position of the adjacent voxel within the corresponding tile.
3 . The method of claim 1 wherein the digital flat volume is a one-row tile flat volume.
4 . The method of claim 1 further comprising identifying one or more boundary voxels in at least one tile that are not updated according to the defined neighborhood.
5 . The method of claim 1 wherein the digital flat volume is a texture of the graphics processing unit.
6 . The method of claim 1 wherein forming the rendered volume image further comprises performing image segmentation on the flat volume image.
7 . The method of claim 1 wherein at least two voxels are updated simultaneously.
8 . The method of claim 1 wherein the digital volume image is a computed tomography image.
9 . The method of claim 1 wherein updating comprises applying a filter to the digital volume image data.
10 . The method of claim 1 wherein updating comprises applying segmentation to the digital volume image data.
11 . A method for processing digital volume image data, the method executed at least in part on a computer that has an associated graphics processing unit, comprising:
ordering the digital volume image data as a stack of two-dimensional image slices, each slice containing a plurality of voxels;
forming a 1:1 mapping of each of the slices, in sequential order, to a corresponding tile in a digital one-row flat volume in a global memory of the graphics processing unit;
defining, for each voxel of a plurality of voxels in the digital flat volume, a neighborhood that comprises, relative to said each voxel, eight adjacent voxels that are within the corresponding tile of said each voxel, nine adjacent voxels that are within the preceding tile, and nine adjacent voxels that are within the next tile;
updating said each voxel according to information stored for the adjacent voxels in the defined neighborhood for said each voxel to form a rendered volume image; and
displaying the updated rendered volume image.
12 . The method of claim 11 wherein defining the neighborhood comprises addressing each adjacent voxel using an x-y offset that gives the position of the adjacent voxel within the corresponding tile.
13 . The method of claim 11 wherein the digital flat volume is a one-row tile flat volume.
14 . The method of claim 11 further comprising identifying one or more boundary voxels in at least one tile that are not updated according to the defined neighborhood.
15 . The method of claim 11 wherein the digital flat volume is a texture of the graphics processing unit.
16 . The method of claim 11 wherein forming the rendered volume image further comprises performing image segmentation on the flat volume image.
17 . The method of claim 11 wherein at least two voxels are updated simultaneously.
18 . The method of claim 11 wherein the digital volume image is a computed tomography image.
19 . The method of claim 11 wherein updating comprises applying a filter to the digital volume image data.