全局匹配优化的超声图像三维重建方法
【专利摘要】本发明公开全局匹配优化的超声图像三维重建方法。该方法包括：读取超声探头实时采集的二维超声图像，获取二维超声图像的大小、间隔两个参数；通过空间定位系统跟踪定位每一幅二维超声图像的空间位置，并获取每一幅二维超声造影图像与超声三维数据的相对位置关系；遍历二维超声图像中每一个像素，根据预设阈值，将像素灰度值赋值给超声三维数据中的体素；遍历超声三维数据，获取体素空缺区域，并检测其边界；计算空缺区域边界上每一个体素的修补权重系数；遍历整幅超声三维数据，寻找与修补权重系数最大的模块最匹配的模块，修补权重系数最大的模块，完成超声三维数据的重建。本发明能够快速精准的重建超声三维数据，为临床病灶诊疗提供理论指导。
【专利说明】全局匹配优化的超声图像三维重建方法
【技术领域】
[0001] 本发明涉及一种全局匹配优化的超声图像三维重建方法，适合于临床超声二维序 列图像的三维重建。
【背景技术】
[0002] 在现代医学影像学中，超声成像与CT、X射线、核医学、磁共振成像并驾齐驱，互为 补充。而超声成像以其强度低，对人体无损伤，无痛苦，显示方法多样而著称，尤其对人体软 组织的探测和心血管脏器的血流动力学观察有其独到之处，对孕妇胎儿的发育及健康情况 的检查更是具有不可替代的作用。超声辅助诊断和治疗的发展过程经历了从一维到二维， 从抽象到直观的过程。在B模式下，医生已经可以从人体内部器官的图像上清楚地看到器 官的病变情况，为诊断治疗提供依据。然而，二维B模式存在着其固有的不足：（I)B模式超 声图像是二维的，而实际人体器官组织是三维实体，因此，诊断医生必须在大脑里对多幅二 维图像进行整合，形成实际三维图像进行诊断。这样就使得诊断过程效率低，并且诊断结果 依靠诊断医生的临床经验，具有很强的主观性，且易造成误诊。（2) -些疾病的诊断常常需 要测量器官或病变的体积。例如，在前列腺癌检查时，除了要检测前列腺抗体的含量，还要 检测前列腺的大小，根据前列腺大小对应的抗体含量进行综合诊断。二维图像中对体积的 测量往往是利用平面图上的长轴与短轴作为参数的估计模型进行估计，测量精度低，变化 大，且主观性强。（3)基于二维的超声诊断不具有可重复性。在一次诊断中，无法重现同一 位置的图像，即使凑巧换能器正好在同一位置，人也无法断定是否在同一位置，因此在重复 检查或追踪疗效时，如果间隔时间短或无大的明显治疗效果，则无法与先前诊断进行比较， 也就无法对病情的发展状况或疗效进行有效的估计。（4)由于人体结构的固有限制，用二维 超声诊断时，一些部位的切面图片无法看到，对准确诊断也是一个无法克服的障碍。而在三 维诊断系统中，利用人体器官或病理组织的二维图像重建成三维图像，诊断医生可以在本 机或其它计算机上处理或测量这个三维图像。而且，可以无限制地显示器官组织的任意方 向和部位的二维切面图像。对每一个器官病理部位也可以进行标记，就可以很容易地与后 继治疗或复查的三维图像相比较，从而很方便地得知病情的发展或治疗的效果。
【发明内容】
[0003] 本发明提出了一种全局匹配优化的超声图像三维重建方法，能够在整幅超声三维 数据寻找最优匹配模块修补空缺体素，为超声图像疾病诊断提供了一种便利的工具。
[0004] 该全局匹配优化的超声图像三维重建方法，包括以下步骤：
[0005] 第一步：读取超声探头实时采集的二维超声图像，获取二维超声图像的大小、间隔 两个参数；
[0006] 第二步：通过空间定位来跟踪定位每一幅二维超声图像的空间位置，并获取每一 幅二维超声造影图像与超声三维数据的相对位置关系；
[0007] 第三步：遍历二维超声图像中每一个像素，将像素灰度值赋值给超声三维数据中 的体素；
[0008] 第四步：遍历超声三维数据，获取体素空缺区域，并检测其边界；
[0009] 第五步：计算空缺区域边界上每一个体素的修补权重系数；
[0010] 第六步：遍历整幅超声三维数据，寻找与修补权重系数最大的模块最匹配的模块， 修补权重系数最大的模块，完成超声三维数据的重建。
【专利附图】
【附图说明】
[0011] 图1是本发明所提出的工作流程图；
[0012] 图2是本发明所提出的空缺区域修补方法示意图；
[0013] 图3是本发明所提出的体素赋值方法示意图；
[0014] 图4是本发明的用户操作流程图。
【具体实施方式】
[0015] 关于本发明的优点与精神可以通过以下的发明详述及附图得到进一步的了解。 [0016]附图1为重建流程图，所述超声三维数据重建包括以下几个步骤：
[0017] 步骤S101，读取超声探头实时采集的二维超声图像，获取二维超声图像的大小、间 隔两个参数。
[0018] 定义获取得到的二维超声图像的大小为II1Xn2，图像间距为[&1，a2]。
[0019] 步骤S102,通过空间定位系统跟踪定位每一幅二维超声图像的空间位置，并获取 每一幅二维超声造影图像与超声三维数据的相对位置关系。
[0020] 在超声三维数据重建系统中存在三个坐标系，分别是超声三维数据所在的三维坐 标系0-XYZ，二维超声切片图像所在的二维坐标系0' -UV，以及定位跟踪系统所在的三维坐 标系Os-XsYsZ s。对于坐标系0' -UV与Os-XsYsZs，可以根据超声图像采集系统提供的超声图 像信息完成两个坐标系之间的相互转换，超声图像信息包括超声图像大小、图像间隔。定 义第j幅二维超声切片图像上一个像素的坐标为（Up%)，该幅超声切片图像的大小为 (N ju, &)，图像间隔为（S, Sp)，而由空间定位系统得到的该幅超声切片图像的空间坐标为 (Xs,，Y5,，Z S,)。如果坐标系o' -UV与Os-XsYsZs定义如图3所示，则像素点（Uj i, Vji)在坐标系 Os-Ws中的坐标(X'W计算公式为：
[0021]
【权利要求】
1. 全局匹配优化的超声图像三维重建方法，其特征在于，包括以下步骤： 第一步：读取超声探头实时采集的二维超声图像，获取二维超声图像的大小、间隔两个 参数； 第二步：通过空间定位来跟踪定位每一幅二维超声图像的空间位置，并获取每一幅二 维超声造影图像与超声三维数据的相对位置关系； 第三步：遍历二维超声图像中每一个像素，将像素灰度值赋值给超声三维数据中的体 素； 第四步：遍历超声三维数据，获取体素空缺区域，并检测其边界； 第五步：计算空缺区域边界上每一个体素的修补权重系数； 第六步：遍历整幅超声三维数据，寻找与修补权重系数最大的模块最匹配的模块，修补 权重系数最大的模块，完成超声三维数据的重建。
2. 如权利要求1所述的全局匹配优化的超声图像三维重建方法，其特征在于，第四步 中，根据canny算子检测得到超声三维数据中空缺区域的边界。
3. 如权利要求1或2所述的全局匹配优化的超声图像三维重建方法，其特征在于，第五 步中所述的体素的修补权重系数由该体素所在模块的灰度变化强度决定。
4. 如权利要求2所述的全局匹配优化的超声图像三维重建方法，其特征在于，Canny算 子检测的具体方法如下： Stepl :用高斯滤波器对图像滤波，以消除图像中的噪声； Step2 :对滤波后图像中的每个像素，计算其梯度幅值和方向； Step3 :对梯度幅值进行非极大值抑制，图像上（i，j)点处幅值A(i，j)的大小不能确定 该点是否为边缘点，需要采用非极大值抑制确定边缘点，即细化幅值图像中的屋脊带，保留 幅值局部变化最大的点； Step4 :用双阈值算法检测和连接边缘，对非极大值抑制幅值进行阈值化，得到边缘阵 列图像。
【文档编号】A61B8/00GK104306021SQ201410543394
【公开日】2015年1月28日 申请日期:2014年10月15日 优先权日:2014年10月15日
【发明者】杨健, 丛伟建, 刘越, 王涌天, 艾丹妮 申请人:北京理工大学