CN115797485A - 图像去伪影方法、系统、电子设备及存储介质 - Google Patents

图像去伪影方法、系统、电子设备及存储介质 Download PDF

Info

Publication number
CN115797485A
CN115797485A CN202211319983.6A CN202211319983A CN115797485A CN 115797485 A CN115797485 A CN 115797485A CN 202211319983 A CN202211319983 A CN 202211319983A CN 115797485 A CN115797485 A CN 115797485A
Authority
CN
China
Prior art keywords
detector
data
ratio
slit region
image
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.)
Pending
Application number
CN202211319983.6A
Other languages
English (en)
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.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare 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 Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN202211319983.6A priority Critical patent/CN115797485A/zh
Publication of CN115797485A publication Critical patent/CN115797485A/zh
Priority to US18/383,910 priority patent/US20240138798A1/en
Pending legal-status Critical Current

Links

Images

Classifications

    • 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/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
    • 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/037—Emission tomography
    • 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/48—Diagnostic techniques
    • A61B6/483—Diagnostic techniques involving scattered radiation
    • 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/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • 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
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Toxicology (AREA)
  • Pulmonology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明公开了一种图像去伪影方法、系统、电子设备及存储介质。其中，图像去伪影方法包括：获取对目标对象进行计算机断层扫描的过程中，探测器的原始扫描数据；采用散射参考信号对包含散射信号的原始扫描数据进行伪影校正；其中，散射参考信号根据进行计算机断层扫描的过程中，无狭缝区域的探测器的能量数据，和有狭缝区域的探测器的能量数据得到。本发明采用根据无狭缝区域的探测器的能量数据，和有狭缝区域的探测器的能量数据得到的散射参考信号，对包含散射信号的原始扫描数据进行伪影校正，进而去除医学图像的环形伪影，达到更好的显像效果。

Description

图像去伪影方法、系统、电子设备及存储介质
技术领域
本发明涉及计算机断层图像技术领域，尤其涉及一种图像去伪影方法、系统、电子设备及存储介质。
背景技术
现在的计算机断层图像设备是通过透射线(线束)构建图像的，透射线通过X射线产生，但是X射线通过人体或物体后，会产生透射线和散射线，散射线对重构图像是有害的，需要消除掉。探测器能将照射在其上的X射线转换成电信号，是装置中最重要的部分。探测器将射线管发出的X射线束导向探测器模块的相应像素点上。探测器以模块的形式组装在一起，因此，探测器模块之间不可避免的存在一些缝隙，如图1所示。有些X射线可能经过缝隙进入探测器最边缘的像素点(即图1中的深色方块)时发生散射现象，这样会导致探测器边缘像素点接收到比实际更多的X射线的接收信号。在最终重建的图像上，探测器边缘的像素点就很容易产生伪影，影响最终重建的图像。现有技术中，为了去除伪影通常会在硬件设备上进行改进，例如可以在图1的缝隙上方增加栅格，减少X射线的散射信号。但是这类方法一方面增加了硬件的复杂度，提高了产品的制作成本，另一方面也无法根本性地解决计算机断层图像设备扫描到的图像会产生伪影的问题。
发明内容
本发明要解决的技术问题是为了克服现有技术中因探测器模块之间存在的缝隙，X射线在探测器边缘产生散射时，探测器会接收到比实际更多的X射线，导致最终重建的图像上产生伪影缺陷，提供一种图像去伪影方法、系统、电子设备及存储介质。
本发明是通过下述技术方案来解决上述技术问题：
第一方面，本发明提供一种图像去伪影方法，所述图像去伪影方法包括：
获取对目标对象进行计算机断层扫描的过程中，探测器的原始扫描数据；
采用散射参考信号对包含散射信号的所述原始扫描数据进行伪影校正；其中，所述散射参考信号根据进行计算机断层扫描的过程中，无狭缝区域的所述探测器的能量数据，和有狭缝区域的所述探测器的能量数据得到。
优选地，所述无狭缝区域的所述探测器的能量数据包括无模体无狭缝区域的所述探测器的能量强度值，和有模体无狭缝区域的所述探测器的能量强度值；
所述有狭缝区域的所述探测器的能量数据包括无模体有狭缝区域的所述探测器的能量强度值，和有模体有狭缝区域的所述探测器的能量强度值；
所述模体为进行计算机断层扫描的过程中模拟所述目标对象的模体。
优选地，所述采用散射参考信号对包含散射信号的所述原始扫描数据进行伪影校正的步骤包括：
对所述原始扫描数据进行散射强度评估得到原始散射信号；
将所述原始散射信号输入多项式拟合结果得到散射信号；其中，所述多项式拟合结果的系数根据所述散射参考信号确定；
将所述原始扫描数据去除所述散射信号，得到伪影校正后的原始扫描数据。
优选地，所述多项式拟合结果的系数为第一比值与第二比值的差值；
其中，第一比值为所述有模体无狭缝区域的所述探测器的能量强度值与无模体无狭缝区域的所述探测器的能量强度值的比值；
第二比值为所述有模体有狭缝区域的所述探测器的能量强度值与无模体有狭缝区域的所述探测器的能量强度值；
或，
所述多项式结果的系数为第一比值与第二比值的差值的平滑差值；所述平滑差值为对第一比值与第二比值的差值进行平滑处理后得到的差值。
优选地，所述图像去伪影方法还包括：
采用投影参考信号对包含投影信号的所述原始扫描数据进行伪影校正；
其中，所述投影参考信号根据进行计算机断层扫描的过程中，无模体有狭缝区域的所述探测器的能量强度值和有模体有狭缝区域的所述探测器的能量强度值得到。
优选地，所述采用投影参考信号对包含投影信号的所述原始扫描数据进行伪影校正的步骤包括：
将所述原始扫描数据转换成投影数据；
将所述投影数据输入多项式拟合结果得到投影校正数据；其中多项式拟合结果的系数根据投影参考信号确定。
优选地，所述多项式拟合结果的系数为第三比值的对数值，其中，所述第三比值为所述无模体有狭缝区域的所述探测器的能量强度值和所述有模体有狭缝区域的所述探测器的能量强度值的比值；
或，
所述多项式拟合结果的系数为第三比值的对数值和平滑对数值；所述平滑对数值为对第三比值的对数值进行平滑处理的对数值。
第二方面，本发明提供一种图像去伪影系统，所述图像去伪影系统包括：
获取模块，用于获取对目标对象进行计算机断层扫描的过程中，探测器的原始扫描数据；
校正模块，用于采用散射参考信号对包含散射信号的所述原始扫描数据进行伪影校正；其中，所述散射参考信号根据进行计算机断层扫描的过程中，无狭缝区域的所述探测器的能量数据，和有狭缝区域的所述探测器的能量数据得到。
第三方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述的图像去伪影方法。
第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述的图像去伪影方法。
本发明的积极进步效果在于：
本发明采用根据无狭缝区域的探测器的能量数据，和有狭缝区域的探测器的能量数据得到的散射参考信号，对包含散射信号的原始扫描数据进行伪影校正，在实践中，可以通过伪影校正后的原始扫描数据重建图像，进而得到去除伪影的医学影像，达到更好的显像效果。
附图说明
图1为X射线通过探测器模块之间的缝隙产生伪影的原理示意图；
图2为本发明的实施例1提供的一种图像去伪影方法的第一流程图；
图3为本发明的实施例1提供的一种图像去伪影方法的第二流程图；
图4为本发明的实施例1提供的一种图像去伪影方法的第三流程图；
图5为本发明的实施例1提供的一种图像去伪影方法的原始散射信号强度示意图；
图6为本发明的实施例1提供的一种图像去伪影方法的拟合后的散射信号强度示意图；
图7为本发明的实施例1提供的一种图像去伪影方法的原始散射信号强度和拟合后的散射信号强度的对比图；
图8为本发明的实施例2提供的一种图像去伪影系统的第一结构图；
图9为本发明的实施例2提供的一种图像去伪影系统的第二结构图；
图10为本发明的实施例3提供的一种电子设备的结构示意图。
具体实施方式
下面通过实施例的方式进一步说明本发明，但并不因此将本发明限制在所述的实施例范围之中。
实施例1
本实施例提供一种图像去伪影方法，其可应用于CT设备(Computed Tomography，即电子计算机断层扫描)和PET设备(Positron Emission Computed Tomography，正电子发射型计算机断层显像)等计算机断层图像设备构建的图像。参见图2，图像去伪影方法包括：
S1、获取对目标对象进行计算机断层扫描的过程中，探测器的原始扫描数据。其中，目标对象为生物组织，例如可以通过CT设备或PET设备对人的头部进行计算机断层扫描，目标对象即为人的头部。
S2、采用散射参考信号对包含散射信号的原始扫描数据进行伪影校正。
其中，散射参考信号根据进行计算机断层扫描的过程中，无狭缝区域的探测器的能量数据，和有狭缝区域的探测器的能量数据得到。其中，能量数据包括散射信号、投影信号等能量数据。
当获取狭缝区域的探测器的能量数据时，探测器中除了狭缝区域之外的其他区域被遮挡以使狭缝区域的原始扫描数据包含的散射信号的含量小于含量阈值。将狭缝区域设置于一块挡板上，狭缝区域由一条有固定预设宽度、狭窄细长的缝孔组成，将有狭缝区域的挡板放置在探测器上，以接收原始扫描数据的固定含量的散射信号。含量阈值根据狭缝宽度和原始扫描数据包含的散射信号总含量确定。例如：狭缝区域只允许照射20个像素点(该像素点为探测器像素点，用于接受X射线，输出响应数据)的散射信号。在实际操作过程中，需要移动狭缝区域，以使所有探测器像素点都能被X射线照射到。
在一个可选的实施方式中，无狭缝区域的探测器的能量数据包括无模体无狭缝区域的探测器的能量强度值，和有模体无狭缝区域的探测器的能量强度值。
有狭缝区域的探测器的能量数据包括无模体有狭缝区域的探测器的能量强度值，和有模体有狭缝区域的探测器的能量强度值。
其中，模体为进行计算机断层扫描的过程中模拟目标对象的模体。其中，模体可采用特氟龙、水等材质的模体，模体厚度大于或者等于3毫米，以模仿人体组织的厚度。
根据伪影校正后的原始扫描数据进行图像重建，进而就能得到目标对象的去伪影医学图像。
在一个可选的实施方式中，可采用滤波反投影方法(FBP)对伪影校正后的原始扫描数据进行图像重建，以获得目标对象的去伪影医学图像。
在本实施例中，采用根据无狭缝区域的探测器的能量数据，和有狭缝区域的探测器的能量数据得到的散射参考信号，对包含散射信号的原始扫描数据进行伪影校正，进而去除医学图像的环形伪影，达到更好的显像效果。
在一个可选的实施方式中，参见图3，步骤S2包括：
S21、对原始扫描数据进行散射强度评估得到原始散射信号。
散射强度评估是通过算法、模型、神经网络等方式对原始扫描数据信号进行分析，以得到原始扫描数据中的散射信号即原始散射信号。散射强度评估方式包括卷积、进行蒙特卡罗算法和深度学习等方式。
S22、将原始散射信号输入多项式拟合结果得到散射信号。
其中，多项式拟合是用一个多项式展开去拟合整个原始散射信号，以得到客观的散射信号。经过多项式拟合后的散射信号更贴近原始扫描数据中实际包含的散射信号。
其中，多项式拟合结果的系数根据散射参考信号确定。
在一个可选的实施方式中，多项式拟合结果的系数为第一比值与第二比值的差值。其中，第一比值为有模体无狭缝区域的探测器的能量强度值与无模体无狭缝区域的探测器的能量强度值的比值；第二比值为有模体有狭缝区域的探测器的能量强度值与无模体有狭缝区域的探测器的能量强度值。
其中，模体可采用特氟龙、水等材质的模体，模体的厚度大于或者等于3毫米，以模仿人体组织的厚度。
在一个可选的实施方式中，多项式结果的系数为第一比值与第二比值的差值的平滑差值；平滑差值为对第一比值与第二比值的差值进行平滑处理后得到的差值。
平滑处理方式包括一阶指数平滑处理方式、二阶指数平滑处理方式和三阶指数平滑处理方式等。
S23、将原始扫描数据去除散射信号，得到伪影校正后的原始扫描数据。
在本实施方式中，通过散射强度评估和多项式拟合算法，最终达到去除原始扫描数据的散射信号的效果，实现对原始扫描数据进行伪影校正，简化了伪影校正的方式，降低了维运成本，达到了更好的显像效果。
在一个可选的实施方式中，图像去伪影方法还包括：
S3、采用投影参考信号对包含投影信号的原始扫描数据进行伪影校正。
其中，投影参考信号根据进行计算机断层扫描的过程中，无模体有狭缝区域的探测器的能量强度值和有模体有狭缝区域的探测器的能量强度值得到。其中，模体为进行计算机断层扫描的过程中模拟目标对象的模体。其中，模体可采用特氟龙、水等材质的模体，模体的厚度大于或者等于3毫米，以模仿人体组织的厚度。
在一个可选的实施方式中，参见图4，步骤S3包括：
S31、将原始扫描数据转换成投影数据。
通常利用log运算(对数运算)，将从探测器采集到的原始扫描数据转换成投影数据。投影数据即探测器采集的原始数据排列成以探测器通道为横轴、以扫描视野为纵轴的二维矩阵，其本质上是医学成像上各点所形成的曲线缠绕叠加。
S32、将投影数据输入多项式拟合结果得到投影校正数据。
其中，多项式拟合结果的系数根据投影参考信号确定。
在一个可选的实施方式中，多项式拟合结果的系数为第三比值的对数值，其中，第三比值为无模体有狭缝区域的探测器的能量强度值和有模体有狭缝区域的探测器的能量强度值的比值。其中，模体可采用特氟龙、水等材质的模体，模体的厚度大于或者等于3毫米，以模仿人体组织的厚度。
在一个可选的实施方式中，多项式拟合结果的系数为第三比值的对数值和平滑对数值；平滑对数值为对第三比值的对数值进行平滑处理的对数值。
平滑处理方式包括一阶指数平滑处理方式、二阶指数平滑处理方式和三阶指数平滑处理方式等。
在本实施方式中，通过投影参考信号，进一步对原始扫描数据的投影信号进行伪影校正，从而更优化了的医学影像的显像效果。
下面介绍一种图像去伪影方法的示例的具体步骤：
S41、对原始扫描数据进行空气校正。
空气校正即在扫描设备的扫描范围内不放置任何物体进行一系列扫描得到的一组数据，然后从原始扫描数据中减去只扫描空气所得的参考值，以得到更精确的原始扫描数据。
S42、对空气校正后的原始扫描数据执行log运算(对数运算)，以得到投影数据。
S43、对投影数据输入多项式拟合结果A得到投影校正数据。
其中，多项式拟合结果A的系数根据投影参考信号确定。在本示例中，多项式拟合结果A的系数为对数值的平滑对数值，平滑对数值为对第三比值的对数值进行平滑处理的对数值。其中，第三比值为无模体有狭缝区域的探测器的能量强度值和有模体有狭缝区域的探测器的能量强度值的比值。其中，模体可采用特氟龙、水等材质的模体，模体厚度大于或者等于3毫米，以模仿人体组织的厚度。
S44、对投影校正数据进行反log运算(反对数运算)，得到投影伪影校正后的原始扫描数据，对投影伪影校正后的原始扫描数据进行散射强度评估，得到原始散射信号。散射评估方式包括卷积、进行蒙特卡罗算法和深度学习等方式。图5为进行散射强度评估后，得到原始散射信号曲线。图5的横坐标表征探测器通道的编号；图5的纵坐标表征散射信号的强度值。
S45、将原始散射信号输入多项式拟合结果B得到散射信号。
图6中的曲线部分为原始散射信号(同图5)，曲线突出部分为拟合得到的散射信号，图6的横坐标表征探测器通道的编号；图6的纵坐标表征散射信号的强度值。图6相对于图5更贴近实际获取到的散射信号的强度值，曲线突出部分产生就是因为有些X射线可能经过缝隙进入探测器最边缘的像素点(即图1中的深色方块)时发生散射现象，这样会导致探测器边缘像素点接收到比实际更多的X射线的接收信号，所以会在原本平滑的曲线上产生多个突出尖峰。这也是在最终重建的图像上，探测器边缘的像素点产生伪影的原因。
图7为原始散射信号和多项式拟合后得到的散射信号的对比图。图7的横坐标表征探测器通道的编号；图7的纵坐标表征散射信号的强度值。
其中，在本示例中，多项式拟合结果B的系数为第一比值与第二比值的差值的平滑差值。其中，第一比值为有模体无狭缝区域的探测器的能量强度值与无模体无狭缝区域的探测器的能量强度值的比值；第二比值为有模体有狭缝区域的探测器的能量强度值与无模体有狭缝区域的探测器的能量强度值。平滑差值为对第一比值与第二比值的差值进行平滑处理后得到的差值。其中，模体可采用特氟龙、水等材质的模体，模体的厚度大于或者等于3毫米，以模仿人体组织的厚度。
S46、将上述原始扫描数据去除散射信号，进一步进行伪影校正。
S47、采用滤波反投影方法(FBP)对经过上述伪影校正操作后的原始扫描数据进行图像重建，以获得目标对象的去伪影医学图像。
实施例2
本实施例提供一种图像去伪影系统，参见图8，图像去伪影系统包括：
获取模块1，用于获取对目标对象进行计算机断层扫描的过程中，探测器的原始扫描数据。
校正模块2，用于采用散射参考信号对包含散射信号的原始扫描数据进行伪影校正；其中，散射参考信号根据进行计算机断层扫描的过程中，无狭缝区域的探测器的能量数据，和有狭缝区域的探测器的能量数据得到。
在一个可选的实施方式中，无狭缝区域的探测器的能量数据包括无模体无狭缝区域的探测器的能量强度值，和有模体无狭缝区域的探测器的能量强度值。
有狭缝区域的探测器的能量数据包括无模体有狭缝区域的探测器的能量强度值，和有模体有狭缝区域的探测器的能量强度值。
模体为进行计算机断层扫描的过程中模拟目标对象的模体。其中，模体可采用特氟龙、水等材质的模体，模体的厚度大于或者等于3毫米，以模仿人体组织的厚度。
在一个可选的实施方式中，参见图9，图像去伪影系统还包括：
重建模块3，用于采用滤波反投影方法(FBP)对伪影校正后的原始扫描数据进行图像重建，以获得目标对象的去伪影医学图像。
在本实施例中，图像去伪影系统采用根据无狭缝区域的探测器的能量数据，和有狭缝区域的探测器的能量数据得到的散射参考信号，对包含散射信号的原始扫描数据进行伪影校正，进而去除医学图像的环形伪影，达到更好的显像效果。
在一个可选的实施方式中，参见图9，图像去伪影系统还包括：
拟合模块4，用于对原始扫描数据进行散射强度评估得到原始散射信号；还用于将原始散射信号输入多项式拟合结果得到散射信号。
其中，多项式拟合结果的系数根据散射参考信号确定。
在一个可选的实施方式中，多项式拟合结果的系数为第一比值与第二比值的差值。其中，第一比值为有模体无狭缝区域的探测器的能量强度值与无模体无狭缝区域的探测器的能量强度值的比值；第二比值为有模体有狭缝区域的探测器的能量强度值与无模体有狭缝区域的探测器的能量强度值。其中，模体可采用特氟龙、水等材质的模体，模体的厚度大于或者等于3毫米，以模仿人体组织的厚度。
在一个可选的实施方式中，多项式结果的系数为第一比值与第二比值的差值的平滑差值；平滑差值为对第一比值与第二比值的差值进行平滑处理后得到的差值。
平滑处理方式包括一阶指数平滑处理方式、二阶指数平滑处理方式和三阶指数平滑处理方式等。
校正模块2，还用于将原始扫描数据去除散射信号，得到伪影校正后的原始扫描数据。
在本实施方式中，图像去伪影系统通过散射强度评估和多项式拟合算法，最终达到去除原始扫描数据的散射信号的效果，实现对原始扫描数据进行伪影校正，简化了伪影校正的方式，降低了维运成本，达到了更好的显像效果。
在一个可选的实施方式中，校正模块2，还用于采用投影参考信号对包含投影信号的原始扫描数据进行伪影校正。
其中，投影参考信号根据进行计算机断层扫描的过程中，无模体有狭缝区域的探测器的能量强度值和有模体有狭缝区域的探测器的能量强度值得到。其中，模体为进行计算机断层扫描的过程中模拟目标对象的模体。模体可采用特氟龙、水等材质的模体，模体的厚度大于或者等于3毫米，以模仿人体组织的厚度。
在一个可选的实施方式中，拟合模块4，还用于将原始扫描数据转换成投影数据；还用于将投影数据输入多项式拟合结果得到投影校正数据。
其中，多项式拟合结果的系数根据投影参考信号确定。
在一个可选的实施方式中，多项式拟合结果的系数为第三比值的对数值，其中，第三比值为无模体有狭缝区域的探测器的能量强度值和有模体有狭缝区域的探测器的能量强度值的比值。其中，模体可采用特氟龙、水等材质的模体，模体的厚度大于或者等于3毫米，以模仿人体组织的厚度。
在一个可选的实施方式中，多项式拟合结果的系数为第三比值的对数值和平滑对数值；平滑对数值为对第三比值的对数值进行平滑处理的对数值。
平滑处理方式包括一阶指数平滑处理方式、二阶指数平滑处理方式和三阶指数平滑处理方式等。
在本实施方式中，图像去伪影系统通过投影参考信号，进一步对原始扫描数据的投影信号进行伪影校正，从而更优化了的医学影像的显像效果。
实施例3
本实施例提供了一种电子设备，图10为该电子设备的模块示意图。电子设备包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行程序时实现实施例1的图像去伪影方法。图10显示的电子设备30仅仅是一个示例，不应对本发明实施例的功能和使用范围带来任何限制。
如图10所示，电子设备30可以以通用计算设备的形式表现，例如其可以为服务器设备。电子设备30的组件可以包括但不限于：上述至少一个处理器31、上述至少一个存储器32、连接不同系统组件(包括存储器32和处理器31)的总线33。
总线33包括数据总线、地址总线和控制总线。
存储器32可以包括易失性存储器，例如随机存取存储器(RAM)321和/或高速缓存存储器322，还可以进一步包括只读存储器(ROM)323。
存储器32还可以包括具有一组(至少一个)程序模块324的程序/实用工具325，这样的程序模块324包括但不限于：操作系统、一个或者多个应用程序、其它程序模块以及程序数据，这些示例中的每一个或某种组合中可能包括网络环境的实现。
处理器31通过运行存储在存储器32中的计算机程序，从而执行各种功能应用以及数据处理，例如本发明实施例1的图像去伪影方法。
电子设备30也可以与一个或多个外部设备34(例如键盘、指向设备等)通信。这种通信可以通过输入/输出(I/O)接口35进行。并且，模型生成的设备30还可以通过网络适配器36与一个或者多个网络(例如局域网(LAN)，广域网(WAN)和/或公共网络，例如因特网)通信。如图10所示，网络适配器36通过总线33与模型生成的设备30的其它模块通信。应当明白，尽管图中未示出，可以结合模型生成的设备30使用其它硬件和/或软件模块，包括但不限于：微代码、设备驱动器、冗余处理器、外部磁盘驱动阵列、RAID(磁盘阵列)系统、磁带驱动器以及数据备份存储系统等。
应当注意，尽管在上文详细描述中提及了电子设备的若干单元/模块或子单元/模块，但是这种划分仅仅是示例性的并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多单元/模块的特征和功能可以在一个单元/模块中具体化。反之，上文描述的一个单元/模块的特征和功能可以进一步划分为由多个单元/模块来具体化。
实施例4
本实施例提供了一种计算机可读存储介质，其上存储有计算机程序，程序被处理器执行时实现实施例1的图像去伪影方法。
其中，可读存储介质可以采用的更具体可以包括但不限于：便携式盘、硬盘、随机存取存储器、只读存储器、可擦拭可编程只读存储器、光存储器件、磁存储器件或上述的任意合适的组合。
在可能的实施方式中，本发明还可以实现为一种程序产品的形式，其包括程序代码，当程序产品在终端设备上运行时，程序代码用于使终端设备执行实现实施例1的图像去伪影方法。
其中，可以以一种或多种程序设计语言的任意组合来编写用于执行本发明的程序代码，程序代码可以完全地在用户设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行、部分在用户设备上部分在远程设备上执行或完全在远程设备上执行。
虽然以上描述了本发明的具体实施方式，但是本领域的技术人员应当理解，这仅是举例说明，本发明的保护范围是由所附权利要求书限定的。本领域的技术人员在不背离本发明的原理和实质的前提下，可以对这些实施方式做出多种变更或修改，但这些变更和修改均落入本发明的保护范围。

Claims (10)

1.一种图像去伪影方法，其特征在于，所述图像去伪影方法包括：
获取对目标对象进行计算机断层扫描的过程中，探测器的原始扫描数据；
采用散射参考信号对包含散射信号的所述原始扫描数据进行伪影校正；其中，所述散射参考信号根据进行计算机断层扫描的过程中，无狭缝区域的所述探测器的能量数据，和有狭缝区域的所述探测器的能量数据得到。
2.如权利要求1所述的图像去伪影方法，其特征在于，所述无狭缝区域的所述探测器的能量数据包括无模体无狭缝区域的所述探测器的能量强度值，和有模体无狭缝区域的所述探测器的能量强度值；
所述有狭缝区域的所述探测器的能量数据包括无模体有狭缝区域的所述探测器的能量强度值，和有模体有狭缝区域的所述探测器的能量强度值；
所述模体为进行计算机断层扫描的过程中模拟所述目标对象的模体。
3.如权利要求2所述的图像去伪影方法，其特征在于，所述采用散射参考信号对包含散射信号的所述原始扫描数据进行伪影校正的步骤包括：
对所述原始扫描数据进行散射强度评估得到原始散射信号；
将所述原始散射信号输入多项式拟合结果得到散射信号；其中，所述多项式拟合结果的系数根据所述散射参考信号确定；
将所述原始扫描数据去除所述散射信号，得到伪影校正后的原始扫描数据。
4.如权利要求3所述的图像去伪影方法，其特征在于，所述多项式拟合结果的系数为第一比值与第二比值的差值；
其中，第一比值为所述有模体无狭缝区域的所述探测器的能量强度值与无模体无狭缝区域的所述探测器的能量强度值的比值；
第二比值为所述有模体有狭缝区域的所述探测器的能量强度值与无模体有狭缝区域的所述探测器的能量强度值；
或，
所述多项式结果的系数为第一比值与第二比值的差值的平滑差值；所述平滑差值为对第一比值与第二比值的差值进行平滑处理后得到的差值。
5.如权利要求1-4任一项所述的图像去伪影方法，其特征在于，所述图像去伪影方法还包括：
采用投影参考信号对包含投影信号的所述原始扫描数据进行伪影校正；
其中，所述投影参考信号根据进行计算机断层扫描的过程中，无模体有狭缝区域的所述探测器的能量强度值和有模体有狭缝区域的所述探测器的能量强度值得到。
6.如权利要求5所述的图像去伪影方法，其特征在于，所述采用投影参考信号对包含投影信号的所述原始扫描数据进行伪影校正的步骤包括：
将所述原始扫描数据转换成投影数据；
将所述投影数据输入多项式拟合结果得到投影校正数据；其中多项式拟合结果的系数根据投影参考信号确定。
7.如权利要求6所述的图像去伪影方法，其特征在于，所述多项式拟合结果的系数为第三比值的对数值，其中，所述第三比值为所述无模体有狭缝区域的所述探测器的能量强度值和所述有模体有狭缝区域的所述探测器的能量强度值的比值；
或，
所述多项式拟合结果的系数为第三比值的对数值的平滑对数值；所述平滑对数值为对第三比值的对数值进行平滑处理的对数值。
8.一种图像去伪影系统，其特征在于，所述图像去伪影系统包括：
获取模块，用于获取对目标对象进行计算机断层扫描的过程中，探测器的原始扫描数据；
校正模块，用于采用散射参考信号对包含散射信号的所述原始扫描数据进行伪影校正；其中，所述散射参考信号根据进行计算机断层扫描的过程中，无狭缝区域的所述探测器的能量数据，和有狭缝区域的所述探测器的能量数据得到。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1-7中任一项所述的图像去伪影方法。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1-7中任一项所述的图像去伪影方法。
CN202211319983.6A 2022-10-26 2022-10-26 图像去伪影方法、系统、电子设备及存储介质 Pending CN115797485A (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202211319983.6A CN115797485A (zh) 2022-10-26 2022-10-26 图像去伪影方法、系统、电子设备及存储介质
US18/383,910 US20240138798A1 (en) 2022-10-26 2023-10-26 Method for removing artifact in image, electronic device, and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211319983.6A CN115797485A (zh) 2022-10-26 2022-10-26 图像去伪影方法、系统、电子设备及存储介质

Publications (1)

Publication Number Publication Date
CN115797485A true CN115797485A (zh) 2023-03-14

Family

ID=85433928

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211319983.6A Pending CN115797485A (zh) 2022-10-26 2022-10-26 图像去伪影方法、系统、电子设备及存储介质

Country Status (2)

Country Link
US (1) US20240138798A1 (zh)
CN (1) CN115797485A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116563413A (zh) * 2023-07-04 2023-08-08 中国科学院深圳先进技术研究院 校正参数确定方法、图像重建方法、装置、设备及介质

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118593003B (zh) * 2024-07-31 2024-11-01 上海西门子医疗器械有限公司 医用x射线成像的散射校正方法及装置、存储介质

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116563413A (zh) * 2023-07-04 2023-08-08 中国科学院深圳先进技术研究院 校正参数确定方法、图像重建方法、装置、设备及介质
CN116563413B (zh) * 2023-07-04 2024-04-09 中国科学院深圳先进技术研究院 校正参数确定方法、图像重建方法、装置、设备及介质

Also Published As

Publication number Publication date
US20240138798A1 (en) 2024-05-02

Similar Documents

Publication Publication Date Title
US9265475B2 (en) Methods and apparatus for scatter correction for CBCT system and cone-beam image reconstruction
CN112348936B (zh) 一种基于深度学习的低剂量锥束ct图像重建方法
US9036771B2 (en) System and method for denoising medical images adaptive to local noise
CN115797485A (zh) 图像去伪影方法、系统、电子设备及存储介质
Niu et al. Shading correction for on‐board cone‐beam CT in radiation therapy using planning MDCT images
US7187794B2 (en) Noise treatment of low-dose computed tomography projections and images
US8532350B2 (en) Dose reduction and image enhancement in tomography through the utilization of the object's surroundings as dynamic constraints
US7747057B2 (en) Methods and apparatus for BIS correction
KR20190060054A (ko) 치과용 콘-빔 ｃｔ의 산란선 보정방법 및 보정장치
US9589373B2 (en) Monte carlo modeling of field angle-dependent spectra for radiographic imaging systems
EP4235581A1 (en) Medical image processing method, medical image processing apparatus, program, method for producing a trained machine learning model, apparatus for producing a trained machine learning model and computer-readable storage medium
JPWO2015020072A1 (ja) Ｘ線ｃｔ装置および補正処理装置
US8938104B2 (en) Systems and methods for adaptive filtering
Fan et al. Image‐domain shading correction for cone‐beam CT without prior patient information
CN107106114B (zh) 运算装置、x射线ct装置及图像重构方法
CN113096210A (zh) 图像重建方法、装置、电子设备和存储介质
EP3404618B1 (en) Poly-energetic reconstruction method for metal artifacts reduction
Wang et al. Locally linear transform based three‐dimensional gradient‐norm minimization for spectral CT reconstruction
US20230394630A1 (en) Method and system for training and tuning neural network models for denoising
US20230274473A1 (en) Artificial intelligence based 3d reconstruction
Elhamiasl et al. Low-dose CT simulation from an available higher dose CT scan
Li et al. An empirical data inconsistency metric (DIM) driven CT image reconstruction method
US20240312076A1 (en) Predicting scattered signal of x-ray, and correcting scattered beam
CN111091516B (zh) 一种基于人工智能的抗散射光栅方法及装置
Elhamiasl et al. The Effect of CT Dose Reduction on Proton Therapy Dose Calculation

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