CN110796605B - 图像处理方法和装置、存储介质、处理器和终端 - Google Patents

图像处理方法和装置、存储介质、处理器和终端 Download PDF

Info

Publication number
CN110796605B
CN110796605B CN201810864673.XA CN201810864673A CN110796605B CN 110796605 B CN110796605 B CN 110796605B CN 201810864673 A CN201810864673 A CN 201810864673A CN 110796605 B CN110796605 B CN 110796605B
Authority
CN
China
Prior art keywords
data
image
artifact
channels
average value
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
CN201810864673.XA
Other languages
English (en)
Other versions
CN110796605A (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.)
Siemens Shanghai Medical Equipment Ltd
Original Assignee
Siemens Shanghai Medical Equipment 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 Siemens Shanghai Medical Equipment Ltd filed Critical Siemens Shanghai Medical Equipment Ltd
Priority to CN201810864673.XA priority Critical patent/CN110796605B/zh
Publication of CN110796605A publication Critical patent/CN110796605A/zh
Application granted granted Critical
Publication of CN110796605B publication Critical patent/CN110796605B/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
    • G06T5/00—Image enhancement or restoration
    • G06T5/73—Deblurring; Sharpening
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/20—Image enhancement or restoration using local operators
    • 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)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本申请涉及图像处理方法和装置、存储介质、处理器和终端。图像处理方法包括：获取定位像的数据；确定定位像的全部通道中的边缘通道是否被物体遮挡以生成确定结果；基于确定结果生成表示定位像中的伪影的伪影数据；以及从定位像的数据中去除伪影数据以生成去除伪影的图像数据。本申请的技术方案容易实施，且计算复杂度低，能够集成到当前的定位像重建技术中，实现了优化CT扫描图像，提升医疗用户体验的技术效果。

Description

图像处理方法和装置、存储介质、处理器和终端
技术领域
本申请涉及医疗成像领域。具体地，本申请涉及图像处理方法和装置、存储介质、处理器和终端。
背景技术
采用滚珠轴承的X-射线管组件(XTA)具有沿着Z方向的显著的热膨胀效应，为了补偿该热膨胀效应需要应用准直器实时控制技术。准直器的设计需要考虑焦点热运动(focusthermal movement)，CT系统(computer tomography system，计算机断层扫描系统)能够检测该运动并相应移动准直器。为了减少CT扫描中的定位像辐射剂量(topo dose)，准直器在设计时并未考虑焦点热运动。当发生焦点热运动时，检测该运动造成的误差，调节准直器进行移动的方法称为准直器控制(z-control)。但是在许多情况下，在开始进行定位像的扫描时，准直器没有处于需要的位置，这会造成重建的定位像中产生轻微横纹的伪影，给诊断过程带来不利影响。
虽然准直器控制能够在CT扫描的一些情况下产生有效作用，但对于定位像来说，在很多情况下并不是很有效。
发明内容
本申请实施例提供了去除定位像中伪影的方法、装置、存储介质和处理器，以至少解决现有技术中不能有效消除定位像中的伪影的问题。
根据本申请实施例的一个方面，提供了图像处理方法，包括：获取定位像的数据；确定定位像的全部通道中的边缘通道是否被物体遮挡以生成确定结果；基于确定结果生成表示定位像中的伪影的伪影数据；以及从定位像的数据中去除伪影数据以生成去除伪影的图像数据。
以这样的方式，以图像处理的方式处理定位像，计算出伪影并将伪影从定位像中消除，提升了定位像的准确性，改善了用户体验。
根据本申请的示例性实施例，基于确定结果生成表示定位像中的伪影的伪影数据包括：如果确定结果为边缘通道未被物体遮挡，则计算边缘通道的数据的平均值，并根据边缘通道的数据的平均值生成伪影数据。
以这样的方式，在边缘通道未被遮挡的情况下生成伪影数据。
根据本申请的示例性实施例，通过下式生成伪影数据： 其中，Sj表示伪影数据，T表示边缘通道的通道数量，Ai，j表示定位像的数据，i＝1：N，j＝1：M，N表示全部通道的通道数量，M表示重建后的总图像长度，max(Ai，j)表示定位像的通道的数据中的最大值，min(Ai，j)表示定位像的通道的数据中的最小值。
以这样的方式，提供生成伪影数据的具体方式，以根据边缘通道的数量和定位像的数据生成伪影数据。
根据本申请的示例性实施例，基于确定结果生成表示定位像中的伪影的伪影数据包括：如果确定结果为边缘通道被物体遮挡，则沿通道排列方向计算定位像的数据的平均值；对定位像的数据的平均值进行平滑处理以生成平滑值；以及根据平滑值与定位像的数据的平均值的差值生成伪影数据。
以这样的方式，在边缘通道被物体遮挡的情况下生成伪影数据。
根据本申请的示例性实施例，通过下式计算定位像的数据的平均值：其中，Pj表示定位像的数据的平均值，Ai，j表示定位像的数据，i＝1：N，j＝1：M，N表示全部通道的通道数量，M表示重建后的总图像长度。
以这样的方式，提供计算定位像的数据的平均值的具体方式，以根据定位像的数据的全部通道计算定位像数据的平均值。
根据本申请的示例性实施例，对定位像的数据的平均值进行平滑处理包括：对定位像的数据的平均值应用平滑滤波器，通过下式生成平滑值：其中，/>表示平滑值，h表示平滑滤波器。
以这样的方式，生成平滑的定位像的数据。
根据本申请的示例性实施例，通过下式生成伪影数据：其中，Sj表示伪影数据。
以这样的方式，根据平滑的定位像的数据和定位像的数据的平均值计算出表示伪影的数据。
根据本申请实施例的另一方面，还提供了图像处理装置，包括：接收单元，用于获取定位像的数据；遮挡确定单元，用于确定定位像的全部通道中的边缘通道是否被物体遮挡以生成确定结果；伪影确定单元，用于基于确定结果生成表示定位像中的伪影的伪影数据；以及图像生成单元，用于从定位像的数据中去除伪影数据以生成去除伪影的图像数据。
以这样的方式，以图像处理的方式处理定位像，计算出伪影并将伪影从定位像中消除，提升了定位像的准确性，改善了用户体验。
根据本申请的示例性实施例，基于确定结果生成表示定位像中的伪影的伪影数据包括：如果确定结果为边缘通道未被物体遮挡，则通过伪影确定单元计算边缘通道的数据的平均值，并根据边缘通道的数据的平均值生成伪影数据。
以这样的方式，在边缘通道未被遮挡的情况下生成伪影数据。
根据本申请的示例性实施例，通过下式生成伪影数据： 其中，Sj表示伪影数据，T表示边缘通道的通道数量，Ai，j表示定位像的数据，i＝1：N，j＝1：M，N表示全部通道的通道数量，M表示重建后的总图像长度，max(Ai，j)表示定位像的通道的数据中的最大值，min(Ai，j)表示定位像的通道的数据中的最小值。
以这样的方式，提供生成伪影数据的具体方式，以根据边缘通道的数量和定位像的数据生成伪影数据。
根据本申请的示例性实施例，基于确定结果生成表示定位像中的伪影的伪影数据包括：如果确定结果为边缘通道被物体遮挡，则通过伪影确定单元沿通道排列方向计算定位像的数据的平均值；对定位像的数据的平均值进行平滑处理以生成平滑值；以及根据平滑值与定位像的数据的平均值的差值生成伪影数据。
以这样的方式，在边缘通道被物体遮挡的情况下生成伪影数据。
根据本申请的示例性实施例，通过下式计算定位像的数据的平均值：其中，Pj表示定位像的数据的平均值，Ai，j表示定位像的数据，i＝1：N，j＝1：M，N表示全部通道的通道数量，M表示重建后的总图像长度。
以这样的方式，提供计算定位像的数据的平均值的具体方式，以根据定位像的数据的全部通道计算定位像数据的平均值。
根据本申请的示例性实施例，对定位像的数据的平均值进行平滑处理包括：对定位像的数据的平均值应用平滑滤波器，通过下式生成平滑值：其中，/>表示平滑值，h表示平滑滤波器。
以这样的方式，生成平滑的定位像的数据。
根据本申请的示例性实施例，通过下式生成伪影数据：其中，Sj表示伪影数据。
以这样的方式，根据平滑的定位像的数据和定位像的数据的平均值计算出表示伪影的数据。
根据本申请实施例的另一方面，还提供了存储介质，存储介质包括存储的计算机程序，其中，计算机程序运行时控制存储介质所在的设备执行上述图像处理方法。
根据本申请实施例的另一方面，还提供了处理器，处理器运行计算机程序，其中，计算机程序运行时执行上述图像处理方法。
根据本申请实施例的另一方面，还提供了终端，包括：一个或多个处理器，存储器，以及一个或多个计算机程序，其中，一个或多个计算机程序被存储在存储器中，并且被配置为由一个或多个处理器执行，一个或多个计算机程序执行上述图像处理方法。
根据本申请实施例的另一方面，还提供了计算机程序产品，计算机程序产品被有形地存储在计算机可读介质上并且包括计算机可执行指令，计算机可执行指令在被执行时使至少一个处理器执行上述图像处理方法。
根据本申请实施例的图像处理技术能够实现在存储介质、处理器、终端和计算机程序产品中。以这样的方式，以图像处理的方式处理定位像，计算出伪影并将伪影从定位像中消除，提升了定位像的准确性，改善了用户体验。
在本申请实施例中，提供了以图像处理的方式计算定位像中的伪影并从定位像中减去伪影的技术方案，以至少解决定位像中存在伪影的技术问题。本申请的技术方案容易实施，且计算复杂度低，能够集成到当前的定位像重建技术中，实现了优化CT扫描图像，提升医疗用户体验的技术效果。
附图说明
此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：
图1是根据本申请实施例的图像处理方法的流程图；
图2是根据本申请实施例的在边缘通道未被遮挡的情况下边缘通道的数据的平均值的示意图；；
图3是根据本申请实施例的在边缘通道未被遮挡的情况下存在伪影的定位像的示图；
图4是根据本申请实施例的在边缘通道未被遮挡的情况下去除伪影的定位像的示图；
图5是根据本申请实施例的在边缘通道被遮挡的情况下沿着通道排列方向的定位像的数据的平均值的示意图；
图6是根据本申请实施例的在边缘通道被遮挡的情况下沿着通道排列方向的定位像的数据的平均值被平滑处理的示意图；
图7是根据本申请实施例的在边缘通道被遮挡的情况下存在伪影的定位像的示图；
图8是根据本申请实施例的在边缘通道被遮挡的情况下去除伪影的定位像的示图；
图9是根据本申请实施例的图像处理装置的示意图。
附图标号说明：
S102，获取定位像的数据；
S104，确定定位像的全部通道中的边缘通道是否被物体遮挡以生成确定结果；
S106，基于确定结果生成表示定位像中的伪影的伪影数据；
S108，从定位像的数据中去除伪影数据以生成去除伪影的图像数据；
1，图像处理装置；
101，接收单元；
103，遮挡确定单元；
105，伪影确定单元；
107，图像生成单元。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分的实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
需要说明的是，本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本申请的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或模块或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或模块或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或模块或单元。
根据本申请实施例，提供了图像处理方法。图1是根据本申请实施例的图像处理方法的流程图。如图1所示，根据本申请实施例的图像处理方法包括：步骤S102，获取定位像的数据；步骤104，确定定位像的全部通道中的边缘通道是否被物体遮挡以生成确定结果；步骤S106，基于确定结果生成表示定位像中的伪影的伪影数据；以及步骤S108，从定位像的数据中去除伪影数据以生成去除伪影的图像数据。
根据本申请实施例，在计算机断层扫描中，扫描得到定位像，对定位像的数据进行图像处理，以去除定位像中由于准直器误差造成的伪影，例如准直器对X射线的遮挡造成的横纹。具体地，在根据本申请实施例的图像处理方法中，首先获取定位像的数据。定位像的数据包括多个通道的数据。X射线管组件排列在第一方向，对应于沿第一方向排列的多个通道，因此，通道的排列方向与第一方向相同。例如，在扫描时，X射线管组件沿着与第一方向垂直的第二方向进行扫描，则得到宽为沿第一方向的X射线管组件对应的通道宽度，长为沿第二方向的扫描长度的定位像(或CT扫描图像)。其中，每个通道的数据是根据X射线的扫描结果得到CT扫描对应的CT值。全部多个通道的数据中，针对边缘通道被遮挡或者未被遮挡的两种情况进行不同的下一步的处理。在CT扫描过程中，通常将被扫描的物体(例如人体)置于定位像的画面中央，而在扫描的通道排列方向的两侧不被被扫描物体遮挡，这些在此情况下不被遮挡的部分对应于边缘通道，边缘通道在扫描过程中获得对应于基本为零的CT值的数据，因为扫描的X射线没有被物体吸收。在一些特殊情况下，CT扫描也可能将被扫描物体布置为遮挡这些边缘通道(例如病人身体较宽，或者对特定部位的扫描)，在扫描过程中，这些边缘通道会获得对应于扫描到的物体的CT值的数据。根据本申请实施例，针对边缘通道是否被遮挡提供两种图像处理方法来确定准直器造成的伪影。因此，根据本申请实施例，判断全部通道中的边缘通道是否被物体遮挡。针对边缘通道被遮挡或者未被遮挡的确定结果根据定位像的数据进行相应运算，以图像处理的方式获得表示定位相中的伪影的伪影数据。在获得伪影数据后，以图像处理的方式从定位像的数据中去除伪影数据，得到去除伪影的图像的图像数据。以图像处理的方式处理定位像，计算出伪影并将伪影从定位像中消除，提升了定位像的准确性，改善了用户体验。
根据本申请的示例性实施例，基于确定结果生成表示定位像中的伪影的伪影数据包括：如果确定结果为边缘通道未被物体遮挡，则计算边缘通道的数据的平均值，并根据边缘通道的数据的平均值生成伪影数据。具体地，针对边缘通道被遮挡或者未被遮挡的两种情况中确定边缘通道未被遮挡的情况进行如下处理。针对边缘通道的数据进行平均值运算，对于边缘通道未被遮挡的情况，如果不存在伪影，那么边缘通道中的数据应是X射线仅穿过空气而产生的接近于零值的CT值，若定位像中存在伪影，则由于伪影是准直器对X射线的遮挡产生的，那么伪影在定位相中会产生突变的CT值数据，边缘通道的数据因此不是接近于零值的CT值而是遮挡对应的CT值。对边缘通道的数据进行平均运算后，获得平均的边缘数据，这一数据可以作为定位像中的伪影数据，以在边缘通道未被遮挡的情况下生成伪影数据。
图2是根据本申请实施例的在边缘通道未被遮挡的情况下边缘通道的数据的平均值的示意图。如图2所示，横轴为扫描行数，纵轴为对应于CT值的HU值。在存在伪影的定位像的数据中，边缘通道具有负值的HU值，这是伪影造成的。在对这些数据进行平均后，得到代表平均的伪影数据的一条曲线。
根据本申请的示例性实施例，通过下式生成伪影数据： 其中，Sj表示伪影数据，T表示边缘通道的通道数量，Ai，j表示定位像的数据，i＝1：N，j＝1：M，N表示全部通道的通道数量，M表示重建后的总图像长度，max(Ai，j)表示定位像的通道的数据中的最大值，min(Ai，j)表示定位像的通道的数据中的最小值。根据本申请的示例性实施例，提供生成伪影数据的具体方式。具体地，从定位像的图像的边缘通道中，去除最大值和最小值。可以理解，为了去除误差，将定位像的数据中最大值和最小值去除，然后再进行平均值计算。例如，上式针对图像为N×M大小的定位像，Ai，j作为定位像的数据是第i个通道对应的第j个扫描行的CT值。对于T数量的边缘通道，从这些边缘通道数据的总和中去除被视为误差值的最大值和最小值，然后对剩余的边缘通道的数据求平均值，得到平均的边缘通道的数据，这个数据可以作为伪影数据。对于每一个扫描行，得到与这一行对应的伪影数据。以这样的方式在边缘通道未被遮挡的情况下根据边缘通道的数据以图像处理的方式计算得出伪影数据。从定位像的数据中去除伪影数据，具体地，从每一扫描行的定位像数据中去除对应的一行的伪影数据，生成的图像中被去除了伪影，提高了图像质量。
图3是根据本申请实施例的在边缘通道未被遮挡的情况下存在伪影的定位像的示图。如图3所示，能够观察到图中沿纵向的扫描方向每隔一些扫描行或像素存在横纹的伪像，这些伪像在医生或其他用户进行医疗的过程中造成不利影响。
图4是根据本申请实施例的在边缘通道未被遮挡的情况下去除伪影的定位像的示图。如图4所示，在从定位像的数据中如上去除伪影数据后，得到保持了有用的图像数据但消除了伪影的定位像。
根据本申请的示例性实施例，基于确定结果生成表示定位像中的伪影的伪影数据包括：如果确定结果为边缘通道被物体遮挡，则沿通道排列方向计算定位像的数据的平均值；对定位像的数据的平均值进行平滑处理以生成平滑值；以及根据平滑值与定位像的数据的平均值的差值生成伪影数据。在一些扫描情形中，若定位像中，边缘通道被物体遮挡，则以如下的另外的方式计算伪影数据。在CT扫描中，每一个扫描切片对应于一个扫描行(通道排列方向)通道在扫描时获得的CT值。若这一个扫描切片是正常的扫描切片，则其中的数据对应于X射线穿过被扫描物体(例如人体)而得的对应CT值。由于被扫描物体的连续性，每一个扫描切片与相邻扫描切片的数据的平均值是相同或者接近的。而伪影对应于准直器造成的遮挡，存在伪影的扫描切片的数据相对正常扫描数据是突变值，因此，若某一扫描切片对应的数据的平均值明显不同于前一扫描切片或后一扫描切片，说明图像中的对应位置存在伪影。对平均后的定位像的数据进行平滑处理，可得到对应于连续的被扫描物体的CT值的估计值，根据该估计值和定位像数据的平均值求两者的差值，即从定位像数据中去除对应于被扫描物体的CT值，得到的差值即对应于伪影的伪影数据。
根据本申请的示例性实施例，通过下式计算定位像的数据的平均值：其中，Pj表示定位像的数据的平均值，Ai，j表示定位像的数据，i＝1：N，j＝1：M，N表示全部通道的通道数量，M表示重建后的总图像长度。例如，上式针对图像为N×M大小的定位像，Ai，j作为定位像的数据是第i个通道对应的第j个像素的CT值。图5是根据本申请实施例的在边缘通道被遮挡的情况下沿着通道排列方向的定位像的数据的平均值的示意图。如图5所示，横轴为扫描行数，纵轴为对应于CT值的HU值。例如，对与第1个扫描行，对N个通道的数据的CT值求平均数，接着对下一个扫描行，对N个通道的数据的CT值求平均数，直至第M个扫描行，对N个通道的数据的CT值求平均数。求得每一个扫描行的CT值的平均数，其中，平均数的突变说明了在对应的扫描行存在伪影数据。
以这样的方式，提供计算定位像的数据的平均值的具体方式，以根据定位像的数据的全部通道计算定位像数据的平均值。
根据本申请的示例性实施例，对定位像的数据的平均值进行平滑处理包括：对定位像的数据的平均值应用平滑滤波器，通过下式生成平滑值：其中，/>表示平滑值，h表示平滑滤波器。图6是根据本申请实施例的在边缘通道被遮挡的情况下沿着通道排列方向的定位像的数据的平均值被平滑处理的示意图。图6所示，横轴为扫描行数，纵轴为对应于CT值的HU值。平滑滤波器h可以根据实际情形进行设置。例如，若检测发现每隔多个扫描行出现伪影，则根据间隔的扫描行的数量设置采用的平滑滤波器h。平滑处理后的数据代表了被扫描物体的CT值的估计值。根据平滑处理后的数据和每一行定位像数据的平均值的差值，可以得出代表伪影的伪影数据。根据本申请的示例性实施例，通过下式生成伪影数据：/>其中，Sj表示伪影数据。以这样的方式，根据平滑的定位像的数据和定位像的数据的平均值计算出表示伪影的数据。
图7是根据本申请实施例的在边缘通道被遮挡的情况下存在伪影的定位像的示图。如图7所示，能够观察到图中沿纵向的扫描方向每隔一些扫描行或像素存在横纹的伪像，这些伪像在医生或其他用户进行医疗的过程中造成不利影响。
图8是根据本申请实施例的在边缘通道被遮挡的情况下去除伪影的定位像的示图。如图8所示，在从定位像的数据中如上去除伪影数据后，得到保持了有用的图像数据但消除了伪影的定位像。
根据本申请实施例，还提供了图像处理装置。图9是根据本申请实施例的图像处理装置的示意图。如图9所示，图像处理装置1包括：接收单元101，用于获取定位像的数据；遮挡确定单元103，用于确定定位像的全部通道中的边缘通道是否被物体遮挡以生成确定结果；伪影确定单元105，用于基于确定结果生成表示定位像中的伪影的伪影数据；以及图像生成单元107，用于从定位像的数据中去除伪影数据以生成去除伪影的图像数据。根据本申请实施例的图像处理装置执行如上所述的图像处理方法，在此不再赘述。以这样的方式，以图像处理的方式处理定位像，计算出伪影并将伪影从定位像中消除，提升了定位像的准确性，改善了用户体验。
根据本申请的示例性实施例，基于确定结果生成表示定位像中的伪影的伪影数据包括：如果确定结果为边缘通道未被物体遮挡，则通过伪影确定单元计算边缘通道的数据的平均值，并根据边缘通道的数据的平均值生成伪影数据。以这样的方式，在边缘通道未被遮挡的情况下生成伪影数据。
根据本申请的示例性实施例，通过下式生成伪影数据： 其中，Sj表示伪影数据，T表示边缘通道的通道数量，Ai，j表示定位像的数据，i＝1：N，j＝1：M，N表示全部通道的通道数量，M表示重建后的总图像长度，max(Ai，j)表示定位像的通道的数据中的最大值，min(Ai，j)表示定位像的通道的数据中的最小值。以这样的方式，提供生成伪影数据的具体方式，以根据边缘通道的数量和定位像的数据生成伪影数据。
根据本申请的示例性实施例，基于确定结果生成表示定位像中的伪影的伪影数据包括：如果确定结果为边缘通道被物体遮挡，则通过伪影确定单元沿通道排列方向计算定位像的数据的平均值；对定位像的数据的平均值进行平滑处理以生成平滑值；以及根据平滑值与定位像的数据的平均值的差值生成伪影数据。以这样的方式，在边缘通道被物体遮挡的情况下生成伪影数据。
根据本申请的示例性实施例，通过下式计算定位像的数据的平均值：其中，Pj表示定位像的数据的平均值，Ai，j表示定位像的数据，i＝1：N，j＝1：M，N表示全部通道的通道数量，M表示重建后的总图像长度。以这样的方式，提供计算定位像的数据的平均值的具体方式，以根据定位像的数据的全部通道计算定位像数据的平均值。
根据本申请的示例性实施例，对定位像的数据的平均值进行平滑处理包括：对定位像的数据的平均值应用平滑滤波器，通过下式生成平滑值：其中，/>表示平滑值，h表示平滑滤波器。以这样的方式，生成平滑的定位像的数据。
根据本申请的示例性实施例，通过下式生成伪影数据：其中，Sj表示伪影数据。以这样的方式，根据平滑的定位像的数据和定位像的数据的平均值计算出表示伪影的数据。
如上的图像处理过程由根据本申请实施例的图像处理装置执行，其具体过程与上述的图像处理方法相同，在此不再赘述。
根据本申请实施例的另一方面，还提供了存储介质，存储介质包括存储的计算机程序，其中，计算机程序运行时控制存储介质所在的设备执行上述图像处理方法。
根据本申请实施例的另一方面，还提供了处理器，处理器运行计算机程序，其中，计算机程序运行时执行上述图像处理方法。
根据本申请实施例的另一方面，还提供了终端，包括：一个或多个处理器，存储器，以及一个或多个计算机程序，其中，一个或多个计算机程序被存储在存储器中，并且被配置为由一个或多个处理器执行，一个或多个计算机程序执行上述图像处理方法。
根据本申请实施例的另一方面，还提供了计算机程序产品，计算机程序产品被有形地存储在计算机可读介质上并且包括计算机可执行指令，计算机可执行指令在被执行时使至少一个处理器执行上述图像处理方法。
根据本申请实施例的图像处理技术能够实现在存储介质、处理器、终端和计算机程序产品中。以这样的方式，以图像处理的方式处理定位像，计算出伪影并将伪影从定位像中消除，提升了定位像的准确性，改善了用户体验。
在本申请的上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
在本申请所提供的几个实施例中，应该理解到，所揭露的技术内容，可通过其它的方式实现。其中，以上所描述的装置实施例仅仅是示意性的，例如所述单元或模块的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或模块或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，模块或单元的间接耦合或通信连接，可以是电性或其它的形式。
所述作为分离部件说明的单元或模块可以是或者也可以不是物理上分开的，作为单元或模块显示的部件可以是或者也可以不是物理单元或模块，即可以位于一个地方，或者也可以分布到多个网络单元或模块上。可以根据实际的需要选择其中的部分或者全部单元或模块来实现本实施例方案的目的。
另外，在本申请各个实施例中的各功能单元或模块可以集成在一个处理单元或模块中，也可以是各个单元或模块单独物理存在，也可以两个或两个以上单元或模块集成在一个单元或模块中。上述集成的单元或模块既可以采用硬件的形式实现，也可以采用软件功能单元或模块的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅是本申请的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本申请原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本申请的保护范围。

Claims (10)

1.图像处理方法，其特征在于，包括：
获取定位像的数据；
确定所述定位像的全部通道中的边缘通道是否被物体遮挡以生成确定结果；
基于所述确定结果生成表示所述定位像中的伪影的伪影数据；以及
从所述定位像的数据中去除所述伪影数据以生成去除所述伪影的图像数据；
基于所述确定结果生成表示所述定位像中的伪影的伪影数据包括：
如果所述确定结果为所述边缘通道未被物体遮挡，则计算所述边缘通道的数据的平均值，并根据所述边缘通道的数据的平均值生成所述伪影数据，通过下式生成所述伪影数据：
其中，Sj表示所述伪影数据，T表示所述边缘通道的通道数量，Ai,j表示所述定位像的数据，i＝1:N,j＝1:M，N表示所述全部通道的通道数量，M表示重建后的总图像长度，max(Ai，j)表示所述定位像的通道的数据中的最大值，min(Ai，j)表示所述定位像的通道的数据中的最小值；
基于所述确定结果生成表示所述定位像中的伪影的伪影数据包括：
如果所述确定结果为所述边缘通道被物体遮挡，则沿通道排列方向计算所述定位像的数据的平均值；对所述定位像的数据的平均值进行平滑处理以生成平滑值；以及根据所述平滑值与所述定位像的数据的平均值的差值生成所述伪影数据，通过下式计算所述定位像的数据的平均值：
其中，Pj表示所述定位像的数据的平均值，Ai,j表示所述定位像的数据，i＝1:N,j＝1:M，N表示所述全部通道的通道数量，M表示重建后的总图像长度。
2.根据权利要求1所述的图像处理方法，其特征在于，对所述定位像的数据的平均值进行平滑处理包括：
对所述定位像的数据的平均值应用平滑滤波器，通过下式生成所述平滑值：
其中，表示所述平滑值，h表示所述平滑滤波器。
3.根据权利要求2所述的图像处理方法，其特征在于，通过下式生成所述伪影数据：
其中，Sj表示所述伪影数据。
4.图像处理装置，其特征在于，包括：
接收单元，用于获取定位像的数据；
遮挡确定单元，用于确定所述定位像的全部通道中的边缘通道是否被物体遮挡以生成确定结果；
伪影确定单元，用于基于所述确定结果生成表示所述定位像中的伪影的伪影数据；以及
图像生成单元，用于从所述定位像的数据中去除所述伪影数据以生成去除所述伪影的图像数据；
基于所述确定结果生成表示所述定位像中的伪影的伪影数据包括：
如果所述确定结果为所述边缘通道未被物体遮挡，则通过所述伪影确定单元计算所述边缘通道的数据的平均值，并根据所述边缘通道的数据的平均值生成所述伪影数据，通过下式生成所述伪影数据：
其中，Sj表示所述伪影数据，T表示所述边缘通道的通道数量，Ai,j表示所述定位像的数据，i＝1:N,j＝1:M，N表示所述全部通道的通道数量，M表示重建后的总图像长度，max(Ai，j)表示所述定位像的通道的数据中的最大值，min(Ai，j)表示所述定位像的通道的数据中的最小值；
基于所述确定结果生成表示所述定位像中的伪影的伪影数据包括：
如果所述确定结果为所述边缘通道被物体遮挡，则通过所述伪影确定单元沿通道排列方向计算所述定位像的数据的平均值；对所述定位像的数据的平均值进行平滑处理以生成平滑值；以及根据所述平滑值与所述定位像的数据的平均值的差值生成所述伪影数据，通过下式计算所述定位像的数据的平均值：
其中，Pj表示所述定位像的数据的平均值，Ai,j表示所述定位像的数据，i＝1:N,j＝1:M，N表示所述全部通道的通道数量，M表示重建后的总图像长度。
5.根据权利要求4图像处理装置，其特征在于，对所述定位像的数据的平均值进行平滑处理包括：
对所述定位像的数据的平均值应用平滑滤波器，通过下式生成所述平滑值：
其中，表示所述平滑值，h表示所述平滑滤波器。
6.根据权利要求5所述的图像处理装置，其特征在于，通过下式生成所述伪影数据：
其中，Sj表示所述伪影数据。
7.存储介质，其特征在于，所述存储介质包括存储的计算机程序，其中，所述计算机程序运行时控制所述存储介质所在的设备执行权利要求1至3中任一项所述的图像处理方法。
8.处理器，其特征在于，所述处理器运行计算机程序，其中，所述计算机程序运行时执行权利要求1至3中任一项所述的图像处理方法。
9.终端，其特征在于，包括：一个或多个处理器，存储器，以及一个或多个计算机程序，其中，所述一个或多个计算机程序被存储在所述存储器中，并且被配置为由所述一个或多个处理器执行，所述一个或多个计算机程序执行权利要求1至3中任一项所述的图像处理方法。
10.计算机程序产品，其特征在于，所述计算机程序产品被有形地存储在计算机可读介质上并且包括计算机可执行指令，所述计算机可执行指令在被执行时使至少一个处理器执行根据权利要求1至3中任一项所述的图像处理方法。
CN201810864673.XA 2018-08-01 2018-08-01 图像处理方法和装置、存储介质、处理器和终端 Active CN110796605B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810864673.XA CN110796605B (zh) 2018-08-01 2018-08-01 图像处理方法和装置、存储介质、处理器和终端

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810864673.XA CN110796605B (zh) 2018-08-01 2018-08-01 图像处理方法和装置、存储介质、处理器和终端

Publications (2)

Publication Number Publication Date
CN110796605A CN110796605A (zh) 2020-02-14
CN110796605B true CN110796605B (zh) 2024-04-23

Family

ID=69425040

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810864673.XA Active CN110796605B (zh) 2018-08-01 2018-08-01 图像处理方法和装置、存储介质、处理器和终端

Country Status (1)

Country Link
CN (1) CN110796605B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108352078A (zh) * 2015-09-15 2018-07-31 上海联影医疗科技有限公司 图像重建系统和方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8244057B2 (en) * 2007-06-06 2012-08-14 Microsoft Corporation Removal of image artifacts from sensor dust
CN103961120B (zh) * 2013-01-31 2018-06-08 Ge医疗系统环球技术有限公司 Ct设备及其使用的图像处理方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108352078A (zh) * 2015-09-15 2018-07-31 上海联影医疗科技有限公司 图像重建系统和方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
何卫红 ; 方向军 ; 彭建春 ; 邓承健 ; 范锟 ; .螺旋CT图像常见伪影的分析和处理.中南医学科学杂志.2011,(06),全文. *

Also Published As

Publication number Publication date
CN110796605A (zh) 2020-02-14

Similar Documents

Publication Publication Date Title
US10874367B2 (en) Angiography
US7379575B2 (en) Method for post- reconstructive correction of images of a computer tomograph
US10083526B2 (en) Radiation tomographic imaging method, apparatus, and program
US7822172B2 (en) Method for hardening correction in medical imaging
JP5028528B2 (ja) Ｘ線ｃｔ装置
EP2691932B1 (en) Contrast-dependent resolution image
US20130089252A1 (en) Method and system for noise reduction in low dose computed tomography
US20060285737A1 (en) Image-based artifact reduction in PET/CT imaging
KR20230153347A (ko) 좁은 시야 ｘ선 이미징 시스템 및 방법
US20140270452A1 (en) Image data processing
US8768045B2 (en) Method for acquiring a 3D image dataset freed of traces of a metal object
JP6293713B2 (ja) 画像処理装置、放射線断層撮影装置並びにプログラム
KR20120138451A (ko) Ｘ선 단층촬영 시스템 및 이를 이용한 산란 보정 방법
US10134157B2 (en) Image generating apparatus, radiation tomography imaging apparatus, and image generating method and program
CN110796605B (zh) 图像处理方法和装置、存储介质、处理器和终端
Al-Antari et al. Denoising images of dual energy X-ray absorptiometry using non-local means filters
US11213260B2 (en) Method and apparatus for correcting cone-beam artifact in cone-beam computed tomography image, and cone-beam computed tomography apparatus including the same
US6931094B2 (en) Methods and systems for smoothing
US6009140A (en) Stair-case suppression for computed tomograph imaging
Gao et al. Optimization of system parameters for modulator design in x-ray scatter correction using primary modulation
US20120177173A1 (en) Method and apparatus for reducing imaging artifacts
CN112446931A (zh) 一种重建数据处理方法、装置、医学成像系统及存储介质
Lin et al. Quantification of radiographic image quality based on patient anatomical contrast-to-noise ratio: a preliminary study with chest images
CN107341836B (zh) 一种ct螺旋扫描图像重建方法及装置
JP7403994B2 (ja) 医用画像処理装置および医用画像処理方法

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