CN109671128A - 图像重建过程中的数据处理、图像重建方法和装置 - Google Patents

图像重建过程中的数据处理、图像重建方法和装置 Download PDF

Info

Publication number
CN109671128A
CN109671128A CN201811491822.9A CN201811491822A CN109671128A CN 109671128 A CN109671128 A CN 109671128A CN 201811491822 A CN201811491822 A CN 201811491822A CN 109671128 A CN109671128 A CN 109671128A
Authority
CN
China
Prior art keywords
steel ball
radiographic source
coordinate
detector
ray
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
Application number
CN201811491822.9A
Other languages
English (en)
Other versions
CN109671128B (zh
Inventor
李翰威
吴书裕
齐宏亮
陈宇思
骆毅斌
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.)
Zhuhai Ruitong Medical Equipment Co.,Ltd.
Original Assignee
Guangzhou Huarui Technology Co Ltd
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 Guangzhou Huarui Technology Co Ltd filed Critical Guangzhou Huarui Technology Co Ltd
Priority to CN201811491822.9A priority Critical patent/CN109671128B/zh
Publication of CN109671128A publication Critical patent/CN109671128A/zh
Application granted granted Critical
Publication of CN109671128B publication Critical patent/CN109671128B/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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30—Computing systems specially adapted for manufacturing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Image Processing (AREA)

Abstract

本发明涉及一种图像重建过程中的数据处理、图像重建方法和装置、计算机设备、计算机存储介质。上述图像重建过程中的数据处理方法包括：在射线源发射射线时，获取钢珠在探测器上的钢珠投影坐标；其中，所述钢珠固定设置在校正模体上；根据钢珠在所述校正模体上的钢珠空间坐标以及所述钢珠投影坐标确定DBT系统的投影矩阵；其中，所述投影矩阵为表征所述射线通过所述钢珠照射至所述探测器时，所述钢珠空间坐标与所述钢珠投影坐标之间变换关系的矩阵；根据所述投影矩阵计算所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标；根据所述垂直距离和所述垂点坐标计算所述DBT系统在所述射线源当前发射射线时的贡献值。

Description

图像重建过程中的数据处理、图像重建方法和装置
技术领域
本发明涉及图像处理技术领域，特别是涉及一种图像重建过程中的数据处理、图像重建方法和装置、计算机设备、计算机存储介质。
背景技术
随着计算机技术的不断发展，数字化层析成像技术(Digital Tomosynthesis，DTS)已逐步应用于临床，与CT(computed tomography)相比其特点为高效率、高检出率、低剂量。与传统DR(Digital Radiography)相比，DTS可进行冠状面任意层面和层厚重建，获得单层的影像，消除重叠影的干扰，明确病灶，目前已在乳腺、骨科、泌尿系统和呼吸系统诊断上得到应用。层析成像技术的是原理为体层合成技术(Tomosynthesis)，是一种通过在有限角度范围内采集有限数量的投影数据并进行重建得到一个近似三维断层的成像技术，扫描方式如图1所示(以数字乳腺层析成像Digital Breast Tomosynthesis,DBT为例进行说明)。DTS的采集轨道不局限于圆弧轨道，也可能是直线轨道，其通过球管、病人、探测器的相对运动进而获取不同方向上的投影图像，即可实现对应层面图像的重建。实现DTS必须具有良好的硬件基础，包括对机械安装和球管旋转高精度的要求。在理想的情况下，射线源的中心射线应经过旋转轴并垂直地落在平板探测器的中心位置。但是，在DTS设计和安装调试过程中，由于机械设计和安装水平等原因，导致系统的几何结构和理想的几何结构之间存在一定程度的几何偏移，导致重建图像中产生几何伪影。由于几何伪影降低图像的对比度和分辨率，严重影响重建图像的质量。为了提高重建图像的质量，需要对DTS系统中的几何参数进行校准，进而实现几何伪影的校正。
现有的几何校正方法包括基于图像驱动的几何校正算法和基于校正模体的几何校正算法。基于校正模体的校正方法是现在主流的几何校正方法，多是使用专门设计的校正模体，几何模体的形状各式各样，上述几何模体可以镶嵌着钢珠点。因为钢珠点的衰减系数很大，在图像上的能够清晰地分割并定位出其坐标，通过几何校正算法即可获取系统的几何参数，之后可利用这些参数进行后期的图像重建。
目前，基于校正模体的校正方法都是先扫描校正模体得出几何参数，之后再利用这些几何参数对物体进行重建。在具体的图像重建过程中，射线源在绕轨道旋转时不可避免的会产生随机抖动，如此每个角度下的实际几何参数与之前几何校正获得的几何参数难以保持一致，因而依据上述几何参数所重建的图像存在效果差的问题。
发明内容
基于此，有必要针对依据传统方案所重建的图像存在效果差的技术问题，提供一种图像重建过程中的数据处理、图像重建方法和装置、计算机设备、计算机存储介质。
一种图像重建过程中的数据处理方法，包括：
在射线源发射射线时，获取钢珠在探测器上的钢珠投影坐标；其中，所述钢珠固定设置在校正模体上，所述校正模体固定设置在所述射线源发射射线的一侧；
根据钢珠在所述校正模体上的钢珠空间坐标以及所述钢珠投影坐标确定DBT系统的投影矩阵；其中，所述投影矩阵为表征所述射线通过所述钢珠照射至所述探测器时，所述钢珠空间坐标与所述钢珠投影坐标之间变换关系的矩阵；
根据所述投影矩阵计算所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标；
根据所述垂直距离和所述垂点坐标计算所述DBT系统在所述射线源当前发射射线时的贡献值。
在一个实施例中，所述校正模体为两个底面分别固定设置多个钢珠的圆柱体，所述钢珠分别在圆柱体两个底面呈圆形排列；
所述在射线源发射射线时，获取钢珠在探测器上的钢珠投影坐标的过程之前，还包括：
将所述校正模体固定设置在所述射线源发射射线的一侧，所述校正模体的中心点与所述射线源的准直器中心点在所述探测器的一条垂线上。
在一个实施例中，所述获取钢珠在探测器上的钢珠投影坐标的过程包括：
在探测器获得的投影数据中提取灰度值大于灰度阈值的钢珠投影点，将所述钢珠投影点的坐标确定为钢珠投影坐标。
在一个实施例中，所述根据所述投影矩阵计算所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标的过程包括：
对所述投影矩阵进行RQ分解，获得所述DBT系统的内矩阵，根据所述内矩阵的矩阵元素分别确定所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标。
作为一个实施例，所述垂直距离为：SID＝K11du；
所述垂点坐标为：u0＝K13，v0＝K23；
其中，SID表示垂直距离，K11表示内矩阵K中第一行第一列元素，du表示探测器U轴的像素点尺寸大小，(u0,v0)表示垂点坐标，K13表示内矩阵K中第一行第三列元素，K23表示内矩阵K中第二行第三列元素。
一种图像重建方法，包括如下步骤：
分别在射线源各次发射射线的过程中，根据如上任意一个实施例所述的图像重建过程中的数据处理方法计算所述DBT系统在射线源该次发射射线时的贡献值，并获取所述射线源该次发射射线时，探测器探测得到的投影数据；
根据所述射线源各次发射射线分别产生的贡献值和投影数据进行图像重建。
一种图像重建过程中的数据处理装置，包括：
第一获取模块，用于在射线源发射射线时，获取钢珠在探测器上的钢珠投影坐标；其中，所述钢珠固定设置在校正模体上，所述校正模体固定设置在所述射线源发射射线的一侧；
第一确定模块，用于根据钢珠在所述校正模体上的钢珠空间坐标以及所述钢珠投影坐标确定DBT系统的投影矩阵；其中，所述投影矩阵为表征所述射线通过所述钢珠照射至所述探测器时，所述钢珠空间坐标与所述钢珠投影坐标之间变换关系的矩阵；
第一计算模块，用于根据所述投影矩阵计算所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标；
第二计算模块，用于根据所述垂直距离和所述垂点坐标计算所述DBT系统在所述射线源当前发射射线时的贡献值。
一种图像重建装置，包括：
第三计算模块，用于分别在射线源各次发射射线的过程中，根据如上任意一个实施例所述的图像重建过程中的数据处理装置计算所述DBT系统在射线源该次发射射线时的贡献值，并获取所述射线源该次发射射线时，探测器探测得到的投影数据；
重建模块，用于根据所述射线源各次发射射线分别产生的贡献值和投影数据进行图像重建。
一种计算机设备，包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述任一实施例提供的图像重建过程中的数据处理方法或者图像重建方法。
一种计算机存储介质，其上存储有计算机程序，该程序被处理器执行时实现上述任一实施例提供的图像重建过程中的数据处理方法或者图像重建方法。
上述图像重建过程中的数据处理方法、装置、计算机设备和计算机存储介质，可以在射线源发射射线时，获取钢珠在探测器上的钢珠投影坐标，以根据钢珠在校正模体上的钢珠空间坐标以及上述钢珠投影坐标确定DBT系统的投影矩阵，计算所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标，以根据上述垂直距离和垂点坐标计算所述DBT系统在所述射线源当前发射射线时的贡献值，使所确定的贡献值以射线源发射产生相应投影数据的射线时相对于探测器的垂直距离和垂点坐标为依据，具有较高的准确性，提高了依据上述贡献值进行相应图像重建的效果。
附图说明
图1为传统方案中DBT的成像扫描方式示意图；
图2为一个实施例的图像重建过程中的数据处理方法流程图；
图3为一个实施例的几何映射模型示意图；
图4为一个实施例的校正模体示意图；
图5为一个实施例的射线源与探测器之间的结构示意图；
图6为一个实施例的重建图像示意图和灰度直方图；
图7为一个实施例的图像重建过程中的数据处理装置结构示意图；
图8为一个实施例的重建图像示意图；
图9为一个实施例的计算机设备内部结构示意图。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步的详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本发明，并不限定本发明的保护范围。
需要说明的是，本发明实施例所涉及的术语“第一\第二\第三”仅仅是区别类似的对象，不代表针对对象的特定排序，可以理解地，“第一\第二\第三”在允许的情况下可以互换特定的顺序或先后次序。应该理解“第一\第二\第三”区分的对象在适当情况下可以互换，以使这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。
本发明实施例的术语“包括”和“具有”以及它们任何变形，意图在于覆盖不排他的包含。例如包含了一系列步骤或模块的过程、方法、装置、产品或设备没有限定于已列出的步骤或模块，而是可选地还包括没有列出的步骤或模块，或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤或模块。
在本文中提及“实施例”意味着，结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例，也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是，本文所描述的实施例可以与其它实施例相结合。
在本文中提及的“多个”是指两个或两个以上。“和/或”，描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。
参考图2所示，图2为一个实施例的图像重建过程中的数据处理方法流程图，包括：
S10，在射线源发射射线时，获取钢珠在探测器上的钢珠投影坐标；其中，所述钢珠固定设置在校正模体上，所述校正模体固定设置在所述射线源发射射线的一侧；
上述步骤S10可以在射线源发射射线时，获取成像对象(压迫对象)和校正模体的投影图像，上述投影图像包括钢珠在探测器上的钢珠投影坐标。
上述射线源为DBT系统的射线源，射线源在可以在探测器所在的方向旋转，每隔设定角度发射一次射线，使探测器可以获得校正模体以及压迫板所压迫的压迫对象(如乳腺组织)受射线照射后投影至探测器上的投影数据。上述设定角度分别可以依据重建图像的所需精度确定，比如设定角度可以设置为3°。
在射线源发射射线时，探测器可以获得校正模体和压迫对象对应的投影图像，投影图像所包括的数据为相应的投影数据。具体地，上述投影图像上的高亮点为钢珠对应的钢珠投影点，上述钢珠投影点的坐标为钢珠投影坐标。
上述校正模体可以为固定设置钢珠的圆柱体，上述钢珠可以分别设置在圆柱体的两个底面。上述圆柱体上钢珠中心点的坐标为钢珠空间坐标，钢珠投影坐标为钢珠中心点对应的投影坐标。
S20，根据钢珠在所述校正模体上的钢珠空间坐标以及所述钢珠投影坐标确定DBT系统的投影矩阵；其中，所述投影矩阵为表征所述射线通过所述钢珠照射至所述探测器时，所述钢珠空间坐标与所述钢珠投影坐标之间变换关系的矩阵；
上述投影矩阵为表征所述射线通过所述钢珠照射至所述探测器时，所述钢珠空间坐标与所述钢珠投影坐标之间变换关系的矩阵，投影矩阵也可以表征上述射线经过的空间点与该空间点在探测器上的投影点之间的变换关系。比如参考图3所示，射线源S发射的射线通过空间中的一点(x,y,z)在探测器上形成投影(u,v)，上述投影矩阵可以表征上述空间点(x,y,z)与相应投影点(u,v)之间的变换关系。扫描对象(包括校正模体和压迫对象)的投影过程中，被射线照射的空间点与相应空间点的投影点之间的对应关系为：
其中，w表示尺度因子，(u,v,1)表示图像坐标系下的投影点的齐次坐标，即为经过像素离散化后的投影点的坐标，P表示投影矩阵，上述P可以为3×4的透视投影矩阵，(x,y,z,1)表示模体坐标系下的空间点的坐标。上述坐标可以为相应点的齐次坐标系中的坐标，上述齐次坐标系可以将N维向量用N+1维向量表示，以描述投影几何，可将旋转平移整合成一个矩阵进行计算。
在扫描过程中，校正模体上设置有多个钢珠，各个钢珠分别具有相应的钢珠空间坐标和钢珠投影坐标，上述投影矩阵可以依据校正模体上述各个钢珠的钢珠空间坐标，以及相应的钢珠投影坐标确定。
S30，根据所述投影矩阵计算所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标；
S40，根据所述垂直距离和所述垂点坐标计算所述DBT系统在所述射线源当前发射射线时的贡献值，以同步完成DBT图像重建。
上述垂直距离为射线源发射产生相应投影数据的射线时，相对于探测器的垂直距离；上述垂点坐标为射线源发射产生相应投影数据的射线时，在探测器上的垂点坐标；依据上述垂直距离和垂点坐标计算得到的贡献值具有较高的准确性。依据上述贡献值进行相应的图像重建，保证了图像重建的效果。
本实施例提供的图像重建过程中的数据处理方法，可以在射线源发射射线时，获取钢珠在探测器上的钢珠投影坐标，以根据钢珠在校正模体上的钢珠空间坐标以及上述钢珠投影坐标确定DBT系统的投影矩阵，计算所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标，以根据上述垂直距离和垂点坐标计算所述DBT系统在所述射线源当前发射射线时的贡献值，使所确定的贡献值以射线源发射产生相应投影数据的射线时相对于探测器的垂直距离和垂点坐标为依据，具有较高的准确性，提高了依据上述贡献值进行相应图像重建的效果。
在一个实施例中，所述校正模体为两个底面分别固定设置多个钢珠的圆柱体，所述钢珠分别在圆柱体两个底面呈圆形排列；
所述在射线源发射射线时，获取钢珠在探测器上的钢珠投影坐标的过程之前，还包括：
将所述校正模体固定设置在所述射线源发射射线的一侧，所述校正模体的中心点与所述射线源的准直器中心点在所述探测器的一条垂线上。
本实施例将上述校正模体固定设置在射线源发射射线的一侧(若射线源向下发射射线，则将校正模体设置在射线源下方)，这样，射线源从可以照射到压迫对象的任意角度发射射线，提高了相应DBT系统扫描过程中的便利性；校正模体的中心点与所述射线源的准直器中心点在所述探测器的一条垂线上，即保证射线源的中心射线经过模体中心，上述校正模体各个底面钢珠排列形成的圆为同心圆，以简化依据各个钢珠的钢珠空间坐标以及相应的钢珠投影坐标计算相应投影矩阵过程中的计算量。
作为一个实施例，上述校正模体的尺寸可以依据压迫板的尺寸、射线源与压迫板之间的距离等因素设置，具体地，参考图4所示，上述圆柱体形的校正模体的高可以为2cm(厘米)，直径可以为4cm，校正模体两个底面的钢珠点成同心圆排列，上底面可以设置8个钢珠点，下底面可以设置7个钢珠点，共15个钢珠点。为了将校正模体上钢珠点与其投影点一一对应以便于提取几何参数，可以设计第一个钢珠点直径大小为1.5mm，剩余的钢珠点直径大小为1mm，这样就可以根据钢珠点投影的大小以及钢珠点排布的规律将模板上的钢珠点与其投影点一一对应起来。由于校正模体和乳腺等压迫对象同步成像，校正模体可以选用低密度的材料，如有机玻璃，这样可忽略模体的衰减，即认为每个角度下投影图像是乳腺和钢珠点的共同衰减。
若射线源向下发射射线，将校正模体固定设于射线源下表面和射线源一起转动扫描成像物；将设计好尺寸的模体固定于射线源下的下表面，尽量保证准直器中心和模体的中心保持在一条垂线上，即保证射线源的中心射线经过模体中心。一般来说，准直器的外形为立方体，在校正模体的固定过程中，只需利用工程技术克服模体本身的重力作用不使其在扫描过程中脱落移动即可，比如可选用强力胶固定相应的校正模体。在校正模体固定好之后，射线源与探测器之间的结构示意图可以参考图5所示，上述射线源可以在能照射到压迫对象的任意角度进行扫描。由于校正模体与扫描物不接触，故扫描物(压迫板以及相应的压迫对象等)的摆位可以和常规的DBT系统保持一致，不需要做更改。此时得到的多幅投影图像均为校正模体和乳腺共同衰减的图像，经负对数变换后所得到的图像可以参考图6(a)所示。
在一个实施例中，所述获取钢珠在探测器上的钢珠投影坐标的过程包括：
在探测器获得的投影数据中提取灰度值大于灰度阈值的钢珠投影点，将所述钢珠投影点的坐标确定为钢珠投影坐标。
上述灰度阈值的确定过程可以包括：
根据包括钢珠的投影数据对应的重建图像生成灰度直方图，将所述灰度直方图最后一个波峰对应的灰度值确定为灰度阈值。
由于钢珠点是高衰减物质，经负对数变换后其在投影图像上可以表现为高亮点，参考图6(b)所示，钢珠点在灰度值直方图上对应最后一个波峰附近的区域，如图6(b)所示。如图6(a)表示投影图像对应的重建图像，图6(b)为图像(a)对应的灰度直方图，横轴代表灰度值，纵轴代表灰度值对应的像素数量；可以通过搜索纵轴对应的一阶导数的左正右负的零点确定直方图的波峰。本实施例可以检测上述灰度直方图的最后一个波峰，从而确定灰度阈值，将投影图像中灰度值大于灰度阈值的点确定为钢珠投影点，进而分割出钢珠投影点，确定相应的钢珠投影坐标。
在实际射线照射过程中，射线穿过钢珠球心的路径是最长的，根据Lambert-Beer定律射线经过球心的路径衰减最多，所以钢珠点球心的投影应该是钢珠点投影的最暗区域。由于实际的钢珠点投影分布比较均匀，最暗区域为好几个像素相连，无法直接精准地求出钢珠点球心的坐标。而钢珠点投影的质心则可以精准地分割出来，考虑到钢珠点是球形，其在投影图像上大致呈圆形排列。即便有些钢珠点距离中心射线比较远，其投影为椭圆，由于钢珠点直径较小，而射线源到探测器的距离很大，利用质心确定的坐标和真实的钢珠点球心投影的坐标也相差不大，因而本实施例确定的钢珠投影坐标为钢珠点球心投影的坐标。
在一个实施例中，所述根据所述投影矩阵计算所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标的过程包括：
对所述投影矩阵进行RQ分解，获得所述DBT系统的内矩阵，根据所述内矩阵的矩阵元素分别确定所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标。
作为一个实施例，所述垂直距离为：SID＝K11du；
所述垂点坐标为：u0＝K13，v0＝K23；
其中，SID表示垂直距离，K11表示内矩阵K中第一行第一列元素，du表示探测器U轴的像素点尺寸大小，(u0,v0)表示垂点坐标，K13表示内矩阵K中第一行第三列元素，K23表示内矩阵K中第二行第三列元素。
本实施例可以对射线源发射相应射线时，射线源到探测器的垂直距离以及射线源在探测器上的垂点坐标进行准确计算，进一步保证了后续得到的射线源发射相应射线时对应的贡献值的准确性。
作为一个实施例，上述垂直距离和垂点坐标分别对应的确定方式的推导过程可以包括：
在DBT系统中，射线源以锥形束的形状向探测器发射射线，对于空间上的某个钢珠点(x,y,z)，经过射线的投射，可以在探测器上得到其对应的投影坐标，即钢珠点球心投影的坐标(u,v)，几何映射模型可以参考图3所示。由于已知校正模体的尺寸和钢珠点相对位置关系，若校正模体的几何中心，即圆柱体的几何中心为校正模体坐标系的坐标原点，通过DLT方法证明，只要校正模体上钢珠点的相对关系已知，无论校正模体坐标系如何建立，总可以求解出系统的几何参数信息。将投影点和标记点的成像几何关系在齐次坐标下写成透视投影矩阵的形式，利用投影矩阵P建立变换关系可以参考如下公式：
齐次坐标系是将N维向量用N+1维向量表示描述投影几何，可将旋转平移整合成一个矩阵进行计算。其中(u,v,1)表示图像坐标系下的投影点的齐次坐标，也即为经过像素离散化后的投影点的坐标。(x,y,z,1)表示模体坐标系下的空间点的坐标，P表示3×4的透视投影矩阵，w表示尺度因子。且P＝K[R|t]，其中K为系统的几何参数矩阵，也称为内矩阵。
其中，SID表示射线源到探测器的垂直距离，du表示像素点在探测器U轴的像素点尺寸大小，dv表示像素点在探测器V轴的像素点尺寸大小，(u0,v0)表示射线源在探测器上的垂点。对于DBT系统来说，K矩阵完整代表了用于重建的系统的几何参数。[R|t]矩阵是由R和t合成的矩阵，其中R表示投影物体坐标系到实际探测器坐标系的旋转变换，t表示投影物体坐标系到实际探测器坐标系的平移变换，[R|t]表示投影物体坐标系到实际探测器坐标系的旋转平移变换，也称为外矩阵。由于任意一个旋转操作，总可以分解为在XYZ轴三个方向上的旋转分量的乘积，用公式表示为：
R＝RxRyRz (3)
其中：
则：
t矩阵可表示为XYZ三个方向上的平移作用叠加，
对于几何校正来说，可以根据探测器的投影平面上钢珠点的钢珠投影坐标和钢珠点的钢珠空间坐标，则可求解出P矩阵，对于一个钢珠点(xi,yi,zi)来说，根据公式(1)有：
将其展开写成方程的形式，可以得到：
消除加权因子wi产生两个方程：
将上述两个方程格式化为矩阵形式，可得：
Mip＝0 (7)
Mi是个2×12的矩阵：
p是包含投影矩阵P的12个元素的列向量：
p＝[P11 P12 P13 P14 P21 P22 P23 P24 P31 P32 P33 P34]T (9)
由此可得，p向量有12个未知量，若需要求得p，则必须要n个钢珠点形成多个Mi矩阵，按行共同堆积形成M。其中，n的取值可以为28。相应地M是一个56×12的矩阵，且满足：
Mp＝0 (10)
对M进行奇异值分解(Singular Value Decomposition，SVD)可以有效地求解p：
M＝UDVT (11)
矩阵V的最后一列就是p的解，进而可得投影矩阵P的各个构成部分。P的前三行三列P3×3具有如下特征：P3×3＝K3×3R，即P3×3可分解为K矩阵的前三行三列乘以R矩阵，由于R矩阵是正交矩阵，K矩阵是上三角矩阵，可以通过对P矩阵RQ分解同时获得K和R矩阵。得到K矩阵之后，根据公式(2)则可得到：
由于尺度因子的存在，在进行RQ分解之前，可以对P矩阵进行尺度变换。对于旋转矩阵R来说，有由于K矩阵的第三行第三列值恒为1，则有这一约束可以消去尺度因子的影响，从而获得DBT系统真实对应的投影矩阵P。这样SID和(u0,v0)便可以得到求解。
参考图7，图7所示为一个实施例的图像重建过程中的数据处理装置结构示意图，包括：
第一获取模块10，用于在射线源发射射线时，获取钢珠在探测器上的钢珠投影坐标；其中，所述钢珠固定设置在校正模体上，所述校正模体固定设置在所述射线源发射射线的一侧；
第一确定模块20，用于根据钢珠在所述校正模体上的钢珠空间坐标以及所述钢珠投影坐标确定DBT系统的投影矩阵；其中，所述投影矩阵为表征所述射线通过所述钢珠照射至所述探测器时，所述钢珠空间坐标与所述钢珠投影坐标之间变换关系的矩阵；
第一计算模块30，用于根据所述投影矩阵计算所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标；
第二计算模块40，用于根据所述垂直距离和所述垂点坐标计算所述DBT系统在所述射线源当前发射射线时的贡献值。
在一个实施例中，所述校正模体为两个底面分别固定设置多个钢珠的圆柱体，所述钢珠分别在圆柱体两个底面呈圆形排列；
所述图像重建过程中的数据处理装置，还包括：
设置模块，用于将所述校正模体固定设置在所述射线源发射射线的一侧，所述校正模体的中心点与所述射线源的准直器中心点在所述探测器的一条垂线上。
在一个实施例中，所述第一获取模块进一步用于：
在探测器获得的投影数据中提取灰度值大于灰度阈值的钢珠投影点，将所述钢珠投影点的坐标确定为钢珠投影坐标。
在一个实施例中，所述第一计算模块进一步用于：
对所述投影矩阵进行RQ分解，获得所述DBT系统的内矩阵，根据所述内矩阵的矩阵元素分别确定所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标。
作为一个实施例，所述垂直距离为：SID＝K11du；
所述垂点坐标为：u0＝K13，v0＝K23；
其中，SID表示垂直距离，K11表示内矩阵K中第一行第一列元素，du表示探测器U轴的像素点尺寸大小，(u0,v0)表示垂点坐标，K13表示内矩阵K中第一行第三列元素，K23表示内矩阵K中第二行第三列元素。
关于图像重建过程中的数据处理装置的具体限定可以参见上文中对于图像重建过程中的数据处理方法的限定，在此不再赘述。上述图像重建过程中的数据处理装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
一种图像重建方法，包括如下步骤：
分别在射线源各次发射射线的过程中，根据上述任意一实施例所述的图像重建过程中的数据处理方法计算所述DBT系统在射线源该次发射射线时的贡献值，并获取所述射线源该次发射射线时，探测器探测得到的投影数据；
根据所述射线源各次发射射线分别产生的贡献值和投影数据进行图像重建。
上述射线源为DBT系统的射线源，射线源可以旋转，每隔设定角度发射一次射线，使探测器可以获得压迫板一面的校正模体以及压迫板另一面的压迫对象(如乳腺组织)受射线照射后投影至探测器上的投影数据。上述射线源的旋转角度范围和设定角度分别可以依据重建图像的所需精度确定，具体地，上述射线源的旋转角度范围可以设置为-12.5°至12.5°，上述设定角度可以设置为3°。
本实施例提供的图像重建方法以更为准确的贡献值为依据，有效提高的相应的重建效果。
在一个实施例中，上述根据所述射线源各次发射射线分别产生的贡献值和投影数据进行图像重建的过程包括：
将所述射线源各次发射射线分别产生的贡献值和投影数据代入重建迭代公式进行图像重建；所述重建迭代公式包括：
式中，表示第n次迭代后体素j的取值，J表示体素总数，的初始值取0，I表示探测器所包括的探测器单元个数，在每个探测器单元对应一条射线时，I也可以表示射线的数量。Aij表示压迫对象第j个体素对第i个探测器单元的射线的贡献值，上述贡献值可以通过计算X射线与体素间相交形成的体积或者X射线与体素的相交截取的长度等计算，具体可以根据射线源发射相应射线时相对于探测器的垂直距离和垂点坐标计算。yi表示第i个探测器单元的观测值，即投影数据中的第i个像素值，n代表迭代次数，λ表示迭代重建的松弛因子，用以调节迭代过程中偏差修正的程度，λ可影响重建的收敛速度和求解精度，具体的，λ取值越大，则收敛速度越快，但有可能引入噪声，λ取值越小，其收敛速度越慢，但迭代修正越精细，重建效果一般会越好，λ的取值范围可以为大于0且小于2，比如λ取0.8。
可选地，上述重建过程可以包括同一角度下所有射线对某个体素的修正值加权累加后才更新一次图像，因而可以减小同一投影下某一射线对体素进行修正时引入的误差，提高算法质量。经过多次迭代直至收敛，即可得到待重建的物体图像相应的重建操作可以在显存为12G的CUDA平台上进行并行计算。
作为一个实施例，在分别针对任意轨迹用常规成像方法、任意轨迹采用本发明进行重建成像方法、以及精准轨迹采用本发明提供的图像重建方法进行仿真的过程中，射线源依次向相应的校正模体和压迫对象发射9次射线，探测器可以分别获得9个投影数据；在射线源发射多次射线的过程中，探测器可以保持不动，相应扫描过程的扫描时间可以为7s(秒)；射线源到探测器的距离SID为66cm，旋转中心到探测器的距离(Axis-to-Image-Distance,AID)为4cm。探测器分辨率可以为2400×3000，像素大小为0.1mm。三者分别重建得到的图像可以参考图8所示，图8中，(a)表示任意轨迹用常规成像方法的重建结果，这里的任意轨迹是指在精准圆轨道上对射线源轨迹在上下左右前后6个方向加上±2cm范围内的随机扰动，常规成像即先扫描校正模体得出几何参数之后再扫描物体进行成像，然后利用之前校正模体计算出来的几何参数重建。(b)为任意轨迹采用本发明的方法重建结果，(c)为精准轨迹的重建结果。可以看到，针对常规的几何校正方法，即先扫描校正模体再进行扫描物的成像，然后利用之前校正模体计算出来的几何参数重建，这样重建出来的图像很模糊，与精准轨迹重建的图像相差较大，影响诊断结果。而本发明中的针对任意轨迹的重建结果和精准轨迹重建的图像相差无几，可以实现了射线源任意轨迹的成像。
上述图像重建方法，将校正模体固定于球管(射线源)下与物体(乳腺等压迫对象)一起成像，进行实时的几何校正，保证扫描物体时的相关几何参数与之前几何校正获得的几何参数保持相同，由此提高图像重建的质量。由于钢珠点的线性衰减系数要比人体任何组织都要大，在投影图像上可利用图像阈值分割出钢珠点，之后对分割后的钢珠点求解质心坐标作为钢珠点球心的投影坐标。然后利用钢珠点在模体上坐标和投影的坐标建立对应关系，可以运用DLT变换，可获得实时的几何参数信息，大大降低了对射线源旋转的要求，可以实现任意的轨迹旋转。同时，由于校正模体放置于射线源附近，则模体离成像物体距离较远，可以保证在DTS在有限角度内物体的投影不会有缺失的情况出现，即在有限角度内，物体的投影不会总被覆盖。
一种图像重建装置，包括：
第三计算模块，用于分别在射线源各次发射射线的过程中，根据如上任意一个实施例所述的图像重建过程中的数据处理装置计算所述DBT系统在射线源该次发射射线时的贡献值，并获取所述射线源该次发射射线时，探测器探测得到的投影数据；
重建模块，用于根据所述射线源各次发射射线分别产生的贡献值和投影数据进行图像重建。
关于图像重建装置的具体限定可以参见上文中对于图像重建方法的限定，在此不再赘述。上述图像重建装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是服务器，其内部结构图可以如图9所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的数据库用于存储相关投影数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种图像重建过程中的数据处理方法或者图像重建方法。
本领域技术人员可以理解，图9中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
基于如上所述的示例，在一个实施例中还提供一种计算机设备，该计算机设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其中，处理器执行所述程序时实现如上述各实施例中的任意一种图像重建过程中的数据处理方法或者图像重建方法。
上述计算机设备，通过所述处理器上运行的计算机程序，实现了图像重建效果的提升。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一非易失性的计算机可读取存储介质中，如本发明实施例中，该程序可存储于计算机系统的存储介质中，并被该计算机系统中的至少一个处理器执行，以实现包括如上述图像重建过程中的数据处理方法或者图像重建方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存储记忆体(Random Access Memory，RAM)等。
据此，在一个实施例中还提供一种计算机存储介质，其上存储有计算机程序，其中，该程序被处理器执行时实现如上述各实施例中的任意一种图像重建过程中的数据处理方法或者图像重建方法。
上述计算机存储介质，通过其存储的计算机程序，使所确定的贡献值以射线源发射产生相应投影数据的射线时相对于探测器的垂直距离和垂点坐标为依据，具有较高的准确性，提高了依据上述贡献值进行相应图像重建的效果。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种图像重建过程中的数据处理方法，其特征在于，包括如下步骤：
在射线源发射射线时，获取钢珠在探测器上的钢珠投影坐标；其中，所述钢珠固定设置在校正模体上，所述校正模体固定设置在所述射线源发射射线的一侧；
根据钢珠在所述校正模体上的钢珠空间坐标以及所述钢珠投影坐标确定DBT系统的投影矩阵；其中，所述投影矩阵为表征所述射线通过所述钢珠照射至所述探测器时，所述钢珠空间坐标与所述钢珠投影坐标之间变换关系的矩阵；
根据所述投影矩阵计算所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标；
根据所述垂直距离和所述垂点坐标计算所述DBT系统在所述射线源当前发射射线时的贡献值。
2.根据权利要求1所述的图像重建过程中的数据处理方法，其特征在于，所述校正模体为两个底面分别固定设置多个钢珠的圆柱体，所述钢珠分别在圆柱体两个底面呈圆形排列；
所述在射线源发射射线时，获取钢珠在探测器上的钢珠投影坐标的过程之前，还包括：
将所述校正模体固定设置在所述射线源发射射线的一侧，所述校正模体的中心点与所述射线源的准直器中心点在所述探测器的一条垂线上。
3.根据权利要求1所述的图像重建过程中的数据处理方法，其特征在于，所述获取钢珠在探测器上的钢珠投影坐标的过程包括：
在探测器获得的投影数据中提取灰度值大于灰度阈值的钢珠投影点，将所述钢珠投影点的坐标确定为钢珠投影坐标。
4.根据权利要求1至3任一项所述的图像重建过程中的数据处理方法，其特征在于，所述根据所述投影矩阵计算所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标的过程包括：
对所述投影矩阵进行RQ分解，获得所述DBT系统的内矩阵，根据所述内矩阵的矩阵元素分别确定所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标。
5.根据权利要求4所述的图像重建过程中的数据处理方法，其特征在于，所述垂直距离为：SID＝K11du；
所述垂点坐标为：u0＝K13，v0＝K23；
其中，SID表示垂直距离，K11表示内矩阵K中第一行第一列元素，du表示探测器U轴的像素点尺寸大小，(u0,v0)表示垂点坐标，K13表示内矩阵K中第一行第三列元素，K23表示内矩阵K中第二行第三列元素。
6.一种图像重建方法，其特征在于，包括如下步骤：
分别在射线源各次发射射线的过程中，根据权利要求1至5任一项所述的图像重建过程中的数据处理方法计算所述DBT系统在射线源该次发射射线时的贡献值，并获取所述射线源该次发射射线时，探测器探测得到的投影数据；
根据所述射线源各次发射射线分别产生的贡献值和投影数据进行图像重建。
7.一种图像重建过程中的数据处理装置，其特征在于，包括：
第一获取模块，用于在射线源发射射线时，获取钢珠在探测器上的钢珠投影坐标；其中，所述钢珠固定设置在校正模体上，所述校正模体固定设置在所述射线源发射射线的一侧；
第一确定模块，用于根据钢珠在所述校正模体上的钢珠空间坐标以及所述钢珠投影坐标确定DBT系统的投影矩阵；其中，所述投影矩阵为表征所述射线通过所述钢珠照射至所述探测器时，所述钢珠空间坐标与所述钢珠投影坐标之间变换关系的矩阵；
第一计算模块，用于根据所述投影矩阵计算所述射线源到探测器的垂直距离以及所述射线源在探测器上的垂点坐标；
第二计算模块，用于根据所述垂直距离和所述垂点坐标计算所述DBT系统在所述射线源当前发射射线时的贡献值。
8.一种图像重建装置，其特征在于，包括：
第三计算模块，用于分别在射线源各次发射射线的过程中，根据权利要求7所述的图像重建过程中的数据处理装置计算所述DBT系统在射线源该次发射射线时的贡献值，并获取所述射线源该次发射射线时，探测器探测得到的投影数据；
重建模块，用于根据所述射线源各次发射射线分别产生的贡献值和投影数据进行图像重建。
9.一种计算机设备，包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至5任意一项所述的图像重建过程中的数据处理方法或者权利要求6所述的图像重建方法。
10.一种计算机存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1至5任意一项所述的图像重建过程中的数据处理方法或者权利要求6所述的图像重建方法。
CN201811491822.9A 2018-12-07 2018-12-07 图像重建过程中的数据处理、图像重建方法和装置 Active CN109671128B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811491822.9A CN109671128B (zh) 2018-12-07 2018-12-07 图像重建过程中的数据处理、图像重建方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811491822.9A CN109671128B (zh) 2018-12-07 2018-12-07 图像重建过程中的数据处理、图像重建方法和装置

