Rapid and robust 3D/3D registration technique

• US 7,724,943 B2
• Filed: 04/19/2005
• Issued: 05/25/2010
• Est. Priority Date: 04/21/2004
• Status: Active Grant

• Current Assignee: Siemens Medical Solutions USA Incorporated (Siemens AG)
• Original Assignee: Siemens Medical Solutions USA Incorporated (Siemens AG)
• Inventors: Xu, Chenyang, Sauer, Frank, Khamene, Ali
• Primary Examiner(s): Bali; Vikkram
• Assistant Examiner(s): ENTEZARI, MICHELLE M
• Application Number:
  US11/109,119
• Publication Number:
  US 20050249398A1
• Time in Patent Office: 1,862 Days
• Field of Search: None
• US Class Current: 382/154
• CPC Class Codes: A61P 9/00 Drugs for disorders of the cardiovascular system
• CPC Class Codes: G06T 2207/30004 Biomedical image processing
• CPC Class Codes: G06T 7/38 Registration of image sequences

Abstract

Exemplary methods are provided. In one exemplary method, a lower-dimension signature is generated for each of a plurality of data sets of a given dimension. Registration is performed on the lower-dimension signatures. In another exemplary method, a two-dimensional signature is generated for each of a plurality of three-dimensional volumes. Registration is performed on the two-dimensional signatures.

Claims

1. A method of registering a first volume with a second and different second volume, the method comprising:

generating by a processor three two-dimensional (2D) signatures for the first volume where each corresponds to a different one of three orthogonal directions;

generating by the processor three 2D signatures for the second volume where each corresponds to a different one of the three orthogonal directions;

registering by the processor the 2D signatures of the first volume with the 2D signatures of the second volume for each same orthogonal direction to generate three transformations; and

registering by the processor the first volume with the second volume using the three generated transformations,

wherein each pixel position of each 2D signature includes a minimum voxel intensity and maximum voxel intensity of a ray intersecting the corresponding 3D volume through the corresponding orthogonal direction and each 2D signature further includes rotation information and translation information of the corresponding 3D volume.

2. The method of claim 1, wherein generating a two-dimensional signature comprises:

generating a planar projective of the corresponding volume at the corresponding orthogonal direction.

3. The method of claim 1, wherein generating a two-dimensional signature comprises:

generating a multi-planar reconstruction of the corresponding volume at the corresponding orthogonal direction.

4. The method of claim 1, wherein generating a two dimensional signature comprises:

generating a non-planar reconstruction of the corresponding volume at the corresponding orthogonal direction.

5. The method of claim 1, wherein the two-dimensional signature are generated using a graphics processing unit (GPU).

6. The method of claim 5, wherein the two-dimensional signatures are generated by steering parallel rays in a transformation space of search parameters using a GPU-accelerated technique.

7. The method of claim 5, wherein the two-dimensional signatures are generated by steering parallel rays in a transformation space of search parameters using an orthographic volume rendering technique on the GPU.

8. The method of claim 1, further comprising:

segmenting the two-dimensional signatures for isolating an area of interest; and

performing selective matching in the area of interest.

9. The method of claim 1, further comprising:

segmenting the two-dimensional signatures for isolating a plurality of areas; and

performing piecewise rigid matching on the plurality of areas.

10. The method of claim 1, further comprising:

segmenting the two-dimensional signatures for isolating a plurality of areas; and

performing non-rigid matching on the plurality of areas.

11. The method of claim 1, further comprising:

simultaneously optimizing search parameters for each of the two-dimensional signatures.

12. The method of claims 11, further comprising:

iteratively optimizing the search parameters for each of the two-dimensional signatures.

13. The method of claims 1, wherein registering a 2D signatures of a first volume with a 2D signatures of a second volume, comprises:

calculating a similarity measure on the two-dimensional signatures.

14. A method of registering a first volume with a second and different volume, the method comprising:

generating by a processor three two-dimensional (2D) signatures for the first volume where each corresponds to a different one of three orthogonal directions;

generating by the processor three 2D signatures for the second volume where each corresponds to a different one of the three orthogonal directions;

registering by the processor the 2D signatures of the first volume with the 2D signatures of the second volume for each same orthogonal direction to generate three transformations; and

registering by the processor the first volume with the second volume using the three generated transformations,

wherein each pixel position of each 2D signature includes a minimum voxel intensity and maximum voxel intensity of a ray intersecting the corresponding 3D volume through the corresponding orthogonal direction.

15. A method of registering a first volume with a second and different volume, the method comprising:

generating by the processor three two-dimensional (2D) signatures for the first volume where each corresponds to a different one of three orthogonal directions;

generating by the processor three 2D signatures for the second volume where each corresponds to a different one of the three orthogonal directions;

registering by the processor the 2D signatures of the first volume with the 2D signatures of the second volume for each same orthogonal direction to generate three transformations; and

registering by the processor the first volume with second volume using the three generated transformations,

wherein each pixel position of each 2D signature includes a voxel intensity of a ray intersecting the corresponding 3D volume through the corresponding orthogonal direction.