CN104517315A - 基于交互式区域生长法的双侧输尿管重建方法与系统 - Google Patents

基于交互式区域生长法的双侧输尿管重建方法与系统 Download PDF

Info

Publication number
CN104517315A
CN104517315A CN201410841890.9A CN201410841890A CN104517315A CN 104517315 A CN104517315 A CN 104517315A CN 201410841890 A CN201410841890 A CN 201410841890A CN 104517315 A CN104517315 A CN 104517315A
Authority
CN
China
Prior art keywords
image
ureter
interpolation
segmentation
reconstructing
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.)
Pending
Application number
CN201410841890.9A
Other languages
English (en)
Inventor
彭勇
金军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201410841890.9A priority Critical patent/CN104517315A/zh
Publication of CN104517315A publication Critical patent/CN104517315A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/08Volume rendering

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Graphics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了基于交互式区域生长法的双侧输尿管同时重建方法与系统。本方法基于MATLAB进行三维重建。其主要步骤为:首先,根据输尿管CT序列图,确定输尿管部位的灰度特性;对含有噪声的CT序列图像进行去噪、增强处理;采用交互式方法确定第一幅图像的种子点,使用二次分割法分割图像,分割出输尿管区域;利用种子区域迁移技术,进行序列分割;为了使重建表面更光滑,采用数据插值处理;最后,采用MATLAB体绘制技术,设置好颜色、光线。实验结果表明,采用本发明的方法能准确的重建出双侧输尿管,且有较快的重建速度。

Description

基于交互式区域生长法的双侧输尿管重建方法与系统
技术领域
本发明属于图像处理技术领域,涉及医学图像处理方法,可用于输尿管CT序列图像重构。
背景技术
基于交互式区域生长法的双侧输尿管同时重建方法与系统是为了解决医学CT图中的双侧输尿管腔道曲面重建并同时成像的问题。其实质是将一组输尿管的CT断层图像,通过本算法转化为具有三维立体效果的输尿管图像。双侧输尿管腔道曲面重建同时成像技术是诊断输尿管疾病的科学基础,也是输尿管疾病诊断治疗需要解决的关键技术之一。
如今,输尿管疾病目前已经成为威胁人类健康的主要病种之一。输尿管上接肾盂,下连膀胱,是一对细长的管道。它是泌尿系统的重要组成部分,常见有输尿管狭窄、输尿管积石、扩张等疾病。而目前以CT为代表的医学影像设备中取得的医学图像多为二维断层图像,对于病变部位的大小、形状、位置、以及周围组织的关系、患者病情严重情况的判断,主要依靠的是医生的临床经验。如果通过一组平行的二维切片重建出三维图像,那么就可以帮助医生对病变体或感兴趣区域有更加直观的认识,从而大大提高了临床诊断的准确性。医学图像三维重构作为计算机辅助诊断技术的一种,是辅助医生对病情判断的有效工具,在外科手术、远程医疗、放射治疗、虚拟内窥镜、医学教学等方面都有重要的应用价值。
针对输尿管部位三维重构,有许多方法被提出:
1.Jean-Francois于2006年提出基于韩国可视人体的三维泌尿生殖道重构。他根据韩国可视人计划中已分割好的泌尿生殖道图片,设计一个三维重构系统。
2.Bagley DH于1998年提出针对超声图像的人体输尿管三维重构。通过采集超声波单元发出的横截面图像,将其重构为三维模型。
而上述重构方中,没有具体的针对输尿管CT图,从序列图像分割到三维重构显示的方法。
发明内容
本发明针对已有技术的不足,提出基于交互式区域生长法的双侧输尿管重建方法与系统,使两侧输尿管腔道同时成像,提高医生诊断效率。
为实现上述目的,本发明包括如下步骤:
(1)选取10幅输尿管CT图,根据CT图中输尿管部位的灰度直方图,得到输尿管部位的灰度值范围,使输尿管部位灰度值g(x,y)满足|g(x,y)-center|<thre,center是灰度均值,thre是阈值;
(2)输入输尿管CT序列图,对每一幅图进行图像增强,采用各向异性滤波方法,对图像进行去噪、增强边缘处理;
(3)采用交互式方法,在序列图第一幅图像的输尿管部位上,人为指定若干个种子点;
(4)使用区域生长法,利用阈值thre、灰度均值center、灰度值T、梯度g等特征进行图像初分割;
(5)使用canny算子得到的边缘,结合前一步的梯度信息,进行二次区域生长,得到最终分割图像(输尿管部位);
(6)计算分割图像的有效像素点cont,自适应选择腐蚀算子a的大小,对图像进行形态学腐蚀处理;
(7)得到分割图像的腐蚀图后,将此图作为下一幅图像的种子点区域。判断有没有到最后一幅图像。如果是最后一幅图像,则跳到步骤(8);否则,执行步骤(4);
(8)得到序列分割图,对分割图进行采样插值处理;
(9)对插值后的序列分割图,进行MATLAB三维体绘制技术重构。
本发明与现有的技术相比具有以下优点:
1.本发明提出一套完整的双侧输尿管曲面重建同时成像系统,输尿管双侧同时成像,极大地方便了医生诊断病情。
2.本发明采用体绘制技术重建,处理三维采样信号方便,便于显示物体的内部结构。
3.本发明对重建环境进行设置,渲染重建结果,使其更接近真实输尿管。
4.本发明在MATLAB平台上实现三维重建,可充分发挥MATLAB的强大功能,而且有助于后续输尿管疾病的检测。
附图说明
图1是本发明的总流程图;
图2是用本发明对打造影剂输尿管的重构效果图;
图3是用本发明对未打造影剂输尿管的重构效果图。
具体实施方法
参照附图1,本发明的实现步骤如下:
步骤1.预处理数据,系统初始化
1a)选取10幅输尿管CT图,根据CT图中输尿管部位的灰度直方图,得到输尿管部位的灰度值范围,使输尿管部位灰度值g(x,y)满足|g(x,y)-center|<thre,center是灰度均值,thre是阈值;
1b)输入输尿管CT序列图,对每一幅图进行图像增强,采用各向异性滤波方法,对图像进行去噪、增强边缘处理;
1c)采用交互式方法,在序列图第一幅图像的输尿管部位上,人为指定若干个种子点。
步骤2.序列图像分割
2a)将CT图转换为梯度图,设高于某个值的梯度为边界;
2b)将前一步得到的种子点作为初始种子点;
2c)根据设定的条件,进行区域生长,得到分割图;
2d)使用canny算子得到边界,结合梯度得到的边界,根据生长规则,再次进行区域生长,得到最终分割图;
2e)计算分割图像的有效像素点cont,根据cont大小自适应选择腐蚀算子的大小,对图像进行形态学腐蚀处理;
2f)得到分割图像的腐蚀图后,将此图作为下一幅图像的种子点区域。判断有没有到最后一幅图像。如果是最后一幅图像,则进行图像插值、三维重构;否则,继续上述步骤,分割图像。
步骤3.图像采样插值、三维重构双侧输尿管
3a)得到序列分割图,对分割图进行采样插值处理。具体过程如下:
3a1)设定采样率和插值层数,根据采样率得到新的分割序列图;
3a2)利用数学形态学的方法,得到被插值图像区域的轮廓;
3a3)求出插值图像在上一层被插值图像和下一层被插值图像中所有对应的点,得到插值图像;
3a4)将插值图像和被插值图像融合,构成完整的分割序列图。
3b)对插值后的序列分割图,进行MATLAB三维体绘制技术重构。具体过程如下:
3b1)获取分割图的三维数据体;
3b2)根据实际情况,对数据集进行预处理,减少数据量,平滑数据;
3b3)计算数据集在显示平面的累计投影;
3b4)构造三维体重建碎片;
3b5)设置图像的颜色、阴影、及显示效果。
1)实验条件
本文中所使用的数据均为第四军医大学提供的两套CT人体输尿管CT扫描数据。一套为未打造影剂的输尿管CT扫描数据,共639张,数据格式为DICOM格式,分割出的输尿管部分为的298张输尿管CT数据,图像大小为512×512像素;另一套为打造影剂的输尿管CT扫描数据,共438张,数据格式为DICOM格式,分割出的输尿管部分为的230张输尿管CT数据,图像大小为512×512像素。
2)实验结果
重建出的双侧输尿管三维图像如图2、图3所示。从图2、图3可以看出:本发明利用MATLAB进行双侧输尿管三维重建的方法,能比较逼真地重建出的三维输尿管道模型,很好地反映了二维CT图像的立体效果,并可以应用于其它组织器官的三维重建,尤其是管腔重建。本研究解决了利用CT图像进行双侧输尿管腔道曲面重建并同时成像的问题,也为后续帮助医生进行更有效的治疗,提高医疗诊断的准确性打下了基础。

Claims (3)

1.基于交互式区域生长法的双侧输尿管腔道曲面同时重建方法与系统,包括如下步骤:
(1)选取10幅输尿管CT图,根据CT图中输尿管部位的灰度直方图,得到输尿管部位的灰度范围,使输尿管部位灰度g(x,y)满足|g(x,y)-center|<thre,center是灰度均值,thre是阈值;
(2)输入输尿管CT序列图,对每一幅图进行图像增强,采用各向异性滤波方法,对图像进行去噪、增强边缘处理;
(3)采用交互式方法,在序列图第一幅图像的输尿管部位上,人为指定若干个种子点;
(4)使用区域生长法,利用阈值thre、灰度均值center、灰度值T、梯度g等特征进行图像初分割;
(5)使用canny算子得到的边缘,结合前一步的梯度信息,进行二次区域生长,得到最终分割图像(输尿管部位);
(6)计算分割图像的有效像素点cont,自适应选择腐蚀算子a的大小,对图像进行形态学腐蚀处理;
(7)得到分割图像的腐蚀图后,将此图作为下一幅图像的种子点区域。判断有没有到最后一幅图像。如果是最后一幅图像,则跳到步骤(8);否则,执行步骤(4);
(8)得到序列分割图,对分割图进行采样插值处理;
(9)对插值后的序列分割图,进行MATLAB三维体绘制技术重构。
2.根据权利要求1所述的基于交互式区域生长法的双侧输尿管腔道曲面同时重建方法与系统,其中所诉步骤(4)采用区域生长法分割输尿管部分的CT图像,其具体步骤如下:
2a)将CT图转换为梯度图,设高于某个值的梯度为边界;
2b)将前一步得到的种子点作为初始种子点;
2c)根据设定的条件,进行区域生长,得到分割图。
3.根据权利要求1所述的基于交互式区域生长法的双侧输尿管腔道曲面同时重建方法与系统,其中所诉步骤(8)对分割序列图进行采样插值处理,其具体步骤如下:
3a)设定采样率和插值层数,根据采样率得到新的分割序列图;
3b)利用数学形态学的方法,得到被插值图像区域的轮廓;
3c)求出插值图像在上一层被插值图像和下一层被插值图像中所有对应的点,得到插值图像;
3d)将插值图像和被插值图像融合,构成完整的分割序列图。
CN201410841890.9A 2014-12-26 2014-12-26 基于交互式区域生长法的双侧输尿管重建方法与系统 Pending CN104517315A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410841890.9A CN104517315A (zh) 2014-12-26 2014-12-26 基于交互式区域生长法的双侧输尿管重建方法与系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410841890.9A CN104517315A (zh) 2014-12-26 2014-12-26 基于交互式区域生长法的双侧输尿管重建方法与系统

Publications (1)

Publication Number Publication Date
CN104517315A true CN104517315A (zh) 2015-04-15

Family

ID=52792557

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410841890.9A Pending CN104517315A (zh) 2014-12-26 2014-12-26 基于交互式区域生长法的双侧输尿管重建方法与系统

Country Status (1)

