CN106651753A - 提高ct图像显示效果的方法及装置 - Google Patents
提高ct图像显示效果的方法及装置 Download PDFInfo
- Publication number
- CN106651753A CN106651753A CN201610859978.2A CN201610859978A CN106651753A CN 106651753 A CN106651753 A CN 106651753A CN 201610859978 A CN201610859978 A CN 201610859978A CN 106651753 A CN106651753 A CN 106651753A
- Authority
- CN
- China
- Prior art keywords
- image
- img
- convolution kernel
- ker
- frequency
- 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
Links
- 230000000694 effects Effects 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000002591 computed tomography Methods 0.000 claims abstract description 130
- 239000011159 matrix material Substances 0.000 claims abstract description 33
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 8
- 230000000007 visual effect Effects 0.000 claims description 5
- 230000006870 function Effects 0.000 description 16
- 238000005070 sampling Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 210000000988 bone and bone Anatomy 0.000 description 4
- 210000004872 soft tissue Anatomy 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
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/08—Projecting images onto non-planar surfaces, e.g. geodetic screens
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/92—Dynamic range modification of images or parts thereof based on global image properties
-
- 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本申请提供一种提高CT图像显示效果的方法及装置,该方法包括:对被检测体进行CT扫描,获取被测体的CT投影数据;根据预设的用于重建CT图像的CT图像矩阵的大小Sizeimg和建像视野FOV,计算图像矩阵所支持的图像截止频率Fimg;读入预设的卷积核和所述卷积核的截止频率Fker;当Fker<Fimg时,将所述卷积核作为重建卷积核;当Fker≥Fimg时,调节Fker至预设经验值k乘以Fimg的大小,并截取所述卷积核中频率小于等于调节后Fker的部分作为重建卷积核;根据CT投影数据和重建卷积核,利用卷积反投影算法重建CT图像。本申请能够避免重建的CT图像中出现混叠伪影,提高CT图像的显示效果。
Description
技术领域
本申请涉及医学影像处理技术领域,尤其涉及一种提高CT图像显示效果的方法及装置。
背景技术
CT(Computed Tomography,即电子计算机断层扫描)图像可由CT机采集到的投影数据与预先设定的卷积核反投影计算获得,使用不同的卷积核会得到不同的建像效果(例如突出显示骨组织、软组织等),以满足不同的临床需求。然而,CT图像显示效果的好坏,与CT图像矩阵所支持的截止频率Fimg和图像重建所使用的卷积核截止频率Fker有关。
通常,CT机的采样频率Fs是固定的,即CT机中检测器的采样频率的大小一般是固定的,假设某次扫描数据的最高频率为Fc,根据采样定理,只有当Fs≥2Fc时,重建出的CT图像才会无混叠。而某次扫描数据的最高频率Fc与CT图像矩阵所支持的截止频率Fimg和图像重建所使用的卷积核截止频率Fker有关,当最高频率Fc不满足采样定理时,重建图像中会有混叠伪影。
目前,在利用卷积核和投影数据重建CT图像时,不会对卷积核的截止频率Fker进行调整,当卷积核的截止频率Fker过大时,会造成CT图像的混叠伪影,重建的CT图像质量较差,可能导致病灶的误判断。
发明内容
有鉴于此,本申请提供一种提高CT图像显示效果的方法及装置,以解决目前利用反投影算法重建的CT图像存在混叠伪影的问题。
具体地,本申请是通过如下技术方案实现的:
一种提高CT图像显示效果的方法,包括:
对被检测体进行CT扫描,获取被测体的CT投影数据;
根据预设的用于重建CT图像的CT图像矩阵的大小Sizeimg和建像视野FOV,计算图像矩阵所支持的图像截止频率Fimg;
读入预设的卷积核和所述卷积核的截止频率Fker,所述卷积核包括一一对应的卷积核数据幅值和频率;
当Fker<Fimg时,将所述卷积核作为重建卷积核;
当Fker≥Fimg时,调节Fker至预设经验值k乘以Fimg的大小,并截取所述卷积核中频率小于等于调节后Fker的部分作为重建卷积核,其中,k*Fimg<Fker;
根据所述被测体的CT投影数据和所述重建卷积核,利用卷积反投影算法重建CT图像。
根据本申请提出的一种具体实施方式中,所述调节Fker至预设经验值k乘以Fimg的大小包括:
利用Fimg和预设的窗函数,截取Fker至k乘以Fimg的大小。
根据本申请提出的一种具体实施方式中,所述窗函数是矩形窗或sin窗。
根据本申请提出的一种具体实施方式中,所述Sizeimg为所述CT图像矩阵的行数或列数,所述建像视野为所述CT图像建像视野的直径大小。
根据本申请的一种提高CT图像显示效果的方法,本申请还提供一种提高CT图像显示效果的装置。
一种提高CT图像显示效果的装置,包括:
数据获取模块,根据对被测体进行CT扫描,获取被测体的CT投影数据;
图像频率计算模块,根据预设的用于重建CT图像的CT图像矩阵的大小Sizeimg和建像视野FOV,计算图像矩阵所支持的图像截止频率Fimg;
曲线预设模块,读入预设的卷积核和所述卷积核的截止频率Fker,所述卷积核包括一一对应的卷积核数据幅值和频率;
曲线调节模块,当Fker<Fimg时,将所述卷积核作为重建卷积核;当Fker≥Fimg时,调节Fker至预设经验值k乘以Fimg的大小,并截取所述卷积核中频率小于等于调节后Fker的部分作为重建卷积核,其中,k*Fimg<Fker;
图像重建模块,根据所述被测体的CT投影数据和所述重建卷积核,利用卷积反投影算法重建CT图像。
根据本申请提出的一种具体实施方式中,所述曲线调节模块进一步包括:
利用Fimg和预设的窗函数,截取Fker至k乘以Fimg的大小。
根据本申请提出的一种具体实施方式中,所述窗函数是矩形窗或sin窗。
根据本申请提出的一种具体实施方式中,所述Sizeimg为所述CT图像矩阵的行数或列数,所述建像视野为所述CT图像建像视野的直径大小。
本申请的有益效果:根据CT图像矩阵所支持的图像截止频率对重建CT图像的卷积核截止频率进行调整,使得调整后的卷积核截止频率与图像截止频率相匹配,能够有效抑制卷积核中的高频成分,从而避免重建的CT图像中出现混叠伪影,提高CT图像的显示效果。
附图说明
图1是本申请提供的一种提高CT图像显示效果的方法流程图;
图2是本申请提供的卷积核曲线示意图;
图3是对图2的卷积核曲线调整后的示意图;
图4是本申请提供的一种提高CT图像显示效果的装置的功能模块结构原理示意图;
图5是本申请提供的一种CT机的硬件结构图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。
在医学临床上,需要使用CT机的检测器采集被测体内部组织(例如骨组织、软组织等)的投影数据,再利用投影数据和卷积核反投影计算获得CT图像,最终通过CT图像判断被测体的病灶信息。
现有技术中,在利用卷积核和投影数据重建CT图像时,通常不会对卷积核的截止频率Fker进行调整,当卷积核的截止频率Fker过大时,会造成CT图像的混叠伪影,重建的CT图像质量较差,可能导致病灶的误判断。
为解决上述问题,本申请提供了一种提高CT图像显示效果的方法和装置。其中,该提高CT图像显示效果的方法和装置适用于CT机,以消除CT机采集的病灶图像(CT图像)的混叠伪影,提高CT机采集的病灶图像的质量。
如图1所示,该方法包括:
S101,对被测体进行CT扫描,获取被测体的CT投影数据。
本申请中,CT机中的检测器和射线源相对位置固定不变,检测器和射线源在旋转机构上绕旋转中心进行螺旋旋转,进而相对被测体发生位置改变,因而检测器能够检测到被测体不同位置处的CT投影数据(即CT正投影数据)。
S102,根据预设的用于重建CT图像的CT图像矩阵的大小Sizeimg和建像视野FOV,计算CT图像矩阵所支持的图像截止频率Fimg(单位:lp/cm,即线对/厘米)。
本申请中,Sizeimg为预设重建CT图像的CT图像矩阵的行数或列数,例如,预设重建CT图像需要获得m行、n列的CT图像矩阵,则Sizeimg=m或Sizeimg=n。优选地,m=n,例如m=n=512,或者,m=n=768。
建像视野FOV为该CT图像建像视野的直径大小,例如,在获取CT图像时,建像视野为一个直径D为250mm的圆,则FOV=D。
根据Sizeimg和FOV,得到该CT图像矩阵所支持的图像截止频率Fimg的计算公式如下:
在一可行实施方式中,假设CT机监测的对象为骨组织,重建的CT图像矩阵为由骨组织的CT值组成的768行、768列的矩阵重建CT图像时的建像的视野为直径250mm的圆,则由公式(1)计算出Fimg等于15.36(单位:lp/cm,即线对/厘米)。
S103,读入预设的卷积核和该卷积核的截止频率Fker(单位:lp/cm,即线对/厘米),其中,卷积核包括卷积核数据幅值(单位为1)和频率(单位:lp/cm,即线对/厘米)。其中卷积核数据幅值相当于在重建CT图像时,各频率需要相乘的系数。
参见图2,为本申请卷积核的曲线示意图,该卷积核曲线的纵坐标为卷积核数据幅值,横坐标为与卷积核数据幅值相对应的频率,由该卷积核曲线可获得截止频率Fker。
S104,比较Fker与Fimg的大小,若Fker<Fimg,则进入步骤S105;反之,若Fker≥Fimg,则进入步骤S106。
由于CT机的采样频率Fs是固定的,其某次扫描数据的最高频率Fc在符合采样定理Fs≥2Fc时,重建出的CT图像才会无混叠。而某次扫描数据的最高频率Fc与CT图像矩阵所支持的截止频率Fimg和图像重建所使用的卷积核截止频率Fker有关。
其中,Fc=max(Fimg,Fker) (1)
本实施例中,Fimg需满足:Fs≥2Fimg (2)
具体地,根据用于重建CT图像的CT图像矩阵的大小Sizeimg和建像视野FOV将Fimg的设计成满足Fs≥2Fimg的要求。
当Fker<Fimg时,说明CT图像重建所使用的卷积核的频率较低,保持Fker不变,经该卷积核调整后的被测体的CT投影数据的频率也较低。
由公式(1)可以得出当Fker<Fimg时,Fc=Fimg,再根据公式(2),得出Fc满足采样定理,采用该卷积核作为重建卷积核,重建的CT图像中不存在混叠伪影。故在Fker<Fimg时,无需对Fker的大小进行调整,即可重建出不存在混叠伪影的CT图像。
当Fker≥Fimg时,说明CT图像重建所使用的卷积核存在高频成分,且与Fimg相比,Fker越大,卷积核的高频成分越高。经该卷积核调整后的被测体的CT投影数据会保留较多的高频成分。
由公式(1)可以得出当Fker≥Fimg时,Fc=Fker,再根据公式(2),得出Fc不满足采样定理,若不对Fker的大小进行调整,重建的CT图像会存在混叠伪影。因而,在当Fker≥Fimg时,需对Fker的大小进行调整,Fc=max(Fimg,调整后的Fker),从而使得Fc满足采样定理,并将调整Fker后的卷积核作为重建卷积核,使得经卷积核调整后的CT投影数据中不满足采样定理的高频成分被有效抑制,从而保证重建出的CT图像不存在混叠伪影。
S105,将该卷积核作为重建卷积核。
本申请中,当判断Fker<Fimg时,保持步骤S103读入的卷积核和该卷积核的Fker均不变,即保持图2所示的卷积核曲线不变,将该卷积核作为重建卷积核,以保证重建图像达到Fker所支持的线束效果。
S106,调节Fker至预设经验值k乘以Fimg的大小,并截取卷积核中频率小于等于调节后Fker的部分作为重建卷积核;其中,k*Fimg<Fker。
参见图3,为对图2中读入的预设卷积核曲线进行调整后的卷积核曲线示意图,在该例子中,假设原卷积核的Fker为12.5,Fimg为10,在本申请中,调节Fker为8,并截取原卷积核曲线中频率为0-8的卷积核数据幅值,将该截取的部分作为重建卷积核,以消除重建的CT图像中的混叠伪影,且达到Fimg所支持的显示效果。
可选地,k=0.8,则调节后的Fker=0.8*Fimg。
该预设经验值k的获取方式如下:利用Fimg和预设的窗(Window)函数,截取Fker至k乘以Fimg的大小,k为经验值,且k<Fker/Fimg。
其中窗函数可根据实际情况选择,例如,将窗函数选择为矩形窗或sin窗等函数等。Fimg为需要重建的CT图像矩阵所支持的图像截止频率大小。
窗函数的频率截断长度设置为k*Fimg的大小,根据k*Fimg的大小对Fker进行截取。可选地,从原卷积核曲线的原点位置开始,将Fker截取至k*Fimg的大小。
具体地,预设函数f(Fimg,Window),该函数f(Fimg,Window)根据Fimg和窗函数计算得到截取后的频率k*Fimg。
本实施例中,将k*Fimg作为窗函数的频率截断长度,将窗函数与原卷积核曲线进行相乘,将原卷积核曲线的频率从0截断至k*Fimg的大小,获得频率从0截断至k*Fimg的卷积核数据幅值,从而获得新的卷积核作为重建卷积核。即当Fker≥Fimg时,利用函数f(Fimg,Window)将Fker截取到k*Fimg的大小,并将卷积核数据幅值截取到k*Fimg对应的位置。
S107,根据被测体的CT投影数据和重建卷积核,利用卷积反投影算法重建CT图像。
其中,若Fker<Fimg,则重建卷积核为步骤S103读入的卷积核,可保证被测体重建的CT图像达到该卷积核截止频率Fker所支持的显示效果。
若Fker≥Fimg,则重建卷积核为步骤S105截取的原卷积核中的频率小于等于调节后Fker的部分,即可防止由于Fker过大而导致重建的CT图像的混叠伪影,又可保证重建的CT图像达到图像截止频率所支持的显示效果。根据k和Fimg调节Fker后,重建的CT图像能够有效滤除混叠伪影,增强CT图像的显示质量。
本申请的提高CT图像显示效果的方法是根据CT图像矩阵所支持的图像截止频率对重建CT图像的卷积核截止频率进行调整,使得调整后的卷积核截止频率与图像截止频率相匹配,能够有效抑制卷积核中的高频成分,从而避免CT图像中出现混叠伪影,提高CT图像的显示效果。
参见图4,为本申请提供一种提高CT图像显示效果的装置的功能模块结构原理示意图,该提高CT图像显示效果的装置与上述提高CT图像显示效果的方法相对应,可参照上述提高CT图像显示效果的方法的实施例来理解或解释该提高CT图像显示效果的装置的内容。
根据应用的场景不同,所述装置有可能是通过软件实现的业务逻辑,也可能是硬件或者软硬件结合的设备。下面以软件实现为例介绍本申请装置。软件作为一个逻辑意义上的装置,是通过其所在设备的处理器将非易失性存储器中对应的计算机程序指令读取到内存中运行形成的。如图5所示,为一个例子中,本申请软件装置所在的CT机的硬件结构图。所述CT机除了包括有处理器、内存、IO接口、以及非易失性存储器之外,可能还包括其他硬件,对此不再赘述。
从功能模块来讲,本申请的提高CT图像显示效果的装置包括数据获取模块101、图像频率计算模块102、曲线预设模块103、曲线调节模块104和图像重建模块105。其中,
数据获取模块101,用以根据被测体进行CT扫描,获取CT投影数据。
图像频率计算模块102,用于根据预设的用于重建CT图像的CT图像矩阵的大小Sizeimg和建像视野FOV,计算图像矩阵所支持的图像截止频率Fimg;
其中,Sizeimg为CT图像矩阵的行数m或者列数n,FOV为CT图像建像视野的直径大小。
图像矩阵所支持的图像截止频率Fimg的计算公式如下:
曲线预设模块103,用于读入预设的卷积核和所述卷积核的截止频率Fker,其中,卷积核包括对应的卷积核数据幅值和频率;
曲线调节模块104,用于在Fker<Fimg时,保持曲线预设模块103读入预设的卷积核不变,并将曲线预设模块103读入预设的卷积核设置为重建卷积核;在Fker≥Fimg时,调节Fker至预设经验值k乘以Fimg的大小,并截取卷积核中频率小于等于调节后Fker的部分作为重建卷积核;其中,k*Fimg<Fker。
k<Fker/Fimg,以保证Fker与CT图像矩阵所支持的图像截止频率Fimg相匹配,能够有效抑制卷积核中存在的高频成分,从而起到抗混叠的作用,提高CT图像的显示效果。
可选地,曲线调节模块104进一步包括:利用Fimg和预设的窗函数,截取Fker至k乘以Fimg的大小。其中,窗函数是矩形窗或sin窗等函数。
图像重建模块105,用于根据CT投影数据和重建卷积核,利用卷积反投影算法重建CT图像。
本申请的提高CT图像显示效果的装置是根据CT图像矩阵所支持的图像截止频率Fimg对重建CT图像的卷积核截止频率Fker进行调整,使得调整后的卷积核截止频率Fker与图像截止频率Fimg相匹配,能够有效抑制卷积核中的高频成分,从而避免CT图像中出现混叠伪影,提高CT图像的显示效果。
以上所述仅为本申请的较佳实施例而已,并不用以限制本申请,凡在本申请的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本申请保护的范围之内。
Claims (8)
1.一种提高CT图像显示效果的方法,其特征在于,包括:
对被检测体进行CT扫描,获取被测体的CT投影数据;
根据预设的用于重建CT图像的CT图像矩阵的大小Sizeimg和建像视野FOV,计算图像矩阵所支持的图像截止频率Fimg;
读入预设的卷积核和所述卷积核的截止频率Fker,所述卷积核包括对应的卷积核数据幅值和频率;
当Fker<Fimg时,将所述卷积核作为重建卷积核;
当Fker≥Fimg时,调节Fker至预设经验值k乘以Fimg的大小,并截取所述卷积核中频率小于等于调节后Fker的部分作为重建卷积核,其中,k*Fimg<Fker;
根据所述被测体的CT投影数据和所述重建卷积核,利用卷积反投影算法重建CT图像。
2.如权利要求1所述的提高CT图像显示效果的方法,其特征在于,所述调节Fker至预设经验值k乘以Fimg的大小包括:
利用Fimg和预设的窗函数,截取Fker至k乘以Fimg的大小。
3.如权利要求2所述的提高CT图像显示效果的方法,其特征在于,所述窗函数是矩形窗或sin窗。
4.如权利要求1所述的提高CT图像显示效果的方法,其特征在于,所述Sizeimg为所述CT图像矩阵的行数或列数,所述建像视野为所述CT图像建像视野的直径大小。
5.一种提高CT图像显示效果的装置,其特征在于,包括:
数据获取模块,根据对被测体进行CT扫描,获取被测体的CT投影数据;
图像频率计算模块,根据预设的用于重建CT图像的CT图像矩阵的大小Sizeimg和建像视野FOV,计算图像矩阵所支持的图像截止频率Fimg;
曲线预设模块,读入预设的卷积核和所述卷积核的截止频率Fker,所述卷积核包括对应的卷积核数据幅值和频率;
曲线调节模块,当Fker<Fimg时,将所述卷积核作为重建卷积核;当Fker≥Fimg时,调节Fker至预设经验值k乘以Fimg的大小,并截取所述卷积核中频率小于等于调节后Fker的部分作为重建卷积核,其中,k*Fimg<Fker;
图像重建模块,根据所述被测体的CT投影数据和所述重建卷积核,利用卷积反投影算法重建CT图像。
6.如权利要求5所述的提高CT图像显示效果的装置,其特征在于,所述曲线调节模块进一步包括:
利用Fimg和预设的窗函数,截取Fker至k乘以Fimg的大小。
7.如权利要求6所述的提高CT图像显示效果的装置,其特征在于,所述窗函数是矩形窗或sin窗。
8.如权利要求5所述的提高CT图像显示效果的装置,其特征在于,所述Sizeimg为所述CT图像矩阵的行数或列数,所述建像视野为所述CT图像建像视野的直径大小。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610859978.2A CN106651753B (zh) | 2016-09-28 | 2016-09-28 | 提高ct图像显示效果的方法及装置 |
US15/719,051 US10181207B2 (en) | 2016-09-28 | 2017-09-28 | Reconstructing image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610859978.2A CN106651753B (zh) | 2016-09-28 | 2016-09-28 | 提高ct图像显示效果的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106651753A true CN106651753A (zh) | 2017-05-10 |
CN106651753B CN106651753B (zh) | 2020-03-17 |
Family
ID=58853462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610859978.2A Active CN106651753B (zh) | 2016-09-28 | 2016-09-28 | 提高ct图像显示效果的方法及装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US10181207B2 (zh) |
CN (1) | CN106651753B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107146211A (zh) * | 2017-06-08 | 2017-09-08 | 山东师范大学 | 基于线扩散函数和双边滤波的视网膜血管图像降噪方法 |
CN110717954A (zh) * | 2019-09-27 | 2020-01-21 | 东软医疗系统股份有限公司 | 图像重建方法及装置、电子设备、存储介质 |
CN111462273A (zh) * | 2020-05-14 | 2020-07-28 | 沈阳先进医疗设备技术孵化中心有限公司 | 图像处理方法、装置、ct设备及ct系统 |
CN116091486A (zh) * | 2023-03-01 | 2023-05-09 | 合肥联宝信息技术有限公司 | 表面缺陷检测方法、装置、电子设备及存储介质 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110430357B (zh) * | 2019-03-26 | 2021-01-29 | 华为技术有限公司 | 一种图像拍摄方法与电子设备 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5795455A (en) * | 1992-07-03 | 1998-08-18 | Institut National De La Sante Et De La Recherche Medicale (Inserm) | Method for determining the nucleotide size of DNA fragments |
US6115423A (en) * | 1997-11-19 | 2000-09-05 | Rockwell Science Center, Llc | Image coding for liquid crystal displays |
US20040252806A1 (en) * | 2003-06-16 | 2004-12-16 | Kabushiki Kaisha Toshiba | Cone-beam reconstruction apparatus and computed tomography apparatus |
US20050036582A1 (en) * | 2001-08-30 | 2005-02-17 | Takayuki Nagaoka | X-ray ct apparatus, x-ray ct apparatus afterglow correction method, and afterglow correction program |
US20050254617A1 (en) * | 2004-05-12 | 2005-11-17 | Ge Medical Systems Global Technology Company, Llc | Image reconstruction method and X-ray CT system |
CN102419866A (zh) * | 2010-09-27 | 2012-04-18 | 北京农业智能装备技术研究中心 | 用于ct图像断层重建的滤波函数建立方法及断层重建方法 |
CN102567958A (zh) * | 2010-12-31 | 2012-07-11 | 东软飞利浦医疗设备系统有限责任公司 | 一种去除ct图像中金属伪影的图像后处理方法 |
CN103654818A (zh) * | 2012-09-04 | 2014-03-26 | 沈阳东软医疗系统有限公司 | 一种实现多能量ct扫描重建的方法及装置 |
-
2016
- 2016-09-28 CN CN201610859978.2A patent/CN106651753B/zh active Active
-
2017
- 2017-09-28 US US15/719,051 patent/US10181207B2/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5795455A (en) * | 1992-07-03 | 1998-08-18 | Institut National De La Sante Et De La Recherche Medicale (Inserm) | Method for determining the nucleotide size of DNA fragments |
US6115423A (en) * | 1997-11-19 | 2000-09-05 | Rockwell Science Center, Llc | Image coding for liquid crystal displays |
US20050036582A1 (en) * | 2001-08-30 | 2005-02-17 | Takayuki Nagaoka | X-ray ct apparatus, x-ray ct apparatus afterglow correction method, and afterglow correction program |
US20040252806A1 (en) * | 2003-06-16 | 2004-12-16 | Kabushiki Kaisha Toshiba | Cone-beam reconstruction apparatus and computed tomography apparatus |
US20050254617A1 (en) * | 2004-05-12 | 2005-11-17 | Ge Medical Systems Global Technology Company, Llc | Image reconstruction method and X-ray CT system |
CN102419866A (zh) * | 2010-09-27 | 2012-04-18 | 北京农业智能装备技术研究中心 | 用于ct图像断层重建的滤波函数建立方法及断层重建方法 |
CN102567958A (zh) * | 2010-12-31 | 2012-07-11 | 东软飞利浦医疗设备系统有限责任公司 | 一种去除ct图像中金属伪影的图像后处理方法 |
CN103654818A (zh) * | 2012-09-04 | 2014-03-26 | 沈阳东软医疗系统有限公司 | 一种实现多能量ct扫描重建的方法及装置 |
Non-Patent Citations (1)
Title |
---|
傅健等: "一种用于CT空间分辨率增强的图像重建方法", 《航空动力学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107146211A (zh) * | 2017-06-08 | 2017-09-08 | 山东师范大学 | 基于线扩散函数和双边滤波的视网膜血管图像降噪方法 |
CN110717954A (zh) * | 2019-09-27 | 2020-01-21 | 东软医疗系统股份有限公司 | 图像重建方法及装置、电子设备、存储介质 |
CN110717954B (zh) * | 2019-09-27 | 2023-09-26 | 东软医疗系统股份有限公司 | 图像重建方法及装置、电子设备、存储介质 |
CN111462273A (zh) * | 2020-05-14 | 2020-07-28 | 沈阳先进医疗设备技术孵化中心有限公司 | 图像处理方法、装置、ct设备及ct系统 |
CN111462273B (zh) * | 2020-05-14 | 2024-03-08 | 东软医疗系统股份有限公司 | 图像处理方法、装置、ct设备及ct系统 |
CN116091486A (zh) * | 2023-03-01 | 2023-05-09 | 合肥联宝信息技术有限公司 | 表面缺陷检测方法、装置、电子设备及存储介质 |
CN116091486B (zh) * | 2023-03-01 | 2024-02-06 | 合肥联宝信息技术有限公司 | 表面缺陷检测方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106651753B (zh) | 2020-03-17 |
US20180089865A1 (en) | 2018-03-29 |
US10181207B2 (en) | 2019-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106651753B (zh) | 提高ct图像显示效果的方法及装置 | |
US7676073B2 (en) | System and method for reducing circular artifacts in tomographic imaging | |
JP4820582B2 (ja) | ヘリカルマルチスライスctのための回復ノイズを伴うヘリカルウィンドミルアーチファクトを低減する方法 | |
US7283605B2 (en) | Methods and apparatus for scatter correction | |
US20140161340A1 (en) | Iterative reconstruction | |
US20110286651A1 (en) | Projection-Space Denoising with Bilateral Filtering in Computed Tomography | |
JP2006015136A (ja) | 断層合成造影における直接再生方法及び装置 | |
CN103136773A (zh) | 一种稀疏角度x射线ct成像方法 | |
US9483851B2 (en) | Systems and methods for filtering for image generation | |
CN105528800A (zh) | 一种计算机断层成像伪影校正方法及装置 | |
US10722204B2 (en) | System and method for computer tomography | |
JP2014054553A (ja) | 風車アーチファクト削減用の半径内画像依存検出器列フィルタリング方法、x線ct装置及びプログラム | |
EP2803320A1 (en) | X-ray computed tomography device (x-ray ct device) using successive approximation | |
US20080002806A1 (en) | X-ray ct apparatus and x-ray ct imaging method | |
US9196064B2 (en) | Image reconstruction method and device for tilted helical scan | |
US10810766B2 (en) | Image processing apparatus, image processing method, and image processing program | |
CN111325678B (zh) | 三维体图像数据的伪影校正方法 | |
US9208540B2 (en) | Image producing method, image producing apparatus and radiation tomographic imaging apparatus, and program | |
CN110458908B (zh) | 基于有限角度迭代重建超视野ct图像的方法 | |
US20120163557A1 (en) | Methods and phantoms for calibrating an imaging system | |
JP2011239888A (ja) | X線コンピュータ断層撮影装置 | |
CN112204607B (zh) | 用于x射线成像的散射校正 | |
JP2006239118A (ja) | X線ct装置 | |
CN114298920A (zh) | 一种超视野ct图像重建模型训练和超视野ct图像重建方法 | |
Sato et al. | Aliased noise in X-ray CT images and band-limiting processing as a preventive measure |
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 |
Address after: 110167 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Applicant after: DongSoft Medical System Co., Ltd. Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16 Applicant before: Dongruan Medical Systems Co., Ltd., Shenyang |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |