CN111524067B - 一种图像处理方法、装置和设备 - Google Patents

一种图像处理方法、装置和设备 Download PDF

Info

Publication number
CN111524067B
CN111524067B CN202010252169.1A CN202010252169A CN111524067B CN 111524067 B CN111524067 B CN 111524067B CN 202010252169 A CN202010252169 A CN 202010252169A CN 111524067 B CN111524067 B CN 111524067B
Authority
CN
China
Prior art keywords
image
frame image
target
splicing
target frame
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
CN202010252169.1A
Other languages
English (en)
Other versions
CN111524067A (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.)
Neusoft Medical Systems Co Ltd
Beijing Neusoft Medical Equipment Co Ltd
Original Assignee
Neusoft Medical Systems Co Ltd
Beijing Neusoft Medical Equipment 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 Neusoft Medical Systems Co Ltd, Beijing Neusoft Medical Equipment Co Ltd filed Critical Neusoft Medical Systems Co Ltd
Priority to CN202010252169.1A priority Critical patent/CN111524067B/zh
Publication of CN111524067A publication Critical patent/CN111524067A/zh
Application granted granted Critical
Publication of CN111524067B publication Critical patent/CN111524067B/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
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/38—Registration of image sequences
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00—Indexing scheme for image data processing or generation, in general
    • G06T2200/32—Indexing scheme for image data processing or generation, in general involving image mosaicing
    • 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/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular
    • 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
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D10/00—Energy efficient computing, e.g. low power processors, power management or thermal management

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Abstract

本申请提供的图像处理方法，所述图像包括基准图像序列和至少一个待对齐图像序列；所述方法包括：将基准图像序列中相邻的两帧基准帧图像配准对齐后进行拼接，得到第一拼接图像；针对待对齐图像序列中的每帧目标帧图像，将该目标帧图像与该目标帧图像对应的基准帧图像进行配准对齐，得到该目标帧图像的第一配准位移；根据各目标帧图像的第一配准位移、以及该目标帧图像对应的基准帧图像在第一拼接图像中的位置信息，对待对齐图像序列进行拼接，得到第二拼接图像；根据预设算法对第一拼接图像和第二拼接图像进行处理，得到目标图像。本申请提供的图像处理方法、装置和设备，可减少目标图像的伪影。

Description

一种图像处理方法、装置和设备
技术领域
本申请涉及图像处理领域，尤其涉及一种图像处理方法、装置和设备。
背景技术
数字减影血管造影技术(digital subtraction angiogram，简称DSA)是指利用射线对器官两次成像(第一次成像结果称为蒙片，第二成像结果称为造影片)，进而将两次成像结果相减得到只含有血管的目标图像的技术。在利用射线对器官成像时，对于超过探测器视野的器官(例如，脊柱、上下肢等)，两次成像均需要移动机架，对器官进行多段次拍摄才能获得完整血管图。
对器官进行多段次拍摄时，两次成像均涉及位置误差，例如，机器的运动导致的误差和病人身体的挪动导致的误差，这些位置误差使得目标图像存在伪影。因此，如何校正这些误差成为当前亟待解决的问题。
发明内容
有鉴于此，本申请提供一种图像处理方法、装置和设备，以减少目标图像的伪影。
本申请第一方面提供一种图像处理装置，所述图像包括基准图像序列和至少一个待对齐图像序列；所述基准图象序列和所述待对齐图像序列为对同一对象进行指定段次拍摄得到的图象序列；所述方法包括：
将所述基准图像序列中相邻的两帧基准帧图像配准对齐后进行拼接，得到第一拼接图像；
针对待对齐图像序列中的每帧目标帧图像，将该目标帧图像与该目标帧图像对应的基准帧图像进行配准对齐，得到该目标帧图像的第一配准位移；
根据各目标帧图像的第一配准位移、以及该目标帧图像对应的基准帧图像在第一拼接图像中的位置信息，对所述待对齐图像序列进行拼接，得到第二拼接图像；
根据预设算法对所述第一拼接图像和所述第二拼接图像进行处理，得到目标图像。
本申请第二方面提供一种图像处理装置，所述图像包括基准图像序列和至少一个待对齐图像序列；所述基准图象序列和所述待对齐图像序列为对同一对象进行指定段次拍摄得到的图象序列；所述装置第一拼接模块、对齐模块、第二拼接模块和处理模块，其中：
所述第一拼接模块，用于将所述基准图像序列中相邻的两帧基准帧图像配准对齐后进行拼接，得到第一拼接图像；
所述对齐模块，用于针对待对齐图像序列中的每帧目标帧图像，将该目标帧图像与该目标帧图像对应的基准帧图像进行配准对齐，得到该目标帧图像的第一配准位移；
所述第二拼接模块，用于根据各目标帧图像的第一配准位移、以及该目标帧图像对应的基准帧图像在第一拼接图像中的位置信息，对所述待对齐图像序列进行拼接，得到第二拼接图像；
所述处理模块，用于根据预设算法对所述第一拼接图像和所述第二拼接图像进行处理，得到目标图像。
本申请第三方面提供一种计算机可读存储介质，其上存储有计算机程序，所述程序被处理器执行时实现本申请提供的任一图像处理方法的步骤。
本申请第四方面提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现本申请提供的任一图像处理方法的步骤。
本申请提供的图像处理方法、装置和设备，通过将所述基准图像序列中相邻的两帧基准帧图像配准对齐后进行拼接，并针对待对齐图像序列中的每帧目标帧图像，将该目标帧图像与该目标帧图像对应的基准帧图像进行配准对齐，得到该目标帧图像的第一配准位移，进而根据各目标帧图像的第一配准位移、以及该目标帧图像对应的基准帧图像在第一拼接图像中的位置信息，对所述待对齐图像序列进行拼接，得到第二拼接图像，这样，不仅可克服拼接区域的伪影，还能保证待对齐图像序列中的各帧目标帧图像与其对应的基准帧图像是对齐的，使得目标图像伪影少。
附图说明
图1为本申请提供的图像处理方法实施例一的流程图；
图2为本申请一示例性实施例示出的图像序列的示意图；
图3为本申请另一示例性实施例示出的图像序列的示意图；
图4为对图3所示图像序列进行处理后得到的目标图像的示意图；
图5为本申请提供的图像处理方法实施例二的流程图；
图6为本申请一示例性实施例示出的图像处理方法的实现流程图；
图7为利用图6所示方法对图像进行处理的示意图；
图8为本申请一示例性实施例示出的图像处理装置所在图像处理设备的硬件结构图；
图9为本申请提供的图像处理装置实施例一的结构示意图。
具体实施方式
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。
在本申请使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本申请。在本申请和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解，尽管在本申请可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本申请范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
本申请提供一种图像处理方法、装置和设备，以减少目标图像的伪影。
下面给出几个具体的实施例，用以详细介绍本申请的技术方案，下面这几个具体的实施例可以相互结合，对于相同或相似的概念或过程可能在某些实施例中不再赘述。
图1为本申请提供的图像处理方法实施例一的流程图。请参照图1，本实施例提供的图像处理方法，所述图像包括基准图像序列和至少一个待对齐图像序列；所述基准图象序列和所述待对齐图像序列为对同一对象进行指定段次拍摄得到的图象序列；所述方法可以包括：
S101、将所述基准图像序列中相邻的两帧基准帧图像配准对齐后进行拼接，得到第一拼接图像。
具体的，基准图像序列可以为蒙片序列，至少一个待对齐图像序列可以包括造影片序列。下面以基准图像序列为蒙片序列，待对齐图像序列为造影片序列为例进行说明。
图2为本申请一示例性实施例示出的图像序列的示意图。参照图2，图像序列是不同时间、不同方位对目标依序连续获取的系列图像。例如，本申请中，一个图像序列可以是对下肢进行4段拍摄得到的系列图像(此时，该图像序列包括4帧图像)。此外，参照图2，一个图像序列中各帧图像的大小相同，各帧图像之间存在重叠区域，拼接时要克服重叠区域的伪影。下面将给出具体的实施例，用以详细说明该步骤的具体实现过程，此处不再赘述。
例如，一实施例中，在拼接时，可以先开辟一画布，然后配准第n(n从0到Phase-1(拍摄段次)帧蒙片和第n+1帧蒙片的重叠区域，得到配准位移；进而将n帧蒙片的非重叠区域拷贝进画布，并根据配准位移对重叠区域进行融合，以及将第n+1帧的非重叠区域拷贝进画布。
S102、针对待对齐图像序列中的每帧目标帧图像，将该目标帧图像与该目标帧图像对应的基准帧图像进行配准对齐，得到该目标帧图像的第一配准位移。
具体的，可以采用相关的配准对齐方法进行配准对齐。例如，可以采用钢体配准、弹性配准和基于特征点匹配的配准方法进行配准对齐。
S103、根据各目标帧图像的第一配准位移、以及该目标帧图像对应的基准帧图像在第一拼接图像中的位置信息，对所述待对齐图像序列进行拼接，得到第二拼接图像。
具体的，该步骤的具体实现过程，可以包括：
(1)针对所述待对齐图像序列中的每帧目标帧图像，根据该目标帧图像的第一配准位移，以及该目标帧图像对应的基准帧图像在第一拼接图像中的位置信息，确定该目标帧图像在拼接画布中的位置信息。
需要说明的是，目标帧图像对应的基准帧图像指的是基准图像序列中帧号与该目标帧图像帧号相同的基准帧图像。例如，第1帧目标帧图像对应的基准帧图像指的是第1帧基准帧图像。
此外，基准帧图像序列拼接好后，基准帧图像在第一拼接图像中的位置即确定。基准帧图像在第一拼接图像中的位置信息可以包括基准帧图像左上角第一个像素点在第一拼接图像中的行号和列号。需要说明的是，对于帧号大于或等于1的基准帧图像，其左上角第一个像素点为开始融合的第一个像素点。
此外，目标帧图像在拼接画布中的位置信息可以包括目标帧图像左上角第一个像素点在拼接画布中的行号和列号。
进一步地，目标帧图像在拼接画布中的位置信息可以按照如下公式计算：
cx＝mx+dix，cy＝my+diy
其中，mx，my基准帧图像在第一拼接图像中的位置信息；
cx，cy为目标帧图像在拼接画布中的位置信息；
dix、diy为第一配准位移。
(2)利用各目标帧图像在拼接画布的位置信息对各目标帧图像进行拼接，得到第二拼接图像。
具体实现时，可根据各目标帧图像在拼接画布中的位置信息，将各目标帧图像填充到拼接画布中，并对拼接画布中相邻两目标帧图像的重叠区域进行融合处理，得到所述第二拼接图像。
下面以对第一帧目标帧图像和第二帧目标帧图像进行拼接为例，详细说明该步骤的具体实现过程：
例如，一实施例中，可以先根据第一帧目标帧图像的位置信息和第二帧目标帧图像的位置信息，确定两者之间的重叠区域，例如，一实施例中，第一帧目标帧图像的位置信息为(1,1)，第二帧目标帧图像的位置信息为(421,1)，图像大小为500*320，此时，确定两者之间的重叠区域大小为80*320(第一帧目标帧图像最下方的80行和第二帧目标帧图像最上方的80行为两者的重叠区域)。
进一步地，可根据图像大小和重叠区域大小确定拼接图像大小，进而将拼接图像划分为上部区域、中间区域和下部区域。例如，本例中，确定拼接图像大小为(500+500-80)*320，并且划分得到的上部区域为(1～420)*320，中间区域为(421～500)*320，下部区域为(501～920)*320。(需要说明的是，结合上面的例子，拼接图像第一点在画布的位置为(1,1))。
进一步地，将第一帧目标帧图像的非重叠区域填充到上部区域，将对第一帧目标帧图像和第二帧目标帧图像的重叠区域融合处理后的融合图像填充到中间区域，将第二帧目标帧图像的非重叠区域填充到下部区域，得到拼接图像。
下面详细介绍对第一帧目标帧图像和第二帧目标帧图像的重叠区域进行融合处理的过程，具体的，该过程可以包括：
(1)确定第一帧目标帧图像中的重叠区域中的各像素点的权重；
下面将对如何确定一像素点的权重进行详细介绍，此处不再赘述。
(2)依据该像素点的权重，该像素点在第一帧目标帧图像的像素值，以及该像素点在第二帧目标帧图像中的对应像素点的像素值，对第一帧目标帧图像和第二帧目标帧图像的重叠区域进行融合处理。
具体的，本步骤中，按照如下公式进行融合处理：
I(i，j)＝w*I1(i，j)+(1-w)*I2(i1，j1)
其中，I(i，j)为融合图像中像素点(i，j)的像素值；
I1(i，j)为第一帧目标帧图像中像素点(i，j)的像素值；
I2(i1，j1)为第二帧目标帧图像中像素点(i1，j1)的像素值。
需要说明的是，像素点(i1，j1)为第二帧目标帧图像中与第一帧目标帧图像中像素点(i，j)对应的像素点。当将第一帧目标帧图像和第二帧目标图像按位置信息放置在拼接画布中后，两者之间的重叠区域即确定，两者之间像素点之间的对应关系即确定。
S104、根据预设算法对所述第一拼接图像和所述第二拼接图像进行处理，得到目标图像。
本例中，当基准图像序列为蒙片序列，待拼接图像序列为造影片序列时，此时，预设算法可以为相减算法。
图3为本申请另一示例性实施例示出的图像序列的示意图。图4为对图3所示图像序列进行处理后得到的目标图像的示意图。请参照图3和图4，例如，在本例中，先对下肢进行4段拍摄，得到蒙片序列(图3中的M1～M4)，进而注射显影剂，对下肢再进行4段拍摄，得到造影片序列(图3中的C1～C4)。对M1～M4进行拼接，得到如图4中的A图，对C1～C4进行拼接，得到图4中的B图，进而基于A图和B图，得到C图(目标图像，血管减影图像)。从图4中可以看出，通过本申请提供的方法对图像进行处理后，目标图像中的伪影较少。
本实施例提供的方法，通过将所述基准图像序列中相邻的两帧基准帧图像配准对齐后进行拼接，并针对待对齐图像序列中的每帧目标帧图像，将该目标帧图像与该目标帧图像对应的基准帧图像进行配准对齐，得到该目标帧图像的第一配准位移，进而根据各目标帧图像的第一配准位移、以及该目标帧图像对应的基准帧图像在第一拼接图像中的位置信息，对所述待对齐图像序列进行拼接，得到第二拼接图像，这样，不仅可克服拼接区域的伪影，还能保证待对齐图像序列中的各帧目标帧图像与其对应的基准帧图像是对齐的，使得目标图像伪影少。
图5为本申请提供的图像处理方法实施例二的流程图。请参照图5，在上述实施例的基础上，步骤S101，可以包括：
S201、针对所述基准图像序列中相邻的第一基准帧图像和第二基准帧图像，确定所述第一基准帧图像中的第一重叠区域和所述第二基准帧图像中的第二重叠区域。
具体的，有关该步骤的具体实现过程和实现原理可以参见相关技术中的描述，此处不再赘述。
S202、将所述第一重叠区域和所述第二重叠区域进行配准对齐，得到两者之间的第二配准位移。
具体的，该步骤中，可以基于刚体配准、弹性配准和基于特征点匹配的方法，确定第二配准位移。例如，一实施中，可以根据第一重叠区域的灰度梯度图像，选取特征点，进而利用特征点匹配的方法，用爬山法找到每个特征点的位移，进而得到第二配准位移。
S203、利用所述第二配准位移对所述第一基准帧图像和所述第二基准帧图像进行拼接。
所述第一基准帧图像的拍摄时间早于所述第二基准帧图像，该步骤的具体实现过程，可以包括：
(1)针对所述第一重叠区域中的每个参考像素点，依据所述第二配准位移，确定所述第二重叠区域中与该参考像素点对应的目标像素点。
例如，一实施例中，第二配准位移为dx，dy，其中，dx为水平方向的平移量，dy为竖直方向的平移量。此时，针对第一重叠区域中的参考像素点(i，j)，确定第二重叠区域中与该参考像素点对应的目标像素点为(i+dx，j+dy)。
(2)针对每个参考像素点，确定该参考像素点的权重。
一实施例中，可根据该参考像素点距离第一重叠区域上下边界的距离，确定该参考像素点的权重。例如，本例中，可以按照如下公式确定：
W＝d2/(d2+d1)
其中，d2为该参考像素点距离第一重叠区域下边界的距离；d1为该参考像素点距离第一重叠区域上边界的距离。
从上述公式可以看出，参考像素点距离融合边界越远，权重越大。
(3)根据该所述参考像素点的像素值、该参考像素点的权重和所述目标像素点的像素值，对所述第一重叠区域和所述第二重叠进行融合，得到融合图像；
具体的，融合图像中，各像素点的像素值可按照如下公式计算：
I(i，j)＝w*I1(i，j)+(1-w)*I2(i+dx，j+dy)
需要说明的是，在融合的过程中，对于超出图像边界的像素点，可以进行补零处理。
进一步地，由于X光采集的图像都存在灰度台阶，为了让这部分像素不参与融合运算，一实施例中，设置了reject参数(reject＝round(30/1956*row，其中，row为图像高度。)，此时，参考像素点的权重按照如下公式计算：
W＝d2/(d2+d1-reject)
进一步地，此时，在融合时，针对前reject行，融合时只参考第一重叠区域。例如，一实施例中，重叠区域的大小为80*320；此时，针对前30行，融合时，只参考第一重叠区域，针对后50行，融合时参考第一重叠区域和第二重叠区域。即此时，按照如下公式进行融合：
I(i，j)＝I1(i，j)，其中，(i，j)为融合图像前30行的像素点。
I(i，j)＝w*I1(i，j)+(1-w)*I2(i+dx，j+dy)，其中，(i，j)为融合图像剩余的像素点。
(4)将所述第一基准帧图像中的非重叠区域，所述融合图像和第二基准帧图像中的非重叠区域进行拼接，得到所述第一拼接图像。
具体的，可将第一基准帧图像中的非重叠区域，所述融合图像和第二基准帧图像中的非重叠区域依次填充进画布，得到第一拼接图像。
可选地，在本申请一可能的实现方式中，再将所述第一重叠区域和所述第二重叠区域进行配准对齐之前，所述方法还包括：
根据第零帧基准帧图像的目标第一重叠区域的亮度均值，对所述基准图像序列中各基准帧图像的亮度进行调整，以使各基准帧图像的亮度均值等于所述目标第一重叠区域的亮度均值；
在对待对齐图像序列进行拼接之前，所述方法还包括：
对所述待对齐图像序列中各目标帧图像的亮度均值进行调整，以使各目标帧图像的亮度均值等于所述目标第一重叠区域的亮度均值。
例如，一实施例中，第零帧基准帧图像的第一重叠区域的亮度均值为70，此时，就对基准帧图像序列中各基准帧图像的亮度进行调整，以使各基准帧图像的亮度均值等于所述第一重叠区域的亮度均值。例如，第2帧基准帧图像的亮度均值为75，此时，就将第2帧基准帧图像中各像素点的像素值减去5，以使其亮度均值等于70。
再例如，对待对齐图像序列中各目标帧图像的亮度均值进行调整，使得该目标帧图像的亮度均值等于70。
本实施例提供的方法，通过进行亮度调整，可以图像序列中多帧图像亮度不一致的问题，得到亮度一致的目标图像。
图6为本申请一示例性实施例示出图像处理方法的实现流程图。图7为利用图6所示方法对图像进行处理的示意图。请参照图6和图7，具体可以包括：
S601、将第0帧蒙片图像放置在Mask画布上。
例如，一实施例中，将第0帧蒙片图像放置在(1,1)的位置，即将第0帧蒙片图像左上角第一点放置在(1,1)的位置。本例中，图像大小为100*100。
S602、配准第0帧蒙片图像和第0帧造影片图像，得到第一配准位移；以及根据第0帧蒙片图像在Mask画布的位置信息和第一配准位移，确定第0帧目标帧图像在Contrast画布上的位置信息，并根据确定的位置信息将第0帧目标帧图像放置在Contrast画布上。
例如，一实施例中，第一配准位移为dix＝0，diy＝1，此时，就将第0帧造影片放置在(2,1)的位置。
S603、令fn＝1，重复执行以下步骤，直至fn大于(phase-1)。
S604、配准第fn帧蒙片图像与第fn-1帧蒙片图像的重合区域，得到第二配准位移；以及根据第二配准位移对第fn帧蒙片图像与第fn-1帧蒙片图像进行拼接。
本步骤的具体实现原理和实现过程可以参见前面实施例中的描述，此处不再赘述。
下面以配准第0帧蒙片图像与第1帧蒙片图像的重合区域为例进行说明。例如，一实施例中，确定第0帧蒙片图像的第一重合区域为下方20行的区域，第1帧蒙片图像的第二重叠区域为上方20行的区域，经过配准，确定两者之间的第二配准位移为(1，2)，进一步地，本例中，reject参数为5，即第1帧图像中的上面5行不参与融合，此时，融合后，第一帧蒙片图像在第一拼接图像中的位置信息为(86，3)。其中，86＝(100-20)+1+reject，3＝1+2。
S605、配准第fn帧造影片图像与第fn帧蒙片图像，得到第一配准位移；以及根据第fn帧蒙片图像在Mast画布中的位置信息和第一配准位移确定第fn帧造影片图像在Contrast的位置信息，并根据确定的位置信息将第fn帧造影片图放置contrast画布上，且在放置时，对第fn帧造影片图像与第fn-1帧造影片图像的重叠区域进行融合。
下面以第一帧蒙片图像和第一帧造影片图像为例，详细说明书本步骤的具体实现过程。具体的，例如，一实施例中，第二配准位移为(0，1)。结合上面的例子，第一帧蒙片图像的位置信息为(86，3)，此时，确定第一帧造影片图像的位置信息为(87，3)。
这样，第0帧造影片图像的位置信息为(2，1)，第一帧蒙片图像的位置信息为(86，3)，此时，即可基于位置信息将这两帧造影片放置在画布上，并对两者之间的重合区域进行融合处理。其中，基于两者的位置信息，可知，第0帧造影片图像的下方14行与第1帧造影片图像的上方14行重合，且在融合时，第0帧造影片图像的第j+1(j大于等于2)列与第1帧造影片图像的第j列重合。
此外，本例中，还可以记录各重叠区域的起始行号和终止行号，以将其展示给用户，供手动移动拼接图像以便改善拼接效果。
与前述图像处理方法的实施例相对应，本申请还提供了图像处理装置的实施例。
本申请图像处理装置的实施例可以应用在图像处理设备上。装置实施例可以通过软件实现，也可以通过硬件或者软硬件结合的方式实现。以软件实现为例，作为一个逻辑意义上的装置，是通过其所在计算机设备的处理器将存储器中对应的计算机程序指令读取到内存中运行形成的。从硬件层面而言，如图8所示，为本申请一示例性实施例示出的图像处理装置所在图像处理设备的硬件结构图，除了图8所示的存储器810、处理器820、内存830和网络接口840之外，实施例中装置所在的图像处理设备通常根据该图像处理装置的实际功能，还可以包括其他硬件，对此不再赘述。
图9为本申请提供的图像处理装置实施例一的结构示意图。请参照图9，本实施例提供的图像处理装置，所述图像包括基准图像序列和至少一个待对齐图像序列；所述基准图象序列和所述待对齐图像序列为对同一对象进行指定段次拍摄得到的图象序列；所述装置包括第一拼接模块910、对齐模块920、第二拼接模块930和处理模块940，其中：
所述第一拼接模块910，用于将所述基准图像序列中相邻的两帧基准帧图像配准对齐后进行拼接，得到第一拼接图像；
所述对齐模块920，用于针对待对齐图像序列中的每帧目标帧图像，将该目标帧图像与该目标帧图像对应的基准帧图像进行配准对齐，得到该目标帧图像的第一配准位移；
所述第二拼接模块930，用于根据各目标帧图像的第一配准位移、以及该目标帧图像对应的基准帧图像在第一拼接图像中的位置信息，对所述待对齐图像序列进行拼接，得到第二拼接图像；
所述处理模块940，用于根据预设算法对所述第一拼接图像和所述第二拼接图像进行处理，得到目标图像。
本实施例的装置，可用于执行图1所示方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
本申请还提供一种计算机可读存储介质，其上存储有计算机程序，所述程序被处理器执行时实现本申请提供的任一图像处理方法的步骤。
具体的，适合于存储计算机程序指令和数据的计算机可读介质包括所有形式的非易失性存储器、媒介和存储器设备，例如包括半导体存储器设备(例如EPROM、EEPROM和闪存设备)、磁盘(例如内部硬盘或可移动盘)、磁光盘以及CD ROM和DVD-ROM盘。
请继续参照图8，本申请还提供一种图像处理设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现本申请提供的任一图像处理方法的步骤。
以上所述仅为本申请的较佳实施例而已，并不用以限制本申请，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本申请保护的范围之内。

Claims (6)

1.一种图像处理方法，其特征在于，所述图像包括基准图像序列和至少一个待对齐图像序列；所述基准图像序列和所述待对齐图像序列为对同一对象进行指定段次拍摄得到的图像序列；所述方法包括：
对所述基准图像序列中相邻的两帧基准帧图像，将第一基准帧图像中的第一重叠区域与第二基准帧图像中的第二重叠区域进行配准对齐，得到两者之间的第二配准位移，并利用所述第二配准位移对所述第一基准帧图像和所述第二基准帧像进行拼接，得到第一拼接图像；
针对待对齐图像序列中的每帧目标帧图像，将该目标帧图像与该目标帧图像对应的基准帧图像进行配准对齐，得到该目标帧图像的第一配准位移；
根据各目标帧图像的第一配准位移、以及该目标帧图像对应的基准帧图像在第一拼接图像中的位置信息，确定该目标帧图像在拼接画布中的位置信息；
根据各目标帧图像在拼接画布中的位置信息，将各目标帧图像填充到拼接画布中，并对拼接画布中相邻两目标帧图像的重叠区域进行融合处理，得到第二拼接图像；
根据预设算法对所述第一拼接图像和所述第二拼接图像进行处理，得到目标图像。
2.根据权利要求1所述的方法，其特征在于，所述第一基准帧图像的拍摄时间早于所述第二基准帧图像，所述利用所述第二配准位移对所述第一基准帧图像和所述第二基准帧图像进行拼接，包括：
针对所述第一重叠区域中的每个参考像素点，依据所述第二配准位移，确定所述第二重叠区域中与该参考像素点对应的目标像素点；
针对每个参考像素点，确定该参考像素点的权重；
根据该所述参考像素点的像素值、该参考像素点的权重和所述目标像素点的像素值，对所述第一重叠区域和所述第二重叠进行融合，得到融合图像；
将所述第一基准帧图像中的非重叠区域，所述融合图像和所述第二基准帧图像中的非重叠区域进行拼接，得到所述第一拼接图像。
3.根据权利要求1所述的方法，其特征在于，在将所述第一重叠区域和所述第二重叠区域进行配准对齐之前，所述方法还包括：
根据第零帧基准帧图像的目标第一重叠区域的亮度均值，对所述基准图像序列中各基准帧图像的亮度进行调整，以使各基准帧图像的亮度均值等于所述目标第一重叠区域的亮度均值；
所述对所述待对齐图像序列进行拼接之前，所述方法还包括：
对所述待对齐图像序列中各目标帧图像的亮度均值进行调整，以使各目标帧图像的亮度均值等于所述目标第一重叠区域的亮度均值。
4.一种图像处理装置，其特征在于，所述图像包括基准图像序列和至少一个待对齐图像序列；所述基准图像序列和所述待对齐图像序列为对同一对象进行指定段次拍摄得到的图像序列；所述装置包括第一拼接模块、对齐模块、第二拼接模块和处理模块，其中：
所述第一拼接模块，用于对所述基准图像序列中相邻的两帧基准帧图像，将第一基准帧图像中的第一重叠区域与第二基准帧图像中的第二重叠区域进行配准对齐，得到两者之间的第二配准位移，并利用所述第二配准位移对所述第一基准帧图像和所述第二基准帧像进行拼接，得到第一拼接图像；
所述对齐模块，用于针对待对齐图像序列中的每帧目标帧图像，将该目标帧图像与该目标帧图像对应的基准帧图像进行配准对齐，得到该目标帧图像的第一配准位移；
所述第二拼接模块，用于根据各目标帧图像的第一配准位移、以及该目标帧图像对应的基准帧图像在第一拼接图像中的位置信息，确定该目标帧图像在拼接画布中的位置信息；根据各目标帧图像在拼接画布中的位置信息，将各目标帧图像填充到拼接画布中，并对拼接画布中相邻两目标帧图像的重叠区域进行融合处理，得到第二拼接图像；
所述处理模块，用于根据预设算法对所述第一拼接图像和所述第二拼接图像进行处理，得到目标图像。
5.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述程序被处理器执行时实现本申请提供的任一图像处理方法的步骤。
6.一种图像处理设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现本申请提供的任一图像处理方法的步骤。
CN202010252169.1A 2020-04-01 2020-04-01 一种图像处理方法、装置和设备 Active CN111524067B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010252169.1A CN111524067B (zh) 2020-04-01 2020-04-01 一种图像处理方法、装置和设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010252169.1A CN111524067B (zh) 2020-04-01 2020-04-01 一种图像处理方法、装置和设备

Publications (2)

Publication Number Publication Date
CN111524067A CN111524067A (zh) 2020-08-11
CN111524067B true CN111524067B (zh) 2023-09-12

Family

ID=71902534

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010252169.1A Active CN111524067B (zh) 2020-04-01 2020-04-01 一种图像处理方法、装置和设备

Country Status (1)

Country Link
CN (1) CN111524067B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112991178B (zh) * 2021-03-24 2022-08-02 北京百度网讯科技有限公司 图像拼接方法、装置、设备和介质
CN114245006B (zh) * 2021-11-30 2023-05-23 联想(北京)有限公司 一种处理方法、装置及系统
CN114287963A (zh) * 2021-12-30 2022-04-08 重庆海扶医疗科技股份有限公司 图像处理方法、装置、电子设备及计算机可读介质
CN115517692B (zh) * 2022-11-25 2023-03-14 苏州一目万相科技有限公司 成像系统及其控制方法和控制装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101086766A (zh) * 2007-06-04 2007-12-12 北京航空航天大学 基于波形匹配的指纹图像帧序列的拼接方法
CN101901481A (zh) * 2010-08-11 2010-12-01 深圳市蓝韵实业有限公司 一种图像拼接方法
CN104680486A (zh) * 2013-11-29 2015-06-03 上海联影医疗科技有限公司 一种非刚体的配准方法
CN106875354A (zh) * 2017-01-20 2017-06-20 北京东软医疗设备有限公司 图像去噪方法、装置及设备
CN106886977A (zh) * 2017-02-08 2017-06-23 徐州工程学院 一种多图自动配准及融合拼接方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9531962B2 (en) * 2014-07-24 2016-12-27 Samsung Electronics Co., Ltd. Image set alignment and combination processing
KR102103252B1 (ko) * 2015-12-28 2020-04-22 후아웨이 테크놀러지 컴퍼니 리미티드 이미지 융합 방법 및 장치, 및 단말 디바이스
AU2016266019A1 (en) * 2016-11-30 2018-06-14 Canon Kabushiki Kaisha Image registration method
JP6723173B2 (ja) * 2017-02-10 2020-07-15 富士フイルム株式会社 画像処理装置、方法およびプログラム

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101086766A (zh) * 2007-06-04 2007-12-12 北京航空航天大学 基于波形匹配的指纹图像帧序列的拼接方法
CN101901481A (zh) * 2010-08-11 2010-12-01 深圳市蓝韵实业有限公司 一种图像拼接方法
CN104680486A (zh) * 2013-11-29 2015-06-03 上海联影医疗科技有限公司 一种非刚体的配准方法
CN106875354A (zh) * 2017-01-20 2017-06-20 北京东软医疗设备有限公司 图像去噪方法、装置及设备
CN106886977A (zh) * 2017-02-08 2017-06-23 徐州工程学院 一种多图自动配准及融合拼接方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈蔓 等.基于SIFT字典学习的引导滤波多聚焦图像融合.《哈尔滨工业大学学报》.2018,第50卷(第11期),第59-66页. *

Also Published As

Publication number Publication date
CN111524067A (zh) 2020-08-11

Similar Documents

Publication Publication Date Title
CN111524067B (zh) 一种图像处理方法、装置和设备
CN106056537B (zh) 一种医学图像拼接方法及装置
CN104268846B (zh) 图像拼接方法及装置
CN105389788B (zh) Pet多床位图像的重建方法及装置、合并方法及装置
JPS63149772A (ja) 物体の連続２画像から減算像を決定する方法および装置
CN104992408B (zh) 用于用户终端的全景图像生成方法和装置
JP2005100407A (ja) 複数のソース画像からパノラマ画像を作成するシステム及び方法
US8433114B2 (en) Method and system for elastic composition of medical imaging volumes
US20110058722A1 (en) Using non-attenuation corrected pet emission images to compensate for incomplete anatomic images
AU2005280315B2 (en) Method and system for motion correction in a sequence of images
US20060116570A1 (en) Reducing effects of rotational motion
CN108074270B (zh) 一种pet衰减校正方法和装置
CN107456236B (zh) 一种数据处理方法及医疗扫描系统
US7911208B2 (en) Methods for rectification of B0 inhomogeneity effects in magnetic resonance images
CN111583120B (zh) 图像拼接方法、装置、设备和存储介质
JP2002092589A (ja) 画像の位置合わせ方法
EA024855B1 (ru) Способ получения субтракционного ангиографического изображения
US8068665B2 (en) 3D-image processing apparatus, 3D-image processing method, storage medium, and program
EP0919858B1 (en) Method for reconstructing a radiation image of a body from partial radiation images
CN107622475A (zh) 图像拼接中的灰度校正方法
WO2015075718A1 (en) Method and system for tracking a region in a video image
CN108648145B (zh) 图像拼接方法及装置
JP6329262B2 (ja) 対象物の全体像の作成方法及び顕微鏡画像の作成方法
KR20120036908A (ko) 스테레오 화상 촬영장치 및 그 방법
US8111068B2 (en) Reducing effects of rotational motion

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