CN103258335A - 基于普通c臂机的脊柱横断面图像合成方法 - Google Patents
基于普通c臂机的脊柱横断面图像合成方法 Download PDFInfo
- Publication number
- CN103258335A CN103258335A CN2012100387782A CN201210038778A CN103258335A CN 103258335 A CN103258335 A CN 103258335A CN 2012100387782 A CN2012100387782 A CN 2012100387782A CN 201210038778 A CN201210038778 A CN 201210038778A CN 103258335 A CN103258335 A CN 103258335A
- Authority
- CN
- China
- Prior art keywords
- fan
- image
- arm
- virtual
- light 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
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
一种基于普通C臂机的脊柱横断面图像合成方法,包括以下步骤:步骤一、对普通C臂机获得的2D锥束投影图像进行预处理,得到1D扇束投影数据;步骤二、对步骤一得到的扇束进行虚拟旋转再处理,计算得到虚拟扇束投影数据;步骤三、对2D断层图像合成方法中的相关参数进行调整,采用调整后的方法基于虚拟扇束投影数据合成脊柱横断面图像。本发明合成方法,原理简单,执行方便,除C臂机本身外无需购置或定制其他额外装置,可用于合成在临床诊断、手术规划、引导手术操作方面都有重要意义的脊柱横断面图像。
Description
技术领域
本发明涉及计算机成像技术,特别涉及一种基于普通C臂机的脊柱横断面图像合成方法。
背景技术
在医院外科临床应用中,尤其是脊柱外科领域,C臂机是一种常用的设备。借助C臂机获取的图像,医生可对患者脊柱的病损状况进行判断,还可对患者进行手术规划和植入锥弓根螺钉等手术操作。但是,通过普通C臂机仅可获取患者的投影图像,而临床上对病情的判断及引导手术操作意义重大的脊柱横断面图像却不能通过C臂机直接获取。
计算机断层成像技术利用X线对患者进行扫描获取投影图像序列,利用特定的算法以投影图像序列为基础合成人体横断面图像。经对现有技术的文献检索发现,Feldkamp等人提出的FDK算法(Practical cone-beam reconstruction,J.Opt.Soc.Am.A1,612-619,1984)适用基于扇束和锥束投影的断层图像合成,但其要求具备360°范围内全扫描路径投影数据,而普通C臂旋转扫描角度范围一般不超过180°,不能直接采用FDK算法。Noo和Kudo等人提出的超短扫描路径断层图像合成算法(Image reconstruction from fan-beam projections on less than ashort scan,J.Phys.Med.Biol.47,2525-2546,2002;New Super-Short-ScanAlgorithms for Fan-Beam and Cone-Beam Reconstruction,IEEE Nuclear Science andMedical Imaging Symposium,Record,Norfolk,USA,2003,2:902-906)允许扫描范围低于180°,但其要求设备必须满足等中心旋转扫描条件,而普通C臂锥束中轴线一般不通过C臂旋转中心,其并不满足等中心旋转扫描条件,也不能直接采用这些算法。
发明内容
本发明的目的,就在于克服上述已有技术的不足,提供一种基于普通C臂机的脊柱横断面图像合成方法,以便依托目前国内大多数医院都具备的普通C臂机合成在临床上意义重大的横断面图像。
为实现上述目的,本发明采用了以下技术方案:一种基于普通C臂机的脊柱横断面图像合成方法,包括以下步骤:
步骤一、对普通C臂机获得的2D锥束投影图像进行预处理,得到1D扇束投影数据;
步骤二、对步骤一得到的扇束进行虚拟旋转再处理,计算得到虚拟扇束投影数据;
步骤三、对2D断层图像合成方法中的相关参数进行调整,采用调整后的方法基于虚拟扇束投影数据合成脊柱横断面图像。
步骤一中所述的预处理,是指首先对2D锥束投影图像进行中值滤波除噪,然后提取出图像中间的3列或2列像素的灰度值,对提取的像素的灰度值取均值构成1D扇束投影数据;
所述1D扇束投影与普通C臂机获得的2D锥束投影的关系为:扇束光源点与C臂X线锥束光源点重合,扇束中轴线与C臂X线锥束中轴线重合,扇束光源点到其成像直线之间的距离等于C臂X线锥束光源点到其成像平面之间的距离,扇束成像直线位于C臂X线锥束成像平面上,扇束的扇形角等于C臂X线锥束的锥角γm。
步骤二中所述的再处理,是指对步骤一所得扇束绕其光源点进行虚拟旋转得到虚拟扇束,通过计算后得到对应的虚拟扇束投影数据;
所述虚拟扇束与扇束的关系为:虚拟扇束光源点到其成像直线之间的距离等于扇束光源点到其成像直线之间的距离D,虚拟扇束光源点到其成像直线的垂线通过C臂旋转中心;对虚拟扇束投影和扇束投影分别建立笛卡尔坐标系,令两个坐标系的原点分别位于光源点在各自成像直线上的垂足,且令虚拟扇束投影坐标轴正向与扇束投影坐标轴正向所夹角度为θ,则两个坐标系下投影像素坐标及灰度值对应关系为:
上式中,u为扇束投影在其笛卡尔坐标系中的像素坐标,uv为虚拟扇束投影在其笛卡尔坐标系中的像素坐标,gc(λ,u)为扇束投影在对应坐标处的像素灰度值,gv(λ,uv)为虚拟扇束投影在对应坐标处的像素灰度值,λ为C臂旋转角度。
步骤三中所述的对2D断层图像合成方法中的相关参数进行调整,具体包括:将光源点沿圆弧路径运动的半径R调整为光源点到C臂旋转中心的距离,将对扇束投影数据进行滤波处理的区间下限ul调整为-D.tan(γm-θ),上限uu调整为D.tan(γm+θ),上述各式中,γm为扇束的扇形角,D为扇束光源点到其成像直线之间的距离,θ为虚拟扇束投影坐标轴正向与扇束投影坐标轴正向所夹角度。
所述提取出图像中间的3列或2列像素的灰度值,是指当图像像素列数n为奇数时,提取出图像中间的3列像素的灰度值,当图像像素列数n为偶数时,提取出图像中间的2列像素的灰度值。
本发明的有益效果是,采用该方法,依托目前国内大多数医院都具备的普通非等中心旋转扫描C臂机,可合成在临床诊断、手术规划、引导手术操作方面都有重要意义的脊柱横断面图像。
本发明基于普通C臂机的脊柱横断面图像合成方法原理简单,执行方便,除C臂机本身外无需购置或定制其他额外装置。
附图说明
图1为本发明基于普通C臂机的脊柱横断面图像合成方法的流程图;
图2为本发明步骤一中锥束投影与扇束投影的关系示意图;
图3为本发明步骤二中扇束与虚拟扇束的关系示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
参见图1,本发明基于普通C臂机的脊柱横断面图像合成方法包含三个步骤:
第1步,对普通C臂2D锥束投影图像进行预处理,得到1D扇束投影数据。首先对C臂在每个旋转角度λ下获取的2D锥束投影图像gm(λ,u,v)进行中值滤波除噪(u,v分别为图像中像素的行、列坐标),得到除噪后投影图像ge(λ,u,v)。然后再从图像ge(λ,u,v)中提取出其中间3列像素的灰度值(若图像像素列数n为奇数):gc1(λ,u)=ge(λ,u,v(((n+1)/2)-1),gc2(λ,u)=ge(λ,u,v((n+1)/2)),gc3(λ,u)=ge(λ,u,v(((n+1)/2)+1),或中间2列像素的灰度值(若图像像素列数n为偶数):gcl(λ,u)=ge(λ,u,v(n/2)),gc2(λ,u)=ge(λ,u,v((n/2)+1))。最后对中间3列(或2列)像素的灰度值取均值构成1D扇束投影数据gc(λ,u)。配合参见图2,所得扇束光源点21与C臂X线锥束1的光源点重合,所得扇束中轴线23与C臂X线锥束中轴线重合,所得扇束光源点21到其成像直线25之间的距离等于C臂X线锥束光源点到其成像平面15之间的距离,所得扇束成像直线25位于C臂X线锥束成像平面15上,所得扇束的扇形角22等于C臂X线锥束的锥角,设其为γm。
第2步,配合参见图2和图3,对扇束2进行虚拟旋转再处理,得到虚拟扇束3,并计算得到虚拟扇束投影数据。虚拟扇束光源点与扇束光源点21重合,虚拟扇束光源点到其成像直线35之间的距离等于扇束光源点21到其成像直线25之间的距离,设其为D,虚拟扇束光源点到其成像直线35的垂线33通过C臂旋转中心4。设θ为虚拟扇束投影成像直线35之垂线与扇束投影成像直线25之垂线所夹锐角。对虚拟扇束投影建立笛卡尔坐标系,令原点位于其光源点在其成像直线35上的垂足34。对扇束投影建立笛卡尔坐标系,令原点位于其光源点21在其成像直线25上的垂足24。令虚拟扇束3之投影坐标轴正向与扇束2之投影坐标轴正向所夹角度为θ,令两个坐标系下投影像素坐标及灰度值对应关系为:
上式中,u和uv分别为扇束投影及虚拟扇束投影在各自坐标系中的像素坐标,g(u)和gv(uv)分别为扇束投影及虚拟扇束投影在对应坐标处的像素灰度值,λ为C臂旋转角度。
第3步,对2D断层图像合成方法中的相关参数进行调整,采用调整后的方法基于虚拟扇束投影数据合成横断面图像。配合参见图2和图3,将光源点沿圆弧路径运动的半径R调整为光源点21到C臂旋转中心4之间的距离,将扇束投影数据滤波处理区间的下限ul调整为-D.tan(γm-θ),上限uu调整为D.tan(γm+θ)。采用公式 基于虚拟扇束投影数据合成横断面图像。和分别表示横断面图像像素坐标及像素灰度值。γm的含义见步骤1,D和θ的含义见步骤2。其他中间计算公式具体如下:
d=10°
以上所述仅为本发明的较佳实施例而已,并非用来限定本发明的实施范围。任何所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰,因此本发明的保护范围应当视权利要求书所界定范围为准。
Claims (5)
1.一种基于普通C臂机的脊柱横断面图像合成方法,其特征在于:包括以下步骤:
步骤一、对普通C臂机获得的2D锥束投影图像进行预处理,得到1D扇束投影数据;
步骤二、对步骤一得到的扇束进行虚拟旋转再处理,计算得到虚拟扇束投影数据;
步骤三、对2D断层图像合成方法中的相关参数进行调整,采用调整后的方法基于虚拟扇束投影数据合成脊柱横断面图像。
2.如权利要求1所述的基于普通C臂机的脊柱横断面图像合成方法,其特征在于:步骤一中所述的预处理,是指首先对2D锥束投影图像进行中值滤波除噪,然后提取出图像中间的3列或2列像素的灰度值,对提取的像素的灰度值取均值构成1D扇束投影数据;
所述1D扇束投影与普通C臂机获得的2D锥束投影的关系为:扇束光源点与C臂X线锥束光源点重合,扇束中轴线与C臂X线锥束中轴线重合,扇束光源点到其成像直线之间的距离等于C臂X线锥束光源点到其成像平面之间的距离,扇束成像直线位于C臂X线锥束成像平面上,扇束的扇形角等于C臂X线锥束的锥角γm。
3.如权利要求1或2所述的基于普通C臂机的脊柱横断面图像合成方法,其特征在于:步骤二中所述的再处理,是指对步骤一所得扇束绕其光源点进行虚拟旋转得到虚拟扇束,通过计算后得到对应的虚拟扇束投影数据;
所述虚拟扇束与扇束的关系为:虚拟扇束光源点到其成像直线之间的距离等于扇束光源点到其成像直线之间的距离D,虚拟扇束光源点到其成像直线的垂线通过C臂旋转中心;对虚拟扇束投影和扇束投影分别建立笛卡尔坐标系,令两个坐标系的原点分别位于光源点在各自成像直线上的垂足,且令虚拟扇束投影坐标轴正向与扇束投影坐标轴正向所夹角度为θ,则两个坐标系下投影像素坐标及灰度值对应关系为:
上式中,u为扇束投影在其笛卡尔坐标系中的像素坐标,uv为虚拟扇束投影在其笛卡尔坐标系中的像素坐标,gc (λ,u)为扇束投影在对应坐标处的像素灰度值,gv(λ,uv)为虚拟扇束投影在对应坐标处的像素灰度值,λ为C臂旋转角度。
4.如权利要求1所述的基于普通C臂机的脊柱横断面图像合成方法,其特征在于:步骤三中所述的对2D断层图像合成方法中的相关参数进行调整,具体包括:将光源点沿圆弧路径运动的半径R调整为光源点到C臂旋转中心的距离,将对扇束投影数据进行滤波处理的区间下限ul调整为-D.tan(γm-θ),上限uu调整为D.tan(γm+θ),上述各式中,γm为扇束的扇形角,D为扇束光源点到其成像直线之间的距离,θ为虚拟扇束投影坐标轴正向与扇束投影坐标轴正向所夹角度。
5.如权利要求2所述的基于普通C臂机的脊柱横断面图像合成方法,其特征在于:所述提取出图像中间的3列或2列像素的灰度值,是指当图像像素列数n为奇数时,提取出图像中间的3列像素的灰度值,当图像像素列数n为偶数时,提取出图像中间的2列像素的灰度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210038778.2A CN103258335B (zh) | 2012-02-20 | 2012-02-20 | 基于普通c臂机的脊柱横断面图像合成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210038778.2A CN103258335B (zh) | 2012-02-20 | 2012-02-20 | 基于普通c臂机的脊柱横断面图像合成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103258335A true CN103258335A (zh) | 2013-08-21 |
CN103258335B CN103258335B (zh) | 2015-12-16 |
Family
ID=48962228
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210038778.2A Expired - Fee Related CN103258335B (zh) | 2012-02-20 | 2012-02-20 | 基于普通c臂机的脊柱横断面图像合成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103258335B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4984159A (en) * | 1988-08-09 | 1991-01-08 | General Electric Company | Method and apparatus for estimating elliptical body contours in fan beam computed tomographic systems |
US6196715B1 (en) * | 1959-04-28 | 2001-03-06 | Kabushiki Kaisha Toshiba | X-ray diagnostic system preferable to two dimensional x-ray detection |
CN1637781A (zh) * | 2003-12-26 | 2005-07-13 | 佳能株式会社 | 图像处理装置、图像处理方法以及图像处理系统 |
CN1642479A (zh) * | 2002-03-12 | 2005-07-20 | 爱克斯康特公司 | 电离辐射的以扫描为基础的探测中的曝光控制 |
CN101842052A (zh) * | 2007-07-19 | 2010-09-22 | 北卡罗来纳大学查珀尔希尔分校 | 固定x射线数字化乳房断层合成系统和相关方法 |
-
2012
- 2012-02-20 CN CN201210038778.2A patent/CN103258335B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6196715B1 (en) * | 1959-04-28 | 2001-03-06 | Kabushiki Kaisha Toshiba | X-ray diagnostic system preferable to two dimensional x-ray detection |
US4984159A (en) * | 1988-08-09 | 1991-01-08 | General Electric Company | Method and apparatus for estimating elliptical body contours in fan beam computed tomographic systems |
CN1642479A (zh) * | 2002-03-12 | 2005-07-20 | 爱克斯康特公司 | 电离辐射的以扫描为基础的探测中的曝光控制 |
CN1637781A (zh) * | 2003-12-26 | 2005-07-13 | 佳能株式会社 | 图像处理装置、图像处理方法以及图像处理系统 |
CN101842052A (zh) * | 2007-07-19 | 2010-09-22 | 北卡罗来纳大学查珀尔希尔分校 | 固定x射线数字化乳房断层合成系统和相关方法 |
Non-Patent Citations (3)
Title |
---|
张利军等: "基于C型臂的Tomosynthesis快速重建方法", 《北京航空航天大学学报》 * |
闫士举等: "C形臂X线投影图像3D建模及其应用", 《中国医学物理学杂志》 * |
闫士举等: "C形臂X线投影图像3D建模及其应用", 《中国医学物理学杂志》, vol. 27, no. 5, 30 September 2010 (2010-09-30), pages 2109 - 2112 * |
Also Published As
Publication number | Publication date |
---|---|
CN103258335B (zh) | 2015-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Flohr et al. | Multi–detector row CT systems and image-reconstruction techniques | |
US6298110B1 (en) | Cone beam volume CT angiography imaging system and method | |
JP2020527087A (ja) | 超音波拓本技術に基づく脊椎画像生成システム及び脊柱手術用のナビゲーション・位置確認システム | |
von Atzigen et al. | HoloYolo: A proof‐of‐concept study for marker‐less surgical navigation of spinal rod implants with augmented reality and on‐device machine learning | |
JP2017532165A (ja) | 脊椎不安定性を測定して評価するためのシステム及び方法 | |
US20060109952A1 (en) | Fan-beam and cone-beam image reconstruction using filtered backprojection of differentiated projection data | |
CN103679768B (zh) | 以加权反投影重建ct图像数据的方法、计算单元和系统 | |
CN105640575A (zh) | 医学图像生成装置、方法和程序 | |
JP2013244211A (ja) | 医用画像処理装置、医用画像処理方法及び制御プログラム | |
US20160371860A1 (en) | System and method for generating a ct slice image | |
US20120141008A1 (en) | Method and device for processing 3-d image data of a skull | |
CN105894508A (zh) | 一种医学图像自动定位质量的评估方法 | |
CN105488826A (zh) | 一种基于fbp的能谱ct迭代成像方法和成像系统 | |
CN103366397A (zh) | 基于c形臂2d投影图像的脊柱3d模型构建方法 | |
WO2013132934A1 (ja) | Ct画像生成装置及び方法、ct画像生成システム | |
JP2011031039A (ja) | 傾斜画像スキャン方法、再構成方法及びその装置 | |
JP2001238877A (ja) | X線画像を形成するための方法及び装置 | |
JP5177606B1 (ja) | 3次元超音波画像作成方法およびプログラム | |
CN109350059B (zh) | 用于肘部自动对准的组合的转向引擎和界标引擎 | |
CN106934810A (zh) | 一种脊椎矫正装置 | |
TWI610661B (zh) | 骨結構檢測方法與系統 | |
CN106913353A (zh) | 用于大体积拍摄的断层成像设备和方法 | |
KR20110136487A (ko) | 의료용 융합영상 획득장치 및 획득방법 | |
Yu et al. | Axially extended-volume C-arm CT using a reverse helical trajectory in the interventional room | |
US20110109630A1 (en) | Preparation and Presentation of Patient-Individual Panoramic Visualizations |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151216 Termination date: 20190220 |