The present disclosure relates to systems and methods for segmentation of a lung. More particularly, the present disclosure relates to systems and methods that defines the borders of the lung based on a three-dimensional (3D) model generated based on CT scan image data of a patient's chest.
A patient's lungs are located within a thoracic cavity, which is separated from the abdominal cavity by the muscular diaphragm located at the base of the lungs. Further, the lungs are surrounded by a double-walled sac called the pleura (visceral pleurae and parietal pleurae) and a pleural fluid between the pleura and the lungs. The pleural fluid enables the lungs to expand and contract without adhering the pleura.
Visualization techniques related to visualizing the lungs have been developed so as to help clinicians to perform diagnoses, surgeries and/or other treatment procedures. Visualization is especially important for identifying a location of a diseased region having corresponding symptoms. Further, when treating the diseased region, additional emphasis is given to identification of a correct location so that a proper procedure is performed at the correct location. Furthermore, visualization of borders of the lungs are critical because diseased regions and locations for surgical operations should be within the borders of the lungs. Thus, visualization of pleura is important because the pleura define the borders of the lungs.
Provided in accordance with the present disclosure is a method for determining pleura boundaries of a lung.
In an aspect of the present disclosure, a segmentation method comprises 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.
In another aspect of the present disclosure, the seed voxel is in a portion of the first slice image from the set of binarized 2D slice images corresponding to a trachea of the lung.
In yet another aspect of the present disclosure, the threshold value is greater than or equal to an intensity of the seed voxel.
In a further aspect of the present disclosure, the acquired image data is stored in the digital imaging and communications in medicine (DICOM) image format.
In another aspect of the present disclosure, the image data is acquired via a network device.
In a further aspect of the present disclosure, 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.
In another aspect of the present disclosure, the threshold value causes a high intensity area to appear around the seed voxel in the set of 2D slice images.
In a further aspect of the present disclosure, 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.
In another aspect of the present disclosure, filtering out connected components of the lung includes detecting a connected component 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 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 connected component when it is determined that an area of the second connected component is greater than or equal to the predetermined value.
In a further aspect of the present disclosure, a connected component is an enclosed area with high intensity.
In another aspect of the present disclosure, the connected component is a blood vessel or an airway.
In a further aspect of the present disclosure, 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.
In another aspect of the present disclosure, the three independent directions are axial, coronal, and sagittal directions.
In a further aspect of the present disclosure, each voxel of the set of binarized 2D slice images has either high or low intensity.
In another aspect of the present disclosure, the image data is acquired from 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.
In another aspect of the present disclosure, a system for determining pleura of a lung comprises 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 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, 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.
Any of the above aspects and embodiments of the present disclosure may be combined without departing from the scope of the present disclosure.
Objects and features of the presently disclosed systems and methods will become apparent to those of ordinary skill in the art when descriptions of various embodiments are read with reference to the accompanying drawings, of which:
The present disclosure is related to systems and methods for segmenting image data of a patient's chest to identify the pleural boundaries of the patient's lungs.
Segmenting image data of a patient's chest to identify the pleural boundaries of the patient's lungs may be a necessary component of an ELECTROMAGNETIC NAVIGATION BRONCHOSCOPY® (ENB) procedure using an electromagnetic navigation (EMN) system. An ENB procedure generally involves at least two phases: (1) planning a pathway to a target located within, or adjacent to, the patient's lungs; and (2) navigating a probe to the target along the planned pathway. An example of the planning software described herein can be found in U.S. patent application Ser. Nos. 13/838,805, 13/838,997, and 13/839,224, all of which are filed by Covidien LP on Mar. 15, 2013 and entitled “Pathway Planning System and Method,” all of which are incorporated herein by reference. An example of the planning software can be found in commonly assigned U.S. Provision Patent Application No. 62/020,240 entitled “SYSTEM AND METHOD FOR NAVIGATING WITHIN THE LUNG” the entire contents of which are incorporated herein by reference.
Prior to the planning phase, the patient's lungs are imaged by, for example, a computed tomography (CT) scan, although additional applicable methods of imaging will be known to those skilled in the art. The image data assembled during the CT scan may then be stored in, for example, the Digital Imaging and Communications in Medicine (DICOM) format, although additional applicable formats will be known to those skilled in the art. The CT scan image data may then be loaded into a planning software application (“application”) to be processed for generating a 3D model which may be used during the planning phase of the ENB procedure. Segmenting image data of the patient's chest may be a component of the process of generating the 3D model, or may be performed separately.
The application may use the CT scan image data to generate a 3D model of the patient's lungs. The 3D model may include, among other things, a model airway tree corresponding to the actual airways of the patient's lungs, and showing the various passages, branches, and bifurcations of the patient's actual airway tree. While the CT scan image data may have gaps, omissions, and/or other imperfections included in the image data, the 3D model is a smooth representation of the patient's airways, with any such gaps, omissions, and/or imperfections in the CT scan image data filled in or corrected. After it is generated, the 3D model may be presented in various views. Although the present disclosure will be described in terms of specific illustrative embodiments, it will be readily apparent to those skilled in this art that various modifications, rearrangements and substitutions may be made without departing from the spirit of the present disclosure. The scope of the present disclosure is defined by the claims appended hereto.
The CT scan image data may be stored in the memory 120 in the DICOM format. A set of 2D slice images may be generated based on the CT scan image data. In an aspect, the 2D slice images may be generated to depict the axial, coronal, and sagittal views of the patient at a given location. For example, at each intersecting point of the chest, there may be three different 2D slice images generated in the three independent directions. These 2D slice images may be loaded into application 122 which is executed by processors 110 to reformat the 2D slice images for display. For example, application 122 may convert a color space of the 2D slice images to another color space suitable for display and perform imaging processes, e.g., scale, rotation, translation, or projection, to display the 2D slice images as intended. Based on at least three different 2D slice images, a voxel may be identified. The 2D slice images may be processed using a region growing algorithm to distinguish different types of tissue and materials to identify the pleural boundaries of the lungs. A variety of region growing algorithms are known in the art. Memory 120 may also store applications that may generate the 3D model of the chest or that may identify the pleural boundaries of the lungs.
In one embodiment, segmentation may be used to distinguish an inside area of the lungs from an outside area of the lungs in the set of 2D slice images. Segmentation may include binarization where the 2D slice images are displayed with only two different intensity values, and filtering processes as will be described in greater detail below.
As a first step of the present disclosure, a set of 2D slice images is generated from CT scan image data of the patient's chest. Once the set of 2D slice images is generated, a seed voxel 220 must be selected.
The trachea has a substantially constant diameter along its length. In particular, the upper part of the trachea extends substantially linearly from the larynx and pharynx and behind the sternum. By using these properties of the trachea, the trachea may be found in the 2D slice images. Further details of finding the trachea to determine the starting seed voxel 220 is described in commonly assigned U.S. Provisional Patent Application No. 62/020,253 entitled “Trachea Marking”, filed on Jul. 2, 2014 by Lachmanovich et al.; and U.S. Provisional Patent Application No. 62/020,257 entitled “Automatic Detection Of Human Lung Trachea”, filed on Jul. 2, 2014 by Markov et al.
After the starting seed voxel 220 is selected either manually or automatically, a region growing algorithm is used to binarize the 2D slice images. Based on the region growing algorithm, every voxel in the set of 2D slice images is checked to determine whether a Hounsfield value assigned to each voxel is less than a threshold value and whether each voxel is connected to the starting voxel 220. When it is determined that a value assigned to a voxel has a Hounsfield value less than the threshold value and is connected to the starting seed voxel 220, the Hounsfield value of the voxel is set to one or the maximum value. Otherwise, the Hounsfield value of the voxel is set to zero or the minimum value. As part of the region growing algorithm, the threshold is selected with a high enough value to cause leakage in the lung, and thus fill the lungs with intensity values leaked from the airways.
After every voxel in the set of 2D slice images is set to the maximum or minimum value, the 2D slice images will have only 2 colors of pixels. The result is a set of 2D slice images where the voxels/pixels having the maximum Hounsfield value would appear white, and the voxels/pixels having the minimum Hounsfield value would appear black. As shown in
Since the lung has a porous structure, most of voxels of the binarized set of 2D slice images representing the lung and the airways are shown in black. Some airway tissue (e.g., 355 in
The segmented, binarized, and inversed set of 2D slice images may contain large black regions, which cover the lung. In an aspect, the threshold value may be adjusted to make the large black regions as large as possible to cover the pleural boundaries of the lungs with sufficient accuracy. By decreasing the threshold value, voxels representing walls of small airways or blood vessels appear in white and at the same time, the large black regions become smaller. On the other hand, when the threshold value is increased, the walls of small airways and blood vessels may not appear in white and the large black regions become larger. Thus, the threshold values may be adjusted to a value so that walls of airways and blood vessels, having a certain size, appear in in white and so that the large black regions grow as large as the pleura of the lungs.
The binarized, segmented and inversed set of 2D slice images may be viewed in the axial, coronal, and sagittal directions.
Filtering out connected components starts from detecting connected components in each 2D slice image. Connected components are displayed as white regions in the large black regions, which represent the lung. When the connected components are detected, an area of each connected component is calculated. In a case where a calculated area is less than a threshold value, the corresponding connected component is filtered out, meaning that pixels in the corresponding connected component are reassigned to zero. In other words, the corresponding connected component is merged into the lung area or the large black regions. After all connected components in the 2D slice image have been filtered out, a filtered out 2D slice image 405 along the axial direction is obtained. The filtered out 2D slice image 405 includes large black regions 420, whose boundaries define the pleural boundary of the lungs.
In step 510, the CT scan image data is processed to generate a set of 2D slice images of the patient's chest. The set of 2D slice images may include the lung, the pleura, the trachea, the heart, the stomach, and any other organs in the chest. A voxel is defined as the intersection of three 2D slice images viewed in the three directions, e.g., the axial, sagittal, and coronal directions. Each voxel has a value representing three pixels, each of which is from a corresponding 2D slice image among the three 2D slice images viewed in the three directions. Segmentation of the set of 2D slice images of the chest starts by determining a seed voxel 220 in the set of 2D slice images, and distinguishes a region including the starting seed voxel from an area not including the starting seed voxel. In step 515, the starting seed voxel is determined automatically or manually. The starting seed voxel may be located at an inside area of the trachea or of the lung. When the trachea is used to determine the starting seed voxel, methods for automatically finding a trachea can be found in commonly assigned U.S. Provisional Patent Application No. 62/020,253 entitled “Trachea Marking”, filed on Jul. 2, 2014 by Lachmanovich et al.; and U.S. Provisional Patent Application No. 62/020,261 entitled “Automatic Detection Of Human Lung Trachea”, filed on Jul. 2, 2014 by Markov et al., the entire contents of both of which is incorporated by reference. After finding the trachea, any voxel inside of the trachea may be determined as the starting seed voxel.
In step 520, a region growing algorithm is applied to process the set of 2D slice images. The region growing algorithm is used to connect more voxels to the seed voxel to grow the region containing the starting seed voxel. This region may be used to determine the borders of the lung at the end of the segmentation method 500.
By applying the region growing algorithm, each voxel of the set of 2D slice images is binarized, that is, it has its Hounsfield value compared to a threshold and reassigned to two colors, black and white. Black voxels represent the lung region and white voxels represent the non-lung region. After the step 520, a new binarized set of 2D slice images is generated. Details of the step 520 is further described with reference to
In step 525, all 2D slice images are processed to filter out connected components therein. Connected components may represent blood vessels and walls of airways having a Hounsfield value greater than the threshold identified above. Connected components are displayed as white areas within the lung area in 2D slice images. By filtering out or removing the white areas having a size smaller than a predetermined threshold from within the lung area, an area covering the lung may be displayed only in black and the 2D slice images are prepared for determining the pleural boundaries of the lungs. Detailed descriptions for the filtering out process are provided with reference to
After the filtering out process, in step 540, the processed set of 2D slice images may be used to identify the pleural boundaries of the lungs. The segmentation method is ended after step 540. By performing the binarization and connected component analysis on the 2D slice images with a minimum size limitation, the true boundaries, particularly the pleural boundaries are clearly defined for the clinician. Alternatively, instead of using the original CT scan image data, the 2D slice images used for process 500 may be generated from a different dataset. An example of another potential dataset to be used for process 500 is a set of 2D slice images generated based on a 3D model which was generated based on the original CT scan image data. Such a 3D model may present a smoother view of the patient's lungs with any gaps, omissions, and/or imperfections in the CT scan image data filled in or corrected.
When it is determined that the intensity of the voxel is greater than or equal to the threshold value, step 524 follows, and otherwise step 522 follows. In step 522, connectedness between the voxel and the starting seed voxel is determined. Here, connectedness ensures that the voxel is in the same region as the starting seed voxel, or in other words, the voxel is in the lung region. This is done by considering the surrounding voxels and determining whether they satisfy the same threshold criteria. This process is continued back to the starting seed voxel. Stated another way, when there is a pathway from the voxel to the starting seed voxel and the pathway is made up of a series of adjacent voxels satisfying the threshold criteria, then the voxel is considered connected to the starting seed voxel. When it is determined that the voxel is connected to the starting seed voxel, the intensity value of the voxel is assigned to the maximum value in step 523. The result is a binary mask of voxel intensities, with all connected voxels satisfying the threshold criteria having one assigned intensity value and all others having another intensity value.
In an aspect, voxels, which are considered as seed voxel candidates, directly adjacent to a seed voxel may be compared with the threshold value recursively without considering connectedness. This recursion method starts by comparing voxels adjacent to the starting seed voxel and recursively proceeds in all directions. In this case, checking for connectedness is unnecessary.
When it is determined that the intensity of the voxel is not less than the threshold value in step 521, or that the voxel is not connected to the starting seed point in step 522, the intensity of the voxel is assigned to the minimum value in step 523. In steps 523 and 524, the set of 2D slice images are segmented, resulting in voxels having only the minimum value or the maximum value.
Step 525 follows after the steps 523 and 524. In step 525, it is determined whether there remain voxels which have not been processed through the region growing algorithm. When there is such a voxel, steps 521-525 are repeated until there are no more unprocessed voxels. The result of steps 521-525 is a binarized set of 2D slice images which only includes black voxels for non-lung regions and white voxels for lung regions. In step 526, the assigned value of each voxel is inversed from the maximum to the minimum or from the minimum to the maximum so that the lung regions have black voxels and the non-lung regions have white voxels. In an aspect, the threshold value used in step 521 may be adjusted to cause sufficient leakage (e.g., a high intensity region) in the lung region such that the high intensity region may cover the lung and the borders of the lung. In another aspect, the threshold value is greater than or equal to the intensity value of the seed voxel.
When it is determined that a connected component is detected in step 527, an area of the detected connected component is calculated in step 528. The area may be a number of pixels in a white area representing the connected component.
In step 529, the calculated area is compared with a threshold value. When it is determined that the calculated area is greater than or equal to the threshold value, it goes back to step 526. That means the connected components may not be a blood vessel or a wall of airways. When it is determined that the calculated area is less than the threshold value, in step 530, the connected component is filtered out or removed. In other words, pixels in the connected components are assigned to the minimum value. In result, the connected component is merged into the large black regions covering the lung and the pleura.
In step 531, it is determined whether there are more 2D slice images to be processed. When it is determined that there is an unprocessed 2D image, steps 526-530 are repeated to filter out connected components in the unprocessed 2D slice image. Otherwise, the filtering out step 525 is complete.
In an aspect, the filtering out step 540 may process 2D slice images taken along one direction and process 2D slice images taken along another direction. For example, step 540 processes 2D slice images taken along the axial direction first and other 2D slice images taken along the coronal direction and the sagittal direction in order. In another aspect, step 540 may process 2D slice images one by one without any order.
By binarizing the set of 2D slice images and filtering out connected components in the 2D slice images, the resulting set of 2D slice images, which is the segmented set of 2D slice images, may be used to determine the pleural boundaries of the lungs, which may be used during either the planning or navigation phases of an ENB procedure. For example, this information may be useful during a planning phase to identify when a suspected target is located outside of the lungs. Alternatively, during a procedure, this information can be used to notify the clinician when a tool is approaching the pleural boundaries. In both instances awareness of the location of the pleural boundaries assists the clinician in avoiding complications such a pneumothorax and other conditions which may occur as the clinician considers or approaches targets at or near the boundaries of the lungs.
Returning now to
Display 130 may be touch-sensitive and/or voice-activated, enabling display 130 to serve as both an input device and an output device. Graphics processors 140 may be specialized graphics processors which perform image-processing functions, such as processing the CT scan image data to generate the 3D model, and process the 3D model to generate the 2D slice images of the 3D model in the various orientations as described above, as well as the 3D renderings of the 3D model. Graphics processors 140 may further be configured to generate a graphical user interface (GUI) to be displayed on display 130. The GUI may include views showing the 2D image slices, the 3D rendering, among other things. In embodiments, graphics processors 140 may be specialized graphics processors, such as a dedicated graphics processing unit (GPU), which performs only the image processing functions so that the one or more general processors 110 may be available for other functions. The specialized GPU may be a stand-alone dedicated graphics card, or an integrated graphics card.
Network interface 150 enables device 100 to communicate with other devices through a wired and/or wireless network connection. In an embodiment, device 100 may receive the CT scan image data from an imaging device via a network connection. In other embodiments, device 100 may receive the CT scan image data via a storage device, such as a disk or other external storage media known to those skilled in the art.
Input interface 160 is used for inputting data or control information, such as setting values, text information, and/or controlling device 100. Input interface 160 may include a keyboard, mouse, touch sensor, camera, microphone, or other data input devices or sensors used for user interaction known to those skilled in the art.
In addition, reference is made to following commonly assigned applications: U.S. Provisional Patent Application Ser. No. 62/020,220 entitled “Real-Time Automatic Registration Feedback”, filed on Jul. 2, 2014, by Brown et al.; U.S. Provisional Patent Application Ser. No. 62/020,242 entitled “Unified Coordinate System for Multiple CT Scans of Patient Lungs”, filed on Jul. 2, 2014, by Greenburg.; U.S. Provisional Patent Application Ser. No. 62/020,245 entitled “Alignment CT”, filed on Jul. 2, 2014, by Klein et al.; U.S. Provisional Patent Application Ser. No. 62/020,250 entitled “Algorithm for Fluoroscopic Pose Estimation”, filed on Jul. 2, 2014, by Merlet.; U.S. Provisional Patent Application Ser. No. 62/020,258 entitled “Cone View—A Method of Providing Distance and Orientation Feedback While Navigating in 3D”, filed on Jul. 2, 2014, by Lachmanovich et al.; and U.S. Provisional Patent Application Ser. No. 62/020,262 entitled “Dynamic 3D Lung Map View for Tool Navigation Inside the Lung”, filed on Jul. 2, 2014, by Weingarten et al., the entire contents of all of which are hereby incorporated by reference. All of these references are directed to aspects of modifying and manipulating the DICOM images to provide enhanced clarity and performance for analysis, diagnostic, and treatment systems relating to, among other things, lung treatment planning and navigation.
Although embodiments have been described in detail with reference to the accompanying drawings for the purpose of illustration and description, it is to be understood that the inventive processes and apparatus are not to be construed as limited thereby. It will be apparent to those of ordinary skill in the art that various modifications to the foregoing embodiments may be made without departing from the scope of the disclosure.
The present application is a continuation of U.S. patent application Ser. No. 15/808,975, filed on Nov. 10, 2017, which is a continuation of U.S. patent application Ser. No. 14/754,867, now U.S. Patent No. 9,836,848, filed on Jun. 30, 2015, and entitled “SYSTEM AND METHOD FOR SEGMENTATION OF LUNG,” which claims the benefit of and priority to U.S. Provisional Patent Application Serial No. 62/020,261 filed on Jul. 2, 2014, the entire contents of each of which are incorporated herein by reference.
Number | Name | Date | Kind |
---|---|---|---|
5592939 | Martinelli | Jan 1997 | A |
5611025 | Lorensen et al. | Mar 1997 | A |
5676673 | Ferre et al. | Oct 1997 | A |
5697377 | Wittkampf | Dec 1997 | A |
5699799 | Xu et al. | Dec 1997 | A |
5715836 | Kliegis et al. | Feb 1998 | A |
5729129 | Acker | Mar 1998 | A |
5752513 | Acker et al. | May 1998 | A |
5782762 | Vining | Jul 1998 | A |
5881124 | Giger et al. | Mar 1999 | A |
5891030 | Johnson et al. | Apr 1999 | A |
5913820 | Bladen et al. | Jun 1999 | A |
5920319 | Vining et al. | Jul 1999 | A |
5967980 | Ferre et al. | Oct 1999 | A |
5971767 | Kaufman et al. | Oct 1999 | A |
5987960 | Messner et al. | Nov 1999 | A |
6019725 | Vesely et al. | Feb 2000 | A |
6047080 | Chen et al. | Apr 2000 | A |
6083162 | Vining | Jul 2000 | A |
6138045 | Kupinski et al. | Oct 2000 | A |
6151404 | Pieper | Nov 2000 | A |
6167296 | Shahidi | Dec 2000 | A |
6181348 | Geiger | Jan 2001 | B1 |
6201387 | Govari | Mar 2001 | B1 |
6233476 | Strommer et al. | May 2001 | B1 |
6246784 | Summers et al. | Jun 2001 | B1 |
6266551 | Osadchy et al. | Jul 2001 | B1 |
6332089 | Acker et al. | Dec 2001 | B1 |
6346940 | Fukunaga | Feb 2002 | B1 |
6366800 | Vining et al. | Apr 2002 | B1 |
6381485 | Hunter et al. | Apr 2002 | B1 |
6387092 | Burnside et al. | May 2002 | B1 |
6466815 | Saito | Oct 2002 | B1 |
6496188 | Deschamps et al. | Dec 2002 | B1 |
6501848 | Carroll et al. | Dec 2002 | B1 |
6501981 | Schweikard et al. | Dec 2002 | B1 |
6505065 | Yanof et al. | Jan 2003 | B1 |
6522907 | Bladen et al. | Feb 2003 | B1 |
6526162 | Asano et al. | Feb 2003 | B2 |
6535756 | Simon et al. | Mar 2003 | B1 |
6578579 | Burnside et al. | Jun 2003 | B2 |
6584174 | Schubert et al. | Jun 2003 | B2 |
6603868 | Ludwig et al. | Aug 2003 | B1 |
6611793 | Burnside et al. | Aug 2003 | B1 |
6650927 | Keidar | Nov 2003 | B1 |
6651669 | Burnside | Nov 2003 | B1 |
6694163 | Vining | Feb 2004 | B1 |
6757557 | Bladen et al. | Jun 2004 | B1 |
6783523 | Qin et al. | Aug 2004 | B2 |
6792390 | Burnside et al. | Sep 2004 | B1 |
6829379 | Knoplioch et al. | Dec 2004 | B1 |
6850794 | Shahidi | Feb 2005 | B2 |
6892090 | Verard et al. | May 2005 | B2 |
6898263 | Avinash et al. | May 2005 | B2 |
6909913 | Vining | Jun 2005 | B2 |
6920347 | Simon et al. | Jul 2005 | B2 |
6925200 | Wood et al. | Aug 2005 | B2 |
7006677 | Manjeshwar et al. | Feb 2006 | B2 |
7072501 | Wood et al. | Jul 2006 | B2 |
7085400 | Holsing et al. | Aug 2006 | B1 |
7096148 | Anderson et al. | Aug 2006 | B2 |
7149564 | Vining et al. | Dec 2006 | B2 |
7167180 | Shibolet | Jan 2007 | B1 |
7174202 | Bladen et al. | Feb 2007 | B2 |
7179220 | Kukuk | Feb 2007 | B2 |
7236558 | Saito et al. | Jun 2007 | B2 |
7301332 | Govari et al. | Nov 2007 | B2 |
7315639 | Kuhnigk | Jan 2008 | B2 |
7324104 | Bitter et al. | Jan 2008 | B1 |
7336809 | Zeng et al. | Feb 2008 | B2 |
7397937 | Schneider et al. | Jul 2008 | B2 |
7428334 | Schoisswohl et al. | Sep 2008 | B2 |
7452357 | Vlegele et al. | Nov 2008 | B2 |
7505809 | Strommer et al. | Mar 2009 | B2 |
7517320 | Wibowo et al. | Apr 2009 | B2 |
7518619 | Stoval, III et al. | Apr 2009 | B2 |
7630752 | Viswanathan | Dec 2009 | B2 |
7630753 | Simon | Dec 2009 | B2 |
7659912 | Akimoto et al. | Feb 2010 | B2 |
7702153 | Hong et al. | Apr 2010 | B2 |
7751865 | Jascob et al. | Jul 2010 | B2 |
7756316 | Odry et al. | Jul 2010 | B2 |
7788060 | Schneider | Aug 2010 | B2 |
7792565 | Vining | Sep 2010 | B2 |
7805269 | Glossop | Sep 2010 | B2 |
7809176 | Gundel | Oct 2010 | B2 |
7811294 | Strommer et al. | Oct 2010 | B2 |
7822461 | Geiger et al. | Oct 2010 | B2 |
7901348 | Soper et al. | Mar 2011 | B2 |
7907772 | Wang et al. | Mar 2011 | B2 |
7929014 | Akimoto et al. | Apr 2011 | B2 |
7951070 | Ozaki et al. | May 2011 | B2 |
7969142 | Krueger et al. | Jun 2011 | B2 |
7985187 | Wibowo et al. | Jul 2011 | B2 |
8009891 | de Vaan | Aug 2011 | B2 |
8049777 | Akimoto et al. | Nov 2011 | B2 |
8055323 | Sawyer | Nov 2011 | B2 |
8102416 | Ito et al. | Jan 2012 | B2 |
8126241 | Zarkh et al. | Feb 2012 | B2 |
8131344 | Strommer et al. | Mar 2012 | B2 |
8170328 | Masumoto et al. | May 2012 | B2 |
8199981 | Koptenko et al. | Jun 2012 | B2 |
8200314 | Bladen et al. | Jun 2012 | B2 |
8202213 | Ito et al. | Jun 2012 | B2 |
8208708 | Homan et al. | Jun 2012 | B2 |
8219179 | Ganatra et al. | Jul 2012 | B2 |
8257346 | Qin et al. | Sep 2012 | B2 |
8267927 | Dalal et al. | Sep 2012 | B2 |
8290228 | Cohen et al. | Oct 2012 | B2 |
8298135 | Ito et al. | Oct 2012 | B2 |
8385688 | Gong et al. | Feb 2013 | B2 |
8391952 | Anderson | Mar 2013 | B2 |
8417009 | Mizuno | Apr 2013 | B2 |
8494612 | Vetter et al. | Jul 2013 | B2 |
8509877 | Mori et al. | Aug 2013 | B2 |
8672836 | Higgins et al. | Mar 2014 | B2 |
8682045 | Vining et al. | Mar 2014 | B2 |
8696549 | Holsing et al. | Apr 2014 | B2 |
8698806 | Kunert et al. | Apr 2014 | B2 |
8700132 | Ganatra et al. | Apr 2014 | B2 |
8706193 | Govari et al. | Apr 2014 | B2 |
8709034 | Keast et al. | Apr 2014 | B2 |
8730237 | Ruijters et al. | May 2014 | B2 |
8768029 | Helm et al. | Jul 2014 | B2 |
8784400 | Roschak | Jul 2014 | B2 |
8798227 | Tsukagoshi et al. | Aug 2014 | B2 |
8798339 | Mielekamp et al. | Aug 2014 | B2 |
8801601 | Prisco et al. | Aug 2014 | B2 |
8819591 | Wang et al. | Aug 2014 | B2 |
8862204 | Sobe et al. | Oct 2014 | B2 |
9836848 | Markov et al. | Dec 2017 | B2 |
10074185 | Markov et al. | Sep 2018 | B2 |
20020009215 | Armato, III | Jan 2002 | A1 |
20030095692 | Mundy et al. | May 2003 | A1 |
20030099385 | Zeng | May 2003 | A1 |
20050001832 | Shen | Jan 2005 | A1 |
20050286750 | Dehmeshki | Dec 2005 | A1 |
20070127802 | Odry | Jun 2007 | A1 |
20070249931 | Fain | Oct 2007 | A1 |
20080183073 | Higgins et al. | Jul 2008 | A1 |
20090012390 | Pescatore et al. | Jan 2009 | A1 |
20090030306 | Miyoshi et al. | Jan 2009 | A1 |
20100128940 | Buelow et al. | May 2010 | A1 |
20100235352 | Slutsky et al. | Sep 2010 | A1 |
20100310146 | Higgins et al. | Dec 2010 | A1 |
20100312094 | Guttman et al. | Dec 2010 | A1 |
20110237897 | Gilboa | Sep 2011 | A1 |
20110251607 | Kruecker et al. | Oct 2011 | A1 |
20120203065 | Higgins et al. | Aug 2012 | A1 |
20120249546 | Tschirren et al. | Oct 2012 | A1 |
20120280135 | Bal | Nov 2012 | A1 |
20120287238 | Onishi et al. | Nov 2012 | A1 |
20130165854 | Sandhu et al. | Jun 2013 | A1 |
Number | Date | Country |
---|---|---|
2007-0085120 | Aug 2007 | KR |
2011-0018573 | Feb 2011 | KR |
2008017984 | Feb 2008 | WO |
Entry |
---|
Australian Examination Report issued in corresponding Appl. No. AU 2015284218 dated Mar. 7, 2019 (3 pages). |
Extended European Search report for application No. 15814448.5 dated Dec. 6, 2017 (7 pages). |
Brown et al., Method for Segmenting Chest CT Image Data Using an Anatomical Model: Preliminary Rreults, IEEE Transactiosn on Medical Imaging, vol. 16, No. 6, Dec. 1, 1997, pp. 828-839. |
Japanese Office Action issued in Appl. No. JP 2016-575072 dated Nov. 5, 2018, together with English language translation (8 pages). |
PCT Search Report for International Patent Application No. PCT/US2015/038541, dated Oct. 13, 2015 (4 pages). |
Number | Date | Country | |
---|---|---|---|
20190019294 A1 | Jan 2019 | US |
Number | Date | Country | |
---|---|---|---|
62020261 | Jul 2014 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 15808975 | Nov 2017 | US |
Child | 16124371 | US | |
Parent | 14754867 | Jun 2015 | US |
Child | 15808975 | US |