Creating a blood vessel tree from imaging data

• US 7,983,459 B2
• Filed: 11/22/2006
• Issued: 07/19/2011
• Est. Priority Date: 10/25/2006
• Status: Expired due to Fees

• Current Assignee: Rcadia Medical Imaging, Ltd.
• Original Assignee: Rcadia Medical Imaging, Ltd.
• Inventors: Begelman, Grigory, Walach, Eugene, Levanon, Shai, Ohayon, Shay, Goldenberg, Roman
• Primary Examiner(s): Perungavoor; Sath V
• Assistant Examiner(s): HEIDEMANN, JASON E
• Application Number:
  US11/562,875
• Publication Number:
  US 20080101674A1
• Time in Patent Office: 1,700 Days
• Field of Search: None
• US Class Current: 382/128
• CPC Class Codes: G06T 2207/30048 Heart; Cardiac
• CPC Class Codes: G06T 7/0012 Biomedical image inspection
• CPC Class Codes: G06T 7/12 Edge-based segmentation
• CPC Class Codes: G06V 10/255 Detecting or recognising potential candidate objects based on visual cues, e.g. shapes
• CPC Class Codes: G06V 2201/031 of internal organs

Abstract

A method of defining a heart region from imaging data is provided. Received imaging data is projected into a first plane. A first threshold is applied to the first plane of data to eliminate data associated with air. A largest first connected component is identified from the first threshold applied data. A first center of mass of the identified largest first connected component is calculated to define a first coordinate and a second coordinate of the heart region. The received imaging data is projected into a second plane, wherein the second plane is perpendicular to the first plane. A second threshold is applied to the second plane of data to eliminate data associated with air. A largest second connected component is identified from the second threshold applied data. A second center of mass of the identified largest second connected component is calculated to define a third coordinate of the heart region.

Claims

1. A method of creating a blood vessel tree from imaging data, the method comprising:

(a) receiving a blood vessel object identified from imaging data;

(b) selecting a starting point of the received blood vessel object;

(c) calculating a width map for the received blood vessel object;

(d) identifying an exit point for the received blood vessel object;

(e) calculating a distance map relative to the identified exit point for the received candidate vessel object, wherein the distance map is weighted by the width values of the width map;

(f) identifying an endpoint during calculation of the distance map;

(g) backtracking from the identified endpoint to the selected starting point to define a path in a vessel tree;

(h) repeating (g) for each identified endpoint to create a vessel tree for the blood vessel object; and

filtering the identified endpoint if the identified endpoint is not a local maxima based on a distance to the selected starting point.

2. The method of claim 1, wherein the blood vessel object is identified as a right coronary artery or as a left main artery.

3. The method of claim 1, wherein the endpoint is a voxel that did not update any neighbor voxel during calculation of the distance map.

4. The method of claim 1, further comprising:

determining if the path in the vessel tree includes a voxel in a second path;

if the path includes a voxel in the second path, defining a bifurcation vertex; and

connecting the path with the second path.

5. The method of claim 1, further comprising:

receiving imaging data including a heart region, a pulmonary artery region, an ascending aorta region, and a thoraco-abdominal aorta region of a patient;

identifying the heart region from the imaging data;

analyzing the identified heart region to detect a coronary pathology;

identifying an ascending aorta object from the ascending aorta region of the imaging data;

analyzing the identified ascending aorta object to detect an aortic dissection;

identifying an abdominal aorta object from the thoraco-abdominal aorta region;

analyzing the identified abdominal aorta object to detect an aortic dissection;

identifying a left main pulmonary artery object and a right main pulmonary artery object from the pulmonary artery region of the imaging data;

analyzing the identified left main pulmonary artery object and the identified right main pulmonary artery object to detect a pulmonary embolism; and

generating a report including any detected coronary pathology, any detected aortic dissection, and any detected pulmonary embolism.

6. A method of creating a blood vessel tree from imaging data, the method comprising:

(a) receiving a blood vessel object identified from imaging data;

(b) selecting a starting point of the received blood vessel object;

(c) calculating a width map for the received blood vessel object;

(d) identifying an exit point for the received blood vessel object;

(e) calculating a distance map relative to the identified exit point for the received candidate vessel object, wherein the distance map is weighted by the width values of the width map;

(f) identifying an endpoint during calculation of the distance map;

(g) backtracking from the identified endpoint to the selected starting point to define a path in a vessel tree; and

(h) repeating (g) for each identified endpoint to create a vessel tree for the blood vessel object;

determining if the identified endpoint is a disconnected vessel;

if the identified endpoint is a disconnected vessel, estimating a direction of the disconnected vessel at the identified endpoint;

detecting a connected component in a portion of the imaging data in the estimated direction of the disconnected vessel; and

repeating (g) for the detected connected component.

7. The method of claim 6, further comprising removing a path having a length less than a threshold value.

8. The method of claim 7, wherein the threshold is approximately five millimeters.

9. The method of claim 6, further comprising:

identifying an x-junction in the created vessel tree; and

removing a path associated with the identified x-junction.

10. The method of claim 9, wherein identifying the x-junction comprises identifying two bifurcations located within a distance threshold relative to the other.

11. The method of claim 10, wherein the distance threshold is approximately three millimeters.

12. The method of claim 10, wherein identifying the x-junction further comprises identifying two bifurcations intersecting at an angle between approximately 70 degrees and approximately 110 degrees.

13. A method of creating a blood vessel tree from imaging data, the method comprising:

(a) receiving a blood vessel object identified from imaging data;

(b) selecting a starting point of the received blood vessel object;

(c) calculating a width map for the received blood vessel object;

(d) identifying an exit point for the received blood vessel object;

(e) calculating a distance map relative to the identified exit point for the received candidate vessel object, wherein the distance map is weighted by the width values of the width map;

(f) identifying an endpoint during calculation of the distance map;

(g) backtracking from the identified endpoint to the selected starting point to define a path in a vessel tree; and

(h) repeating (g) for each identified endpoint to create a vessel tree for the blood vessel object;

identifying an x-junction in the created vessel tree, wherein identifying the x-junction comprises identifying two bifurcations located within a distance threshold relative to the other;

removing a path associated with the identified x-junction;

identifying a first path of a first of the identified two bifurcations, wherein the first path has a closest direction to an edge arriving from an aorta object;

identifying a second path of the first of the identified two bifurcations, wherein the second path does not have the closest direction to the edge arriving from the aorta object and

selecting the path as the second path.