CN111415356B - 肺炎征象的分割方法、装置、介质及电子设备 - Google Patents
肺炎征象的分割方法、装置、介质及电子设备 Download PDFInfo
- Publication number
- CN111415356B CN111415356B CN202010189089.6A CN202010189089A CN111415356B CN 111415356 B CN111415356 B CN 111415356B CN 202010189089 A CN202010189089 A CN 202010189089A CN 111415356 B CN111415356 B CN 111415356B
- Authority
- CN
- China
- Prior art keywords
- image
- pneumonia
- images
- region
- symptom
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 206010035664 Pneumonia Diseases 0.000 title claims abstract description 151
- 230000011218 segmentation Effects 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims abstract description 58
- 208000024891 symptom Diseases 0.000 title claims description 83
- 210000004072 lung Anatomy 0.000 claims description 143
- 238000003062 neural network model Methods 0.000 claims description 64
- 238000004590 computer program Methods 0.000 claims description 9
- 230000003628 erosive effect Effects 0.000 claims description 8
- 206010025080 Lung consolidation Diseases 0.000 claims description 5
- 230000007797 corrosion Effects 0.000 claims description 4
- 238000005260 corrosion Methods 0.000 claims description 4
- 239000005337 ground glass Substances 0.000 claims description 4
- 238000002591 computed tomography Methods 0.000 description 112
- 210000000988 bone and bone Anatomy 0.000 description 30
- 238000012545 processing Methods 0.000 description 12
- 210000001519 tissue Anatomy 0.000 description 10
- 238000001514 detection method Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 6
- 208000004434 Calcinosis Diseases 0.000 description 5
- 230000002308 calcification Effects 0.000 description 5
- 238000007781 pre-processing Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000012549 training Methods 0.000 description 5
- 241000894006 Bacteria Species 0.000 description 4
- 241000711573 Coronaviridae Species 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 4
- 230000002829 reductive effect Effects 0.000 description 4
- 201000001178 Bacterial Pneumonia Diseases 0.000 description 3
- 206010035737 Pneumonia viral Diseases 0.000 description 3
- 241000700605 Viruses Species 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 208000009421 viral pneumonia Diseases 0.000 description 3
- 206010011224 Cough Diseases 0.000 description 2
- 241000270295 Serpentes Species 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 210000000621 bronchi Anatomy 0.000 description 2
- 210000003123 bronchiole Anatomy 0.000 description 2
- 238000013135 deep learning Methods 0.000 description 2
- 230000010339 dilation Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000036961 partial effect Effects 0.000 description 2
- 231100000915 pathological change Toxicity 0.000 description 2
- 230000036285 pathological change Effects 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 241001678559 COVID-19 virus Species 0.000 description 1
- 208000035473 Communicable disease Diseases 0.000 description 1
- 206010061218 Inflammation Diseases 0.000 description 1
- 206010030113 Oedema Diseases 0.000 description 1
- 241000125945 Protoparvovirus Species 0.000 description 1
- 241000194017 Streptococcus Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000740 bleeding effect Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 239000005338 frosted glass Substances 0.000 description 1
- 125000001475 halogen functional group Chemical group 0.000 description 1
- 230000002949 hemolytic effect Effects 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 208000015181 infectious disease Diseases 0.000 description 1
- 230000004054 inflammatory process Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000017074 necrotic cell death Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 208000011580 syndromic disease Diseases 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
- 241000712461 unidentified influenza virus Species 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
- G06T5/30—Erosion or dilatation, e.g. thinning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/149—Segmentation; Edge detection involving deformable models, e.g. active contour models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20084—Artificial neural networks [ANN]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20112—Image segmentation details
- G06T2207/20132—Image cropping
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30061—Lung
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Software Systems (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Molecular Biology (AREA)
- Computational Linguistics (AREA)
- General Engineering & Computer Science (AREA)
- Biophysics (AREA)
- Mathematical Physics (AREA)
- Biomedical Technology (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Medical Informatics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种肺炎征象的分割方法、分割装置、计算机可读存储介质以及电子设备,通过将CT图像中的肺部区域图像分别输入多个神经网络模型,以此分别多个肺炎征象图像,然后将多个肺炎征象图像进行组合,得到肺炎综合征象图像,利用神经网络模型可以高效且准确的获取大量的待检测CT图像中的肺炎征象图像,从而为后续诊断肺炎提供可靠的数据依据。
Description
技术领域
本申请涉及图像处理领域,具体涉及基于CT图像中肺炎征象的分割方法、分割装置、计算机可读存储介质以及电子设备。
背景技术
电脑断层扫描(Computed Tomography,简称CT)是一种利用数位几何处理后重建的三维放射线医学影像。该技术主要通过单一轴面的X射线旋转照射人体,由于不同的组织对X射线的吸收能力(或称阻射率)不同,可以用电脑的三维技术重建出断层面影像,经由窗宽、窗位处理,可以得到相应组织的断层影像,将断层影像层层堆叠,即可形成立体影像。
通过CT图像可以获知被检测者是否为肺炎患者,特别是应对新型冠状病毒所造成的肺炎患者的检测,CT图像检测是最为重要且最为准确的方法之一。目前在获取CT图像后大多是先分割出肺部区域,然后通过人工去分割肺炎病灶或征象,并根据肺炎病灶或征象判断是否为肺炎以及肺炎的程度,这样的分割方式明显效率太低,特别是针对目前新型冠状病毒所造成的病毒性肺炎,因为其传染性很强,需要对大量疑似人员快速的甄别,因此,在保证准确度的同时也要提高效率,以尽早遏制住病毒的传播。因此,目前亟需高精度且高效率的肺部病灶或征象的分割方法。
发明内容
为了解决上述技术问题,本申请提出了一种肺炎征象的分割方法、分割装置、计算机可读存储介质以及电子设备,通过将CT图像中的肺部区域图像分别输入多个神经网络模型,以此分别多个肺炎征象图像,然后将多个肺炎征象图像进行组合,得到肺炎综合征象图像,利用神经网络模型可以高效且准确的获取大量的待检测CT图像中的肺炎征象图像,从而为后续诊断肺炎提供可靠的数据依据。
根据本申请的一个方面,提供了一种肺炎征象的分割方法,包括:基于CT图像中的肺部区域图像,分别生成多个肺炎征象图像;以及将所述多个肺炎征象图像组合,得到肺炎综合征象图像;其中,分别生成多个肺炎征象图像的方式包括:将所述肺部区域图像分别输入多个神经网络模型,得到所述多个肺炎征象图像。
在一实施例中,所述多个肺炎征象图像包括如下征象图像中的任一种或多种的组合:肺实变图像、磨玻璃影图像、肿块图像、树芽征图像、结节图像、空洞图像、晕征图像。
在一实施例中,所述肺部区域图像包括多层二维图像;所述将所述肺部区域图像分别输入多个神经网络模型,得到所述多个肺炎征象图像包括:逐次将所述多层二维图像中的一部分分别输入所述多个神经网络模型,得到分别对应所述多个肺炎征象图像的多层二维征象图像;以及分别将对应同一肺炎征象图像的所述多层二维征象图像叠加,得到所述多个肺炎征象图像。
在一实施例中,所述肺部区域图像包括多层二维图像;所述将所述肺部区域图像分别输入多个神经网络模型,得到所述多个肺炎征象图像包括:将所述多层二维图像分别输入所述多个神经网络模型,得到分别对应所述多个肺炎征象图像的多层二维征象图像;以及分别将对应同一肺炎征象图像的所述多层二维征象图像叠加,得到所述多个肺炎征象图像。
在一实施例中,在所述得到所述多个肺炎征象图像之后,所述方法还包括:对所述多个肺炎征象图像分别进行腐蚀膨胀操作。
在一实施例中,所述肺部区域图像的获取方式包括:获取所述CT图像中的肋骨区域图像;获取所述CT图像中肺部区域的粗分割图像;以及以所述肋骨区域图像为所述肺部区域的边界,并且以所述粗分割图像为种子区域,以预设的步长向周围扩张至所述肺部区域的边界,得到所述肺部区域图像。
在一实施例中,在所述获取所述CT图像中肺部区域的粗分割图像之后,所述方法还包括:对所述粗分割图像进行腐蚀操作,得到腐蚀后的粗分割图像;所述以所述肋骨区域图像为所述肺部区域的边界,并且以所述粗分割图像为种子区域,以预设的步长向周围扩张至所述肺部区域的边界,得到所述肺部区域图像包括:以所述肋骨区域图像为所述肺部区域的边界,并且以所述腐蚀后的粗分割图像为种子区域扩张,以预设的步长向周围扩张至所述肺部区域的边界,得到所述肺部区域图像。
根据本申请的另一个方面,提供了一种肺炎征象的分割装置,包括:生成模块,用于基于CT图像中的肺部区域图像,分别生成多个肺炎征象图像;以及组合模块,用于将所述多个肺炎征象图像组合,得到肺炎综合征象图像;其中,生成模块进一步配置为:将所述肺部区域图像分别输入多个神经网络模型,得到所述多个肺炎征象图像。
根据本申请的另一个方面,提供了一种计算机可读存储介质,所述存储介质存储有计算机程序,所述计算机程序用于执行上述任一所述的分割方法。
根据本申请的另一个方面,提供了一种电子设备,所述电子设备包括:处理器;用于存储所述处理器可执行指令的存储器;所述处理器,用于执行上述任一所述的分割方法。
本申请提供的一种肺炎征象的分割方法、分割装置、计算机可读存储介质以及电子设备,通过将CT图像中的肺部区域图像分别输入多个神经网络模型,以此分别多个肺炎征象图像,然后将多个肺炎征象图像进行组合,得到肺炎综合征象图像,利用神经网络模型可以高效且准确的获取大量的待检测CT图像中的肺炎征象图像,从而为后续诊断肺炎提供可靠的数据依据。
附图说明
通过结合附图对本申请实施例进行更详细的描述,本申请的上述以及其他目的、特征和优势将变得更加明显。附图用来提供对本申请实施例的进一步理解,并且构成说明书的一部分,与本申请实施例一起用于解释本申请,并不构成对本申请的限制。在附图中,相同的参考标号通常代表相同部件或步骤。
图1是本申请一示例性实施例提供的一种肺炎征象的分割方法的流程示意图。
图2是本申请另一示例性实施例提供的一种肺炎征象的分割方法的流程示意图。
图3是本申请另一示例性实施例提供的一种肺炎征象的分割方法的流程示意图。
图4是本申请一示例性实施例提供的一种肺部区域图像分割方法的流程示意图。
图5是本申请另一示例性实施例提供的一种肺部区域分割方法的流程示意图。
图6是本申请一示例性实施例提供的一种肋骨区域图像的获取方法的流程示意图。
图7是本申请另一示例性实施例提供的一种肺部区域分割方法的流程示意图。
图8是本申请另一示例性实施例提供的一种肺部分区域割方法的流程示意图。
图9是本申请一示例性实施例提供的一种肺炎征象的分割装置的结构示意图。
图10是本申请一示例性实施例提供的电子设备的结构图。
具体实施方式
下面,将参考附图详细地描述根据本申请的示例实施例。显然,所描述的实施例仅仅是本申请的一部分实施例,而不是本申请的全部实施例,应理解,本申请不受这里描述的示例实施例的限制。
申请概述
肺炎是一种因为多种因素引发的肺部的炎症,对于肺炎患者来说在疾病发生的时候往往会伴随剧烈的咳嗽现象。肺炎分为两种,一种是细菌引起的细菌性的肺炎,因为细菌侵袭肺部导致的肺炎,其中最为常见的肺炎球菌,和甲型溶血性链球菌这两种细菌细菌,大部分的细菌性肺炎都是因为这两种细菌导致的;另一种就是病毒引起的病毒性的肺炎,最常见的就是由一级巨细病毒引发的流感病毒,又例如2019新型冠状病毒等,病毒性的肺炎相对来于细菌性的肺炎来说,不仅更严重一些,治疗的难度相对来说也会更高。
如上所述,CT图像检测是获知被检测者是否为肺炎患者最为重要且最为准确的方法之一,特别是应对新型冠状病毒所造成的肺炎患者。通常得到CT图像后需要对其中的肺部区域进行分割,以得到肺部区域图像,然后对肺部区域内的肺炎征象进行分割,最后由医生根据分割得到的肺炎征象来判断被检测者是否患有肺炎,现有的分割肺炎征象的工作通常是由专业的医护人员手动实现,显然这样的效率是较低的,特别是应对传染性很强的新型冠状病毒所造成的肺炎的检测,因为其具有高传染性而导致大量的人有感染隐患,这就需要对大量的人做CT图像检测,而相对有限甚至紧缺的医务人员显然不能满足大量的CT图像检测的需求。
随着图像处理的快速发展,越来越多的医学图像的处理都可以通过计算机进行处理,例如图像分割,以得到诊断所需要的感兴趣区域图像或者基础数据等。然而肺炎征象包括多个,其中有的征象在CT图像中与其他组织的密度(在CT图像中表象为CT值)区别较大,可以较为容易的分割出来,而有的征象在CT图像中与其他组织的密度区别较小,很难通过密度的比较进行分割,或者分割的精度不高;同时,由于肺炎患者常常因为咳嗽和吸气不足而导致图像的运动噪声较大,并且肺部内密度升高,从而导致肺部区域分割时存在较大的难度,同时,肺炎病灶的区域较大会影响到肺部的形状和结构,从而导致肺部的分割难度进一步提高,这样会影响到分割肺炎征象的精度,产生一定的误差,并且该误差的传递和累积会最终影响医生的判断,造成不可估量的后果。
出于解决上述问题,本申请提供的一种肺炎征象的分割方法、分割装置、计算机可读存储介质以及电子设备,通过将CT图像中的肺部区域图像分别输入多个神经网络模型,以此分别多个肺炎征象图像,然后将多个肺炎征象图像进行组合,得到肺炎综合征象图像,利用神经网络模型可以高效且准确的获取大量的待检测CT图像中的肺炎征象图像,从而为后续诊断肺炎提供可靠的数据依据。
示例性方法
图1是本申请一示例性实施例提供的一种肺炎征象的分割方法的流程示意图。如图1所示,该方法包括如下步骤:
步骤110:基于CT图像中的肺部区域图像,分别生成多个肺炎征象图像;其中,分别生成多个肺炎征象图像的方式包括:将肺部区域图像分别输入多个神经网络模型,得到多个肺炎征象图像。
通过拍摄CT图像来获知被检测者肺部区域图像,并且根据该肺部区域图像可以准确的判断该被检测者是否患有肺炎,对于目前严峻的新型冠状病毒性肺炎,要遏制住病毒的蔓延,就需要将患有肺炎或携带该病毒的患者与其他人群隔离开来,如此就需要对疑似人员进行甄别,而作为甄别肺炎患者的有效措施之一的CT图像检测是特别重要的。然而,CT图像的检测是较为复杂的,其效率不高,并且由于肺炎征象包括多个,需要在获知所有的或者多个征象后才能综合判断是否患有肺炎以及患病的程度,而这多个征象中有些征象的获取难度较大,这也进一步增加了CT图像检测的难度。因此,本申请实施例中通过将肺部区域图像分别输入多个神经网络模型,针对不同的肺炎征象训练得到不同的神经网络模型,可以有针对性的得到精度较高的各个肺炎征象图像,从而为后续肺炎的甄别提供了准确的数据依据。
在一实施例中,肺炎征象图像可以包括如下征象图像中的任一种或多种的组合:肺实变图像、磨玻璃影图像、肿块图像、树芽征图像、结节图像、空洞图像、晕征图像。肺实变是指终末细支气管远侧的含气腔隙被病理性液体、细胞、组织所代替,其最主要的特点是病变区致密,血管不能显影。磨玻璃影是指各种原因引起的肺泡充填或间质增厚,导致肺密度轻度增高,但其内仍可见血管纹理的病变阴影。肿块是指团块状密度增高影,其最大直径≧3cm。结节是指结节状密度增高影,其最大直径<3cm。空洞是指肺内病变坏死、液化,经引流支气管排除及气体进入而形成的透亮区,对于感染性疾病,空洞临床上会倾向于化脓性感染。树芽征是指由终末细支气管和肺泡腔内病变形成的小结节影与分支细线影构成的酷似春天的树枝发芽状,多在肺外围支气管末梢呈2-4mm大小结节与树枝状的高密度影。晕征是指结节/空洞周围环绕的类环形玻璃样密度影,通常代表渗出、出血或水肿。可以根据各个征象的不同特征设置不同的神经网络模型来分别分割,以提高整体的分割精度。
在一实施例中,神经网络模型可以是深度学习神经网络模型,优选地,该神经网络模型可以是Unet神经网络模型。在一实施例中,神经网络模型在训练方式可以是:选取已经由专业医务人员分割并标注出肺炎征象的肺部区域图像作为神经网络模型的训练样本来训练该神经网络模型;并且神经网络模型在分割过程中得到的分割结果也可以由第三方检测机构进行验证和修改,并且可以将修改后的结果作为样本再次训练神经网络模型,从而进一步提高神经网络模型的分割精度。应当理解,本申请实施例中的多个神经网络模型可以是一个神经网络模型,也可以是多个神经网络模型,且多个神经网络模型的类型可以相同,也可以不同,本申请实施例对于分割各个肺炎征象的神经网络模型的具体类型不做限定。
步骤120:将多个肺炎征象图像组合,得到肺炎综合征象图像。
通常肺炎的确定是基于多个肺炎征象综合判断的,例如当只存在结节图像而不存在其他肺炎征象时,是不能确定被检测者为肺炎患者的,因此,在分别得到多个肺炎征象图像(通常是在CT图像中标注出各个肺炎征象的区域)后,将多个肺炎征象图像组合到一起,得到肺炎综合征象图像,以方便医务人员或者其他的检测机构根据该肺炎综合征象图像准确判断被检测者是否为肺炎患者。
本申请提出了一种肺炎征象的分割方法,通过将CT图像中的肺部区域图像分别输入多个神经网络模型,以此分别多个肺炎征象图像,然后将多个肺炎征象图像进行组合,得到肺炎综合征象图像,利用神经网络模型可以高效且准确的获取大量的待检测CT图像中的肺炎征象图像,从而为后续诊断肺炎提供可靠的数据依据。
图2是本申请另一示例性实施例提供的一种肺炎征象的分割方法的流程示意图。肺部区域图像包括多层二维图像,如图2所示,步骤110可以具体包括如下子步骤:
步骤111:逐次将多层二维图像中的一部分输入多个神经网络模型,得到分别对应多个肺炎征象图像的多层二维征象图像。
二维图像的分割要比三维图像的分割难度小、速度快,而CT图像是包括多层二维图像的,CT图像中的肺部区域图像也是包括多层二维图像的,因此,为了提高分割的效率,可以将肺部区域图像的多层二维图像分为多个部分,将这多个部分以多次输入的方式输入多个神经网络模型,由多个神经网络模型分别多次分割二维图像来得到对应的多层二维征象图像,从而提高分割效率。应当理解,本申请实施例可以根据神经网络或处理机器的处理能力适当选取单次输入神经网络模型的二维图像的层数,可以是一层,也可以是多层,还可以是将所有的二维图像单次输入神经网络模型,只要所选取的层数不超过神经网络模型或处理机器所能承载的负荷即可,本申请实施例对于单次输入神经网络模型的二维图像的具体层数不做限定。
步骤112:分别将对应同一肺炎征象图像的多层二维征象图像叠加,得到多个肺炎征象图像。
在得到每个肺炎征象的多层二维征象图像之后,将同一肺炎征象图像的多层二维征象图像叠加起来得到该肺炎征象图像。在一实施例中,相邻次输入神经网络模型的部分二维图像之间存在交叉部分。通过设置交叉部分,可以避免边缘二维图像之间的差异较大,而且可以通过交叉部分的定位更好的实现叠加。
图3是本申请另一示例性实施例提供的一种肺炎征象的分割方法的流程示意图。如图3所示,在步骤112之后,上述实施例还可以包括:
步骤113:对多个肺炎征象图像分别进行腐蚀膨胀操作。
腐蚀操作为形态学操作,其具体操作过程为:沿着图像中物体边界移除像素并缩小物体的大小,即缩小物体的边界以去除图像中物体的噪声。膨胀操作也是形态学操作,其具体操作过程与腐蚀操作正好相反,即沿着图像中物体边界增加像素并扩大物体的大小。通过腐蚀膨胀操作,可以有效去除分割过程中产生的噪声,同时由于相邻的二维征象图像之间是存在相互关联的,通过对叠加后的肺炎征象图像进行腐蚀膨胀操作,可以利用相邻的二维征象图像之间的相互关联性,可以去除个别层的分割误差,例如通过上下两层或多层二维征象图像可以调整中间层的二维征象图像,从而提高肺炎征象图像整体的分割精度。
图4是本申请一示例性实施例提供的一种肺部区域图像分割方法的流程示意图。如图4所示,该方法包括如下步骤:
步骤410:获取CT图像中的肋骨区域图像。
由于肋骨是紧紧包裹肺部区域的,并且肺炎患者的肺部形状和部分区域的特性会影响到CT图像的成像,而肺炎患者的肋骨是不会因为患病而发生变化的,因此,可以通过获取CT图像中的肋骨区域来得到肺部区域的外边界,从而提高分割肺部区域的精度。
步骤420:获取CT图像中肺部区域的粗分割图像。
在一实施例中,步骤420的具体实现方式可以为:将CT图像输入神经网络模型,得到肺部区域的粗分割图像。通训练好的神经网络模型,可以由神经网络模型识别CT图像中的肺部区域,其中神经网络模型可以是Unet网络等神经网络模型,并且神经网络模型在训练方式可以是:选取已经由专业医务人员识别并标注出肺部区域的CT图像作为神经网络模型的训练样本来训练该神经网络模型。由于该步骤中仅仅是获取CT图像中肺部区域的粗分割图像,并非肺部区域的准确分割图像,因此,可以选取适当数量的训练样本来训练该神经网络模型,从而提高整个肺部分割的效率。
在另一实施例中,步骤420的具体实现方式还可以为:根据肺部区域的CT值,选取CT值在肺部区域的CT值范围内的区域作为粗分割图像。CT值是测定人体某一局部组织或器官密度大小的一种计量单位,通常称亨氏单位(hounsfield unit,HU),其中,空气的CT值为-1000,致密骨的CT值为+1000。实际上CT值是CT图像中各组织与X线衰减系数相当的对应值,CT值不是绝对不变的数值,它不仅与人体内在因素如呼吸、血流等有关,而且与X线管电压、CT装置、室内温度等外界因素有关。人体组织中除了骨骼以为,其他组织的CT值都在300以下且在-80以上,其中钙化点的CT值为80-300、脂肪的CT值为-20—-80。由于肺部区域内基本为空气,其CT值相对其他组织较低,因此可以设定CT值范围,选取CT值在该CT值范围内的连通区域作为肺部区域的粗分割图像。
步骤430:以肋骨区域图像为肺部区域的边界,并且以粗分割图像为种子区域,以预设的步长向周围扩张至肺部区域的边界,得到肺部区域图像。
在一实施例中,步骤430的具体实现方式可以是:以肋骨区域图像为边界、粗分割图像为感兴趣区域,并通过活动轮廓模型、以预设的步长对CT图像进行分割,以得到肺部区域图像。由于肺部是被肋骨紧紧包裹的,即肋骨区域为肺部区域的外边界,因此,可以将粗分割图像作为种子区域或者感兴趣区域,并且以肋骨区域图像作为边界,通过活动轮廓模型、以预设的步长对CT图像进行分割,即由感兴趣区域开始、以预设的步长向周围扩张,直至扩张到肋骨区域为止,以得到精确分割的肺部区域的图像,为后续的肺叶分割、肺炎判断等提供了准确的基础图像数据。其中,预设的步长可以根据实际需求而调整,为了得到更高的精度,可以适当减小预设的步长。在进一步的实施例中,活动轮廓模型可以包括LevelSet模型或Snake模型。应当理解,本申请实施例可以根据实际应用场景的需求而选取不同的活动轮廓模型,只要所选取的活动轮廓模型可以以肋骨区域图像为边界、粗分割图像为种子区域得到精确的肺部区域的图像即可,本申请实施例对于活动轮廓模型的具体结构不做限定。
图5是本申请另一示例性实施例提供的一种肺部区域分割方法的流程示意图。如图5所示,在步骤420之后,上述实施例还可以包括:
步骤440:对粗分割图像进行腐蚀操作,得到腐蚀后的粗分割图像。
由于肺部区域的粗分割图像通常不是准确的肺部区域图像,例如,粗分割图像中可能会包含肺部区域以外的其他区域图像,该其他区域图像即为粗分割图像中的噪声区域,通过腐蚀操作可以将粗分割图像中的干扰噪声区域去除,以保证腐蚀后的粗分割图像为肺部区域图像的一部分,而不包含肺部区域以外的区域图像。
同时,步骤430调整为:以肋骨区域图像为肺部区域的边界,并且以腐蚀后的粗分割图像为种子区域,以预设的步长向周围扩张至肺部区域的边界,得到肺部区域图像。
由于种子区域只需要是肺部区域的一部分即可,不需要是肺部区域的全部,因此,可以通过腐蚀操作来保证种子区域内只包含肺部区域,从而避免种子区域在扩张的过程中融入更多的非肺部区域,并且种子区域在扩张过程中是不会去除已存在于该种子区域内的区域的,因此,保证种子区域内只包含肺部区域是保证分割肺部区域的精度的前提条件。
图6是本申请一示例性实施例提供的一种肋骨区域图像的获取方法的流程示意图。如图6所示,该获取方法可以包括如下步骤:
步骤411:基于骨头的CT值,获取CT图像中骨头区域图像。
根据CT图像中各个区域的CT值即可获取CT值最大的骨头区域图像。在一实施例中,可以基于骨头的CT值设定第一CT值阈值,获取CT图像中CT值大于或等于第一CT值阈值的连通区域作为骨头区域图像。其中,第一CT值阈值小于骨头的CT值且大于其他组织的CT值,通过设定第一CT值阈值,获取CT图像中CT值大于或等于第一CT值阈值的连通区域,即可得到骨头区域图像。
在一实施例中,在获取CT图像中CT值大于或等于第一CT值阈值的连通区域后,上述方法还可以包括:去除连通区域中面积小于预设的面积阈值的区域。第一CT值阈值设定的过大时可能会遗漏部分骨头区域,而第一CT值阈值设定的过小时,由于钙化点的CT值较大,又可能会出现肺部或者心脏区域存在一定的钙化点成为骨头区域图像中的干扰噪声,因此,需要将其去除。通常钙化点的面积较小,因此,可以通过去除连通区域中面积小于预设的面积阈值的区域来排出钙化点对骨头区域图像的干扰,其中面积阈值可以根据实际应用而预先设定。
在一实施例中,在获取CT图像中CT值大于或等于第一CT值阈值的连通区域后,上述方法还可以包括:去除位于粗分割图像内的连通区域。通过将粗分割图像内连通区域去除,可以避免肺部或者心脏区域内的钙化点对最终的分割结果造成的影响,从而可以提高后续的分割精度。
步骤412:基于肋骨的特性,分割出骨头区域图像中的肋骨区域图像。
在一实施例中,步骤412的实现方式可以是:将标准肋骨图像与骨头区域图像进行比对,选取骨头区域图像中与标准肋骨图像的相似度大于预设相似度的骨头区域作为肋骨区域图像。由于肋骨相对其他骨头比较规则,肋骨通常是规则的排列在肺部外侧且呈现一定的弧形且为左右对称结构,通过多平面重建的三维视图中可以区分出肋骨与其他骨头的区别,特别是在矢状位上,可以清晰获知肋骨的上述特性,即肋骨的排列和形状有其特殊性,因此,可以通过标准的肋骨图像去比对,选取骨头区域图像中与该标准肋骨图像的相似度满足一定要求(大于预设相似度)的骨头区域,以分割出骨头区域图像中的肋骨区域图像。应当理解,本申请实施例也可以根据实际应用场景的需求而选取其他获取肋骨区域图像的方式,例如可以根据肋骨的弧度是否在预设弧度范围内来判断,也可以根据肋骨之间的间隙是否在预设距离范围内(因为肋骨的排列有一定的规则,相邻的肋骨之间各处的垂直距离都在一定的范围内),还可以直接通过神经网络模型来获取肋骨区域图像,只要所选取的获取肋骨区域图像的方式满足精度需求即可,本申请实施例对于获取肋骨区域图像的具体方式不做限定。
图7是本申请另一示例性实施例提供的一种肺部区域分割方法的流程示意图。如图7所示,在步骤420之前,上述实施例还可以包括:
步骤450:对CT图像进行预处理。
在一实施例中,预处理可以包括以下操作中任一项或多项的组合:去除背景、去除白噪声、裁剪图像、变换窗宽和窗位。其中,去除背景的具体实现方式可以是:通过设置CT值范围,获取该CT值范围内的连通区域,并且只保留连通区域内面积最大的连通区域,其他区域均设置为背景区域,从而排除了其他区域的干扰。去除白噪声的具体实现方式可以是:通过高斯滤波器去除拍摄CT图像过程中引起的白噪声。裁剪图像的具体实现方式可以是:将背景去除,只保留有效区域,以降低后续图像处理的复杂度。变换窗宽和窗位的具体实现方式可以是:通过设定窗宽和窗位的值来重点突出感兴趣区域,从而避免不感兴趣的区域对后续处理的干扰,在本申请实施例中,可以选取窗位为-500、窗宽为1500,当然应当理解,窗位和窗宽的设定值可以根据实际情况而调整。
应当理解,步骤450可以设置于步骤410之前,通过预处理将CT图像中的背景和其他干扰因素都排出,可以有效降低后续步骤的复杂程度,提供肺部分割的效率。
图8是本申请另一示例性实施例提供的一种肺部分区域割方法的流程示意图。如图8所示,在步骤430之后,上述实施例还可以包括:
步骤460:对肺部区域图像的边界进行光滑化处理。
由于根据肋骨区域图像只能确定肺部区域的部分边界,并且是以同一预设的步长通过活动轮廓模型扩张,因而得到的肺部区域图像的边界可能不光滑,在得到肺部区域图像后,对该肺部区域图像的边界进行光滑化处理,可以得到更为精确的肺部区域图像。
示例性装置
图9是本申请一示例性实施例提供的一种肺炎征象的分割装置的结构示意图。如图9所示,该肺炎征象的分割装置90包括如下模块:
生成模块91,用于基于CT图像中的肺部区域图像,分别生成多个肺炎征象图像;以及组合模块92,用于将多个肺炎征象图像组合,得到肺炎综合征象图像;其中,生成模块91进一步配置为:将肺部区域图像分别输入多个神经网络模型,得到多个肺炎征象图像。
本申请提出的肺炎征象的分割装置,通过生成模块91将CT图像中的肺部区域图像分别输入多个神经网络模型,以此分别多个肺炎征象图像,然后组合模块92将多个肺炎征象图像进行组合,得到肺炎综合征象图像,利用神经网络模型可以高效且准确的获取大量的待检测CT图像中的肺炎征象图像,从而为后续诊断肺炎提供可靠的数据依据。
在一实施例中,肺炎征象图像可以包括如下征象图像中的任一种或多种的组合:肺实变图像、磨玻璃影图像、肿块图像、树芽征图像、结节图像、空洞图像、晕征图像。在一实施例中,神经网络模型可以是深度学习神经网络模型,优选地,该神经网络模型可以是Unet神经网络模型。
在一实施例中,如图9所示,生成模块91可以包括如下单元:输入单元911,用于逐次将多层二维图像中的一部分输入多个神经网络模型,得到分别对应多个肺炎征象图像的多层二维征象图像;叠加单元912,用于分别将对应同一肺炎征象图像的多层二维征象图像叠加,得到多个肺炎征象图像。
在一实施例中,如图9所示,生成模块91还可以包括:后处理单元913,用于对多个肺炎征象图像分别进行腐蚀膨胀操作。
在一实施例中,如图9所示,分割装置90还可以包括:获取模块93,用于获取CT图像中的肋骨区域图像;粗分割模块94,用于获取CT图像中肺部区域的粗分割图像;以及精分割模块95,用于以肋骨区域图像为肺部区域的边界,并且以粗分割图像为种子区域,以预设的步长向周围扩张至所述肺部区域的边界,得到肺部区域图像。
在一实施例中,粗分割模块94可以进一步配置为:将CT图像输入神经网络模型,得到肺部区域的粗分割图像。
在一实施例中,粗分割模块94可以进一步配置为:根据肺部区域的CT值,选取CT值在肺部区域的CT值范围内的区域作为粗分割图像。
在一实施例中,精分割模块95可以进一步配置为:以肋骨区域图像为边界、粗分割图像为感兴趣区域,并通过活动轮廓模型、以预设的步长对CT图像进行分割,以得到肺部区域图像。其中,活动轮廓模型可以包括LevelSet模型或Snake模型。
在一实施例中,如图9所示,该分割装置90还可以包括:腐蚀模块96,用于对粗分割图像进行腐蚀操作,得到腐蚀后的粗分割图像。且精分割模块95配置为:以肋骨区域图像为肺部区域的边界,并且以腐蚀后的粗分割图像为种子区域,以预设的步长向周围扩张至肺部区域的边界,得到肺部区域图像。
在一实施例中,如图9所示,获取模块93可以包括:骨头获取单元931,用于基于骨头的CT值,获取CT图像中骨头区域图像;肋骨获取单元932,用于基于肋骨的特性,分割出骨头区域图像中的肋骨区域图像。
在一实施例中,骨头获取单元931可以进一步配置为:基于骨头的CT值设定第一CT值阈值,获取CT图像中CT值大于或等于第一CT值阈值的连通区域作为骨头区域图像。
在一实施例中,获取模块93可以进一步配置为:在获取CT图像中CT值大于或等于第一CT值阈值的连通区域后,去除连通区域中面积小于预设的面积阈值的区域。
在一实施例中,获取模块93可以进一步配置为:在获取CT图像中CT值大于或等于第一CT值阈值的连通区域后,去除位于粗分割图像内的连通区域。
在一实施例中,肋骨获取单元932可以进一步配置为:将标准肋骨图像与骨头区域图像进行比对,选取骨头区域图像中与标准肋骨图像的相似度大于预设相似度的骨头区域作为肋骨区域图像。
在一实施例中,如图9所示,该分割装置90还可以包括:预处理模块97,用于对CT图像进行预处理。在一实施例中,预处理可以包括以下操作中任一项或多项的组合:去除背景、去除白噪声、裁剪图像、变换窗宽和窗位。
在一实施例中,如图9所示,该分割装置90还可以包括:光滑化模块98,用于对肺部区域图像的边界进行光滑化处理。
示例性电子设备
下面,参考图10来描述根据本申请实施例的电子设备。该电子设备可以是第一设备和第二设备中的任一个或两者、或与它们独立的单机设备,该单机设备可以与第一设备和第二设备进行通信,以从它们接收所采集到的输入信号。
图10图示了根据本申请实施例的电子设备的框图。
如图10所示,电子设备10包括一个或多个处理器11和存储器12。
处理器11可以是中央处理单元(CPU)或者具有数据处理能力和/或指令执行能力的其他形式的处理单元,并且可以控制电子设备10中的其他组件以执行期望的功能。
存储器12可以包括一个或多个计算机程序产品,所述计算机程序产品可以包括各种形式的计算机可读存储介质,例如易失性存储器和/或非易失性存储器。所述易失性存储器例如可以包括随机存取存储器(RAM)和/或高速缓冲存储器(cache)等。所述非易失性存储器例如可以包括只读存储器(ROM)、硬盘、闪存等。在所述计算机可读存储介质上可以存储一个或多个计算机程序指令,处理器11可以运行所述程序指令,以实现上文所述的本申请的各个实施例的分割方法以及/或者其他期望的功能。在所述计算机可读存储介质中还可以存储诸如输入信号、信号分量、噪声分量等各种内容。
在一个示例中,电子设备10还可以包括:输入装置13和输出装置14,这些组件通过总线系统和/或其他形式的连接机构(未示出)互连。
例如,在该电子设备是第一设备或第二设备时,该输入装置13可以是摄像头,用于捕捉图像的输入信号。在该电子设备是单机设备时,该输入装置13可以是通信网络连接器,用于从第一设备和第二设备接收所采集的输入信号。
此外,该输入设备13还可以包括例如键盘、鼠标等等。
该输出装置14可以向外部输出各种信息,包括确定出的距离信息、方向信息等。该输出设备14可以包括例如显示器、扬声器、打印机、以及通信网络及其所连接的远程输出设备等等。
当然,为了简化,图10中仅示出了该电子设备10中与本申请有关的组件中的一些,省略了诸如总线、输入/输出接口等等的组件。除此之外,根据具体应用情况,电子设备10还可以包括任何其他适当的组件。
示例性计算机程序产品和计算机可读存储介质
除了上述方法和设备以外,本申请的实施例还可以是计算机程序产品,其包括计算机程序指令,所述计算机程序指令在被处理器运行时使得所述处理器执行本说明书上述“示例性方法”部分中描述的根据本申请各种实施例的分割方法中的步骤。
所述计算机程序产品可以以一种或多种程序设计语言的任意组合来编写用于执行本申请实施例操作的程序代码,所述程序设计语言包括面向对象的程序设计语言,诸如Java、C++等,还包括常规的过程式程序设计语言,诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行、部分在用户计算设备上部分在远程计算设备上执行、或者完全在远程计算设备或服务器上执行。
此外,本申请的实施例还可以是计算机可读存储介质,其上存储有计算机程序指令,所述计算机程序指令在被处理器运行时使得所述处理器执行本说明书上述“示例性方法”部分中描述的根据本申请各种实施例的分割方法中的步骤。
所述计算机可读存储介质可以采用一个或多个可读介质的任意组合。可读介质可以是可读信号介质或者可读存储介质。可读存储介质例如可以包括但不限于电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。
以上结合具体实施例描述了本申请的基本原理,但是,需要指出的是,在本申请中提及的优点、优势、效果等仅是示例而非限制,不能认为这些优点、优势、效果等是本申请的各个实施例必须具备的。另外,上述公开的具体细节仅是为了示例的作用和便于理解的作用,而非限制,上述细节并不限制本申请为必须采用上述具体的细节来实现。
本申请中涉及的器件、装置、设备、系统的方框图仅作为例示性的例子并且不意图要求或暗示必须按照方框图示出的方式进行连接、布置、配置。如本领域技术人员将认识到的,可以按任意方式连接、布置、配置这些器件、装置、设备、系统。诸如“包括”、“包含”、“具有”等等的词语是开放性词汇,指“包括但不限于”,且可与其互换使用。这里所使用的词汇“或”和“和”指词汇“和/或”,且可与其互换使用,除非上下文明确指示不是如此。这里所使用的词汇“诸如”指词组“诸如但不限于”,且可与其互换使用。
还需要指出的是,在本申请的装置、设备和方法中,各部件或各步骤是可以分解和/或重新组合的。这些分解和/或重新组合应视为本申请的等效方案。
提供所公开的方面的以上描述以使本领域的任何技术人员能够做出或者使用本申请。对这些方面的各种修改对于本领域技术人员而言是非常显而易见的,并且在此定义的一般原理可以应用于其他方面而不脱离本申请的范围。因此,本申请不意图被限制到在此示出的方面,而是按照与在此公开的原理和新颖的特征一致的最宽范围。
为了例示和描述的目的已经给出了以上描述。此外,此描述不意图将本申请的实施例限制到在此公开的形式。尽管以上已经讨论了多个示例方面和实施例,但是本领域技术人员将认识到其某些变型、修改、改变、添加和子组合。
Claims (9)
1.一种肺炎征象的分割方法,其特征在于,包括:
基于CT图像中的肺部区域图像,分别生成多个肺炎征象图像;以及
将所述多个肺炎征象图像组合,得到肺炎综合征象图像;
其中,分别生成多个肺炎征象图像的方式包括:
将所述肺部区域图像分别输入多个神经网络模型,得到所述多个肺炎征象图像;
所述肺部区域图像的获取方式包括:
获取所述CT图像中的肋骨区域图像;
获取所述CT图像中肺部区域的粗分割图像;以及
以所述肋骨区域图像为所述肺部区域的边界,并且以所述粗分割图像为种子区域,以预设的步长向周围扩张至所述肺部区域的边界,得到所述肺部区域图像。
2.根据权利要求1所述的分割方法,其特征在于,所述多个肺炎征象图像包括如下征象图像中的任一种或多种的组合:
肺实变图像、磨玻璃影图像、肿块图像、树芽征图像、结节图像、空洞图像、晕征图像。
3.根据权利要求1所述的分割方法,其特征在于,所述肺部区域图像包括多层二维图像;所述将所述肺部区域图像分别输入多个神经网络模型,得到所述多个肺炎征象图像包括:
逐次将所述多层二维图像中的一部分分别输入所述多个神经网络模型,得到分别对应所述多个肺炎征象图像的多层二维征象图像;以及
分别将对应同一肺炎征象图像的所述多层二维征象图像叠加,得到所述多个肺炎征象图像。
4.根据权利要求1所述的分割方法,其特征在于,所述肺部区域图像包括多层二维图像;所述将所述肺部区域图像分别输入多个神经网络模型,得到所述多个肺炎征象图像包括:
将所述多层二维图像分别输入所述多个神经网络模型,得到分别对应所述多个肺炎征象图像的多层二维征象图像;以及
分别将对应同一肺炎征象图像的所述多层二维征象图像叠加,得到所述多个肺炎征象图像。
5.根据权利要求3或4所述的分割方法,其特征在于,在所述得到所述多个肺炎征象图像之后,还包括:
对所述多个肺炎征象图像分别进行腐蚀膨胀操作。
6.根据权利要求1所述的分割方法,其特征在于,在所述获取所述CT图像中肺部区域的粗分割图像之后,还包括:
对所述粗分割图像进行腐蚀操作,得到腐蚀后的粗分割图像;
所述以所述肋骨区域图像为所述肺部区域的边界,并且以所述粗分割图像为种子区域,以预设的步长向周围扩张至所述肺部区域的边界,得到所述肺部区域图像包括:
以所述肋骨区域图像为所述肺部区域的边界,并且以所述腐蚀后的粗分割图像为种子区域扩张,以预设的步长向周围扩张至所述肺部区域的边界,得到所述肺部区域图像。
7.一种肺炎征象的分割装置,其特征在于,包括:
获取模块,用于获取CT图像中的肋骨区域图像;
粗分割模块,用于获取所述CT图像中肺部区域的粗分割图像;
精分割模块,用于以所述肋骨区域图像为所述肺部区域的边界,并且以所述粗分割图像为种子区域,以预设的步长向周围扩张至所述肺部区域的边界,得到所述肺部区域图像;
生成模块,用于基于所述CT图像中的所述肺部区域图像,分别生成多个肺炎征象图像;以及
组合模块,用于将所述多个肺炎征象图像组合,得到肺炎综合征象图像;
其中,生成模块进一步配置为:
将所述肺部区域图像分别输入多个神经网络模型,得到所述多个肺炎征象图像。
8.一种计算机可读存储介质,所述存储介质存储有计算机程序,所述计算机程序用于执行上述权利要求1-6任一所述的分割方法。
9.一种电子设备,所述电子设备包括:
处理器;
用于存储所述处理器可执行指令的存储器;
所述处理器,用于执行上述权利要求1-6任一所述的分割方法。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010189089.6A CN111415356B (zh) | 2020-03-17 | 2020-03-17 | 肺炎征象的分割方法、装置、介质及电子设备 |
PCT/CN2021/072515 WO2021184949A1 (zh) | 2020-03-17 | 2021-01-18 | 肺炎征象的分割方法、装置、介质及电子设备 |
EP21771562.2A EP3971830B1 (en) | 2020-03-17 | 2021-01-18 | Pneumonia sign segmentation method and apparatus, medium and electronic device |
JP2021571620A JP7304437B2 (ja) | 2020-03-17 | 2021-01-18 | 肺炎徴候の分割方法、装置、媒体及び電子デバイス |
US17/540,629 US20220092788A1 (en) | 2020-03-17 | 2021-12-02 | Segmentation method for pneumonia sign, medium, and electronic device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010189089.6A CN111415356B (zh) | 2020-03-17 | 2020-03-17 | 肺炎征象的分割方法、装置、介质及电子设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111415356A CN111415356A (zh) | 2020-07-14 |
CN111415356B true CN111415356B (zh) | 2020-12-29 |
Family
ID=71494475
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010189089.6A Active CN111415356B (zh) | 2020-03-17 | 2020-03-17 | 肺炎征象的分割方法、装置、介质及电子设备 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20220092788A1 (zh) |
EP (1) | EP3971830B1 (zh) |
JP (1) | JP7304437B2 (zh) |
CN (1) | CN111415356B (zh) |
WO (1) | WO2021184949A1 (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111415356B (zh) * | 2020-03-17 | 2020-12-29 | 推想医疗科技股份有限公司 | 肺炎征象的分割方法、装置、介质及电子设备 |
WO2022225794A1 (en) * | 2021-04-23 | 2022-10-27 | The Johns Hopkins University | Systems and methods for detecting and characterizing covid-19 |
CN115496703A (zh) * | 2021-06-18 | 2022-12-20 | 富联精密电子(天津)有限公司 | 肺炎区域检测方法及系统 |
CN114240956A (zh) * | 2021-12-23 | 2022-03-25 | 零氪科技(天津)有限公司 | 一种肺部ct影像管状结构分割方法、装置、电子设备 |
CN114723688B (zh) * | 2022-03-24 | 2024-08-16 | 华中科技大学同济医学院附属协和医院 | 一种用于ct影像肺炎识别的影像处理方法和系统 |
WO2024180773A1 (ja) * | 2023-03-02 | 2024-09-06 | オリンパスメディカルシステムズ株式会社 | 画像処理装置、画像処理方法、および、画像処理プログラム |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109727253A (zh) * | 2018-11-14 | 2019-05-07 | 西安大数据与人工智能研究院 | 基于深度卷积神经网络自动分割肺结节的辅助检测方法 |
Family Cites Families (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2427695C (en) * | 2001-04-26 | 2012-12-11 | Teijin Limited | Method for measuring 3d joint structure |
US20060018524A1 (en) * | 2004-07-15 | 2006-01-26 | Uc Tech | Computerized scheme for distinction between benign and malignant nodules in thoracic low-dose CT |
WO2010009735A2 (en) * | 2008-07-23 | 2010-01-28 | Dako Denmark A/S | Combinatorial analysis and repair |
US20140074493A1 (en) * | 2012-01-03 | 2014-03-13 | David Lyle Schneider | Medical protection lockout system for procedures and devices |
CA2924228C (en) * | 2013-10-02 | 2024-01-16 | Boehringer Ingelheim Vetmedica, Inc. | Pcv2 orf2 protein variant and virus like particles composed thereof |
JP6906347B2 (ja) * | 2017-04-10 | 2021-07-21 | 富士フイルム株式会社 | 医用画像分類装置、方法およびプログラム |
CN107220373B (zh) * | 2017-06-19 | 2019-09-17 | 太原理工大学 | 一种基于医学征象和卷积神经网络的肺结节ct图像哈希检索方法 |
CN107392910B (zh) * | 2017-07-06 | 2020-01-07 | 东软医疗系统股份有限公司 | 一种基于ct图像的肺叶分割方法及装置 |
CN107563383A (zh) * | 2017-08-24 | 2018-01-09 | 杭州健培科技有限公司 | 一种医学影像辅助诊断及半监督样本生成系统 |
CN109493347B (zh) * | 2017-09-12 | 2021-03-23 | 深圳科亚医疗科技有限公司 | 在图像中对稀疏分布的对象进行分割的方法和系统 |
CN107845091B (zh) * | 2017-10-31 | 2019-02-15 | 北京青燕祥云科技有限公司 | 基于预筛查的快速标注方法和系统 |
JP2019169049A (ja) | 2018-03-26 | 2019-10-03 | 富士フイルム株式会社 | 医用画像特定装置、方法およびプログラム |
US10304193B1 (en) * | 2018-08-17 | 2019-05-28 | 12 Sigma Technologies | Image segmentation and object detection using fully convolutional neural network |
CN110503640B (zh) * | 2018-08-21 | 2022-03-22 | 深圳科亚医疗科技有限公司 | 对医学图像进行分析的装置、系统及计算机可读介质 |
CN109636802A (zh) * | 2019-01-18 | 2019-04-16 | 天津工业大学 | 基于深度卷积神经网络的肺实质经ct图像分割方法 |
CN110517225B (zh) * | 2019-07-19 | 2023-07-11 | 平安科技(深圳)有限公司 | Ai影像识别方法、装置、设备及存储介质 |
CN110599448B (zh) * | 2019-07-31 | 2022-03-15 | 浙江工业大学 | 基于MaskScoring R-CNN网络的迁移学习肺部病变组织检测系统 |
CN110706200B (zh) * | 2019-09-02 | 2022-08-05 | 杭州深睿博联科技有限公司 | 数据预测的方法及装置 |
CN111415356B (zh) * | 2020-03-17 | 2020-12-29 | 推想医疗科技股份有限公司 | 肺炎征象的分割方法、装置、介质及电子设备 |
CN111415743B (zh) * | 2020-03-17 | 2021-07-09 | 推想医疗科技股份有限公司 | 肺炎分类方法、装置、介质及电子设备 |
-
2020
- 2020-03-17 CN CN202010189089.6A patent/CN111415356B/zh active Active
-
2021
- 2021-01-18 EP EP21771562.2A patent/EP3971830B1/en active Active
- 2021-01-18 JP JP2021571620A patent/JP7304437B2/ja active Active
- 2021-01-18 WO PCT/CN2021/072515 patent/WO2021184949A1/zh unknown
- 2021-12-02 US US17/540,629 patent/US20220092788A1/en not_active Abandoned
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109727253A (zh) * | 2018-11-14 | 2019-05-07 | 西安大数据与人工智能研究院 | 基于深度卷积神经网络自动分割肺结节的辅助检测方法 |
Non-Patent Citations (2)
Title |
---|
P. Annangi et al..A REGION BASED ACTIVE CONTOUR METHOD FOR X-RAY LUNG SEGMENTATION USING PRIOR SHAPE AND LOW LEVEL FEATURES.《2010 IEEE International Symposium on Biomedical Imaging:From Nano to Macro》.2010, * |
孟琭.基于活动轮廓模型的医学图像分割关键技术研究.《中国博士学位论文全文数据库 信息科技辑》.2013,(第3期), * |
Also Published As
Publication number | Publication date |
---|---|
WO2021184949A1 (zh) | 2021-09-23 |
EP3971830A4 (en) | 2022-09-14 |
CN111415356A (zh) | 2020-07-14 |
EP3971830A1 (en) | 2022-03-23 |
JP7304437B2 (ja) | 2023-07-06 |
US20220092788A1 (en) | 2022-03-24 |
EP3971830C0 (en) | 2023-09-20 |
JP2022535811A (ja) | 2022-08-10 |
EP3971830B1 (en) | 2023-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111415356B (zh) | 肺炎征象的分割方法、装置、介质及电子设备 | |
CN111415743B (zh) | 肺炎分类方法、装置、介质及电子设备 | |
Pu et al. | CT based computerized identification and analysis of human airways: a review | |
Hua et al. | Segmentation of pathological and diseased lung tissue in CT images using a graph-search algorithm | |
Pu et al. | A computational geometry approach to automated pulmonary fissure segmentation in CT examinations | |
US8605978B2 (en) | Medical image processing apparatus and method, and computer readable recording medium on which is recorded program for the same | |
Swift et al. | Automatic axis generation for virtual bronchoscopic assessment of major airway obstructions | |
JP2005199057A (ja) | 肺画像データの位置合わせのための方法及び装置 | |
CN113327225B (zh) | 用于提供气道信息的方法 | |
JP2007526033A (ja) | 構造化オブジェクトの画像を登録するための装置及び方法 | |
CN111260669A (zh) | 基于ct图像的肺叶分割方法和装置 | |
JP7101809B2 (ja) | 画像処理装置、画像処理方法、及びプログラム | |
US11475568B2 (en) | Method for controlling display of abnormality in chest x-ray image, storage medium, abnormality display control apparatus, and server apparatus | |
WO2009063363A2 (en) | Computer-aided detection (cad) of a disease | |
WO2022164374A1 (en) | Automated measurement of morphometric and geometric parameters of large vessels in computed tomography pulmonary angiography | |
CN111311612A (zh) | 肺部分割方法、装置、介质及电子设备 | |
Zhai et al. | Automatic quantitative analysis of pulmonary vascular morphology in CT images | |
JP5048233B2 (ja) | Cadシステムにおける解剖学的形状の検出のための方法及びシステム | |
US11935234B2 (en) | Method for detecting abnormality, non-transitory computer-readable recording medium storing program for detecting abnormality, abnormality detection apparatus, server apparatus, and method for processing information | |
JP2020032043A (ja) | 画像処理装置、方法およびプログラム | |
CN111415341A (zh) | 肺炎阶段的评估方法、装置、介质及电子设备 | |
Rebouças Filho et al. | 3D segmentation and visualization of lung and its structures using CT images of the thorax | |
Wang et al. | Accuracy and reliability analysis of a machine learning based segmentation tool for intertrochanteric femoral fracture CT | |
CN105997119A (zh) | 用于气道壁密度和气道壁炎症的同时评估的方法和系统 | |
CN117475250A (zh) | 基于解剖结构数据模拟病理图像 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: Room B401, floor 4, building 1, No. 12, Shangdi Information Road, Haidian District, Beijing 100085 Applicant after: Tuxiang Medical Technology Co., Ltd Address before: Room B401, floor 4, building 1, No. 12, Shangdi Information Road, Haidian District, Beijing 100085 Applicant before: Beijing Tuoxiang Technology Co.,Ltd. |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |