一种二维超声图像与穿刺针共面判断方法技术

本发明专利技术涉及一种二维超声图像与穿刺针共面判断方法，手持超声探头扫描组织内穿刺针，保存组织内穿刺针的超声图像，其特征在于:步骤1：提取超声图像上的成像区域，将彩色图像转换为灰度图像，经阈值变换后转换成二值图像；步骤2：对步骤1获得的二值图像进行边缘检测，通过Radon变换确定穿刺针所在的直线；步骤3：定位穿刺针针尖所在的位置；步骤4：通过判断穿刺针所在直线与超声成像区域边界交点的灰度值是否大于某一阈值，确定穿刺针是否位于成像平面内。

【技术实现步骤摘要】

本专利技术涉及。
技术介绍
图像引导外科系统利用医学影像和计算机图形图像技术，可在术前对患者多模态 图像数据进行三维重建和可视化，获得三维模型，制定合理、定量的手术计划，开展术前模 拟；在术中利用三维空间定位设备进行图像和病人物理空间的配准，把患者的实际体位、手 术器械的实时空间位置映射到患者的三维图像空间，对手术器械在空间中的位置实时采集 并显示，医生通过观察三维图像中手术器械与病变部位的相对位置关系，对病人进行精确 的手术治疗。这种新兴外科手术不仅能有效地解决病症，而且还有不开刀、创伤小、恢复快、 效果好优点，已经广泛应用于临床。虽然图像引导外科手术有诸多优点，但手术中引导手术 工具到达预定目标组织的跟踪定位仍然是其应用的技术难点。目前，术中图像引导方法有 computed tomography (CT)，X光透视成像，核磁共振成像和超声成像。与CT和核磁共振等 成像方式相比较，超声成像具有廉价、实时、无损伤、无辐射和便携等优势，已经作为一种实 时的术中成像技术被引入，用于软组织及血流成像。特别是近年来发展的图像引导介入手 术，超声成像凸显其直观、实时、易于识别等特点，更方便外科医生观察目标部位的解剖特 征，判断术具与目标部位的相对位置关系。目前超声成像系统成像模式包括二维成像和三 维成像。作为最新的成像方法，三维成像虽然能够实时地显示人体内部组织的形态，然而， 三维图像仅能在超声成像仪内部观察，无法实时地输出到外部计算机上用于图像处理或识 另IJ。因此，目前超声引导介入手术中，对介入工具的引导和识别还是通过二维超声图像。 二维超声图像引导穿刺的关键是定位穿刺针在组织内的位置，即在穿刺过程中， 超声探头通过反复扫描穿刺针及周围组织来确认穿刺针的位置。在二维超声图像显示穿刺 针的时候，由于超声图像是二维平面，而穿刺针相当于空间一条线段。因此，为了更直观地 显示穿刺针，需要移动超声探头使穿刺针位于二维超声成像平面内。为实现这一个目的，超 声探头扫描穿刺针的时候有个技巧，即利用穿刺针针尖在超声图像上亮度较高的特点，首 先使穿刺针针尖落入超声成像平面内，保持针尖始终显示的同时，适当移动探头使整个穿 刺针位于成像平面内。目前该方法存在的主要问题有： 一是需要人工判断穿刺针是否在成像平面内。由于超声成像平面存在一定厚度， 同时受到伪像噪声等影响，实现穿刺针与成像平面精确对齐需要精细的操作和丰富超声图 像识读经验。 二是需要人工判断穿刺针在超声图像上的位置。尽管针尖在超声图像上显示为一 个高亮区域，然而针体的亮度低于针尖亮度，周围组织中的亮斑容易湮没部分针体，导致穿 刺针在图像上的位置难于判定，因此，操作人员的主观判断将影响穿刺针图像坐标定位的 正确性。
技术实现思路
本专利技术目的在于提供，能够在超声探头 扫描穿刺针过程中自动判断穿刺针是否位于成像平面内，自动定位针尖图像位置。 实现本专利技术目的技术方案： -种二维超声图像与穿刺针共面判断方法，手持超声探头扫描组织内穿刺针，保 存组织内穿刺针的超声图像，其特征在于： 步骤1 :提取超声图像上的成像区域，将彩色图像转换为灰度图像，经阈值变换后 转换成二值图像； 步骤2 :对步骤1获得的二值图像进行边缘检测，通过Radon变换确定穿刺针所在 的直线； 步骤3 :定位穿刺针针尖所在的位置； 步骤4 :通过判断穿刺针所在直线与超声成像区域边界交点的灰度值是否大于某 一阈值，确定穿刺针是否位于成像平面内。 步骤2中，使用Sobel算子对步骤1获得的二值图像进行边缘检测，具体方法如 下， 首先通过公式（2)和（3)计算预处理后的二值图像f(x，y)在X轴和Y轴方向的 卷积分量A、B，然后利用公式（4)选取A、B的最大值作为检测算子S (X，y)，最后选取阈值 THs，若S(x，y)>Ths，将图像中的对应点设置为1，其他点设置为0,从而得到该图像的所有 边缘， A = [f (x+1, y-1) +2*f (x+1, y) +f (x+1, y+1) ] ~[f (x~l, y~l) +2*f (x-1, y) +f (x~l, y+1)] (2) B = [f (x-1, y-l)+2f (x, y-l)+f (x+1, y-l)]-[f (x~l, y+l)+2*f (x, y+l)+f (x+1, y+1)] (3) S (x, y) = max (A, B) (4) 步骤2中，通过Radon变换确定穿刺针所在的直线，具体方法如下， 从0°到179°每隔Γ利用公式（5)做一次radon变换。 R( ρ，Θ ) = / / Df(x，y) δ ( p -xcos Θ -ysin Θ ) dxdy (5) 在图像平面D上以图像中心为坐标原点，选取某一条积分直线，P为坐标原点到 积分直线的距离，Θ为上述距离所在直线与X轴夹角，直线上点（ X，y)的灰度值为f(X，y)， δ为Dirac delta函数，它使f (X，y)沿直线p = xcos Θ+ySin Θ进行积分，得到的像 R( p，Θ )就是函数f (X，y)的Radon变换，p - Θ空间中的任意一个点就会对应图像空间的 一条直线； 在P-θ空间中定位最亮点即为图像空间中穿刺针所在的直线，然后将p = xcos Θ +ysin Θ利用公式（6)和（7)转换为最终的直线方程y = mx+n, cos Or *-) m =--C6) sin (^- *-) 180 n = -1- (7) _ sinU* A) 180 步骤3中，沿穿刺针所在的直线，遍历直线上所有点的灰度值，选择其中灰度值最 大的像素点即为针尖位置，其图像坐标即针尖坐标。 步骤4中，如果大于该阈值，则认为穿刺针位于成像平面内；如果小于该阈值，则 认为穿刺针没有完全位于成像平面内。 本专利技术具有的有益效果： 本专利技术仅利用二维超声探头扫描已刺入肌肉内的穿刺针，将显示穿刺针的超声图 像经过一系列预处理后形成二值图像，经过边缘检测后，运用Radon变换检测直线来识别 穿刺针位置，并通过比较被识别出直线的最大亮度点来定位针尖图像位置，再通过判断针 尾是否在成像平面来判断穿刺针是否完全位于超声成像平面内，以此引导超声探头移动使 穿刺针为于成像平面内，从而更好地显示穿刺针在组织中的位置，辅助医生完成超声图像 引导的肌肉穿刺过程。本专利技术能够在超声探头扫描穿刺针过程中自动判断穿刺针是否位于 成像平面内，自动定位针尖图像位置，无需操作经验和超声图像识读知识，能够快速实现二 维超声探头追踪定位穿刺针。 本专利技术与通过改造超声探头方法相比较，本专利技术无需改变超声扫描仪的设计，仅 使用一般二维超声扫描仪就可简便、准确地完成穿刺针的图像定位。本专利技术与通过改变穿 刺针结构或加装元件方法相比较，本专利技术针对普通穿刺针在组织内跟踪定位，利用图像识 别方法就可在超声探头扫描穿刺针过程中判断穿刺针与成像平面的位置关系，能够不依赖 操作者的经验。 本专利技术与使用其他辅助定位方法相比较，本专利技术仅利用超声图像识别算法自动定 位穿刺针位置，给出本文档来自技高网...

【技术保护点】
一种二维超声图像与穿刺针共面判断方法，手持超声探头扫描组织内穿刺针，保存组织内穿刺针的超声图像，其特征在于: 步骤1：提取超声图像上的成像区域，将彩色图像转换为灰度图像，经阈值变换后转换成二值图像； 步骤2：对步骤1获得的二值图像进行边缘检测，通过Radon变换确定穿刺针所在的直线； 步骤3：定位穿刺针针尖所在的位置； 步骤4：通过判断穿刺针所在直线与超声成像区域边界交点的灰度值是否大于某一阈值，确定穿刺针是否位于成像平面内。

【技术特征摘要】
1. 一种二维超声图像与穿刺针共面判断方法，手持超声探头扫描组织内穿刺针，保存 组织内穿刺针的超声图像，其特征在于： 步骤1 :提取超声图像上的成像区域，将彩色图像转换为灰度图像，经阈值变换后转换 成二值图像； 步骤2 :对步骤1获得的二值图像进行边缘检测，通过Radon变换确定穿刺针所在的直 线. 步骤3 :定位穿刺针针尖所在的位置； 步骤4 :通过判断穿刺针所在直线与超声成像区域边界交点的灰度值是否大于某一阈 值，确定穿刺针是否位于成像平面内。2. 根据权利要求1所述的二维超声图像与穿刺针共面判断方法，其特征在于：步骤2 中， 使用Sobel算子对步骤1获得的二值图像进行边缘检测，具体方法如下， 首先通过公式（2)和（3)计算预处理后的二值图像f(x，y)在X轴和Y轴方向的卷积 分量A、B，然后利用公式（4)选取A、B的最大值作为检测算子S(x，y)，最后选取阈值THs， 若S(x，y)>Ths，将图像中的对应点设置为1，其他点设置为0,从而得到该图像的所有边缘，3. 根据权利要求2所述的二维超声图像与穿刺针共面判断方法，其特征在于：步骤2 中，通过Radon变换确定穿刺针所在的直线，具...

【专利技术属性】
技术研发人员：栾宽，李金，丛珊，李海波，孙小越，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙江;23