CN1464301A - X射线立体摄影实现图像重建的ct-解析法 - Google Patents
X射线立体摄影实现图像重建的ct-解析法 Download PDFInfo
- Publication number
- CN1464301A CN1464301A CN 02133358 CN02133358A CN1464301A CN 1464301 A CN1464301 A CN 1464301A CN 02133358 CN02133358 CN 02133358 CN 02133358 A CN02133358 A CN 02133358A CN 1464301 A CN1464301 A CN 1464301A
- Authority
- CN
- China
- Prior art keywords
- image
- ray
- coordinate
- radon
- stereoscopic photograph
- 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.)
- Pending
Links
Images
Landscapes
- Image Processing (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明属计算机图像重建和数字影像处理技术领域。X射线立体摄影获得多幅数字影像,用计算机视觉相关匹配和ShuJia-Radon关系转换,使X射线光源坐标、被摄影物体坐标和影像的像片坐标被统一到SR坐标系中,光源与像素格网灰度单元构成方向射线矢量方程,实现Radon原理图像重建。能同时实现三种影像观察功能;二维平片观察、全仿真立体影像观察、断层CT影像观察。简化了精密扫描轨道和电子探测器等。
Description
1、技术领域:
本发明属计算机图像重建、数字图像处理技术领域。
2、技术背景:
1、张祖勋合编.数字摄影测量学.武汉大学出版社,1997.1
2、夏良正主编.数字图像处理.东南大学出版社,1999.8
3、康晓东编著.现代医学影像技术.天津科技翻译出版公司,2000.1
4、余建明主编.数字减影血管造影技术.人民军医出版社,1999.5
3、发明的目的:
用X射线立体摄获得面状数字影像,用解析计算的方法实现Radon原理计算机图像重建,获得断层CT图像,提高计算机图像重建的精度和速度,减小X射线辐射量,X射线立体摄不需要精密的摄影轨道和排列的探测器。在CT-解析法设备中同时实现三种影像观察功能;二维平片观察、全仿真立体影像观察、断层CT影像观察,并输出三种影像。
4、发明内容:
一、CT-解析法的影像摄影条件
按人的立体视觉原理和立体影像相关匹配原理,观察立体影像必须获得二张对同一物体在不同位置拍摄的中心投影重叠像片。当用左右眼分别观看这两张像片影像时,就会在人脑中产生立体视觉。人的立体视觉符合三维立体几何交会原理,能建立起严密的数学关系模型,在立体观察中能采集、测绘出物体影像的三维空间坐标数据,用计算机能对采集的物体的三维空间影像数据进行计算处理。观察到物体内部结构密度形态的全仿真立体影像。
二、X射线立体影像的获得
按现有常规的X光机设备分为小、中、大型X光机,C臂X光机,在这些常规的X光机设备中,X射线的发光源球管F都只设置一个。用小型X光机一次只能拍单张X光片,中、大型X光机、C臂X光机能拍出连续的X光片或数字X光影像,这些X光机设备都能按立体视觉观察要求拍摄出重叠的立体观察影像。如图(1)、(2)、(3)。图(1)是用小型X光机分步移动拍单张X光片的方法获得重叠的立体像对影像。图(2)是用中、大型X光机连续移动拍多张X光片的方法获得多张重叠立体像对影像。图(3)是C臂X光机用连续拍多张X光片的方法获得环形重叠立体像对影像。
专用的计算机X射线数字立体影像观察CT-解析法的设备设有两个X射线的发光源球管F1,F2,如图(4),两个X发光源球管F1,F2交替开闭发射X射线透射物体,能实时连续不断的获得X光重叠立体像对影像。
在检查物体、货物的X光检查过程中,输送带匀速使物体运动,如图(5)X射线球管按一定的时间间隔频率开闭发射X射线,数字化影像转换屏不断得到物体的X射线立体像对影像。
为了减少伴影,得到清晰的X射线影像,应采用微焦点X射线球管。在X射线立体摄影时,在设备上设定已知位置参数点,作为坐标量测的起算定位像点,确定影像的比例尺和坐标起算位置。
三、X射线影像的显示和全仿真立体观察
X射线影像的显示和立体观察是为了使我们能从计算机显示屏上直接观察到物体内部结构密度形态立体影像,如图(20)是一个典型的CT-解析法设备;(图20中;1、X射线球管 2、物体 3、影像增强器 4、CCD 5、A/D 6、计算机7、偏振光立体眼镜)
在CT-解析法设备中,两个X射线球管F1,F2交替透射物体,影像增强器不断得到左右两张立体像片的光学模拟影像,并将物体X光影像亮度增强几千倍。CCD像机从影像增强器拍摄到增强的影像后,将光学影像数字化后传到计算机进行数字影像纹理处理,对影像灰度进行分类,检测锐化边缘轮廓影像,突出兴趣影像,提高影像灰度进行显示。计算机虚拟立体观察专业图像显示卡和图像操作软件以偏振光形式在计算机显示屏上交替显示X射线立体像对影像,供立体影像观察和量测空间坐标。观测X射线立体影像有两种方法;
a、偏振光观察法,观察时须佩戴偏振光立体眼镜,偏振光立体眼镜的偏振频率由显示卡控制的红外发射器同步控制,使观察者在观察显示屏时只能左眼看左X射线影像,右眼看右X射线影像,在观察者大脑中形成立体视觉,可获得很高的观察精度。
b、边光显示器和3D立体显示器,观察立体时不需佩戴立体眼镜,但观察立体影像精度不如偏振光观察法高,价格也贵。
在显示器上观察立体的特点是视场大,对观察大于1∶1的X射线立体影像很有利。
CT-解析法也可以按常规X光平片的形式在计算机显示屏上显示X射线二维影像,供常规平面X光影像观察。
四、透射光立体摄影数学模型
用透射光按立体摄影条件摄影,能建立透射光摄影立体影像严密空间几何交会数学模型。使立体摄影测量的方法和技术手段,可用于透射光摄影立体影像观察和影像处理。
a、X射线透射成像原理
X射线的光源X射线球管近似中心投影的点光源,当X射线透射物体后,物体因内部的密度不同,在X射线成像面P上产生不同深浅的影像灰度,显示物体内部结构密度形状影像。成像原理如图(7),图中F为X射线球管光源,P是X射线成影面,H是光源F到成影面P的距离,o是F中心光线垂直投影到成影面P的交点,称o为主点。在P面上的成像是空间物体A内部结构密度形态的X射线中心投影的二维平面影像。
当在成影面P上加一个透镜S,如用CCD像机对X片P摄影,如图(8),能得到需要大小尺寸的X射线影像像片P′。o′S为像距f,SP为物距H。
b、X射线立体摄影像对原理
为了观察到X射线立体影像,须得到能建立物体立体影像的两张“像对”像片或直接数字化记录下像对影像信息,用计算机虚拟现实可视化立体技术观察X射线立体影像。
如图(9),设有两个X射线球管F1,F2光源,并相隔一段距离M。F1,F2垂直投射在两像片P1,P2上的交点o1,o2称为左主点和右主点,是P1,P2像片上的坐标原点。当X射线透射物体后,物体A,B在P1,P2上分别的成像为a1,b1和a2,b2,物体A,B在P1,P2上的投影点a1与a2,b1与b2的距离之差为投影差,相似于人的生理视差:
Δx1=a1-a2;Δx2=b1-b2;
既;Δp=Δx1-Δx2;
当获得立体像对P1,P2像片后,根据立体观察原理,两眼同时分别观看P1,P2两像片上的像点a1,b1和a2,b2将看到物体A,B的立体视觉影像。
用透镜的方法如图(8),同样可建立图(9)的关系,物体A,B在P1′,P2′上的成像为a1′,b1′和a2′,b2′,从图中看出在P1′,P2′像片上并未改变X射线透射成像光线的几何条件,当立体观察P1′,P2′像片时,同样能看到物体A,B的立体视觉影像,但立体视觉影像被缩小
c、X射线立体摄影数学模型
分别在两张X射线成影像片P1,P2建立像片平面坐标系,如图(10),o1-x1y1和o2-x2y2,坐标原点为o1,o2,而o1,o2是X射线球管光源F1,F2垂直辐射在P1,P2像片面上的交点。
当得到X射线立体观察像对P1,P2像片后,建立物体与X射线成影像点之间的空间几何交会数学关系。如图(11),H为X射线光源F1,F2至成影像片P1,P2的设定常数距离,为摄影高度。M为两X光源F1、F2之间设定的常数距离,为摄影基线。设物体的空间坐标系X,Y轴与成影像片P1的像片坐标x1,y1轴平行,Z轴垂直于像片P1坐标原点o1,物体的空间坐标系原点O与像片P1平面坐标系原点o1的垂直距离为ΔH。
当两X射线光源F1,F2透射物体A后,空间物体A在P1,P2像对像片上的同名影像点分别为a1,a2,其a1在P1片上的坐标为x1,y1,a2在P2片上的坐标为x2,y2。可设定y1=y2。
依据立体观察原理,X坐标轴平行摄影基线M和像片坐标x1,x2轴。求在两像片P1,P2上的同名像点a1,a2的x坐标之差为Δx=x1-x2,Δx称为像点左右视差,是立体视觉的来源。
在图(11)中,从O作辅助线至A″M得出,当已设定X射线光源F1,F2的固定位置和基线距离M,及F1,F2到成影面P1,P2的摄影高度为H时:
∵ ΔF1o1a1′∽ΔF1OA′;ΔF1a1a1′∽ΔF1AA′
X/x1=Z/H;Y/y1=Z/H;
∴ X=(Z/H)x1;Y=(Z/H)y1;(1)式
又∵ ΔF2a2′o2∽ΔF2A′A″;
A′A″/x2=Z/H;
则 A′A″=(Z/H)x2;
∵ A′A″=M-X;
故 M-X=(Z/H)x2;
(因x2本身带正负号,图(5)中x2为负)
M=X-(Z/H)x2;(2)式
将(1)式中X代入(2)式;
M=Z(x1-x2)/H;
得 Z=(MH)/(x1)-x2);
设 Δx=x1-x2;(3)式
∴ Z=(MH)/Δx;(4)式
从(4)式得出物体的空间高度坐标Z与立体像对中视差Δx有关,证明在X射线摄影时,像点的视差Δx是立体视觉的来源。
同时;0<Δx≥M,ΔH=H-Z,这是透射光摄影立体成像的特性。
连立方程(5)式是物体与X射线像对影像点之间空间坐标关系几何交会计算的联立方程式。
连立方程(5)式称为:透射光立体摄影坐标连立方程。
当在X射线立体摄影立体影像观察中,量测并计录到物体X射线像片P1,P2上的同名像点坐标x1,y1,x2后,用(5)式可计算得到物体内部结构密度形态的三维空间坐标。
当量测到物体不同位置,不同高度点的像点平面坐标x1,y1,x2后,可得到;
空间物体绝对高度变化为ΔH;
ΔHi=H-Zi (6)式
空间物体点与点相对高差为ΔZ;
ΔZi=Z0-Zi (7)式
空间物体点与点之间的相对视差为Δp;
Δpi=Δx0-Δxi (8)式在(6)、(7)、(8)式中;
H为X射线光源F至X光成影片P的距离。
Zi为X射线光源F至物体点的距离。
Z0为物体高度计算起始面。
Δx0为物体高度计算起始面的视差。
Δxi为不同高度的物体面视差。
Δpi为高度计算起始面的视差与不同高度面的视差之差。
五、数字立体影像计算机视觉相关匹配
X射线立体摄影获得数字影像用计算机视觉相关匹配得到影像点的像片坐标,建立物体内部立体影像数字模型(digital inside-form model,DIM)DIM。
CT-解析法的重要技术是应用计算机视觉技术代替人眼的立体影像识别与量测,完成影像和信息的自动提取及处理。由计算机自动数字影像相关匹配,找到在两张立体影像上相同的影像点,称同名点,并得到两像点的灰度和像片坐标,完成像点的空间几何交会计算,得到物体内部结构密度形态的三维空间坐标和相应的影像灰度值,建立X射线物体内部影像数字立体模型DIM。
X射线立体像对影像是光学的模拟影像,影像中已记录了物体内部结构密度形态的堆积影像灰度,经数字化后已得到影像的灰度向量阵列。
计算机视觉技术被称为影像匹配或数字相关。影像相关匹配用于建立X射线摄影物体内部影像数字立体模型DIM是关键工作,数字影像相关匹配的基本目的是:确定两张立体像片上的同名像点在各自像片上的像片坐标和灰度值,并求的同名像点之间在x方向的坐标之差Δx=x1-x2,称Δx为左右视差。影像匹配的原理是:用数值计算的方法比较探求立体像对左右像片上的影像相似程度,完成找到在左右像片上的同名影像点(相同影像点),简称同名点,并确定同名点在左右像片上的灰度值g和在像片坐标中的位置x,y。影像匹配是根据一定的准则,比较立体像对左右二张像片上影像相似性来确定是否为同名影像块,从而确定相应的同名像点,这过程是计算机对数字影像进行数值计算的方式完成影像的匹配,同名点用向量表示为;
P左=(x1,y1,g左)T;P右=(x2,y2,g右)T
数字影像匹配的方法有多种。在X射线摄影数字影像匹配中,因所取的数字影像像素大,采用一维特征匹配。一维匹配的原理是;在左像片影像中先确定一个待定点,称为目标点,以待定点为中心选取m×n(也可以m=n)个像素的灰度阵列为目标区域,称为目标窗口,为了在右像片上搜索同名影像点,在右像片上建立一个m×L(L>n)像素的灰度阵列作为搜索区。如图(6);
一维影像匹配是;搜索工作只在沿像片坐标x方向进行,这个方向也正是产生同名像点之间的相对位移的方向,既左右视差Δx产生的方向。以计算相似性测度ρi(i=i0-L/2+n/2......i0+L/2-n/2),(io)为搜索区中心像素,取c为搜索间隔设定常数。当;
ρc=max{ρi|i=i0-L/2+n/2......i0+L/2-n/2}时,(c,i0)为同名像点。
在影像匹配的过程中,确定了立体像片同名点的影像灰度值g,同时也确定了左右二张像片上同名点的像片坐标x1,y1,x2值,设定y1=y2,求Δx=x1-x2,得到影像同名点视差值Δx,利用透射光立体摄影坐标连立方程(5)式计算物体内部结构密度形态的空间坐标;
用该方程计算得到像点所对应于物体的空间位置Xi,Yi,Zi,后,可取同名点的平均灰度值g平i为该点的参考灰度,得到像点对于物体的目标向量;
pi=(Xi,Yi,Zi,g平i)T或表示为;
pi=(K1x1i,K1y1i,K2Δx,g平i)T
K1,K2是X射线摄影光源位置常数,(K1=Z/H,K2=MH。
得到物体点的目标向量pi,也就完成了数字影像的特征相关匹配工作。
X射线物体内部数字立体影像模型DIM,是目标向量p的集合为A;
A={p}={X,Y,Z,g左}={K1x1,K1y1,K2Δx,g平}
数字影像匹配建立X射线数字立体影像模型DIM的过程如图(12);
X射线物体内部数字立体影像模型DIM,是物体内部结构密度形态的目标向量的集合,是一个有准确距离、面积、体积、影像纹理的模型。用计算机技术可很容易的在DIM中实现物体内部结构密度形态3D图像的显示,实现三维物体3D图像的任意缩放、旋转、扭曲、分割、混合、融色等
六、X射线立体摄影实现Radon图像重建
X射线立体摄影影像是物体内部结构形态不同密度重叠在二维像片上的影像,数字立体影像模型DIM中是像点的重叠影像平均灰度,是靠影像的相似性相关匹配得到像点的像片灰度和坐标,既以密度变化的边界灰度确认相似性,但像点上的重叠灰度问题并未解决,对于DIM的灰度只是一个边界灰度参考值,需对DIM中的灰度进行修正,得到物体内部结构密度形态的实际空间灰度分布值。
图(13)是交线状态X射线立体摄影中的主断面;主断面是由F1F2o1o2构成的面。在这个主断面上,像点的像片坐标y=0。
在图(15)中可看出,这是一个Radon图像重建问题,这个图形与断层摄影(扫描)CT图像重建的情况相似。在多幅X射线立体摄影时已获得物体内部密度的重叠灰度g,按Radon图像重建原理;“用无限多个投影可确定二维或三维物体”的方法,能实现对X射线立体摄影图像重建,实现Radon图像重建的条件是;物体A在空间的位置、X射线束光源F的位置、影像点灰度g的位置,应在同一个坐标系
中,如图(14);
X射线立体摄影是在增强器上获得的是一个面的数字影像灰度,影像灰度点的坐标是像片坐标,物体A是空间坐标,物体A与影像灰度点g不在同一个坐标系中,如图(13),物体A与影像灰度点g没有任何关系,不满足Radon原理图像重建的条件。显然要找到X射线立体摄影与Radon原理图像重建之间的关系,建立X射线立体摄影像片坐标系与Radon原理理图像重坐标系之间的变换的方法。X射线立体摄影实现Radon原理图像重建之间的坐标变换称为;Shujia-Radon关系变换,或称SR-
坐标系变换。
在前面我们已在X射线立体摄影像片上建立了像片坐标系o-x,y,得到物体内部结构密度形态影像重叠灰度g的像片坐标g(x,y),只要能将g(x,y)换算到SR-
坐标系中的位置坐标g
中,就可实现Radon图像重建原理。要将像片二维坐标的像点灰度g(x,y)转换到三维空间SR-
坐标g
中,用透射光立体摄影坐标连立方程(5)式,能推导出X射线立体摄影实现Radon原理图像重建的Shujia-Radon关系变换方程;(9)式
式中;Δx=x1-x2
(9)式是X射线立体摄影实现Radon图像重建原理的坐标系转换方程,λ是坐标转换系数,显然求得坐标转换系数λ是实现Radon图像重建原理的第一步。首先要找到能建立起平面像片坐标系o-x,y与空间SR-
坐标系的关系条件。从图(15)看出;平面像片坐标系o-x,y与空间Radon-X,Y,Z坐标系建立的法则不同。将原Radon坐标系X,Y,Z改为SR-
坐标系,使像片坐标系o-x,y和SR-
坐标系同法则。
建立物体A与像点灰度g之间的关系;在物体A上设三个以上在SR-
坐标系中已知的点B、C、D、...N,当X射线透射物体后,如透射已知点B
C
点后,在X射线立体摄影像片P1上有已知点的影像b1和c1,在P2上有已知点的影像b2和c2如图(15)、(16);
图(15)是在平面轨道上X射线立体摄影,像片P1、P2平行。图(16)是在圆形轨道上X射线交线立体摄影,像片P1、P2不平行。
已知X射线光源F1、F2到像片P1、P2的距离同为H,已知F1到F2的距离为基线M,已知点B在像片P1上的影像点b1的像片坐标为xb1、yb1,在像片P2上的影像点b2的像片坐标为xb2、yb2,可求得Δxb=xb1-xb2。用透射光立体摄影坐标连立方程求得已知点B的立体摄影坐标点B(Xb、Yb、Zb);(5)式
因在图(15)、(16)是X射线立体摄影主断面,yb1=yb2=0,连立方程(5)式变为;(5)-1式在(5)-1式中,用已知点B在X射线立体摄影像片P1、P2上的影像b1、b2的像片坐标xb1、xb2,求得已知点B的立体摄影坐标点B(Xb、Zb),已建立起Radon原理SR-
坐标系与像片坐标系o-x,y之间的相对关系。立体摄影坐标Xb、Zb是像片坐标系o-x,y转换到Radon原理SR-
坐标系的过度坐标。
(9)式和(9)-1式建立起平面像片o-x,y坐标系与空间SR-
坐标系的关系,实现了坐标系之间的转换,将物体B、X射线光源F、灰度g统一到了同一个SR-
坐标系中。(9)-1式相似断层摄影CT在扫描时的状态。
从(9)-1式看出,当X射线光源F1、F2到像片P1、P2的距离同为H,将H钢性固定时,对F1到F2的基线距离M要求并不高,若Shujia-Radon关系变换精度达不到要求时,可反求M。这样的设定可自由的进行X射线立体摄影,降低X射线立体摄影时的精度要求。
b1和b2的像片坐标xb1,xb2由计算机视觉立体影像相关匹配求得,并计算出坐标差Δxb。进行Shujia-Radon关系变换,完全由计算机自动完成。
X射线立体摄影是得到-个面的影像,并非如图(15)、(16)时的一个主断面影像,只有在像片坐标(x,y=0)的像点,才与图(15)、(16)相符。g(x,y□0)的像点不但在x方向有射线夹角αx,在y方向也有射线夹角αy,它的计算方法比图(15)、(16)时的Shujia-Radon关系变换、Radon原理图像重建计算复杂的多。如图(17);
从以上的推导和图(15)、(16)看出,只要将立体摄影像片坐标系o-x,y的像点g(x,y)换算到SR-
坐标系位置g
中,就能确定F到g的光束路经,采集光束路经上的影像灰度值g,实现X射线立体摄影Radon原理进行图像重建。
X射线立体摄影是得到一个面的数字影像,X射线是扇形投影,用扇形投影方法,将面的数字影像按像素分成连续排列的格网灰度单元,从X射线光源F到格网灰度单元g构成光束路经,如图(18);
图(18)每个格网灰度单元对应一光束路经,多幅X射线立体摄影数字影像可组成从多方向路经灰度g,用等间隔扇形投影直线排列数据图像重建,能获得高精度的图像重建。
X射线立体摄影Shujia-Radon关系变换,实现Radon原理进行图像重建计算,目前已提出图像重建计算方法有五种;
1、连立方程法;
2、逆投影法;
3、付立叶变换法;
4、卷积法;
5、最小二乘法;
以上X射线立体摄影图像重建的步骤为;
1、自由的X射线立体摄影获得有已知点的立体数字影像。
2、用计算机视觉相关匹配的方法获得平面像片坐标系中已知点的像点坐标g(x,y)。
4、求每个格网灰度单元对应一光束路经,多幅X射线立体摄影数字影像可组成从多方向路经灰度g。
5、用等间隔扇形投影直线排列数据图像重建方法计算。
6、用Radon原理图像重建计算的灰度值对DIM灰度进行修正,得到物体内部结构密度形态的实际灰度值。
X射线立体摄影不需要精密的断层扫描机械设备,像点的灰度g(x,y)的间隔按像素分成连续排列的格网灰度单元任意设定,不需要用单个精密排列的探测器。不是等间隔的单个探测器和等间隔的线束扫描,是用一幅数字影像,在数字影像上自动按要求的像素大小无间隔的采集影像灰度,灰度g(x,y)的位置由计算机视觉相关匹配求得和像素格网灰度单元无间隔的任意设定,无影像间隔误差,提高了图像重建的精度,能做到全仿真图像重建。可快速获得X射线立体摄影在一个面的完整影像灰度,可快速的获得X射线图像重建影像。
数字影像匹配的计算方法很多,可根据匹配的要求和精度选择不同的计算方法。如德国Ackermann教授提出的最小二乘影像匹配(Least squares imagematching),影像匹配可以达到1/10甚至1/100像素的高精度,还能进行多点影像匹配和多片影像匹配,他不但解决影像重叠产生的遮蔽问题,而且对于解决图像重建时自动求得已知点坐标和多幅影像连接计算提供了方法。获得影像灰度的位置坐标,是计算机视觉自动相关匹配解算同名影像多投影光线同时交会的联立方程得到,而不是在断层扫描时,由运动机械确定影像灰度所在的探测器的位置坐标。
由诺贝尔奖获得者豪斯菲尔德发明的层析摄影CT图象重建,是完全用X射线束和电子探测器模拟Radon图像重建原理。在模拟Radon图像重建原理进行X射线束摄影时,必须保证物体A、X射线束光源F和探测器影像灰度g同在一个坐标系中,X射线束光源F和影像灰度g必须是一个个不连续的已知坐标点。当X射线束摄影透过物体A后,由F点和g点的连线构成有影像灰度的路经矢量Fg,由多方向穿过物体A的矢量Fg组成连立方程组求解,求得二维或三维物体中空间分布的不同灰度值,实现Radon原理图像重建。
在断层摄影CT设备中,是用旋转扫描方式进行断层X射线束摄影,获得X射线束穿透物体后的灰度。在圆形机架的轨道中,X射线束光源F和排列的单个有间隔的电子探测器位置是精密设定的,由电子探测器获得准直X射线束影像的灰度值g
由F
到g
构成路经射线矢量Fg,达到Radon原理图像重建的计算条件,这是用模拟Radon原理制造的仪器设备,确定X射线束光源F和电子探测器g的位置需要用精密的运动机械,精密的电子探测器和精密机械设备庞大,扫描速度与机械的运动速度有关。至今,虽然断层摄影CT图象重建设备已发展了n代,但始终没有改变模拟的设计思想,完全用机械断层X射线束扫描和电子探测器模拟Radon图像重建原理,将这种模仿Radon原理实现图像重建的方法称为;CT-模拟法。
CT-解析法是用计算机视觉技术,用解析计算的方法确定物体A、X射线光源F和灰度值g的位置。用自由摄影的X射线立体摄影数字影像,用解析计算Shujia-Radon关系变换的方法达到Radon原理图像重建条件,用计算机视觉影像匹配替代电子探测器确定影像灰度点的位置。立体影像匹配能在立体摄影的两张或多张二维平面影像中同时得到影像点的位置坐标和灰度值,它没有精密的摄影机械扫描设备和精密电子探测器,不是用X射线束获得一个点的影像灰度,而是用X射线光源辐射物体后,在X射线影像增强器上自由获得一个面的影像灰度,X射线辐射量大大降低,影像灰度获得的速度快,影像灰度分辨精度和影像点的位置精度高。将这种用解析计算实现Radon原理图像重建的方法称为;CT-解析法。
CT-解析法的过程是;
自由X射线立体摄影-立体影像相关匹配-Shujia-Radon关系变换-Radon原理图像重建。
七、实现3D图像显示、提取任意断层CT影像
在数字影像立体模型DIM中,物体内部结构密度形态分布的空间坐标X,Y,Z已获得,与之对应的影像灰度值g已获得修正,摄影时的高度H和基线M是已设定的常数。显然要得到物体内部结构密度形态的3D立体图像和断层CT影像是很容易的。
从X射线数字影像立体模型DIM提取出的断层CT影像是一个二维断面影像。计算机对平面影像的处理有很多成熟的计算处理方法,能进行断层CT影像增强、锐化、多幅CT影像拼接、特征影像提取、色彩融合等。
Claims (8)
- 本发明属计算机图像重建和数字图像处理技术领域。自由X射线立体摄影获得多幅数字影像,用X射线立体摄影空间交会求得像点空间坐标,用计算机视觉相关匹配和Shujia-Radon关系转换,使X射线光源坐标、被摄影物体坐标和数字影像的像片坐标被规划到统一的SR坐标系中,实现Radon原理进行图像重建,获得任意方向的断层CT图像。CT-解析法能同时实现三种影像观察功能;二维平片观察、全仿真立体影像观察、断层CT影像观察。CT-解析法不需要精密的X射线摄影轨道和排列的电子探测器。要求保护权利的特征如下;1、用透射光(如X射线)对物体进行立体摄影,获得物体内部密度的影像用于立体影像观察和图像重建,得到物体空间密度的分布状态。
- 2、摄影时X射线光源和影像增强器在直线轨道或圆形轨道上移动,连续的进行X射线立体摄影获得多幅数字平面影像。
- 3、在物体方设置多个固定已知点,X射线摄影后在影像中有多个固定已知点的数字影像。用透射光立体摄影坐标连立方程计算影像点坐标。
- 4、在平面数字影像上建立像片坐标系,用计算机影像相关匹配找到已知点的数字影像灰度和坐标。将数字影像在像片坐标系分成连续排列的像素格网灰度单元。
- 5、用立体摄影坐标连立方程求出已知点坐标和坐标转换系数。将像素格网灰度单元的像片坐标和X射线光源坐标用关系转换规划到设定的已知点所在的图像重建坐标系中。
- 6、由X射线光源与像素格网灰度单元构成方向射线路经矢量,多方向射线路经矢量组成连立方程求解,得到物体空间密度的分布。
- 7、实现三种影像观察功能;二维平片观察、全仿真立体影像观察、断层CT图像观察。
- 8、在具有中心投影性质的光、声、波发射源构成影像,进行图像重建得到物体空间密度的分布中使用1、2、3、4、5、6、7。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 02133358 CN1464301A (zh) | 2002-06-25 | 2002-06-25 | X射线立体摄影实现图像重建的ct-解析法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 02133358 CN1464301A (zh) | 2002-06-25 | 2002-06-25 | X射线立体摄影实现图像重建的ct-解析法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN1464301A true CN1464301A (zh) | 2003-12-31 |
Family
ID=29744109
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 02133358 Pending CN1464301A (zh) | 2002-06-25 | 2002-06-25 | X射线立体摄影实现图像重建的ct-解析法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1464301A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009043224A1 (fr) * | 2007-09-24 | 2009-04-09 | Shu Jia | Procédé de reconstruction d'image utilisant une photographie volumique par rayons x |
CN104244833A (zh) * | 2012-03-26 | 2014-12-24 | 皇家飞利浦有限公司 | X射线焦斑移动的方向控制 |
CN114199905A (zh) * | 2021-12-13 | 2022-03-18 | 中国航发南方工业有限公司 | 一种机匣内部缺陷的空间定位方法及系统 |
-
2002
- 2002-06-25 CN CN 02133358 patent/CN1464301A/zh active Pending
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009043224A1 (fr) * | 2007-09-24 | 2009-04-09 | Shu Jia | Procédé de reconstruction d'image utilisant une photographie volumique par rayons x |
US20100195891A1 (en) * | 2007-09-24 | 2010-08-05 | Jia Shu | Image reconstructing methjod using x-ray volume photography |
CN101126725B (zh) * | 2007-09-24 | 2010-12-15 | 舒嘉 | 采用x射线容积摄影实现图像重建的方法 |
JP2010540893A (ja) * | 2007-09-24 | 2010-12-24 | ジャ シュ | X線容積撮影による画像再構成方法 |
RU2469298C2 (ru) * | 2007-09-24 | 2012-12-10 | Цзя ШУ | Способ восстановления изображения с использованием объемной рентгеновской фотосъемки |
US8687861B2 (en) | 2007-09-24 | 2014-04-01 | Jia Shu | Image reconstructing method using X-ray volume photography |
CN104244833A (zh) * | 2012-03-26 | 2014-12-24 | 皇家飞利浦有限公司 | X射线焦斑移动的方向控制 |
US9615802B2 (en) | 2012-03-26 | 2017-04-11 | Koninklijke Philips N.V. | Direct control of X-ray focal spot movement |
CN104244833B (zh) * | 2012-03-26 | 2017-09-01 | 皇家飞利浦有限公司 | X射线焦斑移动的方向控制 |
CN114199905A (zh) * | 2021-12-13 | 2022-03-18 | 中国航发南方工业有限公司 | 一种机匣内部缺陷的空间定位方法及系统 |
CN114199905B (zh) * | 2021-12-13 | 2024-02-20 | 中国航发南方工业有限公司 | 一种机匣内部缺陷的空间定位方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111060023B (zh) | 一种高精度3d信息采集的设备及方法 | |
CN100346748C (zh) | 图像处理设备和图像处理方法 | |
Farahani et al. | Three-dimensional imaging and scanning: current and future applications for pathology | |
CN102144927B (zh) | 基于运动补偿的ct设备和方法 | |
WO2018018981A1 (zh) | 用于对三维表面的图像进行重建的方法、装置和系统 | |
US20220301195A1 (en) | Methods and systems for imaging a scene, such as a medical scene, and tracking objects within the scene | |
CN207516234U (zh) | 一种光学投影层析成像的装置 | |
CN107155341B (zh) | 三维扫描系统和框架 | |
CN110288642A (zh) | 基于相机阵列的三维物体快速重建方法 | |
CN111076674B (zh) | 一种近距离目标物3d采集设备 | |
CN106572829A (zh) | 医学成像系统的扫描区域的定位的控制 | |
CN109934935B (zh) | 一种服装与人体图像穿帮、匹配迁移的方法与系统 | |
CN103854270A (zh) | 一种ct和mr异机三维图像融合的配准方法及系统 | |
CN1665313A (zh) | 产生图像的装置和方法 | |
EP3382645B1 (en) | Method for generation of a 3d model based on structure from motion and photometric stereo of 2d sparse images | |
CN107221029A (zh) | 一种三维图像重建方法 | |
CN102779354B (zh) | 基于光度立体技术的中医望诊信息表面三维重建方法 | |
AU2018301580B2 (en) | Three-dimensional ultrasound image display method | |
CN111060006A (zh) | 一种基于三维模型的视点规划方法 | |
CN112629412A (zh) | 一种转动式3d智能视觉设备 | |
CN107256575A (zh) | 一种基于双目立体视觉的三维舌像重建方法 | |
CN1600273A (zh) | 断层造影成像检查设备的运行方法和x射线断层造影设备 | |
CN110619601B (zh) | 一种基于三维模型的图像数据集生成方法 | |
CN106934110B (zh) | 一种由聚焦堆栈重建光场的反投影方法和装置 | |
CN1464301A (zh) | X射线立体摄影实现图像重建的ct-解析法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |