System and method for segmentation of lung

• US 9,836,848 B2
• Filed: 06/30/2015
• Issued: 12/05/2017
• Est. Priority Date: 07/02/2014
• Status: Active Grant

• Current Assignee: Covidien LP (Medtronic Plc)
• Original Assignee: Covidien LP (Medtronic Plc)
• Inventors: Markov, Igor A., Kreinin, Yuri
• Primary Examiner(s): Fitzpatrick, Atiba O
• Application Number:
  US14/754,867
• Publication Number:
  US 20160005193A1
• Time in Patent Office: 889 Days
• Field of Search: None
• CPC Class Codes: G06T 2207/10081 Computed x-ray tomography [CT]
• CPC Class Codes: G06T 2207/20101 Interactive definition of point of interest, landmark or seed
• CPC Class Codes: G06T 2207/20156 Automatic seed setting
• CPC Class Codes: G06T 2207/30061 Lung
• CPC Class Codes: G06T 7/0012 Biomedical image inspection
• CPC Class Codes: G06T 7/11 Region-based segmentation
• CPC Class Codes: G06T 7/13 Edge detection
• CPC Class Codes: G06T 7/136 involving thresholding
• CPC Class Codes: G06T 7/155 involving morphological operators
• CPC Class Codes: G06T 7/187 involving region growing; involving region merging; involving connected component labelling
• CPC Class Codes: G16H 30/20 for handling medical images, e.g. DICOM, HL7 or PACS
• CPC Class Codes: G16H 30/40 for processing medical images, e.g. editing

Abstract

Disclosed are systems, devices, and methods for determining pleura boundaries of a lung, an exemplary method comprising acquiring image data from an imaging device, generating a set of two-dimensional (2D) slice images based on the acquired image data, determining, by a processor, a seed voxel in a first slice image from the set of 2D slice images, applying, by the processor, a region growing process to the first slice image from the set of 2D slice images starting with the seed voxel using a threshold value, generating, by the processor, a set of binarized 2D slice images based on the region grown from the seed voxel, filtering out, by the processor, connected components of the lung in each slice image of the set of binarized 2D slice images, and identifying, by the processor, the pleural boundaries of the lung based on the set of binarized 2D slice images.

Claims

1. A segmentation method for determining pleural boundaries of a lung, comprising:

acquiring image data from an imaging device;

generating a set of two-dimensional (2D) slice images based on the acquired image data;

determining, by a processor, a seed voxel in a trachea connected to the lung in a first slice image from the set of 2D slice images;

applying, by the processor, a region growing process to the first slice image from the set of 2D slice images starting with the seed voxel using a threshold value, wherein a region grown by applying the region growing process corresponds to a lung;

generating, by the processor, a set of binarized 2D slice images based on the region grown from the seed voxel;

filtering out, by the processor, connected components of the lung in each slice image of the set of binarized 2D slice images; and

identifying, by the processor, the pleural boundaries of the lung based on the set of binarized 2D slice images.

2. The segmentation method according to claim 1, wherein the threshold value is greater than or equal to an intensity of the seed voxel.

3. The segmentation method according to claim 1, wherein the acquired image data is stored in the digital imaging and communications in medicine (DICOM) image format.

4. The segmentation method according to claim 1, wherein the image data is acquired via a network device.

5. The segmentation method according to claim 1, wherein applying the region growing process includes:

in a case where an intensity of a first voxel in the first slice image from the set of 2D slice images is lower than the predetermined threshold value and the first voxel is connected to the seed voxel, setting the intensity of the first voxel as a maximum value; and

in a case where an intensity of a second voxel in the first slice image from the set of 2D slice images is not lower than the predetermined threshold value or the first voxel is not connected to the seed voxel, setting the intensity of the second voxel as a minimum value.

6. The segmentation method according to claim 5, wherein the threshold value causes a high intensity area to appear around the seed voxel in the set of 2D slice images.

7. The segmentation method according to claim 5, wherein applying the region growing process further includes inversely assigning values of voxels in the set of 2D slice images, from the minimum value to the maximum value and from the maximum value to the minimum value, to obtain the set of binarized 2D slice images.

8. The segmentation method according to claim 1, wherein filtering out connected components of the lung includes:

detecting a plurality of connected components in the set of binarized 2D slice images;

calculating an area of each connected component in the set of binarized 2D slice images;

determining whether the area of each connected component is less than a predetermined value;

assigning the minimum value to pixels of a first connected component of the plurality of connected components when it is determined that an area of the first connected component is less than the predetermined value; and

assigning the maximum value to pixels of a second connected component of the plurality of connected components when it is determined that an area of the second connected component is greater than or equal to the predetermined value.

9. The segmentation method according to claim 8, wherein each connected component is an enclosed area with high intensity.

10. The segmentation method according to claim 8, wherein one or more of the plurality of connected components are blood vessels or airways.

11. The segmentation method according to claim 8, wherein an intersection of three 2D slice images, each of which is from each of three independent directions, identifies a voxel in the set of 2D slice images.

12. The segmentation method according to claim 11, wherein the three independent directions are axial, coronal, and sagittal directions.

13. The segmentation method according to claim 1, wherein each voxel of the set of binarized 2D slice images has either high or low intensity.

14. The segmentation method according to claim 1, wherein the image data is acquired from one of computed tomographic technique, radiography, tomogram produced by a computerized axial tomography scan, magnetic resonance imaging, ultrasonography, contrast imaging, fluoroscopy, nuclear scans, and positron emission tomography.

15. A system for determining pleural boundaries of a lung, the system comprising:

an imaging device configured to image a chest of a patient to obtain image data; and

an image processing device including;

a memory configured to store data and processor-executable instructions; and

a processor configured to execute the processor-executable instructions to;

generate a set of two-dimensional (2D) slice images based on the acquired image data;

determine a seed voxel in a trachea connected to the lung in a first slice image from the set of 2D slice images;

apply a region growing process to the first slice image from the set of 2D slice images starting with the seed voxel using a threshold value, wherein a region grown by applying the region growing process corresponds to a lung;

generate a set of binarized 2D slice images based on the region grown from the seed voxel;

filter out connected components of the lung in each slice image of the set of binarized 2D slice images; and

identify the pleural boundaries of the lung based on the set of binarized 2D slice images.