[发明专利]一种CT成像中投影图背景不一致的校正方法

摘要

本发明涉及一种CT成像中投影图背景不一致的校正方法，具体是一种基于图像灰度信息的CT投影图背景互相关配准法，包括如下步骤：1)在采集样品CT投影图像之前，采集亮场图像一张。确定亮场图像的光斑中心；2)计算亮场图像补偿函数；3)亮场图像照度补偿；4)投影图像照度补偿；5)投影图背景与亮场图背景的实际偏移量计算；6)图像平移配准。本发明基于图像灰度信息，利用互相关原理实现对CT投影图的配准，具有不需要预先知道探测器具体偏转参数，操作简单，计算效率高的优点。

主权项

1.一种CT成像中投影图背景不一致的校正方法，其特征在于：包括以下步骤：1)在采集样品CT投影图像之前，采集一张亮场图像，确定亮场图像的光斑中心O(i0，j0)；所述确定亮场图像的光斑中心O(i0，j0)的具体操作步骤为：将亮场图像每点的灰度数值读出，表示为一个灰度数值矩阵{pi,j}(i＝0,1,2...,imax；j＝0,1,2,...,jmax)，其中i、j分别表示矩阵{pi,j}的列坐标和行坐标,与亮场图像上各像素点坐标一一对应；imax是图像列方向像素点的最大坐标值，jmax是图像行方向像素点的最大坐标；pi,j表示坐标为(i,j)处的灰度数值；从亮场图像中光斑边界清晰的图像区域选择第j1,j2,...,jq行，共q行的灰度数据，对每一行灰度数据求其对列坐标i的差分f′j(i)：其中Δi取自然数，其数值要在保证光斑边界处的f′_j(i)能明显区别于图像其他位置的前提下取最小值；当差分函数取到最小值f′_j(i)_min和最大值f′_j(i)_max时，对应自变量分别记为和根据公式(2)计算得到光斑在亮场图像第j行的中心坐标其中表示取数值的整数部分；计算出不同j值时的算术平均值作为亮场图像在行方向的光斑中心坐标i₀，即：其中,表示取数值的整数部分；采用同样的办法得到亮场图像列方向光斑中心坐标j0，则亮场图像光斑中心O点坐标确定为O(i0，j0)；2)计算亮场图像补偿函数f(r)；所述亮场图像补偿函数的计算方法为：以光斑中心O(i₀，j₀)为圆心，半径为(r+0.5)作圆，其中r＝[1,2,...,min(i₀,j₀)]，单位为一个像素，min(i₀,j₀)表示取i₀,j₀这两个数值中较小的一个数值，r取不同值时，分别求出圆周上所有像素点灰度值的算术平均值，表示半径为(r+0.5)的圆周上所有像素点灰度值的算术平均值；根据不同r值得到的用多项式拟合法得到关于r的函数,表示为f(r)，即为亮场图像补偿函数；3)亮场图像照度补偿；所述亮场图像照度补偿的具体操作步骤为：根据亮场图像补偿函数f(r)对亮场图像中所有像素点进行照度补偿：亮场图像灰度数值矩阵{pi,j}中，坐标为(i,j)的像素点与光源中心O(i0,j0)的距离表示为R，则：将r＝R代入亮场图像补偿函数f(r)中，计算得到亮场图像补偿后在坐标(i,j)处像素点的灰度值ai,j：对亮场图像中所有像素点执行公式(5)的计算，得到补偿后的亮场图像灰度数值矩阵{ai,j}，根据{ai,j}利用计算机编程可输出补偿后的亮场图像；4)投影图像照度补偿；所述投影图像照度补偿的具体操作步骤为：将每张投影图像每点的灰度数值读出，表示为一个灰度数值矩阵其中m、n分别表示矩阵的列坐标和行坐标,与投影图像上各像素点坐标一一对应,m_max是图像列方向像素点的最大坐标值，n_max是图像行方向像素点的最大坐标值，k表示投影图的编号，k_max是投影图总的张数；表示第k张投影图上坐标为(m,n)处的灰度数值；确定每张投影图的光斑中心坐标分别表示第k张投影图光斑中心的行坐标和列坐标；利用亮场图像补偿函数f(r)对每张投影图的灰度数值进行补偿，表示补偿后的第k张投影图的灰度数值矩阵，表示照度补偿后的第k张投影图在列坐标为m，行坐标为n处的灰度数值；根据利用计算机编程可输出照度补偿后的投影图像；5)投影图背景与亮场图背景的实际偏移量计算；所述投影图背景与亮场图背景的实际偏移量计算方法为：以照度补偿后的亮场图像作为模板图像，照度补偿后的第k张投影图像作为待配准图像，在模板图像和待配准图像中的相同位置处选取大小相同的区域1和区域2，这两个区域所包含像素点的灰度值分别作为两组离散序列ξ(e,f)和χk(e,f)，其中e,f分别表示所选图像区域的列坐标和行坐标，e＝0,1,2,...,emax，f＝0,1,2,...,fmax，emax和fmax分别表示所选图像区域中列方向和行方向的最大坐标值，其数值根据所选图像区域大小确定；根据公式(11)计算互相关函数ck(u,v)的数值，公式(11)中，u、v分别表示亮场图像和投影图像在列方向和行方向可能的背景偏移量，u＝0,1,2，...，u_max,v＝0,1,2,...，v_max,其中u_max与v_max分别表示投影图与亮场图像在列方向和行方向可能的最大背景偏移量，其数值根据X射线CT实验中X射线源、样品台、探测器三者之间可能发生的最大偏移量确定，代表序列ξ(e,f)的算术平均值，代表序列χ^k(e,f)的算术平均值；当c^k(u,v)取到最大值时，对应的背景偏移量(分别记为)即为亮场图像和第k张投影图像的实际背景偏移量；采用相同的方法将所有投影图像背景配准，从而得到每张投影图像背景与亮场图像背景的实际偏移量；6)图像平移配准；所述图像平移配准的具体操作步骤为：根据照度补偿后每张投影图背景与亮场图像背景的实际偏移量可以将所有投影图像进行平移变换，实现CT投影图像背景与亮场图像背景的配准；根据背景配准后的第k张投影图灰度数值矩阵利用计算机编程可输出背景配准后的投影图。