通过CT机轴承间隙判别CT机图像空间分辨率的方法技术

本发明专利技术涉及一种医疗器械技术领域，ＣＴ机由于机械重心位置的位移，而轴承的偏移方向和偏移量发生变化，通过偏移方向建立几何模型、进行数学计算，并通过计算机程序计算出旋转中心的偏移量，此偏移量主要通过ＣＴ机轴承的间隙ｄ来推出，采用相应的已知变量参数角速度ω、旋转周期Ｔ、重力加速度ｇ、向心加速度ａ，再通过几何模型分别求出中心坐标偏移量，将偏移量送入计算机产生图像模体生数据，实现图像的重建，再和标准模体比较，给ＣＴ机图像空间分辨率的判定和机器部件轴承的选择提供有力的依据。

【技术实现步骤摘要】

本专利技术涉及一种医疗器械
，特别运用于CT机以判别空间分辨率的方法。
技术介绍
螺旋CT机中使用的轴承由内外两个同心圆组成，中间嵌有滚珠，而滚珠并不是紧紧的贴在轴承的内外壁上，而是有一定的间隙，我们把它叫做轴承间隙用d来表示，正是由于这个间隙的存在导致轴承内外壁移动和挤压，使轴承向不同的方向偏移旋转中心，偏移的大小等于轴承间隙，此间隙是随机不断变化，在理论上，螺旋CT机的旋转中心和整个机器的重心应该是同一个点，而由于机械上的各种原因，这种要求是很难达到的，往往这两个中心是不能重合，按技术要求这两个中心偏移的距离不得超过10mm，当整个CT机旋转时，相当于机器绕旋转中心做圆周运动，一个物体要想实现圆周运动，必须有向心力的存在，此向心力是由重量和轴承内外壁在不同方向挤压产生的压力共同提供的，机器的重心旋转到不同的位置，轴承将向不同的方向偏移以提供足够的向心力。目前，CT机由于机械结构复杂，精度误差等影响，使得旋转中心与机器重心偏移，甚至有时超过了机器允许限度，其结果，使得CT图像空间分辨率降低，严重的影响了图像效果，给精确的诊断带来不便。
技术实现思路
为解决以上CT机图像空间分辨率的判定，纠正诊断结果产生的误差，本专利技术的目的提供一种通过CT机轴承间隙来判别CT机图像空间分辨率的方法，具体采用几何模型，利用已知变量参数，通过数学模型、计算程序以实现用数据判别CT机空间机分辨率的好坏，或者用以选择CT机轴承精度做有力参照。本专利技术技术方案是这样实现的CT机由于机械重心的位移，而轴承的偏移方向和偏移量发生变化，使CT机图像空间的分辨率降低，通过偏移方向建立几何模型、进行数学计算，并通过计算机控制程序就可求得CT机偏移的中心坐标，具体判别CT机图像空间分辨率方法由下列步骤组成步骤1首先给定已知变量参数 重心与旋转中心的距离r；机器旋转一周所需要的时间T；角速度ω；重量加速度g；向心加速度α；轴承间隙大小d；步骤2建立几何模型当机械的重心在y轴时，即和球管的起始方向相同。其算法的数学表达式如下设一个临时变量temp代表轴承偏移方向的物理量ω＝2*π/Tα＝w2rtemp2＝α2+g2-2*α*g*cos()β＝+arcsin(g*sin()/temp)。当知道了某个旋转角度下轴承的偏移方向，计算实际球管位置所采集到的物体实际的几何模型和计算方法，以求得旋转中心的偏移坐标，其计算参数矢量图(如图2所示)。计算方法如下 dx＝dcos(β+pi-pi/2) (3)dy＝dsin(β+pi-pi/2) (4)centerx＝centerx+dx(5)centery＝centery+dy(6)上式中的(3)，(4)，(5)，(6)步可以用平移和旋转公式直接得到centerxT＝centerx*cosα+(centery+d)*sinαcenteryT＝centerx*sinα+(centery+d)*cosα。其中，图中中间的虚线圆是随着机器的旋转，重心旋转的轨迹，在球管旋转过任何角度时，都可以求得轴承偏移方向和y轴所成的角度β；绕坐标原点旋转α角度；配重码盘半径R；centerx、centery分别为原图像的中心坐标；dx、dy是物体的平移分量；centerxT、centeryT分别为经过平移和旋转过的中心坐标。同时，可建立另一个几何模型当重心不在y轴，也就是它和球管的旋转角度不是在同一角度下而是成一定的角度，当设定它和球管的旋转角度成90度的时候，即机器重心在x轴上(如图3所示)。从图中可以看出它和第一种模型的差别就是偏移的角度差了90度，故可如下算出β。其算法的数学表达式如下设一个临时变量temp代表轴承偏移方向的物理量ω＝2*π/Tα＝w2rtemp2＝α2+g2-2*α*g*cos(+pi/2)β＝+pi/2+arcsin(g*sin(+pi/2)/temp)在该模型下，由于球管旋转到角度时，轴承的偏移方向β从上面的模型中可以计算出来，相当于机器的旋转中心沿偏移方向移动轴承间隙d的大小，但是如果在运算的过程中，对球管和监测器进行几何变换会很麻烦，故我们对物体进行相应的逆变换可以得到相同的效果。求得轴承偏移的角度后仍然可以用上面模型的方法求得射线实际采得的投影数据。整个程序的模拟过程是先定义不同半径的圆，然后模拟上述的模型进行正投影得到的数据再进行图像重建，观察其空间分辨率的变化(如图3所示)。其程序计算方法由下列步骤组成。步骤1定义各圆的半径、中心位置和CT值(101)；步骤2判断所使用的几何模型(102)；步骤3轴承向固定方向偏移(103)或轴承偏移与球管位置的夹角固定(104)；步骤4计算向心力和重量的合力(105)，计算向心力和合力的夹角(106)； 计算合力与y轴正方向的夹角(轴承的偏移方向)(107)；步骤5对各个圆位置根据轴承的偏移方向进行相应的移动(108)；步骤6将各个圆沿轴承偏移方向进行平移(109)；步骤7再对它们进行相应角度的旋转(110)；步骤8按理论的投影方式对各个圆进行正投影得到生数据(111)；步骤9对生数据进行滤波反投影建像(112)；步骤10和已知标准模体作对比，做出CT机空间分辨率高低的判断。附图说明图1为本专利技术机械重心在y轴时的计算参数矢量图；图2为本专利技术对物体的几何变换过程计算参数矢量图；图3为本专利技术机械重心在x轴时的计算参数矢量图；图4为本专利技术计算方法流程图；图5为本专利技术计算模型一般差值CT机分辨率模体样图；图6为本专利技术计算模型1/4差值CT机分辨率模体样图；图7为本专利技术计算模型一般差值CT机分辨率模体样图；图8为本专利技术计算模型1/4差值CT机分辨率模体样图。图中1轴承内径，2轴承外径，3配重码盘半径，4轴承偏移方向，5轴承偏移方向与y轴的角度β，6球管旋转过的角度，7重力方向，8向心力方向，9机器重心，10实际球管位置，11理论球管位置，12旋转的角度α。本专利技术的优点通过CT机轴承移动间隙和已知变量参数，用数学模型的推导和计算得出中心坐标偏移量，将偏移量送入计算机就可得出图像的重建再和标准模体比较，就确定CT机空间分辨率的高低，同时和选择CT机轴承精度提供有力的参考。具体实施例方式本专利技术一种通过CT机轴承间隙判别CT机图像空间分辨率的方法，结合实施例加以说明。第一步在程序执行界面的option-setting fov对话框中输入建像视野(fov)、偏移角度(angle)、旋转速度(speed)、轴承间隙(offset)参数；第二步点击程序执行界面上的option-Modelspace产生图像模体生数据； 第三步点击程序执行界面上的option-Reconstruct实现图像的重建；第四步观察所建立的图像看其分辨率能达到多少。不断的改变步骤一中的offset和模型方式的参数，看其图像效果，其图像对比参看附图，从这些测试数据中可以为机械上轴承的设计或选择，提供一定的依据。其中图5为计算模型1，一般插值，左轴承间隙0.1，右轴承间隙1CT机分辨率的模体图样；图6为计算模型1，四分之一插值，左轴承间隙0.1，右轴承间隙1CT机分辨率的模体图样；图7为计算模型1，一般插值，左轴承间隙0.1，右轴承间隙1CT机分辨率的模体图样；本文档来自技高网...

