CN101884546A - 一种抑制c型臂断层成像中轨道振动伪影的装置与方法 - Google Patents
一种抑制c型臂断层成像中轨道振动伪影的装置与方法 Download PDFInfo
- Publication number
- CN101884546A CN101884546A CN 201010222418 CN201010222418A CN101884546A CN 101884546 A CN101884546 A CN 101884546A CN 201010222418 CN201010222418 CN 201010222418 CN 201010222418 A CN201010222418 A CN 201010222418A CN 101884546 A CN101884546 A CN 101884546A
- Authority
- CN
- China
- Prior art keywords
- prime
- type arm
- beta
- actual
- ray source
- 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
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开一种抑制C型臂断层成像中轨道振动伪影的装置与方法,属于医疗成像领域。所要解决的是抑制C型臂断层成像中X射线源轨道振动带来的图像伪影。所述装置在C型臂上安装三个位置传感器,所述方法包括在C型臂进行扫描并获得投影数据的过程中,用三个位置传感器记录每个投影角度下C型臂振动情况的步骤;根据C型臂振动情况,计算出轨道振动伪影抑制方法中的偏移参数(D′,β′,α)的步骤;最后重建出图像。本发明能有效抑制C型臂振动所带来的轨道振动图像伪影,并且成本低廉,容易实现,对CT图像重建时间影响不大。
Description
技术领域
本发明涉及一种抑制X射线计算机断层成像(CT)伪影的装置与方法,更具体地说,它涉及一种基于传统锥束重建算法的C型臂X射线计算机断层成像中抑制轨道振动伪影的装置与方法,属于医疗成像领域。
背景技术
C型臂X射线计算机断层成像(C-arm CT)技术,是运用C型臂平面探测器获得的二维X射线投影数据去重建CT图像,这是一种全新的图像成像技术。C型臂成像系统的主要组成部分是X射线产生部分(X射线源)、探测器部分(平面探测器)和C型臂构架部分。C型臂系统构架包括了一个支架和一个C型臂,C型臂上一端安装着探测器,另一端安装着X射线源。由于C型臂非封闭式的设计,在有限空间内可以进行最大角度范围的扫描,因此,C型臂CT检测角度的选择十分灵活。另外,C型臂CT比传统CT小巧,移动方便。由于上述原因,C型臂CT广泛用于临床诊断和手术评估中。更重要的是,C型臂CT可用于术中检测,进行手术导航,这是一般CT不能做到的。
当前C-arm CT系统中用到的图像重建算法主要是由Feldkamp等人提出的锥束重建算法,也称FDK算法,见参考文献[1]:Feldkamp,Davis,Kress.,Practical cone-beam algorithm,J.Opt.Soc.Am.A,1984,1(6):612-619。这个算法主要包含了数据加权修正、卷积运算和三维滤波反投影三步。FDK算法简单易行,是一种基于圆轨道扫描的近似重建算法——在任何测量分辨率下,重建结果和真实物体都会有一些偏差。但是,如果选择合适的锥角,这种偏差将会非常小。FDK算法由于其简易性,成为实际锥束重建算法中应用最广泛的算法。
C型臂CT在实际应用中,X射线源的扫描轨道会由于机械振动,引力凹陷,轴承不规则性等因素偏离预先设计的轨道,发生轨道振动现象,从而导致在重建后的CT图像中出现伪影,影响医生诊断。
发明内容
本发明的目的是克服现有技术中的不足,提供一种结构简单、成本低廉的装置和易于实现、高效率的方法,用以抑制C型臂CT中出现的轨道振动伪影。
本发明的技术方案如下:在C型臂上安装三个位置传感器,分别装在C型臂圆弧的中间处、C型臂的X射线源处以及平面探测器处,在C型臂进行扫描并获得投影数据的过程中,用三个传感器记录的信息可以得到C型臂实际扫描过程中振动后的实际扫描平面,以及每个实际投影角度β′下X射线源的实际位置S′(x′,y′,z′)。然后,将每个实际投影角度β′下X射线源的实际位置换算成轨道振动伪影抑制方法中的三个偏移参数(D′,β′,α)。最后,再运用轨道振动修正方法,抑制轨道振动带来的伪影,重建出CT图像。具体抑制方法如下:
第一步、C型臂上的X射线源和平面探测器围绕物体的中心O点沿着扫描轨道进行旋转扫描,记录投影数据P(β′,a,b)。
第二步,在实际扫描过程中,根据三个传感器记录的信息确定C型臂的实际扫描平面,以及每个实际投影角度下X射线源的实际位置S′(x′,y′,z′)。
第三步,由每个实际投影角度β′下实际X射线源的位置S’,确定三个偏移参数(D′,β′,α)。
第四步,根据偏移参数D′和β′,对投影数据P(β′,a,b)进行加权处理。
本发明的有益效果是:
(1)该装置与方法能有效抑制扫描过程中,C型臂振动所带来的轨道振动图像伪影;
(2)成本低廉,结构简单,在几乎不改变原有C型臂CT系统加工工艺的前提下,仅需要加入三个位置传感器;
(3)该方法数学形式简单,容易实现,且修正方法引入的冗余计算少,对CT图像重建时间影响不大。
附图说明
图1是在理想的无振动情况下,C型臂CT进行检测时的扫描示意图:
图2是在出现振动情况后,C型臂CT进行检测时的扫描示意图;
图3是图2的俯视图;
图4是图2的侧视图;
图5为本发明的重建结果;
图6为传统锥束重建方法的重建结果;
图7为无振动情况下的重建结果。
具体实施方式
下面结合附图与实施例对本发明作进一步详细描述。
本发明提供的抑制C型臂断层成像中轨道振动伪影的装置结构如图1所示,在C型臂1的两个圆弧自由端分别固定有X射线源5和平面探测器6,C型臂1上安装三个位置传感器,分别为第一传感器2、第二传感器3和第三传感器4,所述的第一传感器2装在C型臂1圆弧的中心处,第二传感器3安装在C型臂1的X射线源5处,第三传感器4安装在平面探测器6处。C型臂1上的X射线源5和平面探测器6围绕物体8的中心O点沿着理想圆形扫描轨道7进行旋转扫描。
基于上述的抑制C型臂断层成像中轨道振动伪影的装置,本发明还提供了一种抑制C型臂断层成像中轨道振动伪影的方法,具体如下:
第一步、C型臂1上的X射线源5和平面探测器6围绕物体8的中心O点进行旋转扫描,记录投影数据P(β′,a,b)。
如图1所示,X射线源5所在位置记为点S,则点S到物体1中心O的线段SO的长度为D。β为投影角度,即在中心平面(XOY)上,直线SO与Y轴所成的角度。所述的中心平面是指C型臂1在理想状态下(C型臂无振动情况)进行旋转扫描的轨道平面。
图1是在无振动情况下,C型臂1围绕物体8的中心O进行的扫描示意图,C型臂1的运动轨道是一个理想圆形。扫描过程就是把每个投影角度β下的投影数据P(β′,a,b)通过平面探测器6记录下来,用于CT图像重建。由于一般情况下,物体8的中心O点是固定的,记为O(x0,y0,z0),也可以把它作为建立空间坐标系的原点。每个投影角度β下的理想X射线源5的位置点S的坐标S(x,y,z)可以在扫描轨道定下后就计算得到。
第二步,在实际扫描过程中,根据三个传感器记录的信息确定C型臂的实际扫描平面,以及每个实际投影角度下X射线源的实际位置S′(x′,y′,z′)。
C型臂1在实际扫描过程中会出现振动情况,如图2所示,此时,实际X射线源的位置S’与无振动情况下X射线源的位置S之间会出现一个偏移,实际探测器的位置与无振动情况下探测器6的位置也会出现一个偏移。由于C型臂1是一个刚体,X射线源5和探测器6在X、Z方向上发生的偏移大小相同,但方向相反;而X射线源5和探测器6在Y方向上发生的偏移大小和方向都相同。此时,用三个传感器2、3、4分别记录下C型臂1圆弧中心处、X射线源处和平面探测器处的实际位置点。
因为C型臂1是刚体,这三个实际位置点即可确定实际扫描过程中(振动情况下)C型臂1的轨道平面。实际X射线源的位置S′的坐标记为S′(x′,y′,z′),而实际X射线源的位置S′在中心平面(XOY)上投影S″的坐标为S″(x′,y′,z)。
第三步,由每个实际投影角度β′下实际X射线源的位置S′,确定三个偏移参数(D′,β′,α)。
在轨道振动伪影抑制方法中,用参数(D′,β′,α)就能表示实际X射线源与它理想位置的偏移。
理想X射线源的位置S(x,y,z)与实际X射线源的位置S′(x′,y′,z′)的空间偏移为(Δx,Δy,Δz)。
实际X射线源的位置S′(x′,y′,z′)到物体中心O(x0,y0,z0)的距离D′可以直接由两点间距离公式计算得到(如图2):
实际投影角度β′可以由直线S″O与Y轴间夹角计算公式得到(如图3):
实际X射线源的位置S′与中心平面(XOY)所成的垂直偏移角度α可以由垂直三角形S′S″O计算得到(如图4):
第四步,根据偏移参数D′和β′,对投影数据P(β′,a,b)进行加权处理:
其中,β′表示实际投影角度,a为探测器上横坐标,b为探测器上纵坐标,并且
再用沿着纵坐标方向的滤波器hP(a)对加权后的投影数据P′(β′,a,b)进行滤波:
其中:f(x,y,z)为重建图像中的每个像素点的衰减系数,重建图像由它组成。式(6)中,U(x,y,β′,α,D′)=D′+xcosαsinβ′-ycosαcosβ′+zsinα,为积分加权因子。
在对本发明的验证实验中,分别采用传统锥束重建方法和采用本发明的轨道振动伪影抑制方法对存在振动情况的投影数据进行重建,图5为本发明的重建结果,图6为传统锥束重建方法的重建结果,图7为无振动情况下的重建结果,可以看出,本发明的重建结果与无振动情况下的重建结果相近,明显优于传统锥束重建方法的重建结果。且本发明方法在具体实现中结构简单,计算方便,相对于传统方法,并未增加太多复杂度。
Claims (3)
1.一种抑制C型臂断层成像中轨道振动伪影的装置,其特征在于:C型臂上安装三个位置传感器,分别为第一传感器、第二传感器和第三传感器,所述的第一传感器装在C型臂圆弧的中心处,第二传感器安装在C型臂的X射线源处,第三传感器安装在平面探测器处。
2.一种抑制C型臂断层成像中轨道振动伪影的方法,其特征在于:
第一步、C型臂上的X射线源和平面探测器围绕物体的中心O点沿着扫描轨道进行旋转扫描,记录投影数据P(β′,a,b);
第二步,在实际扫描过程中,根据三个传感器记录的信息确定C型臂的实际扫描平面,以及每个实际投影角度下X射线源的实际位置S′(x′,y′,z′);
三个传感器分别记录下C型臂圆弧中心处、X射线源处和平面探测器处的实际位置点,即得到C型臂的实际扫描平面;
理想X射线源的位置S(x,y,z)与实际X射线源的位置S′(x′,y′,z′)的空间偏移为(Δx,Δy,Δz);实际X射线源的位置S′的坐标记为S′(x′,y′,z′),而实际X射线源的位置S′在中心平面XOY上投影S″的坐标为S″(x′,y′,z);
第三步,由每个实际投影角度β′下实际X射线源的位置S′,确定三个偏移参数(D′,β′,α);实际X射线源的位置S′(x′,y′,z′)到物体中心O(x0,y0,z0)的距离D′:
实际投影角度β′由直线S″O与Y轴间夹角计算得到:
实际X射线源的位置S′与中心平面XOY所成的垂直偏移角度α由垂直三角形S′S″O计算得到:
第四步,根据偏移参数D′和β′,对投影数据P(β′,a,b)进行加权处理;加权处理后得到:
其中,β′表示实际投影角度,a为探测器上横坐标,
其中:f(x,y,z)为每点的衰减系数,重建图像由每点的衰减系数组成;式(6)中,U(x,y,β′,α,D′)=D′+xcosαsinβ′-ycosαcosβ′+zsinα,为积分加权因子。
3.根据权利要求2所述的一种抑制C型臂断层成像中轨道振动伪影的方法,其特征在于:第四步中加权处理后还要进行滤波的步骤,所述的滤波步骤沿着纵坐标方向的滤波器hP(a)对加权后的投影数据P′(β′,a,b)进行滤波,得到:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010222418 CN101884546B (zh) | 2010-06-30 | 2010-06-30 | 一种抑制c型臂断层成像中轨道振动伪影的装置与方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010222418 CN101884546B (zh) | 2010-06-30 | 2010-06-30 | 一种抑制c型臂断层成像中轨道振动伪影的装置与方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101884546A true CN101884546A (zh) | 2010-11-17 |
CN101884546B CN101884546B (zh) | 2013-02-27 |
Family
ID=43070754
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010222418 Expired - Fee Related CN101884546B (zh) | 2010-06-30 | 2010-06-30 | 一种抑制c型臂断层成像中轨道振动伪影的装置与方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101884546B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102727229A (zh) * | 2011-03-31 | 2012-10-17 | 上海西门子医疗器械有限公司 | Ct设备中投影角的获取方法、装置、及ct设备 |
JP2016523685A (ja) * | 2013-07-15 | 2016-08-12 | シロナ・デンタル・システムズ・ゲゼルシャフト・ミット・ベシュレンクテル・ハフツング | X線画像を作成する方法 |
CN106821403A (zh) * | 2017-01-19 | 2017-06-13 | 深圳先进技术研究院 | C型臂系统旋转角度校准装置及c型臂系统校准方法 |
WO2018133002A1 (zh) * | 2017-01-19 | 2018-07-26 | 深圳先进技术研究院 | C型臂系统旋转角度校准装置和c型臂系统校准方法 |
CN113570710A (zh) * | 2021-08-09 | 2021-10-29 | 上海优医基医疗影像设备有限公司 | 一种具有c型臂抖动校正的三维成像系统及方法 |
CN113812971A (zh) * | 2021-08-27 | 2021-12-21 | 浙江大学 | 一种多自由度四维双能锥束ct成像系统及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020025017A1 (en) * | 1999-06-17 | 2002-02-28 | Stergios Stergiopoulos | Method for tracing organ motion and removing artifacts for computed tomography imaging systems |
CN1883390A (zh) * | 2005-06-23 | 2006-12-27 | Ge医疗系统环球技术有限公司 | X射线ct系统 |
DE102008008053A1 (de) * | 2008-02-08 | 2009-08-27 | Siemens Aktiengesellschaft | Verfahren zur Aufnahme einer Bildsequenz mit einem digitalen Röntgendetektor und Röntgengerät |
-
2010
- 2010-06-30 CN CN 201010222418 patent/CN101884546B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020025017A1 (en) * | 1999-06-17 | 2002-02-28 | Stergios Stergiopoulos | Method for tracing organ motion and removing artifacts for computed tomography imaging systems |
CN1883390A (zh) * | 2005-06-23 | 2006-12-27 | Ge医疗系统环球技术有限公司 | X射线ct系统 |
DE102008008053A1 (de) * | 2008-02-08 | 2009-08-27 | Siemens Aktiengesellschaft | Verfahren zur Aufnahme einer Bildsequenz mit einem digitalen Röntgendetektor und Röntgengerät |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102727229A (zh) * | 2011-03-31 | 2012-10-17 | 上海西门子医疗器械有限公司 | Ct设备中投影角的获取方法、装置、及ct设备 |
CN102727229B (zh) * | 2011-03-31 | 2014-07-30 | 上海西门子医疗器械有限公司 | Ct设备中投影角的获取方法、装置、及ct设备 |
JP2016523685A (ja) * | 2013-07-15 | 2016-08-12 | シロナ・デンタル・システムズ・ゲゼルシャフト・ミット・ベシュレンクテル・ハフツング | X線画像を作成する方法 |
CN106821403A (zh) * | 2017-01-19 | 2017-06-13 | 深圳先进技术研究院 | C型臂系统旋转角度校准装置及c型臂系统校准方法 |
WO2018133002A1 (zh) * | 2017-01-19 | 2018-07-26 | 深圳先进技术研究院 | C型臂系统旋转角度校准装置和c型臂系统校准方法 |
CN113570710A (zh) * | 2021-08-09 | 2021-10-29 | 上海优医基医疗影像设备有限公司 | 一种具有c型臂抖动校正的三维成像系统及方法 |
CN113570710B (zh) * | 2021-08-09 | 2022-08-12 | 上海优医基医疗影像设备有限公司 | 一种具有c型臂抖动校正的三维成像系统及方法 |
CN113812971A (zh) * | 2021-08-27 | 2021-12-21 | 浙江大学 | 一种多自由度四维双能锥束ct成像系统及方法 |
CN113812971B (zh) * | 2021-08-27 | 2023-10-13 | 浙江大学 | 一种多自由度四维双能锥束ct成像系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101884546B (zh) | 2013-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101884546B (zh) | 一种抑制c型臂断层成像中轨道振动伪影的装置与方法 | |
CN102973291B (zh) | C型臂半精确滤波反投影断层成像方法 | |
CN104809750B (zh) | 一种直线扫描ct系统及图像重建方法 | |
JP6109524B2 (ja) | X線コンピュータ断層撮影装置及び画像補正方法 | |
JP4668161B2 (ja) | 直線軌跡走査を採用した画像復元システム及び方法 | |
CN101133962B (zh) | 再现三维立体图像的方法和x射线设备 | |
US8218717B2 (en) | Reconstruction of 3D image datasets from x-ray and cone-beam data | |
CN100570343C (zh) | 大视场三维ct成像方法 | |
JP2014533168A (ja) | 拡張された視野を備えるcアーム型システム | |
CN102028490A (zh) | 从在遍历轨迹时拍摄的x射线投影确定图像的方法和设备 | |
JP2005504571A (ja) | 多機能コーンビーム結像装置とその方法 | |
KR20060041709A (ko) | 화상 재구성 방법 및 x선 ct 장치 | |
CN102005031A (zh) | 一种mri系统中k空间采样数据的运动伪影消除方法及装置 | |
Bronnikov | Virtual alignment of x-ray cone-beam tomography system using two calibration aperture measurements | |
CN109978988A (zh) | 用于重建三维图像数据组的方法、双平面x射线装置 | |
US20130083986A1 (en) | Method and system for reconstruction of tomographic images | |
CN102488528B (zh) | 一种层析成像几何参数的校准方法 | |
Baek et al. | Use of sphere phantoms to measure the 3D MTF of FDK reconstructions | |
JP2007198866A (ja) | 広義サドルコーンビームct装置および3次元再構成法 | |
CN112120722B (zh) | X射线断层合成装置、图像处理装置以及记录介质 | |
JP3913798B2 (ja) | 画像再構成処理装置 | |
CN104793160B (zh) | 一种减少欠采样磁共振成像的频率混迭效应的方法 | |
JP3878788B2 (ja) | トモシンセシスを行う方法及び装置 | |
JP3825492B2 (ja) | 画像再構成処理装置及びx線ct装置 | |
JP2009237704A (ja) | 3次元画像化方法およびx線断層像撮影装置 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130227 Termination date: 20150630 |
|
EXPY | Termination of patent right or utility model |