CN110664438A - 超声扫查轨迹规划方法、装置、存储介质及计算机设备 - Google Patents
超声扫查轨迹规划方法、装置、存储介质及计算机设备 Download PDFInfo
- Publication number
- CN110664438A CN110664438A CN201911007896.5A CN201911007896A CN110664438A CN 110664438 A CN110664438 A CN 110664438A CN 201911007896 A CN201911007896 A CN 201911007896A CN 110664438 A CN110664438 A CN 110664438A
- Authority
- CN
- China
- Prior art keywords
- curve
- track
- points
- point cloud
- scanning
- 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
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0825—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the breast, e.g. mammography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种超声扫查轨迹规划方法,其包括:获取胸部区域的三维点云数据;根据三维点云数据对胸部区域结构进行骨架模型重建,以得到曲线骨架;按照预设的曲线分割条件对曲线骨架中的各条曲线进行分割,并取各条曲线上的所有分割点;根据预设的超声扫查方向从分割点集合中选取多组分割点,并将每一组分割点连接成一条扫查轨迹曲线;从扫查轨迹曲线中提取多个轨迹点,并计算各轨迹点的姿态角。根据胸部区域的三维点云数据计算生成的扫查轨迹,以使超声探头能够根据接触区域的形状调整扫查姿态,保证获取到的每一帧超声图像所涵盖的信息全面、准确。
Description
技术领域
本发明涉及超声扫查技术领域,尤其涉及一种超声扫查轨迹规划方法、装置、存储介质及计算设备。
背景技术
超声超声波检查是超声波检查的一种方式,这种非手术的诊断性检查,对受检者无痛苦、无损伤、无放射性。并且,超声可以清晰地显示内脏器官及器官周围的各种断面图像,由于图像富于实体感,接近于解剖的真实结构,所以应用超声检查可以早期明确诊断。基于超声检查的诸多优势,从专业的医疗疾病诊断到日常化的健康指标评估,其应用范围越来越广。
随着医疗诊断技术的发展以及经济水平的提高,越来越多女性开始关注乳腺健康,而通过超声可以方便快捷地对乳腺生理状况进行初步判断,因此不管是公立医疗机构,还是盈利性健康服务机构,推出了多种乳腺筛查服务,以满足女性用户的需求。但是,现有的乳腺筛查模式,基本上是借助于常规的超声设备,通过操作人员手持超声探头进行乳房区域的扫查,在扫查过程中,探头的扫查轨迹通常是操作人员的主观选择,可能存在未扫描到的位置,并且探头姿态也难以根据乳腺对应区域的乳房形状进行适应性调整,还可能存在超声图像信息缺失的现象,总的来说,仅依靠操作人员的主观操作来实现对乳房的全面扫查,存在诸多弊端,不利于对乳腺生理状况做出准确判断。
发明内容
本发明的主要目的在于提供一种超声扫查轨迹规划方法,旨在解决现有的乳房超声扫查方式难以获取全面、准确的超声图像的技术问题。应用于使用B型超声探头进行乳房自动扫查的设备中。
为实现上述目的,本发明提供一种超声扫查轨迹规划方法,包括:
获取胸部区域的三维点云数据;
根据所述三维点云数据对胸部区域结构进行骨架模型重建,以得到曲线骨架;
按照预设的曲线分割条件对所述曲线骨架中的各条曲线进行分割,并取各条曲线上的所有分割点,所述分割点表示为{Sij,0≤i<A,0≤j<Bi},其中A是所述曲线骨架中曲线的条数,Bi是第i条曲线上的分割点数;
根据预设的超声扫查方向从分割点集合中选取多组分割点,并将每一组分割点连接成一条扫查轨迹曲线;
从所述扫查轨迹曲线中提取多个轨迹点,并计算各轨迹点的姿态角。
优选地,所述根据所述三维点云数据对胸部区域结构进行骨架模型重建,以得到曲线骨架包括:
通过对所述三维点云数据进行横向切片以得到若干段子点云;
使用贝塞尔曲线对每一段所述子点云进行曲线拟合以得到所述曲线骨架。
优选地,所述计算各轨迹点的姿态角包括:
提取所述轨迹点的邻域点集,并通过对所述邻域点集求PCA得到所述轨迹点所在邻域曲面的指向身体外部的单位法向量Vz,并按照以下公式计算所述轨迹点在XY轴的单位方向向量Vx、Vy,其中:
Vy=Vz×[0 0 1]T
Vx=Vy×Vz
将所述轨迹点XYZ坐标轴的单位方向向量转换为欧拉角的表示形式,并提取所述姿态角。
优选地,在所述获取胸部区域的三维点云数据的步骤之后,所述方法还包括:
对所述三维点云数据进行预处理,所述预处理包括点云降采样、点云滤波和点云平滑。
优选地,在所述计算各轨迹点的姿态角的步骤之后,所述方法还包括:
导入超声探头末端的行程极限数据,并根据所述行程极限数据从所述轨迹点中滤除所述超声探头末端不可到达的点。
优选地,在所述导入超声探头末端的行程极限数据,并根据所述行程极限数据从所述轨迹点中滤除所述超声探头末端不可到达的点的步骤之后,所述方法还包括:
对每条扫查轨迹曲线进行平滑滤波。
本发明还提出一种超声扫查轨迹规划装置,包括:
点云获取模块,用于获取胸部区域的三维点云数据;
骨架模型重建模块,用于根据所述三维点云数据对胸部区域结构进行骨架模型重建,以得到曲线骨架;
曲线分割模块,用于按照预设的曲线分割条件对所述曲线骨架中的各条曲线进行分割,并取各条曲线上的所有分割点,所述分割点表示为{Sij,0≤i<A,0≤j<Bi},其中A是所述曲线骨架中曲线的条数,Bi是第i条曲线上的分割点数;
轨迹生成模块,用于根据预设的超声扫查方向从分割点集合中选取多组分割点,并将每一组分割点连接成一条扫查轨迹曲线;
姿态角计算模块,用于从所述扫查轨迹曲线中提取多个轨迹点,并计算各轨迹点的姿态角。
本发明还提出一种计算机程序存储介质,所述计算机程序存储介质中存储有计算机程序代码,该计算机程序代码被处理器执行时实现上述超声扫查轨迹规划方法的步骤。
本发明还提供一种计算机设备,包括处理器、存储器和存储在所述存储器中的计算机程序代码,所述处理器在调用所述计算机程序代码时,实现上述超声扫查轨迹规划方法的步骤。
相较于现有技术,本发明通过采用全自动机械化的扫查方式来对用户乳房区域进行超声扫查,从而对乳腺及其周边器官、组织的生理状况进行全面、准确的判断,以避免因人为操作而导致扫查得到的超声图像覆盖面不全、信息缺失等情况的发生;并且,该全自动机械化的扫查方式在执行超声扫查的过程中,依赖于计算机根据胸部区域的三维点云数据计算生成的扫查轨迹,以使超声探头能够根据接触区域的形状调整扫查姿态,保证获取到的每一帧超声图像所涵盖的信息全面、准确。
附图说明
图1为本发明公开的多个实施例可以在其中实施的示例环境的结构示意图;
图2为本发明的超声扫查轨迹规划方法一实施例的流程示意图;
图3为本发明的超声扫查轨迹规划方法另一实施例的流程示意图;
图4为本发明的超声扫查轨迹规划方法又一实施例的流程示意图;
图5为本发明公开的多个实施例中原始三维点云的示意图;
图6为本发明公开的多个实施例中经过骨架模型重建得到的曲线骨架示意图;
图7为本发明的超声扫查轨迹规划装置一实施例的功能模块示意图;
图8为本发明公开的多个实施例能够在其中实施的计算机设备的结构示意图。
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明提供一种超声扫查轨迹规划方法,其目的在于根据每位用户的胸部区域情况制定最佳的超声扫查轨迹,以使乳房超声扫查设备能够根据该超声扫查轨迹执行扫查过程,进而生成包含信息全面、准确的超声图像。
如图1所示,为了给该方法的实施提供环境基础,前述乳房超声扫查设备主要包括扫查执行机构10、躺床20和点云获取装置30,其中扫查执行机构10包括主机11、与主机11连接的机械臂12以及安装在该机械臂12的执行末端处的超声探头13,本实施例中,主机11具有能够实现通信、数据处理和运动控制功能的相应硬件,主机11还具有用于安装机械臂12的基础架构,比如机械臂12被构造成能够提供三个直线运动自由度和两个以上旋转自由度的多轴结构,从而保证超声探头13能够根据待扫查区域的表面形状作适应性的姿态变换,具体应用时,机械臂12可以是五轴机械臂,或者是六轴机械臂。躺床20可以是固定式的支撑结构,也可以设置成能够提供位置调节的活动结构,比如通过设置升降机构,以调节躺床20的支撑面的高度,又比如通过设置水平移动机构,以调节躺床20的支撑面的水平位置,从而在不需要用户挪动身躯的情况下调节用户的初始位置。点云获取装置30设置在躺床20的上方,为了更加全面地获取三维点云数据,可以按照图1所示结构的指引,配置两套点云获取装置30,此示例中是以用户身躯的纵向为参照方向布置点云获取装置30的,在其它实施方案中,以用户身躯的横向为参照方向布置点云获取装置30也是能够满足要求的,本实施例的点云获取装置30可以是结构光传感器,当然也可以是激光雷达。
至此,已经详细介绍了本发明各个实施例的应用环境和相关设备的硬件结构和功能,并且上述乳房超声扫查设备的结构组成仅为基本功能的示例,并不是对乳房超声扫查设备实现方式的限制。下面,将基于上述应用环境和相关设备,详细介绍超声扫查轨迹规划方法的各个实施例。
如图2所示,在一实施例中,本发明提出的超声扫查轨迹规划方法包括:
步骤S10,获取胸部区域的三维点云数据。
用户的胸部区域(针对女性)作为容易受到自身姿势和外力影响而产生形状变化的部位,为了满足前述乳房超声扫查设备的技术要求,在执行全面的扫查动作前,需要对胸部区域进行束形,比如通过穿上具有一定弹性的束胸背心来调整胸部区域的形状,并保持外形的稳定性。因此,针对每一次的超声扫查过程,一般而言,均需要重新采集三维点云数据。在实际应用时,用户先平躺在躺床20上,并根据实际情况调整位置,直至满足三维点云数据采集和超声扫查的要求,然后通过点云获取装置30采集胸部区域的三维点云数据,考虑到获取的原始点云数据覆盖面较广,需要对原始点云数据进行界限过滤,以简化数据的后期处理难度。通过采集胸部区域的三维点云数据,可以准确地描述胸部区域的三维结构,由此通过后期的扫查轨迹规划算法生成符合实际扫查接触面的超声探头13运动轨迹。
进一步地,在一较佳实施例中,该超声扫查轨迹规划方法还包括:
对三维点云数据进行预处理,该预处理包括点云降采样、点云滤波和点云平滑。
该步骤是在获取到三维点云数据之后执行的,通过对三维点云数据进行预处理操作,可以获得更加符合超声扫查应用场景的点云数据,同时降低数据的复杂程度,提高设备的数据处理效率。具体地,输入的点云比较稠密,全部处理的话耗时较长,因此先对输入点云进行降采样,降低点云的密度,加快处理速度。直观上来说,点云降采样就是对原始点云每间隔一定的空间距离取一个点代表其邻域内的其它点,这样就可以得到一个更稀疏的点云,具体的点云降采样设定标准可以根据点云采集装置的数据采集规格和后期数据处理精度选择,在此不作限制。此外,理论上胸部区域的点云应当构成一个平滑连续的曲面,但由于各种原因会存在一些异常点云(如孤立的几个离散点),通过点云滤波就可以滤除这些异常点云,输出一个更高质量的点云供后续步骤使用。滤波后的点云由于传感器的测量误差,会有不平滑的现象,如水浪般的波纹,因此,进一步对点云进行平滑处理,使点云曲面更加平滑。
步骤与S20,根据三维点云数据对胸部区域结构进行骨架模型重建,以得到曲线骨架。
如图5所示,获取到的三维点云数据的数据量较为庞大,需要对其进行模型的重建,在简化数据的同时,满足扫查轨迹规划算法的应用要求。具体地,根据预设方向对点云进行切片,以人体身躯的方向为参照,主要沿身躯的横向和纵向这两个方向进行切片操作,并且在优选的切片约束条件中,以等间距的方式切片,从而获得一段段等宽的子点云,每段子点云的宽度可以根据实际情况灵活调整。作为一种可能的实施方式,超声探头13采用条形扫查的方式,并且该条形扫查的方向沿身躯的纵向,因此沿身躯的横向进行点云切片,这种扫查方式对运动机构的要求较低,并且能够保证超声图像的质量。
如图3所示,前述步骤S20包括:
步骤S201,通过对三维点云数据进行横向切片以得到若干段子点云;
步骤S202,使用贝塞尔曲线对每一段子点云进行曲线拟合以得到曲线骨架。
如图6所示,重建出的曲线骨架是胸部区域结构的更稳定可靠的表示方式,有利于算法的后期处理。本步骤中,关于贝塞尔曲线的拟合操作可参照现有技术中关于这方面的详细说明,在此不作赘述。
步骤S30,按照预设的曲线分割条件对曲线骨架中的各条曲线进行分割,并取各条曲线上的所有分割点。
在该步骤中,以前述选择的纵向条形扫查方式为例,对横向分布的各条曲线进行等弧长分割,并且根据超声探头13的覆盖面大小设定分割间距,从而保证超声探头13在扫查过程中能够覆盖完整的待扫查区域,同时又能减少重合区。在执行曲线分割的环节,所得到的分割点表示为{Sij,0≤i<A,0≤j<Bi},其中A是曲线骨架中曲线的条数,Bi是第i条曲线上的分割点数,i和j均取正整数,通过对点云进行坐标变换,可以得到各个分割点在超声探头13对应的运动坐标系下的XYZ坐标值,关于点云坐标变换的原理,可参照现有技术的详细说明,在此不作赘述。
步骤S40,根据预设的超声扫查方向从分割点集合中选取多组分割点,并将每一组分割点连接成一条扫查轨迹曲线。
在该步骤中,根据预设的超声扫查方向从分割点集合中选取能够组合成扫查轨迹曲线的多组分割点,以纵向的条形扫查为例,最简单的分组方式是,选取曲线骨架中每条曲线上同序号j的分割点为一组,这样即可得到了一条完整的轨迹{S0j,S1j,S2j,…,SAj}。除了上述作为示例的分割点组合方式,还可以采用其它任意适用的方式进行分割点分组。
步骤S50,从扫查轨迹曲线中提取多个轨迹点,并计算各轨迹点的姿态角。
在该步骤中,作为较佳的实施方式,轨迹点即为前述分割点,这种轨迹点的提取方式能够简化数据处理的过程。当然,除了提取前述的分割点,还可以在相邻分割点之间额外提取一个或多个点作为轨迹点,这里需要结合乳房超声扫查设备的机构运动参数,避免造成数据冗余。以提取的分割点作为轨迹点为例,并且匹配五自由度的乳房超声扫查设备,则需要获取每个轨迹点坐标值和相应的姿态角,将轨迹点表示为Pi=[Xi,Yi,Zi,Ri,Pi],这五个量分别表示Pi的XYZ坐标值和Pi的Roll、Pitch姿态角,其中Pi的XYZ坐标值根据前述的点云数据计算可以得到,因此该步骤主要计算轨迹点的两个姿态角。但是,如果提取得到的轨迹点并非前述分割点,那么还需计算这些未知轨迹点的XYZ坐标值。通过确定每个轨迹点的五个坐标量,根据XYZ坐标值可以在运动控制程序中控制超声探头13运动至待扫查区域的具体位置,而根据Roll、Pitch姿态角可以在运动控制程序中控制超声探头13应当调整到哪种角度姿态,以使探头表面与待扫查区域表面紧密贴合。
如图4所示,在一较佳实施例中,前述计算各轨迹点的姿态角主要采用以下算法,具体步骤包括:
步骤S501,提取轨迹点的邻域点集,并通过对邻域点集求PCA得到轨迹点所在邻域曲面的指向身体外部的单位法向量Vz。
在该步骤中,以轨迹点为中心提取的邻域点集可根据所期望的计算精度选择边界半径,这里不对邻域点集的范围设定作出限制,设定好邻域点集的提取范围后,即可通过对邻域点集求PCA得到轨迹点所在邻域曲面的指向身体外部的单位法向量Vz。
步骤S502,按照公式Vy=Vz×[0 0 1]T、Vx=Vy×Vz计算轨迹点在XY轴的单位方向向量Vx、Vy。
步骤S503,将轨迹点XYZ坐标轴的单位方向向量转换为欧拉角的表示形式,并提取姿态角。
在将单位方向向量转换为欧拉角的表示形式后,实际上可以获取到三个方向的姿态角,具体是提取几个姿态角,可以结合超声探头13所能提供的运动自由度,本实施例以提取Roll、Pitch姿态角作为示例。
此外,在计算各轨迹点的姿态角的步骤之后,考虑到一些轨迹点可能位于超声探头13末端运动范围之外,因此需要对各轨迹点进行校验,具体地,该超声扫查轨迹规划方法还包括:
导入超声探头末端的行程极限数据,并根据行程极限数据从轨迹点中滤除超声探头末端不可到达的点。
通常,可标定出超声探头13的运动极限,并以数据表的形式存储起来备用,根据该数据表来对轨迹点进行校验,可以避免设备在扫查过程中出现意外。同时,滤除了一些轨迹点后,对每条扫查轨迹曲线进行平滑滤波,使得超声探头13在扫查过程中动作更加平滑,减少对人体的局部挤压。
由此可见,本发明的超声扫查轨迹规划方法通过采用全自动机械化的扫查方式来对用户乳房区域进行超声扫查,从而对乳腺及其周边器官、组织的生理状况进行全面、准确的判断,以避免因人为操作而导致扫查得到的超声图像覆盖面不全、信息缺失等情况的发生;并且,该全自动机械化的扫查方式在执行超声扫查的过程中,依赖于计算机根据胸部区域的三维点云数据计算生成的扫查轨迹,以使超声探头能够根据接触区域的形状调整扫查姿态,保证获取到的每一帧超声图像所涵盖的信息全面、准确。
此外,本发明还提供一种超声扫查轨迹规划装置,如图7所示,该超声扫查轨迹规划装置包括:
点云获取模块100,用于获取胸部区域的三维点云数据;
骨架模型重建模块200,用于根据三维点云数据对胸部区域结构进行骨架模型重建,以得到曲线骨架;
曲线分割模块300,用于按照预设的曲线分割条件对曲线骨架中的各条曲线进行分割,并取各条曲线上的所有分割点,分割点表示为{Sij,0≤i<A,0≤j<Bi},其中A是曲线骨架中曲线的条数,Bi是第i条曲线上的分割点数;
轨迹生成模块400,用于根据预设的超声扫查方向从分割点集合中选取多组分割点,并将每一组分割点连接成一条扫查轨迹曲线;
姿态角计算模块500,用于从扫查轨迹曲线中提取多个轨迹点,并计算各轨迹点的姿态角。
上述超声扫查轨迹规划装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中,也可以以软件形式存储于服务器中的存储器中,以便于计算机设备调用执行以上各个模块对应的操作。该计算机设备可以为中央处理单元(CPU)、微计算机设备、单片机等。上述各功能模块所起到的工作原理及起到的作用可参见图2~4中所示的超声扫查轨迹规划方法的实现过程,在此不作赘述。
本发明还提供一种计算机程序存储介质,该计算机程序存储介质中存储有计算机程序代码,该计算机程序代码被处理器执行时实现如下步骤:
获取胸部区域的三维点云数据;
根据三维点云数据对胸部区域结构进行骨架模型重建,以得到曲线骨架;
按照预设的曲线分割条件对曲线骨架中的各条曲线进行分割,并取各条曲线上的所有分割点,分割点表示为{Sij,0≤i<A,0≤j<Bi},其中A是曲线骨架中曲线的条数,Bi是第i条曲线上的分割点数;
根据预设的超声扫查方向从分割点集合中选取多组分割点,并将每一组分割点连接成一条扫查轨迹曲线;
从扫查轨迹曲线中提取多个轨迹点,并计算各轨迹点的姿态角。
该计算机程序被处理器执行时还实现了超声扫查轨迹规划方法的其它步骤,具体可参见包括上述图3、4对应实施例在内的说明,在此不作赘述。
本发明还提供了一种计算机设备,如图8所示,该计算机设备包括处理器40、存储器50和存储在存储器50中的计算机程序代码,处理器40在调用该计算机程序代码时,实现上述各实施例中提供的一种超声扫查轨迹规划方法的步骤。
具体地,该计算机设备可为个人计算机或服务器。该计算机设备包括通过系统总线连接的处理器40、存储器50和通信接口(图未示)。其中,处理器40用于提供计算和控制能力,支撑整个计算机设备的运行。存储器50包括非易失性存储介质和内存储器。非易失性存储介质中存储有操作系统和计算机程序,该计算机程序被处理器40执行时以实现一种超声扫查轨迹规划方法。内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。通信接口用于与外部的服务器或终端通过网络连接通信。
以上仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (9)
1.一种超声扫查轨迹规划方法,其特征在于,包括:
获取胸部区域的三维点云数据;
根据所述三维点云数据对胸部区域结构进行骨架模型重建,以得到曲线骨架;
按照预设的曲线分割条件对所述曲线骨架中的各条曲线进行分割,并取各条曲线上的所有分割点,所述分割点表示为{Sij,0≤i<A,0≤j<Bi},其中A是所述曲线骨架中曲线的条数,Bi是第i条曲线上的分割点数;
根据预设的超声扫查方向从分割点集合中选取多组分割点,并将每一组分割点连接成一条扫查轨迹曲线;
从所述扫查轨迹曲线中提取多个轨迹点,并计算各轨迹点的姿态角。
2.根据权利要求1所述的超声扫查轨迹规划方法,其特征在于,所述根据所述三维点云数据对胸部区域结构进行骨架模型重建,以得到曲线骨架包括:
通过对所述三维点云数据进行横向切片以得到若干段子点云;
使用贝塞尔曲线对每一段所述子点云进行曲线拟合以得到所述曲线骨架。
3.根据权利要求1所述的超声扫查轨迹规划方法,其特征在于,所述计算各轨迹点的姿态角包括:
提取所述轨迹点的邻域点集,并通过对所述邻域点集求PCA得到所述轨迹点所在邻域曲面的指向身体外部的单位法向量Vz,并按照以下公式计算所述轨迹点在XY轴的单位方向向量Vx、Vy,其中:
Vy=Vz×[0 0 1]T
Vx=Vy×Vz
将所述轨迹点XYZ坐标轴的单位方向向量转换为欧拉角的表示形式,并提取所述姿态角。
4.根据权利要求1所述的超声扫查轨迹规划方法,其特征在于,在所述获取胸部区域的三维点云数据的步骤之后,所述方法还包括:
对所述三维点云数据进行预处理,所述预处理包括点云降采样、点云滤波和点云平滑。
5.根据权利要求1所述的超声扫查轨迹规划方法,其特征在于,在所述计算各轨迹点的姿态角的步骤之后,所述方法还包括:
导入超声探头末端的行程极限数据,并根据所述行程极限数据从所述轨迹点中滤除所述超声探头末端不可到达的点。
6.根据权利要求5所述的超声扫查轨迹规划方法,其特征在于,在所述导入超声探头末端的行程极限数据,并根据所述行程极限数据从所述轨迹点中滤除所述超声探头末端不可到达的点的步骤之后,所述方法还包括:
对每条扫查轨迹曲线进行平滑滤波。
7.一种超声扫查轨迹规划装置,其特征在于,包括:
点云获取模块,用于获取胸部区域的三维点云数据;
骨架模型重建模块,用于根据所述三维点云数据对胸部区域结构进行骨架模型重建,以得到曲线骨架;
曲线分割模块,用于按照预设的曲线分割条件对所述曲线骨架中的各条曲线进行分割,并取各条曲线上的所有分割点,所述分割点表示为{Sij,0≤i<A,0≤j<Bi},其中A是所述曲线骨架中曲线的条数,Bi是第i条曲线上的分割点数;
轨迹生成模块,用于根据预设的超声扫查方向从分割点集合中选取多组分割点,并将每一组分割点连接成一条扫查轨迹曲线;
姿态角计算模块,用于从所述扫查轨迹曲线中提取多个轨迹点,并计算各轨迹点的姿态角。
8.一种计算机程序存储介质,其特征在于,所述计算机程序存储介质中存储有计算机程序代码,该计算机程序代码被处理器执行时实现权利要求1至6中任一项所述超声扫查轨迹规划方法的步骤。
9.一种计算机设备,包括处理器、存储器和存储在所述存储器中的计算机程序代码,其特征在于,所述处理器在调用所述计算机程序代码时,实现权利要求1至6中任意一项所述超声扫查轨迹规划方法的步骤。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911007896.5A CN110664438B (zh) | 2019-10-22 | 2019-10-22 | 超声扫查轨迹规划方法、装置、存储介质及计算机设备 |
PCT/CN2020/121233 WO2021078064A1 (zh) | 2019-10-22 | 2020-10-15 | 超声扫查轨迹规划方法、装置、存储介质及计算机设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911007896.5A CN110664438B (zh) | 2019-10-22 | 2019-10-22 | 超声扫查轨迹规划方法、装置、存储介质及计算机设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110664438A true CN110664438A (zh) | 2020-01-10 |
CN110664438B CN110664438B (zh) | 2021-09-10 |
Family
ID=69083693
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911007896.5A Active CN110664438B (zh) | 2019-10-22 | 2019-10-22 | 超声扫查轨迹规划方法、装置、存储介质及计算机设备 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN110664438B (zh) |
WO (1) | WO2021078064A1 (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112075957A (zh) * | 2020-07-27 | 2020-12-15 | 深圳瀚维智能医疗科技有限公司 | 乳腺环扫轨迹规划方法、装置及计算机可读存储介质 |
CN112419415A (zh) * | 2020-12-08 | 2021-02-26 | 浙江德尚韵兴医疗科技有限公司 | 一种基于crs曲线拟合实现位姿规划的超声扫查方法 |
WO2021078066A1 (zh) * | 2019-10-22 | 2021-04-29 | 深圳瀚维智能医疗科技有限公司 | 乳腺超声筛查方法、装置及系统 |
WO2021078064A1 (zh) * | 2019-10-22 | 2021-04-29 | 深圳瀚维智能医疗科技有限公司 | 超声扫查轨迹规划方法、装置、存储介质及计算机设备 |
CN112767309A (zh) * | 2020-12-30 | 2021-05-07 | 无锡祥生医疗科技股份有限公司 | 超声扫查方法、超声设备及系统 |
CN112767415A (zh) * | 2021-01-13 | 2021-05-07 | 深圳瀚维智能医疗科技有限公司 | 胸部扫查区域自动确定方法、装置、设备及存储介质 |
CN112767237A (zh) * | 2020-12-30 | 2021-05-07 | 无锡祥生医疗科技股份有限公司 | 基于点云数据的环形位姿控制方法、装置及超声设备 |
CN116851930A (zh) * | 2023-04-28 | 2023-10-10 | 华南智能机器人创新研究院 | 一种三维五轴曲面激光切割轨迹生成方法及装置 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113515513B (zh) * | 2021-06-30 | 2023-04-21 | 同济大学 | 轨迹矫正方法及装置、点云地图生成方法及装置 |
CN113456106B (zh) * | 2021-08-03 | 2023-11-21 | 无锡祥生医疗科技股份有限公司 | 颈动脉扫查方法、装置和计算机可读存储介质 |
CN114224381B (zh) * | 2021-12-16 | 2024-02-27 | 中国人民解放军联勤保障部队北戴河康复疗养中心 | 一种超声检查中辅助支撑装置及其使用方法 |
CN114587347B (zh) * | 2022-03-25 | 2023-04-28 | 深圳市华屹医疗科技有限公司 | 肺功能检测方法、系统、装置、计算机设备和存储介质 |
CN116363330B (zh) * | 2023-02-27 | 2023-09-19 | 优鹰智能科技(深圳)有限公司 | 基于激光扫描技术的空间重构方法、装置、设备及介质 |
CN116587290A (zh) * | 2023-07-17 | 2023-08-15 | 广州多浦乐电子科技股份有限公司 | 工件超声无损检测扫查轨迹生成方法及检测方法 |
CN116781837B (zh) * | 2023-08-25 | 2023-11-14 | 中南大学 | 一种自动化激光三维扫描系统 |
CN117392163B (zh) * | 2023-12-11 | 2024-04-12 | 泉州华中科技大学智能制造研究院 | 基于3d视觉的鞋领口贴纸轨迹提取方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050107703A1 (en) * | 2003-09-22 | 2005-05-19 | Bullis James K. | Ultrasonic imaging with spot focused waves |
US20090024030A1 (en) * | 2007-07-20 | 2009-01-22 | Martin Lachaine | Methods and systems for guiding the acquisition of ultrasound images |
CN102473307A (zh) * | 2010-03-15 | 2012-05-23 | 松下电器产业株式会社 | 用于轨迹估计的方法和装置以及用于分割的方法 |
CN102743188A (zh) * | 2011-04-22 | 2012-10-24 | 李百祺 | 超声波自动扫描系统及其扫描方法 |
WO2013101988A1 (en) * | 2011-12-29 | 2013-07-04 | Maui Imaging, Inc. | M-mode ultrasound imaging of arbitrary paths |
CN103690191A (zh) * | 2013-12-03 | 2014-04-02 | 华南理工大学 | 一种超声探头智能连续扫描装置及其扫描方法 |
CN103999126A (zh) * | 2011-12-16 | 2014-08-20 | 吕贝克大学 | 用于估计姿态的方法和装置 |
CN104856720A (zh) * | 2015-05-07 | 2015-08-26 | 东北电力大学 | 一种基于rgb-d传感器的机器人辅助超声扫描系统 |
CN106361372A (zh) * | 2016-09-22 | 2017-02-01 | 华南理工大学 | 一种超声探头智能扫描路径规划方法 |
CN107330903A (zh) * | 2017-06-29 | 2017-11-07 | 西安理工大学 | 一种人体点云模型的骨架提取方法 |
CN107953336A (zh) * | 2017-12-27 | 2018-04-24 | 北京理工大学 | 机械手超声无损检测中被测件装卡偏差的修正方法及系统 |
CN109671155A (zh) * | 2018-12-21 | 2019-04-23 | 北京林业大学 | 基于点云数据的表面网格重建方法、系统及相关设备 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102007046700A1 (de) * | 2007-09-28 | 2009-04-16 | Siemens Ag | Ultraschallvorrichtung |
CN103750864B (zh) * | 2014-01-13 | 2015-12-02 | 华南理工大学 | 一种超声弹性成像的扫描装置及其扫描方法 |
EP3574841A1 (en) * | 2018-05-28 | 2019-12-04 | Koninklijke Philips N.V. | Ultrasound probe positioning system |
CN110060337B (zh) * | 2019-04-25 | 2021-03-02 | 飞依诺科技(苏州)有限公司 | 颈动脉超声扫查三维重建方法及系统 |
CN110664438B (zh) * | 2019-10-22 | 2021-09-10 | 深圳瀚维智能医疗科技有限公司 | 超声扫查轨迹规划方法、装置、存储介质及计算机设备 |
-
2019
- 2019-10-22 CN CN201911007896.5A patent/CN110664438B/zh active Active
-
2020
- 2020-10-15 WO PCT/CN2020/121233 patent/WO2021078064A1/zh active Application Filing
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050107703A1 (en) * | 2003-09-22 | 2005-05-19 | Bullis James K. | Ultrasonic imaging with spot focused waves |
US20090024030A1 (en) * | 2007-07-20 | 2009-01-22 | Martin Lachaine | Methods and systems for guiding the acquisition of ultrasound images |
CN102473307A (zh) * | 2010-03-15 | 2012-05-23 | 松下电器产业株式会社 | 用于轨迹估计的方法和装置以及用于分割的方法 |
CN102743188A (zh) * | 2011-04-22 | 2012-10-24 | 李百祺 | 超声波自动扫描系统及其扫描方法 |
CN103999126A (zh) * | 2011-12-16 | 2014-08-20 | 吕贝克大学 | 用于估计姿态的方法和装置 |
WO2013101988A1 (en) * | 2011-12-29 | 2013-07-04 | Maui Imaging, Inc. | M-mode ultrasound imaging of arbitrary paths |
CN103690191A (zh) * | 2013-12-03 | 2014-04-02 | 华南理工大学 | 一种超声探头智能连续扫描装置及其扫描方法 |
CN104856720A (zh) * | 2015-05-07 | 2015-08-26 | 东北电力大学 | 一种基于rgb-d传感器的机器人辅助超声扫描系统 |
CN106361372A (zh) * | 2016-09-22 | 2017-02-01 | 华南理工大学 | 一种超声探头智能扫描路径规划方法 |
CN107330903A (zh) * | 2017-06-29 | 2017-11-07 | 西安理工大学 | 一种人体点云模型的骨架提取方法 |
CN107953336A (zh) * | 2017-12-27 | 2018-04-24 | 北京理工大学 | 机械手超声无损检测中被测件装卡偏差的修正方法及系统 |
CN109671155A (zh) * | 2018-12-21 | 2019-04-23 | 北京林业大学 | 基于点云数据的表面网格重建方法、系统及相关设备 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021078066A1 (zh) * | 2019-10-22 | 2021-04-29 | 深圳瀚维智能医疗科技有限公司 | 乳腺超声筛查方法、装置及系统 |
WO2021078064A1 (zh) * | 2019-10-22 | 2021-04-29 | 深圳瀚维智能医疗科技有限公司 | 超声扫查轨迹规划方法、装置、存储介质及计算机设备 |
CN112075957A (zh) * | 2020-07-27 | 2020-12-15 | 深圳瀚维智能医疗科技有限公司 | 乳腺环扫轨迹规划方法、装置及计算机可读存储介质 |
CN112419415A (zh) * | 2020-12-08 | 2021-02-26 | 浙江德尚韵兴医疗科技有限公司 | 一种基于crs曲线拟合实现位姿规划的超声扫查方法 |
CN112419415B (zh) * | 2020-12-08 | 2022-06-17 | 浙江德尚韵兴医疗科技有限公司 | 一种基于crs曲线拟合实现位姿规划的超声扫查方法 |
CN112767309A (zh) * | 2020-12-30 | 2021-05-07 | 无锡祥生医疗科技股份有限公司 | 超声扫查方法、超声设备及系统 |
CN112767237A (zh) * | 2020-12-30 | 2021-05-07 | 无锡祥生医疗科技股份有限公司 | 基于点云数据的环形位姿控制方法、装置及超声设备 |
CN112767237B (zh) * | 2020-12-30 | 2024-06-25 | 无锡祥生医疗科技股份有限公司 | 基于点云数据的环形位姿控制方法、装置及超声设备 |
CN112767415A (zh) * | 2021-01-13 | 2021-05-07 | 深圳瀚维智能医疗科技有限公司 | 胸部扫查区域自动确定方法、装置、设备及存储介质 |
CN116851930A (zh) * | 2023-04-28 | 2023-10-10 | 华南智能机器人创新研究院 | 一种三维五轴曲面激光切割轨迹生成方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
WO2021078064A1 (zh) | 2021-04-29 |
CN110664438B (zh) | 2021-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110664438B (zh) | 超声扫查轨迹规划方法、装置、存储介质及计算机设备 | |
CN110786887B (zh) | 乳腺超声筛查方法、装置及系统 | |
CN110675398B (zh) | 乳腺超声筛查方法、装置及计算机设备 | |
CN110751719B (zh) | 乳房三维点云重建方法、装置、存储介质及计算机设备 | |
CN110766704B (zh) | 乳房点云分割方法、装置、存储介质及计算机设备 | |
CN111292314B (zh) | 冠脉分段方法、装置、图像处理系统和存储介质 | |
CA3149335C (en) | Data processing method, apparatus, device, and storage medium | |
CN111210431B (zh) | 一种血管分割方法、装置、设备及存储介质 | |
CN112244883B (zh) | 基于ct图像提取左心耳数据参数的方法及系统 | |
CN112998749B (zh) | 一种基于视觉伺服的自动超声检查系统 | |
CN115359063B (zh) | 基于目标器官三维图像的刚性配准方法及相关装置 | |
CN113129418B (zh) | 基于三维影像的目标表面重建方法、装置、设备和介质 | |
CN110310364B (zh) | 一种心脏3d打印系统及其打印方法 | |
CN113507890B (zh) | 基于ct图像的肘关节屈伸三维运动分析方法及装置 | |
CN112669450B (zh) | 人体模型构建方法和个性化人体模型构建方法 | |
CN111161419B (zh) | 三维模型的剖切方法、装置和存储介质 | |
CN116687452B (zh) | 一种早孕期胎儿超声自主扫查方法、系统及设备 | |
CN111652954A (zh) | 一种基于左心室分割图片的左心室容积自动计算方法、装置、计算机设备以及存储介质 | |
JP7165541B2 (ja) | ボリュームデータ処理装置、方法及びプログラム | |
KR102426925B1 (ko) | 3d 시뮬레이션을 활용한 수술로봇의 동작 정보를 획득하는 방법 및 프로그램 | |
CN208524899U (zh) | 新型多功能彩色超声系统 | |
WO2017035522A1 (en) | Determining geometric changes, rotations, and/or translations of the heart based on electrocardiogram measurements | |
CN113057678A (zh) | 一种基于双目视觉与机器人的乳腺超声扫描方法及系统 | |
CN117731324A (zh) | 对超声探头在接触表面上的实时力交互控制的方法和装置 | |
CN118266992A (zh) | 一种3d超声造影灌注多参量功能成像方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |