[发明专利]一种最大密度投影图像数据后处理方法

说明书

技术领域

本发明涉及图像处理技术领域，具体涉及一种最大密度投影图像数据后处理方法。

背景技术

在图像处理技术中，三维图像的建立可以通过对二维图像序列的投影完成，最大密度投影(MIP，Maximum intensity projection)作为一种三维图像的重建技术，一直是MR(Magnetic resonance imaging，磁共振)、(computed tomography，计算机断层成像术)数据体得到血管、骨骼和软组织等三维结构的有效图像重建方法。为了利用计算机快速进行图像处理，满足图像重建的时间要求，人们研究、发展了很多最大密度投影MIP算法，同时由于在生成MR、CT的原始二维图像序列时，是通过数据采样得到二维图像数据的，这些图像数据不可能完全满足奈奎斯特采样定律，因此经过最大密度投影后不可能完全重现三维物体的原貌，这就需要对二维图像数据在最大密度投影前进行插值处理(称为前插值)，或者对得到的最大密度投影三维图像数据进行投影后插值处理(称为后插值)，以弥补原始图像数据的不足。现有技术的前插值和后插值方法很多，前插值方法一般计算量是很大的，加上三维重建的数据量本身是“海量”的，常常直接造成效率的很低，重建速度让人难以接受。早期的CT、MR的机器，因扫描间隔很大，要进行三维重建，保持采样一致性一般要进行很大计算量的前插值处理。但现在多排CT，强场MR所获取的图像序列的扫描间隔可以很小，只要最大密度投影重建图像方法合理，就可以直接重建了。后插值方法虽然效率很高，但现有技术的后插值方法仅仅针对最大密度投影三维图像数据进行插值处理，而没有考虑最大密度投影前二维图像序列各图像切片之间的空间关系，造成最终得到三维图像的组织纹理细节模糊，有的甚至变形。

发明内容

本发明要解决的技术问题是提供一种最大密度投影图像数据后处理方法，克服现有技术的后插值方法得到三维图像的组织纹理细节模糊，甚至变形的缺陷。

本发明为解决上述技术问题所采用的技术方案为：

一种最大密度投影图像数据后处理方法，包括步骤：

A1、在对图像序列进行最大密度投影前，选取各幅图像的同一位置的顶点，沿所述图像序列递增或递减的顺序计算并记录相邻两幅图像的所述顶点在投影平面上的坐标偏移量；

A2、在对所述图像序列进行最大密度投影后，记录投影数据中每个像素点的坐标、灰度值和所在图像的标识及其关系，对所述像素点进行后处理，包括步骤：

B1、扫描所述投影数据中的像素点，根据所述像素点的坐标和所述坐标偏移量计算所述像素点所在图像的相邻图像在所述投影数据中的投影匹配像素点的坐标，进而找到所述投影匹配像素点的灰度值；

B2、计算所述像素点和所述投影匹配像素点之间的距离，若所述距离大于1，则执行步骤B3，否则针对下一所述像素点，返回步骤B1；

B3、计算所述像素点和所述投影匹配像素点之间的像素点的坐标，进而找到其灰度值；

B4、根据公式：

G′des+i=Gsr×n-in+Gdes×in,]]>

计算所述像素点和所述投影匹配像素点之间的像素点的插值灰度值，其中G_sr是所述像素点的灰度值，G_des是所述投影匹配像素点的灰度值，n是所述像素点和所述投影匹配像素点之间的像素点的个数，i是计算时当前像素点的序数，取值范围是[1，n]，G′_dest+i是当前像素点计算得到的所述插值灰度值；

B5、若计算得到的所述插值灰度值大于所述像素点和所述投影匹配像素点之间的像素点的原有灰度值，则用所述插值灰度值替换所述原有灰度值。