【技术保护点】
一种通过ＣＴ机轴承间隙判别ＣＴ机图像空间分辨率的方法，ＣＴ机由于机械重心位移，而轴承的偏移方向和偏移量发生变化，使ＣＴ机图像空间分辨率降低，通过偏移方向建立几何模型、进行数学计算，并通过计算机控制程序就可求得ＣＴ机偏移中心坐标，其特征在于判别ＣＴ机图像空间分辩率方法由下列步骤组成：　　　　步骤１　　　　首先给定已知变量参数：　　　　重心与旋转中心的距离ｒ；　　　　机器旋转一周所需要的时间Ｔ；　　　　角速度ω；　　　　重量加速度ｇ；　　　　向心加速度α；　　　　轴承间隙大小ｄ；　　　　步骤２　　　　建立几何模型：当机械的重心在ｙ轴时，即和球管的起始方向相同；　　　　其算法的数学表达式如下：设一个临时变量ｔｅｍｐ代表轴承偏移方向的物理量　　　　ω＝２＊π／Ｔ　　　　α＝ｗ↑［２］ｒ　　　　ｔｅｍｐ↑［２］＝α↑［２］＋ｇ↑［２］－２＊α＊ｇ＊ｃｏｓ（φ）　　　　β＝φ＋ａｒｃｓｉｎ（ｇ＊ｓｉｎ（φ）／ｔｅｍｐ）　　　　当知道了某个旋转角度下轴承的偏移方向，计算实际球管位置所采集到的物体实际的几何模型和计算方法，以求得旋转中心的偏移坐标，计算方法如下：　　　　ｘ＝＊　　（１）　　　　Ｒ／ｓｉｎα＝ｘ／ｓｉｎ（π－φ）＊α＝ａｒｃｓｉｎ（ｄ　　ｓｉｎ（β－φ）／ｘ）　　（２）　　　　ｄｘ＝ｄ　　ｃｏｓ（β＋ｐｉ－ｐｉ／２）　　（３）　　　　ｄｙ＝ｄ　　ｓｉｎ（β＋ｐｉ－ｐｉ／２）　　（４）　　　　ｃｅｎｔｅｒｘ＝ｃｅｎｔｅｒｘ＋ｄｘ　　（５）　　　　ｃｅｎｔｅｒｙ＝ｃｅｎｔｅｒｙ＋ｄｙ　　（６）　　　　上式中的（３），（４），（５），（６）步可以用平移和旋转公式直接得到：　　　　ｃｅｎｔｅｒｘＴ＝ｃｅｎｔｅｒｘ＊ｃｏｓ　　α＋（ｃｅｎｔｅｒｙ＋ｄ）＊ｓｉｎ　　α　　　　ｃｅｎｔｅｒｙＴ＝ｃｅｎｔｅｒｘ＊ｓｉｎ　　α＋（ｃｅｎｔｅｒｙ＋ｄ）＊ｃｏｓ　　α　　　　ｃｅｎｔｅｒｘＴ、ｃｅｎｔｅｒｙＴ分别为经过平移和旋转过的中心坐标；　　　　整个程序的模拟过程是先定义不同半径的圆，然后模拟上述的模型进行正投影得到的数据再进行图像重建，观察其空间分辨率的变化。...

【技术特征摘要】

【专利技术属性】
技术研发人员：李双学，楼珊珊，刘晋军，佟丽霞，
申请(专利权)人：东软飞利浦医疗设备系统有限责任公司，
类型：发明
国别省市：89[中国|沈阳]