CN102768759A - 一种术中ct图像射束硬化伪影校正方法及装置 - Google Patents

一种术中ct图像射束硬化伪影校正方法及装置 Download PDF

Info

Publication number
CN102768759A
CN102768759A CN2012102295380A CN201210229538A CN102768759A CN 102768759 A CN102768759 A CN 102768759A CN 2012102295380 A CN2012102295380 A CN 2012102295380A CN 201210229538 A CN201210229538 A CN 201210229538A CN 102768759 A CN102768759 A CN 102768759A
Authority
CN
China
Prior art keywords
image
projection
data
separation
original
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.)
Granted
Application number
CN2012102295380A
Other languages
English (en)
Other versions
CN102768759B (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.)
ANKE HIGH-TECH Co Ltd SHENZHEN CITY
Original Assignee
ANKE HIGH-TECH Co Ltd SHENZHEN CITY
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 ANKE HIGH-TECH Co Ltd SHENZHEN CITY filed Critical ANKE HIGH-TECH Co Ltd SHENZHEN CITY
Priority to CN201210229538.0A priority Critical patent/CN102768759B/zh
Publication of CN102768759A publication Critical patent/CN102768759A/zh
Application granted granted Critical
Publication of CN102768759B publication Critical patent/CN102768759B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Abstract

本发明公开一种术中CT图像射束硬化伪影校正方法及装置。本发明从原始CT图像中分离出仅含高密度组织的分离CT图像，并对这两种图像分别进行重投影，然后对重投影数据进行校正，使用校正后的重投影数据进行图像重建得到误差CT图像，利用该误差CT图像对原始CT图像进行校正，最后获得校正CT图像。通过本发明的校正，CT图像的伪影基本被消除或减弱，实现了对原始CT图像校正的目的。

Description

