CN110265120B - 医学图像处理方法、装置、计算机设备和存储介质 - Google Patents

医学图像处理方法、装置、计算机设备和存储介质 Download PDF

Info

Publication number
CN110265120B
CN110265120B CN201910383106.7A CN201910383106A CN110265120B CN 110265120 B CN110265120 B CN 110265120B CN 201910383106 A CN201910383106 A CN 201910383106A CN 110265120 B CN110265120 B CN 110265120B
Authority
CN
China
Prior art keywords
image
data
scanning
processed
distribution information
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
CN201910383106.7A
Other languages
English (en)
Other versions
CN110265120A (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.)
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 CN201910383106.7A priority Critical patent/CN110265120B/zh
Publication of CN110265120A publication Critical patent/CN110265120A/zh
Application granted granted Critical
Publication of CN110265120B publication Critical patent/CN110265120B/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
    • G06T5/00Image enhancement or restoration
    • G06T5/80Geometric correction
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10104Positron emission tomography [PET]

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Epidemiology (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本申请涉及一种医学图像处理方法、装置、计算机设备和存储介质。所述方法包括:获取预扫描数据,基于预扫描数据得到预扫描数据对应预扫描图像的亮度分布信息;获取正式扫描图像数据;根据亮度分布信息处理正式扫描图像数据,得到待处理图像数据;将待处理图像数据进行回顾式处理;重建回顾式处理后的待处理图像数据,得到目标图像。采用本方法能够系统地处理医学图像中亮度不均匀的分布。

Description

医学图像处理方法、装置、计算机设备和存储介质
技术领域
本申请涉及图像处理技术领域,特别是涉及一种医学图像处理方法、装置、 计算机设备和存储介质。
背景技术
医学图像包括CT(Computed Tomography,电子计算机断层)图像、 PET(PositronEmission Computed Tomography,正电子发射型计算机断层)图 像以及磁共振图像等。在磁共振成像中,一般都使用表面线圈实现特定部位的 扫描,但是由于表面线圈存在信号接收的敏感度分布,导致最终重建出的图像 在亮度分布上会有不均匀的现象,并且该问题随着静磁场的增加会变得更加严 重,甚至影响到正常的临床诊断。目前针对该问题的处理方案有很多种,按照 是否依赖系统参数通常可以分为回顾式方案和前瞻式方案两大类。其中,回顾 式方案是一种主要基于图像处理的方案,例如基于同态滤波技术、基于表面拟合、图像分割技术以及直方图分析等。而前瞻式方案则是一种依赖系统成像参 数的方案,例如依赖于表面线圈的信号接收敏感性分布、成像序列等。
然而,回顾式方案需要一个假设的前提条件,即图像的亮度不均匀性是低 频缓变的,因此回顾式方案无法处理存在高频不均匀亮度信号的图像,导致使 用范围受限。而前瞻式方案虽然能够充分利用系统参数,且无论对于高频的还 是低频的亮度不均匀都能很好的处理。但是,也正是由于其对系统参数的依赖, 使得校正图像亮度不均匀的质量取决于预采集的数据是否是合理的亮度分布。 并且,随着静磁场的增加,射频磁场的发射场和接收场会变得越来越不均匀, 因此图像亮度分布会随着静磁场的增加变得越来越不均匀。而经过理论分析和 实验证明,前瞻式方案还严重依赖于静磁场、射频磁场的发射场和接收场,从 而导致前瞻式方案无法处理高场强下的亮度不均匀图像。
发明内容
基于此,有必要针对上述技术问题,提供一种能够系统地处理医学图像亮 度分布不均匀问题的医学图像处理方法、装置、计算机设备和存储介质。
一种医学图像处理方法,所述方法包括:
获取预扫描数据,基于所述预扫描数据得到所述预扫描数据对应预扫描图 像的亮度分布信息;
获取正式扫描图像数据;
根据所述亮度分布信息处理所述正式扫描图像数据,得到待处理图像数据;
将所述待处理图像数据进行回顾式处理;
重建回顾式处理后的待处理图像数据,得到目标图像。
在其中一个实施例中,所述方法还包括:
评估确定所述目标图像是否存在残留的不均匀场;
若不存在残留的不均匀场,则输出所述目标图像;
若存在残留的不均匀场,则获取所述残留不均匀场,根据所述残留不均匀 场校正所述亮度分布信息,根据校正后的所述亮度分布信息处理下一个正式扫 描图像数据。
在其中一个实施例中,所述获取预扫描数据,基于所述预扫描数据得到亮 度分布信息之后,根据所述亮度分布信息处理所述正式扫描图像数据之前,还 包括:
将所述亮度分布信息进行回顾式处理,得到校正后的亮度分布信息。
在其中一个实施例中,所述获取预扫描数据,基于所述预扫描数据得到所 述预扫描数据对应预扫描图像的亮度分布信息的步骤,包括:
接收医学扫描设备发送的预扫描数据,所述预扫描数据包括容积线圈对应 的第一预扫描数据和表面线圈对应的第二预扫描数据;
基于所述第一预扫描数据和所述第二预扫描数据,得到所述第一预扫描数 据和所述第二预扫描数据对应预扫描图像的亮度分布信息。
在其中一个实施例中,所述基于所述第一预扫描数据和所述第二预扫描数 据,得到所述第一预扫描数据和所述第二预扫描数据对应预扫描图像的亮度分 布信息的步骤,包括:
重建所述第一预扫描数据和第二预扫描数据,分别得到第一预扫描图像和 第二预扫描图像;
计算所述第一预扫描图像和第二预扫描图像的商,得到亮度分布信息。
在其中一个实施例中,将所述亮度分布信息与所述正式扫描图像数据进行 点乘计算,得到待处理图像数据。
在其中一个实施例中,所述将所述待处理图像数据进行回顾式处理,并重 建回顾式处理后的待处理图像数据,得到目标图像的步骤,包括:
获取不均匀乘性场和图像噪声;
基于所述不均匀乘性场和所述图像噪声处理所述待处理图像数据,并重建 处理后的待处理图像数据,得到目标图像。
一种医学图像处理装置,其特征在于,所述装置包括:
重建模块,用于获取预扫描数据,基于所述预扫描数据得到所述预扫描数 据对应预扫描图像的亮度分布信息;
获取模块,用于获取正式扫描图像数据;
亮度处理模块,用于根据所述亮度分布信息处理所述正式扫描图像数据, 得到待处理图像数据;
回顾式处理模块,将所述待处理图像数据进行回顾式处理;
所述重建模块,还用于重建回顾式处理后的待处理图像数据,得到目标图 像。
一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序, 其特征在于,所述处理器执行所述计算机程序时实现上述任一项所述医学图像 处理方法的步骤。
一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计 算机程序被处理器执行时实现上述任一项所述医学图像处理方法的步骤。
上述医学图像处理方法、装置、计算机设备和存储介质,获取预扫描数据, 基于得到的预扫描数据得到与线圈相关的亮度分布信息,根据亮度分布信息处 理正式扫描获取到正式扫描图像数据,即存在亮度不均匀的图像数据,从而得 到待处理图像数据。以此保证亮度不均匀图像中大部分的低频亮度不均匀分布 和所有的高频亮度分布不均匀被校正。然后通过回顾式处理进一步处理待处理 图像数据,校正残留的低频亮度分布不均匀。上述方法不仅能够处理高频亮度 分布不均匀,还能校正低频亮度分布不均匀,从而系统地解决了图像亮度分布 不均匀的问题。
附图说明
图1为一个实施例中医学图像处理方法的应用环境图;
图2为一个实施例中医学图像处理方法的流程示意图;
图3为一个实施例中包括评估步骤的医学图像处理方法的流程示意图;
图4为另一个实施例中医学图像处理方法的流程示意图;
图5为一个实施例中医学图像处理装置的结构框图;
图6为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实 施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅 用以解释本申请,并不用于限定本申请。
本申请提供的医学图像处理方法,可以应用于如图1所示的应用环境中。 其中,医学图像处理装置102通过网络与医学扫描设备104进行通信。医学扫 描设备104根据预扫描参数进行预扫描,将预扫描得到的预扫描数据发送给医 学图像处理装置102。医学图像处理装置102获取从医学扫描设备104接收的预 扫描数据,基于预扫描数据得到亮度分布信息。医学扫描设备104根据正式扫 描参数进行正式扫描,将正式扫描得到的正式扫描图像数据发送给医学图像处 理装置102。医学图像处理装置102获取从医学扫描设备104处接收的正式扫描 图像数据,根据亮度分布信息处理正式扫描图像数据,得到待处理图像。医学图像处理装置102将待处理图像数据进行回顾式处理,并重建回顾式处理后的 待处理图像数据,得到目标图像。其中,医学图像处理装置102可以但不限于 是服务器、各种个人计算机、笔记本电脑、智能手机、平板电脑和便携式可穿 戴设备,医学扫描设备104可以为磁共振扫描设备。
在一个实施例中,如图2所示,提供了一种医学图像处理方法,以该方法 应用于图1中的医学图像处理装置102为例进行说明,包括以下步骤:
步骤S202,获取预扫描数据,基于预扫描数据得到预扫描数据对应预扫描 图像的亮度分布信息。
其中,预扫描是指在正式扫描之前预先进行的扫描,预扫描数据则是预扫 描得到的图像数据。即,在正式扫描前对被扫描物体进行快速扫描采集一些图 像数据,例如图像轮廓信息、扫描定位像等低分辨率的医学图像。亮度分布信 息即是预扫描数据对应的预扫描图像的图像亮度分布信息。
具体地,当医学扫描设备接收到预扫描指令以及预扫描参数时,根据预扫 描参数对被扫描物体进行预扫描。例如,用户通过与医学扫描设备连接的医学 图像处理装置的输入装置或者其他医学扫描设备连接的计算机设备,下发预扫 描指令以及手动输入或者通过点击选择对应的预扫描参数。然后,医学扫描设 备将预扫描得到的预扫描数据发送至连接的医学图像处理装置,医学图像处理 装置获取到从医学扫描设备接收过来的预扫描数据后,将预扫描数据进行重建, 得到预扫描数据对应的预扫描图像,从预扫描图像中提取出亮度分布信息,该 亮度分布信息与医学扫描设备的系统参数相关。其中,医学扫描设备可以理解 为扫描器、例如MRI扫描器(磁共振扫描器)等。可选地,医学扫描设备的系 统参数可包括主磁场的场强、主磁场的均匀性、梯度场的均匀性、射频场的场 强、射频接收线圈的灵敏度等中的一种或多种的组合。
步骤S204,获取正式扫描图像数据。
其中,正式扫描即为对被扫描物体进行正式的扫描,一般被扫描物体进行 扫描时,为了得到合乎标准的医学图像,扫描时所用的扫描参数应当符合标准。 正式扫描即是根据正式扫描参数,扫描得到符合标准医学图像的图像数据。也 就是说,正式扫描和预扫描一样有对应的正式扫描参数,正式扫描与预扫描使 用不同的扫描参数。例如为了得到高质量医学图像,正式扫描的时间可以长于 预扫描的时间等。而正式扫描图像数据就是正式扫描得到的图像数据,在本申 请中也可称之为“成像数据”。
具体地,当进行预扫描得到预扫描数据,并且基于预扫描数据得到亮度分 布信息后,即可进入正式扫描。医学扫描设备同样接收到正式扫描指令以及对 应的正式扫描参数后,根据正式扫描参数进行正式扫描,并将得到的正式扫描 参数发送至医学图像处理装置,即医学图像处理装置可获取从医学扫描设备处 接收到的正式扫描图像数据。
步骤S208,根据亮度分布信息处理正式扫描图像数据,得到待处理图像数 据。
其中,待处理图像数据是指经过亮度分布信息处理,已经消除大部分低频 亮度不均匀分布和所有高频亮度分布不均匀后的正式扫描图像数据。
具体地,通过将亮度分布信息与正式扫描图像数据进行点乘计算得到的数 据即为待处理图像数据。点乘计算即是将亮度分布信息与正式扫描图像进行对 应相乘。以MRI为例,由于MRI进行扫描时依靠MRI扫描器中的线圈发射和接 收信号,以及由于线圈存在信号接收的敏感度分布,最终导致最终重建的图像 在亮度上会有不均匀分布的现象,因此,首先进行预扫描,根据预扫描数据得 到能够表达线圈信号亮度分布的亮度分布信息。根据该亮度分布信息处理正式 扫描图像数据,能够消除正式扫描图像数据中大部分低频亮度不均匀分布和所 有高频亮度分布不均匀。也就是说,将经过亮度分布信息处理后得到的待处理 图像数据进行图像重建得到的图像,相比于直接将正式扫描图像数据重建得到 图像来说,该图像中没有高频亮度不均匀分布以及没有大部分低频亮度不均匀, 只存在小部分亮度分布信息无法消除的低频亮度不均匀。
步骤S210,将待处理图像数据进行回顾式处理。
回顾式处理是一种消除医学图像中低频亮度不均匀的处理方案。由于图像 存在亮度分布不均匀除了扫描线圈造成的之外,还有可能是因为外部干扰而导 致的。然而经过亮度分布信息处理只能够消除扫描线圈引起的亮度不均匀,其 他外部干扰所导致的亮度不均匀无法通过亮度分布信息进行消除。因此,将经 过亮度分布信息处理的正式扫描图像数据,即待处理图像数据再次进行回顾式 处理,能够消除亮度分布信息无法消除的低频亮度不均匀分布。
步骤S212,重建回顾式处理后的待处理图像数据,得到目标图像。
其中,重建是指图像重建,即将图像数据经过数字处理获取三维物体的形 状信息的技术。也就是将待处理图像数据进行图像重建后,能够得到被扫描物 体对应的医学图像。目标图像即为经过处理后得到的亮度均匀的图像。具体地, 将经过回顾式处理的待处理图像数据进行图像重建,得到正式扫描对应的医学 图像,而该医学图像为目标图像,即亮度分布均匀的图像。
上述医学图像处理方法,获取预扫描数据,基于得到的预扫描数据得到与 线圈相关的亮度分布信息,根据亮度分布信息处理正式扫描获取到正式扫描图 像数据,即存在亮度不均匀的图像数据,从而得到待处理图像数据。以此保证 亮度不均匀图像中大部分的低频亮度不均匀分布和所有的高频亮度分布不均匀 被校正。然后通过回顾式处理进一步处理待处理图像数据,校正残留的低频亮 度分布不均匀。上述方法不仅能够处理高频亮度分布不均匀,还能校正低频亮 度分布不均匀,从而系统地解决了图像亮度分布不均匀的问题。
在一个实施例中,将待处理图像数据进行回顾式处理,并重建回顾式处理 后的待处理图像数据,得到目标图像具体包括:获取不均匀乘性场和图像噪声。 基于不均匀乘性场和图像噪声处理待处理图像数据,并重建处理后的待处理图 像数据,得到目标图像。
其中,图像噪声是图像数据中不必要或多余的干扰信息。不均匀乘性场为 预设的低频的亮度不均匀分布。具体地,回顾式处理一般是在原有的均匀图像 上叠加一个低频的不均匀乘性场,计算公式如下:
v(x)=u(x)f(x)+n(x)
其中,v(x)为实际扫描得到的不均匀的图像,u(x)是理想的亮度分布均匀 的图像,f(x)是叠加的不均匀乘性场,n(x)是图像噪声。在本实施例中,需要 将实际扫描得到的不均匀的图像进行处理后得到均匀的图像,即计算公式如下:
Figure BDA0002053979780000071
可以理解为,v(x)为正式扫描图像数据经过亮度分布信息处理后的待处理 图像数据。获取预设的不均匀乘性场f(x),并且从待处理图像数据中获取图像 噪声n(x),将减去图像噪声后的待处理图像数据与不均匀的乘性场做除法运算, 得到目标图像。
在一个实施例中,如图3所示,在步骤S212之后,还包括以下步骤:
步骤S214,评估确定目标图像是否存在残留的不均匀场。
其中,残留的不均匀场即步骤S212得到目标图像经过评估后,确定该目标 图像还有残留的亮度不均匀分布。具体地,当得到目标图像后,通过重新计算 该目标图像的不均匀度,确定该目标图像是否还有残留的不均匀场。
具体地,目标图像的均匀度计算可通过如下方法确定:将目标图像划分为 不同的区域,在不同的区域测量信号强度值。其中,被测量的区域优选为能够 较全面的反应图像的均匀度,且未靠近图像边缘的区域。根据测量得到的信号 强度值进行计算得到均匀度的计算结果。例如,可对测量得到的信号强度值进 行统计分析得到它们的标准偏差值和均值。
在一个实施例中,目标图像均匀度计算公式为:
Figure BDA0002053979780000081
其中,UZ表示图像的均匀度,SMAX表示所测区域中信号最大值,SMIX表示 所测区域中信号最小值。
步骤S216,若不存在残留的不均匀场,则输出目标图像。
步骤S218,若存在残留的不均匀场,则获取残留不均匀场,根据残留不均 匀场校正亮度分布信息,根据校正后的亮度分布信息处理下一个正式扫描图像 数据。
具体地,当目标图像经过评估确定不存在残留的不均匀场,即表明基于预 扫描的亮度分布信息以及回顾式处理后能够完全将医学图像的亮度不均匀分布 消除,则可以将该目标图像输出显示到显示装置上,例如显示屏等,便于用户 查看。
而当目标图像经过评估确定该目标图像实际上还存在残留的不均匀场,则 获取该残留的不均匀场的数据。并将评估所得到残留的不均匀场的数据对预扫 描获取到的亮度分布信息进行校正,将校正后的亮度分布信息处理正式扫描图 像数据。可以理解为,当根据残留的不均匀场对亮度分布信息进行校正后,可 以利用该校正后的亮度分布信息重新处理医学扫描设备发送的正式扫描图像数 据,也可以是利用该校正后的亮度分布信息处理下一个正式扫描图像数据。其 中,根据残留的不均匀场对亮度分布信息进行校正,即为将残留的不均匀场与 亮度分布信息进行叠加操作,即残留的不均匀场与亮度分布信息进行点乘计算。 在本实施中,利用校正后的亮度分布信息处理正式扫描图像数据,减少最终得 到的目标图像还存在残留的不均匀场的概率。
在一个实施例中,如图4所示,提供另一种医学图像处理方法,即在步骤 S202之后,步骤S208之前,还包括步骤S206,将亮度分布信息进行回顾式处 理,得到校正后的亮度分布信息。
具体地,当基于预扫描数据获取到亮度分布信息之后,通过对亮度分布信 息进行回顾式处理,从而得到校正后的亮度分布信息。也就是说,通过直接校 正预扫描得到的亮度分布信息来达到最终校正效果。可以理解为,当获取到预 扫描数据后,基于预扫描数据得到预扫描数据对应预扫描图像的亮度分布信息 后,首先通过对亮度分布信息进行回顾式处理,得到校正后的亮度分布信息。 根据校正后的亮度分布信息处理正式扫描获取的正式扫描图像数据,得到待处 理图像数据。再次对待处理图像数据进行回顾式处理,重建进行回顾式处理后 的待处理图像数据得到目标图像。或者,不再对待处理图像数据进行回顾式处 理,而直接重建待处理图像数据得到对应的医学图像作为目标图像。也就是说, 当亮度分布信息进行过回顾式处理后,后续得到的待处理图像数据可以无需进 行回顾式处理直接进行图像重建,也可以再次进行回顾式处理之后在进行图像 重建。
亮度分布信息进行回顾式处理的计算公式如下:
Figure BDA0002053979780000091
其中,q1(x)为预扫描图像的亮度分布信息、f(x)为不均匀乘性场、n(x)为 图像噪声,q2(x)为校正后的亮度分布信息。
在一个实施例中,获取预扫描数据,基于预扫描数据得到预扫描数据对应 预扫描图像的亮度分布信息的步骤,具体包括:接收医学扫描设备发送的预扫 描数据,预扫描数据包括容积线圈对应的第一预扫描数据和表面线圈对应的第 二预扫描数据;基于第一预扫描数据和所述第二预扫描数据,得到第一预扫描 数据和第二预扫描数据对应预扫描图像的亮度分布信息。
具体地,进行预扫描时,使用容积线圈和表面线圈进行扫描。即向医学扫 描设备下发预扫描指令时,同时在各线圈中选定容积线圈和表面线圈。当医学 扫描设备进行扫描时,利用容积线圈和表面线圈进行扫描,并将容器线圈和表 面线圈扫描得到的预扫描数据发送至医学图像处理装置。医学图像处理装置接 收到预扫描数据则包括容积线圈对应的第一预扫描数据和表面线圈对应的第二 预扫描数据,得到第一预扫描数据和第二预扫描数据对应预扫描图像的亮度分 布信息。
在一个实施例中,基于第一预扫描数据和第二预扫描数据,得到第一预扫 描数据和第二预扫描数据对应预扫描图像的亮度分布信息的步骤,具体包括: 重建第一预扫描数据和第二预扫描数据,分别得到第一预扫描图像和第二预扫 描图像。计算第一预扫描图像和第二预扫描图像的商,得到亮度分布信息。
具体地,当得到第一预扫描数据和第二预扫描数据后,分别对第一预扫描 数据和第二预扫描数据进行图像重建,得到第一预扫描数据对应的第一预扫描 图像、第二预扫描数据对应的第二预扫描图像。然后,基于第一预扫描图像和 第二预扫描图像得到亮度分布信息。即,将第一预扫描图像与第二预扫描图像 进行除法运算所得到的图像即为亮度分布信息。由于医学图像由像素组成,即 基于第一预扫描图像中的像素和第二预扫描图像中的像素进行除法运算,所得 的像素商组成的图像即为亮度分布信息。
应该理解的是,虽然图2-4的流程图中的各个步骤按照箭头的指示依次显 示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明 确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺 序执行。而且,图2-4中的至少一部分步骤可以包括多个子步骤或者多个阶段, 这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻 执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它 步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在一个实施例中,如图5所示,提供了一种医学图像处理装置,包括:重 建模块502、获取模块504、亮度处理模块506和回顾式处理模块508,其中:
重建模块502,用于获取预扫描数据,基于预扫描数据得到预扫描数据对应 预扫描图像的亮度分布信息。
获取模块504,用于获取正式扫描图像数据。
亮度处理模块506,用于根据亮度分布信息处理正式扫描图像数据,得到待 处理图像数据。
回顾式处理模块508,用于将待处理图像数据进行回顾式处理。
重建模块502还用于重建回顾式处理后的待处理图像数据,得到目标图像。
在一个实施例中,医学图像处理装置还包括评估模块510,用于评估确定目 标图像是否存在残留的不均匀场;若不存在残留的不均匀场,则输出目标图像; 若存在残留的不均匀场,则获取残留不均匀场,根据残留不均匀场校正亮度分 布信息,根据校正后的亮度分布信息处理下一个正式扫描图像数据。
在一个实施例中,回顾式处理模块508还用于将亮度分布信息进行回顾式 处理,得到校正后的亮度分布信息。
在一个实施例中,重建模块502还用于接收医学扫描设备发送的预扫描数 据,预扫描数据包括容积线圈对应的第一预扫描数据和表面线圈对应的第二预 扫描数据;基于第一预扫描数据和所述第二预扫描数据,得到第一预扫描数据 和第二预扫描数据对应预扫描图像的亮度分布信息。
在一个实施例中,重建模块502还用于重建第一预扫描数据和第二预扫描 数据,分别得到第一预扫描图像和第二预扫描图像。计算第一预扫描图像和第 二预扫描图像的商,得到亮度分布信息。
在一个实施例中,亮度处理模块506还用于将亮度分布信息与正式扫描图 像数据进行点乘计算,得到待处理图像数据。
在一个实施例中,回顾式处理模块508还用于获取不均匀乘性场和图像噪 声。基于不均匀乘性场和图像噪声处理待处理图像数据,并重建处理后的待处 理图像数据,得到目标图像。
关于医学图像处理装置的具体限定可以参见上文中对于医学图像处理方法 的限定,在此不再赘述。上述医学图像处理装置中的各个模块可全部或部分通 过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算 机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以 便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,该计算机设备可以是服务器, 其内部结构图可以如图6所示。该计算机设备包括通过系统总线连接的处理器、 存储器、网络接口、数据库、显示屏和输入装置。其中,该计算机设备的处理 器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、 内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内 存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计 算机设备的数据库用于存储处理数据。该计算机设备的网络接口用于与外部的 终端通过网络连接通信。该计算机程序被处理器执行时以实现一种医学图像处 理方法。
该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏,该计算机 设备的输入装置可以是显示屏上覆盖的触摸层,也可以是计算机设备外壳上设 置的按键、轨迹球或触控板,还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解,图6中示出的结构,仅仅是与本申请方案相关 的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定, 具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件, 或者具有不同的部件布置。
在一个实施例中,提供了一种计算机设备,包括存储器和处理器,存储器 中存储有计算机程序,该处理器执行计算机程序时实现以下步骤:
获取预扫描数据,基于预扫描数据得到预扫描数据对应预扫描图像的亮度 分布信息;
获取正式扫描图像数据;
根据亮度分布信息处理正式扫描图像数据,得到待处理图像数据;
将待处理图像数据进行回顾式处理;
重建回顾式处理后的待处理图像数据,得到目标图像。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:
评估确定目标图像是否存在残留的不均匀场;若不存在残留的不均匀场, 则输出目标图像;若存在残留的不均匀场,则获取残留不均匀场,根据残留不 均匀场校正亮度分布信息,根据校正后的亮度分布信息处理下一个正式扫描图 像数据。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:
将亮度分布信息进行回顾式处理,得到校正后的亮度分布信息。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:
接收医学扫描设备发送的预扫描数据,预扫描数据包括容积线圈对应的第 一预扫描数据和表面线圈对应的第二预扫描数据;基于第一预扫描数据和所述 第二预扫描数据,得到第一预扫描数据和第二预扫描数据对应预扫描图像的亮 度分布信息。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:
重建第一预扫描数据和第二预扫描数据,分别得到第一预扫描图像和第二 预扫描图像。计算第一预扫描图像和第二预扫描图像的商,得到亮度分布信息。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:
将亮度分布信息与正式扫描图像数据进行点乘计算,得到待处理图像数据。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:获取不均匀 乘性场和图像噪声。基于不均匀乘性场和图像噪声处理待处理图像数据,并重 建处理后的待处理图像数据,得到目标图像。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程 序,计算机程序被处理器执行时实现以下步骤:
获取预扫描数据,基于预扫描数据得到预扫描数据对应预扫描图像的亮度 分布信息;
获取正式扫描图像数据;
根据亮度分布信息处理正式扫描图像数据,得到待处理图像数据;
将待处理图像数据进行回顾式处理;
重建回顾式处理后的待处理图像数据,得到目标图像。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:
评估确定目标图像是否存在残留的不均匀场;若不存在残留的不均匀场, 则输出目标图像;若存在残留的不均匀场,则获取残留不均匀场,根据残留不 均匀场校正亮度分布信息,根据校正后的亮度分布信息处理下一个正式扫描图 像数据。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:
将亮度分布信息进行回顾式处理,得到校正后的亮度分布信息。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:
接收医学扫描设备发送的预扫描数据,预扫描数据包括容积线圈对应的第 一预扫描数据和表面线圈对应的第二预扫描数据;基于第一预扫描数据和所述 第二预扫描数据,得到第一预扫描数据和第二预扫描数据对应预扫描图像的亮 度分布信息。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:
重建第一预扫描数据和第二预扫描数据,分别得到第一预扫描图像和第二 预扫描图像。计算第一预扫描图像和第二预扫描图像的商,得到亮度分布信息。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:
将亮度分布信息与正式扫描图像数据进行点乘计算,得到待处理图像数据。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:
获取不均匀乘性场和图像噪声。基于不均匀乘性场和图像噪声处理待处理 图像数据,并重建处理后的待处理图像数据,得到目标图像。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程, 是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于 一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述 各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、 存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。 非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM (EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存 取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形 式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据 率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM (SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM (DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述 实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特 征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细, 但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的 普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改 进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权 利要求为准。

Claims (8)

1.一种医学图像处理方法,所述方法包括:
获取预扫描数据,基于所述预扫描数据得到所述预扫描数据对应预扫描图像的亮度分布信息;
获取正式扫描图像数据;
根据所述亮度分布信息处理所述正式扫描图像数据,得到待处理图像数据,其中,将所述亮度分布信息与所述正式扫描图像数据进行对应相乘,得到待处理图像数据,所述亮度分布信息用于消除所述正式扫描图像数据中扫描线圈引起的亮度不均匀;
将所述待处理图像数据进行回顾式处理,其中,所述回顾式处理用于消除所述待处理图像数据中外部干扰导致的亮度不均匀;
重建回顾式处理后的待处理图像数据,得到目标图像;
所述将所述待处理图像数据进行回顾式处理,重建回顾式处理后的待处理图像数据,得到目标图像的步骤,包括:
获取不均匀乘性场和图像噪声;
基于所述不均匀乘性场和所述图像噪声处理所述待处理图像数据,并重建处理后的待处理图像数据,得到目标图像,计算公式如下:
Figure DEST_PATH_IMAGE001
其中,
Figure 425086DEST_PATH_IMAGE002
为所述待处理图像数据,
Figure DEST_PATH_IMAGE003
为所述不均匀乘性场,
Figure 60598DEST_PATH_IMAGE004
为所述图像噪声,
Figure DEST_PATH_IMAGE005
为所述目标图像;
所述方法还包括:
所述获取预扫描数据,基于所述预扫描数据得到所述预扫描数据对应预扫描图像的亮度分布信息之后,所述根据所述亮度分布信息处理所述正式扫描图像数据,得到待处理图像数据之前,还包括:
将所述亮度分布信息进行回顾式处理,得到校正后的亮度分布信息,计算公式如下:
Figure 389948DEST_PATH_IMAGE006
其中,
Figure DEST_PATH_IMAGE007
为所述预扫描图像的所述亮度分布信息,
Figure 277264DEST_PATH_IMAGE008
为所述校正后的亮度分布信息。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
评估确定所述目标图像是否存在残留的不均匀场;
若不存在残留的不均匀场,则输出所述目标图像;
若存在残留的不均匀场,则获取所述残留不均匀场,根据所述残留不均匀场校正所述亮度分布信息,根据校正后的所述亮度分布信息处理下一个正式扫描图像数据。
3.根据权利要求1所述的方法,其特征在于,所述获取预扫描数据,基于所述预扫描数据得到所述预扫描数据对应预扫描图像的亮度分布信息的步骤,包括:
接收医学扫描设备发送的预扫描数据,所述预扫描数据包括容积线圈对应的第一预扫描数据和表面线圈对应的第二预扫描数据;
基于所述第一预扫描数据和所述第二预扫描数据,得到所述第一预扫描数据和所述第二预扫描数据对应预扫描图像的亮度分布信息。
4.根据权利要求3所述的方法,其特征在于,所述基于所述第一预扫描数据和所述第二预扫描数据,得到所述第一预扫描数据和所述第二预扫描数据对应预扫描图像的亮度分布信息的步骤,包括:
重建所述第一预扫描数据和第二预扫描数据,分别得到第一预扫描图像和第二预扫描图像;
计算所述第一预扫描图像和第二预扫描图像的商,得到亮度分布信息。
5.根据权利要求1所述的方法,其特征在于,所述将所述亮度分布信息与所述正式扫描图像数据进行对应相乘,得到待处理图像数据的步骤,包括:
将所述亮度分布信息与所述正式扫描图像数据进行点乘计算,得到待处理图像数据。
6.一种医学图像处理装置,其特征在于,所述装置包括:
重建模块,用于获取预扫描数据,基于所述预扫描数据得到所述预扫描数据对应预扫描图像的亮度分布信息;
获取模块,用于获取正式扫描图像数据;
亮度处理模块,用于根据所述亮度分布信息处理所述正式扫描图像数据,得到待处理图像数据,其中,将所述亮度分布信息与所述正式扫描图像数据进行对应相乘,得到待处理图像数据,所述亮度分布信息用于消除所述正式扫描图像数据中扫描线圈引起的亮度不均匀;
回顾式处理模块,用于将所述待处理图像数据进行回顾式处理,其中,所述回顾式处理用于消除所述待处理图像数据中外部干扰导致的亮度不均匀;
所述重建模块,还用于重建回顾式处理后的待处理图像数据,得到目标图像;
所述回顾式处理模块,还用于所述将所述待处理图像数据进行回顾式处理,重建回顾式处理后的待处理图像数据,得到目标图像的步骤,包括:获取不均匀乘性场和图像噪声;基于所述不均匀乘性场和所述图像噪声处理所述待处理图像数据,并重建处理后的待处理图像数据,得到目标图像,计算公式如下:
Figure DEST_PATH_IMAGE009
其中,
Figure 8460DEST_PATH_IMAGE002
为所述待处理图像数据,
Figure 470665DEST_PATH_IMAGE003
为所述不均匀乘性场,
Figure 831852DEST_PATH_IMAGE010
为所述图像噪声,
Figure 37705DEST_PATH_IMAGE005
为所述目标图像;
所述回顾式处理模块,还用于所述获取预扫描数据,基于所述预扫描数据得到所述预扫描数据对应预扫描图像的亮度分布信息之后,所述根据所述亮度分布信息处理所述正式扫描图像数据,得到待处理图像数据之前,还包括:将所述亮度分布信息进行回顾式处理,得到校正后的亮度分布信息,计算公式如下:
Figure DEST_PATH_IMAGE011
其中,
Figure 623407DEST_PATH_IMAGE007
为所述预扫描图像的所述亮度分布信息,
Figure 397459DEST_PATH_IMAGE012
为所述校正后的亮度分布信息。
7.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至5中任一项所述方法的步骤。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至5中任一项所述的方法的步骤。
CN201910383106.7A 2019-05-09 2019-05-09 医学图像处理方法、装置、计算机设备和存储介质 Active CN110265120B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910383106.7A CN110265120B (zh) 2019-05-09 2019-05-09 医学图像处理方法、装置、计算机设备和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910383106.7A CN110265120B (zh) 2019-05-09 2019-05-09 医学图像处理方法、装置、计算机设备和存储介质

Publications (2)

Publication Number Publication Date
CN110265120A CN110265120A (zh) 2019-09-20
CN110265120B true CN110265120B (zh) 2021-10-22

Family

ID=67914406

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910383106.7A Active CN110265120B (zh) 2019-05-09 2019-05-09 医学图像处理方法、装置、计算机设备和存储介质

Country Status (1)

Country Link
CN (1) CN110265120B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112986881B (zh) * 2019-12-17 2022-12-20 上海联影医疗科技股份有限公司 磁共振系统校准方法、成像方法以及磁共振系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103424727A (zh) * 2012-05-23 2013-12-04 深圳市贝斯达医疗器械有限公司 一种磁共振图像亮度不均匀修正算法
CN103445803A (zh) * 2013-09-09 2013-12-18 深圳先进技术研究院 基于正弦图的ct系统射束硬化消除方法及其系统
CN104486993A (zh) * 2013-06-10 2015-04-01 株式会社东芝 磁共振成像装置
CN108022229A (zh) * 2016-10-31 2018-05-11 北京东软医疗设备有限公司 一种图像处理方法及装置
KR20190014325A (ko) * 2017-08-02 2019-02-12 재단법인 아산사회복지재단 구조유사도를 이용한 아티펙트 저감방법 및 프로그램 및 의료영상획득장치

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103424727A (zh) * 2012-05-23 2013-12-04 深圳市贝斯达医疗器械有限公司 一种磁共振图像亮度不均匀修正算法
CN104486993A (zh) * 2013-06-10 2015-04-01 株式会社东芝 磁共振成像装置
CN103445803A (zh) * 2013-09-09 2013-12-18 深圳先进技术研究院 基于正弦图的ct系统射束硬化消除方法及其系统
CN108022229A (zh) * 2016-10-31 2018-05-11 北京东软医疗设备有限公司 一种图像处理方法及装置
KR20190014325A (ko) * 2017-08-02 2019-02-12 재단법인 아산사회복지재단 구조유사도를 이용한 아티펙트 저감방법 및 프로그램 및 의료영상획득장치

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于同态滤波方法的磁共振图像中非均匀场校正;施波 等;《江西科学》;20180210;第36卷(第1期);第15-16页 *

Also Published As

Publication number Publication date
CN110265120A (zh) 2019-09-20

Similar Documents

Publication Publication Date Title
CN109697740B (zh) 图像重建方法、装置和计算机设备
CN109308728B (zh) 正电子发射型计算机断层扫描图像处理方法及装置
CN110570483B (zh) 扫描方法、装置、计算机设备和存储介质
US8781187B2 (en) Methods, apparatuses, and computer program products for identifying a region of interest within a mammogram image
CN110261804B (zh) 磁共振扫描方法、装置、计算机设备和存储介质
CN109712694B (zh) 图像处理方法、装置、系统、图像处理工作站
CN112365413A (zh) 图像处理方法、装置、设备、系统及计算机可读存储介质
CN111157935B (zh) 磁共振成像方法、装置、存储介质及计算机设备
CN110265120B (zh) 医学图像处理方法、装置、计算机设备和存储介质
CN110598025A (zh) 阅片方法、系统、装置、设备和存储介质
CN110570488B (zh) 图像校正方法、计算机设备和存储介质
CN110244249B (zh) 磁共振扫描方法、装置、医学扫描设备和存储介质
CN115526948A (zh) 磁共振图像的射频场均匀性校正方法、装置及磁共振系统
CN114494014A (zh) 一种磁共振图像超分辨率重建方法及装置
CN110264435B (zh) 低剂量mip图像的增强方法、装置、计算机设备和存储介质
CN113393427A (zh) 斑块分析方法、装置、计算机设备和存储介质
CN113256750B (zh) 医疗图像风格重建方法、装置、计算机设备和存储介质
KR102429284B1 (ko) 컴바인이 적용되는 자기 공명 영상 처리 장치 및 그 방법
CN116958289A (zh) 图像处理方法、装置、计算机设备和存储介质
CN117635508A (zh) 医学图像的确定方法、装置、设备、存储介质和程序产品
CN116342734B (zh) 基于ct灌注影像的运动校正和伪影筛选方法及计算机设备
CN115861470A (zh) 图像伪影校正方法、装置、设备、存储介质和程序产品
CN113096775A (zh) 图像保存方法、装置、计算机设备和存储介质
CN117635511A (zh) 医学图像处理方法、装置、计算机设备和存储介质
CN114202593A (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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Applicant after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

GR01 Patent grant
GR01 Patent grant