CN111967540B - 基于ct数据库的颌面部骨折的识别方法、装置及终端设备 - Google Patents
基于ct数据库的颌面部骨折的识别方法、装置及终端设备 Download PDFInfo
- Publication number
- CN111967540B CN111967540B CN202011046297.7A CN202011046297A CN111967540B CN 111967540 B CN111967540 B CN 111967540B CN 202011046297 A CN202011046297 A CN 202011046297A CN 111967540 B CN111967540 B CN 111967540B
- Authority
- CN
- China
- Prior art keywords
- fracture
- data
- maxillofacial
- training
- model
- 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
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- 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
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Public Health (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Radiology & Medical Imaging (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Evolutionary Biology (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Quality & Reliability (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Databases & Information Systems (AREA)
- Pathology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及图像处理领域,其实施方式提供了一种基于CT数据库的颌面部骨折的识别方法,所述方法包括:获取不同骨折类型的颌面部骨折的CT数据样本;对所述CT数据样本进行标注处理和分类处理;采用同一骨折类型的经处理的CT数据样本训练骨折判别模型,训练好的骨折判别模型通过识别输入图像块是否包括某一骨折类型对应的图像特征来确定是否存在对应的骨折类型;将输入的颌面部CT数据输入对应的训练好的骨折判别模型,以确定输入的颌面部CT数据是否存在对应的骨折类型。同时还提供了对应的基于CT数据库的颌面部骨折的识别装置,以及一种终端设备。本发明提供的实施方式适用于医学图像中的CT图像中的骨折识别,提升了识别效率。
Description
技术领域
本发明涉及图像处理领域,特别涉及一种基于CT数据库的颌面部骨折的识别方法、一种基于CT数据库的颌面部骨折的识别装置、一种终端设备以及对应的存储介质。
背景技术
颌面部骨折是交通事故、意外伤害、竞技运动等事故中常见的外伤类型。颌面部骨折的解剖结构复杂多样,仅依靠临床症状及体征的诊断难度较大,常需要影像学的辅助诊断。相比于常规的X线片,颌面部CT的三维影像能够更加清晰、直观地表现骨折的位置及移位的方向,近年来被认为是颌面部骨折诊断的“金标准”。但是CT图像中包含大量信息,仅通过临床医生的视觉检查难以全面准确评估诸如骨骼结构、边界和隐匿骨折线等细节。因此,传统诊疗模式在处理以“快速、准确、高效”为特点的大型赛事及突发事件现场的救治存在较大困难。
近年来,深度学习技术已逐步用于医学领域,在检测癌症、白内障、骨折、脑出血等疾病方面均有较好的效果。其中卷积神经网络(CNN)是医学图像诊断方面最尖端的技术,其高精确性和稳定性弥补了人眼诊断漏诊和误诊的不足,在肺结核、肺结节CT图像、乳腺癌、脑部病变及白内障分级等疾病分类的精确度已证明达到了人类专家的级别。
现有的基于CT影像的颌面骨折识别技术都是依靠影像科的专业医师借助相关软件(如mimics research 19.0.)进行人工判断,骨折类型复杂包括上颌骨骨折、颧骨骨折、下颌体骨折、下颌角及升支骨折、牙槽突骨折、颏部骨折、髁突骨折和喙突骨折,导致骨折的人工诊断折效率不高。
发明内容
有鉴于此,本发明旨在提出一种基于CT数据库的颌面部骨折的识别方法、装置及终端设备,以至少部分地解决上述问题。本发明拟基于卷积神经网络算法,对颌面部骨折CT数据模型进行深度学***台,提高疾病诊治的稳定性和反应效率。本发明解决了颌面骨折识别的效率问题,通过自动检测发生骨折的断层位置,来辅助医生提高诊断效率。
在本发明的第一方面,提供了一种基于CT数据库的颌面部骨折的识别方法,所述方法包括:获取不同骨折类型的颌面部骨折的CT数据样本;对所述CT数据样本进行标注处理和分类处理;采用同一骨折类型的经处理的CT数据样本训练骨折判别模型,训练好的骨折判别模型通过识别输入图像块是否包括某一骨折类型对应的图像特征来确定是否存在对应的骨折类型;所述骨折类型按解剖区域进行划分;将输入的颌面部CT数据按解剖结构分解为不同解剖区域的图像块序列;将所述图像块序列中的每个图像块分别输入对应的训练好的骨折判别模型,以确定输入的颌面部CT数据是否存在对应的骨折类型。
可选的,所述将待识别的颌面部CT数据按解剖结构分解为不同解剖区域的图像块序列,包括:将所述颌面部CT数据的CT图像序列输入训练好的骨结构定位模型;生成检测矩形框序列;根据不同解剖区域的矩形框,定位每个解剖区域的断层范围;分别从定位出的断层范围内的CT图像序列中逐层提取出每个解剖区域的图像块序列。
可选的,所述骨结构定位模型为卷积神经网络模型,所述卷积神经网络模型的输入为所述颌面部CT数据的CT图像序列,输出为对应的矩形框检测结果。
可选的,所述训练好的骨结构定位模型通过以下步骤得到:对所述CT数据样本进行所述标注处理后作为第一训练样本;将所述第一训练样本分为训练集和验证集;将所述骨结构定位模型的参数进行预设后,采用所述训练集中的训练样本和梯度下降算法对所述骨结构定位模型进行迭代训练;采用所述验证集中的训练样本确定所述骨结构定位模型的最优参数。
可选的,在所述将所述训练样本分为训练集和验证集之前,所述方法还包括:将所述CT图像序列的亨氏单位值映射至预设亨氏单位值范围。
可选的,所述骨折判别模型为:EfficientNetb3、ResNet、DenseNet、3D-ResNet中的一种。
可选的,所述训练好的骨折判别模型通过以下步骤得到:确定所述骨折判别模型的模型结构;将所述第一训练样本按解剖结构分解为不同解剖区域的图像块序列,每个解剖区域的图像块序列构成一个分区训练样本;每个不同解剖区域的骨折判别模型均采用以下步骤进行训练:将对应解剖区域的分区训练样本分为训练集和验证集;将所述骨折判别模型的参数进行预设后,采用所述训练集中的训练样本和梯度下降算法对所述骨折判别模型进行迭代训练;采用所述验证集中的训练样本确定所述骨折判别模型的最优参数。
在本发明的第二方面,还提供了一种基于CT数据库的颌面部骨折的识别装置,所述装置包括:样本获取模块,用于获取不同骨折类型的颌面部骨折的CT数据样本;样本处理模块,用于对所述CT数据样本进行标注处理和分类处理;特征提取模块,用于从同一骨折类型的经过处理的CT数据样本中提取所述骨折类型的图像特征;以及骨折识别模块,用于将所述图像特征与输入待识别的颌面部CT数据进行比对,以确定输入的颌面部CT数据是否存在与所述图像特征对应的骨折类型。
在本发明第三方面,还提供了一种终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现前述的基于CT数据库的颌面部骨折的识别方法的步骤。
在本发明第四方面,还提供了一种计算机可读存储介质,所述存储介质中存储有指令,当其在计算机上运行时,使得计算机执行前述的基于CT数据库的颌面部骨折的识别方法。
通过本发明提供的上述技术方案,具有以下有益效果:与现有技术相比,本发明的有益效果是:本发明提出的基于螺旋CT数据库的颌面部骨折人工智能诊断方法,基于卷积神经网络算法,对颌面部骨折CT数据模型进行深度学***台,提高了疾病诊治的稳定性和反应效率。解决了传统诊疗依赖于专业医师和特定诊疗场所的局限性,使得在诸如冬奥会等要求快速反应、准确判断、及时就地处理的事故现场的即时诊治成为可能。
本发明的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施方式及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明一种实施方式提供的基于CT数据库的颌面部骨折的识别方法的流程示意图;
图2是本发明一种实施方式提供的基于CT数据库的颌面部骨折的识别装置的结构示意图;
图3是本发明一种实施方式提供的终端设备的示意图。
具体实施方式
需要说明的是,在不冲突的情况下,本发明中的实施方式及实施方式中的特征可以相互组合。
以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本发明,并不用于限制本发明。
图1是本发明一种实施方式提供的基于CT数据库的颌面部骨折的识别方法的流程示意图,如图1所示。提供了一种基于CT数据库的颌面部骨折的识别方法,所述方法包括:获取不同骨折类型的颌面部骨折的CT数据样本;对所述CT数据样本进行标注处理和分类处理;采用同一骨折类型的经处理的CT数据样本训练骨折判别模型,训练好的骨折判别模型通过识别输入图像块是否包括某一骨折类型对应的图像特征来确定是否存在对应的骨折类型。
如此,基于颌面部骨折的CT数据所构成的大数据样本,提取每种骨折类型的图像特征,以此实现识别所述颌面部CT数据中是否存在骨折,以及存在的骨折类型。本发明提供的实施方式可以通过训练好的机器学习的图像处理特性,将CT数据进行按类型分解和判断,具有大数据处理的优点,同时避免了对颌面部CT数据直接识别造成的准确度低的问题。优选的,其中分类处理和比对处理均通过卷积神经网络实现。
具体的,样本数据库的建立。基于已收集的左/右喙突骨折、左/右升支骨折、左/右下颌角骨折、左/右下颌体骨折、左/右髁突骨折、颏部骨折、牙槽突骨折、上颌骨骨折、左/右颧骨骨折患者的CT数据,将患者的CT数据进行标注。获得逐层的mask标注和骨折类型的标注。将所有标注的DICOM文件以及标注表导出。此步骤又可以分为以下三个子步骤:
1-1:纳入标准
半个月内具有颌面部外伤骨折史的18-80岁汉族国人。纳入标准为:1)汉族成年人,年龄18~80岁;2)半个月内具有颌面部外伤骨折史;3)颌面部无严重肿瘤病史;4)无全身骨代谢疾病;5)无颌面部发育畸形;6)无放化疗史。
骨折类型主要包括:左/右喙突骨折、左/右升支骨折、左/右下颌角骨折、左/右下颌体骨折、左/右髁突骨折、颏部骨折、牙槽突骨折、上颌骨骨折、左/右颧骨骨折。
排除标准为:1)存在先天面部不对称,如严重的偏颌、鼻中隔偏曲、小耳畸形;2)有颌面部硬组织手术史;3)陈旧性骨折或青枝骨折;4)全身情况较差,不能耐受转诊或无法坐立。5)孕期1-3个月的妇女。
1-2:获得标注
标注的内容主要包括逐层的mask标注和骨折类型的标注。逐层标注将受试者CT数据的DICOM格式导入Mimics Medical 21.0软件进行分割,将分割出的颧上颌骨及下颌骨mask以DICOM格式导出。并观察含有骨折线的CT断层,标注并导出;骨折类型标注将骨折类型分为14类,对应每个病例,阳性结果的骨折类型标为“1”,阴性结果标为“0”,将所有标注的DICOM文件以及标注表导出。
1-3:导出标注及原始CT的DICOM格式文件。
通过以上步骤形成样本库,并以该样本库中样本形成不同骨折类型的图像特征。基于每种骨折类型的图像特征,比对判断是否存在该骨折类型。
在本发明提供的一种实施方式中,所述骨折类型按解剖区域划分;将输入的颌面部CT数据按解剖结构分解为不同解剖区域的图像块序列;将所述图像块序列中的每个图像块分别输入对应的训练好的骨折判别模型来确定输入的颌面部CT数据是否存在对应的骨折类型。将颌面部CT数据分解为多个解剖结构上的构件,通过识别单个构件是否骨折来识别所述颌面部CT数据中是否存在骨折,以此实现对颌面部整体判断,还能识别出骨折位置。本发明提供的实施方式通过训练好的机器学习的图像处理特性,将CT数据进行按区域分解和判断,具有处理快速和信息损失小的优点,同时避免了对颌面部CT数据整体判断造成的准确度低的问题。本实施方式算法主要包含两个模型:骨结构定位模型和骨折判别模型,定位模型主要负责找到上颌骨和下颌骨的具***置及轮廓,提取出不同的关键部位,然后将提取出的不同部位输入进骨折判别模型进行骨折类型的识别。
在本发明提供的一种实施方式中,所述将待识别的颌面部CT数据按解剖结构分解为不同解剖区域的图像块序列,包括:将所述颌面部CT数据的CT图像序列输入训练好的骨结构定位模型;生成检测矩形框序列;根据不同解剖区域的矩形框,定位每个解剖区域的断层范围;分别从定位出的断层范围内的CT图像序列中逐层提取出每个解剖区域的图像块序列。将CT图像序列依次输入训练好的检测模型,输出检测矩形框序列,根据检测出的不同类别的矩形框,定位8种解剖区域的断层范围,然后8种解剖区域分别从定位出的断层范围内的CT图像序列中逐层提取出来,由于颧骨、下颌体、下颌角及升支、髁突和喙突结构具有左右对称性,所以最终将提取出13个图像块序列。
在本发明提供的一种实施方式中,所述骨结构定位模型为卷积神经网络模型,所述卷积神经网络模型的输入为所述颌面部CT数据的CT图像序列,输出为对应的矩形框检测结果。具体的,本实施方式可采用RetinaNet卷积神经网络结构,模型的输入为单张CT影像序列,输出为对应的矩形框检测结果。其中,检测模型RetinaNet可以用其它模型替换,比如SSD,YOLO,FasterRCNN等目标检测模型。
在本发明提供的一种实施方式中,所述训练好的骨结构定位模型通过以下步骤得到:对所述CT数据样本进行所述标注处理后作为第一训练样本;将所述第一训练样本分为训练集和验证集;将所述骨结构定位模型的参数进行预设后,采用所述训练集中的训练样本和梯度下降算法对所述骨结构定位模型进行迭代训练;采用所述验证集中的训练样本确定所述骨结构定位模型的最优参数。通过以上步骤形成训练样本,然后采用该训练样本对模型进行训练。将训练样本的10%,即50例样本,作为验证集,其余作为训练集;模型参数采用在大型自然图像数据集ImageNet上预训练好的参数,并采用梯度下降算法SGD对模型进行迭代训练,根据在验证集上的mAP值来确定模型的最优参数。
在本发明提供的一种实施方式中,在所述将所述训练样本分为训练集和验证集之前,所述方法还包括:将所述CT图像序列的亨氏单位值映射至预设亨氏单位值范围。在本发明提供的一种实施方式中,所述将CT图像的亨氏单位值值映射至预设范围,包括:将CT图像的亨氏单位值按照线性映射至[0, 255],映射公式如下:
y=(x-xmin)/(xmax-xmin)*255
其中,y表示映射后的值,x表示原始亨氏单位值,xmin表示CT图像中的最小的亨氏单位,xmax表示CT最大的亨氏单位值。本实施方式通过将亨氏单位值压缩为8位,以实现训练图像在处理中的一致性。
在本发明提供的一种实施方式中,所述骨折判别模型为:EfficientNetb3、ResNet、DenseNet、3D-ResNet中的一种。以上模型均为卷积神经网络结构,模型的输入为单张区域图像块,输出为对应的二值判别结果。
在本发明提供的一种实施方式中,所述训练好的骨折判别模型通过以下步骤得到:确定所述骨折判别模型的模型结构;将所述第一训练样本按解剖结构分解为不同解剖区域的图像块序列,每个解剖区域的图像块序列构成一个分区训练样本;每个不同解剖区域的骨折判别模型均采用以下步骤进行训练:将对应解剖区域的分区训练样本分为训练集和验证集;将所述骨折判别模型的参数进行预设后,采用所述训练集中的训练样本和梯度下降算法对所述骨折判别模型进行迭代训练;采用所述验证集中的训练样本确定所述骨折判别模型的最优参数。具体的,骨折检测步骤包括:根据前文提取出来的区域图像块序列和对应的标注进行深度学习判别模型训练,进行各型骨折的自动识别;其中对骨折判别模型进行模型训练包括:将提取出来的区域图像块序列的80%作为训练集,20%作为验证集;模型参数采用在大型自然图像数据集ImageNet上预训练好的参数,并采用梯度下降算法Adam对模型进行迭代训练,根据在验证集上的F1值确定模型最优参数。
具体的各型骨折检测结果生成:将13个区域图像块序列依次输入训练好的判别模型,输出形成二值判别结果序列。对于具有左右对称性的骨结构,判别结果序列可以看成左右2个二值序列(0表示无骨折,1表示有骨折);不具有左右对称性的骨结构,其判别结果可以看成1个二值序列,最终可获得13个二值序列;对于任意一个二值序列,根据骨折连续性的特点并保证高敏感性的前提下制定如下规则:如果二值序列中连续的元素1的数量大于1,则该例样本的该类型骨折识别结果为阳性,并输出对应的阳性层图像,否则为阴性。以及进一步的,由两名颌面外科住院医师对测试组数据进行标注,作为骨折诊断的“金标准”。将测试组数据机器的输出结果,与“金标准”进行比较,计算各部位骨折诊断的灵敏度、特异度。
在本发明提供的一种实施方式中,还提供了基于CT数据库的颌面部骨折的识别装置。图2是本发明一种实施方式提供的基于CT数据库的颌面部骨折的识别装置的结构示意图,如图2所示。所述装置包括:样本获取模块,用于获取不同骨折类型的颌面部骨折的CT数据样本;样本处理模块,用于对所述CT数据样本进行标注处理和分类处理;特征提取模块,用于从同一骨折类型的经过处理的CT数据样本中提取所述骨折类型的图像特征;以及骨折识别模块,用于将所述图像特征与输入待识别的颌面部CT数据进行比对,以确定输入的颌面部CT数据是否存在与所述图像特征对应的骨折类型。
关于基于CT数据库的颌面部骨折的识别装置的具体限定可以参见上文中对于将基于CT数据库的颌面部骨折的识别方法的限定,在此不再赘述。上述装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在本发明的一种实施方式中,还提供了一种终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现前述的基于CT数据库的颌面部骨折的识别方法的步骤。
图3是本发明一种实施方式提供的终端设备的示意图,如图3所示。所述终端设备10可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。终端设备10可包括,但不仅限于,处理器100、存储器101。本领域技术人员可以理解,图3仅仅是终端设备10的示例,并不构成对终端设备10的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件,例如所述终端设备还可以包括输入输出设备、网络接入设备、总线等。
所述处理器100可以是中央处理单元(Central Processing Unit,CPU),还可以是其他通用处理器、数字信号处理器 (Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列 (Field-Programmable Gate Array,FPGA) 或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
所述存储器101可以是终端设备10的内部存储单元,例如终端设备10的硬盘或内存。所述存储器101也可以是终端设备10的外部存储设备,例如所述终端设备10上配备的插接式硬盘,智能存储卡(Smart Media Card, SMC),安全数字(Secure Digital, SD)卡,闪存卡(Flash Card)等。进一步地,所述存储器101还可以既包括终端设备10的内部存储单元也包括外部存储设备。所述存储器101用于存储所述计算机程序102以及终端设备10所需的其他程序和数据。所述存储器101还可以用于暂时地存储已经输出或者将要输出的数据。
本发明提供的实施方式针对CT数据的骨折识别存在的识别复杂和准确率低的问题,提供了一种基于CT数据库的颌面部骨折的识别方法及装置,该方法采用当前的大数据技术及机器学习模型,提高了图像中骨折识别的效率和准确率。本发明的提供的实施方式应用于医学图像处理***。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
存储器可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (8)
1.一种基于CT数据库的颌面部骨折的识别方法,其特征在于,所述方法包括:
获取不同骨折类型的颌面部骨折的CT数据样本;
对所述CT数据样本进行标注处理和分类处理;
采用同一骨折类型的经处理的CT数据样本训练骨折判别模型,训练好的骨折判别模型通过识别输入图像块是否包括某一骨折类型对应的图像特征来确定是否存在对应的骨折类型;所述骨折类型按解剖区域进行划分;
将输入的颌面部CT数据按解剖结构分解为不同解剖区域的图像块序列;
将所述图像块序列中的每个图像块分别输入对应的训练好的骨折判别模型,以确定输入的颌面部CT数据是否存在对应的骨折类型;
所述将输入的颌面部CT数据按解剖结构分解为不同解剖区域的图像块序列,包括:
将所述颌面部CT数据的CT图像序列输入到训练好的骨结构定位模型;
生成检测矩形框序列;
根据不同解剖区域的矩形框,定位每个解剖区域的断层范围;
分别从定位出的断层范围内的CT图像序列中逐层提取出每个解剖区域的图像块序列。
2.根据权利要求1所述的方法,其特征在于,所述训练好的骨结构定位模型为卷积神经网络模型,所述卷积神经网络模型的输入为所述颌面部CT数据的CT图像序列,输出为对应的矩形框检测结果。
3.根据权利要求2所述的方法,其特征在于,所述训练好的骨结构定位模型通过以下步骤得到:
对所述CT数据样本进行所述标注处理后作为第一训练样本;
将所述第一训练样本分为训练集和验证集;
将所述骨结构定位模型的参数进行预设后,采用所述训练集中的训练样本和梯度下降算法对所述骨结构定位模型进行迭代训练;
采用所述验证集中的训练样本确定所述骨结构定位模型的最优参数。
4.根据权利要求3所述的方法,其特征在于,在所述将所述第一训练样本分为训练集和验证集之前,所述方法还包括:
将所述CT图像序列的亨氏单位值映射至预设亨氏单位值范围。
5.根据权利要求1所述的方法,其特征在于,所述骨折判别模型为:EfficientNetb3、ResNet、DenseNet、3D-ResNet中的一种。
6.根据权利要求3所述的方法,其特征在于,所述训练好的骨折判别模型通过以下步骤得到:
确定所述骨折判别模型的模型结构;
将所述第一训练样本按解剖结构分解为不同解剖区域的图像块序列,每个解剖区域的图像块序列构成一个分区训练样本;
每个不同解剖区域的骨折判别模型均采用以下步骤进行训练:
将对应解剖区域的分区训练样本分为训练集和验证集;
将所述骨折判别模型的参数进行预设后,采用所述训练集中的训练样本和梯度下降算法对所述骨折判别模型进行迭代训练;
采用所述验证集中的训练样本确定所述骨折判别模型的最优参数。
7.一种基于CT数据库的颌面部骨折的识别装置,其特征在于,所述装置包括:
样本获取模块,用于获取不同骨折类型的颌面部骨折的CT数据样本;
样本处理模块,用于对所述CT数据样本进行标注处理和分类处理;
特征提取模块,用于从同一骨折类型的经过处理的CT数据样本训练并存储该骨折类型的特征;骨折识别模块,用于将所述图像特征与输入的待识别的颌面部CT数据进行比对,以确定输入的颌面部CT数据是否存在与所述图像特征对应的骨折类型;所述骨折类型按解剖区域进行划分;以及
序列提取模块,用于将输入的颌面部CT数据按解剖结构分解为不同解剖区域的图像块序列;将所述图像块序列中的每个图像块分别输入对应的训练好的骨折判别模型,以确定输入的颌面部CT数据是否存在对应的骨折类型;
其中,所述将输入的颌面部CT数据按解剖结构分解为不同解剖区域的图像块序列,包括:
将所述颌面部CT数据的CT图像序列输入到训练好的骨结构定位模型;
生成检测矩形框序列;
根据不同解剖区域的矩形框,定位每个解剖区域的断层范围;
分别从定位出的断层范围内的CT图像序列中逐层提取出每个解剖区域的图像块序列。
8.一种终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至6中任一项所述的基于CT数据库的颌面部骨折的识别方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011046297.7A CN111967540B (zh) | 2020-09-29 | 2020-09-29 | 基于ct数据库的颌面部骨折的识别方法、装置及终端设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011046297.7A CN111967540B (zh) | 2020-09-29 | 2020-09-29 | 基于ct数据库的颌面部骨折的识别方法、装置及终端设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111967540A CN111967540A (zh) | 2020-11-20 |
CN111967540B true CN111967540B (zh) | 2021-06-08 |
Family
ID=73386810
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011046297.7A Active CN111967540B (zh) | 2020-09-29 | 2020-09-29 | 基于ct数据库的颌面部骨折的识别方法、装置及终端设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111967540B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113298828B (zh) * | 2021-06-11 | 2023-09-22 | 上海交通大学医学院附属第九人民医院 | 一种基于卷积神经网络的颌骨自动分割方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108305248A (zh) * | 2018-01-17 | 2018-07-20 | 慧影医疗科技(北京)有限公司 | 一种骨折识别模型的构建方法及应用 |
CN108520519A (zh) * | 2018-04-11 | 2018-09-11 | 上海联影医疗科技有限公司 | 一种图像处理方法、装置及计算机可读存储介质 |
CN109119140A (zh) * | 2018-08-27 | 2019-01-01 | 北京大学口腔医学院 | 一种用于准确治疗陈旧性颧骨骨折的计算机辅助导航方法 |
CN110826557A (zh) * | 2019-10-25 | 2020-02-21 | 杭州依图医疗技术有限公司 | 一种骨折检出的方法及装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7664298B2 (en) * | 2003-03-25 | 2010-02-16 | Imaging Therapeutics, Inc. | Methods for the compensation of imaging technique in the processing of radiographic images |
CN110310723A (zh) * | 2018-03-20 | 2019-10-08 | 青岛海信医疗设备股份有限公司 | 骨骼图像处理方法、电子设备及存储介质 |
CN110458799A (zh) * | 2019-06-24 | 2019-11-15 | 上海皓桦科技股份有限公司 | 基于肋骨展开图的肋骨骨折自动检测方法 |
CN111667474A (zh) * | 2020-06-08 | 2020-09-15 | 杨天潼 | 骨折识别方法、装置、设备和计算机可读存储介质 |
-
2020
- 2020-09-29 CN CN202011046297.7A patent/CN111967540B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108305248A (zh) * | 2018-01-17 | 2018-07-20 | 慧影医疗科技(北京)有限公司 | 一种骨折识别模型的构建方法及应用 |
CN108520519A (zh) * | 2018-04-11 | 2018-09-11 | 上海联影医疗科技有限公司 | 一种图像处理方法、装置及计算机可读存储介质 |
CN109119140A (zh) * | 2018-08-27 | 2019-01-01 | 北京大学口腔医学院 | 一种用于准确治疗陈旧性颧骨骨折的计算机辅助导航方法 |
CN110826557A (zh) * | 2019-10-25 | 2020-02-21 | 杭州依图医疗技术有限公司 | 一种骨折检出的方法及装置 |
Non-Patent Citations (1)
Title |
---|
A Review of Intelligent Image Processing Method of Pulmonary CT Images;wenjun Tan et al;《2019 IEEE International Conference on Bioinformatics and Biomedicine (BIBM)》;20200206;1641-1648页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111967540A (zh) | 2020-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021082416A1 (zh) | 网络模型训练方法及装置、病灶区域确定方法及装置 | |
CN112150472A (zh) | 基于cbct的三维颌骨图像分割方法、装置及终端设备 | |
CN111967539B (zh) | 基于cbct数据库的颌面部骨折的识别方法、装置及终端设备 | |
JPWO2019146357A1 (ja) | 医療画像処理装置、方法及びプログラム並びに診断支援装置、方法及びプログラム | |
Lang et al. | Automatic localization of landmarks in craniomaxillofacial CBCT images using a local attention-based graph convolution network | |
CN112150473A (zh) | 基于ct的三维颌骨图像分割建模方法、装置及终端设备 | |
WO2022247007A1 (zh) | 医学图像分级方法、装置、电子设备及可读存储介质 | |
Chan et al. | Quasi-conformal statistical shape analysis of hippocampal surfaces for Alzheimer׳ s disease analysis | |
CN111967540B (zh) | 基于ct数据库的颌面部骨折的识别方法、装置及终端设备 | |
Nijiati et al. | Artificial intelligence assisting the early detection of active pulmonary tuberculosis from chest X-rays: A population-based study | |
CN113554640A (zh) | Ai模型的训练方法、使用方法、计算机设备及存储介质 | |
CN113948190A (zh) | X线头颅正位片头影测量标志点自动识别方法及设备 | |
Lang et al. | DLLNet: an attention-based deep learning method for dental landmark localization on high-resolution 3D digital dental models | |
WO2023029348A1 (zh) | 基于人工智能的图像实例标注方法及相关设备 | |
CN117237351A (zh) | 一种超声图像分析方法以及相关装置 | |
CN113989269B (zh) | 一种基于卷积神经网络多尺度特征融合的中医舌图像齿痕自动检测方法 | |
CN113450306B (zh) | 提供骨折检测工具的方法 | |
Tu et al. | Quantitative evaluation of local head malformations from 3 dimensional photography: application to craniosynostosis | |
El-Fegh et al. | Automated 2-D cephalometric analysis of X-ray by image registration approach based on least square approximator | |
Tanikawa et al. | Machine learning for facial recognition in orthodontics | |
Nalla et al. | Influence of Convolutional Neural Network Depth on the Efficacy of Automated Breast Cancer Screening Systems | |
Millan-Arias et al. | General Cephalometric Landmark Detection for Different Source of X-Ray Images | |
Yan et al. | Applying artificial intelligence algorithm in the design of a guide plate for mandibular angle ostectomy | |
Zhang et al. | The Implementation of Facial Symmetry Assessment before and after Orthognathic Surgery Using Transfer Learning | |
Sadr et al. | Deep learning for tooth identification and enumeration in panoramic radiographs |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |