System and method for segmentation of lung

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.

1 - 16 . (canceled) 
     
     
         17 . A system for determining pleural boundaries of a lung, comprising:
 a processor configured to:
 apply a region growing process to a first slice image of a set of 2D slice images to form a grown region, the region growing process starting with a seed voxel within a trachea and using a threshold value, wherein 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 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 threshold value or the first voxel is not connected to the seed voxel, setting the intensity of the second voxel as a minimum value; 
 
 filter connected components of a lung in each slice image of a set of binarized 2D slice images generated from the grown region; and 
 identify the pleural boundary of the lung based on the set of binarized 2D slice images; and 
   a display operably coupled to the processor and configured to display a three-dimensional rendering or a two-dimensional image of the lung.   
     
     
         18 . The system according to  claim 17 , wherein the threshold value is greater than or equal to an intensity of the seed voxel. 
     
     
         19 . The system according to  claim 17 , wherein the set of 2D slice images is stored in digital imaging and communications in medicine (DICOM) image format. 
     
     
         20 . The system according to  claim 17 , wherein the processor is further configured to:
 generate the set of 2D slice images based on acquired image data of a patient's chest; and   determine the seed voxel in the trachea connected to the lung in the first slice image of the set of 2D slice images.   
     
     
         21 . The system according to  claim 17 , wherein the threshold value causes a high intensity area to appear around the seed voxel in the set of 2D slice images. 
     
     
         22 . The system according to  claim 17 , wherein the processor is configured to filter connected components of the lung by:
 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 a 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 a 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.   
     
     
         23 . The system according to  claim 22 , wherein each connected component is an enclosed area with high intensity. 
     
     
         24 . The system according to  claim 22 , wherein one or more of the plurality of connected components are blood vessels or airways. 
     
     
         25 . The system according to  claim 22 , 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. 
     
     
         26 . The system according to  claim 25 , wherein the three independent directions are axial, coronal, and sagittal directions. 
     
     
         27 . The system according to  claim 17 , wherein each voxel of the set of binarized 2D slice images has either high or low intensity. 
     
     
         28 . The system according to  claim 17 , wherein the set of 2D slice images is acquired from at least 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. 
     
     
         29 . A non-transitory computer-readable storage medium storing instructions which, when executed by a processor, cause the processor to:
 apply a region growing process to a first slice image of a set of 2D slice images to form a grown region, the region growing process starting with a seed voxel within a trachea and using a threshold value, wherein 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 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 threshold value or the first voxel is not connected to the seed voxel, setting the intensity of the second voxel as a minimum value; 
   filter connected components of a lung in each slice image of a set of binarized 2D slice images generated from the grown region; and   identify a pleural boundary of the lung based on the set of binarized 2D slice images.   
     
     
         30 . The non-transitory computer-readable storage medium according to  claim 29 , wherein the threshold value is greater than or equal to an intensity of the seed voxel. 
     
     
         31 . The non-transitory computer-readable storage medium according to  claim 29 , wherein the set of 2D slice images is stored in digital imaging and communications in medicine (DICOM) image format. 
     
     
         32 . The non-transitory computer-readable storage medium according to  claim 29 , wherein the instructions, when executed, further cause the processor to:
 generate the set of 2D slice images based on acquired image data of a patient's chest; and   determine the seed voxel in the trachea connected to the lung in the first slice image of the set of 2D slice images.   
     
     
         33 . A system for determining pleural boundaries of a lung, comprising:
 a processor configured to:
 form a grown region by:
 setting an intensity of a first voxel in a first slice image of a set of 2D slice images as a maximum value in a case where an intensity of the first voxel in the first slice image from the set of 2D slice images is lower than a threshold value and the first voxel is connected to a seed voxel within a trachea; and 
 setting the intensity of a second voxel in the first slice image from the set of 2D slice images as a minimum value in a case where an intensity of the second voxel in the first slice image from the set of 2D slice images is not lower than the threshold value or the first voxel is not connected to the seed voxel; 
 
 filter connected components of a lung in each slice image of a set of binarized 2D slice images generated from the grown region; and 
 identify the pleural boundary of the lung based on the set of binarized 2D slice images. 
   
     
     
         34 . The system according to  claim 33 , wherein the threshold value is greater than or equal to an intensity of the seed voxel. 
     
     
         35 . The system according to  claim 33 , wherein the processor is further configured to:
 generate the set of 2D slice images based on acquired image data of a patient's chest; and   determine the seed voxel in the trachea connected to the lung in the first slice image of the set of 2D slice images.   
     
     
         36 . The system according to  claim 33 , wherein the processor is further configured to inversely assign 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.