CN110084855B - 一种改进cbct几何参数标定方法 - Google Patents
一种改进cbct几何参数标定方法 Download PDFInfo
- Publication number
- CN110084855B CN110084855B CN201910318861.7A CN201910318861A CN110084855B CN 110084855 B CN110084855 B CN 110084855B CN 201910318861 A CN201910318861 A CN 201910318861A CN 110084855 B CN110084855 B CN 110084855B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- ray source
- phantom
- ellipse
- point
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; Lesion
Abstract
本发明公开了一种改进的CBCT几何参数标定方法,通过空间几何关系找到不同坐标系下的变换关系,以参数估计的方法提高标定参数的精确值。解析几何和特征点迭代优化相结合的方式精确计算标定参数,迭代估计方法的引入在提高精度的前提下也克服现有算法在标定过程对标定体模位置的要求,达到忽略体模位置而精确标定参数的目的。
Description
技术领域
本发明涉及医学物理技术领域,具体涉及一种改进CBCT几何参数标定方法。
背景技术
CBCT图像引导系统主要用于辅助引导质子肿瘤治疗,通过与计划CT图像进行配准获取患者摆位误差。因此CT图像和患者的位置坐标需要精确定位。
现有的标定算法大多假设标定体模旋转,而投影装置固定不变。但是对于CBCT装置此类方法在实际应用中存在先天上的缺陷。符合此类要求的CBCT标定算法,通过设计特定的标定体模,计算各种情况下参数变化,但是对于实际场景考虑不够全面,部分参数计算存在误差。
发明内容
为了解决上述的技术问题,本发明的目的在于提供一种改进CBCT几何参数标定方法。
本发明所要解决的技术问题为:
(1)如何进行部分参数的优化。
(2)如何适配实际应用场景。
本发明的目的可以通过以下技术方案实现:一种改进CBCT几何参数标定方法,该方法的步骤如下:
1)按要求设计体模,体模中钢球位置误差在±0.001毫米以内;
2)360度全扇扫描,按顺时针或逆时针方向依次采集体模的投影图像;
3)对投影图像进行特征点检测与追踪,每一帧图像中检测的特征点形成上下两个椭圆,按顺时针方向对特征点坐标进行重新排序,再以差帧的方法追踪特征点,保证特征点一一对应,并在每副图像中利用检测到的特征点坐标拟合成两个椭圆,椭圆方程为:
a(x-u0)2+b(y-v0)2+2c(x-u0)(y-v0)=1;
其中,a、b、c为椭圆参数,u0、v0为椭圆中心坐标;
4)根据投影图像建立空间坐标系;
6)利用已知参数和初值进一步优化估计得到精确解;
7)利用在步骤5)、6)中计算的参数关系,计算世界坐标系和投影图像之间的变换关系。
进一步的,所述步骤5)的具体计算步骤如下:
5-1)面内转角η利用如下公式计算:
5-2)使用以下公式计算θ和φ角度:
其中ak、bk、ck、εk分别为特征点拟合的椭圆参数;其中Pφ=(uφ,vφ)为两个椭圆长轴所在直线交点,与现实平板绕Y轴旋转角度有关;Pθ=(uθ,vθ)为椭圆最两侧点连线构成的交点,与现实平板绕X轴旋转角度有关,其变换矩阵为:
5-3)计算射线源在各个坐标系中的坐标位置,通过以下计算公式可以粗略计算出射线源在各个坐标系中的位置:
进一步的,步骤6)中得到精确解的具体步骤如下:
6-1)特征点位于虚拟坐标系下的坐标可表示为:
6-2)设以体模为基准,图像旋转角度为t,由投影关系可得射线源到体模坐标点的向量和射线源到投影点的坐标的向量对应成比例,得到两个比例关系;每一帧图像由特征点为已知量和5)中得到的θ,φ,η为参数,使用Levenberg-Marquardt非线性优化算法进行参数估计,以步骤5)中计算的结果为初值,t以上一帧计算的结果为初值,估计精确值。
进一步的,步骤7)中计算世界坐标系和投影图像之间的变换关系的具体步骤如下:
7-1)每一帧计算一个射线源位置坐标,在步骤6)中计算的相对体模旋转的前提下,得到射线源在360度旋转扫描下的运动轨迹;
7-2)以第一帧射线源方向为参考方向,射线源轨迹将形成一个偏心圆,圆心即为等中心点;
7-3)利用已知参数,以等中心点为中心构成新的世界坐标系,进行坐标系参数变换。
本发明的有益效果:
本发明针对现有技术进行改进,提高参数计算的准确性。通过特征点的位置关系,采用优化算法的方式,迭代优化部分参数,克服在标定时由于体模摆放位置的不确定而带来的计算误差。
附图说明
下面结合附图对本发明作进一步的说明。
图1是设计标定体模的结构图;
图2是步骤4)中建立的三个几何坐标系;
图3是步骤4)中描述的实际平板相对虚拟平板的旋转关系;
图4是标定体模所产生的投影图像其中的一帧;
图5坐标系变化后射线源和投影点的坐标轨迹。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
请参阅图1-5所示,本实施例提供了本发明的目的可以通过以下技术方案实现:一种改进CBCT几何参数标定方法,该方法的步骤如下:
1)按要求设计体模,体模中钢球位置误差在±0.001毫米以内;首先设计如图1所示的标定体模,由钢球组成的两个圆形环组成,圆环半径为115毫米,封装在一个中空的丙烯酸圆柱体中。每个环包含12个钢球,每个钢球半径为2.5毫米,两环的间距为210毫米。利用激光校准,使体模尽可能置于等中心点,圆环中心连线与旋转轴尽可能平行。
2)360度全扇扫描,按顺时针或逆时针方向依次采集体模的投影图像;对标定体模进行全扇扫描,按顺序采集720帧标定物体的投影图像,在每一帧图像的计算中,建立如图2所示坐标系,平板的偏转角度如图3所示。每一帧图像进行去噪、分割、阈值化等处理后得到类似于图4所示的图像。以圆检测的方法,检测到图像中的圆心作为特征点。
3)对投影图像进行特征点检测与追踪,每一帧图像中检测的特征点形成上下两个椭圆,按顺时针方向对特征点坐标进行重新排序,再以差帧的方法追踪特征点,保证特征点一一对应,并在每副图像中利用检测到的特征点坐标拟合成两个椭圆,椭圆方程为:
a(x-u0)2+b(y-v0)2+2c(x-u0)(y-v0)=1;
其中,a、b、c为椭圆参数,u0、v0为椭圆中心坐标;
4)根据投影图像建立空间坐标系;
5-1)面内转角η利用如下公式计算:
5-2)使用以下公式计算θ和φ角度:
其中ak、bk、ck、εk分别为特征点拟合的椭圆参数;其中Pφ=(uφ,vφ)为两个椭圆长轴所在直线交点,与现实平板绕Y轴旋转角度有关;Pθ=(uθ,vθ)为椭圆最两侧点连线构成的交点,与现实平板绕X轴旋转角度有关,其变换矩阵为:
5-3)计算射线源在各个坐标系中的坐标位置,通过以下计算公式可以粗略计算出射线源在各个坐标系中的位置:
6)利用已知参数和初值进一步优化估计得到精确解;步骤6)中得到精确解的具体步骤如下:
6-1)特征点位于虚拟坐标系下的坐标可表示为:
6-2)设以体模为基准,图像旋转角度为t,由投影关系可得射线源到体模坐标点的向量和射线源到投影点的坐标的向量对应成比例,得到两个比例关系;每一帧图像由特征点为已知量和5)中得到的θ,φ,η为参数,使用Levenberg-Marquardt非线性优化算法进行参数估计,以步骤5)中计算的结果为初值,t以上一帧计算的结果为初值,估计精确值。
7)利用在步骤5)、6)中计算的参数关系,计算世界坐标系和投影图像之间的变换关系。上文计算出所有坐标系之间的变换关系,然后将体模坐标系下的参数转化到以等中心点为中心的世界坐标系下的参数,得到转换坐标系后射线源和等中心点投影点在旋转360度下的运动轨迹如图5所示。
步骤7)中计算世界坐标系和投影图像之间的变换关系的具体步骤如下:
7-1)每一帧计算一个射线源位置坐标,在步骤6)中计算的相对体模旋转的前提下,得到射线源在360度旋转扫描下的运动轨迹;
7-2)以第一帧射线源方向为参考方向,射线源轨迹将形成一个偏心圆,圆心即为等中心点;
7-3)利用已知参数,以等中心点为中心构成新的世界坐标系,进行坐标系参数变换。
以上内容仅仅是对本发明结构所作的举例和说明,所属本技术领域的技术人员对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,只要不偏离发明的结构或者超越本权利要求书所定义的范围,均应属于本发明的保护范围。
Claims (1)
1.一种改进CBCT几何参数标定方法,其特征在于,该方法的步骤如下:
1)按要求设计体模,体模中钢球位置误差在±0.001毫米以内;
2)360度全扇扫描,按顺时针或逆时针方向依次采集体模的投影图像;
3)对投影图像进行特征点检测与追踪,每一帧图像中检测的特征点形成上下两个椭圆,按顺时针方向对特征点坐标进行重新排序,再以差帧的方法追踪特征点,保证特征点一一对应,并在每副图像中利用检测到的特征点坐标拟合成两个椭圆,椭圆方程为:
a(x-u0)2+b(y-v0)2+2c(x-u0)(y-v0)=1;
其中,a、b、c为椭圆参数,u0、v0为椭圆中心坐标;
4)根据投影图像建立空间坐标系;
6)利用已知参数和初值进一步优化估计得到精确解;
7)利用在步骤5)、6)中计算的参数关系,计算世界坐标系和投影图像之间的变换关系;
所述步骤5)的具体计算步骤如下:
5-1)面内转角η利用如下公式计算:
5-2)使用以下公式计算θ和φ角度:
其中ak、bk、ck、εk分别为特征点拟合的椭圆参数;其中Pφ=(uφ,vφ)为两个椭圆长轴所在直线交点,与现实平板绕Y轴旋转角度有关;Pθ=(uθ,vθ)为椭圆最两侧点连线构成的交点,与现实平板绕X轴旋转角度有关,其变换矩阵为:
5-3)计算射线源在各个坐标系中的坐标位置,通过以下计算公式可以粗略计算出射线源在各个坐标系中的位置:
步骤6)中得到精确解的具体步骤如下:
6-1)特征点位于虚拟坐标系下的坐标可表示为:
6-2)设以体模为基准,图像旋转角度为t,由投影关系可得射线源到体模坐标点的向量和射线源到投影点的坐标的向量对应成比例,得到两个比例关系;每一帧图像由特征点为已知量和步骤5)中得到的θ,φ,η为参数,使用Levenberg-Marquardt非线性优化算法进行参数估计,以步骤5)中计算的结果为初值,t以上一帧计算的结果为初值,估计精确值;
步骤7)中计算世界坐标系和投影图像之间的变换关系的具体步骤如下:
7-1)每一帧计算一个射线源位置坐标,在步骤6)中计算的相对体模旋转的前提下,得到射线源在360度旋转扫描下的运动轨迹;
7-2)以第一帧射线源方向为参考方向,射线源轨迹将形成一个偏心圆,圆心即为等中心点;
7-3)利用已知参数,以等中心点为中心构成新的世界坐标系,进行坐标系参数变换。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910318861.7A CN110084855B (zh) | 2019-04-19 | 2019-04-19 | 一种改进cbct几何参数标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910318861.7A CN110084855B (zh) | 2019-04-19 | 2019-04-19 | 一种改进cbct几何参数标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110084855A CN110084855A (zh) | 2019-08-02 |
CN110084855B true CN110084855B (zh) | 2020-12-15 |
Family
ID=67415760
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910318861.7A Active CN110084855B (zh) | 2019-04-19 | 2019-04-19 | 一种改进cbct几何参数标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110084855B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112634353B (zh) * | 2020-12-17 | 2024-03-26 | 广州华端科技有限公司 | Cbct系统几何标定模体自标定方法、装置及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101936720A (zh) * | 2010-07-30 | 2011-01-05 | 北京航空航天大学 | 一种适用于锥束xct系统的探测器扭转角的标定方法 |
CN107592935A (zh) * | 2015-03-20 | 2018-01-16 | 伦斯勒理工学院 | X射线ct的自动系统校准方法 |
CN108982556A (zh) * | 2018-08-22 | 2018-12-11 | 武汉科技大学 | 一种ct参数标定体膜、ct参数标定系统及ct参数标定方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101750021B (zh) * | 2009-12-04 | 2011-05-11 | 深圳先进技术研究院 | Ct系统中几何参数的标定方法、装置 |
CN102743184B (zh) * | 2012-05-14 | 2013-10-16 | 清华大学 | 一种x射线锥束计算机层析成像系统的几何参数标定方法 |
EP2884897A1 (en) * | 2012-08-20 | 2015-06-24 | orangedental GmbH & Co. KG | Geometric characterization and calibration of a cone-beam computer tomography apparatus |
CN103006251B (zh) * | 2012-12-06 | 2015-02-18 | 深圳先进技术研究院 | 用于ct系统中标定几何参数的标定体模、标定装置及标定方法 |
CN104764756B (zh) * | 2015-03-30 | 2017-05-31 | 清华大学 | 锥束ct成像系统的标定方法 |
CN108122203B (zh) * | 2016-11-29 | 2020-04-07 | 上海东软医疗科技有限公司 | 一种几何参数的校正方法、装置、设备及系统 |
CN106994023B (zh) * | 2017-05-27 | 2018-02-23 | 广州华端科技有限公司 | 锥形束计算机断层成像系统的几何参数确定方法 |
CN107684435A (zh) * | 2017-08-16 | 2018-02-13 | 深圳先进技术研究院 | 锥束ct系统几何校准方法及其校准装置 |
-
2019
- 2019-04-19 CN CN201910318861.7A patent/CN110084855B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101936720A (zh) * | 2010-07-30 | 2011-01-05 | 北京航空航天大学 | 一种适用于锥束xct系统的探测器扭转角的标定方法 |
CN107592935A (zh) * | 2015-03-20 | 2018-01-16 | 伦斯勒理工学院 | X射线ct的自动系统校准方法 |
CN108982556A (zh) * | 2018-08-22 | 2018-12-11 | 武汉科技大学 | 一种ct参数标定体膜、ct参数标定系统及ct参数标定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110084855A (zh) | 2019-08-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109118545B (zh) | 基于旋转轴和双目摄像头的三维成像系统标定方法及系统 | |
CN108053450B (zh) | 一种基于多约束的高精度双目相机标定方法 | |
US7844094B2 (en) | Systems and methods for determining geometric parameters of imaging devices | |
CN113409285B (zh) | 一种沉管隧道接头三维变形的监测方法、系统 | |
CN108245788B (zh) | 一种双目测距装置及方法、包括该装置的加速器放疗系统 | |
CN111973212B (zh) | 参数标定方法和参数标定装置 | |
WO2018126335A1 (zh) | 基于小球模体的锥束ct系统几何参数评价及校正方法 | |
KR101320712B1 (ko) | 정렬 마크를 이용한 2축 회전 스테이지의 회전축 정렬 방법 및 그 장치 | |
CN109064516A (zh) | 一种基于绝对二次曲线像的相机自标定方法 | |
CN110503713B (zh) | 一种基于轨迹平面法向量和圆心结合的旋转轴估计方法 | |
JP6746050B2 (ja) | キャリブレーション装置、キャリブレーション方法およびキャリブレーションプログラム | |
CN110084855B (zh) | 一种改进cbct几何参数标定方法 | |
Mathiassen et al. | Real-time biopsy needle tip estimation in 2D ultrasound images | |
CN109445229B (zh) | 一种获取含一阶径向畸变的变焦相机焦距的方法 | |
Cui et al. | The circular mark projection error compensation in camera calibration | |
CN113876346A (zh) | 一种倾斜图像的迭代校正方法 | |
Duan et al. | Knowledge-based self-calibration method of calibration phantom by and for accurate robot-based CT imaging systems | |
TWI639414B (zh) | 電腦斷層掃描影像的對位方法 | |
Lu et al. | On the sensitivity analysis of camera calibration from images of spheres | |
CN115999079A (zh) | 一种调强放疗设备等中心位置校准方法 | |
CN113643428A (zh) | 一种适用于多自由度锥形束ct的全参数几何校准方法 | |
CN109308706B (zh) | 一种通过图像处理得到三维曲面面积的方法 | |
Daponte et al. | A stereo vision method for IMU-based motion tracking systems calibration | |
CN115227397B (zh) | 一种配准板自动对齐方法和装置 | |
Cheng et al. | Camera Calibration Based on Phase Estimation |
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 |