CN108320314B - 一种基于ct横断图像的图像处理方法和装置、显示系统 - Google Patents
一种基于ct横断图像的图像处理方法和装置、显示系统 Download PDFInfo
- Publication number
- CN108320314B CN108320314B CN201711484011.1A CN201711484011A CN108320314B CN 108320314 B CN108320314 B CN 108320314B CN 201711484011 A CN201711484011 A CN 201711484011A CN 108320314 B CN108320314 B CN 108320314B
- Authority
- CN
- China
- Prior art keywords
- rib
- image
- line
- pixel
- calculating
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 10
- 238000000034 method Methods 0.000 claims description 37
- 238000009877 rendering Methods 0.000 claims description 33
- 238000012545 processing Methods 0.000 claims description 26
- 239000013598 vector Substances 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 4
- 230000001360 synchronised effect Effects 0.000 claims description 4
- QTCANKDTWWSCMR-UHFFFAOYSA-N costic aldehyde Natural products C1CCC(=C)C2CC(C(=C)C=O)CCC21C QTCANKDTWWSCMR-UHFFFAOYSA-N 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- ISTFUJWTQAMRGA-UHFFFAOYSA-N iso-beta-costal Natural products C1C(C(=C)C=O)CCC2(C)CCCC(C)=C21 ISTFUJWTQAMRGA-UHFFFAOYSA-N 0.000 claims description 3
- 238000003745 diagnosis Methods 0.000 abstract description 14
- 208000027790 Rib fracture Diseases 0.000 abstract description 5
- 235000012431 wafers Nutrition 0.000 description 20
- 208000010392 Bone Fractures Diseases 0.000 description 16
- 206010017076 Fracture Diseases 0.000 description 15
- 238000004422 calculation algorithm Methods 0.000 description 6
- 210000000038 chest Anatomy 0.000 description 3
- 208000014674 injury Diseases 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000008733 trauma Effects 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 210000001991 scapula Anatomy 0.000 description 1
- 238000002603 single-photon emission computed tomography Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 208000011580 syndromic disease Diseases 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- 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
- A61B6/5223—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data generating planar views from image data, e.g. extracting a coronal view from a 3D image
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Computer Vision & Pattern Recognition (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)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种基于CT横断图像的图像处理方法,该方法包括:获取若干幅CT横断图像;基于所述CT横断图像生成多平面重建图像;根据肋骨轮廓从所述多平面重建图像提取肋骨中心线;基于所述肋骨中心线形成肋骨套管;将所述肋骨套管拉直,并计算所述套管内像素;依据所计算的套管内像素平面化显示肋骨图像。本发明实现了快速而有效的诊断肋骨骨折、减少诊断医生的重复工作强度、减少诊断时间、提高诊断准确率,在实际工作中具有很强的实用性和必要性。
Description
技术领域
本发明涉及医疗图像显示领域,尤其涉及一种基于CT横断图像的图像处理方法和装置、显示系统。
背景技术
随着社会发展,胸部外伤患者增多,CT是诊断肋骨骨折的最重要及最有效的检查手段,目前,传统的肋骨骨折CT诊断模式是以横断图像为主,辅助以三维图像定位。如图1所示,目前利用CT图像诊断的流程如下:a.逐层并完整浏览横断图像,发现骨折及可疑骨折并进行记录;b.在完整浏览横断图像后,转至肋骨VR图像,去掉床板及肩胛骨,对上述发现的骨折及可疑骨折在 VR图像上定位及观察骨折空间形态等;c.进一步旋转VR图像,以查找有无遗漏及可疑处,如果没有发现异常,d.则直接出诊断报告;如果发现异常,e.则将十字标定点拖动至可疑处并返回横断图像观察是否骨折,必要时辅以MPR 或CPR图像,如此反复,直至没有可疑点后,出诊断报告。
随着CT的发展,扫描层厚越来越薄,因此,在当前的诊断模式下,采用常用的1.25mm层厚对一位胸部外伤患者完成一次肋骨诊断,需要观察约300 幅横断图像,同时还需肋骨VR图像反复旋转观察,并返回横断图像以明确。鉴于肋骨数量多,走行不在同一平面且骨折细小难于观察等因素,诊断医生往往需要重复观察,以防漏诊及误诊,导致了诊断效率很低而且机械性重复工作量大。而此类患者往往很急,甚至有可能伴有其他部位损伤,因此,如何在 CT图像上快速而有效的诊断肋骨骨折、减少诊断医生的重复工作强度、减少诊断时间、提高诊断准确率,是需要解决的问题。
发明内容
鉴于上述问题,提出了本发明以便克服上述问题或者至少部分地解决上述问题。
本发明的一个方面,提供了一种基于CT横断图像的图像处理方法,该方法包括:
获取若干幅CT横断图像;
基于所述CT横断图像生成多平面重建图像;
基于肋骨的物理结构的轮廓,从所述多平面重建图像逐层根据正交面确定肋骨中心点并连接成相应肋骨中心线;
基于所述肋骨中心线形成肋骨套管;
将所述肋骨套管拉直,并计算所述套管内像素;
依据所计算的套管内像素平面化显示肋骨图像。
可选的,所述步骤基于所述肋骨中心线形成肋骨套管,具体包括:
计算第一肋骨线、第二肋骨线之间的距离d,所述第一肋骨线和第二肋骨线相邻;
将所述距离的一半作为第二肋骨线对应的肋骨管的半径,生成弯曲肋骨管。
可选的,将所述距离的一半作为第二肋骨线对应的肋骨管的半径,生成弯曲肋骨管,具体包括:
从根部向尖部逐像素遍历每根肋骨线,设为P2;
取相邻肋骨线同步调像素p1,p3;
可选的,将所述肋骨套管拉直,具体包括:
为每根肋骨线上的每个像素对应一个圆盘,圆盘所在平面与肋骨中线垂直,其半径由相邻肋骨线之间的距离确定;
确定圆盘上各像素;
将各像素确定后的圆盘平行摆放,以生成拉直的肋骨管。
可选的,利用体绘制方式确定圆盘上各像素。
可选的,将所述肋骨套管拉直,并计算所述套管内像素,具体包括:
计算单根肋骨从起始点到结束点的实际长度;
根据CT数据中单个体素的尺寸对肋骨进行遍历;
根据遍历过程中计算的圆片的法向量,计算体绘制所需要的光线和方向;
根据光线通过圆片的路径进行体绘制,计算出像素点;
组合所有像素点,形成套管内图像。
本发明还提供一种基于CT横断图像的图像处理装置,该装置包括:
获取单元,用于获取若干幅CT横断图像;
多平面重建单元,用于基于所述CT横断图像生成多平面重建图像;
肋骨中心线提取单元,用于根据肋骨轮廓从所述多平面重建正交图像逐层提取肋骨中心点并连接成中心线;
肋骨套管形成单元,用于基于所述肋骨中心线形成肋骨套管;
肋骨套管数据处理单元,用于将所述肋骨套管拉直,并计算所述套管内像素;
肋骨平面化显示单元,用于依据所计算的套管内像素平面化显示肋骨图像。
可选的,所述肋骨套管形成单元,具体包括:
第一计算模块,用于计算第一肋骨线、第二肋骨线之间的距离d,所述第一肋骨线和第二肋骨线相邻;
弯曲肋骨管形成模块,用于将所述距离的一半作为第二肋骨线对应的肋骨管的半径,生成弯曲肋骨管。
可选的,肋骨套管数据处理单元具体包括:
圆盘确定模块,用于为每根肋骨线上的每个像素对应一个圆盘,圆盘所在平面与肋骨中线垂直,其半径由相邻肋骨线之间的距离确定;
像素确定模块,用于确定圆盘上各像素;
圆盘平行处理模块,用于将各像素确定后的圆盘平行摆放,以生成拉直的肋骨管。
可选的,利用体绘制方式确定圆盘上各像素。
可选的,肋骨套管数据处理单元具体包括:
第二计算模块,用于计算单根肋骨从起始点到结束点的实际长度;
遍历处理模块,用于根据CT数据中单个体素的尺寸对肋骨进行遍历;
体绘制处理模块,用于完成下列步骤:根据遍历过程中计算的圆片的法向量,计算体绘制所需要的光线和方向;根据光线通过圆片的路径进行体绘制,计算出像素点;组合所有像素点,形成套管内图像。
本发明还提供一种肋骨图像显示系统,其包括前面所述的图像处理装置。
本申请实施例中提供的技术方案,至少具有如下技术效果或优点:
利用本发明所述的图像处理方法、装置能够在一个平面上展现各肋骨细节,从而能够帮助快速的确定肋骨骨折位置,避免医生在多画面之间来回切换,能够大大地减轻医生劳动强度,并降低误诊率。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂,以下特举本发明的具体实施方式。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1示出了根据现有技术中的医学影像诊断流程;
图2示出了根据本发明提出的基于CT横断图像的图像处理方法的流程图;
图3示出了肋骨套管的形成过程;
图4示出了一种具体实施方式的体绘制过程;
图5示出了一种体绘制过程中的圆片向量关系;
图6为拉直MIP示意图;
图7为拉直VR示意图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。
本发明从影像诊断和软件实现角度出发,实现的原理为首先要实现肋骨自动提取算法接口,对肋骨进行平铺显示,在平铺显示图中包含单根肋骨CPR (CPR:曲面重建)引导线及骨折标定点,根据该骨折标定点可切换到相应肋骨的VR、MPR显示图像,从而只做一次图像间的切换即可完成一处可疑骨折点的确认。
具体地,本发明提供了一种基于CT横断图像的图像处理方法,如图2所示,该方法包括:
S1.获取若干幅CT横断图像;
S2.基于所述CT横断图像生成多平面重建图像;
S3.根据肋骨轮廓从所述多平面重建正交图像逐层提取肋骨中心点并连接成中心线;
S4.基于所述肋骨中心线形成肋骨套管;
S5.将所述肋骨套管拉直,并计算所述套管内像素;
S5.依据所计算的套管内像素平面化显示肋骨图像。
本发明通过对不在同一平面内的肋骨图像的处理,能够在一个平面上展现各肋骨的图像细节,从而能够帮助医生通过观察一个平面即可初步判断肋骨骨折点,并可进行标记,避免了医生短时间内观察大量地CT图像横断图像。本发明生成的平面化肋骨显示图像中包含肋骨引导线及骨折标定点,如果医生需要进一步观察肋骨更进一步的细节特征,可通过肋骨引导线和骨折标定点快速定位到相应的三维重建VR(容积重建)图像,甚至可多角度旋转地观察三维重建图像,在需要的情况下,也可进一步定位到MPR(多平面重建)图像或者CPR(曲面重建)图像,因为很多时候,凭借对VR图像的多角度观察仍不能确定是否骨折,还需要进一步观察MPR(多平面重建)图像或者CPR(曲面重建)图像,从而帮助医生更准确地诊断病情,能够大大地降低误诊率。
可选的,所述步骤基于所述肋骨中心线形成肋骨套管,如图3所示,具体包括:
计算第一肋骨线11、第二肋骨线l2之间的距离d,所述第一肋骨线和第二肋骨线相邻;
将所述距离的一半作为第二肋骨线l2对应的肋骨管的半径,生成弯曲肋骨管,具体的算法可包括:
从根部向尖部逐像素遍历每根肋骨线,设为P2;
取相邻肋骨线同步调像素p1,p3;
在将肋骨套管拉直后,需要构建肋骨套管内的像素,以能够对肋骨图像进行完整显示,即像素的填充,可采用两种技术手段,一种是利用最大密度投影方式,一种是体绘制方式,都通过具体的算法来实现。最大密度投影方式具体是由一个或多个功能图像(206-PET,SPECT,MRI)组成一路输入信号,在模态是一种分子成像的情况下,可以使用不同的放射示踪剂来成像,将这路输入信号与三维图像(本申请优选了基于CT图像的多平面重建图像)进行对齐,对齐的方式不限,可以通过同步采集,也可能通过配准算法,从而能够生成弯曲的肋骨套管图像(最大密度投影(英文:maximal intensity projection,MIP)的定义,是一种应用广泛的CT及MR图像后处理技术。
MIP运用透视法获得二维图像,即通过计算沿着被扫描物每条射线上所遇到的最大密度像素而产生的。当光纤束通过一段组织的原始图像时,图像中密度最大的像素被保留,并被投影到一个二维平面上,从而形成MIP重建图像)MIP能反应相应像素的X线衰减值,较小的密度变化也能在MIP图像上显示,能很好地显示血管的狭窄、扩张、充盈缺损及区分血管壁上的钙化与血管腔内的对比剂。因为算法的关系,最大密度投影对于骨骼的走形表现较好,但对骨裂表现与体绘制方式相比较,稍欠一筹。因为体绘制方式会强化裂缝区域,因此本申请优选的方式是采用体绘制的方式。体绘制(体绘制的定义:将所有体数据(这里指CT数据)中的数据信息同时展现在二维图片上的技术,称之为体绘制技术。利用体绘制技术,可以在一幅图像中显示多种物质的综合分布情况,并且可以通过不透明度的控制,反应等值面及数据内部的情况。)的流程包括:
计算单根肋骨从起始点到结束点的实际长度;
根据CT数据中单个体素的尺寸对肋骨进行遍历;
根据遍历过程中计算的圆片的法向量,计算体绘制所需要的光线和方向;
根据光线通过圆片的路径进行体绘制,计算出像素点;
组合所有像素点,形成套管内图像。
下面说明采用体绘制的一种具体实施方式。
如图4所示,首先计算出单根肋骨从起始点到结束点的实际长度,可记为SingleBoneLength,根据CT数据单个体素的尺寸确定遍历步距和所述长度对整根肋骨进行遍历。CT数据单个体素的尺寸有的是0.5×0.5×0.8,单位mm,我们取最小数值0.5作为遍历步距,以保证最小误差,遍历的长度就是所述肋骨从起始点到结束点的实际长度。
遍历过程中,每步均可计算出一圆片的法向量,所述圆片可理解为垂直于肋骨中心线的一个小截面。计算过程如下:取得当前圆片的中心坐标,设为P1,取得下一个圆片的中心坐标,设为P2,则当前圆片的法向量为N1=P2-P1,因为相邻两圆片的法向量,可能会变化较大,过度不平滑,甚至产生跳变,可取前后n个法向量做平均,作为当前法向量,最末端的圆片,因为没有后续圆片,因此就不存在后续法向量,可直接沿用前面计算出的法向量。
根据圆片的法量,计算出体绘制所需要的光线的方向,具体过程可如下:定义一个向上的向量N0,通过两向量相乘的方式得到光线的向量NL=N0× N1;N1在圆片内的向量Nx=N1×NL,向量Nx也处于N1,N0所在的平面,特别需要注意的是如图5所示,Nx,NL均处于圆片所在的平面内。
圆片的尺寸根据上下肋骨中线的距离来确定,优选是上下肋骨中线的距离的一半。
根据光线(平行光线)通过圆片的过程,得到光线的起点和终点,从而取到光线通过圆片的路径,根据该路径进行体绘制,计算出结果像素点,如图6 所示,然后组合所有像素点,形成肋骨套管图像。
形成肋骨套管图像后,将所述肋骨套管拉直,具体包括:
为每根肋骨线上的每个像素对应一个圆盘,圆盘所在平面与肋骨中线垂直,其半径由相邻肋骨线之间的距离确定;
确定圆盘上各像素;
将各像素确定后的圆盘平行摆放,以生成拉直的肋骨套管。拉直的肋骨套管的图像可如图6或者图7所示。
本发明还提供一种基于CT横断图像的图像处理装置,该装置包括:
获取单元,用于获取若干幅CT横断图像;
多平面重建单元,用于基于所述CT横断图像生成多平面重建图像;
肋骨中心线提取单元,用于根据肋骨轮廓从所述多平面重建图像提取肋骨中心线;
肋骨套管形成单元,用于基于所述肋骨中心线形成肋骨套管;
肋骨套管数据处理单元,用于将所述肋骨套管拉直,并计算所述套管内像素;
肋骨平面化显示单元,用于依据所计算的套管内像素平面化显示肋骨图像。
该图像处理装置通过在通用计算机上运行计算机程序来完成,
可选的,所述肋骨套管形成单元,具体包括:
第一计算模块,用于计算第一肋骨线、第二肋骨线之间的距离d,所述第一肋骨线和第二肋骨线相邻;
弯曲肋骨管形成模块,用于将所述距离的一半作为第二肋骨线对应的肋骨管的半径,生成弯曲肋骨管。
可选的,肋骨套管数据处理单元具体包括:
圆盘确定模块,用于为每根肋骨线上的每个像素对应一个圆盘,圆盘所在平面与肋骨中线垂直,其半径由相邻肋骨线之间的距离确定;
像素确定模块,用于确定圆盘上各像素;
圆盘平行处理模块,用于将各像素确定后的圆盘平行摆放,以生成拉直的肋骨管。
可选的,利用体绘制方式确定圆盘上各像素。
可选的,肋骨套管数据处理单元具体包括:
第二计算模块,用于计算单根肋骨从起始点到结束点的实际长度;
遍历处理模块,用于根据CT数据中单个体素的尺寸对肋骨进行遍历;
体绘制处理模块,用于完成下列步骤:根据遍历过程中计算的圆片的法向量,计算体绘制所需要的光线和方向;根据光线通过圆片的路径进行体绘制,计算出像素点;组合所有像素点,形成套管内图像。
本发明还提供一种肋骨图像显示系统,其包括前面所述的图像处理装置。
本申请实施例中提供的技术方案,至少具有如下技术效果或优点:
利用本发明所述的图像处理方法、装置能够在一个平面上展现各肋骨细节,从而能够帮助快速的确定肋骨骨折位置,避免医生在多画面之间来回切换,能够大大地减轻医生劳动强度,并降低误诊率。
在此提供的算法和显示不与任何特定计算机、虚拟系统或者其它设备固有相关。各种通用系统也可以与基于在此的示教一起使用。根据上面的描述,构造这类系统所要求的结构是显而易见的。此外,本发明也不针对任何特定编程语言。应当明白,可以利用各种编程语言实现在此描述的本发明的内容,并且上面对特定语言所做的描述是为了披露本发明的最佳实施方式。
在此处所提供的说明书中,说明了大量具体细节。然而,能够理解,本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中,并未详细示出公知的方法、结构和技术,以便不模糊对本说明书的理解。
类似地,应当理解,为了精简本发明并帮助理解各个发明方面中的一个或多个,在上面对本发明的示例性实施例的描述中,本发明的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而,并不应将该公开的方法解释成反映如下意图:即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说,如下面的权利要求书所反映的那样,发明方面在于少于前面公开的单个实施例的所有特征。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其中每个权利要求本身都作为本发明的单独实施例。
本领域那些技术人员可以理解,可以对实施例中的设备中的模块进行自适应性地改变并且把它们设置在与该实施例不同的一个或多个设备中。可以把实施例中的模块或单元或组件组合成一个模块或单元或组件,以及此外可以把它们分成多个子模块或子单元或子组件。除了这样的特征和/或过程或者单元中的至少一些是相互排斥之外,可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述,本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
此外,本领域的技术人员能够理解,尽管在此的一些实施例包括其它实施例中所包括的某些特征而不是其它特征,但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如,在下面的权利要求书中,所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
本发明的各个部件实施例可以在一个或者多个处理器上运行的软件模块实现,或者以它们的组合实现。本领域的技术人员应当理解,可以在实践中使用微处理器或者数字信号处理器(DSP)来实现根据本发明实施例的网关、代理服务器、系统中的一些或者全部部件的一些或者全部功能。本发明还可以实现为用于执行这里所描述的方法的一部分或者全部的设备或者装置程序(例如,计算机程序和计算机程序产品)。这样的实现本发明的程序可以存储在计算机可读介质上,或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到,或者在载体信号上提供,或者以任何其他形式提供。
应该注意的是上述实施例对本发明进行说明而不是对本发明进行限制,并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中,不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的元件或步骤。本发明可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中,这些装置中的若干个可以是通过同一个硬件项来具体体现。
Claims (9)
1.一种基于CT横断图像的图像处理方法,其特征在于,该方法包括:
获取若干幅CT横断图像;
基于所述CT横断图像生成多平面重建图像;
从所述多平面重建图像逐层根据正交面确定肋骨中心点并连接成相应肋骨中心线;
基于所述肋骨中心线形成肋骨套管;
将所述肋骨套管拉直,并计算所述套管内像素,其中,具体包括:
为每根肋骨线上的每个像素对应一个圆盘,圆盘所在平面与肋骨中线垂直,其半径由相邻肋骨线之间的距离确定;确定圆盘上各像素;将各像素确定后的圆盘平行摆放,以生成拉直的肋骨管;
计算单根肋骨从起始点到结束点的实际长度;根据CT数据中单个体素的尺寸对肋骨进行遍历;根据遍历过程中计算的圆片的法向量,计算体绘制所需要的光线和方向;根据光线通过圆片的路径进行体绘制,计算出像素点;组合所有像素点,形成套管内图像;
依据所计算的套管内像素平面化显示肋骨图像。
2.根据权利要求1所述的方法,其特征还在于,基于所述肋骨中心线形成肋骨套管,具体包括:
计算第一肋骨线、第二肋骨线之间的距离d,所述第一肋骨线和第二肋骨线相邻;
将所述距离的一半作为第二肋骨线对应的肋骨管的半径,生成弯曲肋骨管。
4.根据权利要求1所述的方法,其特征还在于,利用体绘制方式确定圆盘上各像素。
5.一种基于CT横断图像的图像处理装置,其特征在于,该装置包括:
获取单元,用于获取若干幅CT横断图像;
多平面重建单元,用于基于所述CT横断图像生成多平面重建图像;
肋骨中心线提取单元,用于从所述多平面重建图像逐层根据正交面确定肋骨中心点并连接成相应肋骨中心线;
肋骨套管形成单元,用于基于所述肋骨中心线形成肋骨套管;
肋骨套管数据处理单元,用于将所述肋骨套管拉直,并计算所述套管内像素,具体包括:
圆盘确定模块,用于为每根肋骨线上的每个像素对应一个圆盘,圆盘所在平面与肋骨中线垂直,其半径由相邻肋骨线之间的距离确定;像素确定模块,用于确定圆盘上各像素;圆盘平行处理模块,用于将各像素确定后的圆盘平行摆放,以生成拉直的肋骨管;
第二计算模块,用于计算单根肋骨从起始点到结束点的实际长度;遍历处理模块,用于根据CT数据中单个体素的尺寸对肋骨进行遍历;体绘制处理模块,用于完成下列步骤:根据遍历过程中计算的圆片的法向量,计算体绘制所需要的光线和方向;根据光线通过圆片的路径进行体绘制,计算出像素点;组合所有像素点,形成套管内图像;
肋骨平面化显示单元,用于依据所计算的套管内像素平面化显示肋骨图像。
6.根据权利要求5所述的装置,其特征还在于,所述肋骨套管形成单元,具体包括:
第一计算模块,用于计算第一肋骨线、第二肋骨线之间的距离d,所述第一肋骨线和第二肋骨线相邻;
弯曲肋骨管形成模块,用于将所述距离的一半作为第二肋骨线对应的肋骨管的半径,生成弯曲肋骨管。
7.根据权利要求5所述的装置,其特征还在于,利用体绘制方式确定圆盘上各像素。
8.一种肋骨图像显示系统,其包括权利要求5-7任一项所述的图像处理装置。
9.一种存储介质,其存储有计算机程序,所述计算机程序被执行以实现权利要求1-4任一项所述方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711484011.1A CN108320314B (zh) | 2017-12-29 | 2017-12-29 | 一种基于ct横断图像的图像处理方法和装置、显示系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711484011.1A CN108320314B (zh) | 2017-12-29 | 2017-12-29 | 一种基于ct横断图像的图像处理方法和装置、显示系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108320314A CN108320314A (zh) | 2018-07-24 |
CN108320314B true CN108320314B (zh) | 2021-07-09 |
Family
ID=62893584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711484011.1A Active CN108320314B (zh) | 2017-12-29 | 2017-12-29 | 一种基于ct横断图像的图像处理方法和装置、显示系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108320314B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111968728B (zh) * | 2019-05-20 | 2024-03-08 | 杭州依图医疗技术有限公司 | 一种图像的处理方法和处理设备 |
CN111968069B (zh) * | 2019-05-20 | 2023-06-27 | 杭州依图医疗技术有限公司 | 一种图像的处理方法、处理装置和计算机可读存储介质 |
CN110807770A (zh) * | 2019-10-30 | 2020-02-18 | 杭州依图医疗技术有限公司 | 医学影像的处理、识别、显示方法及存储介质 |
CN111243082A (zh) * | 2020-01-08 | 2020-06-05 | 浙江大学 | 获得数字影像重建图像的方法、系统、装置及存储介质 |
CN112185519A (zh) * | 2020-09-01 | 2021-01-05 | 杭州依图医疗技术有限公司 | 一种医学图像处理方法、处理装置和医学图像的显示方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101923713A (zh) * | 2010-08-04 | 2010-12-22 | 中国科学院自动化研究所 | 一种提取冠状动脉血管中心线的方法 |
CN101987020A (zh) * | 2009-08-04 | 2011-03-23 | Ge医疗系统环球技术有限公司 | 倾斜图像扫描方法和重建方法及装置 |
CN103679642A (zh) * | 2012-09-26 | 2014-03-26 | 上海联影医疗科技有限公司 | 一种ct图像金属伪影校正方法、装置及ct设备 |
CN107106112A (zh) * | 2015-02-12 | 2017-08-29 | 株式会社日立制作所 | X射线ct装置、图像处理装置及图像重构方法 |
CN107194909A (zh) * | 2016-03-14 | 2017-09-22 | 东芝医疗系统株式会社 | 医用图像处理装置及医用图像处理程序 |
CN107204045A (zh) * | 2017-05-11 | 2017-09-26 | 青岛大学附属医院 | 基于ct图像的虚拟内窥镜系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8423124B2 (en) * | 2007-05-18 | 2013-04-16 | Siemens Aktiengesellschaft | Method and system for spine visualization in 3D medical images |
CN101452577B (zh) * | 2008-11-26 | 2010-12-29 | 沈阳东软医疗系统有限公司 | 一种肋骨自动标定的方法及装置 |
US9020233B2 (en) * | 2011-09-19 | 2015-04-28 | Siemens Aktiengesellschaft | Method and system for up-vector detection for ribs in computed tomography volumes |
US10991069B2 (en) * | 2014-10-08 | 2021-04-27 | Samsung Electronics Co., Ltd. | Method and apparatus for registration of medical images |
CN105825532A (zh) * | 2016-03-23 | 2016-08-03 | 长沙博为软件技术股份有限公司 | 一种优化的医学影像多平面重建方法 |
CN107154038B (zh) * | 2016-04-22 | 2021-05-14 | 孔德兴 | 一种肋骨可视化的肋骨骨折辅助诊断方法 |
-
2017
- 2017-12-29 CN CN201711484011.1A patent/CN108320314B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101987020A (zh) * | 2009-08-04 | 2011-03-23 | Ge医疗系统环球技术有限公司 | 倾斜图像扫描方法和重建方法及装置 |
CN101923713A (zh) * | 2010-08-04 | 2010-12-22 | 中国科学院自动化研究所 | 一种提取冠状动脉血管中心线的方法 |
CN103679642A (zh) * | 2012-09-26 | 2014-03-26 | 上海联影医疗科技有限公司 | 一种ct图像金属伪影校正方法、装置及ct设备 |
CN107106112A (zh) * | 2015-02-12 | 2017-08-29 | 株式会社日立制作所 | X射线ct装置、图像处理装置及图像重构方法 |
CN107194909A (zh) * | 2016-03-14 | 2017-09-22 | 东芝医疗系统株式会社 | 医用图像处理装置及医用图像处理程序 |
CN107204045A (zh) * | 2017-05-11 | 2017-09-26 | 青岛大学附属医院 | 基于ct图像的虚拟内窥镜系统 |
Non-Patent Citations (1)
Title |
---|
多层螺旋CT不同重建方法在诊断肋骨骨折中的对比研究;李云卿等;《中国CT和MRI杂志》;20051031;第3卷(第4期);第25-27页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108320314A (zh) | 2018-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108320314B (zh) | 一种基于ct横断图像的图像处理方法和装置、显示系统 | |
US8625869B2 (en) | Visualization of medical image data with localized enhancement | |
US9471987B2 (en) | Automatic planning for medical imaging | |
US9401047B2 (en) | Enhanced visualization of medical image data | |
US10497123B1 (en) | Isolation of aneurysm and parent vessel in volumetric image data | |
US9478048B2 (en) | Prior image based three dimensional imaging | |
US9424680B2 (en) | Image data reformatting | |
US20130190602A1 (en) | 2d3d registration for mr-x ray fusion utilizing one acquisition of mr data | |
JP5415068B2 (ja) | カーブした細長い構造の切断面の視覚化 | |
US10460508B2 (en) | Visualization with anatomical intelligence | |
JP2010528750A (ja) | 管状構造の検査 | |
US9691157B2 (en) | Visualization of anatomical labels | |
WO2012061452A1 (en) | Automatic image-based calculation of a geometric feature | |
US9019272B2 (en) | Curved planar reformation | |
CN115082348A (zh) | 图像重建方法及装置 | |
US8938107B2 (en) | System and method for automatic segmentation of organs on MR images using a combined organ and bone atlas | |
CN102573640B (zh) | 医用图像处理装置及医用图像处理方法 | |
US20120078101A1 (en) | Ultrasound system for displaying slice of object and method thereof | |
US20180092608A1 (en) | Reconstruction of Flow Data | |
US9786069B2 (en) | Refined reconstruction of time-varying data | |
US9585569B2 (en) | Virtual endoscopic projection image generating device, method and program | |
Kiraly et al. | A novel visualization method for the ribs within chest volume data |
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 |