Publications (2)

Publication Number Publication Date
CN109671128A true CN109671128A (zh) 2019-04-23
CN109671128B CN109671128B (zh) 2023-06-30

Family

ID=66144662

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811491822.9A Active CN109671128B (zh) 2018-12-07 2018-12-07 图像重建过程中的数据处理、图像重建方法和装置

Country Status (1)

Country Link
CN (1) CN109671128B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110327064A (zh) * 2019-06-19 2019-10-15 东软医疗系统股份有限公司 源像距检测方法、装置及设备
CN111223159A (zh) * 2019-12-31 2020-06-02 广州华端科技有限公司 锥形束成像方法、装置、计算机设备和存储介质
TWI708217B (zh) * 2019-10-30 2020-10-21 行政院原子能委員會核能研究所 用於雙軸數位斷層合成造影系統的幾何校正方法及其系統
CN113395762A (zh) * 2021-04-18 2021-09-14 湖南财政经济学院 超宽带定位网络中位置校正方法及装置
CN113884519A (zh) * 2021-09-29 2022-01-04 清华大学 自导航x射线成像系统及成像方法
CN113963056A (zh) * 2021-09-07 2022-01-21 于留青 Ct图像重建方法、装置、电子设备以及存储介质
CN114324423A (zh) * 2021-12-24 2022-04-12 武汉联影生命科学仪器有限公司 焦点校正方法、装置、系统、计算机设备和可读存储介质
CN115844429A (zh) * 2022-12-05 2023-03-28 有方（合肥）医疗科技有限公司 成像数据校准方法、电子设备、存储介质及程序产品

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130243293A1 (en) * 2012-03-14 2013-09-19 Siemens Aktiengesellschaft Method for reducing motion artifacts
CN104665862A (zh) * 2015-02-16 2015-06-03 清华大学 在cbct中消除几何伪影的方法以及使用该方法的cbct系统
CN106056646A (zh) * 2016-05-24 2016-10-26 广州华端科技有限公司 Ct图像重建方法和系统
CN107233105A (zh) * 2017-05-24 2017-10-10 深圳先进技术研究院 一种用于ct图像重建的修正方法及修正系统
CN107545551A (zh) * 2017-09-07 2018-01-05 广州华端科技有限公司 数字乳腺体层合成图像的重建方法和系统
CN108122203A (zh) * 2016-11-29 2018-06-05 上海东软医疗科技有限公司 一种几何参数的校正方法、装置、设备及系统
CN108201447A (zh) * 2017-11-21 2018-06-26 深圳先进技术研究院 一种静态ct系统几何参数校正方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130243293A1 (en) * 2012-03-14 2013-09-19 Siemens Aktiengesellschaft Method for reducing motion artifacts
CN104665862A (zh) * 2015-02-16 2015-06-03 清华大学 在cbct中消除几何伪影的方法以及使用该方法的cbct系统
CN106056646A (zh) * 2016-05-24 2016-10-26 广州华端科技有限公司 Ct图像重建方法和系统
CN108122203A (zh) * 2016-11-29 2018-06-05 上海东软医疗科技有限公司 一种几何参数的校正方法、装置、设备及系统
CN107233105A (zh) * 2017-05-24 2017-10-10 深圳先进技术研究院 一种用于ct图像重建的修正方法及修正系统
CN107545551A (zh) * 2017-09-07 2018-01-05 广州华端科技有限公司 数字乳腺体层合成图像的重建方法和系统
CN108201447A (zh) * 2017-11-21 2018-06-26 深圳先进技术研究院 一种静态ct系统几何参数校正方法

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110327064A (zh) * 2019-06-19 2019-10-15 东软医疗系统股份有限公司 源像距检测方法、装置及设备
TWI708217B (zh) * 2019-10-30 2020-10-21 行政院原子能委員會核能研究所 用於雙軸數位斷層合成造影系統的幾何校正方法及其系統
CN111223159A (zh) * 2019-12-31 2020-06-02 广州华端科技有限公司 锥形束成像方法、装置、计算机设备和存储介质
CN113395762A (zh) * 2021-04-18 2021-09-14 湖南财政经济学院 超宽带定位网络中位置校正方法及装置
CN113963056A (zh) * 2021-09-07 2022-01-21 于留青 Ct图像重建方法、装置、电子设备以及存储介质
CN113963056B (zh) * 2021-09-07 2022-08-26 于留青 Ct图像重建方法、装置、电子设备以及存储介质
CN113884519A (zh) * 2021-09-29 2022-01-04 清华大学 自导航x射线成像系统及成像方法
CN113884519B (zh) * 2021-09-29 2022-07-12 清华大学 自导航x射线成像系统及成像方法
CN114324423A (zh) * 2021-12-24 2022-04-12 武汉联影生命科学仪器有限公司 焦点校正方法、装置、系统、计算机设备和可读存储介质
CN115844429A (zh) * 2022-12-05 2023-03-28 有方（合肥）医疗科技有限公司 成像数据校准方法、电子设备、存储介质及程序产品
CN115844429B (zh) * 2022-12-05 2024-06-21 有方（合肥）医疗科技有限公司 成像数据校准方法、电子设备、存储介质及程序产品

