CN105094725A

CN105094725A - 图像显示方法

Info

Publication number: CN105094725A
Application number: CN201410202944.7A
Authority: CN
Inventors: 张丽; 陈志强; 黄清萍; 孙运达; 唐智; 沈乐; 唐虎; 任乾鲁; 金鑫
Original assignee: Nuctech Co Ltd
Current assignee: Nuctech Co Ltd
Priority date: 2014-05-14
Filing date: 2014-05-14
Publication date: 2015-11-25
Anticipated expiration: 2034-05-14
Also published as: JP2018128458A; JP2019144250A; EP3021283B1; EP3021283A2; HK1217796A1; WO2015172726A1; JP2015227872A; CN105094725B; US9805504B2; US20150332498A1; RU2599596C1; JP6770121B2; EP3021283A3

Abstract

公开了一种CT系统的图像显示方法。在该方法中，对被检查物体进行CT扫描，得到CT投影数据，按照预定的间隔对CT投影数据进行整理，以一固定角度作为起始角度并以360度为间隔，从整理后的CT投影数据中抽取基本数据，基于所抽取的基本数据形成DR图像，从CT投影数据重建被检查物体的三维图像，以及在屏幕上同时显示DR图像和重建的三维图像。该方案对CT数据进行处理，得到DR数据。得到DR数据后，直接利用DR数据处理算法，得到DR图像。这样能够让判图员利用已有的DR图像判图经验，更加准确，更加快速地检查旅客携带的物品。

Description

图像显示方法

技术领域

本申请涉及计算机断层成像(CT)技术，具体涉及一种CT系统中的图像显示方法。

背景技术

美国的9·11事件发生以后，航空领域的安全检查工作越来越受到重视。基于CT技术的安检产品有着不可比拟的先进性，性能指标优势明显。美国TSA认为，CT设备及其检测方法是实现机场安全检查的重要途径。普通的物品机，即X光机透视成像能够得到清晰的DR图像，这样的DR(DigitalRadiography)图像已经在用户的心目中形成了固定的印象。用户已经非常熟悉DR图像了。在DR系统中，X射线源和探测器是固定的，是不旋转的。这些DR图像通常是顶视角图像或者底视角图像，非常直观，方便看图。读图和判图所形成的丰富经验，都是基于这些直观的DR图像。一位优秀的判图员，需要阅览过数万幅图像以上的实际DR图像。往往需要90个工作日以上的时间，才能培养出一位优秀的判图员。目前，已经有为数众多的优秀的判图员在医疗、安检等行业工作着。没有DR图像时，这些经验丰富的老判图员无法进行正确的判图。

现有的一些CT系统虽然有现成的三维数据，但是没有DR图像。当判图员一下子看到三维图像时，觉得不习惯。例如，由于观察角度不合适，违禁品(例如手枪)在三维图像中会是一个小的聚团，形态不清楚。判图员突然看到三维图像时，搞不清楚怎么手枪会突然成了一个小聚团。判图员单独观察三维图像，仍然不清楚旅客携带的物品究竟是什么。因此，即使已经存在三维图像时，判图员仍然希望能够将三维图像与以前曾经很熟悉的DR图像进行比照，对比这二种图像，并且利用已有的DR图像的丰富的判图经验。即，判图员希望这样布置：左边放一个显示器，显示DR图像，右边放一个显示器，显示CT图像。利用了已有的DR图像判图经验，众多判图员就可以更快速，更方便，更准确地知道旅客携带的物品究竟是什么。这样就可以节省了大量的人力，节省了时间。

此外，对于已有DR系统的CT设备，用户选择了特定位置，要求DR-CT配准，显示出特定位置的切片，可以依靠硬件装置的皮带编码，得知特定位置是行李包裹的头部，中部或者尾部。但是，由于皮带编码在传输过程中丢失或缺失，或者有时不能匹配，导致该方法无法准确进行配准。迫切需要有新的有效的配准方法。

发明内容

考虑到现有技术中的一个或多个问题，提出了一种CT系统的图像显示方法，能够在CT系统不具备DR成像设备的情况下，与被检查物体的三维图像一起显示DR图像。

在本发明的一个方面，提出了一种CT系统的图像显示方法，包括步骤：对被检查物体进行CT扫描，得到CT投影数据；按照预定的间隔对CT投影数据进行整理；以一固定角度作为起始角度并以360度为间隔，从整理后的CT投影数据中抽取基本数据；基于所抽取的基本数据形成DR图像；从CT投影数据重建被检查物体的三维图像；以及在屏幕上同时显示DR图像和重建的三维图像。

根据一些实施例，所述的图像显示方法还包括步骤：以所述固定角度加上180度作为起始角度并以360度为间隔，从整理后的CT投影数据中抽取补充数据；将所抽取的补充数据存储在矩阵中；对所述矩阵做左右镜像处理；以行为单位，将镜像处理后的矩阵中的补充数据间隔插入所述基本数据中，得到补充后的基本数据；基于补充后的基本数据形成DR图像。

根据一些实施例，所述的图像显示方法还包括步骤：以另一固定角度作为起始角度并以360度为间隔，从整理后的CT投影数据中抽取基本数据；基于所抽取的基本数据形成另一DR图像；以及在屏幕上同时显示所述DR图像、所述另一DR图像和重建的三维图像。

根据一些实施例，按照预定的间隔对CT投影数据进行整理的步骤包括：当采样频率高于每预定间隔投影数据1次时，按照欠采样规则，将CT投影数据调整为每预定间隔投影数据1次；当采样频率低于每预定间隔投影数据1次时，按照线性插值，将CT投影数据调整为每预定间隔投影数据1次。

根据一些实施例，当所述基本数据为高能和低能数据时，按照物质识别算法对所述基本数据进行处理得到灰度数据和材料数据。

根据一些实施例，将DR图像插值成为用户要求的分辨率，在屏幕上显示。

根据一些实施例，重建的三维图像的数据是以下之一：物理密度、电子密度数据、原子序数数据、高能数据、低能数据、灰度数据、材料数据。

在本发明的另一方面，提出了一种CT系统的图像显示方法，包括步骤：对被检查物体进行CT扫描，得到CT投影数据；从CT投影数据重建被检查物体的三维图像；按照一视角对被检查物体的三维图像数据H(x,y,z)进行投影，得到该视角下的DR数据，其中三维图像数据H(x,y,z)各维的尺寸为X×Y×Z，维度x从1到X变化，该变化为水平面内垂直于皮带运动的方向，维度y从1到Y变化，该变化为竖直向上的方向，维度z从1到Z变化，该变化为水平面内沿皮带运动的方向；基于该视角下的DR数据形成DR图像；在屏幕上同时显示DR图像和三维图像。

根据一些实施例，所述的图像显示方法还包括步骤：按照另一视角对被检查物体的三维图像数据H(x,y,z)进行投影，得到该视角下的DR数据；基于该另一视角下的DR数据形成DR图像；在屏幕上同时显示两幅DR图像和三维图像。

根据一些实施例，所述视角和所述另一视角是从水平面内垂直于皮带运动的方向、竖直方向和水平面内沿皮带运动的方向中选择的。

根据一些实施例，将三维数据H(x,y,z)沿着维度x的方向投影，得到二维数据I(y,z)，投影公式如下：

I (y, z) = Σ_{x = 1}^{X} H (x, y, z)

该二维数据为侧视角的DR数据，数据维度尺寸为Y×Z。

根据一些实施例，将三维数据H(x,y,z)沿着维度y的方向投影，得到二维数据J(x,z)，投影公式如下：

J (x, z) = Σ_{y = 1}^{Y} H (x, y, z)

该二维数据为底视角或顶视角的DR数据，数据维度尺寸为X×Z。

根据一些实施例，将三维数据H(x,y,z)沿着维度z即皮带运动的方向投影，得到二维数据K(x,y)，投影公式如下：

K (x, y) = Σ_{z = 1}^{Z} H (x, y, z)

该二维数据为前视角的DR数据，数据维度尺寸为X×Y。

根据上述方案，对CT投影数据进行处理，得到DR数据，得到DR数据后，利用DR数据处理算法，得到DR图像。这样能够让判图员利用已有的DR图像判图经验，更加准确，更加快速地检查旅客携带的物品。

在本发明的又一方面，提出了一种在DR-CT系统中进行图像配准的方法，所述DR-CT成像系统包括DR设备和CT设备，所述方法包括步骤：用CT设备对被检查物体进行扫描以得到CT投影数据；从所述CT投影数据重建被检查物体的三维图像；用DR设备对被检查物体进行扫描以得到DR图像；在屏幕上显示所述DR图像和所述三维图像；响应于用户对DR图像中的特定位置的选择，从DR图像的数据中抽取相应的列；从所述三维图像中得到与所述DR图像相同视角下的辅助DR数据；计算所述DR图像中的所述列与辅助DR数据中的各列之间的相关性；将辅助DR数据中与所述DR图像中的所述列之间相关性最大的那一列对应的切片图像与所述DR图像同时显示在屏幕上。

根据一些实施例，从所述三维图像中得到与所述DR图像相同视角下的辅助DR数据的步骤包括：沿着维度y方向对被检查物体的三维图像数据H(x,y,z)进行投影，得到该视角下的DR数据，其中三维图像数据H(x,y,z)各维的尺寸为X×Y×Z，维度x从1到X变化，该变化为水平面内垂直于皮带运动的方向，维度y从1到Y变化，该变化为竖直向上的方向，维度z从1到Z变化，该变化为水平面内沿皮带运动的方向。

J (x, z) = Σ_{y = 1}^{Y} H (x, y, z)

在本发明的再一方面，提出了一种在DR-CT系统中进行图像配准的方法，所述DR-CT成像系统包括DR设备和CT设备，所述方法包括步骤：用CT设备对被检查物体进行扫描以得到CT投影数据；用DR设备对被检查物体进行扫描以得到DR图像；在屏幕上显示所述DR图像；响应于用户对DR图像中的特定位置的选择，从DR图像的数据中抽取相应的列；从所述CT投影数据中得到与所述DR图像相同视角下的辅助DR数据；计算所述DR图像中的所述列与辅助DR数据中的各列之间的相关性；从所述CT投影数据重建被检查物体的三维图像；以及将辅助DR数据中与所述DR图像中的所述列之间相关性最大的那一列对应的切片图像与所述DR图像同时显示在屏幕上。

根据一些实施例，从所述CT投影数据中得到与所述DR图像相同视角下的辅助DR数据的步骤包括：按照预定的间隔对CT投影数据进行整理；以一固定角度作为起始角度并以360度为间隔，从整理后的CT投影数据中抽取得到所述辅助数据。

根据一些实施例，所述的方法还包括步骤：以所述固定角度加上180度作为起始角度并以360度为间隔，从整理后的CT投影数据中抽取补充数据；将所抽取的补充数据存储在矩阵中；对所述矩阵做左右镜像处理；以行为单位，将镜像处理后的矩阵中的补充数据间隔插入所述辅助DR数据中，得到补充后的辅助DR数据。

附图说明

为了更好地理解本发明，将根据以下附图对本发明进行详细描述：

图1示出了根据本发明实施例的CT设备的结构示意图；

图2示出了如图1所示的计算机数据处理器的结构框图；

图3示出了根据本发明实施方式的控制器的结构框图；

图4是描述根据本发明一个实施例的在CT系统中显示图像的方法的流程图；

图5是描述整理CT投影数据的过程的示意图；

图6示出了从CT投影数据得到的DR图像的示意图；

图7是描述根据本发明另一实施例的在CT系统中显示图像的方法的流程图；

图8是描述根据本发明又一实施例的在CT系统中配准DR图像和CT图像的方法的流程图；以及

图9是描述根据本发明再一实施例的在CT系统中配准DR图像和CT图像的方法的流程图。

具体实施方式

下面将详细描述本发明的具体实施例，应当注意，这里描述的实施例只用于举例说明，并不用于限制本发明。在以下描述中，为了提供对本发明的透彻理解，阐述了大量特定细节。然而，对于本领域普通技术人员显而易见的是：不必采用这些特定细节来实行本发明。在其他实例中，为了避免混淆本发明，未具体描述公知的结构、材料或方法。

在整个说明书中，对“一个实施例”、“实施例”、“一个示例”或“示例”的提及意味着：结合该实施例或示例描述的特定特征、结构或特性被包含在本发明至少一个实施例中。因此，在整个说明书的各个地方出现的短语“在一个实施例中”、“在实施例中”、“一个示例”或“示例”不一定都指同一实施例或示例。此外，可以以任何适当的组合和/或子组合将特定的特征、结构或特性组合在一个或多个实施例或示例中。此外，本领域普通技术人员应当理解，这里使用的术语“和/或”包括一个或多个相关列出的项目的任何和所有组合。

针对现有技术中一些CT设备不具备DR成像模块但是在判图过程中判图员又需要显示DR图像的问题，本发明的一些实施例提出了从CT投影数据中得到被检查物体的DR图像，与重建的三维图像一起在屏幕上显示，方便判图员利用已有的判图经验来进行安全检查。根据上述方案，能够更加准确，更加快速地检查旅客携带的物品。在一些实施例中，从重建的三维图像的数据中抽取特定视角下的投影数据，将这些投影数据组合成为该视角下的DR数据。在其他实施例中，对重建的三维图像的数据按照不同的视角进行投影，得到该视角下的DR数据，形成DR图像。

在另外的实施例中，针对一些配备了DR设备的CT设备中由于皮带编码缺失或者丢失造成的无法对DR图像和CT图像进行配准的问题，提出了从CT投影数据中或者重建的三维图像的数据中得到特定视角下的DR图像，并且将DR设备得到的DR图像中的所选列与从三维图像的数据中得到的DR图像进行相关性分析，确定相关性最大的那一列所在的切片与DR图像一起显示出来。

图1是根据本发明实施方式的CT设备的结构示意图。如图1所示，根据本实施方式的CT设备包括：机架20、承载机构40、控制器50、计算机数据处理器60等。机架20包括发出检查用X射线的射线源10，诸如X光机，以及探测和采集装置30。承载机构40承载被检查行李70穿过机架20的射线源10与探测和采集装置30之间的扫描区域，同时机架20围绕被检查行李70的前进方向转动，从而由射线源10发出的射线能够透过被检查行李70，对被检查行李70进行CT扫描。

探测和采集装置30例如是具有整体模块结构的探测器及数据采集器，例如平板探测器，用于探测透射被检液态物品的射线，获得模拟信号，并且将模拟信号转换成数字信号，从而输出被检查行李70针对X射线的投影数据。控制器50用于控制整个系统的各个部分同步工作。计算机数据处理器60用来处理由数据采集器采集的数据，对数据进行处理并重建，输出结果。

如图1所示，射线源10置于可放置被检物体的一侧，探测和采集装置30置于被检查行李70的另一侧，包括探测器和数据采集器，用于获取被检查行李70的多角度投影数据。数据采集器中包括数据放大成形电路，它可工作于(电流)积分方式或脉冲(计数)方式。探测和采集装置30的数据输出电缆与控制器50和计算机数据处理器60连接，根据触发命令将采集的数据存储在计算机数据处理器60中。

图2示出了如图1所示的计算机数据处理器60的结构框图。如图2所示，数据采集器所采集的数据通过接口单元68和总线64存储在存储器61中。只读存储器(ROM)62中存储有计算机数据处理器的配置信息以及程序。随机存取存储器(RAM)63用于在处理器66工作过程中暂存各种数据。另外，存储器61中还存储有用于进行数据处理的计算机程序。内部总线64连接上述的存储器61、只读存储器62、随机存取存储器63、输入装置65、处理器66、显示装置67和接口单元68。

在用户通过诸如键盘和鼠标之类的输入装置65输入的操作命令后，计算机程序的指令代码命令处理器66执行预定的数据处理算法，在得到数据处理结果之后，将其显示在诸如LCD显示器之类的显示装置67上，或者直接以诸如打印之类硬拷贝的形式输出处理结果。

图3示出了根据本发明实施方式的控制器的结构框图。如图3所示，控制器50包括：控制单元51，根据来自计算机60的指令，来控制射线源10、承载机构40和探测和采集装置30；触发信号产生单元52，用于在控制单元的控制下产生用来触发射线源10、探测和采集装置30以及承载机构40的动作的触发命令；第一驱动设备53，它在根据触发信号产生单元52在控制单元51的控制下产生的触发命令驱动承载机构40传送被检查行李70；第二驱动设备54，它根据触发信号产生单元52在控制单元51的控制下产生的触发命令机架20旋转。探测和采集装置30获得的投影数据存储在计算机60中进行CT断层图像重建，从而获得被检查行李70的断层图像数据。然后计算机60例如通过执行软件来从断层图像数据得到被检查行李70的至少一个视角下的DR图像，与重建的三维图像一起显示，方便判图员进行安全检查。根据其他实施例，上述的CT成像系统也可以是双能CT系统，也就是机架20的X射线源10能够发出高能和低能两种射线，探测和采集装置30探测到不同能量水平下的投影数据后，由计算机数据处理器60进行双能CT重建，得到被检查行李70的各个断层的等效原子序数和等效电子密度数据。

图4是描述根据本发明一个实施例的在CT系统中显示图像的方法的流程图。如图4所示，在步骤S41，通过CT系统对被检查物体70进行CT扫描，得到CT投影数据。这里的CT投影数据可以是，但不限于是以下的几种：物理密度或电子密度数据；原子序数数据；高能数据；低能数据；灰度数据；材料数据。该实施例中，探测器的排数以1排为例。当探测器为多排时，将多排看做为一个整体，当作是1排，进行数据处理。

在步骤S42，按照预定的间隔对CT投影数据进行整理。例如，CT投影数据以1度为间隔，进行整理。即：每间隔1度，有投影数据1次。在一些实施例中，当采样频率高于每间隔1度投影数据1次时，按照欠采样规则，调整为每间隔1度投影数据1次。当采样频率过低时，按照线性插值，调整为每间隔1度投影数据1次。

例如，可以将CT投影数据，即正弦图，整理排列为矩阵A，其中的元素记为A(m,n)。n表示探测器阵元，n变化为探测器阵列的1～N，N为探测器阵列中阵元的个数。m表示角度，m变化为角度数变化，即矩阵A的第1行为第1度的探测器阵列的1～N数据，第2行为第2度的探测器阵列的1～N数据。所有的角度以360为周期，周而复始，排列情况如下：1度～360度，361度～720度，721度～1080度，1081度～1440度……，如图5所示。

在步骤S43，以一固定角度作为起始角度，并以360度为间隔，从整理后的CT投影数据中抽取基本数据。例如，从CT投影数据中，按照圆周期360度为间隔，以固定的起始角度S，抽取若干个角度的数据，例如：抽取了1度，361度，721度，1081度，1441度。即：假设固定的起始角度为S，则第i次抽取，对于的角度数Angle为

Angle＝S+i×360

其中i＝0,1,2,……(1)

将抽取出的基本数据，重新组合成为矩阵B。

在步骤S44，基于所抽取的基本数据形成DR图像，如图6所示。例如，将重新组合形成的矩阵B，按照已有的DR数据处理算法，处理形成DR图像。当数据为高能和低能数据时，按照物质识别算法，得到灰度数据和材料数据。

在步骤S45，计算机数据处理器60运行软件从CT投影数据重建被检查物体的三维图像。

在步骤S46，在屏幕上同时显示DR图像和重建的三维图像。一些情况下，将DR图像插值成为用户要求的分辨率，在显示器上显示。

在其他实施例中，如果矩阵B的数据量不够，需要增加数据时，按照以下公式，继续抽取。假设固定的起始角度为S，则第i次抽取，对于的角度数Angle为

Angle＝S+i×360+180

其中i＝0,1,2,……(2)将抽取出的数据，重新组合成为矩阵B1。

然后，对矩阵B1做左右镜像处理，即左右对称的元素，相互交换，公式如下：

B1(m,1)<＝＝>B1(m,N)

B1(m,2)<＝＝>B1(m,N-1)(3)

……

接下来，将得到的矩阵B1，以行为单位，间隔插入矩阵B，与矩阵B组合，形成矩阵C。数据维度尺寸为M×N。组合的公式如下：

最后，将矩阵C按照已有的DR数据处理算法，处理形成DR图像。当数据为高能和低能数据时，按照物质识别算法，得到灰度数据和材料数据。

在上述实施例中，固定的起始角度S，决定了DR图像的视角的方向。例如：S＝1时可以得到顶视角；S＝180时可以得到底视角；S＝90时可以得到左侧视角；S＝270时可以得到右侧视角。

以上的实施例中，得到的是一个视角下的DR图像。如果需要的话，可以以与上述的角度S不同的另一固定角度作为起始角度来形成该角度下的DR数据，得到相应的DR图像后与上述的DR图像和重建的三维图像一起显示在屏幕上，供判图员进行安全检查。

此外，一些实施例中，CT系统的探测器可以是一排或者多排，并且可以利用CT数据得到任意一个视角(包括顶视角、底视角、侧视角)的图像。或者可以利用CT数据得到任意二个即双视角的图像。以及可以利用CT数据得到任意多个即多视角的图像。

图7是描述根据本发明另一实施例的在CT系统中显示图像的方法的流程图。如图7所示，在步骤S71，利用CT系统对被检查物体进行CT扫描，得到CT投影数据，然后在步骤S72，通过计算机60从CT投影数据重建被检查物体的三维图像。这里的重建后的三维数据可以是，但不限于是以下的几种：物理密度或电子密度数据；原子序数数据；高能数据；低能数据；灰度数据；材料数据。

在步骤S73，按照一视角对被检查物体的三维图像数据H(x,y,z)进行投影，得到该视角下的DR数据，其中三维图像数据H(x,y,z)各维的尺寸为X×Y×Z，维度x从1到X变化，该变化为水平面内垂直于皮带运动的方向，维度y从1到Y变化，该变化为竖直向上的方向，维度z从1到Z变化，该变化为水平面内沿皮带运动的方向；

例如，按照需要的视角，进行投影，具体分为3种情况：

(1)将三维数据H(x,y,z)沿着维度x的方向投影，得到二维数据I(y,z)，投影公式如下：

I (y, z) = Σ_{x = 1}^{X} H (x, y, z) - - - (5)

该二维数据为侧视角的DR数据。数据维度尺寸为Y×Z。

(2)将三维数据H(x,y,z)沿着维度y的方向投影，得到二维数据J(x,z)，投影公式如下：

J (x, z) = Σ_{y = 1}^{Y} H (x, y, z) - - - (6)

该二维数据为底视角或顶视角的DR数据。数据维度尺寸为X×Z。

(3)将三维数据H(x,y,z)沿着维度z即皮带运动的方向投影，得到二维数据K(x,y)，投影公式如下：

K (x, y) = Σ_{z = 1}^{Z} H (x, y, z) - - - (7)

该二维数据为前视角的DR数据。数据维度尺寸为X×Y。

在步骤S74，基于该视角下的DR数据形成DR图像。例如，得到DR数据后，按照已有的DR数据处理算法，处理形成DR图像。当数据为高能和低能数据时，按照物质识别算法，得到灰度数据和材料数据。

在步骤S75，在屏幕上同时显示DR图像和三维图像。例如将DR图像插值成为用户要求的分辨率，在显示器上显示。

虽然上面描述的是CT系统不具备DR设备的情况，但是在一些具备DR设备的CT系统中，需要配准DR图像与切片图像。例如已有DR系统的CT设备，要求当用户点击DR图像的某个地方后，能够从三维数据中得到准确的切片位置是位于行李包裹的头部，中部或者尾部，从而显示出该切片。三维数据中有若干个切片，为了从三维数据中得到准确的位置来显示切片，可以依靠硬件装置的皮带编码，得知是行李包裹的头部，中部或者尾部。但是，由于皮带编码在传输过程中丢失或缺失，或者有时不能匹配，该方法无法准确进行配准。

针对上述问题，本发明的实施例提出了将DR设备得到的DR图像中的一列与从重建的三维图像中抽取的DR数据中的每一列进行相关性计算，将相关性最大的那一列所对应的切片图像与DR图像一起显示在屏幕上。这里的DR数据可以是，但不限于是以下的几种：物理密度或电子密度数据；原子序数数据；高能数据；低能数据；灰度数据；材料数据等等。

图8是描述根据本发明又一实施例的在CT系统中配准DR图像和CT图像的方法的流程图。

在步骤S81，用CT设备对被检查物体进行扫描以得到CT投影数据。在步骤S82，通过计算机60执行软件来从CT投影数据重建被检查物体的三维图像。

在步骤S82，用DR设备对被检查物体进行扫描以得到DR图像，然后，在步骤S83,在屏幕上显示DR设备得到的DR图像和重建的三维图像，方便判图员进行安全检查。

判图员操作诸如鼠标之类的输入装置65，点击DR图像中的某一列，希望看到与该列对应的切片图像，以便进行更为准确的检查。在步骤S85，响应于用户对DR图像中的特定位置的选择，从DR图像的数据中抽取相应的列。接下来，在步骤S86，计算机60从三维图像中得到与DR图像相同视角下的辅助DR数据。

例如DR系统生成DR数据，用户在该DR数据中用鼠标指定某个位置。该位置所在的DR数据沿探测器阵元方向扩展，记为Array0(n)，其中n从1～N变化。N为探测器阵列中阵元的个数。该数据实际为一个数组，存储了所有探测器阵元的数据。

按照如图7所示的方法，从CT数据中得到DR数据，记为矩阵J，其中的元素为J(x,z)。将矩阵J转置，得到J1，其中的元素为J1(z,x)，将J1按照要求的维度尺寸M×N插值，得到矩阵C，数据维度尺寸为M×N。

在一些实施例中，从三维图像中得到与所述DR图像相同视角下的辅助DR数据的步骤包括：沿着维度y方向对被检查物体的三维图像数据H(x,y,z)进行投影，得到该视角下的DR数据，其中三维图像数据H(x,y,z)各维的尺寸为X×Y×Z，维度x从1到X变化，该变化为水平面内垂直于皮带运动的方向，维度y从1到Y变化，该变化为竖直向上的方向，维度z从1到Z变化，该变化为水平面内沿皮带运动的方向。例如，将三维数据H(x,y,z)沿着维度y的方向投影，参考上述公式(6)得到二维数据J(x,z)。

然后，在步骤S87，计算机60计算所述DR图像中的所述列与辅助DR数据中的各列之间的相关性，并且将辅助DR数据中与所述DR图像中的所述列之间相关性最大的那一列对应的切片图像与所述DR图像同时显示在屏幕上。

例如，取出矩阵C的某个横行，例如m＝1时，n从1～N变化，记入数组Array1(n)。按照下面的公式，计算误差E：

E = Σ_{n = 1}^{N} {((Array 1 (n) - Array (0))}^{2} - - - (8)

重复上述步骤，取出下一个横行，例如m＝2时，n从1～N变化，记入数组Array1(n)，直到矩阵C所有的横行都被遍历。最后，从所有的E中，逐个查找，找到最小值对应的那个m，该位置匹配。取出三维数据H(x,y,z)的切片数据H(x,y,z＝m)与DR图像一起在显示器上显示。

以上描述的是在重建三维图像后得到辅助DR数据进行配准的实施例。类似地，也可以在重建三维图像之前得到辅助DR数据，然后用于图像的配准。图9是描述根据本发明再一实施例的在CT系统中配准DR图像和CT图像的方法的流程图。

在如图9所示的DR-CT系统中进行图像配准的方法中，DR-CT成像系统包括DR设备和CT设备。在步骤S91，用CT设备对被检查物体进行扫描以得到CT投影数据，在步骤S92，用DR设备对被检查物体进行扫描以得到DR图像。然后，在步骤S93，在屏幕上显示所述DR图像。

判图员操作诸如鼠标之类的输入装置65，点击DR图像中的某一列，希望看到与该列对应的切片图像，以便进行更为准确的检查。在步骤S94，响应于用户对DR图像中的特定位置的选择，从DR图像的数据中抽取相应的列。

在步骤S95，从所述CT投影数据中得到与所述DR图像相同视角下的辅助DR数据；

在步骤S96，计算所述DR图像中的所述列与辅助DR数据中的各列之间的相关性；

在步骤S97，从所述CT投影数据重建被检查物体的三维图像；以及

在步骤S98，将辅助DR数据中与所述DR图像中的所述列之间相关性最大的那一列对应的切片图像与所述DR图像同时显示在屏幕上。

在一些实施例中，从CT投影数据中得到与所述DR图像相同视角下的辅助DR数据的步骤包括：按照预定的间隔对CT投影数据进行整理；以一固定角度作为起始角度并以360度为间隔，从整理后的CT投影数据中抽取得到所述辅助数据。

在其他实施例中，上述的方法还包括步骤：以该固定角度加上180度作为起始角度并以360度为间隔，从整理后的CT投影数据中抽取补充数据。将所抽取的补充数据存储在矩阵中；对所述矩阵做左右镜像处理；以行为单位，将镜像处理后的矩阵中的补充数据间隔插入所述辅助DR数据中，得到补充后的辅助DR数据。

以上的详细描述通过使用示意图、流程图和/或示例，已经阐述了图像显示方法和配准图像的方法的众多实施例。在这种示意图、流程图和/或示例包含一个或多个功能和/或操作的情况下，本领域技术人员应理解，这种示意图、流程图或示例中的每一功能和/或操作可以通过各种结构、硬件、软件、固件或实质上它们的任意组合来单独和/或共同实现。在一个实施例中，本发明的实施例所述主题的若干部分可以通过专用集成电路(ASIC)、现场可编程门阵列(FPGA)、数字信号处理器(DSP)、或其他集成格式来实现。然而，本领域技术人员应认识到，这里所公开的实施例的一些方面在整体上或部分地可以等同地实现在集成电路中，实现为在一台或多台计算机上运行的一个或多个计算机程序(例如，实现为在一台或多台计算机系统上运行的一个或多个程序)，实现为在一个或多个处理器上运行的一个或多个程序(例如，实现为在一个或多个微处理器上运行的一个或多个程序)，实现为固件，或者实质上实现为上述方式的任意组合，并且本领域技术人员根据本公开，将具备设计电路和/或写入软件和/或固件代码的能力。此外，本领域技术人员将认识到，本公开所述主题的机制能够作为多种形式的程序产品进行分发，并且无论实际用来执行分发的信号承载介质的具体类型如何，本公开所述主题的示例性实施例均适用。信号承载介质的示例包括但不限于：可记录型介质，如软盘、硬盘驱动器、紧致盘(CD)、数字通用盘(DVD)、数字磁带、计算机存储器等；以及传输型介质，如数字和/或模拟通信介质(例如，光纤光缆、波导、有线通信链路、无线通信链路等)。

虽然已参照几个典型实施例描述了本发明，但应当理解，所用的术语是说明和示例性、而非限制性的术语。由于本发明能够以多种形式具体实施而不脱离发明的精神或实质，所以应当理解，上述实施例不限于任何前述的细节，而应在随附权利要求所限定的精神和范围内广泛地解释，因此落入权利要求或其等效范围内的全部变化和改型都应为随附权利要求所涵盖。

Claims

1.一种CT系统的图像显示方法，包括步骤：

对被检查物体进行CT扫描，得到CT投影数据；

按照预定的间隔对CT投影数据进行整理；

以一固定角度作为起始角度并以360度为间隔，从整理后的CT投影数据中抽取基本数据；

基于所抽取的基本数据形成DR图像；

从CT投影数据重建被检查物体的三维图像；以及

在屏幕上同时显示DR图像和重建的三维图像。

2.如权利要求1所述的图像显示方法，还包括步骤：

以所述固定角度加上180度作为起始角度并以360度为间隔，从整理后的CT投影数据中抽取补充数据；

将所抽取的补充数据存储在矩阵中；

对所述矩阵做左右镜像处理；

以行为单位，将镜像处理后的矩阵中的补充数据间隔插入所述基本数据中，得到补充后的基本数据；

基于补充后的基本数据形成DR图像。

3.如权利要求1所述的图像显示方法，还包括步骤：

以另一固定角度作为起始角度并以360度为间隔，从整理后的CT投影数据中抽取基本数据；

基于所抽取的基本数据形成另一DR图像；以及

在屏幕上同时显示所述DR图像、所述另一DR图像和重建的三维图像。

4.如权利要求1所述的图像显示方法，按照预定的间隔对CT投影数据进行整理的步骤包括：

当采样频率高于每预定间隔投影数据1次时，按照欠采样规则，将CT投影数据调整为每预定间隔投影数据1次；

当采样频率低于每预定间隔投影数据1次时，按照线性插值，将CT投影数据调整为每预定间隔投影数据1次。

5.如权利要求1所述的图像显示方法，其中，当所述基本数据为高能和低能数据时，按照物质识别算法对所述基本数据进行处理得到灰度数据和材料数据。

6.如权利要求1所述的图像显示方法，其中，将DR图像插值成为用户要求的分辨率，在屏幕上显示。

7.如权利要求1所述的图像显示方法，其中重建的三维图像的数据是以下之一：物理密度、电子密度数据、原子序数数据、高能数据、低能数据、灰度数据、材料数据。

8.一种CT系统的图像显示方法，包括步骤：

对被检查物体进行CT扫描，得到CT投影数据；

从CT投影数据重建被检查物体的三维图像；

按照某一视角对被检查物体的三维图像数据H(x,y,z)进行投影，得到该视角下的DR数据，其中三维图像数据H(x,y,z)各维的尺寸为X×Y×Z，维度x从1到X变化，该变化为水平面内垂直于皮带运动的方向，维度y从1到Y变化，该变化为竖直向上的方向，维度z从1到Z变化，该变化为水平面内沿皮带运动的方向；

基于该视角下的DR数据形成DR图像；

在屏幕上同时显示DR图像和三维图像。

9.如权利要求8所述的图像显示方法，还包括步骤：

按照另一视角对被检查物体的三维图像数据H(x,y,z)进行投影，得到该视角下的DR数据；

基于该另一视角下的DR数据形成DR图像；

在屏幕上同时显示两幅DR图像和三维图像。

10.如权利要求8所述的图像显示方法，其中所述视角和所述另一视角是从水平面内垂直于皮带运动的方向、竖直方向和水平面内沿皮带运动的方向中选择的。

11.如权利要求8所述的图像显示方法，其中将三维数据H(x,y,z)沿着维度x的方向投影，得到二维数据I(y,z)，投影公式如下：

I (y, z) = Σ_{x = 1}^{X} H (x, y, z)

该二维数据为侧视角的DR数据，数据维度尺寸为Y×Z。

12.如权利要求8所述的图像显示方法，其中将三维数据H(x,y,z)沿着维度y的方向投影，得到二维数据J(x,z)，投影公式如下：

J (x, z) = Σ_{y = 1}^{Y} H (x, y, z)

13.如权利要求8所述的图像显示方法，其中将三维数据H(x,y,z)沿着维度z即皮带运动的方向投影，得到二维数据K(x,y)，投影公式如下：

K (x, y) = Σ_{z = 1}^{Z} H (x, y, z)

该二维数据为前视角的DR数据，数据维度尺寸为X×Y。

14.如权利要求8所述的图像显示方法，其中重建的三维图像的数据是以下之一：物理密度、电子密度数据、原子序数数据、高能数据、低能数据、灰度数据、材料数据。

15.一种在DR-CT系统中进行图像配准的方法，所述DR-CT成像系统包括DR设备和CT设备，所述方法包括步骤：

用CT设备对被检查物体进行扫描以得到CT投影数据；

从所述CT投影数据重建被检查物体的三维图像；

用DR设备对被检查物体进行扫描以得到DR图像；

在屏幕上显示所述DR图像和所述三维图像；

响应于用户对DR图像中的特定位置的选择，从DR图像的数据中抽取相应的列；

从所述三维图像中得到与所述DR图像相同视角下的辅助DR数据；

计算所述DR图像中的所述列与辅助DR数据中的各列之间的相关性；

将辅助DR数据中与所述DR图像中的所述列之间相关性最大的那一列对应的切片图像与所述DR图像同时显示在屏幕上。

16.如权利要求15所述的方法，其中从所述三维图像中得到与所述DR图像相同视角下的辅助DR数据的步骤包括：

沿着维度y方向对被检查物体的三维图像数据H(x,y,z)进行投影，得到该视角下的DR数据，其中三维图像数据H(x,y,z)各维的尺寸为X×Y×Z，维度x从1到X变化，该变化为水平面内垂直于皮带运动的方向，维度y从1到Y变化，该变化为竖直向上的方向，维度z从1到Z变化，该变化为水平面内沿皮带运动的方向。

17.如权利要求16所述的方法，其中将三维数据H(x,y,z)沿着维度y的方向投影，得到二维数据J(x,z)，投影公式如下：

J (x, z) = Σ_{y = 1}^{Y} H (x, y, z)

18.一种在DR-CT系统中进行图像配准的方法，所述DR-CT成像系统包括DR设备和CT设备，所述方法包括步骤：

用CT设备对被检查物体进行扫描以得到CT投影数据；

用DR设备对被检查物体进行扫描以得到DR图像；

在屏幕上显示所述DR图像；

从所述CT投影数据中得到与所述DR图像相同视角下的辅助DR数据；

从所述CT投影数据重建被检查物体的三维图像；以及

19.如权利要求18所述的方法，其中从所述CT投影数据中得到与所述DR图像相同视角下的辅助DR数据的步骤包括：

按照预定的间隔对CT投影数据进行整理；

以一固定角度作为起始角度并以360度为间隔，从整理后的CT投影数据中抽取得到所述辅助数据。

20.如权利要求19所述的方法，还包括步骤：

将所抽取的补充数据存储在矩阵中；

对所述矩阵做左右镜像处理；

以行为单位，将镜像处理后的矩阵中的补充数据间隔插入所述辅助DR数据中，得到补充后的辅助DR数据。