CN103258336B - 基于扇束虚拟平移的脊柱横断面图像合成方法 - Google Patents
基于扇束虚拟平移的脊柱横断面图像合成方法 Download PDFInfo
- Publication number
- CN103258336B CN103258336B CN201210038779.7A CN201210038779A CN103258336B CN 103258336 B CN103258336 B CN 103258336B CN 201210038779 A CN201210038779 A CN 201210038779A CN 103258336 B CN103258336 B CN 103258336B
- Authority
- CN
- China
- Prior art keywords
- fan
- virtual
- light source
- image
- source point
- 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.)
- Expired - Fee Related
Links
Landscapes
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
Abstract
一种基于扇束虚拟平移的脊柱横断面图像合成方法,包括以下步骤:步骤一、对普通C臂机获得的2D锥束投影图像进行预处理,得到1D扇束投影数据;步骤二、对步骤一得到的扇束进行虚拟平移再处理,计算得到虚拟扇束投影数据;步骤三、对2D断层图像合成方法中的相关参数进行调整,采用调整后的方法基于虚拟扇束投影数据合成脊柱横断面图像。本发明合成方法,原理简单,执行方便,除C臂机本身外无需购置或定制其他额外装置,可用于合成在临床诊断、手术规划、引导手术操作方面都有重要意义的脊柱横断面图像。
Description
技术领域
本发明涉及计算机成像技术,特别涉及一种基于扇束虚拟平移的脊柱横断面图像合成方法。
背景技术
在医院外科临床应用中,尤其是脊柱外科领域,C臂机是一种常用的设备。借助C臂机获取的图像,医生可对患者脊柱的病损状况进行判断,还可对患者进行手术规划和植入锥弓根螺钉等手术操作。但是,通过普通C臂机仅可获取患者的投影图像,而临床上对病情的判断及引导手术操作意义重大的脊柱横断面图像却不能通过C臂机直接获取。
计算机断层成像技术利用X线对患者进行扫描获取投影图像序列,利用特定的算法以投影图像序列为基础合成人体横断面图像。经对现有技术的文献检索发现,Feldkamp等人提出的FDK算法(Practicalcone-beamreconstruction,J.Opt.Soc.Am.A1,612-619,1984)适用基于扇束和锥束投影的断层图像合成,但其要求具备360°范围内全扫描路径投影数据,而普通C臂旋转扫描角度范围一般不超过180°,不能直接采用FDK算法。Noo和Kudo等人提出的超短扫描路径断层图像合成算法(Imagereconstructionfromfan-beamprojectionsonlessthanashortscan,J.Phys.Med.Biol.47,2525-2546,2002;NewSuper-Short-ScanAlgorithmsforFan-BeamandCone-BeamReconstruction,IEEENuclearScienceandMedicalImagingSymposium,Record,Norfolk,USA,2003,2:902-906)允许扫描范围低于180°,但其要求设备必须满足等中心旋转扫描条件,而普通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。
步骤二中所述的再处理,是指对步骤一所得扇束进行虚拟平移得到虚拟扇束,虚拟平移方向与步骤一所得扇束中轴线垂直且朝向C臂旋转中心,通过计算后得到对应的虚拟扇束投影数据;
所述虚拟扇束与扇束的关系为:虚拟扇束光源点到其成像直线之间的距离等于扇束光源点到其成像直线之间的距离D,虚拟扇束光源点到其成像直线的垂线通过C臂旋转中心;设R为扇束光源点到C臂旋转中心之间的距离,扇束光源点和C臂旋转中心之连线与扇束中轴线之间的夹角为θ;对虚拟扇束投影和扇束投影分别建立笛卡尔坐标系,令两个坐标系的原点分别位于光源点在各自成像直线上的垂足,且令虚拟扇束投影坐标轴正向与扇束投影坐标轴正向相同,则两个坐标系下投影像素坐标及灰度值对应关系为:
上式中,u和uv分别为扇束投影和虚拟扇束投影在各自坐标系中的像素坐标,gc(λ,u)和gv(λ,uv)分别为扇束投影及虚拟扇束投影在对应坐标处的像素灰度值,λ为C臂旋转角度。
步骤三中所述的对2D断层图像合成方法中的相关参数进行调整,具体包括:光源点沿圆弧路径运动的半径调整为R.cos(θ),对扇束投影数据进行滤波处理的区间下限ul调整为R.sin(θ)-D.tan(γm),上限uu调整为R.sin(θ)+D.tan(γm),上述各式中,γm为扇束的扇形角,D为扇束光源点到其成像直线之间的距离,θ为扇束光源点和C臂旋转中心之连线与扇束中轴线之间的夹角。
所述提取出图像中间的3列或2列像素的灰度值,是指当图像像素列数n为奇数时,提取出图像中间的3列像素的灰度值,当图像像素列数n为偶数时,提取出图像中间的2列像素的灰度值。
本发明的有益效果是,采用该方法,依托目前国内大多数医院都具备的普通非等中心旋转扫描C臂机,可合成在临床诊断、手术规划、引导手术操作方面都有重要意义的脊柱横断面图像。
本发明基于扇束虚拟平移的脊柱横断面图像合成方法原理简单,执行方便,除C臂机本身外无需购置或定制其他额外装置。
附图说明
图1为本发明基于扇束虚拟平移的脊柱横断面图像合成方法的流程图;
图2为本发明步骤一中锥束投影与扇束投影的关系示意图;
图3为本发明步骤二中扇束与虚拟扇束的关系示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
参见图1,本发明基于扇束虚拟平移的脊柱横断面图像合成方法包含三个步骤:
第1步,对普通C臂2D锥束投影图像进行预处理,得到1D扇束投影数据。首先对C臂在每个旋转角度λ下获取的2D锥束投影图像gm(λ,u,v)进行中值滤波除噪(u,v分别为图像中像素的行、列坐标),得到除噪后投影图像ge(λ,u,v)。然后再从图像ge(λ,u,v)中提取出其中间3列像素的灰度值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)(若图像像素列数n为奇数),或中间2列像素的灰度值gc1(λ,u)=ge(λ,u,v(n/2)),gc2(λ,u)=ge(λ,u,v((n/2)+1))(若图像像素列数n为偶数)。最后对中间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,并计算得到虚拟扇束投影数据。虚拟平移方向与步骤一所得扇束中轴线垂直且朝向C臂旋转中心。虚拟扇束光源点31到其成像直线35之间的距离等于扇束光源点21到其成像直线25之间的距离,设其为D,虚拟扇束光源点31到其成像直线35的垂线通过C臂旋转中心4。设R为扇束光源点21到C臂旋转中心4之间的距离,扇束光源点21和C臂旋转中心4之连线与扇束中轴线23之间的夹角为θ。对虚拟扇束投影建立笛卡尔坐标系,令原点位于其光源点31在其成像直线35上的垂足34。对扇束投影建立笛卡尔坐标系,令原点位于其光源点21在其成像直线25上的垂足24。令虚拟扇束3之投影坐标轴正向与扇束2之投影坐标轴正向相同,令两个坐标系下投影像素坐标及灰度值对应关系为:
上式中,u和uv分别为扇束投影及虚拟扇束投影在各自坐标系中的像素坐标,g(u)和gv(uv)分别为扇束投影及虚拟扇束投影在对应坐标处的像素灰度值,λ为C臂旋转角度。
第3步,对2D断层图像合成方法中的相关参数进行调整,采用调整后的方法基于虚拟扇束投影数据合成横断面图像。配合参见图2和图3,将光源点沿圆弧路径运动的半径调整为R.cos(θ),对将扇束投影数据滤波处理区间的下限ul调整为R.sin(θ)-D.tan(γm),上限uu调整为R.sin(θ)+D.tan(γm)。
采用公式 基于虚拟扇束投影数据合成横断面图像。和分别表示横断面图像像素坐标及像素灰度值。γm的含义见步骤一。D,R和θ的含义见步骤二。其他中间计算公式具体如下:
以上所述仅为本发明的较佳实施例而已,并非用来限定本发明的实施范围。任何所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰,因此本发明的保护范围应当视权利要求书所界定范围为准。
Claims (2)
1.一种基于扇束虚拟平移的脊柱横断面图像合成方法,其特征在于:包括以下步骤:
步骤一、对普通C臂机获得的2D锥束投影图像进行预处理,得到1D扇束投影数据;
步骤二、对步骤一得到的扇束进行虚拟平移再处理,计算得到虚拟扇束投影数据;
步骤三、对2D断层图像合成方法中的相关参数进行调整,采用调整后的方法基于虚拟扇束投影数据合成脊柱横断面图像;
步骤一中所述的预处理,是指首先对2D锥束投影图像进行中值滤波除噪,当图像像素列数n为奇数时,提取出图像中间的3列像素的灰度值,当图像像素列数n为偶数时,提取出图像中间的2列像素的灰度值,对提取的像素的灰度值取均值构成1D扇束投影数据;
所述1D扇束投影与普通C臂机获得的2D锥束投影的关系为:扇束光源点与C臂X线锥束光源点重合,扇束中轴线与C臂X线锥束中轴线重合,扇束光源点到其成像直线之间的距离等于C臂X线锥束光源点到其成像平面之间的距离,扇束成像直线位于C臂X线锥束成像平面上,扇束的扇形角等于C臂X线锥束的锥角γm;
步骤二中所述的再处理,是指对步骤一所得扇束进行虚拟平移得到虚拟扇束,虚拟平移方向与步骤一所得扇束中轴线垂直且朝向C臂旋转中心,通过计算后得到对应的虚拟扇束投影数据;
所述虚拟扇束与扇束的关系为:虚拟扇束光源点到其成像直线之间的距离等于扇束光源点到其成像直线之间的距离D,虚拟扇束光源点到其成像直线的垂线通过C臂旋转中心;设R为扇束光源点到C臂旋转中心之间的距离,扇束光源点和C臂旋转中心之连线与扇束中轴线之间的夹角为θ;对虚拟扇束投影和扇束投影分别建立笛卡尔坐标系,令两个坐标系的原点分别位于虚拟扇束光源点和扇束光源点在各自成像直线上的垂足,且令虚拟扇束投影坐标轴正向与扇束投影坐标轴正向相同,则两个坐标系下投影像素坐标及灰度值对应关系为:
上式中,u和uv分别为扇束投影和虚拟扇束投影在各自坐标系中的像素坐标,gc(λ,u)和gv(λ,uv)分别为扇束投影及虚拟扇束投影在对应坐标处的像素灰度值,λ为C臂旋转角度。
2.如权利要求1所述的基于扇束虚拟平移的脊柱横断面图像合成方法,其特征在于:步骤三中所述的对2D断层图像合成方法中的相关参数进行调整,具体包括:光源点沿圆弧路径运动的半径调整为R.cos(θ),对扇束投影数据进行滤波处理的区间下限ul调整为R.sin(θ)-D.tan(γm),上限uu调整为R.sin(θ)+D.tan(γm),上述各式中,γm为扇束的扇形角,D为扇束光源点到其成像直线之间的距离,θ为扇束光源点和C臂旋转中心之连线与扇束中轴线之间的夹角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210038779.7A CN103258336B (zh) | 2012-02-20 | 2012-02-20 | 基于扇束虚拟平移的脊柱横断面图像合成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210038779.7A CN103258336B (zh) | 2012-02-20 | 2012-02-20 | 基于扇束虚拟平移的脊柱横断面图像合成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103258336A CN103258336A (zh) | 2013-08-21 |
CN103258336B true CN103258336B (zh) | 2015-12-09 |
Family
ID=48962229
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210038779.7A Expired - Fee Related CN103258336B (zh) | 2012-02-20 | 2012-02-20 | 基于扇束虚拟平移的脊柱横断面图像合成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103258336B (zh) |
Citations (4)
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 |
CN1336811A (zh) * | 1997-10-10 | 2002-02-20 | 模拟技术公司 | 计算层析扫描目标探测 |
CN1642479A (zh) * | 2002-03-12 | 2005-07-20 | 爱克斯康特公司 | 电离辐射的以扫描为基础的探测中的曝光控制 |
CN101842052A (zh) * | 2007-07-19 | 2010-09-22 | 北卡罗来纳大学查珀尔希尔分校 | 固定x射线数字化乳房断层合成系统和相关方法 |
-
2012
- 2012-02-20 CN CN201210038779.7A patent/CN103258336B/zh not_active Expired - Fee Related
Patent Citations (4)
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 |
CN1336811A (zh) * | 1997-10-10 | 2002-02-20 | 模拟技术公司 | 计算层析扫描目标探测 |
CN1642479A (zh) * | 2002-03-12 | 2005-07-20 | 爱克斯康特公司 | 电离辐射的以扫描为基础的探测中的曝光控制 |
CN101842052A (zh) * | 2007-07-19 | 2010-09-22 | 北卡罗来纳大学查珀尔希尔分校 | 固定x射线数字化乳房断层合成系统和相关方法 |
Non-Patent Citations (3)
Title |
---|
C形臂X线投影图像3D建模及其应用;闫士举等;《中国医学物理学杂志》;20100930;第27卷(第5期);第2109-2112页 * |
C形臂X线投影图像3D建模及其应用;闫士举等;《中国医学物理学杂志》;20100930;第27卷(第5期);第2109-2112页第1.1小节,第2-3节 * |
基于C型臂的Tomosynthesis快速重建方法;张利军等;《北京航空航天大学学报》;20060930;第32卷(第9期);第1110-1116页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103258336A (zh) | 2013-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7271646B2 (ja) | X線コンピュータ断層撮影装置、スキャン計画設定支援装置、医用画像診断システム、制御方法及び制御プログラム | |
US6298110B1 (en) | Cone beam volume CT angiography imaging system and method | |
CN107595387A (zh) | 一种基于超声拓片技术的脊椎图像生成系统以及脊柱手术导航定位系统 | |
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 | |
CN108601630A (zh) | 在脊柱区域和邻近胸廓、骨盆或头部的区域的外科手术期间的术中图控导航装置 | |
McNamara et al. | A flexible multicamera visual‐tracking system for detecting and correcting motion‐induced artifacts in cardiac SPECT slices | |
CN102456227A (zh) | Ct图像重建方法及装置 | |
CN105640575A (zh) | 医学图像生成装置、方法和程序 | |
CN104323789A (zh) | 一种ct扫描图像重建方法及ct扫描仪 | |
CN104200500A (zh) | 一种心脏图像的重建方法及装置 | |
JP2013244211A (ja) | 医用画像処理装置、医用画像処理方法及び制御プログラム | |
CN103750851A (zh) | 一种移动式x射线诊断装置 | |
CN109360213A (zh) | 一种基于脊柱超声冠状面图像的自动化椎体识别方法 | |
CN105894508A (zh) | 一种医学图像自动定位质量的评估方法 | |
US8934604B2 (en) | Image display apparatus and X-ray diagnostic apparatus | |
CN103366397A (zh) | 基于c形臂2d投影图像的脊柱3d模型构建方法 | |
CN106446515A (zh) | 一种三维医学图像显示方法及装置 | |
Chen et al. | Automatic inference and measurement of 3D carpal bone kinematics from single view fluoroscopic sequences | |
CN103258336B (zh) | 基于扇束虚拟平移的脊柱横断面图像合成方法 | |
CN103258335B (zh) | 基于普通c臂机的脊柱横断面图像合成方法 | |
US20110109630A1 (en) | Preparation and Presentation of Patient-Individual Panoramic Visualizations | |
CN103295263B (zh) | 基于锥束虚拟旋转的脊柱3d建模方法 | |
CN105069823B (zh) | 基于非对称横向双边截断投影数据的扇束ct重建方法 | |
JP4357928B2 (ja) | 脳虚血診断支援装置及びx線コンピュータ断層撮影装置 | |
Chen et al. | Development of Automatic Assessment Framework for Spine Deformity using Freehand 3D Ultrasound Imaging System |
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: 20151209 Termination date: 20190220 |
|
CF01 | Termination of patent right due to non-payment of annual fee |