Nothing Special   »   [go: up one dir, main page]

CN113538251B - 确定医学影像拼接异常的方法及装置 - Google Patents

确定医学影像拼接异常的方法及装置 Download PDF

Info

Publication number
CN113538251B
CN113538251B CN202111083284.1A CN202111083284A CN113538251B CN 113538251 B CN113538251 B CN 113538251B CN 202111083284 A CN202111083284 A CN 202111083284A CN 113538251 B CN113538251 B CN 113538251B
Authority
CN
China
Prior art keywords
image
image sequence
spliced
continuity
sheet
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
CN202111083284.1A
Other languages
English (en)
Other versions
CN113538251A (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.)
Zhejiang Taimei Medical Technology Co Ltd
Original Assignee
Zhejiang Taimei Medical 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 Zhejiang Taimei Medical Technology Co Ltd filed Critical Zhejiang Taimei Medical Technology Co Ltd
Priority to CN202111083284.1A priority Critical patent/CN113538251B/zh
Publication of CN113538251A publication Critical patent/CN113538251A/zh
Application granted granted Critical
Publication of CN113538251B publication Critical patent/CN113538251B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4038Image mosaicing, e.g. composing plane images from plane sub-images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4007Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00ICT specially adapted for the handling or processing of medical images
    • G16H30/20ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Epidemiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • Image Processing (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本申请提供了一种确定医学影像拼接异常的方法,包括:获取在第一分段床位扫描的第一影像序列,并获取在第二分段床位扫描的第二影像序列,其中第一影像序列与第二影像序列为待拼接的影像序列;确定拼接张,拼接张为第二影像序列中的影像,拼接张用于拼接第一影像序列和第二影像序列;分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果;根据连续性分析结果,确定第一影像序列和第二影像序列是否拼接异常。本申请的确定医学影像拼接异常的方法通过判断多床位拼接影像的连续性,实现了对医学影像拼接异常的自动提醒,帮助医生和操作人员快速确认多床位影像拼接的正确性,提高影像处理的效率和质量。

Description

确定医学影像拼接异常的方法及装置
技术领域
本申请涉及医学影像技术领域,并且更为具体地,涉及一种确定医学影像拼接异常的方法及装置。
背景技术
DICOM(医学数字成像和通信,Digital Imaging and Communications inMedicine)是医学图像和相关信息的国际标准(ISO 12052)。它定义了质量能满足临床需要的可用于数据交换的医学图像格式。
在人体影像扫描的过程中,由于人体的轴向长度大于扫描仪的长度,需要逐个床位扫描,因此需要对多个床位扫描的影像进行拼接。现有的多床位影像拼接的方式大多为人工拼接或由扫描设备完成,医生和操作人员的工作量大,在拼接的过程中很有可能因为没有注意影像的拼接位置关系而造成影像的错误拼接。
发明内容
有鉴于此,本申请实施例提供了一种确定医学影像拼接异常的方法及装置,能够针对医学影像多床位拼接过程中出现错误拼接的情况给出提醒。
第一方面,本申请的实施例提供了一种确定医学影像拼接异常的方法,包括:获取在第一分段床位扫描的第一影像序列,并获取在第二分段床位扫描的第二影像序列,其中第一影像序列与第二影像序列为待拼接的影像序列;确定拼接张,拼接张为第二影像序列中的影像,拼接张用于拼接第一影像序列和第二影像序列;分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果;根据连续性分析结果,确定第一影像序列和第二影像序列是否拼接异常。
在本申请某些实施例中,确定拼接张包括:若第一影像序列与第二影像序列无交叠,选取第二影像序列中距离第一影像序列最近的第二影像作为拼接张;若第一影像序列与第二影像序列有交叠,选取第二影像序列中与第一影像序列无交叠的影像中的距离第一影像序列最近的第二影像作为拼接张。
在本申请某些实施例中,分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,包括:根据拼接张到第一影像序列中与拼接张相邻的第一影像的第一距离,以及拼接张到第二影像序列中与拼接张相邻的第二影像的第二距离,确定拼接张的距离差异;若距离差异大于第一预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
在本申请某些实施例中,分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,包括:采用插值法针对拼接张所处的位置进行拟合,得到拟合拼接张;确定拼接张和拟合拼接张的第一图像差异;若第一图像差异大于第二预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
在本申请某些实施例中,采用插值法针对拼接张所处的位置进行拟合,得到拟合拼接张,包括:根据cubic算法将与拼接张相邻的n张第一影像和与拼接张相邻的n张第二影像进行拟合,得到拟合拼接张,其中n为大于0的整数。
在本申请某些实施例中,分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,还包括:根据拼接张到第一影像序列中与拼接张相邻的第一影像的第一距离,以及拼接张到第二影像序列中与拼接张相邻的第二影像的第二距离,确定拼接张的距离差异;将距离差异和第一图像差异进行加权,得到第一连续性差异;若第一连续性差异大于第三预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
在本申请某些实施例中,分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,包括:根据cubic算法确定拼接张与第一影像序列中与拼接张相邻的至少一个第一影像的第二图像差异;根据cubic算法确定拼接张与第二影像序列中与拼接张相邻的至少一个第二影像的第三图像差异;根据第二图像差异和第三图像差异,确定第二连续性差异;若第二连续性差异大于预设第四阈值,确定第一影像序列和第二影像序列之间不具备连续性。
第二方面,本申请的实施例提供了一种确定医学影像拼接异常的装置,包括:获取模块,用于获取在第一分段床位扫描的第一影像序列,并获取在第二分段床位扫描的第二影像序列,其中第一影像序列与第二影像序列为待拼接的影像序列;第一确定模块,用于确定拼接张,拼接张为第二影像序列中的影像,拼接张用于拼接第一影像序列和第二影像序列;分析模块,用于分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果;第二确定模块,用于根据连续性分析结果,确定第一影像序列和第二影像序列是否拼接异常。
第三方面,本申请的实施例提供了一种计算机可读存储介质,存储介质存储有计算机程序,计算机程序用于执行上述第一方面的确定医学影像拼接异常的方法。
第四方面,本申请的实施例提供了一种电子设备,包括:处理器;用于存储处理器可执行指令的存储器,其中,处理器用于执行上述第一方面的确定医学影像拼接异常的方法。
本申请的确定医学影像拼接异常的方法通过判断多床位拼接影像的连续性,实现了对医学影像拼接异常的自动提醒,帮助医生和操作人员快速确认多床位影像拼接的正确性,提高影像处理的效率和质量。
附图说明
图1为本申请一实施例所提供的确定医学影像拼接异常的一种实施环境的示意图。
图2为本申请一实施例提供的确定医学影像拼接异常的方法的流程示意图。
图3为本申请一实施例提供的确定拼接张的流程示意图。
图4为本申请另一实施例提供的分析连续性得到连续性分析结果的流程示意图。
图5为本申请另一实施例提供的分析连续性得到连续性分析结果的流程示意图。
图6为本申请另一实施例提供的分析连续性得到连续性分析结果的流程示意图。
图7为本申请另一实施例提供的分析连续性得到连续性分析结果的流程示意图。
图8为本申请另一实施例提供的分析连续性得到连续性分析结果的流程示意图。
图9a所示为本申请一示例性实施例提供的医学影像多床位扫描过程的示意图。
图9b所示为本申请一示例性实施例提供的具备连续性的影像的示意图。
图9c所示为本申请一示例性实施例提供的多床位影像不同交叠情况的示意图。
图9d所示为本申请一示例性实施例提供的不同交叠情况下确定的拼接张的示意图。
图10为本申请一实施例提供的确定医学影像拼接异常的装置的结构示意图。
图11所示为本申请一示例性实施例提供的用于执行确定医学影像拼接异常的方法的电子设备的框图。
具体实施方式
下面将参照附图更详细地描述本申请的实施例。虽然附图中显示了本申请的某些实施例,然而应当理解的是,本申请可以通过各种形式来实现,而且不应该被解释为限于这里阐述的实施例,相反提供这些实施例是为了更加透彻和完整地理解本申请。应当理解的是,本申请的附图及实施例仅用于示例性作用,并非用于限制本申请的保护范围。
本申请使用的术语“包括”及其变形是开放性包括,即“包括但不限于”。术语“根据”是“至少部分地根据”。术语“一个实施例”表示“至少一个实施例”;术语“另一实施例”表示“至少一个另外的实施例”。其它术语的相关定义将在下文描述中给出。
本申请的实施例可以用于医学影像阅片系统(Medical Imaging ReadingSystem,MIRS)中。在影像图片管理方面,医学影像阅片系统支持多中心的影像上传,支持影像查询,并对上传影像进行审核和质控管理。在阅片管理方面,支持阅片流程设计、多层级阅片的分配、跟踪和查询,并且支持多重阅片。在整个业务流程中,对影像的上传审核阅片进行智能的统计管理,实时跟进影像状态和阅片进度。
医学影像,也可以称为医学图像,可以为电子计算机断层成像(ComputedTomography,CT)、核磁共振成像(Magnetic Resonance Imaging,MRI)、正电子发射型计算机断层显像(Positron Emission Computed Tomography,PET)等医学影像,本申请实施例对此不作具体限定。
DICOM图像文件结构主要由文件头和数据集合两大部分组成。文件头用于表示该文件是否为DICOM图像文件。数据集合包括:受试者(例如,病人)的姓名、成像方式、图像大小之类的与成像实例相关的数据,以及成像实例的图像像素数据。本申请实施例并不限定医学影像的具体形式,可以是原始医学图像,也可以是经过预处理后的医学图像,还可以是原始医学图像的一部分。
患者在做扫描检查时,一般是患者仰卧、侧卧或者俯卧于床板上,通过扫描仪的进出孔径进行横断位扫描。图9a所示为本申请一示例性实施例提供的医学影像多床位扫描过程的示意图。如图9a所示,由于扫描仪的孔径轴向长度往往小于人体的轴向长度,采样区域无法一次性覆盖到全部的人体轴向范围,因此需要采用多床位扫描的方式。
对患者进行多床位扫描后,需要对各个床位的扫描影像进行拼接。申请人在研究多床位影像拼接的过程中发现,现有的多床位影像拼接的方式大多为人工拼接或由扫描设备完成,例如医生或操作人员将影像的状位信息和序列信息等复杂的参数输入给设备,再由设备进行拼接。对于这种拼接方式,一方面,医生和操作人员的工作量大,一般不会仔细排查生成的拼接影像,在拼接的过程中很有可能因为没有注意影像的拼接位置关系而造成影像的错误拼接;另一方面,在拼接的过程中,并没有根据影像的切层连续性进行分析,在拼接异常时无法给出拼接异常提示。因此,拼接影像的过程中经常出现错误拼接的情况。
为了解决上述问题,本申请提供了一种确定医学影像拼接异常的方法。
图1是本申请实施例所提供的一种实施环境的示意图。该实施环境包括CT扫描仪110和计算机设备120。
计算机设备120可以从CT扫描仪110处获取医学影像。例如,计算机设备120可以通过有线网络或无线网络与CT扫描仪110进行通信。
CT扫描仪110用于对人体组织进行X线扫描,得到人体组织的CT医学影像。在一实施例中,通过CT扫描仪110分别对人体的头部、躯干部和腿部进行多床位扫描,并进行拼接,可以得到拼接后的人体CT医学影像。
计算机设备120可以是通用型计算机或者由专用的集成电路组成的计算机装置等,本申请实施例对此不做限定。例如,计算机设备120可以是平板电脑等移动终端设备,或者也可以是个人计算机(Personal Computer,PC),比如膝上型便携计算机和台式计算机等等。
本领域技术人员可以知晓,上述CT扫描仪110的数量可以一个或多个,其类型可以相同或者不同。比如上述CT扫描仪110可以为一个,或者上述CT扫描仪110为几十个或几百个,或者更多数量。此外,上述计算机设备120的数量可以一个或多个,其类型可以相同或者不同。比如上述计算机设备120可以为一个,或者上述计算机设备120为几十个或几百个,或者更多数量。本申请实施例对CT扫描仪110的数量和设备类型以及计算机设备120的数量和设备类型不加以限定。
在一些可选的实施例中,计算机设备120从CT扫描仪110处获取医学影像,获取在第一分段床位扫描的第一影像序列,并获取在第二分段床位扫描的第二影像序列,确定拼接张,分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,并根据连续性分析结果,确定第一影像序列和第二影像序列是否拼接异常。
图2示意性示出了本发明一实施例提供的确定医学影像拼接异常的方法的流程示意图。图2所述的方法由计算设备(例如,服务器)来执行,但本申请实施例不以此为限。服务器可以是一台服务器,或者由若干台服务器组成,或者是一个虚拟化平台,或者是一个云计算服务中心,本申请实施例对此不作限定。如图2所示,该方法包括如下内容:
S210:获取在第一分段床位扫描的第一影像序列,并获取在第二分段床位扫描的第二影像序列,其中第一影像序列与第二影像序列为待拼接的影像序列。
由于扫描仪的孔径轴向长度小于人体长度,在对人体进行扫描时需要逐床分段扫描,如图9a所示。例如,第一分段床位可以为某一患者头部所在的位置,第二分段床位可以为该患者的胸部所在的位置,相应地,第一影像序列为针对该患者头部区域扫描的多张第一影像,第二影像序列为针对该患者胸部区域扫描的多张第二影像,在完成扫描后需要将第一影像序列和第二影像序列进行拼接,即第一影像序列与第二影像序列为待拼接的影像序列。
本申请的第一影像序列和第二影像序列也可以为其他人体部位的分段扫描影像,例如第一影像序列为胸部的分段扫描影像、第二影像序列为头部的扫描影像,或第一影像序列为腿部的分段扫描影像、第二影像序列为躯干的扫描影像,本申请对第一影像序列和第二影像序列所扫描的人体部位不做限定。
S220:确定拼接张,拼接张为第二影像序列中的影像,拼接张用于拼接第一影像序列和第二影像序列。
由于本申请的第一影像序列和第二影像序列仅为区分不同分段床位扫描的影像,因此也可以将第一影像序列中的影像作为拼接张,用于拼接第一影像序列和第二影像序列。
在一实施例中,确定拼接张包括:若第一影像序列与第二影像序列无交叠,选取第二影像序列中距离第一影像序列最近的第二影像作为拼接张;若第一影像序列与第二影像序列有交叠,选取第二影像序列中与第一影像序列无交叠的影像中的距离第一影像序列最近的第二影像作为拼接张。
对于拼接正确的医学影像,后一个床位的扫描范围应和前一个床位匹配上,即位置无缝对接。但实际扫描过程中可能存在两相邻床位有交叠的情况或者存在间隙的情况。在这种情况,通过选择合适的拼接张,能够舍弃掉一些空余的图像,从而能够减少后续处理的计算量。
图9c所示为本申请一示例性实施例提供的多床位影像不同交叠情况的示意图,图9d所示为本申请一示例性实施例提供的不同交叠情况下确定的拼接张的示意图。
如图9c和图9d所示,本实施例可以根据不同交叠情况确定拼接张:1)当第一影像序列与第二影像序列无交叠且无间隙时,选取第二影像序列中距离第一影像序列最近的第二影像作为拼接张;2)当第一影像序列与第二影像序列无交叠且有间隙时,同样选取第二影像序列中距离第一影像序列最近的第二影像作为拼接张;3)当第一影像序列与第二影像序列有交叠时,选取第二影像序列中与第一影像序列无交叠的影像中的距离第一影像序列最近的第二影像作为拼接张。影像间的距离可以根据各个影像的切片坐标(slicelocation)进行计算。
本实施例在无交叠和有交叠两种情况下对拼接张进行选取,使拼接张的选取更全面。
S230:分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果。
图9b所示为本申请一示例性实施例提供的具备连续性的影像的示意图。由图9b可见,正确拼接时的人体影像具备连续性,相邻影像的差异较小,因此可以通过分析拼接张及与拼接张相邻的影像的连续性,得到连续性分析结果。
在一实施例中,分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,包括:根据拼接张到第一影像序列中与拼接张相邻的第一影像的第一距离,以及拼接张到第二影像序列中与拼接张相邻的第二影像的第二距离,确定拼接张的距离差异;若距离差异大于第一预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
具体地,设拼接张到第一影像序列中与拼接张相邻的第一影像的第一距离为L,拼接张到第二影像序列中与拼接张相邻的第二影像的第二距离为l,拼接张的距离差异T1的计算公式如下:
Figure 426000DEST_PATH_IMAGE001
拼接张的距离差异T1以L和l的比例来进行评估,距离差异T1的大小影响相邻影像的连续性分析。当L小于2l时,距离差异T1对相邻影像连续性分析几乎无影响,认为可以忽略距离差异T1对相邻影像连续性分析的影响;当L大于4l时,距离差异T1对相邻影像连续性分析的影响较大。因此,可以设置第一预设阈值为4l,若距离差异T1大于第一预设阈值,则确定第一影像序列和第二影像序列之间不具备连续性。
本实施例将影像间的距离作为分析影像连续性的因素,使连续性分析结果更加全面、准确。
在一实施例中,分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,包括:采用插值法针对拼接张所处的位置进行拟合,得到拟合拼接张;确定拼接张和拟合拼接张的第一图像差异;若第一图像差异大于第二预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
连续性分析也可以采用插值法针对拼接张所处的位置进行拟合,得到拟合拼接张,如果第一影像序列和第二影像序列具备连续性,则拟合拼接张与拼接张的第一图像差异不应过大。因此,可以预设第二预设阈值,当第一图像差异大于第二预设阈值时,确定第一影像序列和第二影像序列之间不具备连续性。第一图像差异可以根据拼接张的DICOM值、第一影像序列中的距离拼接张最近的至少一个第一影像的DICOM值以及第二影像序列中的距离拼接张最近的至少一个第二影像的DICOM值进行计算,例如确定拟合拼接张的各像素的DICOM值的可能区间,根据拼接张的各像素的DICOM值超出可能区间的像素的比例计算第一图像差异。在CT影像中DICOM值为CT值,在PET影像中DICOM值为PET值。插值法可以为cubic算法,也可以为双线性插值法或最近邻插值法。
本实施例将拼接张和拟合拼接张之间的差异作为分析影像连续性的因素,使连续性分析结果更加全面、准确。
在一实施例中,采用插值法针对拼接张所处的位置进行拟合,得到拟合拼接张,包括:根据cubic算法将与拼接张相邻的n张第一影像和与拼接张相邻的n张第二影像进行拟合,得到拟合拼接张,其中n为大于0的整数。
以n=2的情况对本实施例进行说明。设与拼接张相邻的两张第一影像的编号为k-2,k-1,拼接张的编号为k,与拼接张相邻的两张第二影像的编号为k+1,k+2,拟合拼接张的编号为k'。
对拟合拼接张k'的各像素通过cubic算法将前两张和后两张进行拟合,算法为
Figure 20930DEST_PATH_IMAGE002
其中,m和n为拟合取值的方向,如果均为0,在该方向为正向拟合,如果有一个不为0则在该方向为斜向拟合,Vk'mn(i,j)为拟合拼接张第(i,j)个像素的拟合DICOM值,将通过其附近对角线的区域进行计算评估其可能值。因为m、n各有3种取值,因此一共可以得到Vk'mn(i,j)的9个可能值,根据9个Vk'mn(i,j)可能值建立DICOM值的可能区间,例如选取9个Vk'mn(i,j)可能值中的极值 max[Vk'mn(i,j)]和min[Vk'mn(i,j)],得到极值区间。根据分布的可能性以正态分布计算,可根据拼接张k的第(i,j)个像素的DICOM值Vkmn(i,j)落在可能区间内的比例的可能性得到损失函数,按照如下公式计算:
Figure 688672DEST_PATH_IMAGE003
对拼接张k和拟合拼接张k'的所有像素点的loss值进行加权平均,得到第一图像差异LOSS,即:
Figure 436179DEST_PATH_IMAGE004
本实施例采用cubic算法,能够准确计算出拼接张和拟合拼接张之间的差异。
在一实施例中,分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,还包括:根据拼接张到第一影像序列中与拼接张相邻的第一影像的第一距离,以及拼接张到第二影像序列中与拼接张相邻的第二影像的第二距离,确定拼接张的距离差异;将距离差异和第一图像差异进行加权,得到第一连续性差异;若第一连续性差异大于第三预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
连续性分析结果也可以综合考虑距离差异T1距离以及第一图像差异LOSS进行判断,例如可以根据影像的不同模态赋予距离差异T1距离和第一图像差异LOSS不同的权重,得到第一连续性差异ERRORCOST,算法为:
Figure 23018DEST_PATH_IMAGE005
其中,α和β为权重系数,根据影像的模态不同而有差别,距离差异T1和第一图像差异LOSS的计算方法可参见上述实施例,在此不再赘述。当第一连续性差异ERRORCOST大于第三预设阈值时,确定第一影像序列和第二影像序列之间不具备连续性。
本实施例综合考虑距离差异和第一图像差异的权重,使连续性分析结果更加全面、准确。
在一实施例中,分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,包括:根据cubic算法确定拼接张与第一影像序列中与拼接张相邻的至少一个第一影像的第二图像差异;根据cubic算法确定拼接张与第二影像序列中与拼接张相邻的至少一个第二影像的第三图像差异;根据第二图像差异和第三图像差异,确定第二连续性差异;若第二连续性差异大于预设第四阈值,确定第一影像序列和第二影像序列之间不具备连续性。
以与拼接张前后各相邻两张的情况对本实施例进行说明。设与拼接张相邻的两个第一影像的编号为k-2,k-1,拼接张的编号为k,与拼接张相邻的两个第二影像的编号为k+1,k+2,则有:
根据cubic算法计算第二图像差异Vk'mn(i,j)between的值:
Figure 511768DEST_PATH_IMAGE006
其中,Vk'mn(i,j)between为拼接张k与第一影像序列中与拼接张相邻的两张第一影像(即k-2、k-1)的第二图像差异。
根据cubic算法计算第三图像差异Vk'mn(i,j)in的值:
Figure 717622DEST_PATH_IMAGE007
其中,Vk'mn(i,j)in为拼接张k与第二影像序列中与拼接张相邻的两张第二影像(即k+1、k+2)的第三图像差异。
如果第一影像序列与第一影像序列为正常拼接,第二图像差异Vk'mn(i,j)between和第三图像差异Vk'mn(i,j)in应该不会有较大差别。第二连续性差异DIFF的计算方法为:
首先根据第二图像差异Vk'mn(i,j)between和第三图像差异Vk'mn(i,j)in计算difference值:
Figure 821100DEST_PATH_IMAGE008
再对各像素点的difference值进行加权计算,得到第二连续性差异DIFF,即:
Figure 844420DEST_PATH_IMAGE009
可以针对不同模态的影像设定不同的第四预设阈值,如果第二连续性差异DIFF的值大于预设第四阈值,则确定第一影像序列和第二影像序列之间不具备连续性。
本实施例将同床位相邻影像之间的连续性差异和不同床位相邻影像之间的连续性差异的差异作为分析影像连续性的因素,使连续性分析结果更加全面、准确。
S240:根据连续性分析结果,确定第一影像序列和第二影像序列是否拼接异常。
如果确定第一影像序列和第二影像序列之间不具备连续性,则确定第一影像序列和第二影像序列拼接异常,发出拼接异常提醒。
本申请的确定医学影像拼接异常的方法通过判断多床位拼接影像的连续性,实现了对医学影像拼接异常的自动提醒,帮助医生和操作人员快速确认多床位影像拼接的正确性,提高影像处理的效率和质量。
图3为本申请一实施例提供的确定拼接张的流程示意图。
S310:若第一影像序列与第二影像序列无交叠,选取第二影像序列中距离第一影像序列最近的第二影像作为拼接张。
S320:若第一影像序列与第二影像序列有交叠,选取第二影像序列中与第一影像序列无交叠的影像中的距离第一影像序列最近的第二影像作为拼接张。
关于S310~S320的具体内容,可以参考上述实施例中的描述,为避免重复,此处不再赘述。
本实施例在无交叠和有交叠两种情况下对拼接张进行选取,使拼接张的选取更全面。
图4为本申请另一实施例提供的分析连续性得到连续性分析结果的流程示意图。
S410:根据拼接张到第一影像序列中与拼接张相邻的第一影像的第一距离,以及拼接张到第二影像序列中与拼接张相邻的第二影像的第二距离,确定拼接张的距离差异。
S420:若距离差异大于第一预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
关于S410~S420的具体内容,可以参考上述实施例中的描述,为避免重复,此处不再赘述。
本实施例将影像间的距离作为分析影像连续性的因素,使连续性分析结果更加全面、准确。
图5为本申请另一实施例提供的分析连续性得到连续性分析结果的流程示意图。
S510:采用插值法针对拼接张所处的位置进行拟合,得到拟合拼接张;
S520:确定拼接张和拟合拼接张的第一图像差异;
S530:若第一图像差异大于第二预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
关于S510~S530的具体内容,可以参考上述实施例中的描述,为避免重复,此处不再赘述。
本实施例将拼接张和拟合拼接张之间的差异作为分析影像连续性的因素,使连续性分析结果更加全面、准确。
图6为本申请另一实施例提供的分析连续性得到连续性分析结果的流程示意图。
S610:根据cubic算法将与拼接张相邻的n张第一影像和与拼接张相邻的n张第二影像进行拟合,得到拟合拼接张,其中n为大于0的整数。
S620:确定拼接张和拟合拼接张的第一图像差异;
S630:若第一图像差异大于第二预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
关于S610~S630的具体内容,可以参考上述实施例中的描述,为避免重复,此处不再赘述。
本实施例采用cubic算法,能够准确计算出拼接张和拟合拼接张之间的差异。
图7为本申请另一实施例提供的分析连续性得到连续性分析结果的流程示意图。
S710:根据拼接张到第一影像序列中与拼接张相邻的第一影像的第一距离,以及拼接张到第二影像序列中与拼接张相邻的第二影像的第二距离,确定拼接张的距离差异;
S720:将距离差异和第一图像差异进行加权,得到第一连续性差异;
S730:若第一连续性差异大于第三预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
关于S710~S730的具体内容,可以参考上述实施例中的描述,为避免重复,此处不再赘述。
本实施例综合考虑距离差异和第一图像差异的权重,使连续性分析结果更加全面、准确。
图8为本申请另一实施例提供的分析连续性得到连续性分析结果的流程示意图。
S810:根据cubic算法确定拼接张与第一影像序列中与拼接张相邻的至少一个第一影像的第二图像差异;
S820:根据cubic算法确定拼接张与第二影像序列中与拼接张相邻的至少一个第二影像的第三图像差异;
S830:根据第二图像差异和第三图像差异,确定第二连续性差异;
S840:若第二连续性差异大于预设第四阈值,确定第一影像序列和第二影像序列之间不具备连续性。
关于S810~S840的具体内容,可以参考上述实施例中的描述,为避免重复,此处不再赘述。
本实施例将同床位相邻影像之间的连续性差异和不同床位相邻影像之间的连续性差异的差异作为分析影像连续性的因素,使连续性分析结果更加全面、准确。
图10为本申请一实施例提供的确定医学影像拼接异常的装置的结构示意图,包括:
获取模块1010,用于获取在第一分段床位扫描的第一影像序列,并获取在第二分段床位扫描的第二影像序列,其中第一影像序列与第二影像序列为待拼接的影像序列;
第一确定模块1020,用于确定拼接张,拼接张为第二影像序列中的影像,拼接张用于拼接第一影像序列和第二影像序列;
分析模块1030,用于分析第一影像序列中的距离拼接张最近的至少一个第一影像和第二影像序列中的距离拼接张最近的至少一个第二影像的连续性,得到连续性分析结果;
第二确定模块1040,用于根据连续性分析结果,确定第一影像序列和第二影像序列是否拼接异常。
本申请的确定医学影像拼接异常的装置通过判断多床位拼接影像的连续性,实现了对医学影像拼接异常的自动提醒,帮助医生和操作人员快速确认多床位影像拼接的正确性,提高影像处理的效率和质量。
根据本申请的实施例,第一确定模块1020包括:若第一影像序列与第二影像序列无交叠,选取第二影像序列中距离第一影像序列最近的第二影像作为拼接张;若第一影像序列与第二影像序列有交叠,选取第二影像序列中与第一影像序列无交叠的影像中的距离第一影像序列最近的第二影像作为拼接张。
根据本申请的实施例,分析模块1030根据拼接张到第一影像序列中与拼接张相邻的第一影像的第一距离,以及拼接张到第二影像序列中与拼接张相邻的第二影像的第二距离,确定拼接张的距离差异;若距离差异大于第一预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
根据本申请的实施例,分析模块1030采用插值法针对拼接张所处的位置进行拟合,得到拟合拼接张;确定拼接张和拟合拼接张的第一图像差异;若第一图像差异大于第二预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
根据本申请的实施例,分析模块1030采用插值法针对拼接张所处的位置进行拟合,得到拟合拼接张,包括:根据cubic算法将与拼接张相邻的n张第一影像和与拼接张相邻的n张第二影像进行拟合,得到拟合拼接张,其中n为大于0的整数。
根据本申请的实施例,分析模块1030根据拼接张到第一影像序列中与拼接张相邻的第一影像的第一距离,以及拼接张到第二影像序列中与拼接张相邻的第二影像的第二距离,确定拼接张的距离差异;将距离差异和第一图像差异进行加权,得到第一连续性差异;若第一连续性差异大于第三预设阈值,确定第一影像序列和第二影像序列之间不具备连续性。
根据本申请的实施例,分析模块1030根据cubic算法确定拼接张与第一影像序列中与拼接张相邻的至少一个第一影像的第二图像差异;根据cubic算法确定拼接张与第二影像序列中与拼接张相邻的至少一个第二影像的第三图像差异;根据第二图像差异和第三图像差异,确定第二连续性差异;若第二连续性差异大于预设第四阈值,确定第一影像序列和第二影像序列之间不具备连续性。
关于确定医学影像拼接异常的装置的具体限定可以参见上文中关于确定医学影像拼接异常的方法的限定,在此不再赘述。
图11所示为本申请一示例性实施例提供的用于执行确定医学影像拼接异常的方法的电子设备的框图,包括处理器1110和存储器1120。
存储器1120用于存储所述处理器可执行指令。处理器用于运行可执行指令以执行上述实施例中任一项所述的确定医学影像拼接异常的方法。
本申请还提供一种计算机可读存储介质,存储介质存储有计算机程序,计算机程序用于执行上述实施例中任一项所述的确定医学影像拼接异常的方法。
本申请的确定医学影像拼接异常的方法和装置,通过判断多床位拼接影像的连续性,实现了对医学影像拼接异常的自动提醒,帮助医生和操作人员快速确认多床位影像拼接的正确性,提高影像处理的效率和质量。
在上述实施例中,可以全部或部分地通过软件、硬件、固件或者其他任意组合来实现。当使用软件实现时,可以全部或部分地以计算机程序产品的形式实现。所述计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行所述计算机程序指令时,全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中,或者从一个计算机可读存储介质向另一个计算机可读存储介质传输,例如,所述计算机指令可以从一个网站站点、计算机、服务器端或数据中心通过有线(例如同轴电缆、光纤、数字用户线(digital subscriber line,DSL))或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器端或数据中心进行传输。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器端、数据中心等数据存储设备。所述可用介质可以是磁性介质(例如,软盘、硬盘、磁带)、光介质(例如数字视频光盘(digital video disc,DVD))、或者半导体介质(例如固态硬盘(solid statedisk,SSD))等。
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本申请的范围。
在本申请所提供的几个实施例中,应该理解到,所揭露的系统、装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。

Claims (9)

1.一种确定医学影像拼接异常的方法,其特征在于,包括:
获取在第一分段床位扫描的第一影像序列,并获取在第二分段床位扫描的第二影像序列,其中所述第一影像序列与所述第二影像序列为待拼接的影像序列;
确定拼接张,所述拼接张为所述第二影像序列中的影像,所述拼接张用于拼接所述第一影像序列和所述第二影像序列;
分析所述第一影像序列中的距离所述拼接张最近的至少一个第一影像和所述第二影像序列中的距离所述拼接张最近的至少一个第二影像的连续性,得到连续性分析结果;
根据所述连续性分析结果,确定所述第一影像序列和所述第二影像序列是否拼接异常,其中,所述确定拼接张,包括:
若所述第一影像序列与所述第二影像序列无交叠,选取所述第二影像序列中距离所述第一影像序列最近的第二影像作为所述拼接张;
若所述第一影像序列与所述第二影像序列有交叠,选取所述第二影像序列中与所述第一影像序列无交叠的影像中的距离所述第一影像序列最近的第二影像作为所述拼接张。
2.根据权利要求1所述的方法,其特征在于,所述分析所述第一影像序列中的距离所述拼接张最近的至少一个第一影像和所述第二影像序列中的距离所述拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,包括:
根据所述拼接张到所述第一影像序列中与所述拼接张相邻的第一影像的第一距离,以及所述拼接张到所述第二影像序列中与所述拼接张相邻的第二影像的第二距离,确定所述拼接张的距离差异;
若所述距离差异大于第一预设阈值,确定所述第一影像序列和所述第二影像序列之间不具备连续性。
3.根据权利要求1所述的方法,其特征在于,所述分析所述第一影像序列中的距离所述拼接张最近的至少一个第一影像和所述第二影像序列中的距离所述拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,包括:
采用插值法针对所述拼接张所处的位置进行拟合,得到拟合拼接张;
确定所述拼接张和所述拟合拼接张的第一图像差异;
若所述第一图像差异大于第二预设阈值,确定所述第一影像序列和所述第二影像序列之间不具备连续性。
4.根据权利要求3所述的方法,其特征在于,所述采用插值法针对所述拼接张所处的位置进行拟合,得到拟合拼接张,包括:
根据cubic算法将与所述拼接张相邻的n张第一影像和与所述拼接张相邻的n张第二影像进行拟合,得到拟合拼接张,其中n为大于0的整数。
5.根据权利要求3所述的方法,其特征在于,所述分析所述第一影像序列中的距离所述拼接张最近的至少一个第一影像和所述第二影像序列中的距离所述拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,还包括:
根据所述拼接张到所述第一影像序列中与所述拼接张相邻的第一影像的第一距离,以及所述拼接张到所述第二影像序列中与所述拼接张相邻的第二影像的第二距离,确定所述拼接张的距离差异;
将所述距离差异和所述第一图像差异进行加权,得到第一连续性差异;
若所述第一连续性差异大于第三预设阈值,确定所述第一影像序列和所述第二影像序列之间不具备连续性。
6.根据权利要求1所述的方法,其特征在于,所述分析所述第一影像序列中的距离所述拼接张最近的至少一个第一影像和所述第二影像序列中的距离所述拼接张最近的至少一个第二影像的连续性,得到连续性分析结果,包括:
根据cubic算法确定所述拼接张与所述第一影像序列中与所述拼接张相邻的至少一个第一影像的第二图像差异;
根据cubic算法确定所述拼接张与所述第二影像序列中与所述拼接张相邻的至少一个第二影像的第三图像差异;
根据所述第二图像差异和所述第三图像差异,确定第二连续性差异;
若所述第二连续性差异大于预设第四阈值,确定所述第一影像序列和所述第二影像序列之间不具备连续性。
7.一种确定医学影像拼接异常的装置,其特征在于,包括:
获取模块,用于获取在第一分段床位扫描的第一影像序列,并获取在第二分段床位扫描的第二影像序列,其中所述第一影像序列与所述第二影像序列为待拼接的影像序列;
第一确定模块,用于确定拼接张,所述拼接张为所述第二影像序列中的影像,所述拼接张用于拼接所述第一影像序列和所述第二影像序列;
分析模块,用于分析所述第一影像序列中的距离所述拼接张最近的至少一个第一影像和所述第二影像序列中的距离所述拼接张最近的至少一个第二影像的连续性,得到连续性分析结果;
第二确定模块,用于根据所述连续性分析结果,确定所述第一影像序列和所述第二影像序列是否拼接异常,其中,所述第一确定模块包括:
若所述第一影像序列与所述第二影像序列无交叠,选取所述第二影像序列中距离所述第一影像序列最近的第二影像作为所述拼接张;
若所述第一影像序列与所述第二影像序列有交叠,选取所述第二影像序列中与所述第一影像序列无交叠的影像中的距离所述第一影像序列最近的第二影像作为所述拼接张。
8.一种计算机可读存储介质,所述存储介质存储有计算机程序,所述计算机程序用于执行上述权利要求1至6中任一项所述的确定医学影像拼接异常的方法。
9.一种电子设备,包括:
处理器;
用于存储所述处理器可执行指令的存储器,
其中,所述处理器用于执行上述权利要求1至6中任一项所述的确定医学影像拼接异常的方法。
CN202111083284.1A 2021-09-16 2021-09-16 确定医学影像拼接异常的方法及装置 Active CN113538251B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111083284.1A CN113538251B (zh) 2021-09-16 2021-09-16 确定医学影像拼接异常的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111083284.1A CN113538251B (zh) 2021-09-16 2021-09-16 确定医学影像拼接异常的方法及装置

