[发明专利]一种CT扫描图像旋转误差校正方法

权利要求书

1.一种CT扫描图像旋转误差校正方法，其特征在于包括以下步骤：

（1）选取不同驱替时刻位于相同扫描位置的两张CT扫描图像；

（2）选取阈值t₁、t₂，，将图像进行二值化，选取图像中灰度值介于t₁和t₂之间的像素在二值图中赋值为1，其余像素赋值为0，选取岩石基质作为后续校正方法的参照物；

（3）将得到的二值图像由直角坐标图像转换为极坐标图像，转换公式为：

θ=arctan(yx)r=x2+y2]]>

其中，(x,y)为直角坐标系坐标，(θ,r)为极坐标系坐标，图像灰度值计算采用最邻近插值算法；

（4）统计极坐标图像中每个岩石基质所占像素的坐标，计算该岩石基质的质心坐标，计算公式为：

θ′=Σi=1nθi/nr′=Σi=1nri/n]]>

其中，n为某一岩石基质区域内的像素点数目，θ_i、r_i为第i个像素点的横坐标和纵坐标；

（5）对比不同驱替时刻极坐标图像中相同岩石基质质心横坐标的平移量，反求两张原始CT扫描图像的相对旋转角度，所述横坐标的平移量计算公式为：

Δθj=θj1-θj2,θj1-θj2≥0θj1-θj2+2πθj1-θj2<0,]]>

Δθ=Σj=1mΔθj/m,]]>

其中，m为极坐标图像中岩石基质的数目，θ_j1为第1张极坐标图像中第j个岩石基质的横坐标，θ_j2为第2张极坐标图像中第j个岩石基质的横坐标，△θ_j为不同驱替时刻极坐标图像中第j个岩石基质的质心横坐标平移量，△θ为不同驱替时刻两张极坐标图像的相对平移量，即不同驱替时刻两张CT扫描图像的相对旋转角度；

（6）根据前面反求的旋转角度，对原始CT扫描图像进行旋转操作，旋转公式为：

x′=xcosα-ysinαy′=xsinα+ycosα]]>

其中，(x′,y′)为(x,y)进行旋转变换的结果，图像灰度值计算采用双线性插值算法。

2.根据权利要求1所述的CT扫描图像旋转误差校正方法，其特征在于，所选取的阈值t₁<t₂。