CN107154038B - 一种肋骨可视化的肋骨骨折辅助诊断方法 - Google Patents
一种肋骨可视化的肋骨骨折辅助诊断方法 Download PDFInfo
- Publication number
- CN107154038B CN107154038B CN201710273035.6A CN201710273035A CN107154038B CN 107154038 B CN107154038 B CN 107154038B CN 201710273035 A CN201710273035 A CN 201710273035A CN 107154038 B CN107154038 B CN 107154038B
- Authority
- CN
- China
- Prior art keywords
- rib
- point
- matrix
- ribs
- points
- 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 67
- 238000003745 diagnosis Methods 0.000 title claims abstract description 19
- 238000012800 visualization Methods 0.000 title claims abstract description 11
- 208000027790 Rib fracture Diseases 0.000 title claims abstract description 10
- 210000001562 sternum Anatomy 0.000 claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims description 46
- 230000008569 process Effects 0.000 claims description 37
- 239000013598 vector Substances 0.000 claims description 23
- 210000000038 chest Anatomy 0.000 claims description 19
- 230000009466 transformation Effects 0.000 claims description 18
- 238000002591 computed tomography Methods 0.000 claims description 14
- 238000012937 correction Methods 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 9
- 230000011218 segmentation Effects 0.000 claims description 9
- 230000008602 contraction Effects 0.000 claims description 8
- 210000000988 bone and bone Anatomy 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 238000000844 transformation Methods 0.000 claims description 5
- 238000013519 translation Methods 0.000 claims description 4
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000003902 lesion Effects 0.000 claims description 2
- 230000003014 reinforcing effect Effects 0.000 claims description 2
- 208000010392 Bone Fractures Diseases 0.000 abstract description 8
- 206010017076 Fracture Diseases 0.000 abstract description 5
- 238000003759 clinical diagnosis Methods 0.000 abstract description 3
- 238000005452 bending Methods 0.000 abstract description 2
- 230000006872 improvement Effects 0.000 abstract description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 abstract 1
- 210000000056 organ Anatomy 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- 210000004204 blood vessel Anatomy 0.000 description 3
- 101100272279 Beauveria bassiana Beas gene Proteins 0.000 description 1
- 208000004434 Calcinosis Diseases 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002308 calcification Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 210000003734 kidney Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000007794 visualization technique Methods 0.000 description 1
Images
Classifications
-
- 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
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本发明涉及图像处理技术,旨在提供一种肋骨可视化的肋骨骨折辅助诊断方法。该种肋骨可视化的肋骨骨折辅助诊断方法包括过程:获取图片;提取肋骨;人工修正,去掉脊柱、胸骨;展开肋骨。本发明通过把肋骨展开投影到平面上,消除了肋骨弯曲遮挡对诊断的影响,相对与传统临床诊断方法,在平面上观察提高率诊断的直观性,本发明可以作为诊断微小骨折骨裂的辅助方法。
Description
技术领域
本发明是关于图像处理技术领域,特别涉及一种肋骨可视化的肋骨骨折辅助诊断方法。
背景技术
计算机断层扫描(CT)是诊断肋骨骨折的主要方法,但肋骨特殊的组织结构导致肋骨遮挡因而微小骨折及骨裂容易漏诊,精确的诊断需要经验的积累,常规的二维方法需要手动追踪多个切片,传统的三维可视化方法在诊断时需要对图像进行旋转观察从而避免肋骨遮挡的影响。常用重建方法有多层面重组法、最大密度投影法及表面阴影显示法。多层面重组法需要做一系列图像且非目标器官存在形变;最大密度投影法低密度组织器官易被遗漏且噪声影响大;表面阴影显示法容易产生伪影且看不出密度信息。
发明内容
本发明的主要目的在于克服现有技术中的不足,提供一种通过把肋骨投影到二维平面上,消除在临床诊断时肋骨遮挡导致微小骨折及骨裂的骨折辅助诊断方法。为解决上述技术问题,本发明的解决方案是:
提供一种肋骨可视化的肋骨骨折辅助诊断方法,包括下述过程:
一、获取图片;
二、提取肋骨;
三、人工修正,去掉脊柱、胸骨;
四、展开肋骨;
所述过程一是指:取得原始的人体胸部CT扫描图像;
所述过程二是指:将人体胸部CT扫描图像转化为256阶灰度图像,然后用增强滤波器对肋骨进行增强处理,调整像素值大小取一个范围提取出整个胸部骨骼,整个胸部骨骼包括肋骨、脊柱及胸骨,按照肋骨空间位置大小从上向下标记左右肋骨分支;
所述过程三是指:创建三维窗口,用户通过鼠标右键点击划圈标示出不属于肋骨的部分,计算机根据用户的标示将圈内部分(脊柱、胸骨)去除;
所述过程四是指:泊松曲面重建肋骨轮廓,拉普拉斯网格收缩提取肋骨骨架;应用最近点迭代矫正骨架,对准骨架位置;应用曲面重建算法(Curved Planar Reformation,CPR)对每一只肋骨进行展开,得到最终的肋骨可视化展开结果,用于辅助医生进行肋骨病变诊断。
在本发明中,所述过程二具体包括下述步骤:
步骤b:然后使用一个增强滤波算子对肋骨区域进行强化,从图像中逐点提取平滑的黑塞矩阵:
其中,所述I是指每点的二阶导数矩阵;所述Iij是指二阶导数,i、j分别取x、y或者z;
步骤c:黑塞矩阵I通过高斯核的二阶偏导数与原图像进行卷积求得:
其中,Gσ(v)=exp(-vtv/σ2),v∈R3为高斯核;所述Δ为拉普拉斯算符;所述为步骤a中的256色灰度图像;所述σ是指当前高斯核的标准差;所述t是指向量的转置;所述exp是指自然指数函数;所述*指的是函数间的卷积;所述R3为三维欧氏空间;所述v为欧氏空间中的一个三维向量;
步骤d:求出黑塞矩阵I的特征值λ1,λ2,λ3,且|λ1|<|λ2|<|λ3|;
计算每点的增强值:
其中,所述V0(I)用于表示像素值类似血管(vesselness)的程度;所述α、β、c都分别大于0(均为可调整经验参数,可根据实际灰度图像自行调节);
步骤e:通过调整像素值范围提取出胸部骨骼,包含肋骨、脊柱及胸骨;
步骤f:按照肋骨空间位置大小从上到下对左右肋骨进行标记。
在本发明中,所述过程四中,泊松曲面重建肋骨轮廓具体为:
输入经过过程三人工修正后提取的肋骨分割结果的表面点集,记为点集合S,该点集合满足任意元素s∈S,且任意元素s都具有两个已知参数:元素s的三维空间中对应的位置s.p、元素s的内法线方向同时,记Vec为三维欧式空间上的向量场,该向量场满足对所有的s∈S,向量场的值在点s.p上的值为而向量场在其余点的值为0:
在本发明中,所述过程四中,拉普拉斯网格收缩提取肋骨骨架是指:在提取骨架过程中保持了细节形变,通过估计肋骨上点云邻域建立散乱点之间的连接关系,然后构建拉普拉斯加权矩阵,将所有点作为全局的位置约束,进行点云收缩,对收缩后的点云进行聚类,连接成线得到肋骨骨架;具体为:
对于任意一个表面重建结果G=(V,E);其中,V为该表面的顶点集合(泊松重建所提取的等值面中的点);E是该表面重建结果中边的集合(表面网格中边的集合),是V中点的连线;假设G中的顶点数,即集合V的元素个数为n,记V={V1,...,Vn},Vi∈R3;定义并计算该表面重建结果的拉普拉斯矩阵为L,对于该表面重建结果G(V,E),记αij和βij为边(i,j)∈E的对角,其中(i,j)∈E代表从Vi到Vj的连线包含在E中,i,j分别为V中点的索引;定义角度权重:
wi,j=cotαi,j+cotβi,j
通过如下方式计算局部拉普拉斯矩阵L={Li,j}n×n
上式中,(i,j)中的i和j分别指代当前边的两个顶点的索引;k指代在边(i,k)中不同于当前索引i的那个顶点的索引;Li,j是指拉普拉斯矩阵第i行第j列的元素;
拉普拉斯网格收缩提取肋骨骨架利用一个迭代过程实现,记上一步骤中得到的肋骨泊松表面重建结果为GP=G1=(V1,E),其中G1,V1和E分别代表算法初始时的表面重建结果、与之对应的初始顶点集合以及边的集合,V1=VP,E=EP;设初始权重矩阵和
首先计算对应的拉普拉斯矩阵L1,在第t次迭代过程中,通过已知的第t步的结果其变量分别为第t步的表面重建结果,权重,以及拉普拉斯矩阵,算法重复以下步骤,并更新结果到第t+1次迭代中,直到当前迭代生成的表面重建结果体积为0,迭代过程具体为:
1)更新第t+1步的表面重建结果Gt+1(Vt+1,E),其中Vt+1是第t+1次更新后的顶点集合,通过求解最小二乘:[WL tLt WH t]TVt+1=[0 WH tVt]T;
其中,大写的T代表向量的转置,Vt为第t次迭代生成的顶点集合,WL t以及WH t为第t次迭代的权重矩阵;
3)通过表面重建结果Gt+1(Vt+1,E),计算第t+1步的拉普拉斯矩阵Lt+1;
迭代过程结束后,最后生成的表面重建结果为G∞(V∞,EP),并且获得该最后收敛的表面重建结果的骨架V∞,V∞为该表面的顶点集合。
在本发明中,所述过程四中,应用最近点迭代矫正骨架具体是指:对给定两组不同的三维数据点集,其中一个为源数据点集(待配准点集)P={P1,..,PN},另一个为目标点集Q={Q1,...,QM},其中M≥N,ICP试图找出一种刚性变换,将点集P与目标Q进行匹配,即ICP通过迭代的方法来寻找这样的最佳刚性匹配;
在第k次迭代过程中,ICP算法找到与源数据点集P中的每个点对应的目标点集Q中的点,并且计算最佳的刚性变换,使得下述距离最小:
具体来说ICP主要包括下述步骤:
(5)如果dk+1不小于给定的阈值δ,返回到(1)执行,直到dk+1<δ或迭代次数大于预设的最大迭代次数为止;
因源数据点集P为肋骨骨架内的骨骼点V∞,而点集Q为与之对应的肋骨二值分割结果(看作点集),通过上述应用最近点迭代矫正骨架的迭代算法,即实现将肋骨骨架与分割结果对齐。
在本发明中,所述过程四中,CPR方法展开肋骨具体包括下述步骤:
(1)应用得到的肋骨骨架(即通过应用最近点迭代矫正后的骨架),均匀选取上面的3D采样点;
(2)在每个采样点上计算基于相邻点和前一点的正交观测方向,在每个点附近采样一个平面;
(3)在每个点上重复上述过程,取一个切片,重复第二步;
(4)把采样的切片相互堆叠在一起形成一个直卷;
(5)重复每根肋骨从而把肋骨展开。
与现有技术相比,本发明的有益效果是:
本发明通过把肋骨展开投影到平而上,消除了肋骨弯曲遮挡对诊断的影响,相对与传统临床诊断方法,在平面上观察提高率诊断的直观性,本发明可以作为诊断微小骨折骨裂的辅助方法。
附图说明
图1为肋骨二值图。
图2为修正后去掉脊柱和胸骨。
图3为肋骨XZ方向展开图。
图4为肋骨YZ方向展开图。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述:
一种新的肋骨可视化肋骨骨折辅助诊断方法,包括下述过程:
一、获取图片;
二、提取肋骨;
三、人工修正,去掉脊柱、胸骨;
四、展开肋骨。
所述过程一具体包括:通过和医院联系得到人体胸部CT扫描图片。
所述过程二具体包括:
对腹部CT图像进行处理,是基于骨骼的密度比较大,在CT扫描中CT值比较大,表现为骨骼与其他组织器官相比亮度高,因而有明显区别,通过选取适当的窗宽窗位使用多尺度的增强滤波算子提取出肋骨,在保留骨骼区域的同时避免了其他组织和器官(如血管、钙化点、肾脏)的干扰。具体步骤是:
步骤b:然后使用一个增强滤波算子对目标区域进项强化,从图像中逐点提取平滑的黑塞矩阵
其中,所述I是指每点的二阶导数矩阵;所述Iij是指二阶导数,i、j分别取x、y或者z;
步骤c:黑塞矩阵I通过高斯核的二阶偏导数与原图像进行卷积求得:
其中,Gσ(v)=exp(-vtv/σ2),v∈R3为高斯核;所述△为拉普拉斯算符;所述为步骤a中的256色灰度图像;所述σ是指当前高斯核的标准差;所述t是指向量的转置;所述exp是指自然指数函数;所述*指的是函数间的卷积,所述R3为三维欧氏空间,所述v为欧氏空间中的一个三维向量;
步骤d:求出黑塞矩阵I的特征值λ1,λ2,λ3,且|λ1|<|λ2|<|λ3|;
计算每点的增强值:
其中,所述V0(I)用于表示像素值类似血管(vesselness)的程度;所述α、β、c都分别大于0(均为可调整经验参数,可根据实际灰度图像自行调节);
步骤e:通过调整像素值范围提取出胸部骨骼,包含胸骨、脊柱;
步骤f:按照肋骨空间位置大小从上到下对左右肋骨进行标记。
所述过程三具体包括:
创建三维窗口,用户通过鼠标右键点击划圈标示出不属于肋骨的部分,计算机根据用户的标示将圈内部分(脊柱、胸骨)去除。
所述过程四具体包括:
泊松曲面重建肋骨轮廓是通过估计模型的指示函数和提取等值面,对表面重建一个无缝的三角逼近;向量场V为空间M(在模型表面具有高分辨率,远离表面地方具有低分辨率)中的函数线性总和,构造并解泊松方程,提取求取的指示函数的等值面(它是一种全局的方法,把有法向点集的表面重建转化为一个空间的泊松问题,从而把表面重建问题表示为泊松方程的解)。
输入经过过程三人工修正后提取的肋骨分割结果的表面点集,记为点集合S,该点集合满足任意元素s∈S都具有两个已知参数:该元素s的三维空间中对应的位置s.p,以及该元素的内法线方向同时,记Vec为三维欧式空间上的向量场,该向量场满足对所有的s∈S,向量场的值在点s.p上的值为而向量场在其余点的值为0。泊松重建希望寻找三维空间上的一个标量函数χ,使得式了最小。其中符号||·||代表欧氏距离范数。泊松重建转而求解然后将作为上述χ的近似解。上式中△代表拉普拉斯算符,为三维空间上的一个标量函数。最后,泊松重建通过求取同属一个水平集的点集:获得表面重建结果的顶点集合VP,并对VP构造相应八叉树并获得对应表面重建结果GP。上式中VP代表待求水平集上的点集,q代表三维欧氏空间中的点,其中|·|代表集合的元素个数。GP代表泊松重建表面网格结果。记GP=(VP,EP),其中EP是表面网格中边的集合,EP中任意元素VP中两点间的连线,下同。
在提取骨架过程中保持了细节形变,通过估计肋骨上点云邻域建立散乱点之间的连接关系,然后构建拉普拉斯加权矩阵,将所有点作为全局的位置约束,进行点云收缩,对收缩后的点云进行聚类,连接成线得到肋骨骨架。
对任意一个表面重建结果G=(V,E),其中V为该表面的顶点集合。而E是该表面重建结果中边的集合,为V中点的连线。假设G中的顶点数,即集合V的元素个数为n。在此情况下,记V={V1,...,Vn},Vi∈R3。在这样的定义下,我们定义并计算该表面重建结果的拉普拉斯矩阵为L。对于该表面G(V,E),记αij和βij为边(i,j)∈E的对角,其中(i,j)∈E代表从Vi到Vj的连线包含在E中,i,j分别为V中点的索引。定义角度权重:
wi,j=cotαi,j+cotβi,j
我们通过如下方式计算局部拉普拉斯矩阵L={Li,j}n×n
上式中,(i,j)中的i和j分别指代当前边的两个顶点的索引。k指代在边(i,k)中不同于当前索引i的那个顶点的索引。Li,j是指拉普拉斯矩阵第i行第j列的元素。
拉普拉斯网格收缩提取肋骨骨架利用一个迭代过程实现,我们记上一步骤中得到的肋骨泊松表面重建结果为GP=G1=(V1,E),其中G1,V1和E分别代表算法初始时的表面重建结果、与之对应的初始顶点集合以及边的集合,V1=VP,E=EP。设初始权重矩阵和算法首先计算对应的拉普拉斯矩阵L1。在第t次迭代过程中,通过已知的第t步的结果其变量分别为第t步的表面重建结果,权重,以及拉普拉斯矩阵,算法重复以下步骤,并更新结果到第t+1次迭代中,直到当前迭代生成的表面重建结果体积为0:
1)更新第t+1步的表面重建结果Gt+1(Vt+1,E),其中Vt+1是第t+1次更新后的顶点集合,通过求解最小二乘:[WL tLt WH t]TVt+1=[0 WH tVt]T。其中大写的T代表向量的转置,Vt为第t次迭代生成的顶点集合,WL t以及WH t为第t次迭代的权重矩阵。
3)通过表面重建结果Gt+1(Vt+1,E)计算第t+1步的图像拉普拉斯矩阵Lt+1。
记最后生成的表面重建结果为G∞(V∞,EP),则该表面的骨架即为V∞,其中V∞为该表面的顶点集合。
在本专利中,算法将第三步生成的表面重建结果GP=(VP,EP)作为初始值即G1=(V1,E)带入上述算法,并且获得该表面重建结果的骨架。这个骨架是一系列点的集合,该集合即为最后收敛的表面重建结果G∞(V∞,EP)中的V∞。
应用最近点迭代矫正骨架具体是指:对给定两组不同的三维数据点集,其中一个为源数据点集(待配准点集)P={P1,...,PN},另一个为目标点集Q={Q1,...,QM},其中M≥N,ICP试图找出一种刚性变换,将点集P与目标Q进行匹配,即ICP通过迭代的方法来寻找这样的最佳刚性匹配)。
在第k次迭代过程中,ICP算法找到与源点集P中的每个点对应的目标点集Q中的点,并且计算最佳的刚性变换,使得下述距离最小:
其中,分别为第k步迭代的点集P,Q中的对应点,T=(R;t),R∈SO(3),t∈R3为刚性变换,Γ代表所有刚性变换的集合;其中SO(3)代表三阶正交群,||·||代表欧氏度量,符号代表函数复合。(所有参数及符号定义下同)
具体来说ICP主要包括下述步骤:
(5)如果dk+1不小于给定的阈值δ,返回到(1)执行,直到dk+1<δ或迭代次数大于预设的最大迭代次数为止;
因点集P为肋骨骨架内的骨骼点V∞,而点集Q为与之对应的肋骨二值分割结果(看作点集),通过该算法,即实现将肋骨骨架与分割结果对齐。
利用CPR方法展开肋骨:
(1)应用得到的肋骨骨架(即通过应用最近点迭代矫正后的骨架),均匀选取上面的3D采样点;
(2)在每个采样点上计算基于相邻点和前一点的正交观测方向,在每个点附近采样一个平面;
(3)在每个点上重复上述过程,取一个切片,重复第二步;
(4)把采样的切片相互堆叠在一起形成一个直卷;
(5)重复每根肋骨从而把肋骨展开。图3和图4是过程四得到的肋骨XZ和YZ方向的展开图,可以从图上观察到骨断裂处。
最后,需要注意的是,以上列举的仅是本发明的具体实施例。显然,本发明不限于以上实施例,还可以有很多变形。本领域的普通技术人员能从本发明公开的内容中直接导出或联想到的所有变形,均应认为是本发明的保护范围。
Claims (2)
1.一种肋骨可视化的肋骨骨折辅助诊断方法,其特征在于,包括下述过程:
一、获取图片;
二、提取肋骨;
三、人工修正,去掉脊柱、胸骨;
四、展开肋骨;
所述过程一是指:取得原始的人体胸部CT扫描图像;
所述过程二是指:将人体胸部CT扫描图像转化为256阶灰度图像,然后用增强滤波器对肋骨进行增强处理,调整像素值大小取一个范围提取出整个胸部骨骼,整个胸部骨骼包括肋骨、脊柱及胸骨,按照肋骨空间位置大小从上向下标记左右肋骨分支;
过程二具体包括下述步骤:
步骤b:然后使用一个增强滤波算子对肋骨区域进行强化,从图像中逐点提取平滑的黑塞矩阵:
其中,所述I是指每点的二阶导数矩阵;
步骤c:黑塞矩阵I通过高斯核的二阶偏导数与原图像进行卷积求得:
其中,Gσ(v)=exp(-vTv/σ2),v∈R3为高斯核;所述Δ为拉普拉斯算符;所述为步骤a中的256色灰度图像;所述σ是指当前高斯核的标准差;所述T是指向量的转置;所述exp是指自然指数函数;所述*指的是函数间的卷积;所述R3为三维欧氏空间;所述v为欧氏空间中的一个三维向量;
步骤d:求出黑塞矩阵I的特征值λ1,λ2,λ3,且|λ1|<|λ2|<|λ3|;
计算每点的增强值:
其中,所述V0(I)用于表示像素值类似管状结构的程度;所述α、β、c都分别大于0;
步骤e:通过调整像素值范围提取出胸部骨骼,包含肋骨、脊柱及胸骨;
步骤f:按照肋骨空间位置大小从上到下对左右肋骨进行标记;
所述过程三是指:创建三维窗口,用户通过鼠标右键点击划圈标示出不属于肋骨的部分,计算机根据用户的标示将圈内部分去除;
所述过程四是指:泊松曲面重建肋骨轮廓,拉普拉斯网格收缩提取肋骨骨架;应用最近点迭代矫正骨架,对准骨架位置;应用曲面重建算法对每一只肋骨进行展开,得到最终的肋骨可视化展开结果,用于辅助医生进行肋骨病变诊断;其中,
所述泊松曲面重建肋骨轮廓具体为:
输入经过过程三人工修正后提取的肋骨分割结果的表面点集,记为点集合S,该点集合满足任意元素s∈S,且任意元素s都具有两个已知参数:元素s的三维空间中对应的位置s.p、元素s的内法线方向同时,记Vec为三维欧式空间上的向量场,该向量场满足对所有的s∈S,向量场的值在点s.p上的值为而向量场在其余点的值为0;
所述拉普拉斯网格收缩提取肋骨骨架是指:在提取骨架过程中保持了细节形变,通过估计肋骨上点云邻域建立散乱点之间的连接关系,然后构建拉普拉斯加权矩阵,将所有点作为全局的位置约束,进行点云收缩,对收缩后的点云进行聚类,连接成线得到肋骨骨架;具体为:
对于任意一个表面重建结果G=(V,E);其中,V为该表面的顶点集合;E是该表面重建结果中边的集合,是V中点的连线;假设G中的顶点数,即集合V的元素个数为n,记V={V1,...,Vn},Vi∈R3;定义并计算该表面重建结果的拉普拉斯矩阵为L,对于该表面重建结果G(V,E),记αij和βij为边(i,j)∈E的对角,其中(i,j)∈E代表从Vi到Vj的连线包含在E中,i,j分别为V中点的索引;定义角度权重:
wi,j=cotαi,j+cotβi,j
通过如下方式计算局部拉普拉斯矩阵L={Li,j}n×n
上式中,(i,j)中的i和j分别指代当前边的两个顶点的索引;k指代在边(i,k)中不同于当前索引i的那个顶点的索引;Li,j是指拉普拉斯矩阵第i行第j列的元素;
拉普拉斯网格收缩提取肋骨骨架利用一个迭代过程实现,记上一步骤中得到的肋骨泊松表面重建结果为GP=G1=(V1,E),其中G1,V1和E分别代表算法初始时的表面重建结果、与之对应的初始顶点集合以及边的集合,V1=VP,E=EP;设初始权重矩阵和
首先计算对应的拉普拉斯矩阵L1,在第t次迭代过程中,通过已知的第t步的结果其变量分别为第t步的表面重建结果,权重,以及拉普拉斯矩阵,算法重复以下步骤,并更新结果到第t+1次迭代中,直到当前迭代生成的表面重建结果体积为0,迭代过程具体为:
1)更新第t+1步的表面重建结果Gt+1(Vt+1,E),其中Vt+1是第t+1次更新后的顶点集合,通过求解最小二乘:[WL tLt WH t]TVt+1=[0 WH tVt]T;
其中,大写的T代表向量的转置,Vt为第t次迭代生成的顶点集合,WL t以及WH t为第t次迭代的权重矩阵;
3)通过表面重建结果Gt+1(Vt+1,E),计算第t+1步的拉普拉斯矩阵Lt+1;
迭代过程结束后,最后生成的表面重建结果为G∞(V∞,EP),并且获得该最后收敛的表面重建结果的骨架V∞,V∞为该表面的顶点集合;
所述应用最近点迭代矫正骨架具体是指:对给定两组不同的三维数据点集,其中一个为源数据点集P={P1,...,PN},另一个为目标点集Q={Q1,...,QM},其中M≥N,ICP试图找出一种刚性变换,将点集P与目标Q进行匹配,即ICP通过迭代的方法来寻找这样的最佳刚性匹配;
在第m次迭代过程中,ICP算法找到与源数据点集P中的每个点对应的目标点集Q中的点,并且计算最佳的刚性变换T,使得下述距离最小:
具体来说ICP主要包括下述步骤:
(5)如果dm+1不小于给定的阈值δ,返回到(1)执行,直到dm+1<δ或迭代次数大于预设的最大迭代次数为止;
因源数据点集P为肋骨骨架内的骨骼点V∞,而点集Q为与之对应的肋骨二值分割结果,通过上述应用最近点迭代矫正骨架的迭代算法,即实现将肋骨骨架与分割结果对齐。
2.根据权利要求1所述的一种肋骨可视化的肋骨骨折辅助诊断方法,其特征在于,所述过程四中,曲面重建算法展开肋骨具体包括下述步骤:
(1)应用得到的肋骨骨架,均匀选取上面的3D采样点;
(2)在每个采样点上计算基于相邻点和前一点的正交观测方向,在每个点附近采样一个平面;
(3)在每个点上重复上述过程,取一个切片,重复第二步;
(4)把采样的切片相互堆叠在一起形成一个直卷;
(5)重复每根肋骨从而把肋骨展开。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610263681 | 2016-04-22 | ||
CN201610263681X | 2016-04-22 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107154038A CN107154038A (zh) | 2017-09-12 |
CN107154038B true CN107154038B (zh) | 2021-05-14 |
Family
ID=59793820
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710273035.6A Active CN107154038B (zh) | 2016-04-22 | 2017-04-21 | 一种肋骨可视化的肋骨骨折辅助诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107154038B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108320314B (zh) * | 2017-12-29 | 2021-07-09 | 北京优视魔方科技有限公司 | 一种基于ct横断图像的图像处理方法和装置、显示*** |
CN108446714B (zh) * | 2018-02-06 | 2021-09-24 | 山东科技大学 | 一种多工况下的非马尔科夫退化***剩余寿命预测方法 |
CN108491770B (zh) * | 2018-03-08 | 2023-05-30 | 李书纲 | 一种基于骨折影像的数据处理方法 |
CN108596877B (zh) * | 2018-03-28 | 2022-02-08 | 苏州科技城医院 | 肋骨ct数据分析*** |
CN109124662B (zh) * | 2018-07-13 | 2021-12-03 | 上海皓桦科技股份有限公司 | 肋骨中心线检测装置及方法 |
CN109035141B (zh) * | 2018-07-13 | 2023-03-24 | 上海皓桦科技股份有限公司 | 肋骨展开装置及方法 |
CN111968728B (zh) * | 2019-05-20 | 2024-03-08 | 杭州依图医疗技术有限公司 | 一种图像的处理方法和处理设备 |
CN110458799A (zh) * | 2019-06-24 | 2019-11-15 | 上海皓桦科技股份有限公司 | 基于肋骨展开图的肋骨骨折自动检测方法 |
CN111311626A (zh) * | 2020-05-11 | 2020-06-19 | 南京安科医疗科技有限公司 | 基于ct图像的颅骨骨折自动检测方法及电子介质 |
CN112349391A (zh) * | 2020-11-10 | 2021-02-09 | 山东大学齐鲁医院(青岛) | 一种优化肋骨自动标号方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101452577A (zh) * | 2008-11-26 | 2009-06-10 | 沈阳东软医疗***有限公司 | 一种肋骨自动标定的方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140071125A1 (en) * | 2012-09-11 | 2014-03-13 | The Johns Hopkins University | Patient-Specific Segmentation, Analysis, and Modeling from 3-Dimensional Ultrasound Image Data |
-
2017
- 2017-04-21 CN CN201710273035.6A patent/CN107154038B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101452577A (zh) * | 2008-11-26 | 2009-06-10 | 沈阳东软医疗***有限公司 | 一种肋骨自动标定的方法及装置 |
Non-Patent Citations (6)
Title |
---|
CPR - Curved Planar Reformation;Armin Kanitsar,et.al;《IEEE Visualization 2002》;20021231;第37-44页 * |
Skeleton Extraction by Mesh Contraction;,Au, O. K. C., et.al;《ACM transactions on graphics》;20080831;第27卷(第3期);第1-10页 * |
The ribs unfolded - a CT visualization algorithm for fast detection of rib fractures: effect on sensitivity and specificity in trauma patients;Helmut Ringl. et. al;《European Radiololgy》;20150214;第1865-1874页 * |
三维医学图像刚体配准;郭变芳;《优秀硕士论文全文数据库》;20111231(第14期);第5-6、30-31页 * |
基于CT图像的血管三维分割研究与应用;马骁;《优秀硕士论文全文数据库》;20130630(第6期);第12-17页 * |
基于泊松方程的三维表面重建算法的研究;张凯;《优秀硕士论文全文数据库》;20150331(第3期);第26-32页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107154038A (zh) | 2017-09-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107154038B (zh) | 一种肋骨可视化的肋骨骨折辅助诊断方法 | |
EP2916738B1 (en) | Lung, lobe, and fissure imaging systems and methods | |
CN110570515B (zh) | 一种利用ct图像进行人体骨骼三维建模的方法 | |
CN105719278B (zh) | 一种基于统计形变模型的医学图像分割方法 | |
EP3547207A1 (en) | Blood vessel extraction method and system | |
Baiker et al. | Atlas-based whole-body segmentation of mice from low-contrast Micro-CT data | |
US8055046B2 (en) | Shape reconstruction using X-ray images | |
JP6570145B2 (ja) | 画像を処理する方法、プログラム、代替的な投影を構築する方法および装置 | |
JP4885138B2 (ja) | 一連の画像における動き修正のための方法およびシステム | |
EP2056255B1 (en) | Method for reconstruction of a three-dimensional model of an osteo-articular structure | |
KR102527440B1 (ko) | 화상 해석 방법, 세그먼테이션 방법, 골밀도 측정 방법, 학습 모델 작성 방법 및 화상 작성 장치 | |
CN102651145B (zh) | 股骨三维模型可视化方法 | |
KR20150045885A (ko) | 초음파 및 ct 영상의 정합에 관한 시스템 및 작동 방법 | |
Chintalapani et al. | Statistical atlases of bone anatomy: construction, iterative improvement and validation | |
CN109509193B (zh) | 一种基于高精度配准的肝脏ct图谱分割方法及*** | |
CN105139377A (zh) | 一种腹部ct序列图像肝脏的快速鲁棒自动分割方法 | |
Chartrand et al. | Semi-automated liver CT segmentation using Laplacian meshes | |
JP6458166B2 (ja) | 医用画像処理方法及び装置及びシステム及びプログラム | |
CN114842154B (zh) | 一种基于二维x射线图像重建三维影像的方法和*** | |
CN115830016B (zh) | 医学图像配准模型训练方法及设备 | |
Hong et al. | Automatic lung nodule matching on sequential CT images | |
CN104933672B (zh) | 基于快速凸优化算法配准三维ct与超声肝脏图像的方法 | |
EP3389006B1 (en) | Rib unfolding from magnetic resonance images | |
Abdolali et al. | Mandibular canal segmentation using 3D Active Appearance Models and shape context registration | |
CN114298986A (zh) | 一种基于多视点无序x光片的胸腔骨骼三维构建方法及*** |
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 |