CN105869130B - Ct系统几何校正效果的检验方法及装置 - Google Patents
Ct系统几何校正效果的检验方法及装置 Download PDFInfo
- Publication number
- CN105869130B CN105869130B CN201610251571.1A CN201610251571A CN105869130B CN 105869130 B CN105869130 B CN 105869130B CN 201610251571 A CN201610251571 A CN 201610251571A CN 105869130 B CN105869130 B CN 105869130B
- Authority
- CN
- China
- Prior art keywords
- circle
- point
- mtd
- evaluation
- image
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 225
- 230000000694 effects Effects 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000007689 inspection Methods 0.000 title claims abstract description 18
- 238000011156 evaluation Methods 0.000 claims abstract description 159
- 239000000463 material Substances 0.000 claims description 17
- 238000001514 detection method Methods 0.000 claims description 12
- 230000005484 gravity Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000012795 verification Methods 0.000 claims description 5
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 229910000831 Steel Inorganic materials 0.000 description 33
- 239000010959 steel Substances 0.000 description 33
- 230000008569 process Effects 0.000 description 14
- 238000004422 calculation algorithm Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 4
- 239000004926 polymethyl methacrylate Substances 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 229920003229 poly(methyl methacrylate) Polymers 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 229920005479 Lucite® Polymers 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005693 optoelectronics Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G06T5/80—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/10—Geometric effects
- G06T15/20—Perspective computation
- G06T15/205—Image-based rendering
-
- 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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
- G06T2207/10124—Digitally reconstructed radiograph [DRR]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2215/00—Indexing scheme for image rendering
- G06T2215/06—Curved planar reformation of 3D line structures
Abstract
本发明涉及CT系统几何校正效果的检验方法及装置。所述方法包括：扫描预设的评价模体，得到所述评价模体的投影图像；根据几何校正参数以及所述评价模体的投影图像，重建所述评价模体的三维图像；从所述三维图像中分割出待处理图像，检测所述待处理图像中圆的边缘；统计出所述待处理图像中包含的圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和；计算圆内点CT均值和圆外点CT均值以及两者差的绝对值，作为评价指数；比较各组几何校正参数对应的评价指数，得出各组几何校正参数的几何校正效果的好坏。本发明能够提高评价CT图像的几何校正效果的准确度。
Description
技术领域
本发明涉及计算机断层摄影技术(Computed Tomography,CT)领域，特别是涉及CT系统几何校正效果的检验方法及装置。
背景技术
CT是用X射线束对人体某部一定厚度的层面进行扫描，由探测器接收透过该层面的X射线，转变为可见光后，由光电转换变为电信号，再经模拟/数字转换器转为数字，输入计算机处理。CT图像形成的处理有如对选定层面分成若干个体积相同的长方体，称之为体素(voxel)。扫描所得信息经计算而获得每个体素的X射线衰减系数(或吸收系数)，再排列成数字矩阵。经数字/模拟转换器把数字矩阵中的每个数字转为由黑到白不等灰度的小方块，即像素(pixel)，并按矩阵排列，即构成CT图像，所以CT图像是采用重建算法得到的重建图像。
在CT系统中，通常采用重建算法对投影数据进行重建，但是重建算法的重建效果对CT系统的几何模型要求比较高。例如，在圆轨道的CT扫描系统中，要求的几何模型包括：X线源、被测物体中心、探测器阵列中心在一条直线上。然而，在CT系统经人工安装定位后，由于机械精度达不到要求，导致实际的系统几何模型和重建算法要求的几何模型之间存在一定的几何偏差，由该几何偏差导致的重建伪影称为几何伪影。几何伪影严重影响CT图像的质量。因此CT系统需要进行几何校正，以及对几何校正的效果进行检验，为进一步提高CT图像质量提供依据。
几何校正的目的是减少重建图像的几何伪影，目前已提出了多种几何校正方法，以及对应的几何校正效果的评价方法。如一种几何校正方法：引入一种几何校正模体，该校正模体是由上下两层平行的有机玻璃板组成，上层有机玻璃板均匀镶嵌着24个钢珠，下层有机玻璃板均匀镶嵌20个钢珠。首先，把校正模体放在探测器上并保证该校正模体的投影完全在FOV(Field of View，视野)内。其次，获得每个角度下钢珠投影图像并从投影图像中分割出每个钢珠的投影坐标。最后，根据每个钢珠的空间坐标和对应的投影坐标，计算出每个投影角度下的映射矩阵和几何参数，把计算出的几何参数应用到重建算法中，由此得到几何校正后的重建图像。
但是，由于几何校正参数不能直接反映几何校正效果，因此现有的对几何校正效果的评价方法，都是根据有经验的相关人员观察CT图像得出结论。这种校验方式往往受到个人经验值以及人眼观察细致度等因素的影响，难以对几何校正效果进行准确的检验。
发明内容
基于此，本发明提供CT系统几何校正效果的检验方法及装置，能够提高评价CT图像的几何校正效果的准确度。
本发明一方面提供CT系统几何校正效果的检验方法，包括：
扫描预设的校正模体，得到所述校正模体的投影图像，根据所述校正模体及其投影图像得出CT系统的几何校正参数；
扫描预设的评价模体，得到所述评价模体的投影图像；根据几何校正参数以及所述评价模体的投影图像，重建所述评价模体的三维图像；所述评价模体由正方体和球体组成，所述球体球心和正方体重心重合，且球体半径小于正方体的内切球体半径；球体物质密度不同于正方体的物质密度；
从所述三维图像中分割出球体的图像，检测图像中的球心位置，获取球心位置所在图层作为所述评价模体的待处理图像，并检测所述待处理图像中圆的边缘；
根据圆的边缘统计出所述待处理图像中的圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和；
计算圆内点CT均值、圆外点CT均值，得出两者的差的绝对值，作为评价对应几何校正参数的几何校正效果的评价指数；
比较各组几何校正参数对应的评价指数，根据所述评价指数的大小得出各组几何校正参数的几何校正效果的好坏。
本发明另一方面提供CT系统几何校正效果的检验装置，包括：
校正参数计算模块，用于扫描预设的校正模体，得到所述校正模体的投影图像，根据所述校正模体及其投影图像得出CT系统的几何校正参数；
校正重建模块，用于扫描预设的评价模体，得到所述评价模体的投影图像；根据各几何校正参数以及所述评价模体的投影图像，重建所述评价模体的三维图像；所述评价模体由正方体和球体组成，所述球体球心和正方体重心重合，且球体半径小于正方体的内切球体半径；球体物质密度不同于正方体的物质密度；
重建图像预处理模块，用于从所述三维图像中分割出球体的图像，检测图像中的球心位置，获取球心位置所在图层作为所述评价模体的待处理图像，检测所述待处理图像中圆的边缘；
图像点识别模块，用于根据圆的边缘统计出所述待处理图像中的圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和；
评价指数计算模块，用于计算圆内点CT均值、圆外点CT均值，得出两者的差的绝对值，作为评价对应几何校正参数的几何校正效果的评价指数；
比较及校验模块，用于比较各组几何校正参数对应的评价指数，根据所述评价指数的大小得出各组几何校正参数的几何校正效果的好坏。
上述技术方案，通过扫描预设的评价模体，得到所述评价模体的投影图像；根据几何校正参数以及所述评价模体的投影图像，重建所述评价模体的三维图像；所述评价模体由正方体和位于所述正方体内部的球体组成，所述球体球心和正方体重心重合；球体物质密度不同于正方体的物质密度；从所述三维图像中分割出球体的图像，检测图像中的球心位置，获取球心位置所在图层作为所述评价模体的待处理图像；检测所述待处理图像中圆的边缘，根据圆的边缘统计出所述待处理图像中的圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和；计算圆内点CT均值、圆外点CT均值，得出两者的差的绝对值，作为评价对应几何校正参数的几何校正效果的评价指数；比较各组几何校正参数对应的评价指数，根据所述评价指数的大小得出各组几何校正参数的几何校正效果的好坏。
本发明上述技术方案，根据几何校正参数重建评价模体的CT图像，并计算出评价指数对几何校正效果进行校验，所述评价指数能够客观反映CT图像的几何校正效果，为评价几何校正效果提供统一的标准，避免了人为评估几何校正效果时的主观因素，有利于提高评价CT图像的几何校正效果的准确度。
附图说明
图1为一实施例的CT系统几何校正效果的检验方法的示意性流程图；
图2为一实施例中评价模体的示意结构图；
图3为一实施例中待处理图像的边缘检测结果图；
图4为一实施例中待处理图像逐点判断密度区域示意图；
图5为一实施例的评价指数与参数R的关系图；
图6为一实施例的评价指数与r/R比值的关系图；
图7为一实施例的评价指数与几何校正效果示意图；
图8为一实施例的评价指数准确度的展示图；
图9为一实施例的CT系统几何校正效果的检验装置的示意性结构图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
图1为一实施例的CT系统几何校正效果的检验方法的示意性流程图；如图1所示，本实施例中的CT系统几何校正效果的检验方法包括步骤：
S11，扫描预设的校正模体，得到所述校正模体的投影图像，根据所述校正模体及其投影图像得出CT系统的几何校正参数；并模拟不同几何校正效果下的几何校正参数；
优选的，本实施例中的校正模体为：表面等螺距等角度的镶嵌多个标记点的空心圆柱体。例如：该校正模体是在空心圆柱表面等螺距等角度的镶嵌钢珠，总共18个钢珠，钢珠直径是1mm。对应的，步骤S11中根据所述校正模体及其投影图像得出CT系统的几何校正参数的方式包括：获得每个角度下钢珠的投影图像，并从钢珠的投影图像中检测出钢珠的投影坐标，再根据每个钢珠在CT系统中的空间坐标及其对应的投影坐标，得出每个投影角度下所述CT系统的几何校正参数。具体可为：获得每个角度下钢珠的投影图像，并从钢珠的投影图像中检测出钢珠的投影坐标；由钢珠的空间坐标值与投影坐标值的关系，计算出每个投影角度下所述CT系统的几何校正参数。本实施例中，通过在钢珠的空间坐标值加上一个随机数，以模拟CT系统不同几何校正效果下的几何校正参数。
S12，扫描预设的评价模体，得到所述评价模体的投影图像；根据几何校正参数以及所述评价模体的投影图像，重建所述评价模体的三维图像；所述评价模体由正方体和球体组成，所述球体球心和正方体重心重合，且球体半径小于正方体的内切球体半径；球体物质密度不同于正方体的物质密度；
优选的，如图2所示，本实施例中所述评价模体由正方体和位于所述正方体内部的球体组成，所述球体球心和正方体重心重合，且球体半径小于正方体的内切球体半径；球体物质密度不同于正方体的物质密度。例如：球体内物质为空气，球体与正方体之间填充高密度树脂材料。
本步骤中将几何校正参数加入CT系统的重建算法中，结合所述评价模体的投影图像重建所述评价模体的三维图像。即步骤S11中计算出的各组几何参数分别应用到CT系统的重建算法中，得到与各个几何参数对应的几何校正后的重建图像(所述评价模体的三维图像)。
S13，从所述三维图像中分割出球体的图像，检测图像中的球心位置，获取球心位置所在图层作为所述评价模体的待处理图像，并检测所述待处理图像中圆的边缘；
优选的，检测所述待处理图像中圆的边缘的过程包括：
首先，计算所述待处理图像的梯度，得到待处理图像中圆的模糊边缘(如图3所示)，计算公式包括：
其中，代表图像f的梯度，符号*表示卷积，g表示Sobel模板，其定义为：
然后，获取设定的图像梯度阈值，根据所述图像梯度阈值去除所述模糊边缘中的噪声点，得到待处理图像中圆的边缘。
S14，根据圆的边缘统计出所述待处理图像中的圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和；
首先，根据圆的边缘确定出所述待处理图像中圆的圆心位置和半径；
然后，识别所述待处理图像中的圆内点和圆外点，识别方法包括：
将满足条件（xi-a）2＋（yi-b）2≤r2的点识别为圆内点，将满足条件r2≤(xi-a)2＋(yi-b)2≤(r＋R)2的点识别为圆外点；其中，(xi,yi)为所述待处理图像中第i个点的坐标，(a,b)为所述待处理图像中圆的圆心坐标，r为所述待处理图像中圆的半径，R为预设的参数，其取值范围为大于0且小于正方体内切球半径与球体半径的差值；
最后，根据识别结果统计出圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和。
S15，计算圆内点CT均值、圆外点CT均值，得出两者的差的绝对值，作为评价对应几何校正参数的几何校正效果的评价指数；
例如：如果圆内点总数为A，圆内点的CT值总和为M，圆外点(即图4中圆环内的点)总数为B，圆外点的CT值总和为N，可得到对应的评价指数F为：
S16，比较各组何校正参数对应的评价指数，根据所述评价指数的大小得出各组几何校正参数的几何校正效果的好坏。
重复上述步骤S12-S16可得不同几何校正参数的几何校正效果对应的评价指数，评价指数对应的重建图像如图7所示，可知所述评价指数与几何校正效果之间存在正相关的关系，即评价指数越高，对应的CT系统的几何校正效果越好；反之，几何校正效果越差。
优选的，步骤S11中扫描校正模体的方式可为：射线源和探测器绕所述校正模体旋转一周，等间隔采集所述校正模体的投影图像；所述校正模体位于旋转中心并且所述校正模体的投影完全在所述CT系统的FOV内。具体如：把校正模体放在旋转中心并保证校正模体的投影完全在FOV内；设置CT系统扫描参数，管电压设为85kV，管电流设为6mA，扫描间隔角度为1°，模式为快速扫描。获得每个角度下校正模体的投影图像。
优选的，步骤S12中扫描评价模体的方式包括：射线源和探测器绕所述评价模体旋转一周，等间隔采集所述评价模体的投影图像；所述评价模体位于旋转中心并且所述评价模体的投影完全在所述CT系统的FOV内。例如：把评价模体放入CT系统，采用快速扫描模式，射线源和探测器绕评价模体旋转一周，每旋转1°采集一张评价模体的投影图像，共采集360张投影图像。
为了便于理解本发明实施例的CT系统几何校正效果的检验方法，下面通过一个具体应用示例进行进一步的说明。使用的CT扫描器为牙科CT系统。所述CT系统几何校正效果的检验方法过程如下：
过程(1)，首先设置几何校正模体，该校正模体是在空心圆柱表面等螺距等角度的镶嵌钢珠，总共18个钢珠。校正模体的尺寸是10*10*20cm3，钢珠直径是1mm。其次，把校正模体放在旋转中心并保证校正模体的投影完全在FOV内。再次，设置CT系统扫描参数，管电压设为85kV，管电流设为6mA，扫描间隔角度为1°，模式为快速扫描。获得每个角度下钢珠投影图像并从投影图像中检测出钢珠的投影坐标。最后，由钢珠的空间坐标值与投影坐标值的关系，计算出CT系统的几何校正参数。优选的，通过在钢珠的空间坐标值加上一个随机数，以模拟CT系统不同几何校正效果下的几何校正参数。
过程(2)，设计并扫描评价模体。评价模体由边长为96mm的正方体有机玻璃外壳和半径为18mm的球体组成，球体球心和正方体的重心重合，球体与正方体之间填充高密度树脂。把评价模体放入CT系统，采用快速扫描模式，射线源和探测器绕评价模体旋转一周，每1°采集一张评价模体的投影图像，共采集360张投影图像。
过程(3)，利用几何校正参数和采集的评价模体的投影，将几何校正参数加入所述CT系统的重建算法，得到重建出评价模体的三维图像。
过程(4)，首先，选择三维图像中心层附近的图像，设定阈值，通过阈值分割出球体的图像，计算球体的图像中球心位置，获取球心位置所在图层作为所述评价模体的待处理图像。
其次，求得待处理图像的梯度图，设置阈值，对梯度图进行分割得到梯度图的逻辑图像，即待处理图像中的圆的模糊边缘(如图3所示)。具体计算公式如下：
其中，其中代表图像f的梯度；符号*表示卷积，g表示Sobel模板，其定义如下：
通过以上公式可以得到待处理图像的梯度。设置图像梯度阈值T＝50，去除模糊边缘中对圆的边缘贡献小的点，这样，得到待处理图像中圆的边缘。
接着，通过待处理图像中圆的边缘求得待处理图像的圆心和半径。
最后，通过比较待处理图像上任一点到圆心的距离与圆的半径的大小，判断该点是在低密度区(圆内)还是高密度区(圆外)，统计出所述待处理图像中包含的圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和。过程包括：
（xi-a）2＋(yi-b)2≤r2 (1)
r2≤(xi-a)2+(yi-b)2≤(r＋R)2,0≤R≤30mm (2)
若待处理图像中的点满足条件(1)，则该点为圆内点；若满足条件(2)，则该点为圆外点(即图4中圆环内的点)。
例如，若圆内点总数为A，圆内点的CT值总和为M，圆外点总数为B，圆外点的CT值总和为N，可得到对应的评价指数F为：
过程(5)，调整球体的半径r和参数R的比值，确定最优比的范围。首先，设定低密度区圆的半径r＝18mm，对所述参数R以0.3mm为间隔在0-30mm范围内自动取值，按照上述过程(4)求得不同几何校正效果对应的评价指数曲线F。评价指数F与参数R的关系如图5所示，图中横轴代表R的取值，纵轴代表评价指数F的取值，不同曲线代表不同的几何校正效果。其结果如表1所示。
表1：
由表1可知，所述参数R的取值为(1mm～5mm)时，评价指数能较好的反应不同几何校正效果。
将r/R的比值作为横轴，评价指数F为纵轴，不同曲线代表不同的几何校正效果，即可得到r/R与评价指数F的关系，如图6。其结果如表2所示。
表2：
由表2可知，球体半径r与所述参数R的比值r/R为4～6时，评价指数能较好的反应不同几何校正效果。
过程(6)，用评价指数评价几何校正效果。为了模拟不同几何校正效果对应的评价指数，取R＝3mm，r＝18mm。首先，在过程(1)中计算几何校正参数时，在钢珠的空间坐标值加一个随机数，就可以得到不同几何校正效果对应的几何参数。经过过程(2)-(4)可得不同几何校正效果对应的评价指数，各评价指数对应的重建图像如图7所示。评价指数与几何校正效果之间有正相关的关系，即评价指数越高，CT系统的几何校正效果越好；反之，评价指数越低，CT系统的几何校正效果越差。
过程(7)，与评价图像空间分辨率(即观察线对数目)的方法作比较，展示评价指数能够提高评价CT图像的几何校正效果的准确度。为了模拟细微差别的几何校正效果对应的评价指数，取R＝3mm，r＝18mm。首先，在过程(1)中计算几何校正参数时，在钢珠的空间坐标值加一个限制在0.05mm内的随机数，就可以得到细微差别的几何校正效果对应的几何校正参数。经过过程(2)-(4)可得不同几何校正效果对应的评价指数，同时对线对分辨率模体进行扫描重建，各评价指数对应的线对模体重建图像如图8所示。评价指数可以把细微差别的几何校正效果分为四个等级，但线对模体显示的均为5lp/cm，说明评价指数能够提高评价CT图像的几何校正效果的准确度。
可见，本发明实施例提供了CT系统几何校正效果的检验方法，提出评价指数，该指数能够客观的反应出几何校正效果，为评价几何校正效果提供统一的标准，因此避免了由于主观因素导致评价几何校正效果出现误差的问题。并且所需的几何校正模体的结构简单，制造成本低廉，易于实现。
需要说明的是，对于前述的各方法实施例，为了简便描述，将其都表述为一系列的动作组合，但是本领域技术人员应该知悉，本发明并不受所描述的动作顺序的限制，因为依据本发明，某些步骤可以采用其它顺序或者同时进行。
基于与上述实施例中的CT系统几何校正效果的检验方法相同的思想，本发明还提供CT系统几何校正效果的检验装置，该装置可用于执行上述CT系统几何校正效果的检验方法。为了便于说明，CT系统几何校正效果的检验装置实施例的结构示意图中，仅仅示出了与本发明实施例相关的部分，本领域技术人员可以理解，图示结构并不构成对装置的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
图9为本发明一实施例的CT系统几何校正效果的检验装置的示意性结构图；如图9所示，本实施例的CT系统几何校正效果的检验装置包括：校正参数计算模块710、校正重建模块720、重建图像预处理模块730、重建图像预处理模块740、评价指数计算模块750以及比较及校验模块760，各模块详述如下：
上述校正参数计算模块710，用于扫描预设的校正模体，得到所述校正模体的投影图像，根据所述校正模体及其投影图像得出CT系统的几何校正参数；
上述校正重建模块720，用于扫描预设的评价模体，得到所述评价模体的投影图像；根据各几何校正参数以及所述评价模体的投影图像，重建所述评价模体的三维图像；所述评价模体由正方体和位于所述正方体内部的球体组成，所述球体球心和正方体重心重合，且球体半径小于正方体的内切球体半径；球体的物质密度不同于正方体的物质密度；
上述重建图像预处理模块730，用于从所述三维图像中分割出球体的图像，检测图像中的球心位置，获取球心位置所在图层作为所述评价模体的待处理图像，检测所述待处理图像中圆的边缘；
上述图像点识别模块740，用于根据圆的边缘统计出所述待处理图像中的圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和；
上述评价指数计算模块750，用于计算圆内点CT均值、圆外点CT均值，得出两者的差的绝对值，作为评价对应几何校正参数的几何校正效果的评价指数；
上述比较及校验模块760，用于比较各组几何校正参数对应的评价指数，根据所述评价指数的大小得出各组几何校正参数的几何校正效果的好坏。优选的，所述图像点识别模块740包括:
圆心及半径子模块，用于根据圆的边缘确定出所述待处理图像中圆的圆心位置和半径；识别子模块，用于识别所述待处理图像中的圆内点和圆外点，包括：将满足条件(xi-a)2+(yi-b)2≤r2的点识别为圆内点，将满足条件r2≤(xi-a)2+(yi-b)2≤(r+R)2的点识别为圆外点；以及，统计子模块，用于根据识别结果统计出圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和。其中，(xi,yi)为所述待处理图像中第i个点的坐标，(a,b)为所述待处理图像中圆的圆心坐标，r为所述待处理图像中圆的半径，R为预设的参数，其取值范围为大于0且小于正方体内切球半径与球体半径的差。
优选的，所述重建图像预处理模块730包括：边缘检测子模块，用于计算所述待处理图像的梯度，得到待处理图像中圆的模糊边缘，计算公式包括：
其中，代表图像f的梯度，符号*表示卷积，g表示Sobel模板，其定义为：
获取设定的图像梯度阈值，根据所述图像梯度阈值去除所述模糊边缘中的噪声点，得到待处理图像中圆的边缘。
优选的，所述校正参数计算模块710包括：校正模体扫描子模块，用于检测到所述校正模体位于旋转中心并且所述校正模体的投影完全在所述CT系统的FOV内，扫描预设的校正模体，使射线源和探测器绕所述校正模体旋转一周，等间隔采集所述校正模体的投影图像.
所述校正重建模块720包括：评价模体子模块，用于检测到所述评价模体位于旋转中心并且所述评价模体的投影完全在所述CT系统的FOV内，扫描预设的评价模体，使射线源和探测器绕所述评价模体旋转一周，等间隔采集所述评价模体的投影图像。
进一步的，所述校正模体为：表面等螺距等角度的镶嵌多个标记点(例如钢珠)的空心圆柱体。对应的，所述校正参数计算模块710还包括：参数计算子模块，用于获得每个角度下钢珠的投影图像，并从钢珠的投影图像中检测出钢珠的投影坐标；获取每个钢珠在CT系统中的空间坐标值，用一个随机数对每个钢珠的空间坐标值进行校正；根据每个钢珠对应的投影坐标及其校正后的空间坐标，得出每个投影角度下所述CT系统的几何校正参数。
需要说明的是，上述示例的CT系统几何校正效果的检验装置的实施方式中，各模块/单元之间的信息交互、执行过程等内容，由于与本发明前述方法实施例基于同一构思，其带来的技术效果与本发明前述方法实施例相同，具体内容可参见本发明方法实施例中的叙述，此处不再赘述。
此外，上述示例的CT系统几何校正效果的检验装置的实施方式中，各功能模块的逻辑划分仅是举例说明，实际应用中可以根据需要，例如出于相应硬件的配置要求或者软件的实现的便利考虑，将上述功能分配由不同的功能模块完成，即将所述CT系统几何校正效果的检验装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。其中各功能模既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。
一种终端，包括内存和处理器，所述内存中储存有指令，所述指令被所述处理器执行时，可使得所述处理器执行以下步骤：
扫描预设的校正模体，得到所述校正模体的投影图像，根据所述校正模体及其投影图像得出CT系统的几何校正参数；
扫描预设的评价模体，得到所述评价模体的投影图像；根据各几何校正参数以及所述评价模体的投影图像，重建所述评价模体的三维图像；所述评价模体由正方体球体组成，所述球体球心和正方体重心重合，且球体半径小于正方体的内切球体半径；球体内的填充物质密度低于球体与正方体之间的填充物质密度；
从所述三维图像中分割出球体的图像，检测图像中的球心位置，获取球心位置所在图层作为所述评价模体的待处理图像，并检测所述待处理图像中圆的边缘；
根据圆的边缘统计出所述待处理图像中的圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和；
计算圆内点CT均值、圆外点CT均值，得出两者的差的绝对值，作为评价对应几何校正参数的几何校正效果的评价指数；
比较各几何校正参数对应的评价指数，根据所述评价指数的大小得出各几何校正参数的几何校正效果的高低。
本领域普通技术人员可以理解，实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，作为独立的产品销售或使用。所述程序在执行时，可执行如上述各方法的实施例的全部或部分步骤。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-OnlyMemory，ROM)或随机存储记忆体(Random Access Memory，RAM)等。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其它实施例的相关描述。
以上所述实施例仅表达了本发明的几种实施方式，不能理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。
Claims (8)
1.一种CT系统几何校正效果的检验方法，其特征在于，包括：
扫描预设的校正模体，得到所述校正模体的投影图像，根据所述校正模体及其投影图像得出CT系统的几何校正参数；
扫描预设的评价模体，得到所述评价模体的投影图像；根据几何校正参数以及所述评价模体的投影图像，重建所述评价模体的三维图像；所述评价模体由正方体和球体组成，所述球体球心和正方体重心重合，且球体半径小于正方体的内切球体半径；球体的物质密度不同于正方体的物质密度；
从所述三维图像中分割出球体的图像，检测图像中的球心位置，获取球心位置所在图层作为所述评价模体的待处理图像，并检测所述待处理图像中圆的边缘；
根据圆的边缘统计出所述待处理图像中的圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和；
计算圆内点CT均值、圆外点CT均值，得出两者的差的绝对值，作为评价对应几何校正参数的几何校正效果的评价指数；
比较各组几何校正参数对应的评价指数，根据所述评价指数的大小得出各组几何校正参数的几何校正效果的好坏；
所述扫描预设的校正模体，得到所述校正模体的投影图像，包括：
检测到所述校正模体位于旋转中心并且所述校正模体的投影完全在所述CT系统的FOV内，扫描预设的校正模体，使射线源和探测器绕所述校正模体旋转一周，等间隔采集所述校正模体的投影图像；
所述扫描预设的评价模体，得到所述评价模体的投影图像，包括：
检测到所述评价模体位于旋转中心并且所述评价模体的投影完全在所述CT系统的FOV内，扫描预设的评价模体，使射线源和探测器绕所述评价模体旋转一周，等间隔采集所述评价模体的投影图像；
所述校正模体为：表面等螺距等角度的镶嵌多个标记点的空心圆柱体。
2.根据权利要求1所述的CT系统几何校正效果的检验方法，其特征在于，所述根据圆的边缘统计出所述待处理图像中的圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和，包括：
根据圆的边缘确定出所述待处理图像中圆的圆心位置和半径；
识别所述待处理图像中的圆内点和圆外点，包括：
将满足条件(xi-a)2+(yi-b)2≤r2的点识别为圆内点，
将满足条件r2≤(xi-a)2+(yi-b)2≤(r+R)2的点识别为圆外点；
根据识别结果统计出圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和；
其中，(xi,yi)为所述待处理图像中第i个点的坐标，(a,b)为所述待处理图像中圆的圆心坐标，r为所述待处理图像中圆的半径，R为预设的参数，其取值范围为大于0且小于正方体内切球半径与球体半径的差。
3.根据权利要求1所述的CT系统几何校正效果的检验方法，其特征在于，所述检测所述待处理图像中圆的边缘，包括：
计算所述待处理图像的梯度，得到待处理图像中圆的模糊边缘，计算公式包括：
|▽f|2＝(g*f)2+(gT*f)2
其中，|▽f|代表图像f的梯度，符号*表示卷积，g表示Sobel模板，其定义为：
<mrow>
<mi>g</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
根据设定的图像梯度阈值去除所述模糊边缘中的噪声点，得到待处理图像中圆的边缘。
4.根据权利要求1所述的CT系统几何校正效果的检验方法，其特征在于，所述根据所述校正模体及其投影图像得出CT系统的几何校正参数，包括：
获得每个角度下标记点的投影图像，并从标记点的投影图像中检测出标记点的投影坐标；由标记点的空间坐标值与投影坐标值的关系，计算出每个投影角度下所述CT系统的几何校正参数。
5.一种CT系统几何校正效果的检验装置，其特征在于，包括：
校正参数计算模块，用于扫描预设的校正模体，得到所述校正模体的投影图像，根据所述校正模体及其投影图像得出CT系统的几何校正参数；
校正重建模块，用于扫描预设的评价模体，得到所述评价模体的投影图像；根据几何校正参数以及所述评价模体的投影图像，重建所述评价模体的三维图像；所述评价模体由正方体和球体组成，所述球体球心和正方体重心重合，且球体半径小于正方体的内切球体半径；球体的物质密度不同于正方体的物质密度；
重建图像预处理模块，用于从所述三维图像中分割出球体的图像，检测图像中的球心位置，获取球心位置所在图层作为所述评价模体的待处理图像，检测所述待处理图像中圆的边缘；
图像点识别模块，用于根据圆的边缘统计出所述待处理图像中的圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和；
评价指数计算模块，用于计算圆内点CT均值、圆外点CT均值，得出两者的差的绝对值，作为评价对应几何校正参数的几何校正效果的评价指数；
比较及校验模块，用于比较各组几何校正参数对应的评价指数，根据所述评价指数的大小得出各组几何校正参数的几何校正效果的好坏；
所述校正参数计算模块包括：
校正模体扫描子模块，用于检测到所述校正模体位于旋转中心并且所述校正模体的投影完全在所述CT系统的FOV内，扫描预设的校正模体，使射线源和探测器绕所述校正模体旋转一周，等间隔采集所述校正模体的投影图像；
所述校正重建模块包括：
评价模体子模块，用于检测到所述评价模体位于旋转中心并且所述评价模体的投影完全在所述CT系统的FOV内，扫描预设的评价模体，使射线源和探测器绕所述评价模体旋转一周，等间隔采集所述评价模体的投影图像；
所述校正模体为：表面等螺距等角度的镶嵌多个标记点的空心圆柱体。
6.根据权利要求5所述的CT系统几何校正效果的检验装置，其特征在于，所述图像点识别模块包括:
圆心及半径子模块，用于根据圆的边缘确定出所述待处理图像中圆的圆心位置和半径；
识别子模块，用于识别所述待处理图像中的圆内点和圆外点，包括：
将满足条件(xi-a)2+(yi-b)2≤r2的点识别为圆内点，
将满足条件r2≤(xi-a)2+(yi-b)2≤(r+R)2的点识别为圆外点；
统计子模块，用于根据识别结果统计出圆内点数、圆内点CT值总和、圆外点数以及圆外点CT值总和；
其中，(xi,yi)为所述待处理图像中第i个点的坐标，(a,b)为所述待处理图像中圆的圆心坐标，r为所述待处理图像中圆的半径，R为预设的参数，其取值范围为大于0且小于正方体内切球半径与球体半径的差。
7.根据权利要求5所述的CT系统几何校正效果的检验装置，其特征在于，所述重建图像预处理模块包括：
边缘检测子模块，用于计算所述待处理图像的梯度，得到待处理图像中圆的模糊边缘，计算公式包括：
|▽f|2＝(g*f)2+(gT*f)2
其中，|▽f|代表图像f的梯度，符号*表示卷积，g表示Sobel模板，其定义为：
<mrow>
<mi>g</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
获取设定的图像梯度阈值，根据所述图像梯度阈值去除所述模糊边缘中的噪声点，得到待处理图像中圆的边缘。
8.根据权利要求5所述的CT系统几何校正效果的检验装置，其特征在于，所述校正参数计算模块还包括：
参数计算子模块，用于获得每个角度下标记点的投影图像，并从标记点的投影图像中检测出标记点的投影坐标；由标记点的空间坐标值与投影坐标值的关系，计算出CT系统的几何参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610251571.1A CN105869130B (zh) | 2016-04-20 | 2016-04-20 | Ct系统几何校正效果的检验方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610251571.1A CN105869130B (zh) | 2016-04-20 | 2016-04-20 | Ct系统几何校正效果的检验方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105869130A CN105869130A (zh) | 2016-08-17 |
| CN105869130B true CN105869130B (zh) | 2018-02-23 |
Family
ID=56632743
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610251571.1A Active CN105869130B (zh) | 2016-04-20 | 2016-04-20 | Ct系统几何校正效果的检验方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105869130B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108122203B (zh) * | 2016-11-29 | 2020-04-07 | 上海东软医疗科技有限公司 | 一种几何参数的校正方法、装置、设备及系统 |
| CN112529820A (zh) * | 2020-12-29 | 2021-03-19 | 广州华端科技有限公司 | Cbct系统半扇束几何校正方法、装置、设备和介质 |
| CN113643428A (zh) * | 2021-08-17 | 2021-11-12 | 北京唯迈医疗设备有限公司 | 一种适用于多自由度锥形束ct的全参数几何校准方法 |
| CN116359257A (zh) * | 2021-12-27 | 2023-06-30 | 同方威视技术股份有限公司 | 标定组件、标定模体及标定方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103716602A (zh) * | 2013-12-24 | 2014-04-09 | 北京淳中视讯科技有限公司 | 投影图像的几何校正方法、装置及系统 |
| CN103961125A (zh) * | 2013-01-31 | 2014-08-06 | 东北大学 | 一种用于锥束ct的ct值校正方法 |
| CN105023251A (zh) * | 2015-07-07 | 2015-11-04 | 苏州海斯菲德信息科技有限公司 | 一种高分辨锥束ct系统的几何校正方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006030501A1 (ja) * | 2004-09-15 | 2006-03-23 | Mitsubishi Denki Kabushiki Kaisha | 画像投影システム及び画像幾何補正装置 |
-
2016
- 2016-04-20 CN CN201610251571.1A patent/CN105869130B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103961125A (zh) * | 2013-01-31 | 2014-08-06 | 东北大学 | 一种用于锥束ct的ct值校正方法 |
| CN103716602A (zh) * | 2013-12-24 | 2014-04-09 | 北京淳中视讯科技有限公司 | 投影图像的几何校正方法、装置及系统 |
| CN105023251A (zh) * | 2015-07-07 | 2015-11-04 | 苏州海斯菲德信息科技有限公司 | 一种高分辨锥束ct系统的几何校正方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105869130A (zh) | 2016-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105931202B (zh) | 几何校正模体的校正方法和系统 | |
| CN104107065B (zh) | 3d图像集在不同空间之间的最佳变换 | |
| CN105869130B (zh) | Ct系统几何校正效果的检验方法及装置 | |
| CN100587391C (zh) | 一种适用于2d-ct扫描系统的投影旋转中心测量方法 | |
| Stevens et al. | Alignment of a volumetric tomography system | |
| US8000435B2 (en) | Method and system for error compensation | |
| US20090202127A1 (en) | Method And System For Error Compensation | |
| CN103961125B (zh) | 一种用于锥束ct的ct值校正方法 | |
| US7672424B2 (en) | Image reconstruction with voxel dependent interpolation | |
| US8075183B2 (en) | Method of evaluating the resolution of a volumetric imaging system and image phantom used during the resolution evaluation | |
| CN106667512A (zh) | X射线成像设备的几何校正方法、乳腺断层成像设备 | |
| CN108511043B (zh) | 基于数值模拟的x-ct虚拟数据采集及图像重建方法及系统 | |
| US11504085B2 (en) | Method for calibrating defective channels of a CT device | |
| CN206315096U (zh) | Ct性能检测模体 | |
| CN107948637A (zh) | 一种模体图像质量检测方法和系统 | |
| WO2023124452A1 (zh) | 标定组件、标定模体及标定方法 | |
| CN106572832A (zh) | 数据处理方法、数据处理装置以及x射线ct装置 | |
| CN109146987B (zh) | 一种基于gpu的快速锥束计算机断层成像重建方法 | |
| CN102488528B (zh) | 一种层析成像几何参数的校准方法 | |
| US7769217B2 (en) | Fast iterative 3D PET image reconstruction using a set of 2D linogram transformations | |
| CN108109185A (zh) | 一种生成用于消除ct伪影的校正系数的方法，以及一种基于校正系数消除ct伪影的方法 | |
| CN107157505A (zh) | 一种提高锥形束spect成像分辨率的方法 | |
| CN107080552B (zh) | 一种测试多种医学x射线设备成像质量的体模及其系统、方法 | |
| CN109919868A (zh) | 一种锥束ct射束硬化曲线侦测及投影加权校正方法 | |
| US7058156B2 (en) | Imaging method for a multi-slice spiral CT scan with 3D reconstruction, and a computed tomography unit for carrying out this method |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220913 Address after: Room 901, China Railway South Headquarters Building, No. 3333, Central Road, Azure Coast Community, Yuehai Street, Nanshan District, Shenzhen, Guangdong 518000 Patentee after: Shenzhen Aomei Future Medical Technology Co.,Ltd. Address before: 10th Floor, Qixi Building, No. 286, Science Avenue, Luogang District, Guangzhou City, Guangdong Province, 510530 Patentee before: GUANGZHOU HUADUAN TECHNOLOGY Co.,Ltd. |