CN113804105B - 一种加速器射野几何与透视成像几何耦合校准装置及方法 - Google Patents
一种加速器射野几何与透视成像几何耦合校准装置及方法 Download PDFInfo
- Publication number
- CN113804105B CN113804105B CN202110942168.4A CN202110942168A CN113804105B CN 113804105 B CN113804105 B CN 113804105B CN 202110942168 A CN202110942168 A CN 202110942168A CN 113804105 B CN113804105 B CN 113804105B
- Authority
- CN
- China
- Prior art keywords
- geometric
- accelerator
- die body
- perspective
- perspective imaging
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000008878 coupling Effects 0.000 title claims abstract description 19
- 238000010168 coupling process Methods 0.000 title claims abstract description 19
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 19
- 239000002184 metal Substances 0.000 claims abstract description 33
- 239000011159 matrix material Substances 0.000 claims abstract description 30
- 239000012780 transparent material Substances 0.000 claims abstract description 3
- 230000005855 radiation Effects 0.000 claims description 10
- 239000000758 substrate Substances 0.000 claims 4
- VVQNEPGJFQJSBK-UHFFFAOYSA-N Methyl methacrylate Chemical compound COC(=O)C(C)=C VVQNEPGJFQJSBK-UHFFFAOYSA-N 0.000 claims 1
- 229920005372 Plexiglas® Polymers 0.000 claims 1
- 239000011521 glass Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010330 laser marking Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- 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]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
本发明公开了一种加速器射野几何与透视成像几何耦合校准装置,包括透明材料制成的基体和嵌在基体内部的金属球;任意四个球体不共面,且一个金属球球心位于基体的几何中心;基体表面设有两组标记线。本发明还公开一种基于上述装置的加速器射野几何与透视成像几何耦合校准方法,将校准装置置于加速器射野范围和透视成像装置成像范围内,根据不同角度的影像中金属球在画面中的位置,加速器射野中心点和透视成像装置的几何信息获得透视成像装置相对于加速器射野中心点的耦合几何信息,本发明可实时获得相对于加速器射野的透视装置几何参数,避免了基于激光灯模体摆位的误差传递,提升了用于加速器图像引导的透视装置的几何精度,整体校正效率高。
Description
技术领域
本发明涉及加速器的校准装置及方法,尤其涉及一种加速器射野几何与透视成像几何耦合校准装置及方法。
背景技术
在CBCT(锥形束CT)、透视图像引导摆位等透视图像的应用中,对于透视成像装置的几何参数有精度要求,在获取透视图像的同时需要得到形成透视图的投影透视成像装置的几何参数,这些几何参数指透视投影矩阵或与之等价的透视源坐标、成像平板坐标与方向等参数,这些几何参数需要以加速器的射野等中心为基准获得。
目前该领域的应用中,在获得透视装置的几何参数时,选取基准一般使用预先校准的纵横两组激光灯来标记射野坐标系下各坐标轴的方向,通过激光灯交点标记射野等中心位置。但这样的方法受激光灯校准的误差影响,并随着加速器的使用误差会发生递增。
发明内容
发明目的:本发明的目的是提供一种减小误差影响的加速器校准装置及方法。
技术方案:本发明的一种加速器射野几何与透视成像几何耦合校准装置,其特征在于,所述装置包括透明材料制成的基体和嵌在基体内部的若干金属球。金属球数量至少为透视成像装置的数量的八倍,金属球直径为3-10mm;任意四个球体不共面,其中一个金属球的球心位于基体的几何中心。
基体由有机玻璃制成,下端面为平面,其表面设有两组标记线,每组标记线包括一条竖直标记线和一条水平标记线,通过两组标记线中心点的两组水平面与竖直面交于基体的几何中心,通过两组标记线的水平面与竖直面的两条交线相互垂直。
本发明所述的使用上述校准装置的加速器射野几何与透视成像几何耦合校准方法,包括如下步骤:
(1)将校准装置放于平台上,调整其位置确保校准装置同时处于加速器射野范围和透视成像装置成像范围内;
(2)旋转加速器机头,获得多于2个角度的加速器射野影像,根据这些影像中金属球在画面中的位置于金属球在模体坐标系下的坐标,获得加速器成像的几何信息,包括加速器机头、射野成像平板的坐标和角度;
(3)填充射野影像中金属投影区域并卷积获得射野中心投影坐标,并根据该角度的射野成像装置几何信息计算得射野该角度的三维射野中心射线;由多个角度的射野中心射线拟合得到其三维交点,即模体坐标系下的射野等中心点(x0,y0,z0);
(4)不移动模体,启动透视成像装置,获得透视影像,从透视影像中获得金属球二维投影坐标,与金属球三维坐标联立方程计算获得透视成像装置在模体坐标系下的几何投影矩阵表示的几何信息,矩阵为3x4的形式,联立方程如下:
其中pij指投影矩阵第i行第j列的元素的值,(xi,yi,zi)指第i个金属球的三维坐标,(ui,vi,wi)为三维坐标的二维投影加权坐标,对应的二维投影坐标计算方法为:
(5)根据步骤(3)获得的模体坐标系下的加速器射野中心点坐标形成模体坐标系原点到加速器中心点的平移向量(-x0,-y0,-z0),代入步骤(4)获得的模体坐标系下的几何投影矩阵,最终通过下式形成最终的耦合加速器坐标系下的几何投影矩阵P:
模体坐标系按如下原则形成:以通过第一组标记线的水平面与竖直面交线为x轴,通过第二组标记线的水平面与竖直面的交线为y轴,根据右手定则由x轴和y轴形成z轴,以模体几何中心为坐标原点;记录模体坐标系下各个金属球的位置。
有益效果:与现有技术相比,本发明具有如下显著优点:实现了加速器射野与透视成像耦合的几何校正方法,可实时获得相对于加速器射野的透视装置几何参数,避免了基于激光灯模体摆位的误差传递,提升了用于加速器图像引导的透视装置的几何精度,并且可将激光校准与透视成像几何校准流程合并,提升整体校正效率。
附图说明
图1为本发明的流程图;
图2为本发明加速器射野与透视成像装置示意图;
图3为本发明加速器多角度成像示意图。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
实施例1的加速器射野几何与透视成像几何耦合校准装置由30个直径3mm的金属球嵌入有机玻璃制成的基体中制成,30个金属球各向非共面,基体下端面为平面,以平稳放置在台面;其中任意4个球体不共面,且有一个金属球球心位于基体的几何中心点;模体装置有两组标记线,每组标记线由一条水平标记线和一条竖直标记线组成,通过两组标记线中心点的两组水平面与竖直面交于基体的几何中心,通过两组标记线的水平面与竖直面的两条交线相互垂直。
使用上述装置的加速器射野几何与透视成像几何耦合校准方法包括如下步骤:
(1)将模体放置于加速器下,并调整其位置,确保模体同时处于加速器射野范围和透视成像装置成像范围内,对于已有激光灯的加速器治疗室,可直接利用激光灯放置于治疗床上,将两组激光灯的扫描线分别配准模体侧面的两组横竖标记线,使激光扫描线交点与模体标记线交点重合;这样的情况下,可认为激光坐标系与模体坐标系一致,可利用后续射野中心几何直接校准激光标记***;
(2)旋转加速器机头,获得多于2个角度的加速器射野影像,根据这些影像中金属球在画面中的位置于金属球在模体坐标系下的坐标,获得加速器成像的几何信息;如图3所示获得2个角度的加速器射野影像后计算模体坐标系下射野等中心点的坐标,具体步骤如下:
(21)这些成射野影像的角度需不平行且不近似平行,以便提升后续射野中心射线交点计算精度;
(22)根据这些影像中金属球在画面中的位置于金属球在模体坐标系下的坐标,获得加速器机头、射野成像平板的坐标、角度等完备射野成像装置几何信息,每个角度生成一组这样的信息。这些信息包括:射野源点坐标、射野影像行列方向、射野影像基准点坐标,该几何参数即等价于三维坐标经过射野成像***投影到射野影像上的二维坐标的3x4投影矩阵,该过程通过计算机自动优化,以获得能够精确拟合模体内金属点的投影点在图像坐标系下的坐标和其在模体坐标系下的三维坐标的投影矩阵;
(23)通过填充射野影像中金属投影区域并卷积获得射野中心投影坐标,填充射野影像,通过射野影像内高于一定阈值的像素的值与其图像坐标的加权计算获得射野中心的图像坐标,结合投影矩阵可计算得到投影点坐标,等价于投影中心射线的几何信息;
(24)如图3所示,由多个角度的射野中心射线拟合得到其三维交点,即模体坐标系下的射野等中心点,拟合过程为找到与各投影中心射线距离最小的点为射野等中心点;
(3)填充射野影像中金属投影区域并卷积获得射野中心投影坐标,并根据该角度的射野成像装置几何信息计算得射野该角度的三维射野中心射线;由多个角度的射野中心射线拟合得到其三维交点,即模体坐标系下的射野等中心点;
(4)不移动模体,启动透视成像装置,获得透视影像,并计算获得透视成像装置在模体坐标系下的几何信息,该步骤与(22)类似;
(5)根据步骤(3)、(4)获得的加速器射野等中心点和透视成像装置的几何信息获得透视成像装置相对于加速器射野等中心点的耦合几何信息。
Claims (7)
1.一种加速器射野几何与透视成像几何耦合校准装置的校准方法,其特征在于,所述装置包括透明材料制成的基体和嵌在基体内部的若干金属球,所述金属球直径相同,任意四个球体不共面,其中一个金属球的球心位于基体的几何中心,所述方法包括如下步骤:
(1)将校准装置放于平台上,调整其位置确保校准装置同时处于加速器射野范围和透视成像装置成像范围内;
(2)旋转加速器机头,获得多于2个角度的加速器射野影像,根据这些影像中金属球在画面中的位置于金属球在模体坐标系下的坐标,获得加速器成像的几何信息;
(3)填充射野影像中金属投影区域并卷积获得射野中心投影坐标,并根据该角度的射野成像装置几何信息计算得射野该角度的三维射野中心射线;由多个角度的射野中心射线拟合得到其三维交点,即模体坐标系下的射野等中心点;
(4)不移动模体,启动透视成像装置,获得透视影像,从透视影像中获得金属球二维投影坐标,与金属球三维坐标联立方程计算获得透视成像装置在模体坐标系下的几何投影矩阵表示的几何信息,矩阵为3x4的形式,联立方程如下:
其中pij指投影矩阵第i行第j列的元素的值,(xi,yi,zi)指第i个金属球的三维坐标,(ui,vi,wi)为三维坐标的二维投影加权坐标,对应的二维投影坐标计算方法为:
(5)根据步骤(3)获得的模体坐标系下的加速器射野中心点坐标形成模体坐标系原点到加速器中心点的平移向量(-x0,-y0,-z0),代入步骤(4)获得的模体坐标系下的几何投影矩阵,最终通过下式形成最终的耦合加速器坐标系下的几何投影矩阵P:
2.根据权利要求1所述的加速器射野几何与透视成像几何耦合校准装置的校准方法,其特征在于,所述模体坐标系按如下原则形成:以通过第一组标记线的水平面与竖直面交线为x轴,通过第二组标记线的水平面与竖直面的交线为y轴,根据右手定则由x轴和y轴形成z轴,以模体几何中心为坐标原点;记录模体坐标系下各个金属球的位置。
3.根据权利要求2所述的加速器射野几何与透视成像几何耦合校准装置的校准方法,其特征在于,所述步骤(2)中的几何信息包括加速器机头、射野成像平板的坐标和角度。
4.根据权利要求1所述的加速器射野几何与透视成像几何耦合校准装置的校准方法,其特征在于,所述基体表面设有两组标记线,每组标记线包括一条竖直标记线和一条水平标记线,通过两组标记线中心点的两组水平面与竖直面交于基体的几何中心,通过两组标记线的水平面与竖直面的两条交线相互垂直。
5.根据权利要求1所述的加速器射野几何与透视成像几何耦合校准装置的校准方法,其特征在于,所述基体由有机玻璃制成。
6.根据权利要求1所述的加速器射野几何与透视成像几何耦合校准装置的校准方法,其特征在于,所述基体下端面为平面。
7.根据权利要求4所述的加速器射野几何与透视成像几何耦合校准装置的校准方法,其特征在于,所述金属球直径为3-10mm,所述金属球数量至少为用于进行校准的透视成像装置的数量的八倍。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110942168.4A CN113804105B (zh) | 2021-08-17 | 2021-08-17 | 一种加速器射野几何与透视成像几何耦合校准装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110942168.4A CN113804105B (zh) | 2021-08-17 | 2021-08-17 | 一种加速器射野几何与透视成像几何耦合校准装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113804105A CN113804105A (zh) | 2021-12-17 |
CN113804105B true CN113804105B (zh) | 2024-04-05 |
Family
ID=78893661
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110942168.4A Active CN113804105B (zh) | 2021-08-17 | 2021-08-17 | 一种加速器射野几何与透视成像几何耦合校准装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113804105B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116850484B (zh) * | 2023-08-17 | 2024-03-26 | 迈胜医疗设备有限公司 | 图像引导***及校准装置、位置校准方法、放射治疗设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN200951265Y (zh) * | 2006-08-11 | 2007-09-26 | 天津医科大学附属肿瘤医院 | 放疗设备的射野测量及等中心检测装置 |
CN104783824A (zh) * | 2014-01-20 | 2015-07-22 | 上海联影医疗科技有限公司 | X射线成像***的校正方法 |
CN211669376U (zh) * | 2020-02-25 | 2020-10-13 | 闫冰 | 一种用于加速器日检的模体 |
CN112837246A (zh) * | 2021-03-25 | 2021-05-25 | 中科超精(南京)科技有限公司 | 一种图像引导***几何校正方法 |
CN112837358A (zh) * | 2021-03-25 | 2021-05-25 | 中科超精(南京)科技有限公司 | 一种多透视图耦合的台面物体定位方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005015125A1 (en) * | 2003-08-08 | 2005-02-17 | University Health Network | Method and system for calibrating a source and detector instrument |
US8007173B2 (en) * | 2009-10-14 | 2011-08-30 | Siemens Medical Solutions Usa, Inc. | Calibration of imaging geometry parameters |
-
2021
- 2021-08-17 CN CN202110942168.4A patent/CN113804105B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN200951265Y (zh) * | 2006-08-11 | 2007-09-26 | 天津医科大学附属肿瘤医院 | 放疗设备的射野测量及等中心检测装置 |
CN104783824A (zh) * | 2014-01-20 | 2015-07-22 | 上海联影医疗科技有限公司 | X射线成像***的校正方法 |
CN211669376U (zh) * | 2020-02-25 | 2020-10-13 | 闫冰 | 一种用于加速器日检的模体 |
CN112837246A (zh) * | 2021-03-25 | 2021-05-25 | 中科超精(南京)科技有限公司 | 一种图像引导***几何校正方法 |
CN112837358A (zh) * | 2021-03-25 | 2021-05-25 | 中科超精(南京)科技有限公司 | 一种多透视图耦合的台面物体定位方法 |
Non-Patent Citations (1)
Title |
---|
利用电子射野影像***快速调整立体定向放疗***等中心精度;康德华;邓小武;黄劭敏;周颖娟;;中国肿瘤(第08期);681-684 * |
Also Published As
Publication number | Publication date |
---|---|
CN113804105A (zh) | 2021-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11633629B2 (en) | Method of calibration of a stereoscopic camera system for use with a radio therapy treatment apparatus | |
US6594623B1 (en) | Determining three-dimensional orientation of objects | |
CN111514476B (zh) | 一种用于x射线图像引导***中的校准方法 | |
CN103735282B (zh) | 一种锥束ct***探测器几何校正装置及其校正方法 | |
WO2018218611A1 (zh) | 锥形束计算机断层成像***的几何参数确定方法 | |
US20190001156A1 (en) | Determination of geometrical information about a medical treatment arrangement comprising a rotatable treatment radiation source unit | |
CN107595388B (zh) | 一种基于反光球标记点的近红外双目视觉立体匹配方法 | |
CN109925052B (zh) | 靶点路径的确定方法、装置和***、可读存储介质 | |
CN1612713A (zh) | 用于外部波束放射治疗的记录,引导和目标命中的装置和方法 | |
CN112932665B (zh) | 一种用于x光成像的双位置手术定位标尺 | |
CN113804105B (zh) | 一种加速器射野几何与透视成像几何耦合校准装置及方法 | |
CN108022267A (zh) | 一种摄像机标定装置及标定方法 | |
EP3294137A1 (en) | A monitoring system | |
US20220203131A1 (en) | Radiotherapy system, and device and method for verifying same | |
JP2016221156A (ja) | 放射線治療装置較正用ファントムおよびこれを用いた放射線治療装置の位置合わせ方法 | |
CN114404041B (zh) | C型臂成像参数校准***及方法 | |
CN1853737A (zh) | 一种放射治疗的定位方法 | |
CN112603346B (zh) | 一种基于标记物成像的探测器偏转校正方法 | |
CN114596222A (zh) | 适用于一般轨迹锥束ct***几何校正的模体与标定方法 | |
CN117351091A (zh) | 一种相机阵列标定装置及其使用方法 | |
CN210228333U (zh) | 定位*** | |
CN115836875A (zh) | 校正方法及*** | |
CN213491591U (zh) | 一种用于x射线图像引导***中的校准模体 | |
CN113876346B (zh) | 一种倾斜图像的迭代校正方法 | |
CN113808211A (zh) | 一种用于图像引导***的几何校正模体及制造方法 |
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 |