CN106372035B - 一种点扩散函数估计方法及系统 - Google Patents
一种点扩散函数估计方法及系统 Download PDFInfo
- Publication number
- CN106372035B CN106372035B CN201510437211.6A CN201510437211A CN106372035B CN 106372035 B CN106372035 B CN 106372035B CN 201510437211 A CN201510437211 A CN 201510437211A CN 106372035 B CN106372035 B CN 106372035B
- Authority
- CN
- China
- Prior art keywords
- spread function
- rectangular block
- standard deviation
- point spread
- line
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000003384 imaging method Methods 0.000 claims abstract description 33
- 239000000463 material Substances 0.000 claims abstract description 32
- 238000009825 accumulation Methods 0.000 claims description 19
- 230000008859 change Effects 0.000 claims description 19
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 14
- 229910052742 iron Inorganic materials 0.000 claims description 7
- 238000003672 processing method Methods 0.000 claims description 5
- 229920000573 polyethylene Polymers 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 3
- 239000004698 Polyethylene Substances 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 229910002804 graphite Inorganic materials 0.000 claims description 3
- 239000010439 graphite Substances 0.000 claims description 3
- 239000011133 lead Substances 0.000 claims description 3
- -1 polyethylene Polymers 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 10
- 238000013316 zoning Methods 0.000 description 5
- 241000208340 Araliaceae Species 0.000 description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 235000008434 ginseng Nutrition 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000007630 basic procedure Methods 0.000 description 1
- 238000013398 bayesian method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010921 in-depth analysis Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
- G06T7/41—Analysis of texture based on statistical description of texture
- G06T7/44—Analysis of texture based on statistical description of texture using image operators, e.g. filters, edge density metrics or local histograms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/77—Determining position or orientation of objects or cameras using statistical methods
-
- 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/10116—X-ray image
-
- 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/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
-
- 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/20—Special algorithmic details
- G06T2207/20076—Probabilistic image processing
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30168—Image quality inspection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Probability & Statistics with Applications (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Image Processing (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明提供一种点扩散函数估计方法和系统,所述方法包括:对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像,所述多个矩形块的制作材料不同和/或质量厚度不同,其中扫描成像时射线的入射角度与扫描方向垂直,且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面;根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数,并分别得到所述线扩散函数的标准差参数;将两个方向的标准差参数结合,得到二维点扩散函数的参数,以估计所述点扩散函数。本发明自动化程度高,执行速度快。
Description
技术领域
本发明涉及透视成像技术领域,尤其涉及一种点扩散函数估计方法及系统。
背景技术
在透视成像系统中,经常会出现由于扫描装置运动和探测器串扰引起的图像模糊情况,对图像视觉效果和检查人员观察图像内容及发现图像中的可疑区域造成一定的影响。系统的图像质量退化可以近似为一个点扩散函数(PSF,Point Spreading Function)对图像的卷积,如果系统点扩散函数已知,可以通过一些图像复原技术对图像进行去模糊处理,从而得到清晰的图像。另外,点扩散函数的定量描述本身可以作为评价系统成像质量的一个重要指标,通过对它的深入分析还可以为系统软件硬件的改进提供指导意见。
为了改善图像质量和分析图像退化原因,需要估计系统的点扩散函数。但在以往的估计方法中,大多数是针对可见光成像系统的,如典型倒谱法、变分贝叶斯法、稀疏约束最优化方法等。这些方法都针对一般图像,没有利用特定装置,估计精度较差,且大多数方法计算速度非常慢。并且针对X-射线成像系统的方法往往需要精确已知测试装置和系统的一些参数,比如缝的宽度、孔的半径、探测器的物理尺寸等等。如果测试装置设计不合理或加工精度不高,容易造成估计不准确,计算过程也比较复杂。
发明内容
本发明提供一种点扩散函数估计方法及系统,以解决现有技术中点扩散函数估计方法自动化程度不高、执行速度慢的技术问题。
为解决上述技术问题,本发明提供一种点扩散函数估计方法,包括:
对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像,所述多个矩形块的制作材料不同和/或质量厚度不同,其中扫描成像时射线的入射角度与扫描方向垂直,且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面;
根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数,并分别得到所述线扩散函数的标准差参数;
将两个方向的标准差参数结合,得到二维点扩散函数的参数,以估计所述点扩散函数。
进一步地,所述矩形块的边长不小于所述扫描装置的探测器尺寸的预定倍数。
进一步地,所述根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数还包括:
根据所述矩形块的个数和排布方式确定各个矩形块之间边界的位置,对两个方向的边界分别计算累积灰度变化过渡曲线,并对所述累积灰度变化过渡曲线进行差分得到所述两个方向的线扩散函数。
进一步地,所述分别得到所述线扩散函数的标准差参数包括:
对所述两个方向的线扩散函数进行高斯拟合别得到所述线扩散函数的标准差参数。
进一步地,在所述分别得到所述线扩散函数的标准差参数之后,所述将两个方向的所述标准差参数结合之前,还包括:
对同一方向的多个线扩散函数的标准差参数进行平均,并结合所述矩形块的材料和质量厚度,获取此方向的平均标准差参数。
另一方面,本发明还提供一种点扩散函数估计系统,包括:
成像单元,用于对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像,所述多个矩形块的制作材料不同和/或质量厚度不同,其中扫描成像时射线的入射角度与扫描方向垂直,且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面;
参数获取单元,用于根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数,并分别得到所述线扩散函数的标准差参数;
估计单元,用于将两个方向的标准差参数结合,得到二维点扩散函数的参数,以估计所述点扩散函数。
进一步地,所述矩形块的边长不小于所述扫描装置的探测器尺寸的预定倍数。
进一步地,所述参数获取单元还用于:
根据所述矩形块的个数和排布方式确定各个矩形块之间边界的位置,对两个方向的边界分别计算累积灰度变化过渡曲线,并对所述累积灰度变化过渡曲线进行差分得到所述两个方向的线扩散函数。
进一步地,所述参数获取单元还用于:
对所述两个方向的线扩散函数进行高斯拟合别得到所述线扩散函数的标准差参数。
进一步地,所述系统还包括:
平均单元,与所述参数获取单元和所述估计单元分别相连,用于对同一方向的多个线扩散函数的标准差参数进行平均,并结合所述矩形块的材料和质量厚度,获取此方向的平均标准差参数。
再一方面,本发明还提供一种用于X射线成像装置的图像处理方法,包括如上任一项所述的点扩散函数估计方法的步骤。
另一方面,本发明还提供一种用于X射线成像装置的图像处理系统,其采用如上所述的用于X射线成像装置的图像处理方法。
再一方面,本发明还提供一种用于X射线成像装置的图像处理系统,包括辐射源、探测器和矩形块以及如权利要求上任一项所述的点扩散函数估计系统。
进一步地,所述矩形块的材料为石墨、铅、铁、聚乙烯中的一种或多种。
进一步地,所述矩形块的质量厚度为2-50g/cm2。
进一步地,所述矩形块的长和/或宽为所述探测器有效收集X射线尺寸的5-100倍。
可见,在本发明所提供的点扩散函数估计方法和系统中,能够在无需已知矩形块质量厚度和大小,也无需已知探测器尺寸的情况下估计扫描装置的点扩散函数,可以提供一种自动、快速的扫描装置的点扩散函数估计方式。本发明利用估计得到的系统点扩散函数的定量描述,可以对扫描装置的图像成像质量给出一个定量的评价,通过图像复原技术以对图像进行去模糊处理,从而提高图像的清晰度,并可根据点扩散函数的参数对扫描装置成像的硬件改进提出相应的建议。本发明适用于各种扫描成像系统的点扩散函数估计尤其适用于可近似为高斯形状、但双向标准差参数不同的点扩散函数的估计。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例点扩散函数估计方法的基本流程示意图;
图2是本发明实施例1点扩散函数估计方法的一个优选实施例流程示意图;
图3是本发明实施例1的矩形块扫描图像示意图;
图4是本发明实施例1中步骤202的具体执行步骤示意图;
图5是本发明实施例1中竖直方向累积灰度变化过渡曲线的计算区域示意图;
图6是本发明实施例1中水平方向累积灰度变化过渡曲线的计算区域示意图;
图7是本发明实施例1中图5所示区域经计算得到的累积灰度变化过渡曲线;
图8是本发明实施例1中图7所示累积灰度变化过渡曲线经差分得到的线扩散函数曲线示意图;
图9是本发明实施例1中图8所示的线扩散函数曲线经去除噪声干扰后的曲线(实线)及高斯拟合结果(虚线);
图10是本发明实施例2的点扩散函数估计方法的一个优选实施例流程示意图;
图11是本发明实施例3中点扩散函数估计系统的基本结构示意图;
图12是本发明实施例6中用于X射线成像装置的图像处理系统的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例首先提供一种点扩散函数估计方法,参见图1,包括:
步骤101:对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像,所述多个矩形块的制作材料不同和/或质量厚度不同,其中扫描成像时射线的入射角度与扫描方向垂直,且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面;
步骤102:根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数,并分别得到所述线扩散函数的标准差参数;
步骤103:将两个方向的标准差参数结合,得到二维点扩散函数的参数,以估计所述点扩散函数。
可选地,所述矩形块的边长不小于所述扫描装置的探测器尺寸的预定倍数,如矩形块的边长为扫描装置的探测器尺寸的5-100倍。其中,探测器尺寸指探测器有效收集X射线的尺寸。
可选地,根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数(LSF,Line Spreading Function)还可以包括:根据所述矩形块的个数和排布方式确定各个矩形块之间边界的位置,对两个方向的边界分别计算累积灰度变化过渡曲线,并对所述累积灰度变化过渡曲线进行差分得到所述两个方向的线扩散函数。
可选地,分别得到所述线扩散函数的标准差参数可以包括:对所述两个方向的线扩散函数进行高斯拟合别得到所述线扩散函数的标准差参数。
可选地,在所述分别得到所述线扩散函数的标准差参数之后,所述将两个方向的所述标准差参数结合之前,还可以包括:对同一方向的多个线扩散函数的标准差参数进行平均,并结合所述矩形块的材料和质量厚度,获取此方向的平均标准差参数。
实施例1:
本发明实施例1提供了利用多块大小基本相同、质量厚度不同的均匀板状材料的矩形块进行点扩散函数估计的方法,具体参见图2,包括:
步骤201:对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像。
本步骤中,利用由不同质量厚度的铁和聚乙烯材料块组成的多个矩形块进行扫描成像,扫描图像如图3所示。其中各材料块大小均为200mm×200mm,第1排为质量厚度分别为6g/cm2、8g/cm2、10g/cm2的铁材料;第2排第1个为质量厚度为4g/cm2的铁材料、第2、3个分别为质量厚度为2g/cm2和4g/cm2的聚乙烯材料;第3排第1个为质量厚度为3g/cm2的铁材料、第2、3个分别为质量厚度为6g/cm2和8g/cm2的聚乙烯材料。将这一经紧密排列固定的材料放置在表面垂直于射线入射方向的扫描装置中进行扫描,获得扫描图像,其中扫描方向与图3中矩形块的水平边方向平行。
步骤202:根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数,并分别得到所述线扩散函数的标准差参数。
本步骤的具体执行流程参见图4。其中,根据已知的上述多个矩形块每行及每列块个数信息,可在步骤401中通过简单计算得到每个矩形块边界的位置。进而,在步骤402中,在每个矩形块边界附近选择以边界中点为中心的一个矩形区域用以计算累积灰度变化过渡曲线。如选择图5所示的上排第一块和第二块之间竖直边界的计算区域,图6所示的中排第二块和下排第二块之间水平边界的计算区域,随后沿边界走向累加得到竖直和水平两个方向分别对应的累积灰度变化过渡曲线。并在步骤403中得到两个方向的线扩散函数。例如,图7表示图5所示的计算区域经纵向累加得到的累积灰度变化过渡曲线。对此曲线可以计算一阶差分从而得到近似导数,并根据灰度变化过渡曲线两端的灰度值的大小,以确定采用向前差分还是向后差分(其中当左边灰度值较小时,采用向前差分,即f(x)-f(x-1);反之,则采用向后差分,即f(x)-f(x+1)),以保证经计算得到的绝大多数差分值为正,从而得到近似线扩散函数曲线,如图8所示。
在步骤403获取两个方向的线扩散函数之后,分别从其最大值向两边搜索,遇到第一个负值或零值后,将其后所有值置零,以得到去除噪声干扰的线扩散函数曲线。其中图8中的线扩散函数曲线去除噪声干扰之后如图9中实线所示,对图9中的实线在步骤404中以高斯形状拟合,计算其均值和标准差,获取线扩散函数的标准差参数,拟合结果如图9中虚线所示。
步骤203:将两个方向的标准差参数结合,得到二维点扩散函数的参数,以估计所述点扩散函数。
本步骤中,将水平和竖直(即x和y)两个方向的标准差参数进行结合,可以估计点扩散函数PSF的分布为:
其中C为归一化系数,σx、σy分别为x方向和y方向的标准差参数。
实施例2:
参见图10,本发明实施例2中的点扩散函数估计的方法与实施例1中基本相同,其中,在步骤1004(对应本发明实施例1步骤203)之前,还包括:
步骤1003:对同一方向的多个线扩散函数的标准差参数进行平均,并结合所述矩形块的材料和质量厚度,获取此方向的平均标准差参数。
本步骤中,对本发明实施例1中所有6个内部竖直块边界和6个内部水平块边界进行前述步骤的线扩散函数计算和标准差参数拟合后,可以分别对两个方向的标准差参数进行平均,得到水平和竖直方向的平均标准差参数。
与此对应地,在步骤1004中,估计点扩散函数所采用的x方向和y方向的标准差参数σx、σy也将相应替换为x方向和y方向的平均标准差参数。
实施例3:
本发明实施例3提供了一种点扩散函数估计系统,参见图11,包括:
成像单元1101,用于对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像,所述多个矩形块的制作材料不同和/或质量厚度不同,其中扫描成像时射线的入射角度与扫描方向垂直,且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面;
参数获取单元1102,用于根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数,并分别得到所述线扩散函数的标准差参数;
估计单元1103,用于将两个方向的标准差参数结合,得到二维点扩散函数的参数,以估计所述点扩散函数。
其中,所述矩形块的边长不小于所述扫描装置的探测器尺寸的预定倍数。
其中,参数获取单元1102还可以用于:根据所述矩形块的个数和排布方式确定各个矩形块之间边界的位置,对两个方向的边界分别计算累积灰度变化过渡曲线,并对所述累积灰度变化过渡曲线进行差分得到所述两个方向的线扩散函数。
其中,参数获取单元1102还可以用于:对所述两个方向的线扩散函数进行高斯拟合别得到所述线扩散函数的标准差参数。
其中,所述系统还可以包括:平均单元,与所述参数获取单元和所述估计单元分别相连,用于对同一方向的多个线扩散函数的标准差参数进行平均,并结合所述矩形块的材料和质量厚度,获取此方向的平均标准差参数。
实施例4:
本发明实施例4通过一种用于X射线成像装置的图像处理方法,其应用了如本发明实施例所述的点扩散函数估计方法,以判断X射线成像装置的成像质量或进行图像复原。
其中,当得到所估计的点扩散函数之后,可以通过预设的像素范围阈值判断成像是否清晰,成像质量是否满足要求。当所成图像清晰时,点扩散函数的像素范围可能较小,如为3×3像素;而当成像质量较差时,点扩散函数的像素范围可能相应较大,如为7×7像素。
另外,当所获取的图像质量不尽如人意时,可以利用所得到的点扩散函数与所获取的图像进行相应运算,以复原得到质量较好的图像,提高成像质量。图像复原的运算方法可以有多种,例如可以为反卷积运算,其通过将点扩散函数与所获取的图像进行反卷积,可以得到较为理想的图像。
实施例5:
本发明实施例5提供一种利用本发明实施例4中的用于X射线成像装置的图像处理方法的用于X射线成像装置的图像处理系统。
实施例6:
本发明实施例6提供一种用于X射线成像装置的图像处理系统,参见图12,包括辐射源1201、探测器1202、矩形块1203以及如本发明实施例3中的点扩散函数估计系统1204。
其中,辐射源1201发射X射线,照射到多个紧密排列的矩形块1203上,并通过探测器1202收集信号,将所收集的信号通过点扩散函数估计系统1204进行分析,以进行点扩散函数估计。
可选地,矩形块的材料可以为石墨、铅、铁、聚乙烯中的一种或多种。
可选地,矩形块的质量厚度范围可以为2-50g/cm2。
可选地,矩形块的长和/或宽可以为探测器有效收集X射线尺寸的5-100倍。
可见,在本发明实施例所提供的点扩散函数估计方法和系统中,至少具有如下有益效果:
(1)使用简单、自动化程度高。仅要求各个材料块为厚度均匀、大小相同的矩形块,但对具体厚度和大小并没有精确要求,降低了加工要求;除扫描图像外,无需人工干预,自动完成区域定位和参数计算等过程;
(2)鲁棒性好、自适应能力强。综合利用多个区域间的阶跃边缘,受噪声影响较小;仅要求各个材料块为厚度均匀、大小基本相同的矩形块,并已知两个方向材料块的数目。但对具体厚度和大小并没有精确要求,降低了对材料厚度、尺寸等物理加工精度的要求;
(3)执行速度快。由于点扩散函数近似为高斯形状,估计过程中不需大量反卷积操作,处理速度得到了明显的提高。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (16)
1.一种点扩散函数估计方法,其特征在于,包括:
对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像,所述多个矩形块的制作材料不同和/或质量厚度不同,其中扫描成像时射线的入射角度与扫描方向垂直,且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面;
根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数,并分别得到所述线扩散函数的标准差参数;
将两个方向的标准差参数结合,得到二维点扩散函数的参数,以估计所述点扩散函数。
2.根据权利要求1所述的点扩散函数估计方法,其特征在于:
所述矩形块的边长不小于所述扫描装置的探测器尺寸的预定倍数。
3.根据权利要求1所述的点扩散函数估计方法,其特征在于,所述根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数还包括:
根据所述矩形块的个数和排布方式确定各个矩形块之间边界的位置,对两个方向的边界分别计算累积灰度变化过渡曲线,并对所述累积灰度变化过渡曲线进行差分得到所述两个方向的线扩散函数。
4.根据权利要求1所述的点扩散函数估计方法,其特征在于,所述分别得到所述线扩散函数的标准差参数包括:
对所述两个方向的线扩散函数进行高斯拟合别得到所述线扩散函数的标准差参数。
5.根据权利要求1所述的点扩散函数估计方法,其特征在于,在所述分别得到所述线扩散函数的标准差参数之后,所述将两个方向的所述标准差参数结合之前,还包括:
对同一方向的多个线扩散函数的标准差参数进行平均,并结合所述矩形块的材料和质量厚度,获取此方向的平均标准差参数。
6.一种点扩散函数估计系统,其特征在于,包括:
成像单元,用于对大小相同且紧密排列的多个矩形块利用扫描装置进行扫描成像,所述多个矩形块的制作材料不同和/或质量厚度不同,其中扫描成像时射线的入射角度与扫描方向垂直,且所述射线的入射角度垂直于所述多个矩形块的紧密排列的表面;
参数获取单元,用于根据扫描图像获取沿所述矩形块的长和宽两个方向的线扩散函数,并分别得到所述线扩散函数的标准差参数;
估计单元,用于将两个方向的标准差参数结合,得到二维点扩散函数的参数,以估计所述点扩散函数。
7.根据权利要求6所述的点扩散函数估计系统,其特征在于:
所述矩形块的边长不小于所述扫描装置的探测器尺寸的预定倍数。
8.根据权利要求6所述的点扩散函数估计系统,其特征在于,所述参数获取单元还用于:
根据所述矩形块的个数和排布方式确定各个矩形块之间边界的位置,对两个方向的边界分别计算累积灰度变化过渡曲线,并对所述累积灰度变化过渡曲线进行差分得到所述两个方向的线扩散函数。
9.根据权利要求6所述的点扩散函数估计系统,其特征在于,所述参数获取单元还用于:
对所述两个方向的线扩散函数进行高斯拟合别得到所述线扩散函数的标准差参数。
10.根据权利要求6所述的点扩散函数估计系统,其特征在于,所述系统还包括:
平均单元,与所述参数获取单元和所述估计单元分别相连,用于对同一方向的多个线扩散函数的标准差参数进行平均,并结合所述矩形块的材料和质量厚度,获取此方向的平均标准差参数。
11.一种用于X射线成像装置的图像处理方法,其特征在于,包括如权利要求1-5中任一项所述的点扩散函数估计方法的步骤。
12.一种用于X射线成像装置的图像处理系统,其特征在于,采用如权利要求11所述的用于X射线成像装置的图像处理方法。
13.一种用于X射线成像装置的图像处理系统,其特征在于,包括辐射源、探测器和矩形块以及如权利要求6-10中任一项所述的点扩散函数估计系统。
14.如权利要求13所述的用于X射线成像装置的图像处理系统,其特征在于:所述矩形块的材料为石墨、铅、铁、聚乙烯中的一种或多种。
15.如权利要求13所述的用于X射线成像装置的图像处理系统,其特征在于:所述矩形块的质量厚度为2-50g/cm2。
16.如权利要求13所述的用于X射线成像装置的图像处理系统,其特征在于:所述矩形块的长和/或宽为所述探测器有效收集X射线尺寸的5-100倍。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510437211.6A CN106372035B (zh) | 2015-07-23 | 2015-07-23 | 一种点扩散函数估计方法及系统 |
PCT/CN2016/090336 WO2017012528A1 (zh) | 2015-07-23 | 2016-07-18 | 一种点扩散函数估计方法及系统 |
US15/217,613 US10152774B2 (en) | 2015-07-23 | 2016-07-22 | Method and system for estimating point spread function |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510437211.6A CN106372035B (zh) | 2015-07-23 | 2015-07-23 | 一种点扩散函数估计方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106372035A CN106372035A (zh) | 2017-02-01 |
CN106372035B true CN106372035B (zh) | 2018-12-18 |
Family
ID=57833710
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510437211.6A Active CN106372035B (zh) | 2015-07-23 | 2015-07-23 | 一种点扩散函数估计方法及系统 |
Country Status (3)
Country | Link |
---|---|
US (1) | US10152774B2 (zh) |
CN (1) | CN106372035B (zh) |
WO (1) | WO2017012528A1 (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106372035B (zh) * | 2015-07-23 | 2018-12-18 | 清华大学 | 一种点扩散函数估计方法及系统 |
EP3680851B1 (en) * | 2019-01-09 | 2021-11-03 | Agfa Nv | Method and system to calculate the point spread function of a digital image detector system based on a mtf modulated quantum-noise measurement |
CN111488701B (zh) * | 2020-06-28 | 2020-10-16 | 中国人民解放军国防科技大学 | 一种基于点扩散函数的计算机辅助装调方法、装置及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101291391A (zh) * | 2007-04-20 | 2008-10-22 | 致伸科技股份有限公司 | 图像处理方法以及相关的部分点扩散函数估测方法 |
CN101354307A (zh) * | 2008-09-22 | 2009-01-28 | 哈尔滨工业大学 | 动态目标调制传递函数测量方法与装置 |
CN102819830A (zh) * | 2012-08-15 | 2012-12-12 | 北京交通大学 | 一种新的基于卡尔曼滤波的点扩散函数估计方法 |
CN103927740A (zh) * | 2014-02-18 | 2014-07-16 | 北京象开科技有限公司 | 获得图像的点扩展函数的方法和装置 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7251306B2 (en) * | 2004-11-17 | 2007-07-31 | General Electric Company | Methods, apparatus, and software to facilitate iterative reconstruction of images |
US8249357B2 (en) * | 2006-10-23 | 2012-08-21 | Ben Gurion University Of The Negev, Research And Development Authority | Blind restoration of images degraded by isotropic blur |
US8326011B2 (en) * | 2008-05-21 | 2012-12-04 | Varian Medical Systems, Inc. | Methods, systems, and computer-program products for estimating scattered radiation in radiographic projections |
WO2010104969A1 (en) * | 2009-03-11 | 2010-09-16 | Zoran Corporation | Estimation of point spread functions from motion-blurred images |
KR102020531B1 (ko) * | 2012-08-20 | 2019-09-10 | 삼성전자주식회사 | 선형 감마선원을 이용하여 고해상도의 pet(양전자 방출 단층 촬영) 영상을 생성하는 방법 및 장치 |
US9785860B2 (en) * | 2014-07-16 | 2017-10-10 | The Cleveland Clinic Foundation | Real-time image enhancement for X-ray imagers |
US9779491B2 (en) * | 2014-08-15 | 2017-10-03 | Nikon Corporation | Algorithm and device for image processing |
CN106372035B (zh) * | 2015-07-23 | 2018-12-18 | 清华大学 | 一种点扩散函数估计方法及系统 |
KR101850871B1 (ko) * | 2015-08-26 | 2018-04-23 | 주식회사 디알텍 | 방사선 영상의 처리방법 및 방사선 촬영시스템 |
KR20170088681A (ko) * | 2016-01-25 | 2017-08-02 | 삼성전자주식회사 | 단층 촬영 장치 및 그에 따른 단층 영상 복원 방법 |
-
2015
- 2015-07-23 CN CN201510437211.6A patent/CN106372035B/zh active Active
-
2016
- 2016-07-18 WO PCT/CN2016/090336 patent/WO2017012528A1/zh active Application Filing
- 2016-07-22 US US15/217,613 patent/US10152774B2/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101291391A (zh) * | 2007-04-20 | 2008-10-22 | 致伸科技股份有限公司 | 图像处理方法以及相关的部分点扩散函数估测方法 |
CN101354307A (zh) * | 2008-09-22 | 2009-01-28 | 哈尔滨工业大学 | 动态目标调制传递函数测量方法与装置 |
CN102819830A (zh) * | 2012-08-15 | 2012-12-12 | 北京交通大学 | 一种新的基于卡尔曼滤波的点扩散函数估计方法 |
CN103927740A (zh) * | 2014-02-18 | 2014-07-16 | 北京象开科技有限公司 | 获得图像的点扩展函数的方法和装置 |
Non-Patent Citations (3)
Title |
---|
SU Jian等.Study of the material photon and electron background and the liquid argon detector veto efficiency of the CDEX-10 experiment.《Chinese Physics C》.2015,第39卷(第3期), * |
基于点扩散函数参数辨识的运动模糊图像的盲恢复研究;廖永忠等;《光学技术》;20140531;第40卷(第3期);全文 * |
基于相位相关的匀速直线运动模糊图像位移参数估计;孙辉等;《中国光学》;20120430;第5卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
WO2017012528A1 (zh) | 2017-01-26 |
US20170024862A1 (en) | 2017-01-26 |
US10152774B2 (en) | 2018-12-11 |
CN106372035A (zh) | 2017-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10539516B2 (en) | X-ray inspection apparatus and X-ray inspection method | |
US6879715B2 (en) | Iterative X-ray scatter correction method and apparatus | |
US8548252B2 (en) | Edge detection | |
EP3558127B1 (en) | Self-calibrating ct detectors, systems and methods for self-calibration | |
CN104161536B (zh) | 一种基于互补光栅的锥束ct散射校正方法及其装置 | |
CN106372035B (zh) | 一种点扩散函数估计方法及系统 | |
CA2515564A1 (en) | Method and apparatus for removing uneven brightness in an image | |
WO2011036961A1 (ja) | X線ct装置 | |
US10282826B2 (en) | Despeckling method for radiographic images | |
WO2013146283A1 (ja) | 画像処理装置及び画像処理方法 | |
CN103886552B (zh) | 一种去除x射线影像滤线栅条纹的方法及系统 | |
KR101146006B1 (ko) | 엑스선 씨티에서의 링 아티팩트 제거방법 및 장치 | |
EP3226205B1 (en) | Image processing method and image processing apparatus | |
CN106651768B (zh) | 图像校正方法及装置、x射线摄影设备 | |
JP6685762B2 (ja) | 画像処理装置、画像処理方法、およびプログラム | |
KR20150089382A (ko) | 디지털 유방 엑스선 촬영용 시스템을 위한 검출양자효율의 검출방법, 검출시스템 및 그 방법을 실행하는 프로그램이 기록된 기록매체 | |
JP2014176565A (ja) | 画像処理装置、放射線撮影装置、画像処理方法、コンピュータプログラム及び記憶媒体 | |
Donovan et al. | Step by step analysis toward optimal MTF algorithm using an edge test device | |
Mertens et al. | A Study of Pb-Rich Dendrites in a Near-Eutectic 63Sn-37Pb Solder Microstructure via Laboratory-Scale Micro X-ray Computed Tomography (μ XCT) | |
CN109813259B (zh) | 高动态x射线成像方法、存储介质和装置 | |
KR101270002B1 (ko) | 임상용 시스템 검출 양자효율 검출 방법 및 그 시스템 | |
JP2009000164A (ja) | 放射線画像処理装置 | |
Yip et al. | Automated scoring method for the CDMAM phantom | |
CN116413290B (zh) | 一种工业ct散射校正方法 | |
JP6570594B2 (ja) | 画像処理装置、放射線撮影装置、画像処理方法、コンピュータプログラム及び記憶媒体 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |