CN104207801A - 一种超声检测图像三维标定方法 - Google Patents
一种超声检测图像三维标定方法 Download PDFInfo
- Publication number
- CN104207801A CN104207801A CN201310221987.5A CN201310221987A CN104207801A CN 104207801 A CN104207801 A CN 104207801A CN 201310221987 A CN201310221987 A CN 201310221987A CN 104207801 A CN104207801 A CN 104207801A
- Authority
- CN
- China
- Prior art keywords
- coordinate
- coordinate system
- template
- intersection point
- cotton thread
- 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
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明涉及一种超声检测图像三维标定方法,该方法通过超声探头上的接收器实时获取超声探头的位姿信息,从而可将模板中任意点(包括棉线交点)坐标转换至接收器坐标系中;通过超声探头对网状模板中棉线交点的扫描,获得棉线交点的超声图像;通过模板中棉线交点构建空间矢量,按照模板中空间矢量间的几何特性计算超声图像中棉线交点成像点的坐标,进而应用同源点最小二乘匹配计算定位接收器坐标系Or与探头扫描图像坐标系Oi间的关系。与现有技术相比,本发明通过空间矢量关系计算图像坐标系中同源点坐标,避免了以往方法中人工拾取图像同源点坐标时人为因素影响,从而具有标定精度与重建精度高、标定速度快等优点。
Description
技术领域
本发明涉及超声标定技术领域,尤其是涉及一种超声检测图像三维标定方法。
背景技术
三维(3-D)超声标定是超声图像3-D可视化不可缺少的关键步骤,因为手持超声探头Freehand采样,得到的是一系列非平行、不规则的二维(2-D)超声图像,只有从超声探头获取的2-D图像信息出发,计算其对应的3-D空间中的几何信息,再由此生成规则的平行图像序列,才能进行后续的三维重建和目标分割。这一定标计算过程称为3-D超声标定,其通常通过在超声探头上固定定位跟踪设备,计算该设备与超声探头扫描平面间的方位关系来实现。该设备可以实时记录超声探头以任意方位扫描目标时的方位信息,因此为扫描得到的一系列非平行超声图像提供了相应的空间方向与位置信息。
目前常用的两种超声探头定位跟踪方法是磁定位与光学定位,本专利采用磁定位方法。而定位跟踪设备一旦固定,其与超声探头扫描平面间的方位关系即已确定,这种转换关系一般不能通过物理方法来测量,但可以表示为定位跟踪设备坐标相对于超声图像平面像素坐标的转换转关系(旋转、平移和比例关系),其数学形式为4×4矩阵。此转换转关系矩阵亦称标定矩阵,以往计算标定矩阵中通过同源点在不同坐标系下的坐标,应用最小二乘匹配算法获得两个坐标系间的转换矩阵。其中同源点在图像坐标系中的坐标依赖人工从图像中拾取,这样不可避免的带来人为因素对标定精度的影响。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷,从同源点在图像中的坐标能够被自动提取的角度进行研究,独特之处在于应用同源点的空间矢量特征关系进行图像坐标的自动计算,从而提供减少算法中人为因素的影响,实现一种标定精度与重建精度高、标定速度快的超声检测图像三维标定方法。
本发明的目的可以通过以下技术方案来实现:
一种超声检测图像三维标定方法,定位接收器与定位发射器连接,该方法通过模板计算固定在超声探头上的定位接收器坐标系Or与探头扫描图像坐标系Oi间的关系,所述的模板上设有两层网状棉线,网状棉线在扫描面上有n个棉线交点,所述的标定方法包括以下步骤:
1)取超声扫描图像中的点I1作为基准点,获取I1在探头扫描图像坐标系上的坐标,I1在定位接收器坐标系中的坐标R1为:
R1=Tr←t·Tt←p·P1
其中,P1为基准点I1在模板坐标系Op中的对应坐标,Tr←t为定位发射器坐标系Ot相对于定位接收器坐标系Or的转换矩阵,Tt←p为模板坐标系Op相对于定位发射器坐标系Ot的转换矩阵;
2)根据模板上棉线交点的位置在扫描图像中拾取点I2,获得向量的方向,I2在模板坐标系Op中的对应坐标为P2,根据P1与P2的距离d12和的方向计算I2在探头扫描图像坐标系上的坐标,I2在定位接收器坐标系中的坐标R2为:
R2=Tr←t·Tt←p·P2;
3)根据模板上棉线交点的位置在扫描图像中拾取点I3,I3在模板坐标系Op中的对应坐标为P3,向量间的夹角为θ,根据和θ确定的方向,根据P1与P3的距离d13和的方向计算I3在探头扫描图像坐标系上的坐标,I3在定位接收器坐标系中的坐标R3为:
R3=Tr←i·Tt←p·P3;
4)重复步骤2)、3),获得扫描面上n个棉线交点在定位接收器坐标系与探头扫描图像坐标系中的坐标,并根据以下公式进行同源点最小二乘原理匹配,获得标定矩阵Tr←t:
其中,Ri为第i个棉线交点在定位接收器坐标系中的坐标,Ii为第i个棉线交点在探头扫描图像坐标系上的坐标,s为比例系数,R为旋转关系矩阵,p为平移向量;
5)采集m幅超声扫描图像,重复步骤1)-4),计算m个标定矩阵Tr←i的平均值作为最终标定结果。
所述的基准点I1为网状棉线中最上层交点的中间一点在超声图像中的成像点。
所述的的方向取决于I2沿基准点I1的方向。
所述的Tt←p通过模板标定方法计算得到,所述的模板标定方法具体为:
模板两侧设有多个通孔,且通孔坐标已知,将定位接收器的探针针尖依次接触通孔,根据针尖相对于其上定位接收器的位置关系及Tr←t获得这些通孔在定位发射器坐标系Ot上的坐标,根据同源点最小二乘原理匹配获得Tt←p。
与现有技术相比,本发明具有以下优点:
1、通过空间矢量的方法计算同源点坐标,是以往研究中没有采用过的方法。以往的方法中通常采用直接拾取图像同源点坐标或应用相似三角形的方法计算同源点坐标。
2、同源点数目少,通过构建空间向量,得到三对同源点即可完成标定算法,标定简单快速。
3、基准点选择在图像中心位置,实验证明此位置受超声束扇形发散效应影响小,可提高标定精度与重建精度。
4、基准点选择在图像偏上位置,实验证明此位置受超声束深度影响小,可提高标定精度与重建精度高。
5、超声图像中另外2个同源点选择在靠近图像中心线的位置,可减少超声束扇形发散效应影响,提高标定精度与重建精度。
6、采用同源点最小二乘原理匹配及多次采集超声图像再求平均值的计算方法,来计算最终的标定矩阵,可以修正单次操作中由于探头对正、图像成像点拾取等因互引入的误差。
附图说明
图1为本发明的原理示意图;
图2为模板扫描面示意图;
图3为超声扫描图像选点示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
一种超声检测图像三维标定方法,采用NDI磁定位方式,定位接收器与定位发射器连接,该方法需辅以模板来完成,通过模板计算固定在超声探头上的定位接收器坐标系Or与探头扫描图像坐标系Oi间的关系,所述的模板上设有两层网状棉线,网状棉线在扫描面上有n(n为10)个棉线交点。
标定过程所需的探头扫描图像坐标系Oi为自定义,定位接收器坐标系Or为NDI公司定义,定位发射器坐标系Ot为NDI公司定义,模板坐标系Op为自定义。
当超声探头如图1所示扫描模板时,取模板上处于图像视野中的某点Pρ=[xkykzk1]T,其成像点在坐标系Oi中的坐标记为Pi=[su·uksv·vk01]T,su、sv为超声图像两个方向的分辨率(比例系数,单位mm/pixel),一般情况下:su=sv,则Op与Oi间的转换关系以齐次坐标变换矩阵表示为:
简写为:Pp=Tp←t·Tt←r·Tr←i·Pi (2)
其中,Tt←r为定位接收器坐标系Or相对于定位发射器坐标系Ot的转换矩阵,由磁定位跟踪装置固有的定位算法计算可得;Tp←t为发射器坐标系Ot相对于模板坐标系Op的转换矩阵,由模板标定方法计算得到。(2)式全部以矩阵形式表示为:
Tp←i=Tp←i·Tt←r·Tr←i (3)
Tr←i=Tt←r -1·Tp←t -1·Tp←i (4)
Tr←i即为标定算法的最终目标。
所述的Tt←p通过模板标定方法计算得到,所述的模板标定方法具体为:
模板两侧设有多个通孔,且通孔坐标已知,将定位接收器的探针针尖依次接触通孔,根据针尖相对于其上定位接收器的位置关系及Tr←t获得这些通孔在定位发射器坐标系Ot上的坐标,根据同源点最小二乘原理匹配获得Tt←p。
如图2所示,由于模板设计的尺寸已知,所以模板上的P1,P2,P2,d12、d13、向量和间的夹角θ是能够计算得到的,均作为已知条件。
如图1-图3所示,该标定方法包括以下步骤:
1)取超声扫描图像中的点I1作为基准点,I1为网状棉线中最上层交点的中间一点,获取I1在探头扫描图像坐标系上的坐标,I1在定位接收器坐标系中的坐标R1为:
R1=Tr←t·Ti←p·P1
其中,P1为基准点I1在模板坐标系Op中的对应坐标,Tr←t为定位发射器坐标系Ot相对于定位接收器坐标系Or的转换矩阵,tt←p为模板坐标系Op相对于定位发射器坐标系Ot的转换矩阵;
2)根据模板上棉线交点的位置在扫描图像中拾取点I2,I2可选临近基准点上层左侧的一点,的方向取决于I2沿基准点I1的方向,可平行于模板坐标系中的x轴方向,I2在模板坐标系Op中的对应坐标为P2,将P1与P2的距离d12转换为像素,计算I1与I2的距离:|I1I2|=s·d12,s为比例系数,单位是:pixel/mm,也就是每毫米代表多少个像素,根据距离和方向计算I2在探头扫描图像坐标系上的坐标,I2在定位接收器坐标系中的坐标R2为:R2=Tr←t·Tt←p·P2;
3)根据模板上棉线交点的位置在扫描图像中拾取点I3,I3可选临近基准点下层左侧的一点,I3在模板坐标系Op中的对应坐标为P3,向量间的夹角为θ,根据和θ确定的方向,根据P1与P3的距离d13和的方向计算I3在探头扫描图像坐标系中的坐标: I3在定位接收器坐标系中的坐标R3为:R3=Tr←t·Tt←p·P3;
4)重复步骤2)、3),获得扫描面上10个棉线交点在定位接收器坐标系与探头扫描图像坐标系中的坐标,并根据以下公式进行同源点最小二乘原理匹配,获得标定矩阵Tr←i:
其中,R1为第i个棉线交点在定位接收器坐标系中的坐标,Ii为第i个棉线交点在探头扫描图像坐标系上的坐标,R为旋转关系矩阵,p为平移向量;
5)采集10幅超声扫描图像,重复步骤1)-4),计算m个标定矩阵Tr←i的平均值作为最终标定结果。
最后的Tr←i是一个4×4矩阵,前三列为旋转关系R,最后一列为平移向量p:
Claims (4)
1.一种超声检测图像三维标定方法,定位接收器与定位发射器连接,其特征在于,该方法通过模板建立同源点空间矢量来计算固定在超声探头上的定位接收器坐标系Or与探头扫描图像坐标系Oi间的关系,所述的模板上设有两层网状棉线,网状棉线在扫描面上有n个棉线交点,所述的标定方法包括以下步骤:
1)取超声扫描图像中的点I1作为基准点,获取I1在探头扫描图像坐标系上的坐标,I1在定位接收器坐标系中的坐标R1为:
R1=Tr←t·Tt←p·P1
其中,P1为基准点I1在模板坐标系Op中的对应坐标,Tr←t为定位发射器坐标系Ot相对于定位接收器坐标系Or的转换矩阵,Tt←p为模板坐标系Op相对于定位发射器坐标系Ot的转换矩阵;
2)根据模板上棉线交点的位置在扫描图像中拾取点I2,获得向量的方向,I2在模板坐标系Op中的对应坐标为P2,根据P1与P2的距离d12和的方向计算I2在探头扫描图像坐标系上的坐标,I2在定位接收器坐标系中的坐标R2为:
R2=Tr←t·Tt←p·P2;
3)根据模板上棉线交点的位置在扫描图像中拾取点I3,I3在模板坐标系Op中的对应坐标为P3,向量间的夹角为θ,根据和θ确定的方向,根据P1与P3的距离d13和的方向计算I3在探头扫描图像坐标系上的坐标,I3在定位接收器坐标系中的坐标R3为:
R3=Tr←t·Tt←p·P3;
4)重复步骤2)、3),获得扫描面上n个棉线交点在定位接收器坐标系与探头扫描图像坐标系中的坐标,并根据以下公式进行同源点最小二乘原理匹配,获得标定矩阵Tr←t:
其中,Ri为第i个棉线交点在定位接收器坐标系中的坐标,Ii为第i个棉线交点在探头扫描图像坐标系上的坐标,s为比例系数,R为旋转关系矩阵,p为平移向量;
5)采集m幅超声扫描图像,重复步骤1)-4),计算m个标定矩阵Tr←t的平均值作为最终标定结果。
2.根据权利要求1所述的一种超声检测图像三维标定方法,其特征在于,所述的基准点I1为网状棉线中最上层交点的中间一点在超声图像中的成像点。
3.根据权利要求1所述的一种超声检测图像三维标定方法,其特征在于,所述的的方向取决于I2沿基准点I1的方向。
4.根据权利要求1所述的一种超声检测图像三维标定方法,其特征在于,所述的Tt←p通过模板标定方法计算得到,所述的模板标定方法具体为:
模板两侧设有多个通孔,且通孔坐标已知,将定位接收器的探针针尖依次接触通孔,根据针尖相对于其上定位接收器的位置关系及Tr←i获得这些通孔在定位发射器坐标系Ot上的坐标,根据同源点最小二乘原理匹配获得Tt←p。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310221987.5A CN104207801B (zh) | 2013-06-05 | 2013-06-05 | 一种超声检测图像三维标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310221987.5A CN104207801B (zh) | 2013-06-05 | 2013-06-05 | 一种超声检测图像三维标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104207801A true CN104207801A (zh) | 2014-12-17 |
CN104207801B CN104207801B (zh) | 2016-08-24 |
Family
ID=52089993
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310221987.5A Active CN104207801B (zh) | 2013-06-05 | 2013-06-05 | 一种超声检测图像三维标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104207801B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105193445A (zh) * | 2015-09-01 | 2015-12-30 | 中国科学院深圳先进技术研究院 | 一种超声探头标定体模、超声探头标定***及其标定方法 |
CN106199507A (zh) * | 2016-06-15 | 2016-12-07 | 北京工业大学 | 基于单程调频连续电磁波的刚体空间位姿直接测量*** |
CN106725595A (zh) * | 2016-12-05 | 2017-05-31 | 华南理工大学 | 一种电磁定位与b超一体化探头标定装置及其标定方法 |
CN110554095A (zh) * | 2019-08-16 | 2019-12-10 | 上海工程技术大学 | 一种三维超声探头标定装置及标定三维超声探头的方法 |
CN110742652A (zh) * | 2019-09-18 | 2020-02-04 | 中国科学院西安光学精密机械研究所 | 一种对消化道肿瘤磁辅助超声图像的三维重建设备及方法 |
CN111449680A (zh) * | 2020-01-14 | 2020-07-28 | 深圳大学 | 一种超声扫描路径的优化方法及超声设备 |
CN111671519A (zh) * | 2020-06-15 | 2020-09-18 | 新博医疗技术有限公司 | 一种超声图像的标定***及标定方法 |
CN113133832A (zh) * | 2021-03-26 | 2021-07-20 | 中国科学院深圳先进技术研究院 | 一种双臂机器人穿刺***标定方法及*** |
CN113827277A (zh) * | 2021-10-21 | 2021-12-24 | 复旦大学 | 一种声致超声成像方法 |
WO2022032964A1 (zh) * | 2020-08-12 | 2022-02-17 | 中国科学院深圳先进技术研究院 | 一种双臂机器人标定方法、***、终端以及存储介质 |
CN114699168A (zh) * | 2022-01-23 | 2022-07-05 | 东南大学 | 一种双平面超声探头的标定***及标定方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103230283A (zh) * | 2013-04-16 | 2013-08-07 | 清华大学 | 一种超声探头成像平面空间位置标定的优化方法 |
CN203324232U (zh) * | 2013-06-05 | 2013-12-04 | 上海工程技术大学 | 一种三维超声定位网状模板 |
CN104215695A (zh) * | 2013-06-05 | 2014-12-17 | 上海工程技术大学 | 一种三维超声定位网状模板 |
-
2013
- 2013-06-05 CN CN201310221987.5A patent/CN104207801B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103230283A (zh) * | 2013-04-16 | 2013-08-07 | 清华大学 | 一种超声探头成像平面空间位置标定的优化方法 |
CN203324232U (zh) * | 2013-06-05 | 2013-12-04 | 上海工程技术大学 | 一种三维超声定位网状模板 |
CN104215695A (zh) * | 2013-06-05 | 2014-12-17 | 上海工程技术大学 | 一种三维超声定位网状模板 |
Non-Patent Citations (1)
Title |
---|
宋章军等: "一种快速、自动的三维超声标定方法", 《北京生物医学工程》, vol. 28, no. 6, 15 December 2009 (2009-12-15) * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017036044A1 (zh) * | 2015-09-01 | 2017-03-09 | 中国科学院深圳先进技术研究院 | 一种超声探头标定体模、超声探头标定***及其标定方法 |
US10555724B2 (en) | 2015-09-01 | 2020-02-11 | Shenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences | Ultrasound probe calibration phantom, ultrasound probe calibration system and calibration method thereof |
CN105193445A (zh) * | 2015-09-01 | 2015-12-30 | 中国科学院深圳先进技术研究院 | 一种超声探头标定体模、超声探头标定***及其标定方法 |
CN106199507A (zh) * | 2016-06-15 | 2016-12-07 | 北京工业大学 | 基于单程调频连续电磁波的刚体空间位姿直接测量*** |
CN106725595A (zh) * | 2016-12-05 | 2017-05-31 | 华南理工大学 | 一种电磁定位与b超一体化探头标定装置及其标定方法 |
CN106725595B (zh) * | 2016-12-05 | 2023-04-21 | 华南理工大学 | 一种电磁定位与b超一体化探头标定装置及其标定方法 |
CN110554095B (zh) * | 2019-08-16 | 2022-07-26 | 上海工程技术大学 | 一种三维超声探头标定装置及标定三维超声探头的方法 |
CN110554095A (zh) * | 2019-08-16 | 2019-12-10 | 上海工程技术大学 | 一种三维超声探头标定装置及标定三维超声探头的方法 |
CN110742652A (zh) * | 2019-09-18 | 2020-02-04 | 中国科学院西安光学精密机械研究所 | 一种对消化道肿瘤磁辅助超声图像的三维重建设备及方法 |
CN111449680A (zh) * | 2020-01-14 | 2020-07-28 | 深圳大学 | 一种超声扫描路径的优化方法及超声设备 |
CN111449680B (zh) * | 2020-01-14 | 2023-07-18 | 深圳大学 | 一种超声扫描路径的优化方法及超声设备 |
CN111671519A (zh) * | 2020-06-15 | 2020-09-18 | 新博医疗技术有限公司 | 一种超声图像的标定***及标定方法 |
WO2022032964A1 (zh) * | 2020-08-12 | 2022-02-17 | 中国科学院深圳先进技术研究院 | 一种双臂机器人标定方法、***、终端以及存储介质 |
CN113133832A (zh) * | 2021-03-26 | 2021-07-20 | 中国科学院深圳先进技术研究院 | 一种双臂机器人穿刺***标定方法及*** |
CN113827277A (zh) * | 2021-10-21 | 2021-12-24 | 复旦大学 | 一种声致超声成像方法 |
CN113827277B (zh) * | 2021-10-21 | 2023-10-03 | 复旦大学 | 一种声致超声成像方法 |
CN114699168A (zh) * | 2022-01-23 | 2022-07-05 | 东南大学 | 一种双平面超声探头的标定***及标定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104207801B (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104207801A (zh) | 一种超声检测图像三维标定方法 | |
CN106595630B (zh) | 一种基于激光导航变电站巡检机器人的建图***与方法 | |
CN100476345C (zh) | 基于双目立体视觉技术的空间圆几何参数的测量方法 | |
CN102589530B (zh) | 基于二维相机和三维相机融合的非合作目标位姿测量方法 | |
CN101569541B (zh) | 三维超声波成像*** | |
CN105157609A (zh) | 基于两组相机的大型零件全局形貌测量方法 | |
CN105067023A (zh) | 一种全景三维激光传感器数据校准方法和装置 | |
CN103793915B (zh) | 神经外科导航中低成本无标记配准***及配准方法 | |
CN105509733A (zh) | 非合作空间圆目标的相对位姿测量方法 | |
CN101504275A (zh) | 一种基于空间无线定位的手持式线激光三维测量*** | |
Lim et al. | Real-time single camera SLAM using fiducial markers | |
CN103759669A (zh) | 一种大型零件的单目视觉测量方法 | |
CN101858730A (zh) | 一种自动测量煤堆体积的方法及其专用装置 | |
CN112825190B (zh) | 一种精度评估方法、***、电子设备和存储介质 | |
CN104990515A (zh) | 大型物体三维形状测量***及其测量方法 | |
CN110065075B (zh) | 一种基于视觉的空间细胞机器人外部状态感知方法 | |
CN106996748A (zh) | 一种基于双目视觉的轮径测量方法 | |
CN104034269A (zh) | 一种单目视觉测量方法与装置 | |
CN104634248A (zh) | 一种双目视觉下的转轴标定方法 | |
CN101194845A (zh) | 基于双x射线机的同步实时动态跟踪立体定位方法及装置 | |
CN103900504A (zh) | 纳米尺度下的实时三维视觉信息反馈方法 | |
CN111126116A (zh) | 无人船河道垃圾识别方法及*** | |
CN101865656B (zh) | 一种使用少数共面点精确定位多摄像头***位姿的方法 | |
CN112168357A (zh) | C臂机空间定位模型构建***及方法 | |
CN102999895B (zh) | 利用两个同心圆线性求解摄像机内参数 |
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 |