一种由非平行断层图像序列重建任意切面的方法技术

本发明专利技术涉及一种图像重建和可视化的方法，尤其涉及一种由非平行断层图像序列重建任意切面的方法。它由原始的断层图像序列直接重建得到任意切面、即一种由非平行断层图像序列重建任意切面的方法。步骤如下：确定插值过程中的距离阈值Ｔ；选择断层图像中的感兴趣区域；确定重建切面的位置；确定断层图像中的重建区域；将断层图像上重建区域内的所有像素点映射到重建切面；计算重建切面上每个像素的灰度值；显示重建得到的切面图像。本发明专利技术方法计算量小、耗时短、精度高、实用性强，且对图像中的噪声能起到一定的抑制作用。

【技术实现步骤摘要】

【技术保护点】
一种由非平行断层图像序列重建任意切面的方法，步骤如下：　Ｓ１）确定插值过程中的距离阈值Ｔ　若重建切面上的某像素周围距离Ｔ的范围内原始断层图像上有灰度值非０的像素，则重建切面的此像素被赋值其灰度值由其周围距离Ｔ范围内的原始断层图像 上所有灰度值非０像素的灰度值平方距离加权得到；否则，重建切面上的此像素的灰度值被赋值为０；　Ｓ２）选择断层图像中的感兴趣区域　在断层图像中与重建有关的信息往往只存在于图像中心的某个区域可在图像中选择一个区域作为图像中的感兴趣区域 ，此区域将适用于断层图像序列中的所有图像，重建过程将只对断层图像中感兴趣区域内的像素进行处理；　Ｓ３）确定重建切面的位置　通过输入重建切面的坐标或者使用图形化的人机交互界面确定重建切面的位置及大小，在存储空间中为重建切面上的每个 像素开辟一个链表，并以此平面为ＸＯＹ平面建立重建空间坐标系；　Ｓ４）确定断层图像中的重建区域　若两者不平行，设每帧断层图像其和重建切面的夹角为α，则断层图像上所有到重建切面的距离小于Ｔ的像素都将位于断层图像与重建切面交线两侧宽度 为Ｗ＝Ｔ／ｓｉｎα的区域内，将此区域和断层图像中感兴趣区域的重叠区域称为重建区域；若两者平行，计算两者的距离，距离小于Ｔ则将断层图像中的感兴趣区域作为重建区域，距离大于Ｔ这帧断层图像中的重建区域为空；　Ｓ５）将断层图像上重建区域内的所 有像素点映射到重建切面　对断层图像上重建区域中的每个像素点，计算其在重建空间坐标系中的坐标（ｘ′，ｙ′，ｚ′）点，以（ｘ′，ｙ′，０）点为圆心，＊＊＊为半径在重建切面上画圆，则重建切面上圆内的所有像素点到此断层图像像素点的距离都小于Ｔ ，将断层图像上此像素点的灰度值及其到重建切面上圆内像素点的距离分别存储到重建切面上对应像素点的链表当中；　Ｓ６）计算重建切面上每个像素的灰度值　对重建切面上的每个像素点，根据存储在其对应链表中的距离小于Ｔ的断层图像上的像素点的灰 度值和两者之间的距离，平方距离加权得到其灰度值，其公式如下：　ｐ↓［ｉ］＊Ｗ↓［ｋ］＊Ｃ↓［ｋ］／＊Ｗ↓［ｋ］　其中Ｗ↓［ｋ］＝１／ｄ↓［ｋ］↑［２］，ｐ↓［ｉ］是重建切面上像素点的灰度值，Ｃ↓［ｋ］和ｄ↓［ｋ］分别是ｐ↓［ｉ］ 对应像素的链表中存储的第ｋ个灰度值和距离，ｎ为ｐ↓［ｉ］对应像素的链表中存储的灰度值的数目，即重建切面上每一个像素点的灰度值都...

【技术特征摘要】

【专利技术属性】
技术研发人员：孙丰荣，王庆浩，王文明，王丽梅，张明强，
申请(专利权)人：山东大学，
类型：发明
国别省市：88[中国|济南]