CN1963869A - 空腔性脏器内壁快速体绘制方法 - Google Patents
空腔性脏器内壁快速体绘制方法 Download PDFInfo
- Publication number
- CN1963869A CN1963869A CNA200610118894XA CN200610118894A CN1963869A CN 1963869 A CN1963869 A CN 1963869A CN A200610118894X A CNA200610118894X A CN A200610118894XA CN 200610118894 A CN200610118894 A CN 200610118894A CN 1963869 A CN1963869 A CN 1963869A
- Authority
- CN
- China
- Prior art keywords
- binaryzation
- volume data
- data
- volume
- inwall
- 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.)
- Granted
Links
Landscapes
- Image Generation (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及的是一种图像处理技术领域的空腔性脏器内壁快速体绘制方法。首先对医学原始图像进行分割，得到脏器的空腔部分的二值化体数据，接着对这一体数据进行三维形态学膨胀运算，然后将膨胀后的二值化体数据与原二值化体数据进行异或操作，得到脏器壁二值化体数据，再利用这一数据在医学原始图像中找出所有前景体素，最后对这些前景体素进行体绘制。本发明较现有技术的改进在于在体绘制时去掉了所有离内壁较远的体素，可以大大减少体绘制的计算量，提高体绘制的速度，并且能保持空腔性脏器内壁的原始形态结构，保证体绘制效果的真实性。
Description
技术领域
本发明涉及的是一种图像处理技术领域的方法，具体涉及一种空腔性脏器内壁快速体绘制方法。
技术背景
目前，利用CT、MRI等医学成像数据，对空腔性脏器内壁进行三维可视化重建主要采用的是虚拟内窥镜技术。这一技术将视角置于脏器内部的空腔内，然后以内部的视角对脏器的内壁进行观察，从而模拟传统光学内窥镜的观察效果。在目前的虚拟内窥镜技术中，对脏器内壁进行重建主要有面绘制和体绘制两种方法。面绘制是基于对分割后的空腔体数据等值面的抽取，只对内壁的表面进行重建，故效果不够真实。体绘制无须经过数据分割，是将视点放在空腔内部，直接对所有体数据进行绘制，这样效果更为真实，然而重建的速度较慢，对每一帧图像的计算时间一般会达到秒级，无法实现实时漫游体绘制。
经对现有技术的文献检索发现，耿国华等在《计算机应用》(2002年第22卷第11期，第54-55页)上发表的“交互式实时虚拟内窥镜***中的关键技术”中提出了一种基于动态精简来提高体绘制速度的方法。采用这种方法，每次进行体绘制时只需对一个四棱锥内的数据进行绘制而无需考虑其它体素。其不足之处在于，这种方法对空腔性脏器内壁进行体绘制时，仍然会有很多离开内壁表面有一定距离的体素被考虑进行了绘制运算。但是，这些体素对当前视角所观察到的图像没有任何影响，在微机平台对256×256×69规模的数据作整体绘制仍需1秒左右，从而影响了体绘制的速度。
发明内容
本发明的目的在于针对现有虚拟内窥镜技术中体绘制速度上的不足，提出一种空腔性脏器内壁快速体绘制方法。本发明利用三维形态学膨胀运算，来最大限度地减少体绘制时所需计算的体素，但同时又不影响体绘制效果。
本发明是通过以下技术方案实现的，本发明首先对医学原始图像进行分割，得到脏器的空腔部分的二值化体数据，接着对这一体数据进行三维形态学膨胀运算，然后利用膨胀后的二值化体数据与原二值化体数据，得到脏器壁二值化体数据，再利用这一数据在医学原始图像中找出所有前景体素，最后对这些前景体素进行体绘制。
所述的医学原始图像，是指利用CT、MRI等医学成像手段，所获得的包含待可视化重建的空腔性脏器的医学图像。
所述的进行分割，得到脏器的空腔部分的二值化体数据，是指，对医学原始图像采用图像分割的方法，得到空腔部分的体素，然后重新建立一个二值化体数据，使得其在空腔内部的值为1，外部的值都为0。
所述的对这一体数据进行三维形态学膨胀运算，是指，对刚刚得到的脏器空腔部分的二值化体数据进行三维数学形态学膨胀运算，得到一个新的膨胀后的二值化体数据。
所述的利用膨胀后的二值化体数据与原二值化体数据，得到脏器壁二值化体数据，是指对原空腔的二值化体数据中的每个体素值与膨胀后的二值化体数据中的每个体素值分别进行异或运算，从而得到一个新的二值化体数据。在这个体数据中，只在膨胀运算对空腔扩大的部分的体素值为1，其余都为0，故可称之为脏器壁二值化体数据。
所述的利用这一数据在医学原始图像中找出所有的前景体素，是指，利用前面得到的脏器壁二值化体数据，给出其中每一个值为1的体素的坐标，然后在医学原始数据中找出所有这些坐标上的点，作为前景体素，并将其余点都看作背景体素，值置为零。
所述的对这些前景体素进行体绘制，是指，不考虑所有值为零的体素，只考虑前述所得到的前景体素，进行体绘制。体绘制的过程中，仍可以进一步加入动态精简等加速的方法。
本发明相对于现有技术的优势在于：可以在不影响内壁绘制效果，保持其真实性准确性的前提下，排除离内壁距离相差超过所预设数值的体素，使得体绘制的计算量可以减少到普通体绘制方法的十分之一或更多，从而加快体绘制的速度达十倍以上。
附图说明
图1空腔形态学膨胀示意图
图2空腔体数据与膨胀后体数据异或示意图
图3体绘制示意图
具体实施方式
下面结合附图对本发明的实施例作详细说明：本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和过程，但本发明的保护范围不限于下述的实施例。
本实施例具体步骤如下：
(1)从医学原始图像中分割出空腔二值化体数据。
(2)对空腔二值化体数据进行形态学膨胀运算。
(3)将膨胀后的体数据与原空腔体数据进行异或运算，得到脏器壁体数据。
(4)在医学原始图像中找出所对应的体素，并将其作为前景，其余点均作为背景置零，不予考虑。
(5)对前景体素进行体绘制。
该实施例采用人体下腹部CT扫描断层数据，要对结肠内壁进行内窥体绘制。整个实施例的实现过程如下：
1.读入人体下腹部CT扫描断层数据，利用区域生长法分割出结肠内部空腔，并令空腔内的体素值为1，其余为0。
2.对结肠空腔数据进行形态学膨胀运算。膨胀的结构元素为球体，球体半径的选取跟所需体绘制效果的精度有关。半径选得越大，体绘制的精度越高，计算量也越大；半径选得越小，体绘制的精度降低，计算量也随之降低。在本实施例中，半径先取为10个体素。膨胀效果的示意图如图1所示。经过膨胀运算，结肠空腔的外缘得到扩大。
3.将膨胀后的体数据同原空腔二值化体数据进行异或运算，得到一个新的结肠壁二值化体数据。图2是这一步骤的示意图。经过这一步操作，结肠空腔内部全都变成零，只有结肠壁部分数据得到保留。
4.对结肠壁二值化体数据，给出其中每一个值为1的体素的坐标，然后在原始CT断层数据中找出所有这些坐标上的点，作为前景体素，并将其余点都看作背景体素，值置为零。
5.不考虑所有背景体素，只对前景体素进行体绘制计算。体绘制的示意图如图3所示。体绘制的视角置于空腔内部，体绘制区域仅为阴影部分即结肠壁的数据。
这样的实施方式，在体绘制时去掉了所有离内壁较远，超过10个体素的数据，这样可以大大减少体绘制的计算量。跟考虑一定区域内的全部空腔外部体素相比，计算量减少到原来的二十分之一左右，体绘制的速度提高为原来的二十倍，并且最后的绘制效果能保持空腔性脏器内壁的原始形态结构，保证真实性。
Claims (7)
1、一种空腔性脏器内壁快速体绘制方法，其特征在于，首先对医学原始图像进行分割，得到脏器的空腔部分的二值化体数据，接着对这一体数据进行三维形态学膨胀运算，然后利用膨胀后的二值化体数据与原二值化体数据，得到脏器壁二值化体数据，再利用这一数据在医学原始图像中找出所有前景体素，最后对这些前景体素进行体绘制。
2、根据权利要求1所述的一种空腔性脏器内壁快速体绘制方法，其特征是，所述的医学原始图像，是指利用医学成像手段，所获得的包含待可视化重建的空腔性脏器的医学图像。
3、根据权利要求1所述的一种空腔性脏器内壁快速体绘制方法，其特征是，所述的进行分割，得到脏器的空腔部分的二值化体数据，是指，对医学原始图像采用图像分割的方法，得到空腔部分的体素，然后重新建立一个二值化体数据，使得其在空腔内部的值为1，外部的值都为0。
4、根据权利要求1所述的一种空腔性脏器内壁快速体绘制方法，其特征是，所述的对这一体数据进行三维形态学膨胀运算，是指，对刚刚得到的脏器空腔部分的二值化体数据进行三维数学形态学膨胀运算，得到一个新的膨胀后的二值化体数据。
5、根据权利要求1所述的一种空腔性脏器内壁快速体绘制方法，其特征是，所述的利用膨胀后的二值化体数据与原二值化体数据，得到脏器壁二值化体数据，是指对原空腔的二值化体数据中的每个体素值与膨胀后的二值化体数据中的每个体素值分别进行异或运算，从而得到一个新的二值化体数据。
6、根据权利要求1所述的一种空腔性脏器内壁快速体绘制方法，其特征是，所述的利用这一数据在医学原始图像中找出所有的前景体素，是指，利用前面得到的脏器壁二值化体数据，给出其中每一个值为1的体素的坐标，然后在医学原始数据中找出所有这些坐标上的点，作为前景体素，并将其余点都看作背景体素，信置为零。
7、根据权利要求1所述的一种空腔性脏器内壁快速体绘制方法，其特征是，所述的对这些前景体素进行体绘制，是指，不考虑所有值为零的体素，只考虑前述所得到的前景体素，进行体绘制。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200610118894A CN100583159C (zh) | 2006-11-30 | 2006-11-30 | 空腔性脏器内壁快速体绘制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200610118894A CN100583159C (zh) | 2006-11-30 | 2006-11-30 | 空腔性脏器内壁快速体绘制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1963869A true CN1963869A (zh) | 2007-05-16 |
| CN100583159C CN100583159C (zh) | 2010-01-20 |
Family
ID=38082909
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN200610118894A Expired - Fee Related CN100583159C (zh) | 2006-11-30 | 2006-11-30 | 空腔性脏器内壁快速体绘制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100583159C (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8064723B2 (en) | 2007-04-19 | 2011-11-22 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Method and apparatus for fast volume rendering of 3D ultrasound image |
| CN101561930B (zh) * | 2009-05-07 | 2012-06-13 | 上海交通大学 | 改良的管道状物体虚拟外翻可视化方法 |
| CN103679808A (zh) * | 2013-12-24 | 2014-03-26 | 通号通信信息集团有限公司 | 利用二维核磁图像重建头部三维模型的方法及*** |
| CN109035261A (zh) * | 2018-08-09 | 2018-12-18 | 北京市商汤科技开发有限公司 | 医疗影像处理方法及装置、电子设备及存储介质 |
-
2006
- 2006-11-30 CN CN200610118894A patent/CN100583159C/zh not_active Expired - Fee Related
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8064723B2 (en) | 2007-04-19 | 2011-11-22 | Shenzhen Mindray Bio-Medical Electronics Co., Ltd. | Method and apparatus for fast volume rendering of 3D ultrasound image |
| CN101561930B (zh) * | 2009-05-07 | 2012-06-13 | 上海交通大学 | 改良的管道状物体虚拟外翻可视化方法 |
| CN103679808A (zh) * | 2013-12-24 | 2014-03-26 | 通号通信信息集团有限公司 | 利用二维核磁图像重建头部三维模型的方法及*** |
| CN109035261A (zh) * | 2018-08-09 | 2018-12-18 | 北京市商汤科技开发有限公司 | 医疗影像处理方法及装置、电子设备及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN100583159C (zh) | 2010-01-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107230206B (zh) | 一种基于多模态数据的超体素序列肺部图像的3d肺结节分割方法 | |
| CN109389590B (zh) | 结肠图像数据处理***及方法 | |
| CN107808156A (zh) | 感兴趣区域提取方法 | |
| CN107590809A (zh) | 肺分割方法及医学成像*** | |
| CN107067398B (zh) | 用于三维医学模型中缺失血管的补全方法及装置 | |
| CN107595390B (zh) | 一种超声影像与ct影像的实时匹配融合方法 | |
| CN107895364B (zh) | 一种用于虚拟手术术前规划的三维重建*** | |
| JP5846841B2 (ja) | 医用画像表示装置 | |
| CN1973298A (zh) | 显示气管支气管树 | |
| JPS63186628A (ja) | ３次元画像処理装置 | |
| CN103236058A (zh) | 获取四维心脏图像感兴趣体积的方法 | |
| CN100583159C (zh) | 空腔性脏器内壁快速体绘制方法 | |
| CN102779354B (zh) | 基于光度立体技术的中医望诊信息表面三维重建方法 | |
| CN105139377A (zh) | 一种腹部ct序列图像肝脏的快速鲁棒自动分割方法 | |
| CN103077557A (zh) | 一种自适应分层次胸部大数据显示的实现方法 | |
| CN110882057A (zh) | 一种退针布针***及其定位方法 | |
| CN113643176A (zh) | 一种肋骨显示方法和装置 | |
| CN107204045A (zh) | 基于ct图像的虚拟内窥镜*** | |
| CN108510580A (zh) | 一种椎骨ct图像三维可视化方法 | |
| CN114445431A (zh) | 一种医学三维影像任意切割的方法与装置 | |
| CN1430185A (zh) | 基于单层表面跟踪的超大规模医学图像表面重建方法 | |
| CN101996415B (zh) | 眼球的三维建模方法 | |
| CN102548480B (zh) | 医用图像显示设备和方法 | |
| US9585569B2 (en) | Virtual endoscopic projection image generating device, method and program | |
| CN105869169B (zh) | 一种肿瘤组织微排列图像的自动分割方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100120 Termination date: 20121130 |