CN114010215B - 用于医学影像辅助诊断支气管扩张的方法及装置 - Google Patents
用于医学影像辅助诊断支气管扩张的方法及装置 Download PDFInfo
- Publication number
- CN114010215B CN114010215B CN202111528793.0A CN202111528793A CN114010215B CN 114010215 B CN114010215 B CN 114010215B CN 202111528793 A CN202111528793 A CN 202111528793A CN 114010215 B CN114010215 B CN 114010215B
- Authority
- CN
- China
- Prior art keywords
- bronchus
- data
- node
- image
- dimensional
- 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
- 238000000034 method Methods 0.000 title claims abstract description 54
- 201000009267 bronchiectasis Diseases 0.000 title claims abstract description 51
- 238000003745 diagnosis Methods 0.000 title claims abstract description 31
- 210000000621 bronchi Anatomy 0.000 claims abstract description 225
- 210000002345 respiratory system Anatomy 0.000 claims abstract description 63
- 238000004364 calculation method Methods 0.000 claims abstract description 42
- 230000004807 localization Effects 0.000 claims abstract description 10
- 210000004072 lung Anatomy 0.000 claims description 34
- 238000003325 tomography Methods 0.000 claims description 34
- 230000002159 abnormal effect Effects 0.000 claims description 22
- 238000004422 calculation algorithm Methods 0.000 claims description 15
- 238000004140 cleaning Methods 0.000 claims description 14
- 238000003709 image segmentation Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 8
- 238000007667 floating Methods 0.000 claims description 8
- 239000000126 substance Substances 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 7
- 241000935061 Larrea Species 0.000 claims description 6
- 238000012216 screening Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 abstract description 16
- 230000002829 reductive effect Effects 0.000 abstract description 5
- 239000003814 drug Substances 0.000 abstract 1
- 238000002591 computed tomography Methods 0.000 description 52
- 210000000038 chest Anatomy 0.000 description 22
- 238000011156 evaluation Methods 0.000 description 11
- 238000004891 communication Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 239000006185 dispersion Substances 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000011002 quantification Methods 0.000 description 3
- 230000000241 respiratory effect Effects 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 210000000115 thoracic cavity Anatomy 0.000 description 3
- 210000001519 tissue Anatomy 0.000 description 3
- 230000006378 damage Effects 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000002372 labelling Methods 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 206010006440 Bronchial obstruction Diseases 0.000 description 1
- 241000711573 Coronaviridae Species 0.000 description 1
- 206010011224 Cough Diseases 0.000 description 1
- 206010016654 Fibrosis Diseases 0.000 description 1
- 208000000616 Hemoptysis Diseases 0.000 description 1
- 206010061218 Inflammation Diseases 0.000 description 1
- 201000005702 Pertussis Diseases 0.000 description 1
- 206010035664 Pneumonia Diseases 0.000 description 1
- 206010036790 Productive cough Diseases 0.000 description 1
- 206010057190 Respiratory tract infections Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 206010006451 bronchitis Diseases 0.000 description 1
- 208000013116 chronic cough Diseases 0.000 description 1
- 230000001684 chronic effect Effects 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 210000004177 elastic tissue Anatomy 0.000 description 1
- 230000004761 fibrosis Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 230000004054 inflammatory process Effects 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003097 mucus Anatomy 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000001717 pathogenic effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 210000003456 pulmonary alveoli Anatomy 0.000 description 1
- 208000008128 pulmonary tuberculosis Diseases 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 208000023504 respiratory system disease Diseases 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 208000024794 sputum Diseases 0.000 description 1
- 210000003802 sputum Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- 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/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- 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
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Public Health (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biophysics (AREA)
- Geometry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Pulmonology (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本申请涉及医学技术领域,特别涉及一种用于医学影像辅助诊断支气管扩张的方法及装置,其中,方法包括:采集人体呼吸系统的二维影像数据,并基于二维影像数据,重建支气管树空间三维数据;以支气管树空间三维数据为基础,识别每根支气管的中心线,并对每根支气管的中心线进行节点划分,计算每根支气管对应的所有节点参数,获得CT影像预设范围内所有支气管直径;对每根支气管进行命名,以对扩张支气管在肺叶中的位置进行靶向定位追踪。由此,实现了影像可及范围支气管直径的精确输出计算,大大提升了传统手动提取方法准确度,显著降低了时间成本。
Description
技术领域
本申请涉及医学技术领域,特别涉及一种用于医学影像辅助诊断支气管扩张的方法及装置。
背景技术
支气管扩张是由于支气管及其周围肺组织慢性化脓性炎症和纤维化,使支气管壁的肌肉和弹性组织破坏,导致支气管壁损伤引起的气道变形及不可逆性扩张。典型的症状有慢性咳嗽、咳大量脓痰和反复咯血,主要致病因素为支气管感染、阻塞和牵拉,是新型冠状病毒肺炎、肺结核、百日咳等呼吸道疾病的常见并发症。
准确的识别支气管扩张,是提高诊疗精确度的前提。现有医学研究发现,胸片评价支气管扩张的敏感性较低,在轻中度病例中,扩张支气管增厚的管壁表现的平行轨道征和环状影以及充满粘液的扩张气道表现的管样密度增高影在胸片上可以显示,但许多支气管扩张症在胸片上难以识别。而薄层胸部CT(Computed Tomography,电子计算机断层扫描)比胸片更敏感,该方法通过应用间隔为1.0mm或更小的高分辨率CT进行螺旋扫描,对于支气管扩张的评价有较高准确度,是目前支气管扩张特征性评价和鉴别评价的参考标准。
支气管扩张的核心影像表征是支气管直径增粗,判断标准主要有三种,分别为:(1)支气管-动脉比值增加:通过测量支气管横断面直径和肺内伴行动脉直径,计算支气管-动脉直径比值,正常支气管-动脉比值平均为0.65~0.7;(2)缺乏支气管逐渐变细的征象:缺乏支气管逐渐变细的表现,这一征象表现为气道从支气管分支点到远端至少50px直径保持不变;(3)支气管轮廓异常:支气管扩张可分为柱状、曲张型、囊状扩张,识别这些轮廓异常改变可以做出评价。由此可发现,准确的测量支气管管径是CT影像评价的关键。
然而,人体呼吸道解剖结构呈树状,由支气管(第1级)至肺泡约有24级分支,各级以二叉树的形式生长,支气管数量庞大,同时,每支支气管沿长度方向的断面直径并不相同,且与受试者个体、呼吸状态及病程发展等因素显著均相关,临床阅片的工作量极大。因此,传统支气管扩张的判别依赖于影像科医生依据经验判断扩张支气管,其量化则采用手动标注,通过在薄层CT上取若干断面手动测量支气管直径,工作量大且精确度欠佳。由于缺乏高精度的全自动支气管直径提取及计算工具,导致医生在支气管扩张程度分析及多次CT支气管扩张病程发展判断时,缺乏可靠的数据支撑。
申请内容
本申请提供一种用于医学影像辅助诊断支气管扩张的方法及装置,实现了影像可及范围支气管直径的自动精确输出计算,大大提升了传统手动提取方法准确度,显著降低了时间成本。
本申请第一方面实施例提供一种用于医学影像辅助诊断支气管扩张的方法,包括以下步骤:
采集人体呼吸系统的二维影像数据,并基于所述二维影像数据,重建支气管树空间三维数据;
以所述支气管树空间三维数据为基础,识别每根支气管的中心线,并对所述每根支气管的中心线进行节点划分,计算所述每根支气管对应的所有节点参数,判断并筛选所述每根支气管对应的有效节点,获得CT影像预设范围内所有支气管直径;
对所述每根支气管进行命名,以对扩张支气管在肺叶中的位置进行靶向定位追踪。
可选地,所述识别每根支气管的中心线,并对所述每根支气管的中心线进行节点划分,计算所述每根支气管对应的所有节点参数,包括:
识别所述每根支气管的中心线的同时,将所述中心线沿支气管生长方向划分成多个计算节点;
以每个计算节点为中心,采用最小二乘法获得所述每个计算节点对应的支气管断面最佳拟合圆直径,并输出各节点空间位置坐标及其对应的支气管断面最佳拟合圆直径的数据集;
查找所述数据集中空间位置坐标为空的无效节点,并依据二叉树算法结构遍历所述数据集,剔除非二叉树的冗余支气管数据,以完成所述数据集的二次数据清洗,并输出有效节点的空间位置坐标及其对应的支气管断面最佳拟合圆直径。
可选地,所述每个支气管有效节点的支气管断面最佳拟合圆直径数据集服从正态分布,其中,所述判断并筛选所述每根支气管对应的有效节点,获得CT影像预设范围内所有支气管直径,包括:
依据正态分布拉依达原则,对所述数据集的异常值进行剔除,计算所述每个支气管的有效节点最佳拟合圆直径的平均值和标准差;
判断所述每个支气管的有效节点的最佳拟合圆直径与所述平均值的偏差是否超过预设标准差,并在超过所述预设标准差时,判定为异常值节点;
剔除所述每个支气管有效节点中的异常值节点后,对所述每个支气管剩余所有节点对应的拟合圆直径进行算数平均,输出所述CT影像预设范围内所有支气管的直径标准数据集。
可选地,所述对所述每根支气管进行命名,以对扩张支气管在肺叶中的位置进行靶向定位追踪,包括:
依据所述每根支气管在肺叶中的位置及级数,对所述每根支气管进行命名,得到所述每根支气管命名序列;
根据所述每根支气管命名序列对所述每根支气管在肺叶中的位置进行定位,以对所述CT影像中的扩张支气管直径进行靶向定位追踪。
可选地,所述采集人体呼吸系统的二维影像数据,并基于所述二维影像数据,重建支气管树空间三维数据,包括:
围绕人体胸部进行断层扫描,获得层厚小于预设阈值的薄层人体胸部二维断层扫描影像;并基于所述人体胸部二维断层扫描影像,获得多个由黑到白不同灰度按矩阵排列的像素点,并获得上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像;
根据所述上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像获取人体呼吸道支气管树的所述二维影像数据;
根据所述二维影像数据分别对所述冠状位、水平位、矢状位断面影像进行逐层细化处理,得到去除背景噪声和像素中的漂浮物质后的每张薄层CT支气管完整精细的图像分割结果;并根据预设呼吸道阈值确定等值面,并利用线性差值逼近所述等值面,由所述图像分割结果生成网格模型数据,得到所述支气管树空间三维数据。
本申请第二方面实施例提供一种用于医学影像辅助诊断支气管扩张的装置,包括:
重建模块,用于采集人体呼吸系统的二维影像数据,并基于所述二维影像数据,重建支气管树空间三维数据
获取模块,用于以所述支气管树空间三维数据为基础,识别每根支气管的中心线,并对所述每根支气管的中心线进行节点划分,计算所述每根支气管对应的所有节点参数,判断并筛选所述每根支气管对应的有效节点,获得CT影像预设范围内所有支气管直径;
定位模块,用于对所述每根支气管进行命名,以根据所述每根支气管命名序列对扩张支气管在肺叶中的位置进行靶向定位追踪。
可选地,所述获取模块,具体用于:
识别所述每根支气管的中心线的同时,将所述中心线沿支气管生长方向划分成多个计算节点;
以每个计算节点为中心,采用最小二乘法获得所述每个计算节点对应的支气管断面最佳拟合圆直径,并输出各计算节点的空间位置坐标及其对应的支气管断面最佳拟合圆直径的数据集;
查找所述数据集中空间位置坐标为空的无效节点,并依据二叉树算法结构遍历所述数据集,剔除非二叉树的冗余支气管数据,以完成所述数据集的二次数据清洗,并输出有效节点的空间位置坐标及其对应的支气管断面最佳拟合圆直径。
可选地,所述每个支气管有效节点的支气管断面最佳拟合圆直径数据集服从正态分布,所述获取模块,具体用于:
依据正态分布拉依达原则,对所述数据集的异常值进行剔除,计算所述每个支气管的有效节点最佳拟合圆直径的平均值和标准差;
判断所述每个支气管的有效节点的最佳拟合圆直径与所述平均值的偏差是否超过预设标准差,并在超过所述预设标准差时,判定为异常值节点;
剔除所述每个支气管有效节点中的异常值节点后,对所述每个支气管剩余所有节点对应的拟合圆直径进行算数平均,输出所述CT影像预设范围内所有支气管的直径标准数据集。
可选地,所述定位模块,具体用于:
依据所述每根支气管在肺叶中的位置及级数,对所述每根支气管进行命名,得到所述每根支气管命名序列;
根据所述每根支气管命名序列对所述每根支气管在肺叶中的位置进行定位,以对所述CT影像中的扩张支气管直径进行靶向定位追踪。
可选地,所述重建模块,具体用于:
围绕人体胸部进行断层扫描,获得层厚小于预设阈值的薄层人体胸部二维断层扫描影像,并基于所述人体胸部二维断层扫描影像,获得多个由黑到白不同灰度按矩阵排列的像素点,并获得上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像;
根据所述上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像获取人体呼吸道支气管树的所述二维影像数据;
根据所述二维影像数据分别对所述冠状位、水平位、矢状位断面影像进行逐层细化处理,得到去除背景噪声和像素中的漂浮物质后的每张薄层CT支气管完整精细的图像分割结果;并根据预设呼吸道阈值确定等值面,并利用线性差值逼近所述等值面,由所述图像分割结果生成网格模型数据,得到所述支气管树空间三维数据。
本申请第三方面实施例提供一种电子设备,包括:存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述程序,以实现如上述实施例所述的用于医学影像辅助诊断支气管扩张的方法。
本申请第四方面实施例提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行,以用于实现上述的用于医学影像辅助诊断支气管扩张的方法。
由此,通过对CT影像获得的大量呼吸道二维断层扫描图像进行形态提取及空间离散,利用数据抓取及清洗算法,实现了影像可及范围支气管直径的自动精确输出计算,大大提升了传统手动提取方法准确度,显著降低了时间成本,在支气管扩张的辅助评价方面具有极高的应用前景。本申请附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本申请的实践了解到。
附图说明
本申请上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本申请实施例提供的一种用于医学影像辅助诊断支气管扩张的方法的流程图;
图2为根据本申请一个实施例的目标个体胸腔区域CT断层扫描示例图;
图3为根据本申请一个实施例的薄层CT图像的支气管阈值分割示意图;
图4为根据本申请一个实施例的基于Marching Cubes算法生成的支气管树空间三维模型的示例图;
图5为根据本申请一个实施例的支气管中心线识别及节点划分示意图;
图6为根据本申请一个实施例的各节点支气管断面最佳拟合圆直径计算的示例图;
图7为根据本申请一个实施例的中心气道部分节点空间位置及断面最佳拟合圆直径计算示意图;
图8为根据本申请一个实施例的扩张支气管直径输出及靶向定位示意图;
图9为根据本申请实施例的用于医学影像辅助诊断支气管扩张的装置的示例图;
图10为根据本申请实施例的电子设备的示例图。
具体实施方式
下面详细描述本申请的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本申请,而不能理解为对本申请的限制。
下面参考附图描述本申请实施例的用于医学影像辅助诊断支气管扩张的方法及装置。
在介绍本申请实施例的用于医学影像辅助诊断支气管扩张的方法之前,先简单介绍下相关技术中的支气管直径的量化方法。
相关技术中,支气管直径的量化采用离体组织切片或人工标记的方法,对于患者CT影像学评价,多用人工标记方法。该方法的核心是影像科医生根据影像判断标准,通过审阅患者胸部薄层CT扫描影像,找到扩张支气管位置,通过CT机配套的后处理软件辅助测量支气管直径,判断是否存在支气管扩张。该方法大大依赖于医生的阅片经验,由于支气管数量极多,人工提取的工作量很大,阅片中易存在漏判、误判的情况,准确度亟待提升。
正是基于上述问题,本申请提供了一种用于医学影像辅助诊断支气管扩张的方法,在该方法中,通过对CT影像获得的大量呼吸道二维断层扫描图像进行形态提取及空间离散,利用数据抓取及清洗算法,实现了影像可及范围支气管直径的自动精确输出计算,大大提升了传统手动提取方法准确度,显著降低了时间成本,在支气管扩张的辅助评价方面具有极高的应用前景。
具体而言,图1为本申请实施例所提供的一种用于医学影像辅助诊断支气管扩张的方法的流程示意图。
如图1所示,该用于医学影像辅助诊断支气管扩张的方法包括以下步骤:
在步骤S101中,采集人体呼吸系统的二维影像数据,并基于二维影像数据,重建支气管树空间三维数据。
可选地,采集人体呼吸系统的二维影像数据,并基于二维影像数据,重建支气管树空间三维数据,包括:围绕人体胸部进行断层扫描,获得层厚小于预设阈值的薄层人体胸部二维断层扫描影像;并基于人体胸部二维断层扫描影像,获得多个由黑到白不同灰度按矩阵排列的像素点,并获得上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像;根据上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像获取人体呼吸道支气管树的二维影像数据;根据二维影像数据分别对冠状位、水平位、矢状位断面影像进行逐层细化处理,得到去除背景噪声和像素中的漂浮物质后的每张薄层CT支气管完整精细的图像分割结果;并根据预设呼吸道阈值确定等值面,并利用线性差值逼近等值面,由图像分割结果生成网格模型数据,得到支气管树空间三维数据。
具体而言,本申请实施例可以通过电子计算机断层扫描技术,利用精确准直的X线束,与灵敏度极高的探测器一同围绕人体胸部进行一个接一个的断面扫描,获得层厚小于1.5mm的薄层人体胸部二维断层扫描影像。依据不同组织对X线的吸收率与透过率的不同,获得一定数目由黑到白不同灰度按矩阵排列的像素点,然后将测量所获取的数据输入电子计算机,对数据进行处理后,显示上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像,快速、便捷地获取人体呼吸道支气管树原始二维影像数据。
进一步地,本申请实施例可以定义呼吸道阈值范围,对冠状位、水平位、矢状位三个视窗断层影像数据人工逐层细化处理,去除背景噪声和像素中的漂浮物质,实现每张薄层CT支气管完整精细的图像分割。本申请实施例可以采用Marching Cubes(MC,移动立方体)算法,假设原始数据为三维离散数据场,通过已定义的阈值确定等值面,利用线性差值逼近等值面,将支气管数据生成网格模型数据,输出支气管树空间三维模型。
在步骤S102中,以支气管树空间三维数据为基础,识别每根支气管的中心线,并对每根支气管的中心线进行节点划分,计算每根支气管对应的所有节点参数,判断并筛选每根支气管对应的有效节点,获得CT影像预设范围内所有支气管直径。
具体地,本申请实施例可以以支气管树空间三维数据为基础,识别每根支气管的中心线及计算对应参数,并且由每根支气管的中心线得到CT影像,识别预设范围内所有支气管直径,同时对每根支气管进行命名,以对扩张支气管进行靶向定位追踪。
可选地,在一些实施例中,识别每根支气管的中心线,并对每根支气管的中心线进行节点划分,计算每根支气管对应的所有节点参数,包括:识别每根支气管的中心线的同时,将中心线沿支气管生长方向划分成多个计算节点;以每个计算节点为中心,采用最小二乘法获得每个计算节点对应的支气管断面最佳拟合圆直径,并输出各节点空间位置坐标及其对应的支气管断面最佳拟合圆直径的数据集;利用计算机编程自动查找数据集中空间位置坐标为空的无效节点,并依据二叉树算法结构遍历数据集,剔除非二叉树的冗余支气管数据,以完成数据集的二次数据清洗,并输出有效节点的空间位置坐标及其对应的支气管断面最佳拟合圆直径。
可选地,在一些实施例中,每个支气管有效节点的支气管断面最佳拟合圆直径数据集服从正态分布,其中,判断并筛选每根支气管对应的有效节点,获得CT影像预设范围内所有支气管直径,包括:依据正态分布拉依达原则,利用计算机编程对数据集的异常值进行剔除,计算每个支气管的有效节点最佳拟合圆直径的平均值和标准差;判断每个支气管的有效节点的最佳拟合圆直径与平均值的偏差是否超过预设标准差(如3倍标准差),并在超过预设标准差时,判定为异常值节点;剔除每个支气管有效节点中的异常值节点后,对每个支气管剩余所有节点对应的拟合圆直径进行算数平均,输出CT影像预设范围内所有支气管的直径标准数据集。由此,完成CT影像预设范围内所有支气管直径的全自动提取及高精度计算,并输出支气管直径标准数据集,可辅助影像科医生对二维CT影像中的扩张支气管直径进行快速量化。其中,预设标准差可以为
具体地,每支支气管沿长度方向截面直径不相等,为了让支气管的直径计算更加准确,以输出的支气管树空间三维模型为基础,识别每一根支气管中心线,将中心线沿支气管生长方向划分为若干个计算节点,相邻节点间距小于1.3mm。以每一个节点为中心,采用最小二乘法获得每个节点对应的支气管断面最佳拟合圆直径,输出各节点空间位置坐标及其对应的支气管断面最佳拟合圆直径的.txt或.csv数据集。
进一步地,基于计算机高级编程语言,建立支气管直径自动计算及存储介质。对.txt或.csv后缀的原始数据集进行初步判断、分项、数据清洗,依据正态分布拉依达原则,对异常值进行剔除,计算每个支气管有效节点最佳拟合圆直径的平均值和标准差;此后,判断每个支气管有效节点的最佳拟合圆直径与平均值的偏差是否超过预设标准差,并在超过预设标准差时,判定为异常值节点;最后,剔除每个支气管有效节点中的异常值节点后,对每个支气管剩余所有节点对应的拟合圆直径进行算数平均,输出CT影像预设范围内所有支气管的直径标准数据集。
在步骤S103中,对每根支气管进行命名,以对扩张支气管在肺叶中的位置进行靶向定位追踪。
可选地,对每根支气管进行命名,以对扩张支气管在肺叶中的位置进行靶向定位追踪,包括:依据每根支气管在肺叶中的位置及级数,对每根支气管进行命名,得到每根支气管命名序列;根据每根支气管命名序列对每根支气管在肺叶中的位置进行定位,以对CT影像中的扩张支气管直径进行靶向定位追踪。由此,可辅助影像科医生对二维CT影像中的扩张支气管直径进行靶向追踪,完成多次CT扫描支气管扩张病程发展分析
由此,通过对人体胸部二维CT断层扫描影像获得的大量无序节点数据集进行算法编程,实现了影像提取范围内所有支气管直径的自动自动输出计算,适用于任意单体呼吸道支气管直径识别及分析,显著提升了传统手动标注的准确性,缩短了人工提取工作量,为支气管扩张的影像学量化及靶向追踪提供了新的技术方法。
为使得本领域技术人员进一步了解本申请实施例的用于医学影像辅助诊断支气管扩张的方法,下面基于具体实施例进行详细阐述。
具体地,本申请基于目标个体CT断层扫描影像,还原前7级呼吸道形态学结构,计算全部支气管直径并输出,辅助影像科医生完成支气管扩张部位的定位及直径量化工作,主要包括以下步骤:
(1)人体呼吸系统二维影像数据采集模块:
采用电子计算机断层扫描(CT)技术对目标个体进行胸腔CT断层扫描,快速、便捷地获得目标个体高精度呼吸系统的冠状位、矢状位、水平位断面影像,如图2(a)-图2(c)所示,图2(a)为目标个体矢状位胸腔区域CT断层扫描图,图2(b)为目标个体冠状位胸腔区域CT断层扫描图,图2(c)为目标个体水平位胸腔区域CT断层扫描图,准确再现不同性别、年龄、健康状况等因素下目标人群呼吸系统解剖学结构差异。其中,黑影表示低吸收区,即低密度区,如含气体多的呼吸道;白影表示高吸收区,即高密度区,如骨骼。
(2)支气管树空间三维数据重建模块:
将步骤(1)中人体呼吸道二维影像数据导入支气管树空间三维数据重建模块,定义呼吸道阈值范围,对冠状位、水平位、矢状位三个视窗断层影像数据人工逐层细化处理,去除背景噪声和像素中的漂浮物质,实现每张薄层CT支气管完整精细的图像分割,如图3所示,其中,图3为薄层CT图像的支气管阈值分割示意图,图3(a)为中心气道(0级);图3(b)为左主和右主支气管(1级);图3(c)2级支气管。本申请实施例采用Marching Cubes(移动立方体)算法,假设原始数据为三维离散数据场,通过已定义的阈值确定等值面,利用线性差值逼近等值面,将支气管数据生成网格模型数据,输出支气管树空间三维模型,基于MarchingCubes算法生成的支气管树空间三维模型可以如图4所示。
(3)支气管树中心线提取及参数计算模块:
以步骤(2)输出的支气管树空间三维模型为基础,识别每一根支气管中心线,将中心线沿支气管生长方向划分为若干个计算节点,相邻节点间距为~1.3mm,如图5所示,图5为支气管中心线识别及节点划分示意图。以每一个节点为中心,采用最小二乘法获得每个节点对应的支气管断面最佳拟合圆直径,如图6所示。
由此,即可输出各节点空间位置坐标及其对应的支气管断面最佳拟合圆直径的.txt或.csv数据集。
(4)支气管直径自动输出计算模块
步骤(3)中输出的.txt或.csv数据集,数据量依据个体呼吸道形态结构差异有所不同,包括了千/万级节点空间位置坐标及其对应的支气管断面最佳拟合圆直径原始数据,如图7所示。基于计算机高级编程语言,利用数据抓取及清洗算法,实现支气管直径自动计算、输出及存储。对.txt或.csv后缀的原始数据集进行初步判断、分项、数据清洗,寻找无效节点并剔除冗余数据信息,此后,进行二次数据清洗,将同一支气管中心线所有节点对应的拟合圆直径进行算数平均,输出CT影像可识别范围内所有支气管直径标准数据集。
也即是说,本申请实施例可以首先进行每根支气管的中心线节点划分,然后计算每个节点的拟合圆直径,从而计算每根支气管所有拟合圆直径的算数平均值-获得支气管直径,计算公式为:
其中,n为支气管根数,Dn为第n根支气管直径,i为第i个节点,Di为第i个节点断面最佳拟合圆直径,m为每根支气管的节点数。
经计算,该算例共识别支气管数量276根,对每一根支气管均进行命名,命名结果可以如表1所示,表1为部分支气管直径自动输出计算表(部分)。
表1
由此,如图8所示,白色虚线为扩张呼吸道-复查曲线,灰色曲线为扩张呼吸道-感染曲线,基于影像科医生对支气管扩张位置的预判,实现了扩张支气管的靶向定位追踪,以及扩张直径的定量化输出,为支气管扩张的影像评价提供了技术支持。
根据本申请实施例提出的用于医学影像辅助诊断支气管扩张的方法,通过对CT影像获得的大量呼吸道二维断层扫描图像进行形态提取及空间离散,利用数据抓取及清洗算法,实现了影像可及范围支气管直径的自动精确输出计算,大大提升了传统手动提取方法准确度,显著降低了时间成本,在支气管扩张的辅助评价方面具有极高的应用前景。
其次参照附图描述根据本申请实施例提出的用于医学影像辅助诊断支气管扩张的装置。
图9是本申请实施例的用于医学影像辅助诊断支气管扩张的装置的方框示意图。
如图9所示,该用于医学影像辅助诊断支气管扩张的装置10包括:重建模块100、获取模块200和定位模块300。
其中,重建模块100用于采集人体呼吸系统的二维影像数据,并基于二维影像数据,重建支气管树空间三维数据
获取模块200用于以支气管树空间三维数据为基础,识别每根支气管的中心线,并对每根支气管的中心线进行节点划分,计算每根支气管对应的所有节点参数,判断并筛选每根支气管对应的有效节点,获得CT影像预设范围内所有支气管直径;
定位模块300用于对每根支气管进行命名,以根据每根支气管命名序列对扩张支气管在肺叶中的位置进行靶向定位追踪。
可选地,获取模块200具体用于:
识别每根支气管的中心线的同时,将中心线沿支气管生长方向划分成多个计算节点;
以每个计算节点为中心,采用最小二乘法获得每个计算节点对应的支气管断面最佳拟合圆直径,并输出各计算节点的空间位置坐标及其对应的支气管断面最佳拟合圆直径的数据集;
查找数据集中空间位置坐标为空的无效节点,并依据二叉树算法结构遍历数据集,剔除非二叉树的冗余支气管数据,以完成数据集的二次数据清洗,并输出有效节点的空间位置坐标及其对应的支气管断面最佳拟合圆直径。
可选地,每个支气管有效节点的支气管断面最佳拟合圆直径数据集服从正态分布,获取模块200具体用于:
依据正态分布拉依达原则,对数据集的异常值进行剔除,计算每个支气管的有效节点最佳拟合圆直径的平均值和标准差;
判断每个支气管的有效节点的最佳拟合圆直径与平均值的偏差是否超过预设标准差,并在超过预设标准差时,判定为异常值节点;
剔除每个支气管有效节点中的异常值节点后,对每个支气管剩余所有节点对应的拟合圆直径进行算数平均,输出CT影像预设范围内所有支气管的直径标准数据集。
可选地,定位模块300具体用于:
依据每根支气管在肺叶中的位置及级数,对每根支气管进行命名,得到每根支气管命名序列;
根据每根支气管命名序列对每根支气管在肺叶中的位置进行定位,以对CT影像中的扩张支气管直径进行靶向定位追踪。
可选地,重建模块100具体用于:
围绕人体胸部进行断层扫描,获得层厚小于预设阈值的薄层人体胸部二维断层扫描影像,并基于人体胸部二维断层扫描影像,获得多个由黑到白不同灰度按矩阵排列的像素点,并获得上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像;
根据上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像获取人体呼吸道支气管树的二维影像数据;
根据二维影像数据分别对冠状位、水平位、矢状位断面影像进行逐层细化处理,得到去除背景噪声和像素中的漂浮物质后的每张薄层CT支气管完整精细的图像分割结果;并根据预设呼吸道阈值确定等值面,并利用线性差值逼近等值面,由图像分割结果生成网格模型数据,得到支气管树空间三维数据。
需要说明的是,前述对用于医学影像辅助诊断支气管扩张的方法实施例的解释说明也适用于该实施例的用于医学影像辅助诊断支气管扩张的装置,此处不再赘述。
根据本申请实施例提出的用于医学影像辅助诊断支气管扩张的装置,通过对CT影像获得的大量呼吸道二维断层扫描图像进行形态提取及空间离散,利用数据抓取及清洗算法,实现了影像可及范围支气管直径的自动精确输出计算,大大提升了传统手动提取方法准确度,显著降低了时间成本,在支气管扩张的辅助评价方面具有极高的应用前景。
图10为本申请实施例提供的电子设备的结构示意图。该电子设备可以包括:
存储器1001、处理器1002及存储在存储器1001上并可在处理器1002上运行的计算机程序。
处理器1002执行程序时实现上述实施例中提供的用于医学影像辅助诊断支气管扩张的方法。
进一步地,电子设备还包括:
通信接口1003,用于存储器1001和处理器1002之间的通信。
存储器1001,用于存放可在处理器1002上运行的计算机程序。
存储器1001可能包含高速RAM存储器,也可能还包括非易失性存储器(non-volatile memory),例如至少一个磁盘存储器。
如果存储器1001、处理器1002和通信接口1003独立实现,则通信接口1003、存储器1001和处理器1002可以通过总线相互连接并完成相互间的通信。总线可以是工业标准体系结构(Industry Standard Architecture,简称为ISA)总线、外部设备互连(PeripheralComponent,简称为PCI)总线或扩展工业标准体系结构(Extended Industry StandardArchitecture,简称为EISA)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示,图10中仅用一条粗线表示,但并不表示仅有一根总线或一种类型的总线。
可选地,在具体实现上,如果存储器1001、处理器1002及通信接口1003,集成在一块芯片上实现,则存储器1001、处理器1002及通信接口1003可以通过内部接口完成相互间的通信。
处理器1002可能是一个中央处理器(Central Processing Unit,简称为CPU),或者是特定集成电路(Application Specific Integrated Circuit,简称为ASIC),或者是被配置成实施本申请实施例的一个或多个集成电路。
本实施例还提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如上的用于医学影像辅助诊断支气管扩张的方法。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本申请的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或N个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本申请的描述中,“N个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为,表示包括一个或更N个用于实现定制逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分,并且本申请的优选实施方式的范围包括另外的实现,其中可以不按所示出或讨论的顺序,包括根据所涉及的功能按基本同时的方式或按相反的顺序,来执行功能,这应被本申请的实施例所属技术领域的技术人员所理解。
在流程图中表示或在此以其他方式描述的逻辑和/或步骤,例如,可以被认为是用于实现逻辑功能的可执行指令的定序列表,可以具体实现在任何计算机可读介质中,以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用,或结合这些指令执行系统、装置或设备而使用。就本说明书而言,"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下:具有一个或N个布线的电连接部(电子装置),便携式计算机盘盒(磁装置),随机存取存储器(RAM),只读存储器(ROM),可擦除可编辑只读存储器(EPROM或闪速存储器),光纤装置,以及便携式光盘只读存储器(CDROM)。另外,计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质,因为可以例如通过对纸或其他介质进行光学扫描,接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序,然后将其存储在计算机存储器中。
应当理解,本申请的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中,N个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。如,如果用硬件来实现和在另一实施方式中一样,可用本领域公知的下列技术中的任一项或他们的组合来实现:具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路,具有合适的组合逻辑门电路的专用集成电路,可编程门阵列(PGA),现场可编程门阵列(FPGA)等。
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,该程序在执行时,包括方法实施例的步骤之一或其组合。
此外,在本申请各个实施例中的各功能单元可以集成在一个处理模块中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。
上述提到的存储介质可以是只读存储器,磁盘或光盘等。尽管上面已经示出和描述了本申请的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本申请的限制,本领域的普通技术人员在本申请的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (8)
1.一种用于医学影像辅助诊断支气管扩张的方法,其特征在于,包括以下步骤:
采集人体呼吸系统的二维影像数据,并基于所述二维影像数据,重建支气管树空间三维数据;
以所述支气管树空间三维数据为基础,识别每根支气管的中心线,并对所述每根支气管的中心线进行节点划分,计算所述每根支气管对应的所有节点参数,判断并筛选所述每根支气管对应的有效节点,获得CT影像预设范围内所有支气管直径;
对所述每根支气管进行命名,以对扩张支气管在肺叶中的位置进行靶向定位追踪;
其中,所述识别每根支气管的中心线,并对所述每根支气管的中心线进行节点划分,计算所述每根支气管对应的所有节点参数,包括:识别所述每根支气管的中心线的同时,将所述中心线沿支气管生长方向划分成多个计算节点;以每个计算节点为中心,采用最小二乘法获得所述每个计算节点对应的支气管断面最佳拟合圆直径,并输出各计算节点的空间位置坐标及其对应的支气管断面最佳拟合圆直径的数据集;查找所述数据集中空间位置坐标为空的无效节点,并依据二叉树算法结构遍历所述数据集,剔除非二叉树的冗余支气管数据,以完成所述数据集的二次数据清洗,并输出有效节点的空间位置坐标及其对应的支气管断面最佳拟合圆直径;
所述对所述每根支气管进行命名,以对扩张支气管在肺叶中的位置进行靶向定位追踪,包括:依据所述每根支气管在肺叶中的位置及级数,对所述每根支气管进行命名,得到所述每根支气管命名序列;根据所述每根支气管命名序列对所述每根支气管在肺叶中的位置进行定位,以对所述CT影像中的扩张支气管直径进行靶向定位追踪。
2.根据权利要求1所述的方法,其特征在于,所述每根 支气管有效节点的支气管断面最佳拟合圆直径数据集服从正态分布,其中,所述判断并筛选所述每根支气管对应的有效节点,获得CT影像预设范围内所有支气管直径,包括:
依据正态分布拉依达原则,对所述数据集的异常值进行剔除,计算所述每个支气管的有效节点最佳拟合圆直径的平均值和标准差;
判断所述每个支气管的有效节点的最佳拟合圆直径与所述平均值的偏差是否超过预设标准差,并在超过所述预设标准差时,判定为异常值节点;
剔除所述每个支气管有效节点中的异常值节点后,对所述每个支气管剩余所有节点对应的拟合圆直径进行算数平均,输出所述CT影像预设范围内所有支气管的直径标准数据集。
3.根据权利要求1所述的方法,其特征在于,所述采集人体呼吸系统的二维影像数据,并基于所述二维影像数据,重建支气管树空间三维数据,包括:
围绕人体胸部进行断层扫描,获得层厚小于预设阈值的薄层人体胸部二维断层扫描影像,并基于所述人体胸部二维断层扫描影像,获得多个由黑到白不同灰度按矩阵排列的像素点,并获得上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像;
根据所述上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像获取人体呼吸道支气管树的所述二维影像数据;
根据所述二维影像数据分别对所述冠状位、水平位、矢状位断面影像进行逐层细化处理,得到去除背景噪声和像素中的漂浮物质后的每张薄层CT支气管完整精细的图像分割结果,并根据预设呼吸道阈值确定等值面,并利用线性差值逼近所述等值面,由所述图像分割结果生成网格模型数据,得到所述支气管树空间三维数据。
4.一种用于医学影像辅助诊断支气管扩张的装置,其特征在于,包括:
重建模块,用于采集人体呼吸系统的二维影像数据,并基于所述二维影像数据,重建支气管树空间三维数据;
获取模块,用于以所述支气管树空间三维数据为基础,识别每根支气管的中心线,并对所述每根支气管的中心线进行节点划分,计算所述每根支气管对应的所有节点参数,判断并筛选所述每根支气管对应的有效节点,获得CT影像预设范围内所有支气管直径;
定位模块,用于对所述每根支气管进行命名,以根据所述每根支气管命名序列对扩张支气管在肺叶中的位置进行靶向定位追踪;
其中,所述获取模块,具体用于:识别所述每根支气管的中心线的同时,将所述中心线沿支气管生长方向划分成多个计算节点;以每个计算节点为中心,采用最小二乘法获得所述每个计算节点对应的支气管断面最佳拟合圆直径,并输出各计算节点的空间位置坐标及其对应的支气管断面最佳拟合圆直径的数据集;查找所述数据集中空间位置坐标为空的无效节点,并依据二叉树算法结构遍历所述数据集,剔除非二叉树的冗余支气管数据,以完成所述数据集的二次数据清洗,并输出有效节点的空间位置坐标及其对应的支气管断面最佳拟合圆直径;
所述定位模块,具体用于:依据所述每根支气管在肺叶中的位置及级数,对所述每根支气管进行命名,得到所述每根支气管命名序列;根据所述每根支气管命名序列对所述每根支气管在肺叶中的位置进行定位,以对所述CT影像中的扩张支气管直径进行靶向定位追踪。
5.根据权利要求4所述的装置,其特征在于,所述每根 支气管有效节点的支气管断面最佳拟合圆直径数据集服从正态分布,所述获取模块,具体用于:
依据正态分布拉依达原则,对所述数据集的异常值进行剔除,计算所述每个支气管的有效节点最佳拟合圆直径的平均值和标准差;
判断所述每个支气管的有效节点的最佳拟合圆直径与所述平均值的偏差是否超过预设标准差,并在超过所述预设标准差时,判定为异常值节点;
剔除所述每个支气管有效节点中的异常值节点后,对所述每个支气管剩余所有节点对应的拟合圆直径进行算数平均,输出所述CT影像预设范围内所有支气管的直径标准数据集。
6.根据权利要求4所述的装置,其特征在于,所述重建模块,具体用于:
围绕人体胸部进行断层扫描,获得层厚小于预设阈值的薄层人体胸部二维断层扫描影像,并基于所述人体胸部二维断层扫描影像,获得多个由黑到白不同灰度按矩阵排列的像素点,并获得上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像;
根据所述上呼吸道、下呼吸道及肺部的冠状位、水平位、矢状位断面影像获取人体呼吸道支气管树的所述二维影像数据;
根据所述二维影像数据分别对所述冠状位、水平位、矢状位断面影像进行逐层细化处理,得到去除背景噪声和像素中的漂浮物质后的每张薄层CT支气管完整精细的图像分割结果;并根据预设呼吸道阈值确定等值面,并利用线性差值逼近所述等值面,由所述图像分割结果生成网格模型数据,得到所述支气管树空间三维数据。
7.一种电子设备,其特征在于,包括:存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述程序,以实现如权利要求1-3任一项所述的用于医学影像辅助诊断支气管扩张的方法。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行,以用于实现如权利要求1-3任一项所述的用于医学影像辅助诊断支气管扩张的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111528793.0A CN114010215B (zh) | 2021-12-14 | 2021-12-14 | 用于医学影像辅助诊断支气管扩张的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111528793.0A CN114010215B (zh) | 2021-12-14 | 2021-12-14 | 用于医学影像辅助诊断支气管扩张的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114010215A CN114010215A (zh) | 2022-02-08 |
CN114010215B true CN114010215B (zh) | 2022-08-05 |
Family
ID=80068676
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111528793.0A Active CN114010215B (zh) | 2021-12-14 | 2021-12-14 | 用于医学影像辅助诊断支气管扩张的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114010215B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115661149B (zh) * | 2022-12-26 | 2023-03-10 | 天津医科大学总医院 | 一种基于肺部组织数据的肺部影像处理系统 |
CN116612891B (zh) * | 2023-07-14 | 2023-09-29 | 营动智能技术(山东)有限公司 | 一种慢性病患者数据处理系统 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7346203B2 (en) * | 2003-11-19 | 2008-03-18 | General Electric Company | Methods and apparatus for processing image data to aid in detecting disease |
WO2006000953A1 (en) * | 2004-06-22 | 2006-01-05 | Koninklijke Philips Electronics N.V. | Displaying a tracheobronchial tree |
CN112006772B (zh) * | 2020-08-04 | 2021-04-02 | 清华大学 | 完整人体体外呼吸道建立方法及系统 |
CN113628346B (zh) * | 2021-06-30 | 2023-08-08 | 上海市胸科医院 | 一种自由浏览vb的方法以及进行标记的方法和系统 |
-
2021
- 2021-12-14 CN CN202111528793.0A patent/CN114010215B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114010215A (zh) | 2022-02-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110969622B (zh) | 用于辅助肺炎诊断的图像处理方法和系统 | |
Yun et al. | Improvement of fully automated airway segmentation on volumetric computed tomographic images using a 2.5 dimensional convolutional neural net | |
Kirişli et al. | Standardized evaluation framework for evaluating coronary artery stenosis detection, stenosis quantification and lumen segmentation algorithms in computed tomography angiography | |
CN114010215B (zh) | 用于医学影像辅助诊断支气管扩张的方法及装置 | |
Paniagua et al. | Clinical application of SPHARM-PDM to quantify temporomandibular joint osteoarthritis | |
CN1663530B (zh) | 用于处理图像数据以帮助检测疾病的方法和装置 | |
US7970196B2 (en) | Automatic detection and quantification of plaque in the coronary arteries of subjects from CT scans | |
JP2008503303A (ja) | 気管気管支樹を表示する方法及び装置 | |
Pu et al. | CT based computerized identification and analysis of human airways: a review | |
Thein et al. | An image preprocessing method for kidney stone segmentation in CT scan images | |
US20080279435A1 (en) | Automated Calcium Scoring of the Aorta | |
EP3971830B1 (en) | Pneumonia sign segmentation method and apparatus, medium and electronic device | |
CN112862833A (zh) | 血管分段方法、电子装置和存储介质 | |
CN115222937A (zh) | 一种脊柱侧弯检测方法及装置 | |
WO2022164374A1 (en) | Automated measurement of morphometric and geometric parameters of large vessels in computed tomography pulmonary angiography | |
Rebouças Filho et al. | 3D segmentation and visualization of lung and its structures using CT images of the thorax | |
Peng et al. | Pulmonary lobe segmentation in CT images based on lung anatomy knowledge | |
Nie et al. | OII-DS: A benchmark Oral Implant Image Dataset for object detection and image classification evaluation | |
CN114203297B (zh) | 用于医学影像的呼吸道疾病随诊辅助方法及装置 | |
Wang et al. | Accuracy and reliability analysis of a machine learning based segmentation tool for intertrochanteric femoral fracture CT | |
Perchet et al. | Advanced navigation tools for virtual bronchoscopy | |
Hu et al. | Slice grouping and aggregation network for auxiliary diagnosis of rib fractures | |
CN111415341A (zh) | 肺炎阶段的评估方法、装置、介质及电子设备 | |
Wang et al. | Thoracic aorta diameter calculation by artificial intelligence can predict the degree of arterial stiffness | |
CN111358484B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |