CN110602990B - 用于锥形束计算机断层扫描的患者移动校正方法 - Google Patents

用于锥形束计算机断层扫描的患者移动校正方法 Download PDF

Info

Publication number
CN110602990B
CN110602990B CN201880028718.8A CN201880028718A CN110602990B CN 110602990 B CN110602990 B CN 110602990B CN 201880028718 A CN201880028718 A CN 201880028718A CN 110602990 B CN110602990 B CN 110602990B
Authority
CN
China
Prior art keywords
ray
projection
geometry
coordinate system
transformation
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
Application number
CN201880028718.8A
Other languages
English (en)
Other versions
CN110602990A (zh
Inventor
M·利尔佳
K·卡尔胡
J·拉海尔马
K·奈豪姆
A·希塔宁
T·米勒
S·克藤恩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Planmeca Oy
Original Assignee
Planmeca Oy
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Planmeca Oy filed Critical Planmeca Oy
Publication of CN110602990A publication Critical patent/CN110602990A/zh
Application granted granted Critical
Publication of CN110602990B publication Critical patent/CN110602990B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4064—Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
    • A61B6/4085—Cone-beams
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46—Arrangements for interfacing with the operator or the patient
    • A61B6/461—Displaying means of special interest
    • A61B6/466—Displaying means of special interest adapted to display 3D data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/44—Constructional features of apparatus for radiation diagnosis
    • A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
    • A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/501—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the head, e.g. neuroimaging or craniography
    • 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/10—Image acquisition modality
    • G06T2207/10116—X-ray image
    • G06T2207/10124—Digitally reconstructed radiograph [DRR]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/424—Iterative

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pulmonology (AREA)
  • Human Computer Interaction (AREA)
  • Quality & Reliability (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Neurology (AREA)
  • Neurosurgery (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

基于初始投影几何形状估计的迭代改进的用于锥形束计算机断层扫描的患者移动校正方法，该初始投影几何形状估计描述了X射线源和检测器在X射线投影图像的采集期间的空间位置和方向。取得的X射线投影图像组被用于在校正迭代期间计算中间CBCT重建，并且通过最大化测量的和重新投影的X射线投影图像的相似度来建立投影图像特定的校正几何变换。在投影图像特定的旋转坐标系中确定校正几何变换，该投影图像特定的旋转坐标系与X射线源和X射线检测器在X射线投影图像的采集期间的空间位置和方向相对应。使用校正的投影几何形状估计来计算出最终的CBCT重建。

Description

用于锥形束计算机断层扫描的患者移动校正方法
技术领域
本发明关注在锥形束计算机断层扫描领域的患者移动校正方法。
背景技术
以下公开涉及用于医疗应用的三维X射线锥形束计算机断层扫描，其中，从不同角度取得的多个X射线投影图像被用于重建患者的解剖结构的3-D截面图像。因为X射线管(X射线源)和传感器(X射线检测器)必须物理地经过与采集角度相对应的空间轨迹，所以X射线投影图像采集的持续时间通常大约是10秒-30秒。成像轨迹通常由旋转和平移机制来实现。
成像轨迹应以足够的精度而为人所知，并且被成像的物体在X射线投影图像采集期间应该保持足够地静止，以便重建的CBCT图像相对于解剖结构是清晰的并且真实的。因为投影图像测量被假定是表示静止物体的共同配准的整合图，这造成几何上一致的测量的集合，该几何上一致的测量的集合可以被用于重建表示该研究的解剖结构的衰减分布。每当这些假定被违反时，由于投影测量变得相互不一致，所以重建图像的精度下降。
最终图像的精度通常取决于在重建过程中做出的假定与实际的物理图像采集过程对应地有多准确。与每个取得的X射线图像相对应的X射线源和检测器的估计的空间位置影响在CBCT重建过程期间的射线路径的计算。由于成像设备的固有的制造和操作公差以及潜在的变形，所以实现的旋转角度和位置倾向于与根据成像轨迹的假定形式的理想值相偏离。然而可以通过使用在一定时间段或操作循环后重复进行的不同的校准方法来解决***偏差。
在医疗CBCT成像中，几何误差的最显著的来源是患者在X射线投影图像的采集期间的潜在移动。即，如果被成像的对象在X射线投影图像的采集期间移动，那么射线测量的有效空间路径变得相互不一致。虽然已知患者不应该移动并且常规地患者被指示不要移动，但是患者在X射线投影的采集期间通常不能保持完全地静止。通常通过支撑患者来解决这个问题。然而，过紧地支撑患者是不方便的且不舒服的。另外，防止所有患者移动将需要使用高度约束的支撑，这在常规成像中是不适用的。
在医疗CBCT成像中，在所有实际测量中上述的两种误差都一定程度地呈现出。在最坏的情况设想下，在放射科医生已经检查了图像重建的质量之后，投影图像测量的得到的几何不一致性可能甚至需要重复扫描。这是不期望的，因为旨在将与X射线图像采集相关的辐射剂量保持尽可能合理地低。
计算的方法已经被开发以解决计算机断层扫描成像中的几何误差的问题。在文献中报导的方法中，通过在固定坐标系中的刚性几何变换的方式，X射线源和X射线检测器的虚拟运动已经被应用于建模和补偿在投影图像采集期间的刚性运动。在与医疗CBCT成像相关的近期的方法中，通过最大化得到的CBCT重建的清晰度来使这样的几何变换最优化。通常这样的校正过程被迭代地执行。
在CBCT成像中，因为在CBCT成像中的固有的几何自由度不能被坐标系分离，所以将固定坐标系应用于患者运动的建模和补偿是不理想的。在CBCT成像设备中，特别的是X射线束发散并且形成金字塔形状的锥形。结果，沿着与X射线源和X射线检测器的中心邻接的等射线的位移将只影响倍率因子，然而沿着X射线检测器的像素阵列的平面内的位移将导致被成像的对象在其投影图像内的最大位移。而且，通过得到的校正的CBCT重建和未校正的CBCT重建的刚性配准的已知方式来防止由所应用的几何校正所引起的净变换在计算上是昂贵的，尤其是如果在几何校正过程中被重复应用
发明内容
在本公开的用于医疗的锥形束计算机断层扫描(CBCT)的患者移动校正过程中，利用附接到X射线源和检测器的空间位置和方向的固有的旋转坐标系的数据驱动算法被使用以建立与测量的X射线投影图像相对应的实际成像几何形状的改进估计。通过计算净变换并且直接从估计的校正几何变换参数中减去净变换来防止固定参考坐标系中的净变换。结果，基于X射线投影图像内容回顾性地估计和改进X射线源和检测器***相对于患者的相对位置，反之亦然。校正过程的目的是通过提高X射线图像测量的几何一致性来提高得到的CBCT重建图像质量，并且进而越来越好地满足断层图像重建的要求。固有坐标系使得能够以一种方式来定义应用的几何自由度，该方式与应用的几何自由度对校正过程以及得到的图像质量两者的相对重要性相对应。
本公开的患者移动校正过程将其输入作为用于计算CBCT图像重建通常所需的数据：X射线投影图像的集合和与X射线源和X射线检测器在X射线投影图像的采集期间的空间位置相对应的3-D投影几何形状的估计。类似于相应的方法，首先使用估计的投影几何形状来计算中间重建。然后，利用与X射线源和检测器在图像采集期间的虚拟移动相对应的变换来建立校正几何变换，该校正几何变换提高每个访问的X射线投影图像与其余的X射线投影图像的几何对应性。特别地，投影图像特定的旋转坐标系被应用以用于确定校正几何变换。校正过程可以被多次迭代，该校正过程包括使用投影几何形状的当前估计的中间重建的计算以及随后的校正变换的最优化。在建立校正几何变换之后，使用X射线投影图像和与校正变换相对应的投影几何形状的最终估计来计算最终CBCT重建。
在重建最终CBCT图像之前，本公开的患者移动校正过程的好处是投影几何一致性的提高，这导致在清晰度、细节水平和对比度方面有较高的图像质量。校正过程的另一好处是通过潜在地防止由于投影几何形状相关的图像质量下降而导致的重新扫描，可以减少对患者产生的辐射剂量。
附图说明
附图提供了本公开的患者移动校正方法的进一步理解。在图中：
图1示出一个典型CBCT成像装置的示例。
图2示出本公开的患者移动校正方法的流程图。
图3示出详述图2的步骤203的流程图。
图4示出在CBCT成像装置的背景下在本公开的患者移动校正方法中应用的旋转坐标系。
具体实施方式
为了计算3-D CBCT重建，使用围绕被成像的解剖结构旋转的X射线源和X射线检测器来测量若干2-D X射线投影图像。每个X射线投影图像的曝光通常发生在毫秒的范围内，并且当使用典型的扫描速度时，X射线源和检测器的连续的移动不会在投影图像中产生任何明显的运动模糊。为了使用诸如FDK算法之类的合适的方法来计算3-D重建，X射线源和检测器在包含成像或投影几何形状的合适的参考坐标系中的3-D空间位置必须可用于重建算法。
成像几何形状的典型的描述由X射线源的焦点和X射线检测器的中心的3-D位置以及足够唯一确定X射线检测器的方向的信息组成。这样的信息可以包括例如参考坐标系中的旋转角度，当这样信息根据预先定义的协定被应用时，将确定检测器的像素阵列的水平轴和垂直轴的方向。投影几何形状的描述通常是基于对成像设备的物理测量以及与给定的成像程序相对应的理想曝光轨迹的认识。另外，周期性校准过程通常被执行以确保投影几何形状的足够的精度。
在重建过程中，通过被成像的解剖结构的X辐射的空间传播被建模。通常，为了简单起见，假定直线传播，并且从源到检测器的X射线束路径被建模成线积分，该线积分基于被包含在投影几何形状描述中的信息来计算。使用合适的投影仪算法和投影几何形状描述来解决X射线路径和用于重建的3-D图像体素阵列的元素的空间重叠。基于空间重叠和X射线投影图像中的测量值，可以使用诸如著名的FDK算法之类的合适的方法来重建与被成像的解剖结构相对应的总X射线衰减分布。
然而，在重建过程中采用的基本假定是被成像的解剖结构在X射线投影图像的采集期间已保持足够地静止。基本原理是X射线投影应该表示静止对象的共同配准测量，然后该静止对象的共同配准测量可以被一致地组合以重建物体的3-D结构。可以由周期性校准过程来补偿投影几何形状中的***误差，该周期性校准过程消除了任何偏离假定的理想的X射线投影图像采集轨迹而产生的影响。
患者在成像期间的移动会产生更困难的问题，这是不可预知的并且本质上是随机的。虽然通过在成像期间支撑患者在一定程度上防止了该问题，但是CBCT成像中的X射线投影图像采集的相对长的持续时间(大约10秒)使得完全地消除患者移动是难以实现的。而且，某些患者(比如小孩和老人)更难以在成像期间保持静止。
万一在X射线投影图像采集期间出现明显的患者移动，那么明显的条纹状或模糊的伪影(artefacts)将继而产生在重建的图像中，这在最坏的情况下可能致使图像对于预期的医疗目是无用的。重复扫描可能会提供更好的结果，但是以额外的辐射剂量为代价。
在用于补偿患者移动的不利影响的本公开的回顾的方法中，首先使用2-D X射线投影图像和估计的成像几何形状来计算中间的CBCT重建。以比当为诊断目的而进行重建时通常使用的分辨率粗略的分辨率来计算中间的CBCT重建是足够的。中间重建用作集合来自所有X射线投影图像的信息，该中间重建具有反映测量的2-D X射线投影图像的相互几何一致性的外观。
投影图像特定的几何形状的最优化是基于测量物理X射线投影图像和对应的形成数字重建射线照片(DRR)的中间CBCT重建的重新投影数据的相似度。基本原理是当估计的投影几何形状是一致时，从定义上尝试满足由从正投影操作的意义上说的投影图像传达的测量的重建图像将产生与测量的数据非常匹配的重新投影。万一几何形状不匹配，中间重建的重新投影将偏离测量的投影。由于中间重建用作所有投影图像的集合，它的重新投影反映所有投影图像的总和以及测量的X射线投影图像的最大相似度，并且从中间重建的意义上说当对应的重新投影几何形状与平均正确投影几何形状匹配时，可以预期正投影操作将被最大化。
在典型的CBCT成像设备中，发射的X射线束发散并且形成金字塔形状的锥形。沿着与X射线源和X射线检测器的中心邻接的等射线的位移将只影响倍率因子，然而沿着X射线检测器的平面的位移将导致被成像的对象在其投影图像内的最大位移。从这个观点来看，本公开的方法采用旋转坐标系，该旋转坐标系在图像采集期间被附接到X射线源和检测器的物理位置。即，两个坐标轴被附接到矩形X射线检测器像素阵列，而剩余的垂直轴被附接到检测器的像素阵列的法线。通过限制每个投影图像沿着这些固定轴的投影几何形状的几何变换，可以分离几何自由度，从几何精度的意义来说几何自由度的重要性有所变化。
当每个X射线投影图像都经过单独的几何变换时，得到的平均变换可能将对重建图像产生净效应。这体现为例如重建的解剖结构相对于未校正的解剖结构的整体净位移或旋转。解剖结构的净变换可能将对重建图像的适用性产生不利影响。因此，还公开了用于净变换的补偿方法。可以通过将与每个变换的投影图像相对应的变换从旋转坐标轴映射到固定坐标轴来估计净变换。例如，通过假定每个投影图像仅沿着其水平轴移动，可以基于投影特定的旋转坐标系的已知的水平轴来计算在固定坐标系中的对应的位移，并且可以提取平均值以表示固定坐标系中的净位移。通过上述映射的逆，净变换的逆可以被映射回到旋转坐标系并且从投影图像特定变换中被减去。结果，固定坐标系中的净变换被消除。
在最优化过程中，由与变换的投影几何形状相对应的正投影图像和原始X射线投影图像的相似度来测量给出的X射线投影图像的校正几何变换的优度。例如由图像的均方差、相关系数或梯度相关系数可以测量重新投影和X射线投影的相似度。然后可以通过找到随着几何变换的参数变化的在正投影图像和X射线投影图像之间的相似度测量的极值来确定每个给出了中间CBCT重建的投影图像的最优(从相似度测量的意义上来说)几何变换。
然后最优化过程包括使用X射线投影几何形状的初始估计来计算中间CBCT重建；访问所有的测量的X射线投影图像或其子集；(对于每个被访问的投影图像)通过找到随着旋转坐标系中执行的几何变换的参数变化的在测量的X射线投影图像和中间CBCT重建的对应的正投影之间的最大相似度来建立投影图像特定的校正变换；估计固定参考坐标系中的净变换并且将对应的变换从旋转坐标系中的变换参数中减去；当估计已获得足够的校正结果时，计算最终的CBCT重建。
在图1中，示出医疗CBCT成像装置100的示例，该医疗CBCT成像装置100包括垂直底座结构101、患者支撑部件107以及臂部103，其中支撑结构102从垂直底座结构101水平地延伸，臂部103支撑支持成像部件的结构，臂部104。针对支撑成像部件的臂部104，布置彼此相距一定距离的X射线源105和X射线图像信息的接收器部件109(X射线检测器)，X射线源105和X射线图像信息的接收器部件109相对于患者支撑部件107布置使得放置在X射线源105和X射线图像信息的接收器部件109之间的成像站108被形成使得由X射线源105生成的束可对准以穿过成像站108朝向X射线图像信息的接收器部件109。支撑成像部件的臂部104被布置成可旋转，并且相对于支撑它的结构103和/或患者支撑站108，臂部104的位置可以被布置成可变化的。该布置包括控制部件，关于该控制部件图1示出与支撑患者支撑部件107的支撑结构102连接放置的控制面板106。成像装置100可以被布置成经由线缆连接到控制器110，该控制器包括布置有用于处理由成像装置产生的图像信息的部件的计算机，以及可以显示图像的显示器111。控制器110还包括至少一个处理器112和至少一个存储器113。至少一个处理器112可以被配置成执行计算程序，并且至少一个存储器113被配置成存储计算机程序和相关数据。控制器110可以是用于实现以下描述的过程的通用计算机或者特别制造的设备。
图2描述了本公开的患者移动校正方法的步骤，该方法可以被用于例如处理由图1的成像布置取得的图像。该方法是基于找到在被附接到X射线源105和X射线检测器109的位置的旋转坐标系中的初始估计的投影几何形状的校正几何变换。
在图2中的方法的步骤200中，取得输入数据，该输入数据包含测量的X射线投影图像和与采集过程相对应的投影几何形状的初始估计。在图1的装置中，当X射线源105和X射线检测器109借助于臂103和104围绕成像站108旋转和平移时，投影几何形状限定了X射线源105和X射线检测器109的物理轨迹。投影几何形状还确定X射线检测器109的估计方向。投影几何形状的形式通常是基于利用例如具有不透射线标志的已知的参考模型的几何形状校准程序。
在步骤201中，使用步骤200中取得的输入数据来计算为了投影最优化方法的目的的中间CBCT重建。可以使用输入数据的降采样版本，因为通常不需要应用与用于诊断目的的CBCT图像中一样的高空间分辨率。中间重建被理解成集合在步骤200的测量过程期间取得的所有可用的物理信息和几何信息。万一几何不一致，这由中间重建反映出，例如由重建细节的模糊性反映出。
在步骤202中，在旋转坐标系中建立投影图像的校正几何变换。几何变换的目的是补偿在投影几何形状的初始估计中的固有几何不一致。下面参考图3说明步骤202的细节，但是更通常地，通过找到最优几何变换来寻找校正几何变换。变换的优度通过为其分配相似度值来定义。通过将中间CBCT图像重建的重新投影的数字重建射线照片(DRR)与对应的测量的X射线投影图像进行比较来计算出相似度值，其中应用的投影几何形状对应于评估的几何变换。较高的相似度值被视为较好的校正几何变换的指示。应用用于几何变换的旋转坐标系使得能够在锥形束投影几何形状中根据几何自由度对问题的重要性来分离几何自由度。
在步骤203中，减去固定坐标系中的净几何变换。通常参考成像装置100的静态组件(比如与成像站108连接而被固定的成像装置的组件)来定义固定坐标系。通过将在步骤202中建立的几何变换从旋转坐标系线性地变换到固定坐标系来计算净几何变换。基于旋转坐标系和固定坐标系的已知的坐标轴容易地获得线性变换。在净变换被建立之后，净变换的逆从固定坐标系被线性地变换到其旋转坐标系。于是净变换的逆与旋转坐标系中的每个投影图像的几何变换参数的集合相对应。将这些值添加到在步骤202中建立的几何参数值中将导致净变换在固定坐标系中被抵消。
在步骤204中，从步骤203得到的变换参数值被应用到初始投影几何形状估计以获得校正的投影几何形状估计。特别地，在前进到步骤204之前，步骤201-203可以被迭代地重复多次。
在步骤205中，使用校正的投影几何形状估计来计算最终的CBCT重建。除了步骤201-204未被应用的情况之外，以由在步骤201-204中获得的校正的投影几何形状估计来替换初始投影几何形状估计的常规的方式来计算最终的CBCT重建。
图3作为示例描述了根据图2的步骤202的细节。在步骤300中，评估的几何变换被应用到旋转坐标系中的投影图像的初始几何形状。再次，坐标系被定义成与X射线源105和X射线检测器109在被考虑的投影图像的物理采集期间的空间位置和方向一致。在旋转参考系中应用变换涉及从通常定义成像几何形状的固定坐标系到执行变换的旋转坐标系的线性映射，接着是从旋转坐标系到固定坐标系的逆线性映射。作为实现的效果的简单示例，为了沿着X射线检测器109的水平轴来应用平移，投影图像的初始投影几何形状被映射到旋转坐标系并且沿着旋转坐标系的该轴以与X射线源和检测器的虚拟移动相对应的给定量被平移，并且然后被映射回到固定坐标系。在固定坐标系中表示变换的投影几何形状，并将获得的变换的投影几何形状用作后续步骤301的输入。
在步骤301中，使用在步骤300中获得的变换的投影几何形状来计算在步骤201中计算的中间CBCT重建的重新投影DRR图像。可以使用诸如Siddon射线投射法之类的标准算法来执行DRR图像的计算。在典型的重新投影算法中，算法的输入由3-D X射线表示的端点和其DRR图像被计算的源图像组成，该源图像包括在相同坐标系中表示的它的空间位置和方向的知识。在所描述的设置中，该坐标系与固定坐标系相对应。在步骤300中应用的几何变换的效果是改变通过中间CBCT重建的每个虚拟X射线路径的端点，这将几何变换的效果传播到获得的DRR图像。
在步骤302中，评估在步骤301中获得的DRR图像与步骤200中获得的X射线投影图像之间的相似度。相似度是基于使用已建立的方法(比如图像的平均平方差或它们的互相关)对图像进行逐点比较。用于评估相似度的具体的测量对所描述的方法并不重要。获得的相似度值被分配给作为步骤300的输入而给出的几何参数。较高的相似度被认为是更合适的几何变换参数的指示。
在步骤303中，通过找到与通过应用步骤300-302获得的最高相似度值相对应的参数来建立最优几何变换参数。在典型的设置中，由诸如著名的Nelder-Mead单纯型算法之类的合适的最小化算法来重复地评估步骤300-302，以建立最优几何变换参数。存储这些参数并且将这些参数分配到特定的投影图像，直到这些参数通过步骤201-203的重复被潜在地改变为止。
图4示出在步骤202中应用的旋转坐标系。在CBCT成像装置中，C-臂部403通常支撑X射线源404和X射线检测器405。C-臂由垂直基座结构401和肩臂部402支撑。在CBCT成像中，特别的是X射线束发散并且形成金字塔形状的锥形。沿着与X射线源404和X射线检测器405的中心邻接的等射线的偏移只会影响倍率因子，然而X射线检测器的像素阵列的平面内的偏移将导致被成像的对象在它的投影图像中的最大位移。从这个观点来看，本公开的方法采用旋转的uvw坐标系406，该旋转的uvw坐标系406在图像采集期间被附接到X射线源和检测器像素阵列的物理位置和方向。通过限制每个投影图像的沿着这些固有轴的投影几何形状的几何变换，可以分离几何自由度，从几何精度的意义上来说该几何自由度重要性会有所变化。计算净变换的固定坐标系由xyz坐标系407表示。
将坐标系的u轴设置为与X射线检测器的法线相一致，并且将v轴设置为与X射线检测器的像素阵列的水平轴相一致，例如，这些自由度是用于刚性变换：
1、沿X射线检测器的像素阵列的法线(u轴)的纵向位移；2、沿X射线检测器的像素阵列的水平轴(v轴)的横向位移；3、沿X射线检测器的像素阵列的垂直轴(w轴)的垂直位移；4、绕u轴的旋转(滚动角)；5、绕v轴的旋转(俯仰角)；6、绕w轴的旋转(偏航角)。
本公开的患者移动校正方法的结果是与X射线投影图像的物理采集和相应地CBCT重建图像相对应的投影几何形状的改进的估计，其中几何不一致的影响被减小。
本公开的方法可以被实现作为在计算设备中执行的计算机软件。该软件在计算机可读介质上被具体化，以便该软件可以被提供到诸如图1的控制器110之类的计算设备。
如上所述，示例实施例的组件可以包括用于保存根据本实施例的教导而被编程的指令并且用于保存本文中描述的数据结构、表、记录和/或其它数据的计算机可读介质或存储器。计算机可读介质可以包括参与将指令提供给处理器以用于执行的任何合适的介质。计算机可读介质的常见形式可以包括，例如，软盘、柔性盘、硬盘、磁带、任何其它合适的磁介质、CD-ROM、CD±R、CD±RW、DVD、DVD-RAM、DVD±RW、DVD±R、HD DVD、HD DVD-R、HD DVD-RW、HD DVD-RAM、蓝光光盘、任何其它合适的光学介质、RAM、PROM、EPROM、闪存-EPROM、任何其它合适的存储器芯片或盒带、载波或计算机可以读取的任何其它合适的介质。
对于本领域的技术人员清楚的是随着技术的进步，可以以各种方式来实现自校准医学成像装置的基本想法。因此，自校准医学成像装置和其实施例不限于上述示例；反而，它们可以在权利要求的范围内变化。

Claims (9)

1.一种用于锥形束计算机断层扫描的患者移动校正方法，包括：
使用包含X射线源和X射线检测器的X射线成像部件来采集被成像的患者的一组X射线投影图像；
定义初始投影几何形状估计，该初始投影几何形状估计描述X射线源和X射线检测器在X射线投影图像的采集期间的空间位置和方向；
使用X射线投影图像和初始投影几何形状估计来计算中间CBCT重建；
基于X射线投影图像、初始投影几何形状估计以及中间CBCT重建的重新投影的数字重建射线照片来确定用于初始投影几何形状估计的投影图像特定的校正几何变换；
并且使用X射线投影图像和与确定的校正几何变换相对应的投影几何形状的校正的估计来计算最终CBCT重建；
其特征在于，初始投影几何形状估计的校正几何变换在投影图像特定的旋转坐标系中被确定，该投影图像特定的旋转坐标系与X射线源和X射线检测器在X射线投影图像的采集期间的空间位置和方向相对应。
2.根据权利要求1所述的方法，其中，通过将两个垂直轴设置为与X射线检测器的像素阵列相一致，并且将第三垂直轴设置为与X射线检测器的像素阵列的法线相一致来定义所述投影图像特定的旋转坐标系。
3.根据权利要求1或2所述的方法，其中，通过将旋转坐标系中的几何变换应用到初始投影几何形状估计来获得重新投影的数字重建射线照片的投影几何形状。
4.根据权利要求1或2所述的方法，其中，通过找到在X射线投影图像与中间CBCT重建的对应的重新投影的数字重建射线照片之间的最大相似度来确定校正几何变换。
5.根据权利要求1或2所述的方法，其中，固定坐标系中的净变换被计算并且从旋转坐标系中确定的校正几何变换中被减去。
6.根据权利要求1或2所述的方法，其中，通过将投影图像特定的几何变换从旋转坐标系映射到固定坐标系并且取投影图像特定的几何变换的平均值来计算固定坐标系中的净变换，并且在从固定坐标系映射到旋转坐标系之后该平均值从几何变换中被减去。
7.根据权利要求1或2所述的方法，其中，中间CBCT重建的计算和确定校正几何变换被多次迭代，伴随着在每次迭代之后将校正的投影几何形状估计用作用于下面的迭代的新的初始投影几何形状估计，并且最终校正的投影几何形状估计与最后应用的迭代的校正的投影几何形状估计相对应。
8.一种计算机可读存储介质，所述计算机可读存储介质存储计算机可执行指令，该计算机可执行指令在由处理器执行时使得处理器执行根据前述权利要求1至7中的任一项所述的方法。
9.一种装置，包括：
至少一个处理器(112)，被配置成执行计算机程序；以及
至少一个存储器(113)，被配置成存储计算机程序和相关数据；
其特征在于，该装置能够连接到医疗成像装置并且被配置成执行根据权利要求1至7中的任一项所述的方法。
CN201880028718.8A 2017-03-17 2018-03-19 用于锥形束计算机断层扫描的患者移动校正方法 Active CN110602990B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FI20175244 2017-03-17
FI20175244A FI20175244L (fi) 2017-03-17 2017-03-17 Itsekalibroiva lääketieteellinen kuvannuslaite
PCT/FI2018/000005 WO2018167355A1 (en) 2017-03-17 2018-03-19 Patient movement correction method for cone-beam computed tomography

Publications (2)

Publication Number Publication Date
CN110602990A CN110602990A (zh) 2019-12-20
CN110602990B true CN110602990B (zh) 2023-09-01

Family

ID=62143215

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201880028718.8A Active CN110602990B (zh) 2017-03-17 2018-03-19 用于锥形束计算机断层扫描的患者移动校正方法

Country Status (9)

Country Link
US (1) US10636179B2 (zh)
EP (1) EP3595528B1 (zh)
JP (1) JP7116078B2 (zh)
KR (1) KR102369653B1 (zh)
CN (1) CN110602990B (zh)
BR (1) BR112019019206A2 (zh)
FI (1) FI20175244L (zh)
RU (1) RU2766743C1 (zh)
WO (1) WO2018167355A1 (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3639752A1 (en) * 2018-10-03 2020-04-22 Canon Medical Systems Corporation Analyzing apparatus and analyzing program
CN110766629B (zh) * 2019-10-17 2022-03-01 广州华端科技有限公司 Cbct***几何校正方法、装置、计算机设备和存储介质
CN111553849B (zh) * 2020-03-28 2023-03-24 中国人民解放军战略支援部队信息工程大学 基于局部特征匹配的锥束ct几何伪影去除方法及装置
CN113520421A (zh) * 2020-05-06 2021-10-22 北京友通上昊科技有限公司 X射线切面成像方法、存储介质及成像***
CN112529820B (zh) * 2020-12-29 2024-06-14 广州开云影像科技有限公司 Cbct***半扇束几何校正方法、装置、设备和介质
US11712215B2 (en) 2021-04-13 2023-08-01 Canon Medical Systems Corporation Devices, systems, and methods for motion-corrected medical imaging
EP4120182A1 (de) * 2021-07-12 2023-01-18 DENTSPLY SIRONA Inc. Selbstkalibrierendes dentales dvt unterstützt durch maschinelles lernen
CN113963056B (zh) * 2021-09-07 2022-08-26 于留青 Ct图像重建方法、装置、电子设备以及存储介质
CN116894798B (zh) * 2023-09-11 2023-12-05 金华飞光科技有限公司 光固化3d打印机的投影畸形矫正方法及***

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5251128A (en) * 1990-11-19 1993-10-05 General Electric Company Motion artifact reduction in projection imaging
CN101111758A (zh) * 2005-02-01 2008-01-23 皇家飞利浦电子股份有限公司 用于x射线投影的校正或扩展的设备和方法
CN102488528A (zh) * 2011-12-07 2012-06-13 华中科技大学 一种层析成像几何参数的校准方法
CN102525501A (zh) * 2010-11-22 2012-07-04 佳能株式会社 图像处理设备和图像处理方法
CN103519833A (zh) * 2013-06-05 2014-01-22 东南大学 一种旋转c型臂x射线机的三维校正重建方法
CN103729827A (zh) * 2012-10-11 2014-04-16 南京普爱射线影像设备有限公司 影像增强器c形臂x射线机三维重建重叠伪影校正方法
CN103735282A (zh) * 2014-01-06 2014-04-23 北京大学 一种锥束ct***探测器几何校正装置及其校正方法
CN104334081A (zh) * 2012-06-05 2015-02-04 皇家飞利浦有限公司 X射线ct成像器的运动层分解校准

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ES2254599T3 (es) * 2002-06-17 2006-06-16 Nucletron B.V. Sistema de planificacion de tratamiento por radiacion en tiempo real.
US6888924B2 (en) * 2003-06-25 2005-05-03 General Electric Company Method, apparatus, and medium for calibration of tomosynthesis system geometry using fiducial markers with non-determined position
US7653229B2 (en) * 2003-12-23 2010-01-26 General Electric Company Methods and apparatus for reconstruction of volume data from projection data
JP4438053B2 (ja) * 2004-05-11 2010-03-24 キヤノン株式会社 放射線撮像装置、画像処理方法及びコンピュータプログラム
DE102004040629B4 (de) * 2004-08-19 2011-05-19 MedCom Gesellschaft für medizinische Bildverarbeitung mbH Einstellen einer Relativposition einer Bestrahlungseinrichtung und eines Patienten
DE102005005087B4 (de) * 2005-02-03 2012-09-27 Siemens Ag Verfahren zur Korrektur nicht reproduzierbarer Geometriefehler eines Röntgen-C-Bogengerätes
US8199168B2 (en) * 2005-11-15 2012-06-12 General Electric Company System and method for 3D graphical prescription of a medical imaging volume
US7894649B2 (en) * 2006-11-02 2011-02-22 Accuray Incorporated Target tracking using direct target registration
CN101553207B (zh) * 2006-11-09 2013-01-16 爱尔康研究有限公司 包含水不溶性聚合物基质的泪点塞
FR2919096A1 (fr) * 2007-07-19 2009-01-23 Gen Electric Procede de correction de recalage d'images radiographiques
EP2295113B1 (en) * 2009-03-10 2020-02-26 Hitachi, Ltd. Controller of radiation therapy equipment
WO2011085448A1 (en) * 2010-01-13 2011-07-21 The Australian National University A computed tomography imaging process and system
JP2015525665A (ja) * 2012-08-20 2015-09-07 オレンジデンタル・ゲゼルシャフト・ミット・ベシュレンクテル・ハフツング・ウント・コムパニー・コマンディットゲゼルシャフトｏｒａｎｇｅｄｅｎｔａｌ ＧｍｂＨ ＆ Ｃｏ． ＫＧ コーンビームコンピュータートモグラフィー装置の幾何学的特徴付け及び較正
WO2015028069A1 (de) * 2013-08-29 2015-03-05 Siemens Aktiengesellschaft Computertomographie mit bewegungskompensation
EP3175790B1 (en) * 2013-11-04 2021-09-08 Ecential Robotics Method for reconstructing a 3d image from 2d x-ray images
US10478148B2 (en) * 2016-02-19 2019-11-19 The Johns Hopkins University Self-calibrating projection geometry for volumetric image reconstruction

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5251128A (en) * 1990-11-19 1993-10-05 General Electric Company Motion artifact reduction in projection imaging
CN101111758A (zh) * 2005-02-01 2008-01-23 皇家飞利浦电子股份有限公司 用于x射线投影的校正或扩展的设备和方法
CN102525501A (zh) * 2010-11-22 2012-07-04 佳能株式会社 图像处理设备和图像处理方法
CN102488528A (zh) * 2011-12-07 2012-06-13 华中科技大学 一种层析成像几何参数的校准方法
CN104334081A (zh) * 2012-06-05 2015-02-04 皇家飞利浦有限公司 X射线ct成像器的运动层分解校准
CN103729827A (zh) * 2012-10-11 2014-04-16 南京普爱射线影像设备有限公司 影像增强器c形臂x射线机三维重建重叠伪影校正方法
CN103519833A (zh) * 2013-06-05 2014-01-22 东南大学 一种旋转c型臂x射线机的三维校正重建方法
CN103735282A (zh) * 2014-01-06 2014-04-23 北京大学 一种锥束ct***探测器几何校正装置及其校正方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
袁翠云 ; 齐宏亮 ; 陈梓嘉 ; 吴书裕 ; 徐圆 ; 周凌宏 ; .基于投影域校正的CT图像环形伪影去除方法.计算机工程与设计.(第03期),全文. *

Also Published As

Publication number Publication date
US20180268574A1 (en) 2018-09-20
WO2018167355A1 (en) 2018-09-20
KR20190126882A (ko) 2019-11-12
US10636179B2 (en) 2020-04-28
EP3595528B1 (en) 2023-07-19
RU2766743C1 (ru) 2022-03-15
EP3595528A1 (en) 2020-01-22
BR112019019206A2 (pt) 2020-04-14
JP2020509887A (ja) 2020-04-02
EP3595528C0 (en) 2023-07-19
JP7116078B2 (ja) 2022-08-09
FI20175244L (fi) 2018-09-18
CN110602990A (zh) 2019-12-20
KR102369653B1 (ko) 2022-03-07

Similar Documents

Publication Publication Date Title
CN110602990B (zh) 用于锥形束计算机断层扫描的患者移动校正方法
JP5906015B2 (ja) 特徴に基づいた２次元／３次元画像のレジストレーション
US7782998B2 (en) Method and apparatus for correcting motion in image reconstruction
JP5572221B2 (ja) Ｘ線画像内のマーカ位置からの患者運動ベクトルの抽出
US8565856B2 (en) Ultrasonic imager for motion measurement in multi-modality emission imaging
JP5600742B2 (ja) 投射のデータを処理するための装置及び方法
EP2036038A1 (en) Method and system for error compensation
JP2006167455A (ja) イメージング・システムの取得時幾何学的構成のための方法及び装置
US20170340287A1 (en) Method And Apparatus For Motion Correction In CT Imaging
JP7436463B2 (ja) デジタル***トモシンセシス撮影装置の自己較正手順
Niebler et al. Projection‐based improvement of 3D reconstructions from motion‐impaired dental cone beam CT data
JP4444100B2 (ja) 多次元構造の解析方法
EP3331449B1 (en) Computed tomography image generation apparatus
RU2789105C2 (ru) Способ автокалибровки устройства для формирования изображений цифровой томографической реконструкции груди
CN113963056B (zh) Ct图像重建方法、装置、电子设备以及存储介质
Wicklein et al. An object-independent measure for improving misalignment correction in c-arm ct
JP2002533146A (ja) 対象物の細部の寸法を得る方法
CN114119801A (zh) 三维数字减影血管造影方法、装置、电子设备及存储介质
WO2023285202A1 (en) Improved metal artifact reduction through motion artifact compensation
RU2021108170A (ru) Способ автокалибровки устройства для формирования изображений цифровой томографической реконструкции груди

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
REG Reference to a national code

Ref country code: HK

Ref legal event code: DE

Ref document number: 40018965

Country of ref document: HK

GR01 Patent grant
GR01 Patent grant