平码五不中公式规律
  • / 10
  • 下载费用:30 金币  

图像金属伪影校正方法.pdf

关 键 ?#21097;?/dt>
图像 金属 校正 方法
  专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
摘要
申请专利号:

CN201611246225.0

申请日:

2016.01.21

公开号:

CN106683143A

公开日:

2017.05.17

当前法律状态:

实审

?#34892;?#24615;:

审中

法?#19978;?#24773;: 实质审查的生效IPC(主分类):G06T 11/00申请日:20160121|||公开
IPC分类号: G06T11/00 主分类号: G06T11/00
申请人: 上海联影医疗科技有限公司
发明人: 王毅
地址: 201807 上海市嘉定区嘉定工业区城北路2258号
优先权:
专利代理机构: 代理人:
PDF完整版下载: PDF下载
法律状态
申请(专利)号:

CN201611246225.0

授权公告号:

|||

法律状态公告日:

2017.06.09|||2017.05.17

法律状态类型:

实质审查的生效|||公开

摘要

本发明公开了一种图像金属伪影校正方法,包括:利用插值法对待校正图像进行金属伪影校正,以获取第一校正后图像;获取待校正图像相对于第一校正后图像的误差图像;根据该误差图像获取第一校正后图像的信息熵,并根据该信息熵确定第一校正后图像引入伪影的程度;根据该第一校正后图像引入伪影的程度调整所述误差图像的权重,并从所述待校正图像中去除调整权重后的误差图像,以获取第二校正后图像;分割出所述待校正图像的高频部分图像及所述第二校正后图像的低频部分图像,并对所述高频部分图像及低频部分图像进行融合,以获得第三校正后图像。本发明所提供的方法及装置可在现有伪影去除方法的基础上,?#34892;б种?#26032;引入伪影,进一步提高图像质量。

权利要求书

1.一种图像金属伪影校正方法,其特征在于,包括:
利用插值法对待校正图像进行金属伪影校正,以获取第一校正后图像;
获取待校正图像相对于第一校正后图像的误差图像;
根据该误差图像获取第一校正后图像的信息熵,并根据该信息熵确定第一校正后图像
引入伪影的程度;
根据该第一校正后图像引入伪影的程度调整所述误差图像的权重,
并从所述待校正图像中去除调整权重后的误差图像,以获取第二校正后图像;
分割出所述待校正图像的高频部分图像及所述第二校正后图像的低频部分图像,并对
所述高频部分图像及低频部分图像进行融合,以获得第三校正后图像。
2.根据权利要求1所述的金属伪影校正方法,其特征在于,所述根据该误差图像获取第
一校正后图像的信息熵包括:
划分所述误差图像与待校正图像每个像素的邻域矩阵;
调整所述误差图像像素邻域矩阵的权重,并根据调整权重后的所述误差图像像素邻域
矩阵与所述待校正图像像素邻域矩阵之差获取所述信息熵。
3.根据权利要求2所述的金属伪影校正方法,其特征在于,获取使所述信息熵最小时相
应的所述误差图像邻域矩阵的权重,将该权重作为所述误差图像的权重。
4.根据权利要求2所述的金属伪影校正方法,其特征在于,根据所述误差图像中金属伪
影去除的程度,确定所述邻域矩阵的尺寸。
5.根据权利要求2所述的伪影校正方法,其特征在于,根据所述待校正图像中的金属图
像形态,确定所述邻域矩阵的尺寸。
6.根据权利要求2所述的伪影校正方法,所述邻域矩阵的尺寸为9-31单位像素。
7.根据权利要求2所述的伪影校正方法,其特征在于,还包括:根据设定的视场对所述
待校正图像、误差图像、第一校正后图像、第二校正后图像中的至少一种进行压缩。

说明书

图像金属伪影校正方法

本申请是申请日为2016年1月21日、申请号为“201610038856.7”、发明名称为“一
种计算机断层成像伪影校正方法及装置”的发明专利申请的分案申请。

【技术领域】

本发明涉及计算机断层成像技术领域,尤其涉及一种计算机断层成像伪影的校正
方法及装置。

【背景技术】

计算机断层成像是用射线对人体的特定部位按一定厚度的层面进行扫描,根据不
同的人体组织对射线的吸收能力不同,对扫描数据利用计算机重建出断层面图像的技术。