一种术中CT图像射束硬化伪影校正方法及装置
技术领域
本发明涉及术中CT图像处理领域，尤其涉及一种术中CT图像射束硬化伪影校正方法及装置。
背景技术
因为术中CT（一种CT医疗诊断设备）具有体积小、质量轻、使用灵活等特点，因而应用越来越广泛，与传统CT一样，术中CT也可对病人进行断层扫描，得到感兴趣区的切片图像，应用于早期诊疗、无创诊断等领域。术中CT的应用范围广，在需要进行CT扫描的场合中都可使用，比如手术室、ICU重症监护室、急诊室等，通过使用术中CT，可以达到提高手术成功率、延长病人寿命的目的。
但是术中CT所得到的图像中广泛存在的一种伪影，其主要是由射束硬化现象引入的，所谓射束硬化现象是指低能量的射线易于被物质吸收衰减的现象，射束硬化还会带来条状伪影，射束硬化对图像质量的影响严重，特别是在对头颅等组织进行扫描时，这种射束硬化伪影更加明显。
因此，现有技术还有待于改进和发展。
发明内容
鉴于上述现有技术的不足，本发明的目的在于提供一种术中CT图像射束硬化伪影校正方法及装置，旨在解决现有技术中术中CT射束硬化现象带来的伪影问题。
本发明的技术方案如下：
一种术中CT图像射束硬化伪影校正方法，其中，包括步骤：
A、通过设置一预定的分割阀值从原始CT图像中分离提取高密度组织区域获得仅包含高密度组织区域的分离CT图像；
B、使用平行光源在预定的角度范围内分别对原始CT图像和分离CT图像进行重投影获得原始CT图像重投影数据和分离CT图像重投影数据；
C、根据所述原始CT图像重投影数据设置一校正区域和重置区域，将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据；
D、对所述分离CT图像重投影校正数据进行CT图像重建获得误差CT图像；
E、利用所述误差CT图像对所述原始CT图像进行校正处理获得校正CT图像。
所述术中CT图像射束硬化伪影校正方法，其中，所述步骤A的分割阀值大于1000Hu，当所述原始CT图像的目标区域的CT值大于所述分割阀值时，则该目标区域为高密度组织区域。
所述术中CT图像射束硬化伪影校正方法，其中，所述步骤B的角度范围为                                                
Figure 802433DEST_PATH_IMAGE001
，平行光源的角度增量
Figure 471312DEST_PATH_IMAGE002
，所述平行光源共有256条射线，所述分离CT图像重投影数据为，
Figure 219005DEST_PATH_IMAGE004
为当前射线入射到原始CT图像的像素，
Figure 213823DEST_PATH_IMAGE006
为平行光源的第i个重投影角度，
Figure 378088DEST_PATH_IMAGE007
为常量。
所述术中CT图像射束硬化伪影校正方法，其中，所述步骤C具体包括：
C1、设置常量
Figure 619714DEST_PATH_IMAGE008
和，当
Figure 131784DEST_PATH_IMAGE010
，则判定该投影数据
Figure 732529DEST_PATH_IMAGE011
为重置区域，当，则判定该投影数据
Figure 540265DEST_PATH_IMAGE011
为校正区域，其中，
Figure 507084DEST_PATH_IMAGE011
为原始CT图像重投影数据在角度s处的投影数据；
C2、将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据。
所述术中CT图像射束硬化伪影校正方法，其中，所述步骤D具体包括：
D1、设置与所述分离CT图像重投影校正数据对应的滤波核；
D2、对所述分离CT图像重投影校正数据和滤波核分别进行快速傅里叶变换获得分离CT图像重投影校正数据频域以及滤波核频域；
D3、将所述分离CT图像重投影校正数据频域以及滤波核频域在频域内相乘获得频域的滤波数据；
D4、对所述滤波数据进行离散傅里叶逆变换获得滤波投影数据；
D5、对所述滤波投影数据进行反投影计算，获得误差CT图像。
所述术中CT图像射束硬化伪影校正方法，其中，所述步骤E具体包括：
将原始CT图像与误差CT图像中的每一像素进行逐行累加，获得校正CT图像。
所述术中CT图像射束硬化伪影校正方法，其中，所述步骤A之前还包括：
利用探测器探测多个视角的投影帧数据，并对所述投影帧数据进行切片重建获取重建切片的原始CT图像。
一种术中CT图像射束硬化伪影校正装置，其中，包括：
分离CT图像获取模块，用于通过设置一预定的分割阀值从原始CT图像中分离提取高密度组织区域获得仅包含高密度组织区域的分离CT图像；
重投影模块，用于使用平行光源在预定的角度范围内分别对原始CT图像和分离CT图像进行重投影获得原始CT图像重投影数据和分离CT图像重投影数据；
重投影数据校正模块，用于根据所述原始CT图像重投影数据设置一校正区域和重置区域，将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据；
误差CT图像重建模块，用于对所述分离CT图像重投影校正数据进行CT图像重建获得误差CT图像；
校正处理模块，用于利用所述误差CT图像对所述原始CT图像进行校正处理获得校正CT图像；
所述分离CT图像获取模块、重投影模块、重投影数据校正模块、误差CT图像重建模块、校正处理模块依次连接。
所述术中CT图像射束硬化伪影校正装置，其中，所述重投影数据校正模块包括：
区域设置单元，用于设置常量
Figure 13152DEST_PATH_IMAGE008
和
Figure 229370DEST_PATH_IMAGE009
，当
Figure 111875DEST_PATH_IMAGE010
，则判定该投影数据
Figure 57834DEST_PATH_IMAGE011
为重置区域，当
Figure 734803DEST_PATH_IMAGE012
，则判定该投影数据
Figure 172738DEST_PATH_IMAGE011
为校正区域，其中，为原始CT图像重投影数据在角度s处的投影数据；
重投影数据校正单元，用于将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据；
所述区域设置单元与重投影数据校正单元相连接。
所述术中CT图像射束硬化伪影校正装置，其中，所述误差CT图像重建模块包括：
滤波核设置单元，用于设置与所述分离CT图像重投影校正数据对应的滤波核；
快速傅里叶变换单元，用于对所述分离CT图像重投影校正数据和滤波核分别进行快速傅里叶变换获得分离CT图像重投影校正数据频域以及滤波核频域；
滤波数据获取单元，用于将所述分离CT图像重投影校正数据频域以及滤波核频域在频域内相乘获得频域的滤波数据；
离散傅里叶变换单元，用于对所述滤波数据进行离散傅里叶逆变换获得滤波投影数据；
误差CT图像获取单元，用于对所述滤波投影数据进行反投影计算，获得误差CT图像；
所述滤波核设置单元、快速傅里叶变换单元、滤波数据获取单元、离散傅里叶变换单元、误差CT图像获取单元依次连接。
有益效果：本发明射束硬化伪影的校正方法及装置，从原始CT图像中分离出仅含高密度组织的分离CT图像，并对这两种图像分别进行重投影，然后对重投影数据进行骨校正，使用校正后的重投影数据进行图像重建得到误差CT图像，利用该误差CT图像对原始CT图像进行校正，最后获得校正CT图像。通过本发明的校正，CT图像的伪影基本被消除或减弱， 实现了对原始CT图像校正的目的，并且本发明降低了数据处理量，加快了图像处理速度，从而提高了校正效率。
附图说明
图1为本发明术中CT图像射束硬化伪影校正方法较佳实施例的流程图。
图2为图1所示方法中CT图像重建的具体流程图。
图3为本发明实施例中原始CT图像的示意图。
图4为采用本发明的校正方法后获得的校正CT图像的示意图。
图5为本发明术中CT图像射束硬化伪影校正装置较佳实施例的结构框图。
图6为图5所示装置中重投影数据校正模块的结构框图。
图7为图5所示装置中误差CT图像重建模块的结构框图。
具体实施方式
本发明提供一种术中CT图像射束硬化伪影校正方法及装置，为使本发明的目的、技术方案及效果更加清楚、明确，以下对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
请参阅图1，图1为本发明术中CT图像射束硬化伪影校正方法较佳实施例的流程图，如图所示，包括步骤：
S101、通过设置一预定的分割阀值从原始CT图像中分离提取高密度组织区域获得仅包含高密度组织区域的分离CT图像；
S102、使用平行光源在预定的角度范围内分别对原始CT图像和分离CT图像进行重投影获得原始CT图像重投影数据和分离CT图像重投影数据；
S103、根据所述原始CT图像重投影数据设置一校正区域和重置区域，将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据；
S104、对所述分离CT图像重投影校正数据进行CT图像重建获得误差CT图像；
S105、利用所述误差CT图像对所述原始CT图像进行校正处理获得校正CT图像。
在步骤S101中，原始CT图像是从对探测器探测到的多个视角的投影帧数据进行切片重建获取的，下文会具体对其进行说明。在获取了原始CT图像后，需要从原始CT图像中获得仅含高密度组织区域的分离CT图像，高密度组织是指组织密度相对较高的物质，例如骨骼、金属等。具体的分离方法为设置一分割阀值
Figure 534766DEST_PATH_IMAGE013
，对于原始CT图像中CT值大于分割阀值
Figure 648216DEST_PATH_IMAGE013
的区域则为高密度组织区域，对于原始CT图像中CT值小于分割阀值
Figure 573446DEST_PATH_IMAGE013
的区域则为低密度组织区域，并将该低密度组织区域的CT值重置为0，如此，即可获得仅含高密度组织区域的分离CT图像。本发明中，对于不同的组织结构，所设定的分割阀值也不同，对于高密度组织，分割阀值一般设置为1000Hu以上。本发明中，CT值的单位为Hu。
假设需要对原始CT图像的头骨进行分离，获得仅含头骨的分离CT图像，用
Figure 391230DEST_PATH_IMAGE014
表示原始CT图像，
Figure 718306DEST_PATH_IMAGE015
表示坐标，其大小为
Figure 2657DEST_PATH_IMAGE016
，其中
Figure 149604DEST_PATH_IMAGE017
，而设置头骨的CT值为1200Hu，那么可设置分割阀值
Figure 646445DEST_PATH_IMAGE018
Hu，通过下面的算法来获得分离得到仅含头骨的分离CT图像
Figure 828027DEST_PATH_IMAGE019
；
FOR i←1:512,j←1:512
DO IF f(i,j)<th
DO b(I,j)=0
                              ELSE b(I,j) = f(i,j)
