CN103033525A - Ct系统和ct图像重建方法 - Google Patents
Ct系统和ct图像重建方法 Download PDFInfo
- Publication number
- CN103033525A CN103033525A CN2011102927738A CN201110292773A CN103033525A CN 103033525 A CN103033525 A CN 103033525A CN 2011102927738 A CN2011102927738 A CN 2011102927738A CN 201110292773 A CN201110292773 A CN 201110292773A CN 103033525 A CN103033525 A CN 103033525A
- Authority
- CN
- China
- Prior art keywords
- camera
- mark group
- examined object
- camera set
- processing unit
- 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
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种CT系统和CT图像重建方法。所述CT系统包括射线源、射线探测器、数据采集与处理单元、主控制器、相机组和标记组以及标记组处理单元,其中,所述标记组被固定在受检物体表面;所述相机组用来获得所述标记组在其上的投影;所述标记组处理单元根据所述投影得出所述标记组的初始位置以及旋转后的位置信息,从而计算得到受检物体的旋转矩阵和位移量。本发明CT图像重建方法包括步骤:将标记组固定在受检物体表面并得到所述标记组的初始位置;在每次受检物体被旋转后获得旋转后标记组的位置信息;接着获得受检物体的旋转矩阵和位移量;根据所述旋转矩阵和位移量来重建CT图像。采用本发明CT系统和CT图像重建方法对机械控制精度要求不高。
Description
技术领域
本发明总体涉及X光计算机断层成像领域,且特别涉及一种CT系统和CT图像重建方法。
背景技术
计算机断层成像(CT)是一种基于X射线扫描数据而得到受检物体图像,进而可以得知受检物体的内部信息的检测技术。该技术已经广泛应用于医疗、机械制造、航空航天、国防军工等领域。使用该技术得到的受检物体的图像具有清晰直观、分辨率高、便于分析和存储等优点。
传统CT的实现是基于X光机和高能加速器等X射线源对受检物体进行照射,并在对应的探测器采集数据,在扫描过程中通过CT系统的旋转架等的旋转平移等而使X射线在不同角度下穿过受检物体,从而得到受检物体不同位置的光强数据(也即扫描数据,又可称为投影数据),利用所得到的扫描数据可以重建出受检物体的断层图像。
传统CT系统需要精确地控制扫描过程中受检物体和加速器以及探测器的相对位置,以满足对断层图像重建所需的数据条件。因此传统CT系统通常包括射线源、射线探测器、数据采集与处理系统、主控制器以及机械控制系统,其中,所述射线探测器探测所述射线源发射的X射线并将其转换成数字信号以传送给数据采集与处理单元;所述射线源和所述数据采集与处理单元以及所述机械控制系统由所述主控制器控制。另外,机械控制系统用来控制受检物体和加速器(即射线源)以及射线探测器之间的相对位置。传统CT系统的机械控制系统通常包括直线导轨、旋转架等。由于该机械控制系统结构复杂,所以不仅使传统CT系统的整体成本高昂,而且为了保证系统的机械精度,安装一套这样的工业CT系统一般需要长时间的调试,对异形大尺寸受检物体还存在很大的机械布局和控制的实现困难。在工业应用中,这使得我们无法实现CT系统快速部署和移动,同时无法适应复杂环境。
另外,随着断层图像重建方法技术的发展,使得对于CT数据构成的要求得以降低。例如基于代数迭代的重建方法结合先验条件可以在几何关系不规律的有限扫描数据情况下实现对受检物体的断层图像的重建。
发明内容
本发明要解决的技术问题是提供一种对机械控制精度要求不高的CT系统和CT图像重建方法。
为了解决上述技术问题,本发明CT系统的技术方案包括射线源、射线探测器、数据采集与处理单元、主控制器,其中,所述射线探测器探测所述射线源发射的X射线并将其转换成数字信号以传送给所述数据采集与处理单元;所述射线源和所述数据采集与处理单元由所述主控制器控制,所述CT系统还包括相机组和标记组以及标记组处理单元,其中:
所述标记组被固定在受检物体表面;
所述相机组用来获得所述标记组在其上的投影;
所述标记组处理单元根据所述投影得出所述标记组的初始位置以及旋转后的位置信息,从而计算得到受检物体的旋转矩阵和位移量。
优选地,所述相机组包括以柱面透镜作为镜头的线阵CCD相机。
优选地,所述相机组包括以球面透镜作为镜头的面阵CCD相机。
其中,所述标记组中的标记点由单色反光材料制成或者具有单色反光材料的涂层。
所述相机组中的相机的镜头上添加有对应颜色的滤光片。
优选地,所述标记组至少包括三个标记点。
所述CT系统可以包括多个相机组且这多个相机组的视野范围能够覆盖在扫描过程中所述标记组中的所有标记点可能出现的区域。
优选地,所述相机组包括至少三个相机且这些相机分布在以被监视区域中心为中心的圆弧上。
进一步地,以所述中心为坐标原点建立坐标系,如果在xy平面内以被监视区域中心为中心的圆弧上分布了三个相机,则左右两个相机的镜头顺着z轴方向摆放而中间的相机的镜头平行于xy平面方向摆放。
相应地,本发明CT图像重建方法包括步骤:
将标记组固定在受检物体表面并通过相机组和标记组处理单元得到所述标记组的初始位置;
在每次受检物体被旋转到某位置后且在该某位置对受检物体进行扫描前通过相机组和标记组处理单元获得旋转后标记组的位置信息;
根据所述标记组的初始位置和旋转后的位置信息获得受检物体的旋转矩阵和位移量;
根据所述旋转矩阵和位移量以及重建参数来重建CT图像。
优选地,所述相机组包括以柱面透镜作为镜头的线阵CCD相机。
优选地,所述相机组包括以球面透镜作为镜头的面阵CCD相机。
其中,所述标记组中的标记点由单色反光材料制成或者具有单色反光材料的涂层。
所述相机组中的相机的镜头上添加有对应颜色的滤光片。
优选地,所述标记组至少包括三个标记点。
进一步地,所述相机组的数目为两个以上且这些相机组的视野范围能够覆盖在扫描过程中所述标记组中的所有标记点可能出现的区域。
另外,所述相机组包括至少三个相机且这些相机分布在以被监视区域中心为中心的圆弧上。
此外,以所述中心为坐标原点建立坐标系,如果在xy平面内以被监视区域中心为中心的圆弧上分布了三个相机,则左右两个相机的镜头顺着z轴方向摆放而中间的相机的镜头平行于xy平面方向摆放。
与现有技术相比,本发明CT系统和CT图像重建方法的有益效果为:
首先,由于本发明CT系统没有采用传统CT系统中的复杂的机械控制系统,而取代之的是相机组和标记组以及标记组处理单元。因此,本发明CT系统在机械控制方面的精度要求不高且容易安装和部署,从而具有更大的灵活性和适用性。
其次, 由于本发明采用将标记组固定在受检物体表面,也就是说标记组能够随着受检物体的旋转而旋转。这样根据标记组初始位置和旋转后的位置信息可以获得受检物体的旋转矩阵和位移量,再根据所述旋转矩阵和位移量以及重建参数就可以重建出CT图像。即使对于大型受检物体,通过本发明的方法也能够得到每次扫描后的精确几何关系从而使得对超大物体的CT检测成为可能。
再者,即使受检物体的运动受到外界因素干扰,仍能通过本发明CT系统和CT图像重建方法得到准确的运动轨迹并实现精确重建,因此本发明CT系统的抗干扰性强。
附图说明
在下文中,将结合下列图形对本发明进行描述,其中,同样的附图标记表明同样的元素,在附图中:
图1示出了依据本发明CT系统的一个实施例的简化结构示意图;
图2示出了依据本发明CT图像重建方法的一个实施例的流程图;
图3为以柱面透镜作为镜头的线阵CCD相机的成像原理示意图;
图4为柱面透镜相对于标记点摆放的位置示意图;
图5示出了标记组的形态及方向向量定义的示意图;
图6为以待重建区域中心位置为坐标原点建立的坐标系的示意图;
图7是依据本发明的一个实施例进行多段直线CT扫描的简化布局;
图8是直线轨迹滤波反投影算法示意图;
图9是采用本发明的技术方案和采用传统方案分别得到的CT重建图像以及原始模型的示意图。
具体实施方式
下列详细说明实际上仅仅是示例性的,并非意图限制本发明或者本发明的应用和使用。如本文所使用的,词“示例性”意指“用作示例、例子或示范”。因此,本文被描述为“示例性”的任何实施例不必被解释为比其它实施例优选或有利。本文所描述的所有实施例都是示例性实施例,其被提供用来使得本领域技术人员能够制作或使用本发明,而并非限制本发明的范围,本发明的范围由权利要求来限定。
如图1所示,其示出了依据本发明的一个实施例的CT系统的简化结构示意图。在图1中,该示例性实施例所图示的CT系统包括射线源1、射线探测器2、数据采集与处理单元3、主控制器4,其中,所述射线探测器2探测所述射线源1发射的X射线并将其转换成数字信号以传送给所述数据采集与处理单元3;所述射线源1和所述数据采集与处理单元3由所述主控制器4控制。另外,本发明CT系统还包括相机组5和标记组7以及标记组处理单元6,其中:
所述标记组7被固定在受检物体表面;所述相机组5用来获得所述标记组7在其上的投影;所述标记组处理单元6根据投影得出所述标记组7的初始位置以及旋转后的位置信息,从而计算得到受检物体的旋转矩阵和位移量。
对于该旋转矩阵和位移量(也可称为平移量),它们是用来在三维解析几何、计算机视觉以及定位领域中描述物体位置和方向信息的常用量。旋转矩阵一般指在空间内某点旋转前后所在坐标系的对应关系。假设有一固定模型,该模型是一系列点坐标a的集合。如果对该模型进行平移则模型内所有点坐标满足a’=a+a0, 其中a0为模型的位移量,a表示平移前各点的坐标,a’为平移后各点的坐标;如果对整个模型进行旋转,则模型内所有点的坐标满足a’=aM,其中M例如可以为3×3的旋转矩阵,根据xyz轴各个方向旋转角度的大小可以唯一确定旋转矩阵中的各元素的数值。综上可知物体在三维空间内的运动引起的坐标变化可以由a’= aM+a0表示。在本文中,旋转矩阵指受检物体初始位置所在坐标系与旋转后所在坐标系的对应关系。
其中,所述相机组5包括以柱面透镜作为镜头的线阵CCD相机。
另外,所述相机组5也可以包括以球面透镜作为镜头的面阵CCD相机。
对于所述标记组7中的标记点可以由单色反光材料制成或者具有单色反光材料的涂层。在所述相机组5中的相机的镜头上添加有对应颜色的滤光片。也就是说,如果标记点由绿色反光材料制成,那么在相机的镜头上应添加绿颜色的滤光片,这样可以滤去其它颜色的干扰。
此外,所述标记组7至少包括三个标记点。
其中,所述CT系统可以包括多个相机组5且这多个相机组5的视野范围能够覆盖在扫描过程中所述标记组7中的所有标记点可能出现的区域。
此外,所述相机组5可以包括至少三个相机且这些相机分布在以被监视区域中心为中心的圆弧上。以所述中心为坐标原点建立坐标系,如果在xy平面内以被监视区域中心为中心的圆弧上分布了三个相机,则左右两个相机的镜头例如顺着z轴方向摆放而中间的相机的镜头平行于xy平面方向摆放。
对于标记组处理单元6,其首先根据标记组7在相机组5上的投影得到标记组7的初始位置(此时,受检物体也处于初始位置,即受检物体未被旋转),然后再根据受检物体旋转后的标记组7在相机组5上的投影得到旋转后标记组7的位置信息。通过标记组7的初始位置和旋转后的位置信息即可得到受检物体的旋转矩阵和位移量,从而根据所得到的旋转矩阵和位移量,采用任何已知的CT图像重建方法就可以得到受检物体的CT重建图像。
下面,介绍一个柱面透镜的成像原理。如图3所示,当标记点71到线阵CCD相机的镜头的距离远远大于焦距时,其在焦平面位置的线阵CCD相机上的投影将是一个点,设该点的位置为z,对于线阵CCD,其输出为与单位时间内成像单元上光通量成正比的离散一维信号,通过寻找该信号峰值的位置便可确定投影点的位置,该投影与柱面透镜光心线和标记点71在同一平面内。对于标记组7中单个标记点71的定位过程,设柱面透镜的焦距为f,线阵CCD相机中心位置为z-0,线阵CCD相机中心位置与柱面透镜光心线所在的平面定义为中心平面,标记点71所在平面与中心平面的夹角为θ,则有: 。设中心平面法向量为n,柱面透镜光心线方向向量q,m=n×q,则标记点所在面的法向量n’为:n’=ncosθ-msinθ;设相机镜头中心位置为C0,标记点71位置为Ci,则所确定的平面方程为:;对于一个相机组5中三个不同的相机所确定出的三个不同的平面,根据三个平面确定一点的原理则有:
其中,C01,C02, C03分别为三个相机的镜头中心位置,n1’, n2’, n3’为每个镜头确定出的该标记点71所在平面的法向量,解此线性方程组即可得到标记点71的位置Ci。同理可以得到标记组7中其它标记点71的位置。
对于标记组7的位置信息,例如此标记组7具有三个标记点71,则其位置可以由四个向量来表示,分别是标记组中心位置C和三个描述标记组朝向的单位向量u,v,w,如图5所示,其示出了标记组7的形态及方向向量定义的示意图。在图5中,标记组7由三个标记点组成,设三个标记点的空间坐标分别为c1,c2,c-3,则C=(c1+c2+c3)/3,u=(c1-C)/|c1-C|,v=(c2-C)/|c2-C|,w=u×v,受检物体移动或旋转后标记组7的位置信息记为C’,u’,v’,w’。同理,C’,u’,v’,w’也可以根据上面描述计算得出。
相应地,如图2所示,其示出了依据本发明CT图像重建方法的一个实施例的流程图,本发明CT图像重建方法包括步骤:
1)将标记组7固定在受检物体表面并通过相机组5和标记组处理单元6得到所述标记组7的初始位置;
2)在每次受检物体被旋转到某位置后且在该某位置对受检物体进行扫描前通过相机组5和标记组处理单元6获得旋转后标记组7的位置信息;
3)根据所述标记组7的初始位置和旋转后的位置信息获得受检物体的旋转矩阵和位移量;
4)根据所述旋转矩阵和位移量以及重建参数来重建CT图像。
从上述可知,本发明CT图像重建方法通过将标记组7固定在受检物体表面并通过相机组5和标记组处理单元6得到标记组7的初始位置;在每次受检物体旋转后,也即标记组7旋转后再获得旋转后标记组7的位置信息。根据标记组7的初始位置和旋转后的位置信息即可得到标记组7的旋转矩阵和位移量,也就是受检物体的旋转矩阵和位移量,这是因为标记组7是固定在受检物体表面上的,这样标记组7就随着受检物体的旋转而旋转,所以通过比较标记组7旋转后的位置信息与标记组7的初始位置之间的旋转关系以及位移关系就可以得到受检物体旋转后的位置与受检物体初始位置之间的旋转关系以及位移关系。接着根据所述旋转矩阵和位移量以及重建参数等来重建CT图像。
其中,所述相机组5包括以柱面透镜作为镜头的线阵CCD相机。所述相机组5还可以包括以球面透镜作为镜头的面阵CCD相机。
其中,所述标记组7中的标记点由单色反光材料制成或者具有单色反光材料的涂层。相应地,在相机组5中的相机的镜头上添加有对应颜色的滤光片,这样可以滤去其它颜色的干扰。
经实验证明,所述标记组7至少包括三个标记点。
优选地,所述相机组5的数目为两个以上且这些相机组的视野范围能够覆盖在扫描过程中所述标记组7中的所有标记点可能出现的区域。
此外,所述相机组5包括至少三个相机且这些相机分布在以被监视区域中心为中心的圆弧上。
以所述中心为坐标原点建立坐标系,如果在xy平面内以被监视区域中心为中心的圆弧上分布了三个相机,则左右两个相机的镜头顺着z轴方向摆放而中间的相机的镜头平行于xy平面方向摆放。
下面以采用以柱面透镜作为镜头的线阵CCD相机组成的相机组5为例来详述本发明CT图像重建方法的技术方案,其中每一个相机组5由三个以柱面透镜作为镜头的线阵CCD相机,采用两个标记组(如图7所示),每个标记组包括三个标记点。
对于每一个线阵CCD相机组5的摆放,如图4所示,在该示例中以被监视区域中心o作为坐标原点建立坐标系,把三个线阵CCD相机51、52、53分别布置在半径为R的圆弧上,其中左右两个线阵CCD相机51和53的镜头顺坐标系的z轴方向摆放,中间的CCD相机的柱面镜头平行于xy平面方向摆放。关于线阵CCD相机51、52和53的视角、到目标区域的距离、在圆弧上的位置等参数可以根据实际系统的尺寸和精度要求而具体设定,例如当相机组5被用来监视1.5m×1.5m的方形区域内标记点71的位置,则可以把视角为30度的镜头摆放在距目标区域中心3m的位置,每个镜头之间在圆弧上相隔45度,来实现对标记点的拍摄以便确定其位置信息。
下面以待重建区域中心位置为坐标原点建立坐标系,如图6所示,把此坐标系定义为系统坐标系。
将标记组5固定在受检物体表面之后,通过相机组5和标记组处理单元6得到所述标记组7的初始位置,即(C,u,v,w),其中C为标记组中心位置坐标,u,v,w为描述标记组朝向的单位向量,即三个标记点的单位方向向量。受检物体移动或旋转后标记组7的位置信息记为(C’,u’,v’,w’)。
标记组方向向量构成的坐标系(图6中以C为原点uvw为坐标轴的坐标系)与系统坐标系(如图6所示)之间的对应关系为矩阵R,该矩阵R满足关系:
受检物体旋转后,也即标记组7旋转后的方向向量所构成的坐标系(图6中以C’为原点u’v’w’为坐标轴的坐标系)与系统坐标系(如图6所示)之间的对应关系为矩阵R’,该矩阵R’满足关系:
则对于受检物体内的任意点,设其初始位置在系统坐标系下的坐标为a,以及在uvw坐标系下的坐标为a0,旋转后在系统坐标系下的坐标为a’,则满足如下公式:
(a-C)R=a0
(a’-C’)R’=a0
通过上述公式可以得到坐标a’和a的关系为:
a’=aRR’-1+C’-CRR’-1,
得到了受检物体的旋转矩阵和位移量之后,就可以根据该旋转矩阵和位移量,采用现有技术已知重建方法来重建CT图像。众所周知,现在存在很多种方法来重建CT图像,例如代数迭代法、滤波反投影法等等。
对于代数迭代方法,其公式如下:
其中f为重建图像的向量表示,g表示射线对应的投影值,H={hij}对应于每条射线与重建物体图像的几何关系,满足,其中每一个元素hij表示第i条射线穿过图像中第j个像素的长度。对于射线源1所发射的每条X射线,我们把其几何关系的定义为由射线源位置p和单位方向向量r来表示,也即每条射线由(p,r)来表示。则hij的数值可由像素j的位置和射线i的参数(p,r)来求解,即由(p,r)来确定直线方程并求其穿过像素j的长度,k表示迭代次数。
设第j条射线的参数为(rj,p)的射线,其中rj=(rx,ry,rz),p=(px,py,pz),由于每个像素的形状为正立方体,设第i个像素所表示的正立方体六个面分别在x0; x1; y0; y1; z0; z1位置,根据正方体的六个面定义该像素,则第i个像素的范围为x0<x<x1,y0<y<y1,z0<z<z1。则第i条射线所在直线方程为:x=px+t×rx; y=py+t×ry; z=pz+t×rz;其中t为直线参数。将z0带入直线方程得t=(z- pz)/rz;x= px+(z- pz)×rx/rz;y= py+(z- pz)×ry/ rz,判断所得的x,y是否满足x0<x<x1,y0<y<y1,如果满足则记ti=t,(i=1,2为符合该条件的t的编号),以此类推分别带入z1, x0,x1,y0,y1并找出符合条件的t1和t2,如果能找到两个符合条件的t,则hij=|t1-t2|,否则hij=0。
p’=p0R’R-1+C-C’RR’-1;
r’=r0R’R-1;
其中p0,r0为初始位置射线源1在系统坐标系中的位置和方向。
Hi=[hi1…hin]则表示第i条X射线的投影与所经过的像素的关系,满足gi=Hif,其中n为图像的总像素数,gi为第i条射线的投影值。代数迭代算法的核心思想就是通过解方程g=Hf中的f来得到断层图像。对于扫描得到的投影数据,按照以上方法对每条数据的射线进行迭代就可以重建出目标断层的图像。
对于滤波反投影方法,通常在机械条件能使扫描轨迹比较规律的情况下采用此方法进行CT断层图像重建。以多段直线CT扫描为例,如图7所示,扫描过程首先是射线源1和射线探测器2先相对于受检物体做直线轨迹扫描,之后再旋转受检物体,在下一个角度下进行直线CT扫描。在本例中配有两个相机组5和5’,这两个相机组5和5’的视野范围的间隔由预计的旋转角度决定。以三段直线CT扫描为例,射线源1的张角为60度,每次直线扫描后把受检物体旋转60度,两个相机组5和5’分别用来监视受检物体表面相隔60度的两个区域。初始标记组7放置在相机组5的视野范围内,旋转60度后标记组7出现在相机组5’的视野范围内,根据两个相机组5和5’确定的标记组7旋转前后的位置差异来得到受检测物体相对于初始位置的旋转矩阵和位移量,由旋转矩阵Ri和位移量Ci表示。
对于多段直线CT扫描,定义每一段扫描得到的投影数据为gi(l,t),如图8所示,l表示射线源和探测器在轨道上的位置,t表示探测器单元在探测器上的位置。对于每一段投影数据gi,重建出对应的图像fi,其公式为:
其中,p(l,t)为扫描所得到的投影数据,D表示射线源到探测器的距离,以重建图像中心为原点建立系统坐标系,z轴垂直于重建图像所在平面,其中a=(x,y,0)表示重建图像中像素的位置,l’=x-t/2-t×y/D表示t位置的探测器单元过(x,y)位置的像素所对应的射线源位置,h(l)为对应的SL或RL滤波函数。根据得到的fi结合所得到的受检物体的旋转矩阵和位移量,在此基础上对每组数据进行图像融合,所得到的重建图像f为:
如图9所示,其图示了受检物体的原始模型与采用本发明和非本发明方法得到的该受检物体的CT重建图像。在该图9中,中间一幅为,采用传统方法得到的CT重建图像;最左边为受检物体的原始模型,最右边为采用本发明CT系统和本发明CT图像重建方法得到的受检物体的CT重建图像。从此图9中可以看出,采用本发明CT图像重建方法得到的受检物体的CT重建图像更接近于该受检物体的原始模型。
尽管在本发明的上述详细说明中已经呈现了至少一个示例性实施例,但是应该理解存在大量的变体。还应当理解一个或多个示例性实施例仅是示例,并且不以任何方式限制本发明的范围、适用性或配置。相反,上述详细说明将为本领域技术人员提供便捷路线图以用于实施本发明的示例性实施例。应当理解在不脱离所附权利要求中所阐述的本发明范围的情况下,可以在示例性实施例中所描述的元素的功能和布置方面进行各种变化。
Claims (18)
1. 一种CT系统,其包括射线源、射线探测器、数据采集与处理单元、主控制器,其中,所述射线探测器探测所述射线源发射的X射线并将其转换成数字信号以传送给所述数据采集与处理单元;所述射线源和所述数据采集与处理单元由所述主控制器控制,其特征在于,所述CT系统还包括相机组和标记组以及标记组处理单元,其中:
所述标记组被固定在受检物体表面;
所述相机组用来获得所述标记组在其上的投影;
所述标记组处理单元根据所述投影得出所述标记组的初始位置以及旋转后的位置信息,从而计算得到受检物体的旋转矩阵和位移量。
2. 如权利要求1所述的CT系统,其中,所述相机组包括以柱面透镜作为镜头的线阵CCD相机。
3. 如权利要求1所述的CT系统,其中,所述相机组包括以球面透镜作为镜头的面阵CCD相机。
4. 如权利要求2或3所述的CT系统,其中,所述标记组中的标记点由单色反光材料制成或者具有单色反光材料的涂层。
5. 如权利要求4所述的CT系统,其中,所述相机组中的相机的镜头上添加有对应颜色的滤光片。
6. 如权利要求5所述的CT系统,其中,所述标记组至少包括三个标记点。
7. 如权利要求6所述的CT系统,其中,所述CT系统可以包括多个相机组且这多个相机组的视野范围能够覆盖在扫描过程中所述标记组中的所有标记点可能出现的区域。
8. 如权利要求7所述的CT系统,其中,所述相机组包括至少三个相机且这些相机分布在以被监视区域中心为中心的圆弧上。
9. 如权利要求8所述的CT系统,其中,以所述中心为坐标原点建立坐标系,如果在xy平面内以被监视区域中心为中心的圆弧上分布了三个相机,则左右两个相机的镜头顺着z轴方向摆放而中间的相机的镜头平行于xy平面方向摆放。
10. 一种CT图像重建方法,其特征在于,包括步骤:
将标记组固定在受检物体表面并通过相机组和标记组处理单元得到所述标记组的初始位置;
在每次受检物体被旋转到某位置后且在该某位置对受检物体进行扫描前通过相机组和标记组处理单元获得旋转后标记组的位置信息;
根据所述标记组的初始位置和旋转后的位置信息获得受检物体的旋转矩阵和位移量;
根据所述旋转矩阵和位移量以及重建参数来重建CT图像。
11. 如权利要求10所述的CT图像重建方法,其中,所述相机组包括以柱面透镜作为镜头的线阵CCD相机。
12. 如权利要求10所述的CT图像重建方法,其中,所述相机组包括以球面透镜作为镜头的面阵CCD相机。
13. 如权利要求10至12所述的CT图像重建方法,其中,所述标记组中的标记点由单色反光材料制成或者具有单色反光材料的涂层。
14. 如权利要求13所述的CT图像重建方法,其中,所述相机组中的相机的镜头上添加有对应颜色的滤光片。
15. 如权利要求14所述的CT图像重建方法,其中,所述标记组至少包括三个标记点。
16. 如权利要求15所述的CT图像重建方法,其中,所述相机组的数目为两个以上且这些相机组的视野范围能够覆盖在扫描过程中所述标记组中的所有标记点可能出现的区域。
17. 如权利要求16所述的CT图像重建方法,其中,所述相机组包括至少三个相机且这些相机分布在以被监视区域中心为中心的圆弧上。
18. 如权利要求17所述的CT图像重建方法,其中,以所述中心为坐标原点建立坐标系,如果在xy平面内以被监视区域中心为中心的圆弧上分布了三个相机,则左右两个相机的镜头顺着z轴方向摆放而中间的相机的镜头平行于xy平面方向摆放。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110292773.8A CN103033525B (zh) | 2011-09-30 | 2011-09-30 | Ct系统和ct图像重建方法 |
HK13109732.9A HK1182446A1 (zh) | 2011-09-30 | 2013-08-20 | 系統和 圖像重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110292773.8A CN103033525B (zh) | 2011-09-30 | 2011-09-30 | Ct系统和ct图像重建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103033525A true CN103033525A (zh) | 2013-04-10 |
CN103033525B CN103033525B (zh) | 2016-03-02 |
Family
ID=48020640
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110292773.8A Active CN103033525B (zh) | 2011-09-30 | 2011-09-30 | Ct系统和ct图像重建方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN103033525B (zh) |
HK (1) | HK1182446A1 (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103776851A (zh) * | 2014-02-07 | 2014-05-07 | 中国科学院高能物理研究所 | 一种智能化ct扫描的方法与装置 |
CN104142353A (zh) * | 2014-06-06 | 2014-11-12 | 清华大学 | 可见光成像辅助定位局部ct扫描方法及系统 |
CN104422704A (zh) * | 2013-08-21 | 2015-03-18 | 同方威视技术股份有限公司 | 对x光能谱ct的能谱信息进行分解的方法和对应的重建方法 |
CN104749197A (zh) * | 2013-12-26 | 2015-07-01 | 清华大学 | Ct系统及其方法 |
CN104757988A (zh) * | 2015-05-04 | 2015-07-08 | 重庆大学 | 一种电子直线扫描微纳焦点ct扫描系统及方法 |
CN107913749A (zh) * | 2017-12-05 | 2018-04-17 | 佛山科学技术学院 | 一种试管结构 |
CN110595999A (zh) * | 2018-05-25 | 2019-12-20 | 上海翌视信息技术有限公司 | 一种图像采集系统 |
CN111707687A (zh) * | 2020-03-26 | 2020-09-25 | 国网河南省电力公司电力科学研究院 | 一种直线轨迹扫描ct检测方法及装置 |
WO2022041092A1 (zh) * | 2020-08-26 | 2022-03-03 | 深圳先进技术研究院 | Ct成像系统及其成像方法 |
CN114152637A (zh) * | 2022-02-07 | 2022-03-08 | 东莞市志橙半导体材料有限公司 | 一种硬质碳化硅材料打孔检测装置与方法 |
CN116196022A (zh) * | 2023-04-28 | 2023-06-02 | 之江实验室 | 扇形x光束穿过介质时的通量分布计算方法和系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1991001514A1 (en) * | 1989-07-14 | 1991-02-07 | Raychem Limited | Laser machining |
US5023894A (en) * | 1989-02-07 | 1991-06-11 | Hamamatsu Photonics K.K. | Data correcting device for use with computed tomography apparatus |
JPH1119080A (ja) * | 1997-07-08 | 1999-01-26 | Shimadzu Corp | X線ct装置 |
US6006126A (en) * | 1991-01-28 | 1999-12-21 | Cosman; Eric R. | System and method for stereotactic registration of image scan data |
US6167295A (en) * | 1991-01-28 | 2000-12-26 | Radionics, Inc. | Optical and computer graphic stereotactic localizer |
EP1642545A1 (en) * | 1997-12-01 | 2006-04-05 | Sherwood Services AG | Surgical positioning system |
CN101776618A (zh) * | 2009-01-08 | 2010-07-14 | 欧姆龙株式会社 | 检查方法、检查装置以及检查用程序 |
WO2011047467A1 (en) * | 2009-10-20 | 2011-04-28 | Imris Inc. | Imaging system using markers |
CN102245106A (zh) * | 2008-12-11 | 2011-11-16 | 皇家飞利浦电子股份有限公司 | 用于产生患者内部和外部图像的系统和方法 |
CN202305443U (zh) * | 2011-09-30 | 2012-07-04 | 清华大学 | Ct系统 |
-
2011
- 2011-09-30 CN CN201110292773.8A patent/CN103033525B/zh active Active
-
2013
- 2013-08-20 HK HK13109732.9A patent/HK1182446A1/zh not_active IP Right Cessation
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5023894A (en) * | 1989-02-07 | 1991-06-11 | Hamamatsu Photonics K.K. | Data correcting device for use with computed tomography apparatus |
WO1991001514A1 (en) * | 1989-07-14 | 1991-02-07 | Raychem Limited | Laser machining |
US6006126A (en) * | 1991-01-28 | 1999-12-21 | Cosman; Eric R. | System and method for stereotactic registration of image scan data |
US6167295A (en) * | 1991-01-28 | 2000-12-26 | Radionics, Inc. | Optical and computer graphic stereotactic localizer |
JPH1119080A (ja) * | 1997-07-08 | 1999-01-26 | Shimadzu Corp | X線ct装置 |
EP1642545A1 (en) * | 1997-12-01 | 2006-04-05 | Sherwood Services AG | Surgical positioning system |
CN102245106A (zh) * | 2008-12-11 | 2011-11-16 | 皇家飞利浦电子股份有限公司 | 用于产生患者内部和外部图像的系统和方法 |
CN101776618A (zh) * | 2009-01-08 | 2010-07-14 | 欧姆龙株式会社 | 检查方法、检查装置以及检查用程序 |
WO2011047467A1 (en) * | 2009-10-20 | 2011-04-28 | Imris Inc. | Imaging system using markers |
CN202305443U (zh) * | 2011-09-30 | 2012-07-04 | 清华大学 | Ct系统 |
Non-Patent Citations (2)
Title |
---|
吕洪波等: "基于虚拟现实技术的机器人外科手术模拟与培训系统研究", 《高技术通讯》, 30 November 2001 (2001-11-30), pages 88 - 92 * |
李天元等: "CT金属物标记定位法在高血压脑出血微创术中的应用", 《中国医药指南》, vol. 6, no. 2, 29 February 2008 (2008-02-29), pages 129 - 130 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104422704A (zh) * | 2013-08-21 | 2015-03-18 | 同方威视技术股份有限公司 | 对x光能谱ct的能谱信息进行分解的方法和对应的重建方法 |
CN104422704B (zh) * | 2013-08-21 | 2017-07-25 | 同方威视技术股份有限公司 | 对x光能谱ct的能谱信息进行分解的方法和对应的重建方法 |
US9786041B2 (en) | 2013-12-26 | 2017-10-10 | Tsinghua University | CT systems and methods thereof |
CN104749197A (zh) * | 2013-12-26 | 2015-07-01 | 清华大学 | Ct系统及其方法 |
CN104749197B (zh) * | 2013-12-26 | 2017-08-11 | 清华大学 | Ct系统及其方法 |
CN103776851A (zh) * | 2014-02-07 | 2014-05-07 | 中国科学院高能物理研究所 | 一种智能化ct扫描的方法与装置 |
CN103776851B (zh) * | 2014-02-07 | 2016-08-17 | 中国科学院高能物理研究所 | 一种智能化ct扫描的方法与装置 |
CN104142353A (zh) * | 2014-06-06 | 2014-11-12 | 清华大学 | 可见光成像辅助定位局部ct扫描方法及系统 |
CN104142353B (zh) * | 2014-06-06 | 2017-01-11 | 清华大学 | 可见光成像辅助定位局部ct扫描方法及系统 |
CN104757988A (zh) * | 2015-05-04 | 2015-07-08 | 重庆大学 | 一种电子直线扫描微纳焦点ct扫描系统及方法 |
CN107913749A (zh) * | 2017-12-05 | 2018-04-17 | 佛山科学技术学院 | 一种试管结构 |
CN110595999A (zh) * | 2018-05-25 | 2019-12-20 | 上海翌视信息技术有限公司 | 一种图像采集系统 |
CN111707687A (zh) * | 2020-03-26 | 2020-09-25 | 国网河南省电力公司电力科学研究院 | 一种直线轨迹扫描ct检测方法及装置 |
WO2022041092A1 (zh) * | 2020-08-26 | 2022-03-03 | 深圳先进技术研究院 | Ct成像系统及其成像方法 |
CN114152637A (zh) * | 2022-02-07 | 2022-03-08 | 东莞市志橙半导体材料有限公司 | 一种硬质碳化硅材料打孔检测装置与方法 |
CN116196022A (zh) * | 2023-04-28 | 2023-06-02 | 之江实验室 | 扇形x光束穿过介质时的通量分布计算方法和系统 |
CN116196022B (zh) * | 2023-04-28 | 2023-08-04 | 之江实验室 | 扇形x光束穿过介质时的通量分布计算方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103033525B (zh) | 2016-03-02 |
HK1182446A1 (zh) | 2013-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103033525B (zh) | Ct系统和ct图像重建方法 | |
CN103759669B (zh) | 一种大型零件的单目视觉测量方法 | |
CN102971604A (zh) | 用于确定车辆车轮对准的系统和相关方法 | |
CN105913417B (zh) | 基于透视投影直线的几何约束位姿方法 | |
CN106780609B (zh) | 视觉定位方法和视觉定位装置 | |
CN104567758B (zh) | 立体成像系统及其方法 | |
CN105043350A (zh) | 一种双目视觉测量方法 | |
CN202305443U (zh) | Ct系统 | |
KR20100087083A (ko) | 물체의 모양을 삼차원 측정하기 위한 시스템 및 방법 | |
CN106408609A (zh) | 一种基于双目视觉的并联机构末端运动位姿检测方法 | |
CN110849331B (zh) | 基于三维点云数据库模型的单目视觉测量与地面试验方法 | |
CN104330074A (zh) | 一种智能测绘平台及其实现方法 | |
CN103196370A (zh) | 一种导管接头空间位姿参数的测量方法和装置 | |
Xia et al. | Global calibration of non-overlapping cameras: State of the art | |
CN104036542A (zh) | 一种基于空间光线聚集性的像面特征点匹配方法 | |
CN108022267A (zh) | 一种摄像机标定装置及标定方法 | |
CN110136047B (zh) | 一种车载单目图像中静止目标三维信息获取方法 | |
CN105574812A (zh) | 多角度三维数据配准方法及装置 | |
Shang et al. | Single-pass inline pipeline 3D reconstruction using depth camera array | |
Mei et al. | Monocular vision for pose estimation in space based on cone projection | |
CN102881040A (zh) | 一种数码相机移动拍摄三维重建方法 | |
Chen et al. | Drogue pose estimation for unmanned aerial vehicle autonomous aerial refueling system based on infrared vision sensor | |
CN117553697A (zh) | 一种基于led的高速摄像测量方法及舱门形变测量系统 | |
CN105931292A (zh) | 一种基于仿射标定的多方向莫尔层析方法 | |
WO2023159180A1 (en) | Single-pass 3d reconstruction of internal surface of pipelines using depth camera array |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 1182446 Country of ref document: HK |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: GR Ref document number: 1182446 Country of ref document: HK |