Publications (2)

Publication Number Publication Date
CN113538251A CN113538251A (zh) 2021-10-22
CN113538251B true CN113538251B (zh) 2021-12-28

Family

ID=78123179

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111083284.1A Active CN113538251B (zh) 2021-09-16 2021-09-16 确定医学影像拼接异常的方法及装置

Country Status (1)

Country Link
CN (1) CN113538251B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111311655A (zh) * 2020-02-13 2020-06-19 东软医疗系统股份有限公司 多模态图像配准方法、装置、电子设备、存储介质
CN111784788A (zh) * 2020-06-04 2020-10-16 深圳深透医疗科技有限公司 一种基于深度学习的pet快速成像方法和系统
CN112907540A (zh) * 2021-02-22 2021-06-04 浙江大华技术股份有限公司 一种拼接异常检测方法、装置、设备及介质
CN113344004A (zh) * 2021-03-26 2021-09-03 慧影医疗科技(北京)有限公司 图像特征生成方法、图像识别方法和装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5380348B2 (ja) * 2010-03-31 2014-01-08 富士フイルム株式会社 内視鏡観察を支援するシステムおよび方法、並びに、装置およびプログラム
CN107067366A (zh) * 2017-04-17 2017-08-18 苏州海斯菲德信息科技有限公司 小动物医学影像多段配准及拼接方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111311655A (zh) * 2020-02-13 2020-06-19 东软医疗系统股份有限公司 多模态图像配准方法、装置、电子设备、存储介质
CN111784788A (zh) * 2020-06-04 2020-10-16 深圳深透医疗科技有限公司 一种基于深度学习的pet快速成像方法和系统
CN112907540A (zh) * 2021-02-22 2021-06-04 浙江大华技术股份有限公司 一种拼接异常检测方法、装置、设备及介质
CN113344004A (zh) * 2021-03-26 2021-09-03 慧影医疗科技(北京)有限公司 图像特征生成方法、图像识别方法和装置

Also Published As

Publication number Publication date
CN113538251A (zh) 2021-10-22

Similar Documents

Publication Publication Date Title
US6549214B1 (en) Automatic scaling of display for image management system
US8041093B2 (en) System and method for definition of DICOM header values
US20140104311A1 (en) Medical image display method using virtual patient model and apparatus thereof
US20100135552A1 (en) Medical Imaging with Accessible Computer Assisted Detection
US10825149B2 (en) Defective pixel correction using adversarial networks
US10600136B2 (en) Identification of medical concepts for imaging protocol selection
US11257211B2 (en) Medical image processing apparatus, medical image processing system, and medical image processing method
CN112150376B (zh) 血管医学影像分析方法、装置、计算机设备和存储介质
US11232604B2 (en) Device for marking image data
US6127669A (en) Computer-aided determination of window and level settings for filmless radiology
US20080292162A1 (en) Method and System for Generating a Collage to Summarize a Medical Dataset
CN113658175B (zh) 一种征象数据的确定方法及装置
US20090080692A1 (en) Size measurement apparatus, image display unit, computer-readable size measurement program storage medium, and computer-readable image display program storage medium
WO2021122342A1 (en) Systems and methods for immediate image quality feedback
CN114864035A (zh) 影像报告生成方法、装置、系统、设备及存储介质
CN113538251B (zh) 确定医学影像拼接异常的方法及装置
WO2023032436A1 (ja) 医用画像処理装置、医用画像処理方法及びプログラム
EP4318488A1 (en) Magnetic resonance image processing apparatus to which slice resolution enhancement is applied, and method therefor
CN116869555A (zh) 扫描协议调节方法、装置以及存储介质
EP3949909A2 (en) Device and application for integrated management of artificial joints
EP4195669A1 (en) Image quality management system
US9734286B2 (en) System and method for splitting DICOM medical image series into framesets
CN113344920B (zh) 一种纠正医学影像摆位的方法和装置
US20060115136A1 (en) Image processing method and apparatus
US20230245760A1 (en) Display apparatus, image management server, and recording medium

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