End IF
如此，即可获得包含了切片组织的原始CT图像
Figure 283279DEST_PATH_IMAGE014
，以及仅包含头骨的分离CT图像
Figure 183102DEST_PATH_IMAGE019
。
S102、获取了原始CT图像和分离CT图像后，需要对他们进行重投影，以获得重投影数据，这种重投影包括对多种扫描方式下的图像进行重投影，例如平行束、扇束扫描等二维扫描方式下的图像，或者锥束扫描、螺旋扫描等立体扫描方式下的图像。而对于重投影则可使用平行光源在预定的角度范围内分别对原始CT图像和分离CT图像进行重投影，以获得原始CT图像重投影数据和分离CT图像重投影数据。
例如对于分辨率的CT图像，平行光源在的角度范围内进行曝光，平行光源的起始重投影角度用
Figure 4931DEST_PATH_IMAGE021
表示，角度增量为
Figure 126470DEST_PATH_IMAGE022
，那么第个重投影角度表示为
Figure 121288DEST_PATH_IMAGE024
，其中
Figure 918343DEST_PATH_IMAGE025
,
Figure 527179DEST_PATH_IMAGE026
。而在位置处，平行光源的相邻射线的间距为
Figure 39249DEST_PATH_IMAGE028
，平行光源中第j条射线的重投影数据可表示为
Figure 7205DEST_PATH_IMAGE029
，其中的
Figure 103337DEST_PATH_IMAGE027
表示重投影角度，l表示射线，i表示第i个重投影角度，j表示第j条射线，那么所有角度下的重投影数据表示为
Figure 814941DEST_PATH_IMAGE030
和
Figure 945708DEST_PATH_IMAGE031
，
Figure 553407DEST_PATH_IMAGE030
表示原始CT图像的重投影数据，
Figure 871255DEST_PATH_IMAGE031
表示分离CT图像的重投影数据，s、t表示坐标。重投影的总的思路是：根据重投影角度确定射线的间隔，然后确定下一个像素坐标，求解当前射线上的线积分，然后对线积分进行对数运算。
现对于重投影的实现方式用一实施例来具体说明，对原始CT图像进行重投影，x和y表示坐标，平行光源曝光的角度范围
Figure 371824DEST_PATH_IMAGE001
，角度增量
Figure 275058DEST_PATH_IMAGE002
，平行光源共有256条分布均匀地射线，在当前角度
Figure 80203DEST_PATH_IMAGE027
下，射线的间距为
Figure 133610DEST_PATH_IMAGE032
，像素的递增分量是
Figure 973390DEST_PATH_IMAGE033
，当前射线通过的入射像素是
Figure 480911DEST_PATH_IMAGE004
，则其下一个像素是
Figure 72430DEST_PATH_IMAGE034
，出射像素是
Figure 32295DEST_PATH_IMAGE035
。如此，
Figure 277332DEST_PATH_IMAGE005
的重投影表示为，这样，所有角度的重投影数据可以表示为
Figure 452278DEST_PATH_IMAGE037
，s表示平行光源的角度，x表示像素。采用上述方法同样可得到分离CT图像重投影数据为。
在步骤S103中，获取了原始CT图像重投影数据和分离CT图像重投影数据之后，需要对其进行校正，具体可根据所述原始CT图像重投影数据设置一校正区域和重置区域，将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据。例如，设置两个常值
Figure 89113DEST_PATH_IMAGE008
和
Figure 824988DEST_PATH_IMAGE009
，若满足
Figure 23888DEST_PATH_IMAGE010
， 
Figure 692767DEST_PATH_IMAGE011
为原始CT图像重投影数据在角度s处的投影数据，则判定此处投影所在的射线不包含高密度组织，并且设置；否则，判定此处投影所在的射线包含高密度组织，并且使用预先定义的预定常量代替该处的投影值
Figure 33935DEST_PATH_IMAGE040
，从而实现对重投影数据的校正获得分离CT图像重投影校正数据。本发明采用上述方法对原始CT图像和分离CT图像进行重投影，然后利用原始CT图像重投影数据对分离CT图像重投影数据进行校正，大大减少了计算量，处理运算过程都大大降低，提高了图像校正的速度和质量。
本步骤可具体细化为以下步骤：
S201、设置常量
Figure 36526DEST_PATH_IMAGE008
和，当，则判定该投影数据为重置区域，当，则判定该投影数据
Figure 353238DEST_PATH_IMAGE011
为校正区域，其中，为原始CT图像重投影数据在角度s处的投影数据，从实际经验得出，
Figure 10802DEST_PATH_IMAGE008
一般需小于1，此时判定结构将更加准确；
S202、将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据。
在步骤S104中，获取了分离CT图像重投影校正数据之后，需要对其进行CT图像重建获取误差CT图像，具体如图2所示，其包括步骤：
S301、设置与所述分离CT图像重投影校正数据对应的滤波核；
S302、对所述分离CT图像重投影校正数据和滤波核分别进行快速傅里叶变换获得分离CT图像重投影校正数据频域以及滤波核频域；
S303、将所述分离CT图像重投影校正数据频域以及滤波核频域在频域内相乘获得频域的滤波数据；
S304、对所述滤波数据进行离散傅里叶逆变换获得滤波投影数据；
S305、对所述滤波投影数据进行反投影计算，获得误差CT图像。
举例说明，假设在角度
Figure 355195DEST_PATH_IMAGE041
处，分离CT图像重投影校正数据表示为
Figure 587593DEST_PATH_IMAGE042
，对
Figure 359240DEST_PATH_IMAGE042
进行快速傅立叶变化，得到分离CT图像重投影校正数据的频域，其表达式为
Figure 309879DEST_PATH_IMAGE043
，其中
Figure 661226DEST_PATH_IMAGE044
，对滤波核
Figure 13710DEST_PATH_IMAGE045
进行快速傅立叶变化，得到滤波核频域的表达式，其中
Figure 394193DEST_PATH_IMAGE047
；在频域内计算
Figure 673864DEST_PATH_IMAGE048
，得到频域的滤波数据；对
Figure 880855DEST_PATH_IMAGE049
进行离散傅立叶逆变换，得到滤波投影数据
Figure 994304DEST_PATH_IMAGE050
；再使用前述的重投影算法进行反投影计算，对当前射线上的所有像素
Figure 919535DEST_PATH_IMAGE051
用进行累加，从而获得误差CT图像。
在步骤S105中，获得误差CT图像后，利用该误差CT图像对原始CT图像进行校正处理得到校正CT图像并存储。其主要是通过对原始CT图像
Figure 958532DEST_PATH_IMAGE014
和误差CT图像中的每一像素进行逐行累加，得到校正CT图像
Figure 992533DEST_PATH_IMAGE054
，用表达式表示为
Figure 174116DEST_PATH_IMAGE055
。例如，对于像素
Figure 629368DEST_PATH_IMAGE056
，首先对
Figure 732453DEST_PATH_IMAGE053
进行插值计算，得到更新后的误差CT图像，
Figure 298563DEST_PATH_IMAGE057
；
第一列像素可以进行如下运算，
Figure 600232DEST_PATH_IMAGE058
；
第二列像素可以进行如下运算，
Figure 226385DEST_PATH_IMAGE059
；
第三列像素可以进行如下运算，
；
这里，
Figure 311202DEST_PATH_IMAGE061
、、、
Figure 342109DEST_PATH_IMAGE064
、
Figure 250022DEST_PATH_IMAGE065
、
Figure 260703DEST_PATH_IMAGE066
、
Figure 228659DEST_PATH_IMAGE067
、
Figure 918267DEST_PATH_IMAGE068
是加权系数，满足
Figure 629871DEST_PATH_IMAGE069
，
Figure 495058DEST_PATH_IMAGE070
。
校正CT图像，通过对原始CT图像进行累加得到，
Figure 633916DEST_PATH_IMAGE071
第一行像素可以进行如下运算，
Figure 686185DEST_PATH_IMAGE072
第二行像素可以进行如下运算，
Figure 201480DEST_PATH_IMAGE073
第三行像素可以进行如下运算，
这里，
Figure 824409DEST_PATH_IMAGE075
、
Figure 895133DEST_PATH_IMAGE076
、
Figure 948539DEST_PATH_IMAGE077
、
Figure 257161DEST_PATH_IMAGE078
、
Figure 737821DEST_PATH_IMAGE079
、
Figure 295841DEST_PATH_IMAGE080
、
Figure 887360DEST_PATH_IMAGE081
、
Figure 581646DEST_PATH_IMAGE082
是加权系数，满足
Figure 92262DEST_PATH_IMAGE083
，
Figure 871999DEST_PATH_IMAGE084
。
通过这些运算，得到校正CT图像，然后使用
Figure 550422DEST_PATH_IMAGE085
转换为16位格式进行存储，其中
Figure 372885DEST_PATH_IMAGE086
运算符表示二进制与运算。图像存储格式需符合DICOM标准（一种图像格式标准），文件开头写入128个空字符(0)，写入字符”DICM”，写入标识扫描机生产商、扫描协议、图像大小等信息的图像信息单元，再写入像素的CT值，然后进行存储。
在本发明中，原始CT图像是通过术中CT机的X射线球管发射X射线，经过准直后形成扇形光束并穿透感兴趣区，衰减后的X射线达到探测器阵列并被采集，X射线球管在切片平面内围绕着被扫描区域进行旋转扫描，在每一个曝光位置，探测器阵列采集到一组投影数据，形成一个视角的投影帧，X射线球管在360度范围内进行曝光形成多帧投影，形成对切片的一次完整扫描，利用上述投影的投影数据进行切片重建即可获得原始CT图像，原始CT图像反映了物质的X射线衰减系数分布函数，也反映了物质的相对密度，在原始CT图像中，X射线衰减系数转化为CT值，用以控制原始CT图像在图像输出设备上的显示。本发明中，到达探测器阵列的X射线强度取决于切片的组织密度，不同的组织密度具有不同的X射线衰减率，探测器阵列中的探测器单元是彼此独立的。
而射束硬化的原因则主要是：X射线球管发射出的X射线是宽能带的，而且相同密度的物质，其X射线衰减系数随着X射线强度而变化，最后导致按照现有重建算法重建的原始CT图像存在射束硬化伪影。本发明就是针对上述问题，提出一种射束硬化伪影的校正方法，从原始CT图像中分离出仅含高密度组织的分离CT图像，并对这两种图像分别进行重投影，然后对重投影数据进行校正，使用校正后的重投影数据进行图像重建得到误差CT图像，利用该误差CT图像对原始CT图像进行校正，最后获得校正CT图像。
本发明的一个具体实施例，采用单排探测器，球管电压为120kV，球管电流80mA，视野为240cm，层厚为5mm，使用头模，机架倾角为20度，球管和探测器共同旋转
Figure 639918DEST_PATH_IMAGE087
，曝光1200次，通过探测器采集1200个位置的投影数据，对这些投影数据进行校正，使用滤波反投影算法进行图像重建获得原始CT图像，如图3所示，通过本发明的校正后获得的校正CT图像如图4所示，从图中可以看出，伪影现象明显减少，图像更加清晰，方便了医生对于患者的情况给出一个准确的判断。
基于上述方法，本发明还提供一种术中CT图像射束硬化伪影校正装置，如图5所示，其包括：
分离CT图像获取模块100，用于通过设置一预定的分割阀值从原始CT图像中分离提取高密度组织区域获得仅包含高密度组织区域的分离CT图像；
重投影模块200，用于使用平行光源在预定的角度范围内分别对原始CT图像和分离CT图像进行重投影获得原始CT图像重投影数据和分离CT图像重投影数据；
重投影数据校正模块300，用于根据所述原始CT图像重投影数据设置一校正区域和重置区域，将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据；
误差CT图像重建模块400，用于对所述分离CT图像重投影校正数据进行CT图像重建获得误差CT图像；
校正处理模块500，用于利用所述误差CT图像对所述原始CT图像进行校正处理获得校正CT图像；
所述分离CT图像获取模块100、重投影模块200、重投影数据校正模块300、误差CT图像重建模块400、校正处理模块500依次连接。
进一步，如图6所示，所述重投影数据校正模块300包括：
区域设置单元310，用于设置常量
Figure 838818DEST_PATH_IMAGE008
和
Figure 507697DEST_PATH_IMAGE009
，当
Figure 94536DEST_PATH_IMAGE010
，则判定该投影数据
Figure 848865DEST_PATH_IMAGE011
为重置区域，当
Figure 585877DEST_PATH_IMAGE012
，则判定该投影数据
Figure 109262DEST_PATH_IMAGE011
为校正区域，其中，
Figure 539107DEST_PATH_IMAGE011
为原始CT图像重投影数据在角度s处的投影数据；
重投影数据校正单元320，用于将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据；
所述区域设置单元310与重投影数据校正单元320连接。
进一步，如图7所示，所述误差CT图像重建模块400包括：
滤波核设置单元410，用于设置与所述分离CT图像重投影校正数据对应的滤波核；
快速傅里叶变换单元420，用于对所述分离CT图像重投影校正数据和滤波核分别进行快速傅里叶变换获得分离CT图像重投影校正数据频域以及滤波核频域；
滤波数据获取单元430，用于将所述分离CT图像重投影校正数据频域以及滤波核频域在频域内相乘获得频域的滤波数据；
离散傅里叶变换单元440，用于对所述滤波数据进行离散傅里叶逆变换获得滤波投影数据；
误差CT图像获取单元450，用于对所述滤波投影数据进行反投影计算，获得误差CT图像；
所述滤波核设置单元410、快速傅里叶变换单元420、滤波数据获取单元430、离散傅里叶变换单元440、误差CT图像获取单元450依次连接。关于上述各模块及单元的功能作用在前面的方法中已有详述，故不再赘述。
综上所述，本发明射束硬化伪影的校正方法及装置，从原始CT图像中分离出仅含高密度组织的分离CT图像，并对这两种图像分别进行重投影，然后对重投影数据进行校正，使用校正后的重投影数据进行图像重建得到误差CT图像，利用该误差CT图像对原始CT图像进行校正，最后获得校正CT图像。通过本发明的校正，CT图像的伪影基本被消除或减弱， 实现了对原始CT图像校正的目的，并且本发明降低了数据处理量，加快了图像处理速度，从而提高了校正效率。
应当理解的是，本发明的应用不限于上述的举例，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，所有这些改进和变换都应属于本发明所附权利要求的保护范围。

Claims (10)

1.一种术中CT图像射束硬化伪影校正方法，其特征在于，包括步骤：
A、通过设置一预定的分割阀值从原始CT图像中分离提取高密度组织区域获得仅包含高密度组织区域的分离CT图像；
B、使用平行光源在预定的角度范围内分别对原始CT图像和分离CT图像进行重投影获得原始CT图像重投影数据和分离CT图像重投影数据；
C、根据所述原始CT图像重投影数据设置一校正区域和重置区域，将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据；
D、对所述分离CT图像重投影校正数据进行CT图像重建获得误差CT图像；
E、利用所述误差CT图像对所述原始CT图像进行校正处理获得校正CT图像。
2.根据权利要求1所述术中CT图像射束硬化伪影校正方法，其特征在于，所述步骤A的分割阀值大于1000Hu，当所述原始CT图像的目标区域的CT值大于所述分割阀值时，则该目标区域为高密度组织区域。
3.根据权利要求1所述术中CT图像射束硬化伪影校正方法，其特征在于，所述步骤B的角度范围为                                               
Figure 2012102295380100001DEST_PATH_IMAGE002
，平行光源的角度增量
Figure 2012102295380100001DEST_PATH_IMAGE004
，所述平行光源共有256条射线，所述分离CT图像重投影数据为
Figure 2012102295380100001DEST_PATH_IMAGE006
，
Figure DEST_PATH_IMAGE008
为当前射线
Figure DEST_PATH_IMAGE010
入射到原始CT图像的像素，
Figure DEST_PATH_IMAGE012
为平行光源的第i个重投影角度，
Figure DEST_PATH_IMAGE014
为常量。
4.根据权利要求3所述术中CT图像射束硬化伪影校正方法，其特征在于，所述步骤C具体包括：
C1、设置常量和
Figure DEST_PATH_IMAGE018
，当
Figure DEST_PATH_IMAGE020
，则判定该投影数据
Figure DEST_PATH_IMAGE022
为重置区域，当
Figure DEST_PATH_IMAGE024
，则判定该投影数据
Figure 660258DEST_PATH_IMAGE022
为校正区域，其中，
Figure 2116DEST_PATH_IMAGE022
为原始CT图像重投影数据在角度s处的投影数据；
C2、将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据。
5.根据权利要求1所述术中CT图像射束硬化伪影校正方法，其特征在于，所述步骤D具体包括：
D1、设置与所述分离CT图像重投影校正数据对应的滤波核；
D2、对所述分离CT图像重投影校正数据和滤波核分别进行快速傅里叶变换获得分离CT图像重投影校正数据频域以及滤波核频域；
D3、将所述分离CT图像重投影校正数据频域以及滤波核频域在频域内相乘获得频域的滤波数据；
D4、对所述滤波数据进行离散傅里叶逆变换获得滤波投影数据；
D5、对所述滤波投影数据进行反投影计算，获得误差CT图像。
6.根据权利要求1所述术中CT图像射束硬化伪影校正方法，其特征在于，所述步骤E具体包括：
将原始CT图像与误差CT图像中的每一像素进行逐行累加，获得校正CT图像。
7.根据权利要求1所述术中CT图像射束硬化伪影校正方法，其特征在于，所述步骤A之前还包括：
利用探测器探测多个视角的投影帧数据，并对所述投影帧数据进行切片重建获取重建切片的原始CT图像。
8.一种术中CT图像射束硬化伪影校正装置，其特征在于，包括：
分离CT图像获取模块，用于通过设置一预定的分割阀值从原始CT图像中分离提取高密度组织区域获得仅包含高密度组织区域的分离CT图像；
重投影模块，用于使用平行光源在预定的角度范围内分别对原始CT图像和分离CT图像进行重投影获得原始CT图像重投影数据和分离CT图像重投影数据；
重投影数据校正模块，用于根据所述原始CT图像重投影数据设置一校正区域和重置区域，将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据；
误差CT图像重建模块，用于对所述分离CT图像重投影校正数据进行CT图像重建获得误差CT图像；
校正处理模块，用于利用所述误差CT图像对所述原始CT图像进行校正处理获得校正CT图像；
所述分离CT图像获取模块、重投影模块、重投影数据校正模块、误差CT图像重建模块、校正处理模块依次连接。
9.根据权利要求8所述术中CT图像射束硬化伪影校正装置，其特征在于，所述重投影数据校正模块包括：
区域设置单元，用于设置常量
Figure 35931DEST_PATH_IMAGE016
和，当
Figure 235892DEST_PATH_IMAGE020
，则判定该投影数据
Figure 561700DEST_PATH_IMAGE022
为重置区域，当
Figure 142198DEST_PATH_IMAGE024
，则判定该投影数据
Figure 267280DEST_PATH_IMAGE022
为校正区域，其中，
Figure 524709DEST_PATH_IMAGE022
为原始CT图像重投影数据在角度s处的投影数据；
重投影数据校正单元，用于将所述分离CT图像重投影数据中的与校正区域相应的区域的CT值替换为预定常量，并将所述分离CT图像重投影数据中的重置区域相应的区域的CT值重置为0，获得分离CT图像重投影校正数据；
所述区域设置单元与重投影数据校正单元相连接。
10.根据权利要求8所述术中CT图像射束硬化伪影校正装置，其特征在于，所述误差CT图像重建模块包括：
滤波核设置单元，用于设置与所述分离CT图像重投影校正数据对应的滤波核；
快速傅里叶变换单元，用于对所述分离CT图像重投影校正数据和滤波核分别进行快速傅里叶变换获得分离CT图像重投影校正数据频域以及滤波核频域；
滤波数据获取单元，用于将所述分离CT图像重投影校正数据频域以及滤波核频域在频域内相乘获得频域的滤波数据；
离散傅里叶变换单元，用于对所述滤波数据进行离散傅里叶逆变换获得滤波投影数据；
误差CT图像获取单元，用于对所述滤波投影数据进行反投影计算，获得误差CT图像；
所述滤波核设置单元、快速傅里叶变换单元、滤波数据获取单元、离散傅里叶变换单元、误差CT图像获取单元依次连接。
CN201210229538.0A 2012-07-04 2012-07-04 一种术中ct图像射束硬化伪影校正方法及装置 Active CN102768759B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210229538.0A CN102768759B (zh) 2012-07-04 2012-07-04 一种术中ct图像射束硬化伪影校正方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210229538.0A CN102768759B (zh) 2012-07-04 2012-07-04 一种术中ct图像射束硬化伪影校正方法及装置

Publications (2)

Publication Number Publication Date
CN102768759A true CN102768759A (zh) 2012-11-07
CN102768759B CN102768759B (zh) 2014-11-26

Family

ID=47096152

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210229538.0A Active CN102768759B (zh) 2012-07-04 2012-07-04 一种术中ct图像射束硬化伪影校正方法及装置

Country Status (1)

Country Link
CN (1) CN102768759B (zh)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104700377A (zh) * 2013-12-06 2015-06-10 Ge医疗系统环球技术有限公司 获得对计算机断层扫描数据进行射束硬化校正的射束硬化校正系数的方法和装置
CN105184835A (zh) * 2015-09-15 2015-12-23 上海联影医疗科技有限公司 乳腺断层图像重建方法和装置
CN105374014A (zh) * 2015-08-25 2016-03-02 上海联影医疗科技有限公司 图像校正方法及装置、医学图像生成方法及装置
CN105374012A (zh) * 2014-08-27 2016-03-02 通用电气公司 用于消除由性能差异的探测器单元所致的条状伪影的方法
CN105528771A (zh) * 2016-01-19 2016-04-27 南京邮电大学 一种使用能量函数方法的锥束ct中杯状伪影的校正方法
CN105989580A (zh) * 2015-02-28 2016-10-05 朱全祥 一种ct成像校正方法
CN106485680A (zh) * 2016-10-13 2017-03-08 上海联影医疗科技有限公司 图像校正方法及装置
CN106551703A (zh) * 2015-09-30 2017-04-05 上海联影医疗科技有限公司 计算机断层成像方法和计算机断层成像系统
CN106551702A (zh) * 2015-09-30 2017-04-05 通用电气公司 Ct图像的射束硬化伪影校正装置及方法
WO2017128892A1 (zh) * 2016-01-30 2017-08-03 上海联影医疗科技有限公司 计算机断层成像伪影校正方法及系统
CN107374657A (zh) * 2017-06-30 2017-11-24 上海联影医疗科技有限公司 对ct扫描数据进行校正的方法及ct扫描系统
CN107430779A (zh) * 2015-03-09 2017-12-01 皇家飞利浦有限公司 多能量（谱）图像数据处理
US9875558B2 (en) 2015-09-15 2018-01-23 Shanghai United Imaging Healthcare Co., Ltd. Image reconstruction system and method
US10586355B2 (en) 2015-09-15 2020-03-10 Shanghai United Imaging Healthcare Co., Ltd. Image reconstruction system and method
US10722204B2 (en) 2015-09-18 2020-07-28 Shanghai United Imaging Healthcare Co., Ltd. System and method for computer tomography
CN112233156A (zh) * 2020-10-14 2021-01-15 首都师范大学 微纳米ct投影数据的中心切片对齐方法
US10977843B2 (en) 2017-06-28 2021-04-13 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for determining parameters for medical image processing
EP4047344A1 (en) 2021-02-19 2022-08-24 Helmholtz-Zentrum Dresden - Rossendorf e.V. A method and device for 3d analyzes of a particulate material

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101698033B1 (ko) * 2016-02-26 2017-01-19 연세대학교 산학협력단 Ct 영상에서의 콘빔 아티팩트 보정 장치 및 방법

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101404088A (zh) * 2008-11-05 2009-04-08 华中科技大学 Ct图像重建的方法及系统
CN102016923A (zh) * 2008-05-06 2011-04-13 皇家飞利浦电子股份有限公司 图像伪影的减少
US20110206258A1 (en) * 2010-02-23 2011-08-25 Shoupu Chen Method and system for cone beam computed tomography high density object artifact reduction

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102016923A (zh) * 2008-05-06 2011-04-13 皇家飞利浦电子股份有限公司 图像伪影的减少
CN101404088A (zh) * 2008-11-05 2009-04-08 华中科技大学 Ct图像重建的方法及系统
US20110206258A1 (en) * 2010-02-23 2011-08-25 Shoupu Chen Method and system for cone beam computed tomography high density object artifact reduction

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李印生等: "基于自适应方向性滤波和非局部均值修补的CT图像金属伪影消除", 《中国生物医学工程学报》 *
魏英等: "CT图像"条状"伪影校正方法研究", 《光学技术》 *

Cited By (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104700377A (zh) * 2013-12-06 2015-06-10 Ge医疗系统环球技术有限公司 获得对计算机断层扫描数据进行射束硬化校正的射束硬化校正系数的方法和装置
CN104700377B (zh) * 2013-12-06 2019-07-30 Ge医疗系统环球技术有限公司 获得对计算机断层扫描数据进行射束硬化校正的射束硬化校正系数的方法和装置
US10217248B2 (en) 2014-08-27 2019-02-26 General Electric Company Method for removing streak from detector cell with performance difference
CN105374012A (zh) * 2014-08-27 2016-03-02 通用电气公司 用于消除由性能差异的探测器单元所致的条状伪影的方法
CN105374012B (zh) * 2014-08-27 2018-11-27 通用电气公司 用于消除由性能差异的探测器单元所致的条状伪影的方法
CN105989580A (zh) * 2015-02-28 2016-10-05 朱全祥 一种ct成像校正方法
CN107430779A (zh) * 2015-03-09 2017-12-01 皇家飞利浦有限公司 多能量（谱）图像数据处理
CN107430779B (zh) * 2015-03-09 2021-08-24 皇家飞利浦有限公司 多能量（谱）图像数据处理
CN105374014B (zh) * 2015-08-25 2018-01-12 上海联影医疗科技有限公司 图像校正方法及装置、医学图像生成方法及装置
CN105374014A (zh) * 2015-08-25 2016-03-02 上海联影医疗科技有限公司 图像校正方法及装置、医学图像生成方法及装置
US10600214B2 (en) 2015-09-15 2020-03-24 Shanghai United Imaging Healthcare Co., Ltd. Image reconstruction system and method
US10586355B2 (en) 2015-09-15 2020-03-10 Shanghai United Imaging Healthcare Co., Ltd. Image reconstruction system and method
US11335041B2 (en) 2015-09-15 2022-05-17 Shanghai United Imaging Healthcare Co., Ltd. Image reconstruction system and method
US9875558B2 (en) 2015-09-15 2018-01-23 Shanghai United Imaging Healthcare Co., Ltd. Image reconstruction system and method
CN105184835A (zh) * 2015-09-15 2015-12-23 上海联影医疗科技有限公司 乳腺断层图像重建方法和装置
CN105184835B (zh) * 2015-09-15 2018-11-06 上海联影医疗科技有限公司 乳腺断层图像重建方法和装置
US11191509B2 (en) 2015-09-18 2021-12-07 Shanghai United Imaging Healthcare Co., Ltd. System and method for computed tomography
US10722204B2 (en) 2015-09-18 2020-07-28 Shanghai United Imaging Healthcare Co., Ltd. System and method for computer tomography
CN106551702A (zh) * 2015-09-30 2017-04-05 通用电气公司 Ct图像的射束硬化伪影校正装置及方法
CN106551702B (zh) * 2015-09-30 2021-05-14 通用电气公司 Ct图像的射束硬化伪影校正装置及方法
CN106551703B (zh) * 2015-09-30 2018-10-30 上海联影医疗科技有限公司 计算机断层成像方法和计算机断层成像系统
CN106551703A (zh) * 2015-09-30 2017-04-05 上海联影医疗科技有限公司 计算机断层成像方法和计算机断层成像系统
CN105528771B (zh) * 2016-01-19 2018-09-18 南京邮电大学 一种使用能量函数方法的锥束ct中杯状伪影的校正方法
CN105528771A (zh) * 2016-01-19 2016-04-27 南京邮电大学 一种使用能量函数方法的锥束ct中杯状伪影的校正方法
US10540766B2 (en) 2016-01-30 2020-01-21 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for artifact correction of computerized tomography
WO2017128892A1 (zh) * 2016-01-30 2017-08-03 上海联影医疗科技有限公司 计算机断层成像伪影校正方法及系统
CN106485680B (zh) * 2016-10-13 2017-11-28 上海联影医疗科技有限公司 图像校正方法及装置
CN106485680A (zh) * 2016-10-13 2017-03-08 上海联影医疗科技有限公司 图像校正方法及装置
US10977843B2 (en) 2017-06-28 2021-04-13 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for determining parameters for medical image processing
US11908046B2 (en) 2017-06-28 2024-02-20 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for determining processing parameter for medical image processing
CN107374657A (zh) * 2017-06-30 2017-11-24 上海联影医疗科技有限公司 对ct扫描数据进行校正的方法及ct扫描系统
CN112233156A (zh) * 2020-10-14 2021-01-15 首都师范大学 微纳米ct投影数据的中心切片对齐方法
CN112233156B (zh) * 2020-10-14 2022-02-15 俐玛精密测量技术(苏州)有限公司 微纳米ct投影数据的中心切片对齐方法
EP4047344A1 (en) 2021-02-19 2022-08-24 Helmholtz-Zentrum Dresden - Rossendorf e.V. A method and device for 3d analyzes of a particulate material

Also Published As

Publication number Publication date
CN102768759B (zh) 2014-11-26

Similar Documents

Publication Publication Date Title
CN102768759B (zh) 一种术中ct图像射束硬化伪影校正方法及装置
US8111894B2 (en) Computer Tomography (CT) C-arm system and method for examination of an object
US8818065B2 (en) Methods and apparatus for scatter correction for CBCT system and cone-beam image reconstruction
JP6492005B2 (ja) Ｘ線ｃｔ装置、再構成演算装置、及び再構成演算方法
US8520974B2 (en) Image processing method, an apparatus therefor and a tomographic apparatus for removing artifacts from a sectional image
EP2283464B1 (en) Method and apparatus for scatter correction
US9008274B2 (en) Systems and methods for selecting image display parameters
CN107260195B (zh) 用于计算机断层扫描成像的伪影去除的系统和方法
Meilinger et al. Metal artifact reduction in cone beam computed tomography using forward projected reconstruction information
US8553959B2 (en) Method and apparatus for correcting multi-modality imaging data
US20060291711A1 (en) Imaging chain for digital tomosynthesis on a flat panel detector
DE10356116A1 (de) Verfahren und Vorrichtung zum Erleichtern eines Verringerns von Artefakten
Zou et al. Imaging cochlear implantation with round window insertion in human temporal bones and cochlear morphological variation using high-resolution cone beam CT
US20120163557A1 (en) Methods and phantoms for calibrating an imaging system
US20170084059A1 (en) Image Generating Apparatus, Radiation Tomography Imaging Apparatus, and Image Generating Method and Program
JP5610474B2 (ja) 画像処理装置およびプログラム並びに画像診断装置
US10552992B2 (en) Poly-energetic reconstruction method for metal artifacts reduction
Seeram et al. Computed tomography: Physical principles, instrumentation, and quality control
JP2010075443A (ja) 断層像処理装置、ｘ線ｃｔ装置およびプログラム
US20120177173A1 (en) Method and apparatus for reducing imaging artifacts
CN103829964A (zh) X射线检测器、准直仪、ct设备及其所用的方法
EP2864964B1 (en) Image reconstruction in interleaved multi-energy imaging
US8620052B2 (en) Projection truncation processing for CBCT
JP2011036684A (ja) コンピュータ支援画像診断システム
Xia et al. Scaling calibration in the ATRACT algorithm

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Intraoperative CT (Computed Tomography) image beam hardening artifact correction method and device

Effective date of registration: 20160301

Granted publication date: 20141126

Pledgee: Shenzhen SME financing Company limited by guarantee

Pledgor: Anke High-Tech Co., Ltd., Shenzhen City

Registration number: 2016990000154

PLDC Enforcement, change and cancellation of contracts on pledge of patent right or utility model
PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20170330

Granted publication date: 20141126

Pledgee: Shenzhen SME financing Company limited by guarantee

Pledgor: Anke High-Tech Co., Ltd., Shenzhen City

Registration number: 2016990000154

PC01 Cancellation of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Intraoperative CT (Computed Tomography) image beam hardening artifact correction method and device

Effective date of registration: 20170419

Granted publication date: 20141126

Pledgee: Shenzhen SME financing Company limited by guarantee

Pledgor: Anke High-Tech Co., Ltd., Shenzhen City

Registration number: 2017990000306

PE01 Entry into force of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20180522

Granted publication date: 20141126

Pledgee: Shenzhen SME financing Company limited by guarantee

Pledgor: Anke High-Tech Co., Ltd., Shenzhen City

Registration number: 2017990000306

PC01 Cancellation of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Intraoperative CT (Computed Tomography) image beam hardening artifact correction method and device

Effective date of registration: 20180626

Granted publication date: 20141126

Pledgee: Shenzhen SME financing Company limited by guarantee

Pledgor: Anke High-Tech Co., Ltd., Shenzhen City

Registration number: 2018990000490

PE01 Entry into force of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20190903

Granted publication date: 20141126

Pledgee: Shenzhen SME financing Company limited by guarantee

Pledgor: Anke High-Tech Co., Ltd., Shenzhen City

Registration number: 2018990000490

PC01 Cancellation of the registration of the contract for pledge of patent right
CP02 Change in the address of a patent holder

Address after: No.2 workshop, Lingya Industrial Park, No.1 Tangtou Road, Tangtou community, Shiyan street, Bao'an District, Shenzhen City, Guangdong Province

Patentee after: SHENZHEN ANKE HIGH-TECH Co.,Ltd.

Address before: 518076, No. 26, Mount Shun Road, Shekou, Guangdong, Shenzhen, Nanshan District

Patentee before: SHENZHEN ANKE HIGH-TECH Co.,Ltd.

CP02 Change in the address of a patent holder