CN101126725B - 采用x射线容积摄影实现图像重建的方法 - Google Patents

采用x射线容积摄影实现图像重建的方法 Download PDF

Info

Publication number
CN101126725B
CN101126725B CN2007100500943A CN200710050094A CN101126725B CN 101126725 B CN101126725 B CN 101126725B CN 2007100500943 A CN2007100500943 A CN 2007100500943A CN 200710050094 A CN200710050094 A CN 200710050094A CN 101126725 B CN101126725 B CN 101126725B
Authority
CN
China
Prior art keywords
ray
image data
coordinate system
battle array
face battle
Prior art date
Application number
CN2007100500943A
Other languages
English (en)
Other versions
CN101126725A (zh
Inventor
舒嘉
Original Assignee
舒嘉
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 舒嘉 filed Critical 舒嘉
Priority to CN2007100500943A priority Critical patent/CN101126725B/zh
Publication of CN101126725A publication Critical patent/CN101126725A/zh
Application granted granted Critical
Publication of CN101126725B publication Critical patent/CN101126725B/zh

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/006Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods

Abstract

本发明公开了采用X射线容积摄影实现立体图像重建的方法,该方法根据数字摄影测量的定位原理,对被关注对象拍摄多张X射线的面阵影像数据,在拍摄X射线像片时,将已知三维位置的多个标志点同时拍摄在X射线像片中,利用标志点在容积摄影立体图像重建直角坐标系中的坐标和标志点在X射线像片坐标系中的坐标建立所述两坐标系的转换关系,通过这个转换关系来确定每条X摄影光束的路径和位置,得到每条X光束通过被关注对象的全部体素,由同一条X光束路径上的体素g和像元G列出Radon方程组求解体素g的灰度值,通过计算机解算出全部体素g的灰度值,由体素g的灰度值重建被关注对象的三维立体容积图像。

Description

采用χ射线容积摄影实现图像重建的方法

技术领域

[0001] 本发明属于图像重建技术领域,特别涉及一种物体内部状态的图像重建方法。 背景技术

[0002] 关于物体内部状态的图像重建,现今主要是以豪斯非尔德(Hoimsfield)发明的 方法进行CT (Computend Temography)扫描,即在精密机械扫描机架上安置X射线源1和 探测器3 (如图1所示),使X射线源1和探测器3围绕物体(被关注对象)4旋转扫描,在 旋转的同时精确的确定每束X光线的发出位置和接收位置,再利用求物体内部密度分布的 Radon计算方法,实现图像重建(参考书:计算机体层成像,美Willi A. Kalender,人民卫生 出版社ISBN7-117-05284-8/R. 5285。实用CT检查技术学,王鸣鹏,科学技术文献出版社, ISBN7-5023-3262-6/R. 656)。

[0003] 自1972年第一台CT机问世以来,CT机已经过了八代的发展(“代”是其发展进 度标志,实质是如何提高扫描速度),但扫描定位原理没有任何变化,一直是采用精密机械 机架滑环或光束钟形器扫描断层数据的模式,这种模拟Radon计算条件的定位、扫描方式 限制了 CT设备的构成和设备部件的充分利用,由此导致:扫描机械和电子设备制造复杂,X 射线照射量大、利用率低,数据采集速度慢、扫描数据精度低,重建图像分辨率低,而且没有 成立真正意义上的三维图像。

发明内容

[0004] 本发明的目的在于克服现有扫描技术的不足,提供一种采用X射线容积摄影实现 图像的重建的方法,此种方法以获得面状X射线摄影影像重建被关注对象的三维立体图 像,对三维立体图像分割任意方向的断层图像,不仅影像分辨率高,X射线照射量少,而且设 备结构简单,制造难度低,投资少。

[0005] 本发明技术方案:根据数字摄影测量的定位原理,对被关注对象拍摄多张X射线 面阵影像数据,在拍摄X射线像片时,将已知三维位置的多个标志点同时拍摄在X射线像片 中,利用标志点在容积摄影立体图像重建直角坐标系O-XYZ中的坐标和标志点在X射线像 片坐标系o-xy中的坐标建立所述两坐标系的转换关系,通过这个转换关系来确定位每条X 射线摄影光束的路径和位置,得到每条X光束通过被关注对象的全部体素,由同一条X光束 路径上的体素g和像元G列出Radon方程组求解体素g的灰度值,通过计算机解算出全部 体素g的灰度值,由体素g的灰度值重建被关注对象的三维立体容积图像。

[0006] 本发明所述采用X射线容积摄影实现图像重建的方法,步骤依次如下:

[0007] (1)鉴定X射线摄影机

[0008] 调整X射线摄影机,使X射线管中心S的垂直光束与X射线数字像片承影面垂直, 然后测定X射线管中心S到X射线数字像片承影面的距离H。

[0009] (2)建立像片坐标系o-xy

[0010] 在X射线数字像片上建立像片坐标系o-xy,像片坐标系原点ο是X射线管中心S的垂直光束在X射线数字像片上的像元点,过像片坐标系原点ο的像元点横向排列方向为 像片坐标系χ轴,过像片坐标系原点ο的像元点纵向排列方向为像片坐标系y轴。

[0011] (3)建立容积摄影立体图像重建直角坐标系O-XYZ

[0012] 对被关注对象至少设置三个标志点并事先测定各标志点之间的三维距离,以任一 标志点作为坐标系原点0,建立容积摄影立体图像重建直角坐标系0-ΧΥΖ,将各标志点在坐 标系O-XYZ中的坐标存入计算机;

[0013] 所述标志点设置在被关注对象上或标志架上,标志架围绕被关注对象设置。

[0014] (4) X射线摄影

[0015] 操作X射线摄影机,使X射线管和X射线数字像片承影面沿弧线或直线轨道运动, 在相隔一定间距的各位点拍摄被关注对象的X射线数字像片,得被关注对象不同位置的多 张X射线数字像片。

[0016] (5)标志点模式识别与标志点影像定位

[0017] 将步骤(4)拍摄的所有X射线数字像片输入计算机,通过计算机对各张X射线数 字像片进行标志点的模式识别与标志点影像定位,得到各标志点在各张像片o-xy坐标系 中的坐标X、y (“标志点的模式识别与标志点影像定位”见张祖勋张剑清“数字摄影测量学”, 武汉大学出版社2001. 5ISBN 7-307-03233-3/P. 7第四章至第六章);

[0018] 所述标志点模式识别是通过做好的标志点形状模板影像和特征提取影像匹配 来完成的,特征提取用各种算子来完成;由于标志点的几何形状各不相同,这些不同形状 的标志点在X像片上具有不同的影像特征,因此特征提取算子可分为点特征提取算子、 线特征提取算子和面特征提取算子;点特征提取算子主要有Moravec算子、Harmah算子 与Forstner算子,线特征提取算子主要有剃度算子、二阶差分算子、高斯一拉普拉斯算子 (LOG算子)、特征分割法、Hough变换,面特征提取主要通过影像区域分割来完成,影像分割 法主要有阀值法、区域生长法、集群分类法。

[0019] 所述标志点影像定位是指精确的确定标志点定位目标处在X射线数字像片坐标 系o-xy中的位置,在提取到的标志点形状区域内,找到标志点的定位目标处,定位目标处 的模式识别计算方法和前述特征提取算子相同。

[0020] (6)计算各张X射线数字像片的方位元素\、Ys、Zs、α、ω、κ

[0021] 用计算机列出各张X射线数字像片上的全部标志点的X射线光束的共线方程式, 对任一标志点(例如标志点Α),其X射线光束共线方程式如下:

x = _H ay{X Α ~ X s)+ bx{YA-Ys)+ cx{Z Α - Zs) α3(Χ Α - X s) + b3(YA - Vs) + c3(Z Α - Zs)

[0022] (1-1) v 二 _H a2{X A- X s)^b2{YA-Ys)+ C2(Za-Zs)

a3{X A - xs)+ b3{YA -Ys)+ A - Z s) J

[0023] 上式中,x、y为标志点A在像片坐标系o-xy中的坐标,H为X射线管中心S到X射 线数字像片承影面的距离,Xa、Ya、Za为标志点A在容积摄影立体图像重建直角坐标系O-XYZ 中的坐标,Xs> Ys、Zs为X射线管中心S在容积摄影立体图像重建直角坐标系O-XYZ中的坐

标,旋转因子为:

[0024]/ a\ a2 Gt3、 f b2 、 =KcI C2 C3, V

cosa 0 - sinaYl 0 0 Ycos/f — sinx* 0)

0

sina 0

0

cosa 人1

sin/f cosa: 0 0 0 1

(1-2)

0 cosa> -sima 0 sino? cos6)

[0025] 在各张X射线数字像片上至少选择三个标志点的X射线光束共线方程式联立解 算,即得出各张X射线数字像片的六个方位元素XS、YS、ZS、α、ω、κ。方位元素中,XS、YS、 Zs为X射线管中心S在容积摄影立体图像重建直角坐标系O-XYZ中的坐标,α、ω、κ为 X射线数字像片在容积摄影立体图像重建直角坐标系O-XYZ中的旋转方向夹角(“旋转因 子”和“共线方程式”原理和算法,见钱曾波“解析摄影测量”测绘出版社15039.新143,第 一章)。

[0026] (7)计算像元G在容积摄影立体图像重建直角坐标系O-XYZ中的坐标

[0027] 依据步骤(6)得到的各张X射线数字像片的方位元素XS、YS、ZS、α、ω、κ禾PX射 线管中心S到X射线数字像片承影面的修正距离(修正距离通过X射线数字像片上 的标志点的X射线光束共线方程式1-1计算)、像元之间的距离1,用计算机计算各张X射线 数字像片上每个像元Gij的容积摄影立体图像重建坐标Xu、Yij. Zij,计算原理和方法如下:

[0028] X射线数字像片上的像元点是按等距离有序排列的,像元G之间的距离1是固定 的。设i、j为像元序号,设定像片坐标系o-xy的坐标轴x、y按像元有序排列的横、纵方向 平行,设定每个像元Gij的像片坐标χ、y是以投影中心S到像片P的垂直距离的交点ο处 起算,S到ο的垂直距离为H。在交点ο处位置的像元点Gtl的X像片坐标为χ = 0、y = 0。 由步骤(6)可知,X射线管中心S在容积摄影立体图像重建直角坐标系O-XYZ中的坐标为: S(Xs、Ys、Zs),引入已计算出的旋转方向夹角元素α、ω、κ,那么就有方程:

[0029] 像元Gtl在容积摄影立体图像重建直角坐标系O-XYZ中的坐标

JTgo =Xs + (-H cosiy sin α) [οο3ο] Ygo =YsjT (-H sin ω)

Zgq = Zs + (-H cos ω cos a)

[0031 ] 像元Gij在容积摄影立体图像重建直角坐标系O-XYZ中的坐标

XGij =Xs+ (-H cosiy sin or) + (j x 0 cos λ: [oo32] YGij =Ys+ {-Hsm ω) + (i χ /) cos κ ZGij =Zs +(-//cosa?cosa)

上式中,i、j为像元序号,1为像元之间的距离,iXl =Xi, iXl =y”将(1_4)

>

(1-3)

>

(1-4)

[0033]

式表示为:

XGij =Xs+ (-H cos ω sin ά) + Xi cos κ [οο34] YGij =Ys ^ i~H sin ω) + yt cos κ ZGij =Zs+ (-H cos ω cos a)

[0035] 用(1-5)式计算得到的像元Gu在容积摄影立体图像重建直角坐标系O-XYZ中的 的坐标为:

r (1-5)[0036] Gij (Xij, Yij, Zij)

[0037] 可以看出,(1-3)式和(1-4)式及(1-5)式表示的是一个将X射线数字像片上的 像元G坐标(x、y)在容积摄影立体图像重建直角坐标系中重排列的过程。

[0038] 通过(1)至(7)步骤,已求得X射线发射中心S和X射线像片像元G在容积摄影 立体图像重建直角坐标系中的坐标,显然已确定了 X射线发射中心S和X射线数字像片上 全部像元的位置,这过程与传统的CT精密机械机架扫描功能相同,只是传统的CT扫描是用 机械的方式获得X射线发射中心S和X射线接收器的位置,本方法是用计算的方式获得X 射线发射中心S和X射线数字像片上全部像元的位置。

[0039] (8)确定每条光束通过被关注对象的全部体素

[0040] 依据步骤(6)得到的各张X射线数字像片的方位元素XS、YS、ZS和步骤(7)得到的 各张X射线数字像片的各像元G的容积摄影立体图像重建坐标Xu、Yij, Zij,通过计算机计 算,找出每条X光束通过被关注对象的全部体素g,计算原理和方法如下:

[0041] 由X射线管中心S和像元Gij两点构成的光束过重建体素gijk方格构成的直线方 程为:

X—Xs Y-Ys Z, — Zs

[0042] _Y=V-V=7-7 ( 1_6)

入 GiJ-入S 1 Glj —1S ^Gjj

[0043] 因此,判定X光束通过体素guk的条件是满足直线方程(1-6)式,找出满足直线方 程(1-6)式上的体素gijk位置坐标,就是光束经过的路径。

[0044] 体素gijk方格的设定大小是:gijk = ItXtXtI,体素gijk方格的八个顶点在容积 (VP)摄影立体图像重建直角坐标系中已知的坐标(见图7)是:

[0045]

g ijk 1 ( X ijk 1 > Y ijk 1 ^ ^ijkO

gijk2 (Xj/fc2、YyA;2^ Z yk

g ijk 8 (XiJkS、Y/A/'8> Z ijk 8 )

[0046] 由(1-6)式计算得到gijk体素坐标值限定在(gijkl.....giJk8)范围内,即认为该X

光束通过该体素。找出每条X光束SGu通过的全部体素guk,就确定了每条光束路径。

[0047] 由于体素gijk方格的八个顶点在容积摄影立体图像重建直角坐标系中的坐标已 知,例如gukl的坐标为(Xi、h、zk、),因此光束路径计算步骤如下:

[0048] ①将Z = Zk代入式(1-6)所得到的关于X、Y的方程就是该X光线在高程Z = Zk 面内的投影方程。

[0049] ②判断该投影方程与体素方格底面四边有无相交,如果没有相交,则该X光线不 经过该体素。如果相交,则求出其交点坐标,设交点坐标为(Spb1)和(a2、b2),需要进一步 判断。

[0050] ③将第二步的交点坐标分别代入式(1-6)求出Z值,判断Z值是否在Z1^PZfta 为正方体体素方格的边长)之间,如果在,则说明该光线经过该体素方格。

[0051] (9)重建三维立体容积图像

[0052] 由每一条X光束路径上的像元G灰度值和体素g列出Radon方程组,通过计算机解算出全部体素g的灰度值(体素灰度值的解算用Radon方法),由体素的灰度值重建被关 注对象的三维立体容积图像(参考书:见数字图像处理导论第七章,赵荣春,西北工业大学 出版社,ISBN7-5612-0737-9/TB. 82)。

[0053] 本发明所述方法有多种用途:

[0054] 1、二维平片观察

[0055] 本发明所述方法步骤(4)X射线摄影,所获得的是被关注对象不同位置的多张X射 线数字像片,通过X射线数字像片,可对被关注对象进行二维平片观察。

[0056] 2、二维数字减影

[0057] X射线增强摄影后,本发明所述方法步骤(4)所得到的X射线数字像片可直接进行 影像灰度反向,与X像片叠加,形成二维数字减影影像。

[0058] 3、实时立体影像观察

[0059] 利用立体观察原理,用两张在不同位置拍摄的有重叠的X射线数字像片,能建立 体像对,对被关注对象进行实时立体观察。

[0060] 4、三维特征影像提取

[0061] 在体素gijk的灰度值重建的三维立体图像中,可提取出特征相同部位的灰度影像 坐标,构成特征影像图观察。

[0062] 5、实时立体定位导航

[0063] 利用提取出的特征相同部位灰度的影像坐标确定被关注对象的空间位置,引导器 械沿确定的路线坐标到达确定的空间位置。

[0064] 6、提取任意断层图像

[0065] 在体素gijk的灰度值重建的三维立体图像中,可提取被关注对象的任意断层图像。

[0066] 本发明具有以下有益效果:

[0067] 1、本发明所述方法用普通的数字X射线机摄影,得到所拍摄物体的二维面状影像 数据,以面状摄影影像重建立体图像,得到整个物体的影像数据,为物体内部状态图像重建 提供了 一种不同构思的技术方案。

[0068] 2、本发明所述方法通过计算机计算确定X光束的路径和位置,与CT扫描用精密机 械滑环确定X光束路径和位置相比,不仅影像分辨率高,X射线照射量少,而且所用设备结 构简单,投资少,耗电少,对工作环境要求低。

[0069] 3、本发明所述方法不仅能重建被关注对象的三维立体图像,而且还具有二维平片 观察、二维数字减影、实时立体影像观察、实时立体定位导航、三维特征影像提取、任意断层 图像提取等用途。

[0070] 4、本发明所述方法可广泛用于被关注对象内部探伤、安全检查、海关检查、被关注 对象内部空间定位导航、医学检测等领域。

附图说明

[0071] 图1是计算机体层成像方法扫描(CT)的原理图。

[0072] 图2是本发明所述采用X射线容积摄影实现图像重建方法的流程图。

[0073] 图3是本发明所述采用X射线容积摄影实现图像重建方法的用途说明图。

[0074] 图4是像片坐标系o-xy的示意图。[0075] 图5是X射线容积摄影时X射线源、被关注对象、标志架、X射线数字像片的位置 关系图,该图还描述了标志点和容积摄影立体图像重建直角坐标系O-XYZ。

[0076] 图6是像片坐标系o-xy与容积摄影立体图像重建直角坐标系O-XYZ的关系图。

[0077] 图7是体素gijk方格与容积摄影立体图像重建直角坐标系O-XYZ的关系图。

[0078] 图8、图9、图10、图11是是多层网架4个不同位置的X射线数字像片,其中,图8 为第1张X射线数字像片,图9为第2张X射线数字像片,图10为第3张X射线数字像片, 图11为第4张X射线数字像片。

[0079] 图中,1一X射线管中心(X射线源)、2—扫描机架、3—探测器、4一被关注对象、5— 图像重建方格、6— X射线数字像片、7—标志架、8—标志点、9一承物床、10—标志点像元、 11 一体素方格。

具体实施方式

[0080] 下面结合附图对本发明所述方法作进一步说明。在下述实施例中,通过X射线容 积摄影对多层网架实现图像重建,步骤如下:

[0081] (1)鉴定X射线摄影机、

[0082] 调整X射线摄影机,使X射线管中心1 (S)的垂直光束与X射线数字像片承影面垂 直,确定X射线管中心1⑶到X射线数字像片承影面的距离H为1000mm(见图4)。

[0083] (2)建立像片坐标系o-xy

[0084] 在X射线数字像片6上建立像片坐标系o-xy,像片坐标系原点ο是X射线管中心 I(S)的垂直光束在X射线数字像片上的像元点,过像片坐标系原点O的像元点横向排列方 向为像片坐标系χ轴,过像片坐标系原点ο的像元点纵向排列方向为像片坐标系y轴,如图 4所示。

[0085] (3)建立容积摄影立体图像重建直角坐标系O-XYZ

[0086] 将用0. 5mm直径铁丝绕制的50 X 50 X 50mm的多层网架设置平面标志架,在标志架 上设置8个标志点A、B、C、D、E、F、G、H并测定各标志点之间的三维距离,以标志点A作为 坐标系原点,建立容积摄影立体图像重建直角坐标系O-XYZ。各标志点在坐标系O-XYZ中的 坐标见表1。

[0087] 表1各标志点在坐标系O-XYZ中的坐标[0088]

[0089] 将各标志点在坐标系O-XYZ中的坐标存入计算机,所用计算机为一般的PC机。

[0090] (4) X射线摄影

[0091] 操作X射线摄影机,使X射线管和X射线数字像片承影面沿直线轨道运动,在相隔 50mm间距的4个位点各拍摄一张X射线数字像片,得多层网架不同位置的4张X射线数字 像片(见图8、图9、图10、图11)。

[0092] (5)标志点模式识别与标志点影像定位

[0093] 将步骤(4)拍摄的4张X射线数字像片输入步骤⑶所述的计算机,通过计算机 用点特征提取算子Moravec匹配计算提取各标志点特征影像,确定各标志点在各张X射线 数字像片上的像元位置,得到各标志点在各张像片o-xy坐标系中的坐标。

[0094] (6)计算各张父射线数字像片的方位元素^、&、0、《、κ

[0095] 将各张X射线数字像片上的标志点在容积摄影立体图像重建直角坐标系O-XYZ中 的坐标和步骤(5)得到的标志点在o-xy坐标系中的坐标分别代入X射线光束共线方程式, 对各张X射线数字像片选择三个标志点的X射线光束共线方程式,通过计算机联立解算,算 出各张X射线数字像片的六个方位元素(计算方法见钱曾波“解析摄影测量”,测绘出版社 15039.新143,第一章),各张X射线数字像片的方位元素见表2。

[0096] 表2方位元素计算结果表

[0097]

[0098] (7)计算像元G在容积摄影立体图像重建直角坐标系O-XYZ中的坐标

[0099] 依据步骤(6)得到的4张X射线数字像片的方位元素XS、YS、ZS、α、ω、κ和X射 线管中心I(S)到X射线数字像片承影面的修正距离H·®*= 1000mm),像元之间的距离 1(1 = 100 μ m),用计算机计算各张X射线数字像片上每个像元Gij的容积摄影立体图像重 建坐标Xu、Yij. Zijo各张X射线数字像片上的标志点像元Gij的容积摄影立体图像重建坐 标 Xij、Yi^Zij 见表 3。

[0100] 表3各标志点像元Gij在坐标系O-XYZ中的坐标

[0102] (8)确定每条光束通过多层网架的全部体素[0103] 依据步骤(6)得到的各张X射线数字像片的方位元素XS、YS、ZS和步骤(7)得到的 各张X射线数字像片的各像元G的容积摄影立体图像重建坐标Xu、Yij, Zij,通过计算机计 算,找出每条X光束通过多层网架的全部体素g。

[0104] (9)重建三维立体图像

[0105] 由每一条X光束路径上的像元G和体素g列出Radon方程组,通过计算机解算出 全部体素的灰度值,由体素的灰度值重建多层网架的立体容积图像。

[0106] 在上述立体影像观察中测量标志点和多层网架特征部位坐标,提取体素g的灰度 影像坐标与实测标志点和多层网架特征部位坐标比较,完全相符,证明标志点模式识别与 标志点影像定位影像匹配正确,相交的光束构成的前方交会计算的同名点对对相交,共线 条件计算正确,X射线投影中心S和像元G的像片坐标构成的光束路径是符合Radon计算 条件的。

Claims (3)

  1. 采用X射线容积摄影实现图像重建的方法,其特征在于根据数字摄影测量的定位原理,对被关注对象拍摄多张X射线面阵影像数据像片,在拍摄X射线面阵影像数据像片时,将已知三维位置的至少三个标志点同时拍摄在X射线面阵影像数据像片中,利用标志点在容积摄影立体图像重建直角坐标系中的坐标和标志点在像片坐标系中的坐标建立所述两坐标系的转换关系,通过这个转换关系来确定每条X射线摄影光束的路径和位置,得到每条X射线摄影光束通过被关注对象的全部体素,由同一条X射线摄影光束路径上的体素和像元列出Radon方程组求解体素的灰度值,通过计算机解算出全部体素的灰度值,由体素的灰度值重建被关注对象的三维立体容积图像。
  2. 2.根据权利要求1所述的采用X射线容积摄影实现图像重建的方法,其特征在于步骤 依次如下:(1)鉴定X射线摄影机调整X射线摄影机,使X射线管中心S(I)的垂直光束与X射线面阵影像数据像片(6) 承影面垂直,然后测定X射线管中心S到X射线面阵影像数据像片承影面的距离H ;(2)建立像片坐标系o-xy在X射线面阵影像数据像片(6)上建立像片坐标系o-xy,像片坐标系原点ο是X射线 管中心S(I)的垂直光束在X射线面阵影像数据像片上的像元点,过像片坐标系原点ο的像 元点横向排列方向为像片坐标系χ轴,过像片坐标系原点ο的像元点纵向排列方向为像片 坐标系y轴;(3)建立容积摄影立体图像重建直角坐标系O-XYZ对被关注对象(4)至少设置三个标志点并测定各标志点(8)之间的三维距离,以任一 标志点作为坐标系原点0,建立容积摄影立体图像重建直角坐标系0-ΧΥΖ,将各标志点在容 积摄影立体图像重建直角坐标系O-XYZ中的坐标存入计算机;(4) X射线摄影操作X射线摄影机,使X射线管和X射线面阵影像数据像片承影面沿弧线或直线轨道 运动,在相隔一定间距的各位点拍摄被关注对象(4)的X射线面阵影像数据像片,得被关注 对象不同位置的多张X射线面阵影像数据像片;(5)标志点模式识别与标志点影像定位将步骤(4)拍摄的所有X射线面阵影像数据像片输入计算机,通过计算机对各张X射 线面阵影像数据像片进行标志点的模式识别与影像定位,得到各标志点在各张像片的像片 坐标系o-xy中的坐标χ、y ;(6)计算各张X射线面阵影像数据像片的方位元素XS、YS、ZS、α、ω、κ用计算机列出各张X射线面阵影像数据像片上的全部标志点的X射线摄影光束的共线 方程式,在各张X射线面阵影像数据像片上至少选择三个标志点的X射线摄影光束共线方 程式联立解算,得到各张X射线面阵影像数据像片的六个方位元素XS、YS、ZS、α、ω、K,所 述方位元素中,Xs> Ys> Zs为X射线管中心S在容积摄影立体图像重建直角坐标系O-XYZ中 的坐标,α、ω、Κ为X射线面阵影像数据像片在容积摄影立体图像重建直角坐标系O-XYZ 中的旋转方向夹角;(7)计算像元在容积摄影立体图像重建直角坐标系O-XYZ中的坐标依据步骤(6)得到的各张X射线面阵影像数据像片的方位元素XS、YS、ZS、α、ω、κ和X射线管中心S到X射线面阵影像数据像片承影面的修正距离、像元之间的距离1,用计 算机计算各张X射线面阵影像数据像片上每个像元的容积摄影立体图像重建直角坐标Xu、 Yij^ Zij,所述修正距离通过X射线面阵影像数据像片上的标志点的X射线摄影光束共线 方程式计算;(8)确定每条X射线摄影光束通过被关注对象的全部体素依据步骤(6)得到的各张X射线面阵影像数据像片的方位元素XS、YS、ZS和步骤(7)得 到的各张X射线面阵影像数据像片的各像元的容积摄影立体图像重建直角坐标Χ『Υ。_、Ζ。_, 通过计算机计算,找出每条X射线摄影光束通过被关注对象的全部体素;(9)重建三维立体容积图像由每一条X射线摄影光束路径上的像元灰度值和体素列出Radon方程组,通过计算机 解算出全部体素的灰度值,由体素的灰度值重建被关注对象的三维立体容积图像。
  3. 3.根据权利要求2所述的采用X射线容积摄影实现图像重建的方法,其特征在于标志 点设置在被关注对象上或标志架上,所述标志架围绕被关注对象设置。
CN2007100500943A 2007-09-24 2007-09-24 采用x射线容积摄影实现图像重建的方法 CN101126725B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2007100500943A CN101126725B (zh) 2007-09-24 2007-09-24 采用x射线容积摄影实现图像重建的方法

Applications Claiming Priority (6)

Application Number Priority Date Filing Date Title
CN2007100500943A CN101126725B (zh) 2007-09-24 2007-09-24 采用x射线容积摄影实现图像重建的方法
US12/679,680 US8687861B2 (en) 2007-09-24 2008-08-04 Image reconstructing method using X-ray volume photography
DE200811002547 DE112008002547T5 (de) 2007-09-24 2008-08-04 Bildrekonstruktionsverfahren mittels Röntgenvolumenphotographie
RU2010144360/28A RU2469298C2 (ru) 2007-09-24 2008-08-04 Способ восстановления изображения с использованием объемной рентгеновской фотосъемки
PCT/CN2008/001417 WO2009043224A1 (en) 2007-09-24 2008-08-04 Image reconstructing method using x-ray volume photography
JP2010525180A JP2010540893A (ja) 2007-09-24 2008-08-04 X線容積撮影による画像再構成方法

Publications (2)

Publication Number Publication Date
CN101126725A CN101126725A (zh) 2008-02-20
CN101126725B true CN101126725B (zh) 2010-12-15

Family

ID=39094803

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2007100500943A CN101126725B (zh) 2007-09-24 2007-09-24 采用x射线容积摄影实现图像重建的方法

Country Status (6)

Country Link
US (1) US8687861B2 (zh)
JP (1) JP2010540893A (zh)
CN (1) CN101126725B (zh)
DE (1) DE112008002547T5 (zh)
RU (1) RU2469298C2 (zh)
WO (1) WO2009043224A1 (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101126725B (zh) 2007-09-24 2010-12-15 舒嘉 采用x射线容积摄影实现图像重建的方法
US9420262B1 (en) 2013-03-15 2016-08-16 Matthew Myers Gyroscopic two dimensions of freedom photography system and method
CN104730091B (zh) * 2015-02-10 2018-01-16 西安交通大学 基于区域分割探测的燃气轮机叶片缺陷提取与分析方法
CN106959308B (zh) * 2017-03-20 2020-04-07 东南大学 一种混凝土结构火灾影响深度检测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1376947A (zh) * 2001-03-26 2002-10-30 舒嘉 X射线立体图像观察系统
CN1464301A (zh) * 2002-06-25 2003-12-31 舒嘉 X射线立体摄影实现图像重建的ct-解析法
DE10244180A1 (de) * 2002-09-23 2004-04-08 Siemens Ag Verfahren zur Bilderstellung in der Computertomographie eines periodisch bewegten Untersuchungsobjektes und CT-Gerät zur Durchführung des Verfahrens
FR2862790A1 (fr) * 2003-11-25 2005-05-27 Ge Medical Syst Sa Procede et dispositif pour l'observation d'un phenomene transitoire a l'aide d'un arceau vasculaire
FR2864669A1 (fr) * 2003-12-24 2005-07-01 Ge Med Sys Global Tech Co Llc Procede pour ameliorer la visualisation d'un vaisseau

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4630203A (en) * 1983-12-27 1986-12-16 Thomas Szirtes Contour radiography: a system for determining 3-dimensional contours of an object from its 2-dimensional images
JPS62198965A (en) * 1986-02-27 1987-09-02 Shigeo Ishii 3-dimensional image processing system
SU1581281A1 (ru) * 1988-02-02 1990-07-30 А.П.Золотухин Способ формировани изображений при компьютерной томографии головы
US6289235B1 (en) 1998-03-05 2001-09-11 Wake Forest University Method and system for creating three-dimensional images using tomosynthetic computed tomography
RU2171630C2 (ru) * 1999-06-18 2001-08-10 Пестряков Андрей Витальевич Способ совмещения трехмерных изображений, полученных с помощью компьютерных томографов, работающих на основе различных физических принципов
DE10141406B4 (de) * 2001-08-23 2008-01-17 Siemens Ag Verfahren und Vorrichtung zum Detektieren der dreidimensionalen Position eines in einem Körperbereich eingeführten Untersuchungsinstruments
DE10201644A1 (de) * 2002-01-17 2003-08-07 Siemens Ag Registrierungsprozedur bei projektiver intraoperativer 3D-Bildgebung
JP3698130B2 (ja) * 2002-07-22 2005-09-21 株式会社日立製作所 Ct装置,ct撮像方法、及びct撮像サービス方法
JP2004065706A (ja) * 2002-08-08 2004-03-04 Ge Medical Systems Global Technology Co Llc 投影データ補正方法、画像生成方法およびx線ct装置
DE10326281B4 (de) * 2003-06-11 2005-06-16 Siemens Ag Verfahren zur Zuordnung von Marken und Verwendungen der Verfahren
JP4495926B2 (ja) 2003-07-01 2010-07-07 株式会社東芝 X線立体再構成処理装置、x線撮影装置、x線立体再構成処理方法及びx線立体撮影補助具
CA2455359C (en) * 2004-01-16 2013-01-08 Geotango International Corp. System, computer program and method for 3d object measurement, modeling and mapping from single imagery
CN101126725B (zh) * 2007-09-24 2010-12-15 舒嘉 采用x射线容积摄影实现图像重建的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1376947A (zh) * 2001-03-26 2002-10-30 舒嘉 X射线立体图像观察系统
CN1464301A (zh) * 2002-06-25 2003-12-31 舒嘉 X射线立体摄影实现图像重建的ct-解析法
DE10244180A1 (de) * 2002-09-23 2004-04-08 Siemens Ag Verfahren zur Bilderstellung in der Computertomographie eines periodisch bewegten Untersuchungsobjektes und CT-Gerät zur Durchführung des Verfahrens
FR2862790A1 (fr) * 2003-11-25 2005-05-27 Ge Medical Syst Sa Procede et dispositif pour l'observation d'un phenomene transitoire a l'aide d'un arceau vasculaire
FR2864669A1 (fr) * 2003-12-24 2005-07-01 Ge Med Sys Global Tech Co Llc Procede pour ameliorer la visualisation d'un vaisseau

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张祖勋 张剑清.高等学校测绘类教材 数字摄影测量学 1.武汉大学出版社,1997,137-138.
张祖勋 张剑清.高等学校测绘类教材 数字摄影测量学 1.武汉大学出版社,1997,137-138. *

Also Published As

Publication number Publication date
US20100195891A1 (en) 2010-08-05
RU2469298C2 (ru) 2012-12-10
CN101126725A (zh) 2008-02-20
DE112008002547T5 (de) 2010-08-05
RU2010144360A (ru) 2012-04-27
JP2010540893A (ja) 2010-12-24
WO2009043224A1 (en) 2009-04-09
US8687861B2 (en) 2014-04-01

Similar Documents

Publication Publication Date Title
Assali et al. Surveying and modeling of rock discontinuities by terrestrial laser scanning and photogrammetry: Semi-automatic approaches for linear outcrop inspection
Silver Determining shape and reflectance using multiple images
AU783085B2 (en) Multi-view image registration
CN100476345C (zh) 基于双目立体视觉技术的空间圆几何参数的测量方法
US20150279087A1 (en) 3d data to 2d and isometric views for layout and creation of documents
CN102003938B (zh) 大型高温锻件热态在位检测方法
Ganapathy Decomposition of transformation matrices for robot vision
DE102008041523B4 (de) Verfahren zur dreidimensionalen Messung und Vorrichtung zur dreidimensionalen Messung
KR100576228B1 (ko) 다시각 x선 입체 영상 방법 및 시스템
Chen et al. Telecentric stereo micro-vision system: Calibration method and experiments
Orteu et al. Multiple-camera instrumentation of a single point incremental forming process pilot for shape and 3D displacement measurements: methodology and results
US20130070048A1 (en) Formation Apparatus Using Digital Image Correlation
Ivanov et al. Computer stereo plotting for 3-D reconstruction of a maize canopy
EP2923174A2 (en) Point-cloud fusion
JP2005530153A (ja) 三次元座標を取得するための立体x線撮像装置
CN104867160B (zh) 一种用于摄像机内外参数标定的方向性标定靶标
CN100430690C (zh) 利用单数码相机自由拍摄进行物体三维测量的方法
CN102663810A (zh) 一种基于相位偏移扫描的三维人脸全自动建模方法
Ai et al. A robust photogrammetric processing method of low-altitude UAV images
Santos et al. Image-based 3D digitizing for plant architecture analysis and phenotyping.
Hernández et al. Autonomous underwater navigation and optical mapping in unknown natural environments
CN104656097B (zh) 基于旋转式二维激光三维重构系统的标定装置及方法
CN102032878B (zh) 基于双目立体视觉测量系统的精确在线测量方法
CN1264062C (zh) 一种多视角x射线立体成像的方法与系统
CN104427325B (zh) 快速集成图像生成方法及与用户交互的裸眼三维显示系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant