[发明专利]一种三维医学图像多尺度异构特征的提取与分类方法

说明书

技术领域

本发明涉及三维医学图像多尺度异构特征的提取与分类方法，具体涉及三维图像曲线、曲面、实体块等结构特征的自动提取与语义分类方法，通过将张量场构建、各向异性小波核函数定义、多尺度分析和自适应聚类进行有机结合，设计了一个统一的算法来实现这些多尺度、多级别异构特征的自动提取。

背景技术

三维医学图像多尺度异构特征的提取与分类方法相关方面为热核理论应用、图像特征提取和基于热扩散的图像处理三个方面。

在热核理论应用方面：

近年来，物理上的热核和热扩散理论在计算机图形图像处理领域受到了众多研究者的青睐，已在二维流形的多尺度特征提取、平滑处理以及形状检索等方面得到了广泛应用。其中，最重要的原因是：热核天然地蕴含了多尺度的性质，并具有等距变换不变性，这对图形或图像的鲁棒特征提取具有重要意义。如：文献提出了一种尺度不变的热核信号计算方法，并将其应用于形状检索，取得了较好的效果。同时，为了提高在大规模网格模型上进行热核计算的时间效率，文献提出了一种多分辨率的计算策略来加速热核信号的计算。

在图像的特征提取方面：

图像的局部特征可被看作是与周围领域具有显著不同、并表征特定物理意义的一种模态，他们一般蕴含丰富的信息并可通过一些微分量来进行验证。其中，现有的以Harris、Hessian、LOG以及DOG等算子为核心的图像特征提取算法的主要思想是：通过对图像的微分属性分析，提取某种物理量的局部极值，并将其作为特征。因此，这些算法一般对尺度变化、噪声和图像形变等因素都比较敏感。为了尽可能的解决这些问题，文献的实验证明：将Harris算子和Hessian算子相结合并辅以尺度选择，是提高图像特征提取鲁棒性的一个可行的策略，而且该方法已被文献推广到了三维情形。应该说，目前2D SIFT算法仍是最为稳定的二维图像特征提取算法，而且它的高维形式的扩展已经开始出现.

在基于热扩散理论的图像处理方面：

在过去的20年中，借助高斯核函数，热扩散理论在图像的多尺度分析领域被广为使用。但是，高斯核一般只能一个很小的邻域内刻画热扩散过程，因为其本质上等价于用各向同性的热传递系数作为权重对原始图像进行卷积处理。自从双边滤波方法被提出来用于边界保持的图像平滑和去噪处理以后，基于各向异性热传递理论的滤波器才开始在二维图像处理领域逐渐流行起来。特别值得指出的是，在二维流形上首先提出的热核信号已被推广用于三维流形，但是这种扩展尚没有考虑三维流形内部材质的差异性。这也是启发本发明直接在三维医学图像上研究各向异性热核的一个主要因素。

发明内容

DTI图像所包含的张量信息对器官结构特征的提取和可视化均具有重要意义，而临床常用的MRI,CT等标量医学图像并不自带张量信息。因此，本部分内容重点研究如何根据三维标量图像的几何信息，来构建描述局部结合结构的张量场。在此基础上，进一步研究如何将二维SIFT特征提取的思想进行扩展，实现三维医学影像数据多尺度几何特征的提取，并研究适用于三维医学数据非刚性等距形变的特征描述子，为配准等其他基于特征的三维医学图像应用提供底层的算法支持。

本发明提出一种三维医学图像多尺度异构特征的提取与分类方法,设计了一种鲁棒的三维图像异构结构特征自动提取、分类、过滤和可视化方法，该方法可同时处理多尺度点特征、多尺度曲面特征、多级别曲线特征和实体块特征，并在异构结构特征提取的基础上，给出了一种自适应的结构特征语义分类方法；基于三维图像局部黑塞矩阵的特征系统定义了一种通用的扩散张量，该张量能够解释图像的局部几何结构，辅助曲线种子点和实体块特征的提取，并可基于扩散张量的距离度量控制线特征的自动生长；基于扩散张量定义了一种数据相关的各向异性小波核函数，在图像的多尺度分析过程中，它可用于图像的各向异性卷积计算和图像的频带分解；基于点到局部拟合曲面的投影距离，提出了一种曲线候选种子点过滤方法，可有效去除由边界区域离散体素的锯齿而导致的伪种子点，同时，给出了基于CUDA的曲线并行生长算法。

附图说明

图1人体头部MRI三维图像异构特征的提取和分类结果。

图2为基于扩散张量和各向异性小波多尺度分析的三维图像异构特征分类算法流程图。

图3为结构张量、扩散张量及其物理意义的示意图。

图4为多尺度曲面特征的提取结果。

图5为多尺度曲面特征的语义分类。

图6为基于MLS局部曲面拟合的种子点过滤方法示意图。

图7为10%噪声扰动的人体头部MRI三维图像的多尺度异构结构特征分类结果。