CN104248454B - 一种二维超声图像与穿刺针共面判断方法 - Google Patents
一种二维超声图像与穿刺针共面判断方法 Download PDFInfo
- Publication number
- CN104248454B CN104248454B CN201410525601.4A CN201410525601A CN104248454B CN 104248454 B CN104248454 B CN 104248454B CN 201410525601 A CN201410525601 A CN 201410525601A CN 104248454 B CN104248454 B CN 104248454B
- Authority
- CN
- China
- Prior art keywords
- puncture needle
- image
- straight line
- place
- coplanar
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000003384 imaging method Methods 0.000 claims abstract description 52
- 238000006243 chemical reaction Methods 0.000 claims abstract description 15
- 229910052704 radon Inorganic materials 0.000 claims abstract description 15
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims abstract description 15
- 238000001514 detection method Methods 0.000 claims abstract description 9
- 239000000284 extract Substances 0.000 claims abstract description 7
- 230000009466 transformation Effects 0.000 claims abstract description 4
- 239000000523 sample Substances 0.000 abstract description 18
- 238000002604 ultrasonography Methods 0.000 abstract description 14
- 210000001519 tissue Anatomy 0.000 description 9
- 238000005481 NMR spectroscopy Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 3
- 238000002675 image-guided surgery Methods 0.000 description 2
- 238000007918 intramuscular administration Methods 0.000 description 2
- 238000002372 labelling Methods 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 101100117236 Drosophila melanogaster speck gene Proteins 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 238000009432 framing Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明涉及一种二维超声图像与穿刺针共面判断方法,手持超声探头扫描组织内穿刺针,保存组织内穿刺针的超声图像,其特征在于:步骤1:提取超声图像上的成像区域,将彩色图像转换为灰度图像,经阈值变换后转换成二值图像;步骤2:对步骤1获得的二值图像进行边缘检测,通过Radon变换确定穿刺针所在的直线;步骤3:定位穿刺针针尖所在的位置;步骤4:通过判断穿刺针所在直线与超声成像区域边界交点的灰度值是否大于某一阈值,确定穿刺针是否位于成像平面内。
Description
技术领域
本发明涉及一种二维超声图像与穿刺针共面判断方法。
背景技术
图像引导外科系统利用医学影像和计算机图形图像技术,可在术前对患者多模态图像数据进行三维重建和可视化,获得三维模型,制定合理、定量的手术计划,开展术前模拟;在术中利用三维空间定位设备进行图像和病人物理空间的配准,把患者的实际体位、手术器械的实时空间位置映射到患者的三维图像空间,对手术器械在空间中的位置实时采集并显示,医生通过观察三维图像中手术器械与病变部位的相对位置关系,对病人进行精确的手术治疗。这种新兴外科手术不仅能有效地解决病症,而且还有不开刀、创伤小、恢复快、效果好优点,已经广泛应用于临床。虽然图像引导外科手术有诸多优点,但手术中引导手术工具到达预定目标组织的跟踪定位仍然是其应用的技术难点。目前,术中图像引导方法有computedtomography(CT),X光透视成像,核磁共振成像和超声成像。与CT和核磁共振等成像方式相比较,超声成像具有廉价、实时、无损伤、无辐射和便携等优势,已经作为一种实时的术中成像技术被引入,用于软组织及血流成像。特别是近年来发展的图像引导介入手术,超声成像凸显其直观、实时、易于识别等特点,更方便外科医生观察目标部位的解剖特征,判断术具与目标部位的相对位置关系。目前超声成像系统成像模式包括二维成像和三维成像。作为最新的成像方法,三维成像虽然能够实时地显示人体内部组织的形态,然而,三维图像仅能在超声成像仪内部观察,无法实时地输出到外部计算机上用于图像处理或识别。因此,目前超声引导介入手术中,对介入工具的引导和识别还是通过二维超声图像。
二维超声图像引导穿刺的关键是定位穿刺针在组织内的位置,即在穿刺过程中,超声探头通过反复扫描穿刺针及周围组织来确认穿刺针的位置。在二维超声图像显示穿刺针的时候,由于超声图像是二维平面,而穿刺针相当于空间一条线段。因此,为了更直观地显示穿刺针,需要移动超声探头使穿刺针位于二维超声成像平面内。为实现这一个目的,超声探头扫描穿刺针的时候有个技巧,即利用穿刺针针尖在超声图像上亮度较高的特点,首先使穿刺针针尖落入超声成像平面内,保持针尖始终显示的同时,适当移动探头使整个穿刺针位于成像平面内。目前该方法存在的主要问题有:
一是需要人工判断穿刺针是否在成像平面内。由于超声成像平面存在一定厚度,同时受到伪像噪声等影响,实现穿刺针与成像平面精确对齐需要精细的操作和丰富超声图像识读经验。
二是需要人工判断穿刺针在超声图像上的位置。尽管针尖在超声图像上显示为一个高亮区域,然而针体的亮度低于针尖亮度,周围组织中的亮斑容易湮没部分针体,导致穿刺针在图像上的位置难于判定,因此,操作人员的主观判断将影响穿刺针图像坐标定位的正确性。
发明内容
本发明目的在于提供一种二维超声图像与穿刺针共面判断方法,能够在超声探头扫描穿刺针过程中自动判断穿刺针是否位于成像平面内,自动定位针尖图像位置。
实现本发明目的技术方案:
一种二维超声图像与穿刺针共面判断方法,手持超声探头扫描组织内穿刺针,保存组织内穿刺针的超声图像,其特征在于:
步骤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-1,y-1)+2*f(x-1,y)+f(x-1,y+1)]
(2)
B=[f(x-1,y-1)+2f(x,y-1)+f(x+1,y-1)]-[f(x-1,y+1)+2*f(x,y+1)+f(x+1,y+1)]
(3)
S(x,y)=max(A,B)
(4)
步骤2中,通过Radon变换确定穿刺针所在的直线,具体方法如下,
从0°到179°每隔1°利用公式(5)做一次radon变换。
R(ρ,θ)=∫∫Df(x,y)δ(ρ-xcosθ-ysinθ)dxdy(5)
在图像平面D上以图像中心为坐标原点,选取某一条积分直线,ρ为坐标原点到积分直线的距离,θ为上述距离所在直线与x轴夹角,直线上点(x,y)的灰度值为f(x,y),δ为Diracdelta函数,它使f(x,y)沿直线ρ=xcosθ+ysinθ进行积分,得到的像R(ρ,θ)就是函数f(x,y)的Radon变换,ρ-θ空间中的任意一个点就会对应图像空间的一条直线;
在ρ-θ空间中定位最亮点即为图像空间中穿刺针所在的直线,然后将ρ=xcosθ+ysinθ利用公式(6)和(7)转换为最终的直线方程y=mx+n,
步骤3中,沿穿刺针所在的直线,遍历直线上所有点的灰度值,选择其中灰度值最大的像素点即为针尖位置,其图像坐标即针尖坐标。
步骤4中,如果大于该阈值,则认为穿刺针位于成像平面内;如果小于该阈值,则认为穿刺针没有完全位于成像平面内。
本发明具有的有益效果:
本发明仅利用二维超声探头扫描已刺入肌肉内的穿刺针,将显示穿刺针的超声图像经过一系列预处理后形成二值图像,经过边缘检测后,运用Radon变换检测直线来识别穿刺针位置,并通过比较被识别出直线的最大亮度点来定位针尖图像位置,再通过判断针尾是否在成像平面来判断穿刺针是否完全位于超声成像平面内,以此引导超声探头移动使穿刺针为于成像平面内,从而更好地显示穿刺针在组织中的位置,辅助医生完成超声图像引导的肌肉穿刺过程。本发明能够在超声探头扫描穿刺针过程中自动判断穿刺针是否位于成像平面内,自动定位针尖图像位置,无需操作经验和超声图像识读知识,能够快速实现二维超声探头追踪定位穿刺针。
本发明与通过改造超声探头方法相比较,本发明无需改变超声扫描仪的设计,仅使用一般二维超声扫描仪就可简便、准确地完成穿刺针的图像定位。本发明与通过改变穿刺针结构或加装元件方法相比较,本发明针对普通穿刺针在组织内跟踪定位,利用图像识别方法就可在超声探头扫描穿刺针过程中判断穿刺针与成像平面的位置关系,能够不依赖操作者的经验。
本发明与使用其他辅助定位方法相比较,本发明仅利用超声图像识别算法自动定位穿刺针位置,给出针尖坐标,无需专业图像识读知识就可实现精确定位。
附图说明
1图为本发明方法流程图。
具体实施方式
如图1所示,本发明二维超声图像与穿刺针共面判断方法包括如下步骤:
步骤1:图像预处理。提取超声图像上的成像区域,将彩色图像转换为灰度图像,经阈值变换后转换成二值图像;
提取超声图像上的成像区域:如扇形探头则提取扇形区域,矩形探头则提取矩形区域。这些区域的坐标由超声成像仪的参数设置决定。
将彩色图像转换为灰度图像:保存的超声图像尽管显示的是灰度物体,但往往图像格式是按照彩色图像存储的,不便于进一步处理,因此需要将提取成像区域后的图像由彩色格式转换为灰度格式。具体是将像素的RGB三个分量通过公式(1)计算其对应的灰度值。
Gray=R*0.299+G*0.587+B*0.114(1)
生成二值图像:对转换后的灰度图像选取适当的阈值,使得大于等于该阈值的像素点置为1,小于该阈值的像素点置为0。
步骤2:对步骤1获得的二值图像进行边缘检测,通过Radon变换确定穿刺针所在的直线;
使用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-1,y-1)+2*f(x-1,y)+f(x-1,y+1)]
(2)
B=[f(x-1,y-1)+2f(x,y-1)+f(x+1,y-1)]-[f(x-1,y+1)+2*f(x,y+1)+f(x+1,y+1)]
(3)
S(x,y)=max(A,B)
(4)
通过Radon变换确定穿刺针所在的直线:
从0°到179°每隔1°利用公式(5)做一次radon变换。
R(ρ,θ)=∫∫Df(x,y)δ(ρ-xcosθ-ysinθ)dxdy(5)
在图像平面D上以图像中心为坐标原点,选取某一条积分直线,ρ为坐标原点到积分直线的距离,θ为上述距离所在直线与x轴夹角,直线上点(x,y)的灰度值为f(x,y),δ为Diracdelta函数,它使f(x,y)沿直线ρ=xcosθ+ysinθ进行积分,得到的像R(ρ,θ)就是函数f(x,y)的Radon变换,ρ-θ空间中的任意一个点就会对应图像空间的一条直线;
在ρ-θ空间中定位最亮点即为图像空间中穿刺针所在的直线,然后将ρ=xcosθ+ysinθ利用公式(6)和(7)转换为最终的直线方程y=mx+n,
步骤3:定位穿刺针针尖所在的位置;
依据穿刺针成像特点,针尖通常亮度较高,易于在组织中发现。因此,步骤2中被识别出的直线上亮度最高的点即为穿刺针针尖。具体方法为沿该直线遍历直线上所有点的灰度值,选择其中灰度值最大的像素点即为针尖位置,其图像坐标也就是针尖坐标。
步骤4:判断穿刺针是否位于成像平面内。
由于两点确定一条直线,针尖已在成像平面内定位,如果图像上也能定位针尾,那么就可以判定穿刺针落入成像平面内。具体方法是判断穿刺针所在直线与超声成像区域(扇形或矩形)边界交点的灰度值是否大于某一阈值。如果大于该阈值,则认为穿刺针位于成像平面内;如果小于该阈值,则认为穿刺针没有完全位于成像平面内,即穿刺针与成像平面成一定角度。
具体实施本发明方法时:
硬件配置环境:包括:(1)一台二维超声扫描仪,配备一个二维超声探头;(2)一根穿刺针,用于被跟踪定位;(3)一台普通计算机,用于处理超声图像。如超声扫描仪不能直接实时输出图像,则需要在计算机上安装视频采集卡。
获取图像:超声探头扫描组织内的穿刺针时,通过移动超声探头借助亮度信息先找到穿刺针尖,保持针尖始终显示并调整探头角度,使得超声图像尽可能完整地显示穿刺针,在此过程中计算机实时获取超声图像。
定位穿刺针:在获取的某一幅穿刺针图像上,利用本发明的图像定位方法判断穿刺针在图像上的位置及针尖的图像坐标,同时依据针尾亮度判断穿刺针是否完全位于成像平面内。
标记超声图像上穿刺针位置:如果穿刺针位于成像平面内,则保存步骤3中得到的针尖及针尾点坐标,并转化为图像坐标。因为Radon变换后得到直线所在的坐标系是以图像中心为坐标原点、向右为X轴正方向、向上为Y轴正方向的坐标系;而原图像的坐标系是以左上角的顶点为坐标原点、向右为X轴正方向、向下为Y轴正方向的坐标系,因此需要最终转化为图像坐标系下的坐标。最后在超声图像上针尖与针尾间画出一段彩色线段,以标记穿刺针位置。
Claims (5)
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,从而得到该图像的所有边缘,
A=[f(x+1,y-1)+2*f(x+1,y)+f(x+1,y+1)]-[f(x-1,y-1)+2*f(x-1,y)+f(x-1,y+1)](2)
B=[f(x-1,y-1)+2f(x,y-1)+f(x+1,y-1)]-[f(x-1,y+1)+2*f(x,y+1)+f(x+1,y+1)](3)
S(x,y)=max(A,B)(4)。
3.根据权利要求2所述的二维超声图像与穿刺针共面判断方法,其特征在于:步骤2中,通过Radon变换确定穿刺针所在的直线,具体方法如下,
从0°到179°每隔1°利用公式(5)做一次radon变换,
R(ρ,θ)=∫∫Df(x,y)δ(ρ-xcosθ-ysinθ)dxdy(5)
在图像平面D上以图像中心为坐标原点,选取某一条积分直线,ρ为坐标原点到积分直线的距离,θ为上述距离所在直线与x轴夹角,直线上点(x,y)的灰度值为f(x,y),δ为Diracdelta函数,它使f(x,y)沿直线ρ=xcosθ+ysinθ进行积分,得到的像R(ρ,θ)就是函数f(x,y)的Radon变换,ρ-θ空间中的任意一个点就会对应图像空间的一条直线;
在ρ-θ空间中定位最亮点即为图像空间中穿刺针所在的直线,然后将ρ=xcosθ+ysinθ利用公式(6)和(7)转换为最终的直线方程y=mx+n,
4.根据权利要求3所述的二维超声图像与穿刺针共面判断方法,其特征在于:步骤3中,沿穿刺针所在的直线,遍历直线上所有点的灰度值,选择其中灰度值最大的像素点即为针尖位置,其图像坐标即针尖坐标。
5.根据权利要求4所述的二维超声图像与穿刺针共面判断方法,其特征在于:步骤4中,如果大于该阈值,则认为穿刺针位于成像平面内;如果小于该阈值,则认为穿刺针没有完全位于成像平面内。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410525601.4A CN104248454B (zh) | 2014-10-09 | 2014-10-09 | 一种二维超声图像与穿刺针共面判断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410525601.4A CN104248454B (zh) | 2014-10-09 | 2014-10-09 | 一种二维超声图像与穿刺针共面判断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104248454A CN104248454A (zh) | 2014-12-31 |
CN104248454B true CN104248454B (zh) | 2016-06-29 |
Family
ID=52183939
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410525601.4A Expired - Fee Related CN104248454B (zh) | 2014-10-09 | 2014-10-09 | 一种二维超声图像与穿刺针共面判断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104248454B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105286960B (zh) * | 2015-12-03 | 2019-01-18 | 上海联影医疗科技有限公司 | 乳腺病灶穿刺位置的定位方法及系统 |
CN106344153B (zh) * | 2016-08-23 | 2019-04-02 | 深圳先进技术研究院 | 一种柔性穿刺针针尖自动跟踪装置及方法 |
CN108335326B (zh) * | 2018-01-31 | 2022-07-26 | 上海联影医疗科技股份有限公司 | 穿刺系统的校正方法、装置及系统 |
US11638569B2 (en) | 2018-06-08 | 2023-05-02 | Rutgers, The State University Of New Jersey | Computer vision systems and methods for real-time needle detection, enhancement and localization in ultrasound |
EP3836843A4 (en) | 2018-08-13 | 2022-02-16 | Rutgers, the State University of New Jersey | ARTIFICIAL VISION SYSTEMS AND METHODS FOR REAL-TIME LOCATION OF NEEDLES IN ULTRASOUND IMAGES |
CN111329585A (zh) * | 2018-12-18 | 2020-06-26 | 深圳迈瑞生物医疗电子股份有限公司 | 超声图像处理方法及超声成像设备 |
CN109907801B (zh) * | 2019-03-08 | 2021-09-10 | 哈尔滨工程大学 | 一种可定位超声引导穿刺方法 |
CN110276772B (zh) * | 2019-05-10 | 2022-10-18 | 深圳大学 | 一种肌肉组织中结构要素的自动定位方法及系统 |
CN111476790A (zh) * | 2020-05-14 | 2020-07-31 | 武汉中旗生物医疗电子有限公司 | 一种超声穿刺中穿刺针增强显示方法及装置 |
CN113592057B (zh) * | 2021-07-22 | 2024-07-30 | 广东环凯生物科技有限公司 | 菌落计数分析方法、装置及可读存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3750972B2 (ja) * | 1998-11-06 | 2006-03-01 | 株式会社東芝 | 3次元超音波診断装置 |
US6438401B1 (en) * | 2000-04-28 | 2002-08-20 | Alpha Intervention Technology, Inc. | Indentification and quantification of needle displacement departures from treatment plan |
JP4388255B2 (ja) * | 2002-05-21 | 2009-12-24 | アロカ株式会社 | 穿刺用超音波探触子 |
US8556815B2 (en) * | 2009-05-20 | 2013-10-15 | Laurent Pelissier | Freehand ultrasound imaging systems and methods for guiding fine elongate instruments |
-
2014
- 2014-10-09 CN CN201410525601.4A patent/CN104248454B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN104248454A (zh) | 2014-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104248454B (zh) | 一种二维超声图像与穿刺针共面判断方法 | |
US11707256B2 (en) | System and method for tracking completeness of co-registered medical image data | |
US11896414B2 (en) | System and method for pose estimation of an imaging device and for determining the location of a medical device with respect to a target | |
CN103371870B (zh) | 一种基于多模影像的外科手术导航系统 | |
RU2748435C2 (ru) | Ультразвуковая система и способ для визуализации ткани груди | |
WO2015161728A1 (zh) | 三维模型的构建方法及装置、图像监控方法及装置 | |
EP2044884B1 (de) | Bestimmung und Erkennung von Lageänderungen von Körperstrukturteilen | |
CN114129240B (zh) | 一种引导信息生成方法、系统、装置及电子设备 | |
CN105678746B (zh) | 一种医学图像中肝脏范围的定位方法及装置 | |
JP5513790B2 (ja) | 超音波診断装置 | |
JP2008073304A (ja) | 超音波乳房診断システム | |
JP2008073305A (ja) | 超音波乳房診断システム | |
CN107835661B (zh) | 超声图像处理系统和方法及其装置、超声诊断装置 | |
CN110400286A (zh) | 一种x射线ct图像中金属针的检测定位方法 | |
CN106157282A (zh) | 图像处理系统及方法 | |
CN103584885B (zh) | 一种基于定位导航穿刺针的自由臂超声标定方法 | |
US10617383B2 (en) | Medical image display apparatus, method, and medium | |
CN106068098B (zh) | 用于超声引导过程的区域可视化 | |
CN105894508A (zh) | 一种医学图像自动定位质量的评估方法 | |
Wijata et al. | Unbiased validation of the algorithms for automatic needle localization in ultrasound-guided breast biopsies | |
JP2018043000A (ja) | Ent画像の位置合わせ | |
JP6875774B1 (ja) | 医療画像用のマーキングシステム及びマーキング支援装置 | |
WO2019130636A1 (ja) | 超音波撮像装置、画像処理装置、及び方法 | |
US20240307030A1 (en) | Medical support device, and operation method and operation program of medical support device | |
Doerfler et al. | Blood vessel detection in navigated ultrasound: An assistance system for liver resections |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160629 |