Also Published As

Publication number Publication date
CN109671128B (zh) 2023-06-30

Similar Documents

Publication Publication Date Title
CN109671128A (zh) 图像重建过程中的数据处理、图像重建方法和装置
Sonke et al. Respiratory correlated cone beam CT
US8615118B2 (en) Techniques for tomographic image by background subtraction
US7711087B2 (en) Patient setup using tomosynthesis techniques
US20110176723A1 (en) Motion Correction in Cone-Beam CT by Tracking Internal and External Markers Using Cone-Beam Projection From a kV On-Board Imager: Four-Dimensional Cone-Beam CT and Tumor Tracking Implications
CN106994023B (zh) 锥形束计算机断层成像系统的几何参数确定方法
CN103329168B (zh) 针对spect/ct系统的迭代锥形射束ct重建的截断补偿
US8705817B2 (en) Measurement of geometric quantities intrinsic to an anatomical system
CN110381839A (zh) 处置计划图像、分次内3D图像和分次内2D x射线图像的图像配准
US20090202127A1 (en) Method And System For Error Compensation
CN107133549B (zh) Ect运动门控信号获取方法及ect图像重建方法
US9795356B2 (en) System and method of X-ray dose distribution for computed tomography based on simulation
Qi et al. Extraction of tumor motion trajectories using PICCS‐4DCBCT: a validation study
WO2010032167A2 (en) Mr segmentation using transmission data in hybrid nuclear/mr imaging
CN108013888A (zh) 用于医学成像设备的散射辐射补偿
CN102016923A (zh) 图像伪影的减少
Lu et al. A comprehensive study on decreasing the kilovoltage cone‐beam CT dose by reducing the projection number
CN106667512A (zh) X射线成像设备的几何校正方法、乳腺断层成像设备
CN110310307A (zh) 在知道计算机断层扫描（ct）重构区域的情况下的可变形图像配准
CN107865658A (zh) 用于修正合成电子密度图的方法和设备
Shi et al. Library based x‐ray scatter correction for dedicated cone beam breast CT
Gagné et al. The impact of tumor motion upon CT image integrity and target delineation: Impact of tumor motion upon CT image integrity
CN109658465A (zh) 图像重建过程中的数据处理、图像重建方法和装置
CN112001979B (zh) 运动伪影处理方法、系统、可读存储介质和设备
CN107192726A (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
TR01 Transfer of patent right

Effective date of registration: 20230809

Address after: No. 501-3, 5th Floor, Building 9, No. 129 Dingxing Road, Tangjiawan Town, High tech Zone, Zhuhai City, Guangdong Province, 519000

Patentee after: Zhuhai Ruitong Medical Equipment Co.,Ltd.

Address before: 10th Floor, Qixi Building, No. 286, Science Avenue, Huangpu District, Guangzhou, Guangdong, 510700

Patentee before: GUANGZHOU HUADUAN TECHNOLOGY Co.,Ltd.

TR01 Transfer of patent right