CN103340644B - 一种图像位置补偿方法及装置 - Google Patents
一种图像位置补偿方法及装置 Download PDFInfo
- Publication number
- CN103340644B CN103340644B CN201310273523.9A CN201310273523A CN103340644B CN 103340644 B CN103340644 B CN 103340644B CN 201310273523 A CN201310273523 A CN 201310273523A CN 103340644 B CN103340644 B CN 103340644B
- Authority
- CN
- China
- Prior art keywords
- image
- scanning bed
- height
- bed
- width
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 87
- 230000000007 visual effect Effects 0.000 claims description 66
- 238000007665 sagging Methods 0.000 claims description 61
- 239000011159 matrix material Substances 0.000 claims description 23
- 238000012545 processing Methods 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims description 10
- 239000000284 extract Substances 0.000 claims description 4
- 238000003745 diagnosis Methods 0.000 abstract description 8
- 239000000463 material Substances 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 238000007689 inspection Methods 0.000 description 8
- 238000012360 testing method Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000004064 recycling Methods 0.000 description 2
- 230000003187 abdominal effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008571 general function Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
Abstract
本发明实施例公开了一种图像位置补偿方法及装置,其中方法具体是测量扫描床在各个z位置的实际高度以及对应的扫描床下垂偏移量;根据所述扫描床下垂偏移量,在后重建过程中,调整每一幅图像的建像中心纵坐标,得到调整后的每一幅重建图像;从所述调整后的每一幅重建图像中获取的每一行的坐标信息和对应的像素值,采用插值的方式得到目标图像每一行的像素值。目标图像就是位置补偿后的图像,通过本发明能够实现对现有的CT设备生成的图像进行补偿,以得到正确的断层图像,为医疗诊断提供更精确的图像信息。
Description
技术领域
本发明涉及医学图像技术领域,特别是涉及一种图像位置补偿方法及装置。
背景技术
X射线计算机断层扫描(Computed Tomography,CT)是利用计算机技术对被测物体断层扫描图像进行重建获得三维断层图像的扫描方式。该扫描方式是通过单一轴面的射线穿透被测物体,根据被测物体各部分对射线的吸收和透过率不同,又计算机采集透过射线并通过三维重构成像。该技术可应用于很多领域,比如:医学检测、工业检测、安保检测、航空运输、运输港湾、大型货物集装箱案件装置等等。
以医学检测为例:CT设备的工作原理具体是:根据人体不同组织对X射线的吸收和透过率的不同,应用灵敏度极高的仪器对人体进行测量,然后将测量所获取的数据输入电子计算机,电子计算机对数据进行处理后,就可摄下人体被检测部位的断面或者立体的图像,发现体内任何部位的细小病变。
在实际检测中,受检者被放置在扫描床上,扫描床从支撑工作台中逐渐伸出,将受检者传送至由射线源和检测器构成的扫描平面内,射线源和检测器围绕受检者旋转,获得受检者的扫描数据。但是由于受到受检者重力的影响,伸出的扫描床伸出的长度不同,受到的重力不同将发生一定的下垂偏移,导致标称的扫描床位置与实际不符。当扫描床位置与实际不符时,CT扫描装置后续重建出的受检者断层图像位置和形态也会与发生偏移,与预期不相符合,进而给医疗诊断造成一定干扰。
基于上述技术问题,目前迫切需要提供一种图像位置补偿方法,以实现使用X射线CT设备的普通功能执行扫描床位置检测方法,通过图像位置补偿方法生成的扫描图像进行补偿得到正确的断层图像,以为医疗诊断提供更精确的图像信息。
发明内容
为了解决上述技术问题,本发明实施例中提供了一种图像位置补偿方法及装置,以实现对现有的CT设备生成的图像进行补偿,以得到正确的断层图像,为医疗诊断提供更精确的图像信息。
本申请实施例公开了如下技术方案:
一种图像位置补偿方法,包括:
测量扫描床在各个z位置的实际高度以及对应的扫描床下垂偏移量;
根据所述扫描床下垂偏移量,在后重建过程中,调整每一幅图像的建像中心纵坐标,得到调整后的每一幅重建图像;
从所述调整后的每一幅重建图像中获取的每一行的坐标信息和对应的像素值,采用插值的方式得到目标图像每一行的像素值。
优选的,所述测量扫描床在z位置的实际高度以及对应的扫描床下垂偏移量,包括:
获取机架的旋转中心高度、每一幅预览图像的图像中心坐标、图像重建视野和图像矩阵大小;
依次搜索每一幅预览图像中像素值大于阈值的行序号,确定出每一预览幅图像中行序号最大的行;
根据预览图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、图像重建视野和所述确定出的行序号最大的行,计算扫描床在每一幅预览图像建立时的实际高度;
依次计算扫描床在每一幅预览图像建立时的实际高度与扫描床的标称床高之差,所述差值作为对应的扫描床下垂偏移量。
优选的,在获取机架的旋转中心高度、预览图像中心坐标和预览图像重建视野和预览图像重建矩阵大小,之后还包括:
根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、预览图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在每一幅图像中的理论像素纵坐标值;
则所述依次搜索每一幅预览图像中像素值大于阈值的行数,确定出行序号最大的行,具体为:
依次根据计算出的扫描床在每一幅预览图像中的理论像素纵坐标值,在行序号为纵坐标值的上下任意多行进行搜索,搜索出存在像素值大于阈值的行数,确定出行序号最大的行。
优选的,所述根据所述扫描床标称床高、预览图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、预览图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在每一幅预览图像中的像素纵坐标值,包括:
获取的所述扫描床标称床高H2、预览图像矩阵大小M×M、机架的旋转中心高度H1、每一幅预览图像的图像中心坐标(cXi,cYi)、图像重建视野FOV;
按照公式:依次计算当扫描床的高度为标称床高时,所述扫描床在每一幅预览图像中的像素纵坐标值其中,公式表示取整运算,表示每个像素所代表的实际宽度,cYi表示每一幅预览图像的图像中心坐标的纵坐标值。
优选的,根据预览图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、预览图像重建视野和所述确定出行序号最大的行,计算扫描床在每一幅预览图像建立时的实际高度,包括:
获取的所述扫描床标称床高H2、预览图像矩阵大小M×M、机架的旋转中心高度H1、每一幅预览图像的图像中心坐标(cXi,cYi)、预览图像重建视野FOV和所述确定出行数目最大的行数Ycouch;
按照公式: 依次计算扫描床在每一幅预览图像建立时的实际高度Hi,其中,为每个像素所代表的实际宽度。
优选的,所述测量扫描床在z位置的实际高度以及对应的扫描床下垂偏移量,包括:
获取机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和图像矩阵大小;
依次搜索平片图像中每一列中的行像素值大于阈值的行序号,确定出每一列中行序号最大的行;
根据平片图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和所述确定出的行序号最大的行,计算扫描床在每一个z位置的实际高度;
依次计算扫描床在每一个z位置的实际高度与扫描床的标称床高之差,所述差值作为对应的扫描床下垂偏移量。
优选的,获取机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和图像矩阵大小,之后还包括:
根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、平片图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在平片图像中的理论像素纵坐标值;
则所述依次搜索平片图像中每一列中的行像素值大于阈值的行序号,确定出每一列中行序号最大的行,具体为:
依次根据计算出的扫描床在平片图像中的理论像素纵坐标值,在行序号为纵坐标值的上下任意多行进行搜索,搜索出存在行像素值大于阈值的行序号,确定出每一列中行序号最大的行。
优选的,所述根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、平片图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在平片图像中的理论像素纵坐标值,包括:
获取的所述扫描床标称床高H2、平片图像矩阵大小N×M、机架的旋转中心高度H1、平片图像的图像中心坐标(cX,cY)、图像重建视野FOV;
按照公式:计算当扫描床的高度为标称床高时,所述扫描床在平片图像中的像素纵坐标值其中,公式表示取整运算,表示每个像素所代表的实际宽度,cY表示平片图像的图像中心坐标的纵坐标值。
优选的,所述根据平片图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和所述确定出的行序号最大的行,计算扫描床在每一个z位置的实际高度,包括:
获取的所述扫描床标称床高H2、平片图像矩阵大小N×M、机架的旋转中心高度H1、平片图像的图像中心坐标(cX,cY)、平片图像重建视野FOV和所述确定出行数目最大的行数Ycouch;
按照公式: 依次计算扫描床在每个z位置Zi处的实际高度Hi,其中,为每个像素所代表的实际高度。
优选的,从所述调整后的每一幅重建图像中获取的每一行的坐标信息和对应的像素值,采用插值的方式得到目标图像每一行的像素值,包括:
读取所述调整后的每一幅重建图像中每一行的坐标(Y(m),Z(m))为:
Z(m)=ZP,其中,m=0,…,M-1表示行数目,ZP表示扫描床在建立第P图像时对应的刻度位置,C(ZP)表示第P图像的调整后图像中心纵坐标;
计算所述当前目标图像每一行的坐标信息(y(n),z(n)),其中,y(n)和z(n)分别代表第n行的y轴和z轴坐标值,n=0,…,M-1表示行数目;
对读取的(Y(m),Z(m))做两个方向的插值处理,得到目标图像坐标为(y(n),z(n))图像像素值。
优选的,所述当前目标图像所在直线通过以下方式得到,包括:
采用最小二乘法对扫描床在每一幅图像建立时所在的刻度位置Zi以及位置Zi对应的下垂偏移量D(Zi)进行曲线拟合,拟合出扫描床偏移所在直线y=kz+b;
根据当前目标图像所在直线与扫描床偏移所在直线的垂直关系,计算出目标图像所在的直线 其中 为扫描床偏移所在直线在z轴上的截距;
按照所述当前目标图像所在的直线 计算出目标图像每一行的坐标信息(y(n),z(n))为:
其中,y(n)和z(n)分别代表第n行的y轴和z轴坐标值,n=0,…,M-1表示行数目。
本发明还提供了一种图像位置补偿装置,包括:
测量单元,用于测量扫描床在各个z位置的实际高度以及对应的扫描床下垂偏移量;
调整单元,用于根据所述扫描床下垂偏移量,在后重建过程中,调整每一幅图像的建像中心纵坐标,得到调整后的每一幅重建图像;
插值单元,用于从所述调整后的每一幅重建图像中获取的每一行的坐标信息和对应的像素值,采用插值的方式得到目标图像每一行的像素值。
由上述实施例可以看出,本发明提供的一种图像位置补偿方法和装置,通过测量扫描床在各个z位置的实际高度以及对应的扫描床下垂偏移量;根据所述扫描床下垂偏移量,在后重建过程中,调整每一幅图像的建像中心纵坐标,得到调整后的每一幅重建图像;从所述调整后的每一幅重建图像中获取的每一行的坐标信息和对应的像素值,采用插值的方式得到目标图像每一行的像素值。由于,扫描床在工作过程中受到受检者重力的影响,在逐渐送入检测设备的过程中会发生不同程度的下垂偏移,导致扫描所得图像并不能真实的反应扫描床下垂偏移的信息,得到的图像反应的是物体斜切面的信息,而实际目标图像是正截面的信息,因此通过利用目标图像与上述重建图像之间的关系,对其进行插值处理最终得到目标图像。可见:通过本发明的图像补偿方法及装置,能够对原始图像进行调整得到正确的断层图像,为医疗诊断提供正确的图像信息。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例一揭示的一种图像位置补偿方法的方法流程图;
图2为本申请实施例二揭示的一种测量扫描床在z位置的实际高度以及对应的扫描床下垂偏移量的方法流程图;
图3为本申请实施例二揭示的另一种测量扫描床在z位置的实际高度以及对应的扫描床下垂偏移量的方法流程图;
图4为本申请实施例三揭示的一种测量扫描床在z位置的实际高度以及对应的扫描床下垂偏移量的方法流程图;
图5为本申请实施例三揭示的一种图像位置补偿装置的装置结构图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图对本发明实施例进行详细描述。
需要说明的是本发明的方法可以用于X射线CT设备中,一般的CT设备包括:扫描床、工作台、射线源和检测器,整个扫描过程为:首先被检测物放在扫描床上;其次,扫描床由工作台中伸出,将被检测物传送至由射线源和检测器构成的扫描平面内;最后,射线源和检测器绕着被检测物旋转,获得被检测物的相关扫描数据,一般情况需要两次建立图像的处理,第一次处理是根据相关扫描数据通过快速建像得到每一幅预览图像,这些图像主要是便于医生能够快速对病灶或者身体部位进行预览或者定位。第二次处理是根据扫描得到的相关扫描数据进行后重建处理,得到每一幅重建图像,这些图像就能够更为细致地反应被检测物的实际情况,将第二次处理过程成为后重建过程。
为了确定后续扫描操作应该扫描的位置范围,在扫描之前CT设备还需要做生成一张平片图像,该图像用于确定后续扫描的位置范围,比如:当前需要对被测试者的腹部进行检测,由于每一个测试者的身高不同、实际躺在扫描床上的位置也不一定完全一致,所以需要建立平片图像,确定后续扫描的z位置范围。
实施例一
请参阅图1,其为本申请实施例一揭示的一种图像补偿方法的方法流程图,该方法包括以下步骤:
S1:测量扫描床在各个z位置的实际高度以及对应的扫描床下垂偏移量;
由于在实际扫描过程中,扫描床是将被检测物传送至由射线源和检测器构成的扫描平面内,受到被检测物的压力,扫描床在不同的刻度位置会有相同或者不同程度的下垂,因此与扫描床的标称床高相比就会产生下垂偏移量。实际应用中任意一个CT设备中的扫描床都有标称床高,一般可以从设备的产品说明书中获取。
在CT设备中扫描床床板上是有刻度位置,也就说,每一个位置都对应着刻度值,比如:当扫描床长度为三米,扫描床的从头到尾依次对应的刻度值为0~300厘米,这些刻度经常用于指出扫描床在某个位置建立的图像,例如当前需要对被检测者进行扫描建立5张图片,在刻度为31~45厘米,每隔1厘米扫描建立一张图像,那扫描结果也就是在刻度位置31厘米、32厘米、33厘米、34厘米、35厘米分别都建立了一张图像,共建立5张图像。在建立这五张图像时,扫描床在对应的刻度位置都发生了下垂,因此可以计算出相对比扫描床标称床高的下垂偏移量。对于整个CT设备而言,射线源和检测器所在的平面看做为XY平面,扫描床所在位置通常情况下成为与XY平面垂直的Z方向,因此扫描床上的刻度位置称为Z位置。
S2:根据所述扫描床下垂偏移量,在后重建过程中,调整每一幅图像的建像中心纵坐标,得到调整后的每一幅重建图像;
在扫描建立图像时,在扫描床的什么刻度位置建立图像,建立多少张图像都可以根据实际需求预先设置,每一张图像都有图像中心坐标,每一张图像的图像中心坐标都是以设备的机架旋转中心为原点的平面坐标系以基准的。
由于扫描床下垂后,对于图像中心坐标而言其横坐标没有发生变化,只是纵坐标发生了变化。在扫描建立图像的确定的图像中心坐标是以扫描床的标称床高为基准设置的,当扫描床向下偏移时,实际该中心坐标也应该向下移动。根据在每一个刻度位置扫描建立图像时扫描床的下垂偏移量,调整图像的图像中心坐标的纵坐标,对应得到重建图像。即,原图像中所有的点都发生了下垂,得到重建图像。调整图像中心坐标的纵坐标具体是:当图像中心坐标的纵坐标为cY,调整后的纵坐标为C(Zi),对应的下垂偏移量为D(Zi),按照公式C(Zi)=cY-D(Zi)进行调整。
S3:从所述调整后的每一幅重建图像中获取的每一行的坐标信息和对应的像素值,采用插值的方式得到目标图像每一行的像素值。
根据与被补偿的图像相邻多张图像的坐标信息和对应的像素值,进行插值处理得到目标图像,所述目标图像就是真正能够反映被测试物的实际信息的图像,也就是被测试物的正截面信息。
通过上述实施例可以看出,本发明的图像位置补偿方法,通过测量扫描床在各个z位置的实际高度以及对应的扫描床下垂偏移量;根据所述扫描床下垂偏移量,在后重建过程中,调整每一幅图像的建像中心纵坐标,得到调整后的每一幅重建图像;从所述调整后的每一幅重建图像中获取的每一行的坐标信息和对应的像素值,采用插值的方式得到目标图像每一行的像素值,也就得到最终的目标图像。这样处理是通过研究分析目标图像与重建图像之间的关系,通过插值计算得到目标图像,这样经过位置补正后的图像是真实反映被检测物的信息的图像,能够为医疗诊断提供正确的图像信息,同时也就提高CT设备测试质量。
实施例二
下面以具体的应用场景对发明的图像位置补偿方法进行解释说明,具体为:CT设备需要扫描被检测者,得到受检者的扫描图像,最初设置需要在扫描床的10个不同的刻度位置,分别扫描建立1张图像,共建立10张图像,通过对预览图像的处理测量扫描床下垂偏移量,具体的方法包括以下步骤:
步骤201:测量扫描床在各个z位置的实际高度以及对应的扫描床下垂偏移量;
每一幅预览图像在最初扫描建立时都对应着扫描床的一个具体的z位置。
优选的,所述测量扫描床在z位置的实际高度以及对应的扫描床下垂偏移量,具体如图2的一种测量扫描床在z位置的实际高度以及对应的扫描床下垂偏移量的方法流程图所示,具体包括步骤A1~A4:
A1:获取机架的旋转中心高度、每一幅预览图像的图像中心坐标、图像重建视野和图像矩阵大小;
A2:依次搜索每一幅预览图像中像素值大于阈值的行序号,确定出每一预览幅图像中行序号最大的行;
A3:根据预览图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、图像重建视野和所述确定出的行序号最大的行,计算扫描床在每一幅预览图像建立时的实际高度;
A4:依次计算扫描床在每一幅预览图像建立时的实际高度与扫描床的标称床高之差,所述差值作为对应的扫描床下垂偏移量。
为了提高测量效率,优选的,上述测量扫描床在每一幅图像建立时的实际高度的方法中的步骤A1之后,还包括:
根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、预览图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在每一幅图像中的理论像素纵坐标值;
则所述依次搜索每一幅预览图像中像素值大于阈值的行数,确定出行序号最大的行,具体为:
依次根据计算出的扫描床在每一幅预览图像中的理论像素纵坐标值,在行序号为纵坐标值的上下任意多行进行搜索,搜索出存在像素值大于阈值的行数,确定出行序号最大的行。具体方法如图3的另一种测量扫描床在z位置的实际高度以及对应的扫描床下垂偏移量的方法流程图所示。
优选的,步骤B2:所述根据所述扫描床标称床高、预览图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、预览图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在每一幅预览图像中的像素纵坐标值,包括:
获取的所述扫描床标称床高H2、预览图像矩阵大小M×M、机架的旋转中心高度H1、每一幅预览图像的图像中心坐标(cXi,cYi)、图像重建视野FOV;
按照公式:依次计算当扫描床的高度为标称床高时,所述扫描床在每一幅预览图像中的像素纵坐标值其中,公式表示取整运算,表示每个像素所代表的实际宽度,cYi表示每一幅预览图像的图像中心坐标的纵坐标值。
优选的,根据预览图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、预览图像重建视野和所述确定出行序号最大的行,计算扫描床在每一幅预览图像建立时的实际高度,包括:
获取的所述扫描床标称床高H2、预览图像矩阵大小M×M、机架的旋转中心高度H1、每一幅预览图像的图像中心坐标(cXi,cYi)、预览图像重建视野FOV和所述确定出行数目最大的行数Ycouch;
按照公式: 依次计算扫描床在每一幅预览图像建立时的实际高度Hi,其中,为每个像素所代表的实际宽度。
步骤202:根据所述扫描床下垂偏移量,在后重建过程中,调整每一幅图像的建像中心纵坐标,得到调整后的每一幅重建图像;
步骤203:从所述调整后的每一幅重建图像中获取的每一行的坐标信息和对应的像素值,采用插值的方式得到目标图像每一行的像素值。
优选的步骤203主要通过以下方式实现,具体包括:
读取所述调整后的每一幅重建图像中每一行的坐标(Y(m),Z(m))为:
Z(m)=ZP,其中,m=0,…,M-1表示行数目,ZP表示扫描床在建立第P图像时对应的刻度位置,C(ZP)表示第P图像的调整后图像中心纵坐标;
计算所述当前目标图像每一行的坐标信息(y(n),z(n)),其中,y(n)和z(n)分别代表第n行的y轴和z轴坐标值,n=0,…,M-1表示行数目;
对读取的(Y(m),Z(m))做两个方向的插值处理,得到目标图像坐标为(y(n),z(n))图像像素值。
优选的,所述当前目标图像所在直线通过以下方式得到,包括:
采用最小二乘法对扫描床在每一幅图像建立时所在的刻度位置Zi以及位置Zi对应的下垂偏移量D(Zi)进行曲线拟合,拟合出扫描床偏移所在直线y=kz+b;
曲线拟合是指选择适当的曲线类型拟合观测数据,并用拟合的曲线方程分析两个变量间的关系。在本发明中利用曲线拟合算法,观测数据具体是指扫描床在每一幅图片建立时所在的刻度位置Zi和扫描床扫描时实际高度,即,曲线方程中的两个变量分别是指扫描建立图片时所在的刻度位置值和实际高度值
根据当前目标图像所在直线与扫描床偏移所在直线的垂直关系,计算出目标图像所在的直线 其中 为扫描床偏移所在直线在z轴上的截距;
按照所述当前目标图像所在的直线 计算出目标图像每一行的坐标信息(y(n),z(n))为:
其中,y(n)和z(n)分别代表第n行的y轴和z轴坐标值,n=0,…,M-1表示行数目。
通过上述实施例可以看出,本发明提供的测量扫描床的方法,可以方便准确地测试出扫描床在各个建立图像位置的实际高度,利用测试的实际高度计算得到的下垂偏移量对每一幅重建图像进行调整,再利用目标图像是正截面与扫描床之间是垂直关系,对所述调整后的重建图像中坐标数据进行插值处理,得到目标图像各个坐标对应的像素值,也就得到最终的目标图像。这样处理使得在后重建过程中就将图像进行补偿,使得处理过程更加方便,这样经过位置补正后的图像是真实反映被检测物的信息的图像,可以校正扫描床下垂带来的图像中心偏移以及图像位置形态不准确性,能够为医疗诊断提供正确的图像信息,同时也就提高CT设备测试质量。
实施例三
本发明实施例三提供了一种通过处理平片图像得到扫描床实际高度以及对应的下垂偏移量的方法,如图4的一种测量扫描床在z位置的实际高度以及对应的扫描床下垂偏移量的方法流程图所示,具体包括以下步骤:
步骤301:获取机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和图像矩阵大小;
步骤302:依次搜索平片图像中每一列中的行像素值大于阈值的行序号,确定出每一列中行序号最大的行;
步骤303:根据平片图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和所述确定出的行序号最大的行,计算扫描床在每一个z位置的实际高度;
优选的,步骤303具体包括:
获取的所述扫描床标称床高H2、平片图像矩阵大小N×M、机架的旋转中心高度H1、平片图像的图像中心坐标(cX,cY)、平片图像重建视野FOV和所述确定出行数目最大的行数Ycouch;
按照公式: 依次计算扫描床在每个z位置Zi处的实际高度Hi,其中,为每个像素所代表的实际高度。
步骤304:依次计算扫描床在每一个z位置的实际高度与扫描床的标称床高之差,所述差值作为对应的扫描床下垂偏移量。
为了提高测量速率,本发明还提供了另外一种测量方法,具体是在上述步骤301:获取机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和图像矩阵大小,之后还包括:
根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、平片图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在平片图像中的理论像素纵坐标值;
则所述依次搜索平片图像中每一列中的行像素值大于阈值的行序号,确定出每一列中行序号最大的行,具体为:
依次根据计算出的扫描床在平片图像中的理论像素纵坐标值,在行序号为纵坐标值的上下任意多行进行搜索,搜索出存在行像素值大于阈值的行序号,确定出每一列中行序号最大的行。
优选的,所述根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、平片图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在平片图像中的理论像素纵坐标值,包括:
获取的所述扫描床标称床高H2、平片图像矩阵大小N×M、机架的旋转中心高度H1、平片图像的图像中心坐标(cX,cY)、图像重建视野FOV;
按照公式:计算当扫描床的高度为标称床高时,所述扫描床在平片图像中的像素纵坐标值其中,公式表示取整运算,表示每个像素所代表的实际宽度,cY表示平片图像的图像中心坐标的纵坐标值。
实施例四
与上述实施例一中的一种图像位置补偿方法相对应,本申请实施例提供了一种图像位置补偿装置。请参阅图5,其为本申请实施例三揭示的一种图像位置补偿装置的装置结构图,该装置包括:测量单元401、调整单元402和插值单元403。下面结合该装置的工作原理进一步介绍其内部结构及其连接关系。
测量单元401,用于测量扫描床在各个z位置的实际高度以及对应的扫描床下垂偏移量;
调整单元402,用于根据所述扫描床下垂偏移量,在后重建过程中,调整每一幅图像的建像中心纵坐标,得到调整后的每一幅重建图像;
插值单元403,用于从所述调整后的每一幅重建图像中获取的每一行的坐标信息和对应的像素值,采用插值的方式得到目标图像每一行的像素值。
优选的,所述测量单元,包括:
第一获取模块,用于获取机架的旋转中心高度、每一幅预览图像的图像中心坐标、图像重建视野和图像矩阵大小;
第一搜索模块,用于依次搜索每一幅预览图像中像素值大于阈值的行序号,确定出每一预览幅图像中行序号最大的行;
第一计算模块,用于根据预览图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、图像重建视野和所述确定出的行序号最大的行,计算扫描床在每一幅预览图像建立时的实际高度;
第一差值模块,用于依次计算扫描床在每一幅预览图像建立时的实际高度与扫描床的标称床高之差,所述差值作为对应的扫描床下垂偏移量。
优选的,所述测量单元,还包括:
第一估算模块,用于根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、预览图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在每一幅图像中的理论像素纵坐标值;
则所述第一搜索模块,具体为:
依次根据计算出的扫描床在每一幅预览图像中的理论像素纵坐标值,在行序号为纵坐标值的上下任意多行进行搜索,搜索出存在像素值大于阈值的行数,确定出行序号最大的行。
优选的,所述第一估算模块,包括:
第一获取子模块,用于获取的所述扫描床标称床高H2、预览图像矩阵大小M×M、机架的旋转中心高度H1、每一幅预览图像的图像中心坐标(cXi,cYi)、图像重建视野FOV;
第一估算子模块,用于按照公式:依次计算当扫描床的高度为标称床高时,所述扫描床在每一幅预览图像中的理论像素纵坐标值其中,公式表示取整运算,表示每个像素所代表的实际宽度,cYi表示每一幅预览图像的图像中心坐标的纵坐标值。
优选的,所述第一计算模块,包括:
第一提取子模块,用于获取的所述扫描床标称床高H2、预览图像矩阵大小M×M、机架的旋转中心高度H1、每一幅预览图像的图像中心坐标(cXi,cYi)、预览图像重建视野FOV和所述确定出行数目最大的行数Ycouch;
第一计算子模块,用于按照公式: 依次计算扫描床在每一幅预览图像建立时的实际高度Hi,其中,为每个像素所代表的实际宽度。
所述测量单元,包括:
第二获取模块,用于获取机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和图像矩阵大小;
第二搜索模块,用于依次搜索平片图像中每一列中的行像素值大于阈值的行序号,确定出每一列中行序号最大的行;
第二计算模块,用于根据平片图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和所述确定出的行序号最大的行,计算扫描床在每一个z位置的实际高度;
第二差值模块,用于依次计算扫描床在每一个z位置的实际高度与扫描床的标称床高之差,所述差值作为对应的扫描床下垂偏移量。
为了提高装置的测量单元的测量精确度,优选的,所述测量单元,还包括:
第二估算模块,用于根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、平片图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在平片图像中的理论像素纵坐标值;
则所述第二搜索模块,具体为:
依次根据计算出的扫描床在平片图像中的理论像素纵坐标值,在行序号为纵坐标值的上下任意多行进行搜索,搜索出存在行像素值大于阈值的行序号,确定出每一列中行序号最大的行。
优选的,所述第二估算模块,包括:
第二获取子模块,用于获取的所述扫描床标称床高H2、平片图像矩阵大小N×M、机架的旋转中心高度H1、平片图像的图像中心坐标(cX,cY)、平片图像重建视野FOV;
第二估算子模块,用于按照公式:计算当扫描床的高度为标称床高时,所述扫描床在平片图像中的像素纵坐标值其中,公式表示取整运算,表示每个像素所代表的实际宽度,cY表示平片图像的图像中心坐标的纵坐标值。
优选的,所述第二计算模块,包括:
第二提取子模块,用于获取的所述扫描床标称床高H2、平片图像矩阵大小N×M、机架的旋转中心高度H1、平片图像的图像中心坐标(cX,cY)、平片图像重建视野FOV和所述确定出行数目最大的行数Ycouch;
第二计算子模块,用于按照公式: 依次计算扫描床在每个z位置Zi处的实际高度Hi,其中,为每个像素所代表的实际高度。
优选的,所述插值单元,包括:
第一坐标模块,用于读取所述调整后的每一幅重建图像中每一行的坐标(Y(m),Z(m))为:
Z(m)=ZP
,其中,m=0,…,M-1表示行数目,ZP表示扫描床在建立第P图像时对应的刻度位置,C(ZP)表示第P图像的调整后图像中心纵坐标;
第二坐标模块,用于计算所述当前目标图像每一行的坐标信息(y(n),z(n)),其中,y(n)和z(n)分别代表第n行的y轴和z轴坐标值,n=0,…,M-1表示行数目;
第一处理模块,用于对读取的(Y(m),Z(m))做两个方向的插值处理,得到目标图像坐标为(y(n),z(n))图像像素值。
优选的,所述第二坐标模块,包括:
拟合子模块,用于采用最小二乘法对扫描床在每一幅图像建立时所在的刻度位置Zi以及位置Zi对应的下垂偏移量D(Zi)进行曲线拟合,拟合出扫描床偏移所在直线y=kz+b;
直线子模块,用于根据当前目标图像所在直线与扫描床偏移所在直线的垂直关系,计算出目标图像所在的直线 其中为扫描床偏移所在直线在z轴上的截距;
坐标子模块,用于按照所述当前目标图像所在的直线
其中,y(n)和z(n)分别代表第n行的y轴和z轴坐标值,n=0,…,M-1表示行数目。
本申请中的上述装置中的各个模块,可以集成在CT扫描设备中,也可以是以独立的设备存在。
通过上述实施例可以看出,本发明提供的图像位置补偿装置,通过测量模块可以方便准确地测试出扫描床在各个建立图像位置的实际高度,利用测试的实际高度计算得到的下垂偏移量对每一幅重建图像进行调整,再利用目标图像是正截面与扫描床之间是垂直关系,对所述调整后的重建图像中坐标数据进行插值处理,得到目标图像各个坐标对应的像素值,也就得到最终的目标图像。这样处理使得在后重建过程中就将图像进行补偿,使得处理过程更加方便,这样经过位置补正后的图像是真实反映被检测物的信息的图像,可以校正扫描床下垂带来的图像中心偏移以及图像位置形态不准确性,能够为医疗诊断提供正确的图像信息,同时也就提高CT设备测试质量。
需要说明的是,本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统或装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
还需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个......”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (20)
1.一种图像位置补偿方法,其特征在于,包括:
测量扫描床在各个z位置的实际高度以及对应的扫描床下垂偏移量;
根据所述扫描床下垂偏移量,在后重建过程中,调整每一幅图像的建像中心纵坐标,得到调整后的每一幅重建图像;
从所述调整后的每一幅重建图像中获取的每一行的坐标信息和对应的像素值,采用插值的方式得到目标图像每一行的像素值,包括:
读取所述调整后的每一幅重建图像中每一行的坐标(Y(m),Z(m))为:
Z(m)=ZP,其中,m=0,…,M-1表示行序号,ZP表示扫描床在建立第P图像时对应的刻度位置,C(ZP)表示第P图像的调整后图像中心纵坐标;表示每个像素所代表的实际宽度;
计算当前目标图像每一行的坐标信息(y(n),z(n)),其中,y(n)和z(n)分别代表第n行的y轴和z轴坐标值,n=0,…,M-1表示行序号;
对读取的(Y(m),Z(m))做两个方向的插值处理,得到目标图像坐标为(y(n),z(n))图像像素值。
2.根据权利要求1所述的方法,其特征在于,所述测量扫描床在z位置的实际高度以及对应的扫描床下垂偏移量,包括:
获取机架的旋转中心高度、每一幅预览图像的图像中心坐标、图像重建视野和图像矩阵大小;
依次搜索每一幅预览图像中像素值大于阈值的行序号,确定出每一预览幅图像中行序号最大的行;
根据预览图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、图像重建视野和所述确定出的行序号最大的行,计算扫描床在每一幅预览图像建立时的实际高度;
依次计算扫描床在每一幅预览图像建立时的实际高度与扫描床的标称床高之差,所述差值作为对应的扫描床下垂偏移量。
3.根据权利要求2所述的方法,其特征在于,在获取机架的旋转中心高度、预览图像中心坐标和预览图像重建视野和预览图像重建矩阵大小,之后还包括:
根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、预览图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在每一幅图像中的理论像素纵坐标值;
则所述依次搜索每一幅预览图像中像素值大于阈值的行数,确定出行序号最大的行,具体为:
依次根据计算出的扫描床在每一幅预览图像中的理论像素纵坐标值,在行序号为纵坐标值的上下任意多行进行搜索,搜索出存在像素值大于阈值的行数,确定出行序号最大的行。
4.根据权利要求3所述的方法,其特征在于,所述根据所述扫描床标称床高、预览图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、预览图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在每一幅预览图像中的像素纵坐标值,包括:
获取的所述扫描床标称床高H2、预览图像矩阵大小M×M、机架的旋转中心高度H1、每一幅预览图像的图像中心坐标(cXi,cYi)、图像重建视野FOV;
按照公式:依次计算当扫描床的高度为标称床高时,所述扫描床在每一幅预览图像中的像素纵坐标值其中,公式表示取整运算,表示每个像素所代表的实际宽度,cYi表示每一幅预览图像的图像中心坐标的纵坐标值。
5.根据权利要求2所述的方法,其特征在于,根据预览图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、预览图像重建视野和所述确定出行序号最大的行,计算扫描床在每一幅预览图像建立时的实际高度,包括:
获取的所述扫描床标称床高H2、预览图像矩阵大小M×M、机架的旋转中心高度H1、每一幅预览图像的图像中心坐标(cXi,cYi)、预览图像重建视野FOV和所述确定出行序号最大的行数Ycouch;
按照公式:依次计算扫描床在每一幅预览图像建立时的实际高度Hi,其中,为每个像素所代表的实际宽度。
6.根据权利要求1所述的方法,其特征在于,所述测量扫描床在z位置的实际高度以及对应的扫描床下垂偏移量,包括:
获取机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和图像矩阵大小;
依次搜索平片图像中每一列中的行像素值大于阈值的行序号,确定出每一列中行序号最大的行;
根据平片图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和所述确定出的行序号最大的行,计算扫描床在每一个z位置的实际高度;
依次计算扫描床在每一个z位置的实际高度与扫描床的标称床高之差,所述差值作为对应的扫描床下垂偏移量。
7.根据权利要求6所述的方法,其特征在于,获取机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和图像矩阵大小,之后还包括:
根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、平片图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在平片图像中的理论像素纵坐标值;
则所述依次搜索平片图像中每一列中的行像素值大于阈值的行序号,确定出每一列中行序号最大的行,具体为:
依次根据计算出的扫描床在平片图像中的理论像素纵坐标值,在行序号为纵坐标值的上下任意多行进行搜索,搜索出存在行像素值大于阈值的行序号,确定出每一列中行序号最大的行。
8.根据权利要求7所述的方法,其特征在于,所述根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、平片图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在平片图像中的理论像素纵坐标值,包括:
获取的所述扫描床标称床高H2、平片图像矩阵大小N×M、机架的旋转中心高度H1、平片图像的图像中心坐标(cX,cY)、图像重建视野FOV;
按照公式:计算当扫描床的高度为标称床高时,所述扫描床在平片图像中的像素纵坐标值其中,公式表示取整运算,表示每个像素所代表的实际宽度,cY表示平片图像的图像中心坐标的纵坐标值。
9.根据权利要求6所述的方法,其特征在于,所述根据平片图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和所述确定出的行序号最大的行,计算扫描床在每一个z位置的实际高度,包括:
获取的所述扫描床标称床高H2、平片图像矩阵大小N×M、机架的旋转中心高度H1、平片图像的图像中心坐标(cX,cY)、平片图像重建视野FOV和所述确定出行序号最大的行数Ycouch;
按照公式:依次计算扫描床在每个z位置Zi处的实际高度Hi,其中,为每个像素所代表的实际宽度。
10.根据权利要求1所述的方法,其特征在于,所述当前目标图像所在直线通过以下方式得到,包括:
采用最小二乘法对扫描床在每一幅图像建立时所在的刻度位置Zi以及位置Zi对应的下垂偏移量D(Zi)进行曲线拟合,拟合出扫描床偏移所在直线y=kz+b;
根据当前目标图像所在直线与扫描床偏移所在直线的垂直关系,计算出目标图像所在的直线其中为扫描床偏移所在直线在z轴上的截距;
按照所述当前目标图像所在的直线计算出目标图像每一行的坐标信息(y(n),z(n))为:
其中,y(n)和z(n)分别代表第n行的y轴和z轴坐标值,n=0,…,M-1表示行序号。
11.一种图像位置补偿装置,其特征在于,包括:
测量单元,用于测量扫描床在各个z位置的实际高度以及对应的扫描床下垂偏移量;
调整单元,用于根据所述扫描床下垂偏移量,在后重建过程中,调整每一幅图像的建像中心纵坐标,得到调整后的每一幅重建图像;
插值单元,用于从所述调整后的每一幅重建图像中获取的每一行的坐标信息和对应的像素值,采用插值的方式得到目标图像每一行的像素值;所述插值单元,包括:
第一坐标模块,用于读取所述调整后的每一幅重建图像中每一行的坐标(Y(m),Z(m))为:
第二坐标模块,用于计算当前目标图像每一行的坐标信息(y(n),z(n)),其中,y(n)和z(n)分别代表第n行的y轴和z轴坐标值,n=0,…,M-1表示行序号;
第一处理模块,用于对读取的(Y(m),Z(m))做两个方向的插值处理,得到目标图像坐标为(y(n),z(n))图像像素值。
12.根据权利要求11所述的装置,其特征在于,所述测量单元,包括:
第一获取模块,用于获取机架的旋转中心高度、每一幅预览图像的图像中心坐标、图像重建视野和图像矩阵大小;
第一搜索模块,用于依次搜索每一幅预览图像中像素值大于阈值的行序号,确定出每一预览幅图像中行序号最大的行;
第一计算模块,用于根据预览图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、图像重建视野和所述确定出的行序号最大的行,计算扫描床在每一幅预览图像建立时的实际高度;
第一差值模块,用于依次计算扫描床在每一幅预览图像建立时的实际高度与扫描床的标称床高之差,所述差值作为对应的扫描床下垂偏移量。
13.根据权利要求12所述的装置,其特征在于,所述测量单元,还包括:
第一估算模块,用于根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、每一幅预览图像的图像中心坐标、预览图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在每一幅图像中的理论像素纵坐标值;
则所述第一搜索模块,具体为:
依次根据计算出的扫描床在每一幅预览图像中的理论像素纵坐标值,在行序号为纵坐标值的上下任意多行进行搜索,搜索出存在像素值大于阈值的行数,确定出行序号最大的行。
14.根据权利要求13所述的装置,其特征在于,所述第一估算模块,包括:
第一获取子模块,用于获取的所述扫描床标称床高H2、预览图像矩阵大小M×M、机架的旋转中心高度H1、每一幅预览图像的图像中心坐标(cXi,cYi)、图像重建视野FOV;
第一估算子模块,用于按照公式:依次计算当扫描床的高度为标称床高时,所述扫描床在每一幅预览图像中的理论像素纵坐标值其中,公式表示取整运算,表示每个像素所代表的实际宽度,cYi表示每一幅预览图像的图像中心坐标的纵坐标值。
15.根据权利要求12所述的装置,其特征在于,所述第一计算模块,包括:
第一提取子模块,用于获取的所述扫描床标称床高H2、预览图像矩阵大小M×M、机架的旋转中心高度H1、每一幅预览图像的图像中心坐标(cXi,cYi)、预览图像重建视野FOV和所述确定出行序号最大的行数Ycouch;
第一计算子模块,用于按照公式:依次计算扫描床在每一幅预览图像建立时的实际高度Hi,其中,为每个像素所代表的实际宽度。
16.根据权利要求11所述的装置,其特征在于,所述测量单元,包括:
第二获取模块,用于获取机架的旋转中心高度、平片图像的图像中心坐标、图像重建视野和图像矩阵大小;
第二搜索模块,用于依次搜索平片图像中每一列中的行像素值大于阈值的行序号,确定出每一列中行序号最大的行;
第二计算模块,用于根据平片图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、平片图像重建视野和所述确定出的行序号最大的行,计算扫描床在每一个z位置的实际高度;
第二差值模块,用于依次计算扫描床在每一个z位置的实际高度与扫描床的标称床高之差,所述差值作为对应的扫描床下垂偏移量。
17.根据权利要求16所述的装置,其特征在于,所述测量单元,还包括:
第二估算模块,用于根据所述扫描床标称床高、图像矩阵大小、机架的旋转中心高度、平片图像的图像中心坐标、平片图像重建视野,计算出当扫描床的高度为标称床高时,所述扫描床在平片图像中的理论像素纵坐标值;
则所述第二搜索模块,具体为:
依次根据计算出的扫描床在平片图像中的理论像素纵坐标值,在行序号为纵坐标值的上下任意多行进行搜索,搜索出存在行像素值大于阈值的行序号,确定出每一列中行序号最大的行。
18.根据权利要求17所述的装置,其特征在于,所述第二估算模块,包括:
第二获取子模块,用于获取的所述扫描床标称床高H2、平片图像矩阵大小N×M、机架的旋转中心高度H1、平片图像的图像中心坐标(cX,cY)、平片图像重建视野FOV;
第二估算子模块,用于按照公式:计算当扫描床的高度为标称床高时,所述扫描床在平片图像中的像素纵坐标值其中,公式表示取整运算,表示每个像素所代表的实际宽度,cY表示平片图像的图像中心坐标的纵坐标值。
19.根据权利要求16所述的装置,其特征在于,所述第二计算模块,包括:
第二提取子模块,用于获取的所述扫描床标称床高H2、平片图像矩阵大小N×M、机架的旋转中心高度H1、平片图像的图像中心坐标(cX,cY)、平片图像重建视野FOV和所述确定出行序号最大的行数Ycouch;
第二计算子模块,用于按照公式:依次计算扫描床在每个z位置Zi处的实际高度Hi,其中,为每个像素所代表的实际宽度。
20.根据权利要求11所述的装置,其特征在于,所述第二坐标模块,包括:
拟合子模块,用于采用最小二乘法对扫描床在每一幅图像建立时所在的刻度位置Zi以及位置Zi对应的下垂偏移量D(Zi)进行曲线拟合,拟合出扫描床偏移所在直线y=kz+b;
直线子模块,用于根据当前目标图像所在直线与扫描床偏移所在直线的垂直关系,计算出目标图像所在的直线其中为扫描床偏移所在直线在z轴上的截距;
坐标子模块,用于按照所述当前目标图像所在的直线
计算出目标图像每一行的坐标信息(y(n),z(n))为:
其中,y(n)和z(n)分别代表第n行的y轴和z轴坐标值,n=0,…,M-1表示行序号。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310273523.9A CN103340644B (zh) | 2013-06-28 | 2013-06-28 | 一种图像位置补偿方法及装置 |
US14/142,945 US9355437B2 (en) | 2013-06-28 | 2013-12-30 | Image position compensation method and apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310273523.9A CN103340644B (zh) | 2013-06-28 | 2013-06-28 | 一种图像位置补偿方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103340644A CN103340644A (zh) | 2013-10-09 |
CN103340644B true CN103340644B (zh) | 2015-07-22 |
Family
ID=49275511
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310273523.9A Active CN103340644B (zh) | 2013-06-28 | 2013-06-28 | 一种图像位置补偿方法及装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US9355437B2 (zh) |
CN (1) | CN103340644B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6072097B2 (ja) * | 2015-01-30 | 2017-02-01 | キヤノン株式会社 | 放射線撮影装置、制御装置、長尺撮影システム、制御方法、及びプログラム |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101879069A (zh) * | 2009-05-08 | 2010-11-10 | 株式会社东芝 | X射线计算机断层摄影装置 |
CN102805643A (zh) * | 2011-05-31 | 2012-12-05 | 上海西门子医疗器械有限公司 | 检查床弯曲的补偿 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7020315B2 (en) * | 2001-01-10 | 2006-03-28 | Elgems Ltd. | Sag correction |
JP2004180846A (ja) * | 2002-12-02 | 2004-07-02 | Shimadzu Corp | X線ct装置 |
US7111985B2 (en) * | 2004-11-04 | 2006-09-26 | General Electric Company | Method and system for measuring table sag |
JP2006288472A (ja) * | 2005-04-06 | 2006-10-26 | Toshiba Corp | X線ct装置及び画像再構成方法 |
US8135198B2 (en) * | 2007-08-08 | 2012-03-13 | Resonant Medical, Inc. | Systems and methods for constructing images |
WO2010150717A1 (ja) * | 2009-06-23 | 2010-12-29 | 株式会社 日立メディコ | X線ct装置 |
-
2013
- 2013-06-28 CN CN201310273523.9A patent/CN103340644B/zh active Active
- 2013-12-30 US US14/142,945 patent/US9355437B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101879069A (zh) * | 2009-05-08 | 2010-11-10 | 株式会社东芝 | X射线计算机断层摄影装置 |
CN102805643A (zh) * | 2011-05-31 | 2012-12-05 | 上海西门子医疗器械有限公司 | 检查床弯曲的补偿 |
Non-Patent Citations (1)
Title |
---|
重力因素对模拟CT定位精度的影响;徐微 等;《实用癌症杂志》;20070131;第22卷(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
US9355437B2 (en) | 2016-05-31 |
CN103340644A (zh) | 2013-10-09 |
US20150003705A1 (en) | 2015-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101473538B1 (ko) | 엑스선 이미지들 내의 마커 위치들로부터의 환자 모션 벡터들의 추출 | |
US7986823B2 (en) | System and method for consistent detection of mid-sagittal planes for magnetic resonance brain scans | |
Stevens et al. | Alignment of a volumetric tomography system | |
US20110123081A1 (en) | Correcting and reconstructing x-ray images using patient motion vectors extracted from marker positions in x-ray images | |
RU2602750C1 (ru) | Способ калибровки компьютерно-томографического изображения, устройство и система компьютерной томографии | |
US9336569B2 (en) | Image compensation method and apparatus | |
CN103385708B (zh) | 一种用于背部躯干的测量与图像处理装置 | |
CN103654833B (zh) | Ct探测器偏转角的确定方法和装置 | |
JP6711410B2 (ja) | 放射線断層撮影装置の撮像倍率校正方法 | |
CN104771189A (zh) | 三维头颅图像对正方法及装置 | |
US20130223703A1 (en) | Medical image processing apparatus | |
CN105319225B (zh) | 一种实现板状样品高分辨率大视野cl成像的扫描方法 | |
CN103860182A (zh) | 一种检查床的位置差异确定方法、系统及医疗设备 | |
CN103340644B (zh) | 一种图像位置补偿方法及装置 | |
EP2745767B1 (en) | Improvements in and relating to ophthalmoscopes | |
CN103720482B (zh) | 影像重建方法与系统及影像建构方法与系统 | |
JP2000046761A (ja) | 物体のイメ―ジング方法及びシステム | |
CN103606144A (zh) | Ct投影旋转中心的确定方法和装置 | |
US7111985B2 (en) | Method and system for measuring table sag | |
CN108364325A (zh) | 规则样品x射线ct投影图位置平移偏差检测及校正方法 | |
CN105321206B (zh) | 一种适用于中子层析成像系统的旋转轴线偏摆角的误差补偿方法 | |
KR101426245B1 (ko) | 3d 단층 촬영장치 | |
CN104395934B (zh) | 交错的多能量成像中的图像重建 | |
CN111184535A (zh) | 手持式无约束扫描无线三维超声实时体素成像系统 | |
US20110103545A1 (en) | Computed tomography method and apparatus for centre-of-rotation determination |
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 | ||
CP03 | Change of name, title or address |
Address after: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Patentee after: DongSoft Medical System Co., Ltd. Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16 Patentee before: Dongruan Medical Systems Co., Ltd., Shenyang |
|
CP03 | Change of name, title or address |