CN105023251B - 一种高分辨锥束ct***的几何校正方法 - Google Patents
一种高分辨锥束ct***的几何校正方法 Download PDFInfo
- Publication number
- CN105023251B CN105023251B CN201510394350.5A CN201510394350A CN105023251B CN 105023251 B CN105023251 B CN 105023251B CN 201510394350 A CN201510394350 A CN 201510394350A CN 105023251 B CN105023251 B CN 105023251B
- Authority
- CN
- China
- Prior art keywords
- projected image
- geometric
- image
- cone
- parameter
- 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
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公布了一种高分辨率锥束CT的几何校正方法,包括方法采用的几何模型、校正模体、目标函数和几何参数求解算法,其特征在于所述的几何模型为六几何参数模型,校正模体为一个金属小球,目标函数的自变量为上述六个几何参数并能够反映重建图像质量的好坏,最终求解几何参数的方法采用的是最优化算法。该方法能够计算出高分辨锥束CT***的几何参数,并且计算得到的参数具有很高的精度,能够满足高分辨锥束CT***对几何校正精度的要求,能够有效的抑制高分辨锥束CT重建图像中的几何伪影,提高重建图像的质量。
Description
技术领域
本发明涉及一种在锥束CT***中对投影图像的处理方法,特别涉及一种在锥束CT***中对于投影图像空间位置偏移的校正方法,属于计算机断层成像的图像处理领域。
背景技术
在锥束CT***中,对于断层图像的重建可以获得被测物体内部图像信息,具有其它无损检测方法无法替代的特点,在工业无损检测领域具有广阔的应用前景。锥束CT***的主要组成部分包括X射线源,载物台及探测器。对于锥束投影一般采用FDK算法及其改进算法。如Practical cone-beam algorithm L.A.Feldkamp,L.C.Davis,and J.W.KressResearch Staff,Ford Motor Company,Dearborn,Michigan 48121,J.Opt.Soc.Am.A/Vol.1,No.6/June 1984中所述,标准的FDK算法要求锥束CT***中的射线源、旋转轴和平板探测器满足一个严格的空间位置关系,这个位置关系描述如下:射线源焦点与探测器中心的连线垂直于探测器平面,且该连线穿过旋转轴,同时要求旋转轴平行于平板探测器的竖直方向的边缘。
在实际的锥束CT***中,射线源、旋转轴和探测器通常无法满足上述空间位置关系,即它们的实际空间位置与理想的空间位置之间存在着一定的误差,如果不纠正这个误差,依然以标准的FDK算法对投影图像进行重建,重建图像就会发生退化,这种退化就称为几何伪影。因而要得到好的重建图像就必须对锥束CT***进行几何校正。而现有的锥束CT几何校正方法由于精度较低,用于高分辨锥束CT时并不能解决重建图像的几何伪影问题,本发明因此而来。
发明内容
针对现有技术存在的问题,本发明提出一种高分辨率锥束CT***的几何校正方法,利用几何参数与小球球心重建图像之间存在的联系,建立了一个能反映该关系的目标函数,解决了现有技术存在的精度不够高的问题。
本发明解决以上技术问题的技术方案:提供一种高分辨锥束CT***中的几何校正方法,其特征在于该方法包括以下步骤为:
根据投影图像的灰度差异选择合适的阈值对选取的不同角度投影图像进行二值化处理;
经对二值化后的投影图像中黑色像素测量得到各个投影图像的投影中心,并对各个投影图像旋转后测量得到旋转后的各个投影图像的投影中心,从而得到投影图像的旋转中心;
采用一个能反映重建体数据质量的目标函数,建立该函数与CT***几何参数之间的数学关系;
采用最优化算法计算目标函数的最大值,并最终得到几何参数的最优解;
用几何参数的最优解对高分辨锥束CT***进行几何校正,就能有效的抑制高分辨CT重建图像中的几何伪影,提高重建图像的质量。
本发明进一步限定的技术方案为:前述的高分辨锥束CT***中的几何校正方法,所述的阈值通过测量被测物体的投影灰度和背景灰度的差异来确定,即记被测物体投影灰度的平均值为A,记背景灰度的平均值为B,采用(A+B)/2作为所述阈值。
前述的高分辨锥束CT***中的几何校正方法,所述的几何参数为一组6参数的向量。
前述的高分辨锥束CT***中的几何校正方法,所述的目标函数为重建体数据的平方和。
进一步的,前述的高分辨锥束CT***中的几何校正方法,所述的最优化算法为模拟退火算法。
本发明的有益效果:本发明技术方案利用了几何参数与小球球心重建图像之间存在的联系,建立了一个能反映该关系的目标函数,并通过计算目标函数的最大值来计算得到几何参数的最优解。试验证明通过该方法计算得到的高分辨CT***的几何参数能够有效的抑制重建图像中的几何伪影,提高重建图像质量。
附图说明
图1为本发明方法采用的锥束CT的几何模型的立体示意图。
图2为本方法采用的锥束CT几何模型的俯视图。
图3为本方法采用的锥束CT几何模型的探测器平面示意图。
图4为本发明方法几何校正后的细竹签的重建图像。
图5为未采用本方法进行几何校正的细竹签的重建图像。
具体实施方式
实施例1
本实施例提供一种高精度的锥束CT几何校正方法,步骤可以分述为:
第一步,将金属小球固定于CT***的载物台上,用30kV/200μA的X射线对金属小球进行扫描,角度间隔为3.6°,扫描范围为360°,得到100张小球的投影图像。
第二步,对投影图像的灰度差异,选择合适的阈值T,对投影图像做二值化处理,设投影图像中在坐标(x,y)处的灰度值为G(x,y),如果G(x,y)≤T,置G(x,y)=0,反之若G(x,y)≥T,置G(x,y)=2N-1,其中N为投影图像的位数,这样通过二值化之后的投影图像中的目标物体为黑色,背景为白色。
第三步,对二值化后的投影图像求中心。由于物体在处理后的投影图像中为黑色,测量得出投影图像中黑色区域的中心即为投影中心。以测量横坐标为例,具体的方法是首先令nXSum=0,nCount=0,其中nXSum代表所有黑色像素的X坐标和,nCount代表黑色像素的个数,其中n表示该数据用int型表示,依次判断二值化之后的图像中的点是否为0,如果为0,nXSum加上该点的横坐标,并且黑色像素的个数nCount加1,整幅图像所有的点判断完之后,nXSum除以nCount即为该幅图像的物体的投影中心。
第四步,用数学公式描述目标函数(重建体数据像素值的平方和)的值与几何参数之间的函数关系,即给出用几何参数来计算相应的目标函数的计算公式。
第五步,选取一组几何参数作为初始值,并采用模拟退火算法迭代地计算目标函数的最大值,直至满足迭代终止条件,得到几何参数的最优解。
第六步,用得到的几何参数的最优解对投影图像进行几何校正。
具体的,本实施例以一组100张12位2048*2048像素的投影图像几何参数计算步骤,步骤如下:
1、针对投影图像的灰度差异,选择合适的阈值,对投影图像做二值化处理。
根据实际情况任意选择一块背景的区域计算出的平均值为3846,而目标物体的灰度值一般不超过2450,根据这两个值选择的阈值为T=3000。设投影图像中在坐标(x,y)处的灰度值为G(x,y),如果G(x,y)≤3000,置G(x,y)=0,反之若G(x,y)≥3000,置G(x,y)=4095。这样通过二值化之后的投影图像中的目标物体为黑色,背景为白色。
2、对二值化后的投影图像求中心。
由于物体在处理后的投影图像中为黑色,只要求出图中黑色区域的中心即为投影中心。具体的方法是首先令nXSum=0,nCount=0,依次判断二值化之后的图像中的点是否为0,如果为0,nXSum加上该点的横坐标,并且黑色像素的个数nCount加1,整幅图像所有的点判断完之后,nXSum除以nCount即为该幅图像的物体的投影中心。
根据锥束CT几何结构的特点,采用如图1中的几何模型和参数(D,θ,η,io,jo)来描述锥束CT的几何结构,并给出目标函数f(D,θ,η,io,jo)的计算公式如下:
αi=arctan(yi/xi)
采用(10000,0,0,0,1024,1024)作为几何参数的初始值,并根据上述的目标函数的计算公式,用模拟退火算法迭代的计算目标函数的最大值,并最终收敛得到几何参数的最优解(12680,5.25,3.38,-0.5375,894.875,1014.3)。
用上述的几何参数的最优解对一组400张12位2048*2048像素的细竹签投影数据进行几何校正,并对校正后的投影数据进行重建,得到校正后的重建图像如图2所示。
通过以上的实验说明采用此方法对高分辨锥束CT进行几何校正能够有效的抑制重建图像中的几何伪影,提高重建图像的质量。图2为采用此方法几何校正后的重建图像,图3为没有采用此方法的重建图像。
上述实例只为说明本发明的技术构思及特点,其目的在于让熟悉此项技术的人是能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡根据本发明精神实质所做的等效变换或修饰,都应涵盖在本发明的保护范围之内。
Claims (4)
1.一种高分辨锥束CT***中的几何校正方法,其特征在于该方法包括以下步骤为:
(1)根据投影图像的灰度差异选择合适的阈值对选取的不同角度投影图像进行二值化处理,所述阈值通过测量被测物体的投影灰度和背景灰度的差异来确定,即记被测物体投影灰度的平均值为A,记背景灰度的平均值为B,采用(A+B)/2作为所述阈值;对投影图像做二值化处理,设投影图像中在坐标(x,y)处的灰度值为G(x,y),如果G(x,y)≤T,置G(x,y)=0,反之若G(x,y)≥T,置G(x,y)=2N-1,其中N为投影图像的位数,这样通过二值化之后的投影图像中的目标物体为黑色,背景为白色;
(2)经对二值化后的投影图像中黑色像素测量得到各个投影图像的投影中心,并对各个投影图像旋转后测量得到旋转后的各个投影图像的投影中心,从而得到投影图像的旋转中心;由于物体在处理后的投影图像中为黑色,测量得出投影图像中黑色区域的中心即为投影中心,测量横坐标时,具体的方法是首先令nXSum=0,nCount=0,其中nXSum代表所有黑色像素的X坐标和,nCount代表黑色像素的个数,其中n表示nXSum、nCount数据用int型表示,依次判断二值化之后的图像中的点是否为0,如果为0,nXSum加上该点的横坐标,并且黑色像素的个数nCount加1,整幅图像所有的点判断完之后,nXSum除以nCount即为该幅图像的物体的投影中心;
(3)采用一个能反映重建体数据质量的目标函数,建立该函数与CT***几何参数之间的数学关系;目标函数的计算公式如下:
αi=arctan(yi/xi)
(4)采用最优化算法计算目标函数的最大值,并最终得到几何参数的最优解;
(5)最后通过几何参数的最优解对高分辨锥束CT***进行几何校正。
2.根据权利要求1所述的高分辨锥束CT***中的几何校正方法,其特征在于:所述的几何参数为一组6参数的向量,其中6个参数分别是D表示的是射线源到旋转轴的距离,坐标点io,jo表示的是射线源沿着穿过它并垂直于z轴的直线在探测器上的投影点;为了能够表达出探测器平面的方向偏转,引入另外三个两两正交的中间变量α,β,ω,因此θ,φ,η为各个维度上的对应角度偏转参数。
3.根据权利要求1所述的高分辨锥束CT***中的几何校正方法,其特征在于;所述的目标函数为重建体数据的平方和。
4.根据权利要求1所述的高分辨锥束CT***中的几何校正方法,其特征在于:所述的最优化算法为模拟退火算法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510394350.5A CN105023251B (zh) | 2015-07-07 | 2015-07-07 | 一种高分辨锥束ct***的几何校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510394350.5A CN105023251B (zh) | 2015-07-07 | 2015-07-07 | 一种高分辨锥束ct***的几何校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105023251A CN105023251A (zh) | 2015-11-04 |
CN105023251B true CN105023251B (zh) | 2018-08-28 |
Family
ID=54413195
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510394350.5A Active CN105023251B (zh) | 2015-07-07 | 2015-07-07 | 一种高分辨锥束ct***的几何校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105023251B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105869130B (zh) * | 2016-04-20 | 2018-02-23 | 广州华端科技有限公司 | Ct***几何校正效果的检验方法及装置 |
WO2018126335A1 (zh) * | 2017-01-03 | 2018-07-12 | 苏州海斯菲德信息科技有限公司 | 基于小球模体的锥束ct***几何参数评价及校正方法 |
CN109556542B (zh) * | 2018-11-14 | 2020-10-20 | 北京卫星制造厂有限公司 | 复杂点阵镂空结构ct尺寸测量方法 |
CN109820531B (zh) * | 2019-01-30 | 2023-04-25 | 东软医疗***股份有限公司 | 一种锥束型ct几何参数的校正方法和装置 |
CN110766629B (zh) * | 2019-10-17 | 2022-03-01 | 广州华端科技有限公司 | Cbct***几何校正方法、装置、计算机设备和存储介质 |
CN110974280A (zh) * | 2019-12-23 | 2020-04-10 | 常州博恩中鼎医疗科技有限公司 | 一种适用于锥束ct的几何参数精确自动化校正方法 |
CN111539878B (zh) * | 2020-03-25 | 2022-11-01 | 中国人民解放军战略支援部队信息工程大学 | 基于镜像投影相似度的锥束ct几何伪影去除方法及装置 |
CN111178457B (zh) * | 2020-04-09 | 2020-07-03 | 南京佗道医疗科技有限公司 | 一种精度校验器及三维影像重建精度评价方法 |
CN112233030B (zh) * | 2020-10-15 | 2023-01-03 | 上海优医基医疗器械技术有限公司 | 一种基于深度学习的ct投影几何校正方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101447085A (zh) * | 2008-12-24 | 2009-06-03 | 苏州和君科技发展有限公司 | Micro CT***中对投影图像旋转中心偏移的校正方法 |
CN103784160A (zh) * | 2014-03-06 | 2014-05-14 | 北京锐视康科技发展有限公司 | 一种锥束ct***几何位置的校正装置及其校正方法 |
CN104504656A (zh) * | 2014-12-10 | 2015-04-08 | 浙江大学 | 锥束ct图像域快速散射修正方法 |
-
2015
- 2015-07-07 CN CN201510394350.5A patent/CN105023251B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101447085A (zh) * | 2008-12-24 | 2009-06-03 | 苏州和君科技发展有限公司 | Micro CT***中对投影图像旋转中心偏移的校正方法 |
CN103784160A (zh) * | 2014-03-06 | 2014-05-14 | 北京锐视康科技发展有限公司 | 一种锥束ct***几何位置的校正装置及其校正方法 |
CN104504656A (zh) * | 2014-12-10 | 2015-04-08 | 浙江大学 | 锥束ct图像域快速散射修正方法 |
Non-Patent Citations (2)
Title |
---|
A wire scanning based method for geometric calibration of high resolution CT system;Ruijie Jiang et al.;《SPIE Medical Imaging,2015》;20150331;论文摘要、第1-3节 * |
Geometric Calibration of a Micro-CT System and Performance for Insect Imaging;Zhanli Hu et al.;《IEEE TRANSACTION ON INFORMATION TECHNOLOGY IN BIOMEDICINE》;20110731;第15卷(第4期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105023251A (zh) | 2015-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105023251B (zh) | 一种高分辨锥束ct***的几何校正方法 | |
CN108053450B (zh) | 一种基于多约束的高精度双目相机标定方法 | |
CN108122203B (zh) | 一种几何参数的校正方法、装置、设备及*** | |
WO2018218611A1 (zh) | 锥形束计算机断层成像***的几何参数确定方法 | |
CN105931202B (zh) | 几何校正模体的校正方法和*** | |
US11478214B2 (en) | Geometric calibration for cone beam CT using line fiducials | |
WO2018126335A1 (zh) | 基于小球模体的锥束ct***几何参数评价及校正方法 | |
Sun et al. | A calibration method for misaligned scanner geometry in cone-beam computed tomography | |
CN111145227B (zh) | 一种地下隧道空间多视点云的可迭代整体配准方法 | |
CN111553849B (zh) | 基于局部特征匹配的锥束ct几何伪影去除方法及装置 | |
CN101126722A (zh) | 基于配准模型仿真的锥束ct射束硬化校正方法 | |
CN101584587B (zh) | 一种ct投影中心的自动校正方法 | |
Ferrucci et al. | Measurement of the X-ray computed tomography instrument geometry by minimization of reprojection errors—Implementation on simulated data | |
CN109685877A (zh) | 一种基于自适应投影图像特征区域匹配的微纳ct焦点漂移校正方法 | |
CN108109185A (zh) | 一种生成用于消除ct伪影的校正系数的方法,以及一种基于校正系数消除ct伪影的方法 | |
CN111553960A (zh) | 一种基于投影均值图像的环状伪影快速校正方法 | |
Yang et al. | Geometry calibration method for a cone‐beam CT system | |
Nguyen et al. | A low-cost geometry calibration procedure for a modular cone-beam X-ray CT system | |
CN107016655A (zh) | 锥束cl几何全参数迭代校正方法 | |
Xu et al. | Simultaneous calibration phantom commission and geometry calibration in cone beam CT | |
CN115266779A (zh) | 一种ct扫描用平板探测器像元响应非线性校正方法 | |
Fleßner et al. | Evaluating and visualizing the quality of surface points determined from computed tomography volume data | |
CN101447085B (zh) | Micro CT ***中对投影图像旋转中心偏移的校正方法 | |
CN113876346A (zh) | 一种倾斜图像的迭代校正方法 | |
CN109919868A (zh) | 一种锥束ct射束硬化曲线侦测及投影加权校正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |