[发明专利]显微CT机用GPU数字图像系统及其控制方法

说明书

技术领域

本发明属于显微CT机成像系统领域，尤其是一种加快并行图像处理速度、从而可在便携式显微CT机上使用的显微CT机用GPU数字图像系统及其控制方法。

背景技术

目前，显微CT机成像系统图像处理速度慢、使用范围受到限制。

发明内容

本发明的目的是提供一种加快并行图像处理速度、从而可在便携式显微CT机上使用的显微CT机用GPU数字图像系统及其控制方法。

本发明的技术方案是：

显微CT机用GPU数字图像系统，其特征是由为成像光源（1）、CMOS或CCD成像器（3）、嵌入式成像和光源控制系统（6）、精密机电耦合XYZ和旋转被测物体成像平台（12）、Linux Os为基础的CPU运算平台（15）组成，

其中，成像光源（1）通过RS232串行线与嵌入式成像和光源控制系统（6）连接；CMOS或CCD成像器（3）的成像聚集导轨（9）与嵌入式成像和光源控制系统（6）中的步进电机驱动单元阵列（10）连接；精密机电耦合XYZ和旋转被测物体成像平台（12）与光源控制系统（6）中的步进电机驱动单元阵列（10）连接；Linux Os为基础的CPU运算平台（15）通过路由器（14）与光源控制系统（6）连接。

显微CT机用GPU数字图像系统的控制方法，其特征是包括下列步骤：

步骤26：主机CPU启动CMOS或CCD成像器，成像器经过控制函数API获取数字图像I Inner(x,y)；

步骤27：检测该数字图像的X方向和Y方向的解析度Rx、Ry；

步骤28：假若GPU的本地内存25为MEM_local大于数字图像I Inner(x,y)的2.5倍，则执行步骤29；否则将数字图像I Inner(x,y)分割成子图像I Inner_sub(x,y)，使得本地内存MEM_local大于子图像的2.5倍，然后执行步骤29；

步骤29：将数字图像I Inner(x,y)全部读入该内存中，

步骤30：对存于TPC内存中的图像或子图像进行二维拉普拉斯-高斯算子卷积：

第一步：首先选取两维卷积用高斯核kernel，G(x,y;x0,y0,S)，如式（1），其中x0,y0设定为0,0点，S根据图像质量选取，常见选取值为[2,3]，x,y的单位为像素（pixel）：

G(x,y;x0,y0,u,S)=

1/(sqrt(2 pi)*S)exp{[(x-x0)^2+(y-y0)^2)]/(2 S*S)}…(1)

第二步：对高斯核进行拉普拉斯算子运算，如式（2）：

d²/dx²+d²/dy²[G(x,y;x0,y0,u,S)]   ...(2)

第三步：使用式（2）二维拉普拉斯-高斯算子对GPU本地内存中的数字图像进行两维卷积运算,即G(z,y;uz,uy,S)*I(z,y)，并将运算结果存储于GPU本地内存中；

步骤31：检查是否对数字图像I Inner(x,y)分割后的每个子图像都进行了读入TPC本地内存，并进行了步骤30第三步的处理，如果没有，则执行步骤29，如果是，则结束。

本发明的效果是：本发明采用基于GPU图形图像处理显卡组成的CPU-GPU系统实现便携式显微CT（Computed Tomography)机使用的并行高速二维拉普拉斯-高斯算子数字图像处理。

所述CPU-GPU系统优化了线程处理阵列（TPC，Thread Processing Clus ter）与被处理图像解析度(resolution)的优化耦合,通过嵌入式实时图像获取系统的软件驱动函数将图像解析参数与CPU-GPU算法结合，优选了图像分割阵列处理，简化了CPU与GPU之间的图像数据交换，将每个线程处理阵列拥有的N(SP)个流处理器（SP，Stream Processor）流多处理器(SM,Streaming Multiprocessor)的运算进行了并行优化处理，解决了GPU没有分支预测机制，因此在分支得到机会执行之前运算被挂起，等待所有的分支路径执行完成的运算瓶颈，加快了并行图像处理速度，从而可在便携式显微CT机上使用。实现对二维拉普拉斯-高斯算子数字图像处理的并行高速处理，加速微小物体三维断层扫描数字成像计算。

下面结合附图和实施例对本发明做进一步的说明。

附图说明

图1是本发明的结构示意图；

图2是图1中显微CT用GPU系统结构示意图；