在计算机断层扫描过程?#26657;?#22240;被扫描物体中含有金属或其它高密度物?#23454;?#33268;重建
后的图像中存在的伪影称为金属伪影。金属伪影的存在会?#26723;?#22270;像质量,并且可能会影响
医生的诊断。因此,在计算机断层扫描的成像?#26657;?#21435;金属伪影(Metal Artifact Reduction,
MAR)即金属伪影校正,所要解决的就是去除由于金属的存在而引入的伪影,恢复被伪影破
坏或掩盖的组织,以便于用户观察。

现有技术?#26657;?#23384;在各种金属伪影的校正方法.这些校正方法大致可以分为迭代重
建方法(迭代法)及投影插值方法(插值法)。但无论现有技术中何种形式的伪影校正,都有
可能引入新的伪影,这部?#20013;?#24341;入的伪影同样会影响图像质量,因而需要对其进行校正。

因此,需要提出一种新的计算机断层成像伪影校正方法,在现有伪影去除方法的
基础上,?#34892;б种?#22240;伪影校正新引入伪影,进一步提高图像质量。

【发明内容】

本发明解决的是现有的计算机断层成像图像伪影校正方法实施过程中出?#20013;?#24341;
入伪影的问题。

为解决上述问题,本发明提出一种计算机断层成像伪影校正方法,包括:接收待校
正图像;对所述待校正图像进行伪影校正,以获取第一校正后图像;获取待校正图像相对于
第一校正后图像的误差图像;根据第一校正后图像引入伪影的程度,调整所述误差图像的
权重,并从所述待校正图像中去除调整权重后的误差图像,以获取第二校正后图像;对所述
待校正图像及第二校正后图像进行频率分割及融合,获得第三校正后图像。

可选地,所述第一校正后图像引入伪影的程度,由所述第一校正后图像的信息熵
进行确定。

可选地,所述信息熵的获取包括:划分所述误差图像与待校正图像每个像素的邻
域矩阵;调整所述误差图像像素邻域矩阵的权重,并根据调整权重后的所述误差图像像素
邻域矩阵与所述待校正图像像素邻域矩阵之差获取所述信息熵。

可选地,获取使所述信息熵最小时相对应的所述误差图像邻域矩阵的权重,将该
权重作为所述误差图像的权重。

可选地,根据所述误差图像中伪影去除的程度,确定所述邻域矩阵的尺寸。

可选地,根据所述待校正图像中的金属图像形态,确定所述邻域矩阵的尺寸。

可选地,所述邻域矩阵的尺寸范围为9-31单位像素。

可选地,还包括根据设定的视场对所述待校正图像、误差图像、第一校正后图像、
第二校正后图像中的至少一种进行压缩。

可选地,所述频率分割及融合包括:分割出所述待校正图像的高频部分图像及所
述第二校正后图像的低频部分图像,并对所述高频部分图像及低频部分图像进行融合。

本发明?#22266;?#20379;一种计算机断层成像伪影校正装置,其特征在于,包括:

输入单元,用于接收待校正图像;第一处理单元,用于对待校正图像进行伪影校
正,以生成第一校正后图像;第二处理单元,用于获取待校正图像相对于第一校正后图像的
误差图像,及根据第一校正后图像引入伪影的程度,调整所述误差图像的权重,并从所述待
校正图像中去除调整权重后的误差图像,以获取第二校正后图像;第三处理单元,用于获得
所述待校正图像的高频部分图像及所述第二校正后图像的低频部分图像,并对所述高频部
分图像及低频部分图像进行融合,获得第三校正后图像。

本发明方案根据原有伪影校正引入伪影的程度,调整误差图像的权重,并在待校
正图像中去除调整权重后的误差图像以形成新的校正后图像,?#34892;б种?#20102;新伪影的产生,
提高了图像质量;同时该方案实现过程对原伪影校正方法及装置?#35272;?#24615;低,因而适用范围
广,可在任何可能会引入新伪影的现有金属伪影校正方法及装置基础上实现;进一步地,频
率分割及融合减少了由于权重系数的差异引起的图像马赛克现象,使得校正后图像更加自
然;进一步地,对图像进行压缩,及根据误差图像中伪影去除的程度对邻域矩阵进行划分,
?#26723;?#20102;?#20302;?#30340;计算量,提升了校正速?#21462;?br />

【附图说明】

图1是本发明的计算机断层成像?#20302;?#30340;结构示意图;

图2是本发明一实施例中伪影校正方法流程示意图;

图3是本发明一实施例中第一校正后图像信息熵的求取流程示意图;

图4是本发明一实施例中伪影校正装置结构示意图。

【具体实施方式】

为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图对本发明
的具体实施方式做详细的说明。

图1是一种计算机断层成像?#20302;?#30340;结构示意图,如图1所示,计算机断层成像?#20302;?br />100包括机架110,所述机架110具有围绕?#20302;持?#32447;旋转的可旋转的部分130。可旋转的部分
130具有相对设置的X射线源131和X射线探测器132的X射线?#20302;场?br />

计算机断层成像?#20302;?00还具有检查床120,在进行检查时,患者在该检查床120上
可以沿着Z轴方向被推入到扫描腔体中。X射线源131绕S轴旋转,探测器132相对于X射线源
131一起运动,以采集投影测量数据,这些数据在之后被用于重建图像。还可以进行螺旋扫
描,在螺旋扫描期间,通过患者沿着S轴的连续运动和X射线源131的同时旋转,X射线源131
相对于患者产生螺旋轨迹。

所述计算机断层成像?#20302;?00还可以包括控?#39057;?#20803;和图像重建单元,所述控?#39057;?br />元用于在扫描过程中根据特定的扫描协议控制计算机断层成像?#20302;?00的各部件。所述图
像重建单元用于根据探测器132采样的待校正数据重建出图像。

以上,仅以示例方式阐释了可使用本发明所提供金属伪影校正方法的计算机断层
成像设备,本领域技术人员理解,如使用X射线的C型臂?#20302;?#31561;设备,或组合式医学成像?#20302;?br />(例如:组合式正电子发射断层成像-计算机断层成像,Positron Emission Tomography-
Computed tomography Tomography,PET-CT),或使用其它类型射线的断层成像设备等,均
可适用本发明所述校正方法及装置,本发明对计算机断层成像设备的类型与结构并不做具
体限定。

当受检对象在上述任一种类的计算机断层成像设备中进行扫描成像时,因金属或
高密度物体的存在会导致伪影的存在,影响图像成像质量及导致用户阅图不便,因而需对
此类伪影进行校正。

图2是本发明一实施例中伪影校正方法流程示意图。参照图2,在该实施例?#26657;?br />

执行步骤S1,接收待校正图像。该待校正图像由上述计算机断层成像设备扫描重
建获得。

执行步骤S2,对待校正图像进行伪影校正,获取第一校正后图像。此步骤?#26657;?#24182;不
限定该伪影校正的具体实现方案。因本发明方案的目的在于解决现有技术中存在的各种金
属伪影的校正方法所存在的新引入伪影缺陷,?#25105;?#19968;种去金属伪影算法或多种去金属伪影
算法的组合(无论迭代法还是插值法,亦或其它种类的伪影校正方法),只要有可能引入新
的伪影,均可应用于此。对待校正图像进行伪影校正后,可获取第一校正后图像。

执行步骤S3,获取待校正图像相对于第一校正后图像的误差图像。在本实施例?#26657;?br />定义Iori为待校正图像,Icorr为第一校正后图像,Ierr为误差图像。则:

Ierr=Iori-Icorr

即,可通过比较待校正图像与第一校正后图像的差异来获取误差图像。

执行步骤S4,根据第一校正后图像引入伪影的程度,调整误差图像的权重,并从待
校正图像中去除调整权重后的误差图像,获取第二校正后图像。第一校正后图像引入伪影
的程度可以通过对误差图像的判断得出,因此,通过对误差图像的每个像素点分配权重(分
配的原则是,对新引入的伪影区域分配少的权重,待校正图像原有的伪影分配多的权重),
再在待校正图像中去除调整权重后的误差图像,即可实现?#31181;?#26032;引入伪影的效果。

在本实施例?#26657;?#31532;一校正后图像新引入伪影的程度,可通过第一校正后图像的信
息熵进行确定。图3示出了第一校正后图像信息熵的求取流程:

首先,执行步骤S401,划分误差图像与待校正图像每个像素的邻域矩阵。具体地,
如对于一个N×N像素的图像I而言,对其图像的每个像素(投影角为p,通道为k)计算一个邻
域矩阵Np,k(I)。经过测试发现,该邻域矩阵的尺寸(MxM)不仅影响着运算速度,而?#19968;?#30452;接
影响到之后步骤中权重的分配,所以邻域大小的划分可选取经验值的方式确定。例如,在本
实施例?#26657;?#37051;域的大小根据误差图像中伪影去除的程?#30830;?#37197;不同的值。一般来说,对于步骤
S2中对待校正图像伪影去除较多的应用场景,邻域划分可以?#28304;螅?#21453;之则可以稍小。进一步
地,邻域的取值范围可在9-31单位像素之间。此外,在对图像边缘像素取邻域矩阵时,可以
在图像边缘填零扩充图像边缘。

此处,作为本实施例的一个变化例,邻域的尺寸还可以根据金属形态信息进行确
定。此处的金属形态信息,指待校正图像中金属物质或较高密度物质的形态信息,其直接影
响到待校正图像中数据破坏的程度:理论上,仅当金属或高密度物体为一规则圆形时,现有
校正方案可通过完全替换原有数据的方式?#34892;?#21435;除金属伪影(即新引入伪影较少)。而实际
中很少有规则的物体,如常见的脊柱钉子等植入物的形状多以?#36824;?#21017;形态呈现。

设金属形态信息可用形态指数来进行表征:在某一断层面?#26657;?#23556;线经过该断层面
一规则的圆形物体,则无论从任一角度进行投影,其面积均是一致的,可将此理想的金属投
影域面积定义为理想面积,并作为判断金属形态的一项参照,其获取公式可为:


该公式?#26657;琒Imetal指的是金属图像中像素个数;Spacechannel为通道分辨?#21097;?br />Spacepixel为像素分辨?#21097;琋view为投影角个数,SIdealPmetal是指规则金属图像(金属图像域像素
个数和实际金属图像的金属像素个数相同)投影后的非零像素个数(即理想面积)。

通过比较金属图像投影数据的理想面积与实际面积可获得金属形态指数。进一步
地,设SPmetal为实际金属图像投影后的非零像素个数(即金属图像实际面积),则金属形态指
数RD的获取公式为:


根据本实施例的一个变化实施例,步骤S2中的伪影校正方法为投影域加权校正:
金属形态信息会影响到步骤S2的投影域加权校正是否进行及加权的强度,进而影响第一校
正后图像的质量,因而,第一校正后图像信息熵求取时邻域的大小,可依据金属形态信息进
行划分。一般来说,若金属图像实际面积与理想面积较为接近(如金属形态指数RD小于或接
近于1),则进行较少强度的加权甚至强度为零的加权,对应选择较大的邻域;若与理想面积
相差较大(如金属形态指数RD为2左右),则进行较高强度的加权校正,对应选择较小的邻
域。

完成划分误差图像与待校正图像每个像素的邻域矩阵后,继续如图3所示,执行步
骤S402,基于划?#20540;?#37051;域矩阵进行图像域加权。具体地,对于每一个误差图像中每个像素的
邻域矩阵进行加权,和待校正图像的邻域矩阵做差得到第一校正后图像的邻域矩阵:

Np,k(Icw)=Np,k(Iori)-wp,k×Np,k(Ierr)

式?#26657;琁cw指第一校正后图像,wp,k为邻域矩阵中每个像素的权重,其大小决定了对
误差图像中伪影的加重或减轻。

执行步骤S403,计算第一校正后图像的信息熵。本实施例中采用信息熵来表征第
一校正后图像中的结构信息。具体地,第一校正后图像的信息熵为:


其?#26657;珽ntropy(Np,k(Icw))表示第一校正后图像像素邻域矩阵的信息熵,P(Np,k
(Icw))表示第一校正后图像像素的邻域矩阵的先验概率函数(该函数可通过直方图方式获
取),i表示邻域矩阵中的像素下标,n表示像素总数。

根据S402及S403中第一校正后图像信息熵的求取过程可知,不同的像素权重wp,k
对应不同的信息熵,通过调整wp,k,可求出使第一校正后图像信息熵Entropy(Np,k(Icw))最
小时的权重即arg w min(Entropy(Np,k(Icw)))。将该权重作为误差图像调整后的权重,并
从待校正图像中去除按该权重进行调整后的误差图像,以得到第二校正后图像。

因通过上述权重调整计算有可能会在第二校正后图像中引入一些马赛克效应,因
此,如图2所示,还需执行步骤S5,对待校正图像及第二校正后图像进行频率分割及融合。

具体地,可通过对待校正图像及第二校正后图像进行频率分割,获取未校正图像
的高频部分图像和校正图像的低频部分图像并对其进行融合以实现,例如可选用高斯低通
滤波器G(σ)通过频域卷积实?#20540;?#36890;滤波,其?#26657;?br />


loriLow=Iori*G(σ)

Icorr2Low=Icorr2*G(σ)

式?#26657;?#20026;一经验值,可选范围为1.2-1.3之间,e表示指数函数的底,x表示距离,
Iori代表待校正图像,IoriLow代表待校正图像的低频部分图像,Icorr2代表第二校正后图
像,Icorr2Low代表第二校正后图像的低频部分图像。当然,?#37096;?#36873;用其它类型的低通滤波
器。

优选地,在待校正图像做滤波之前,可做自适应滤波处理以消除大部?#20540;?#39640;频细
条伪影。得到待校正图像和第二校正后图像的低频部分图像后,可以通过像素相减得到对
应的第二校正后图像的高频部分图像IoriHigh:

IoriHigh=Iori-IoriLow

合并未校正图像的高频图像和校正图像的低频图像,即可得到第三校正后图像。

在上述伪影校正的过程?#26657;?#36824;可根据设定视场对所述待校正图像、误差图像、第一
校正后图像、第二校正后图像中的至少一种进行压缩,以提高计算效率。

本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可
以通过程序来指令相关的硬件来完成,该程序可以存储于计算机可读存储介质?#26657;?#23384;储介
质可以包括但不限于:软盘、光盘、CD-ROM、磁光盘、ROM(只读存储器)、RAM(随机存取存储
器)、EPROM(可擦除可编程只读存储器)、EEPROM(电可擦除可编程只读存储器)、磁卡或光
卡、闪存、或适于存储机器可执行指令的其他类型的介质/机器可读介质。

图4示出了本发明一实施例中计算机断层成像伪影校正装置的结构示意图。包括:

输入单元,用于接收待校正图像;

第一处理单元,用于对待校正图像进行伪影校正,以生成第一校正后图像;

第二处理单元,用于获取待校正图像相对于第一校正后图像的误差图像,及根据
第一校正后图像引入伪影的程度,调整所述误差图像的权重,并从所述待校正图像中去除
调整权重后的误差图像,以获取第二校正后图像;

第三处理单元,用于获得所述待校正图像的高频部分图像及所述第二校正后图像
的低频部分图像,并对所述高频部分图像及低频部分图像进行融合,获得第三校正后图像;
以及输出单元,用于将校正后的图像输出至用户或后续计算设备。

本实施例中的输入单元、第一处理单元及输出单元,可在现有技术中?#25105;?#19968;种伪
影校正装置的基础上实现,因而,本发明的提供的伪影校正装置具备良好的兼容性及较低
的实施成本。

本发明?#26657;?#21508;实施例采用递进式写法,重点描述与前述实施例的不同之处,各实施
例中的相同方法或结构参照前述实施例的相同部分。

本发明虽然已以较佳实施例公开如上,但其并不是用来限定本发明,任?#20266;?#39046;域
技术人员在不脱离本发明的精神和范围内,都可以利用上述揭示的方法?#22270;?#26415;内容对本发
明技术方案做出可能的变动和修改,因此,凡是未脱离本发明技术方案的内容,依据本发明
的技术实质对以上实施例所作的任何简单修改、等同变化及修饰,均属于本发明技术方案
的保护范围。

关于本文
本文标题:图像金属伪影校正方法.pdf
链接地址:http://www.pqiex.tw/p-6079643.html
关于我们 - 网站声明 - 网?#38236;?#22270; - 资源地图 - 友情链接 - 网站客服 - 联系我们

[email protected] 2017-2018 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备17046363号-1 
 


收起
展开
平码五不中公式规律 ak娱乐是什么 杰克棋牌 有山西快乐10分的平台 恐怖嘉年华下载 pt老虎机平台大全 直销超市赚钱吗 天顺娱乐怎么样 天天捕鱼2下载安装 宝龙 单挑机黑红梅方100局