Country Link
CN (1) CN104517315A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105335742A (zh) * 2015-09-02 2016-02-17 清华大学 一种基于鲁棒特征统计的三维区域自适应分割方法
CN106875410A (zh) * 2017-02-09 2017-06-20 深圳市医诺智能科技发展有限公司 图像轮廓自动提取方法及系统
CN111311738A (zh) * 2020-03-04 2020-06-19 杭州市第三人民医院 一种采用影像学的输尿管3d数模建立方法及其数据采集装置
CN112132953A (zh) * 2020-08-24 2020-12-25 扬州哈工科创机器人研究院有限公司 一种ct三维重建方法及装置
CN112258522A (zh) * 2020-10-19 2021-01-22 哈尔滨体育学院 基于二次区域生长的武术竞赛区域分割方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101576997A (zh) * 2009-06-19 2009-11-11 西安电子科技大学 基于二次三维区域生长的腹部器官分割方法
CN103440665A (zh) * 2013-09-13 2013-12-11 重庆大学 膝关节软骨图像自动分割方法
CN103473805A (zh) * 2013-09-18 2013-12-25 内蒙古科技大学 基于改进区域生长算法测量三维重建肝脏模型体积的方法
CN103745470A (zh) * 2014-01-08 2014-04-23 兰州交通大学 基于小波的多边形轮廓演化医学ct图像交互式分割方法
US20140376776A1 (en) * 2007-05-21 2014-12-25 Cornell University. Method for segmenting objects in images

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140376776A1 (en) * 2007-05-21 2014-12-25 Cornell University. Method for segmenting objects in images
CN101576997A (zh) * 2009-06-19 2009-11-11 西安电子科技大学 基于二次三维区域生长的腹部器官分割方法
CN103440665A (zh) * 2013-09-13 2013-12-11 重庆大学 膝关节软骨图像自动分割方法
CN103473805A (zh) * 2013-09-18 2013-12-25 内蒙古科技大学 基于改进区域生长算法测量三维重建肝脏模型体积的方法
CN103745470A (zh) * 2014-01-08 2014-04-23 兰州交通大学 基于小波的多边形轮廓演化医学ct图像交互式分割方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
徐艳蕾: ""基于顺序形态学理论的医学CT图像三维重建方法的研究"", 《中国博士学位论文全文数据库 信息科技辑》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105335742A (zh) * 2015-09-02 2016-02-17 清华大学 一种基于鲁棒特征统计的三维区域自适应分割方法
CN106875410A (zh) * 2017-02-09 2017-06-20 深圳市医诺智能科技发展有限公司 图像轮廓自动提取方法及系统
CN106875410B (zh) * 2017-02-09 2020-07-03 深圳市医诺智能科技发展有限公司 图像轮廓自动提取方法及系统
CN111311738A (zh) * 2020-03-04 2020-06-19 杭州市第三人民医院 一种采用影像学的输尿管3d数模建立方法及其数据采集装置
CN111311738B (zh) * 2020-03-04 2023-08-11 杭州市第三人民医院 一种采用影像学的输尿管3d数模建立方法及其数据采集装置
CN112132953A (zh) * 2020-08-24 2020-12-25 扬州哈工科创机器人研究院有限公司 一种ct三维重建方法及装置
CN112258522A (zh) * 2020-10-19 2021-01-22 哈尔滨体育学院 基于二次区域生长的武术竞赛区域分割方法

Similar Documents

Publication Publication Date Title
CN107230206B (zh) 一种基于多模态数据的超体素序列肺部图像的3d肺结节分割方法
US9659405B2 (en) Image processing method and apparatus
CN108335304B (zh) 一种腹部ct扫描序列图像的主动脉瘤分割方法
CN104517315A (zh) 基于交互式区域生长法的双侧输尿管重建方法与系统
CN107545584A (zh) 医学图像中定位感兴趣区域的方法、装置及其系统
US11468570B2 (en) Method and system for acquiring status of strain and stress of a vessel wall
WO2022105623A1 (zh) 一种基于迁移学习的颅内血管病灶识别方法
CN110751651B (zh) 基于多尺度迁移学习的mri胰腺图像分割方法
CN105912874A (zh) 基于dicom医学图像构建的肝脏三维数据库系统
JP2002504385A (ja) 自動経路立案方法
CN106709920B (zh) 血管提取方法及其装置
CN109166104A (zh) 一种病变检测方法、装置及设备
CN100595791C (zh) 虚拟结肠镜的混合组织图像全自动分割方法
CN111815766A (zh) 基于2d-dsa图像重建血管三维模型处理方法及系统
CN105139377A (zh) 一种腹部ct序列图像肝脏的快速鲁棒自动分割方法
CN104112265B (zh) 结肠图像分割方法及装置
KR102182357B1 (ko) Ct 영상 내 간암 영역 기반으로 3차원 분석을 하는 수술 보조 장치 및 방법
Hu et al. Automated placenta segmentation with a convolutional neural network weighted by acoustic shadow detection
CN117529273A (zh) 从超声图像和相机图像重建3d图像的系统和方法
CN103324934A (zh) 基于平行结构检测与聚类的血管中心线自动提取方法
Dawood et al. The importance of contrast enhancement in medical images analysis and diagnosis
WO2014013285A1 (en) Apparatus and method for determining optimal positions of a hifu probe
CN109785296B (zh) 一种基于cta图像的三维球形指数测定方法
CN104751457A (zh) 一种新的基于变分能量的肝脏分割方法
Zhu et al. 3D automatic MRI level set segmentation of inner ear based on statistical shape models prior

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
DD01 Delivery of document by public notice

Addressee: Peng Yong

Document name: Notification of Publication and of Entering the Substantive Examination Stage of the Application for Invention

DD01 Delivery of document by public notice

Addressee: Jin Jun

Document name: Notification that Application Deemed to be Withdrawn

DD01 Delivery of document by public notice
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20150415

WD01 Invention patent application deemed withdrawn after publication