动态影像序列分割方法
专利汇可以提供动态影像序列分割方法专利检索，专利查询，专利分析的服务。并且一种动态影像序列分割方法，属于 图像处理 领域。方法如下：确定处理对象的边界，得到区分处理对象与图像背景的模板；将动态序列中所有图像的背景分割掉，避免背景因素干扰；构造处理对象区域中所有 像素 点的灰度平均值序列为基准时间强度曲线；利用线性相关系数构造动态序列的综合特征图像；分割综合特征图像得到分割结果图像；从分割结果图像选择感兴趣区域；对分割结果图像中的感兴趣区域的每一 块 划分，以动态序列的第一 帧 图像灰度信息进行填充，然后，该划分的范围内分割该部分图像，得到该划分内的加细分割结果；将加细分割结果进行拼装，得到最终结果图像。本 发明 充分利用动态影像序列中的有效信息，从而实现对动态影像序列准确分割的目的。，下面是动态影像序列分割方法专利的具体信息内容。
1、一种动态影像序列分割方法，其特征在于，方法步骤如下：
步骤一：利用动态序列的第一帧图像的灰度信息，确定处理对象的边界，从 而得到区分处理对象与图像背景的模板；
步骤二：根据步骤一确定的处理对象与背景的位置模板，将动态序列中所有 图像的背景分割掉，避免背景因素在下面的步骤中的干扰；
步骤三：构造步骤二中所确定的处理对象区域中所有像素点的灰度平均值序 列为基准时间强度曲线；
步骤四：利用线性相关系数构造动态序列的综合特征图像；
步骤五：分割步骤四中的所得综合特征图像得到分割结果图像；
步骤六：从步骤五得到的分割结果图像选择感兴趣区域；
步骤七：对分割结果图像中的感兴趣区域的每一块划分，以动态序列的第一 帧图像灰度信息进行填充，然后，该划分的范围内分割该部分图像，得到该划分 内的加细分割结果；
步骤八：将步骤七所得的所有划分的加细分割结果按照步骤六中所在区域的 位置进行拼装，得到最终结果图像。
2、根据权利要求1所述的动态影像序列分割方法，其特征是，所述的步骤四， 利用线性相关系数构造动态序列的综合特征图像具体为：
(1)构造处理对象区域中所有像素点的时间强度曲线与步骤三中构造的基准 时间强度曲线之间的线性相关系数矩阵；
(2)由(1)所得的线性相关系数矩阵构造综合特征图像，该图像综合了动 态序列中有效位置的所有时间强度曲线的特征。
3、根据权利要求1所述的动态影像序列分割方法，其特征是，所述的步骤五， 分割步骤四中的所得综合特征图像得到分割结果图像具体为：
(1)对分割结果图像进行对比度增强操作，获得对应的对比度增强图像
(2)利用分水岭算法分割(1)中所得的对比度增强图像。
4、根据权利要求1所述的动态影像序列分割方法，其特征是，步骤六中，兴 趣区域默认为处理对象所在的整个区域，或者是实际选择的感兴趣区域。
技术领域
本发明涉及一种动态影像序列分割方法，特别是一种利用线性相关系数对动 态影像序列进行分割的方法。属于图像处理领域。
技术背景
动态影像序列是指对于同一功能位置连续采集一段时间段内不同时刻的影像 构成的序列，如核磁共振设备采集的动态增强序列和血管灌注序列。动态影像序 列反映了处理对象在时间段内的影像变化情况。动态序列的基本构成单位是每一 像素点在不同时刻的灰度值构成的时间强度变化曲线。因此，对动态影像序列的 分割有别于对普通图像的分割。动态影像序列的最后分割结果应该不仅反映处理 对象的影像灰度信息，也要能够综合时间反映动态影像序列中全部时间强度变化 曲线中隐藏的功能信息。
经文献检索发现，J.Rogowska，K.Preston，and G.Hunter.在IEEE Transactions on Medical Imaging，(《IEEE医学影像汇刊》)14(3)：480-486， September 1995上发表的“Applications of similarity mapping in dynamic mri” (相似性映射在动态MRI序列中的应用)一文，该文提出利用线性相关系数对动 态MRI序列进行分析，即计算动态影像序列相关感兴趣区域内各点的时间强度变 化曲线与基准时间强度变化曲线之间的线性相关系数，从而产生该区域所有位置 的相关系数值构成的结果图像。但是该项技术有以下两个缺点：(1)线性相关 系数计算结果可用于衡量两支曲线的线性相似程度，但在相关系数计算结果值为 极小的时候，比如，计算A、B时间强度曲线与基准曲线所得线性相关系数的结 果数值均小，表示A、B时间强度曲线均与基准曲线的线性相似度低，但不能说 明A、B时间强度曲线之间的线性相似程度就高。而该方法的结果图像中这种差 异性被忽略了。(2)结果图像综合反映了动态影像序列的时间强度曲线特征，但 忽略了该动态影像序列的原始影像特征，造成有效信息的缺失。
发明内容
本发明是通过以下技术方案实现的，方法步骤如下：
步骤一：利用动态影像序列的第一帧图像的灰度信息，确定处理对象的边界， 从而得到区分处理对象与图像背景的模板。
步骤二：根据步骤一确定的处理对象与背景位置的模板，去掉动态影像序列 中各帧图像的背景部分，这样可以避免背景因素在以下步骤中可能产生的干扰。
步骤三：由步骤二中所确定的处理对象区域中所有像素点的灰度平均值序列 构造基准时间强度曲线。
步骤四：以步骤三所得到的基准时间强度曲线，构造由线性相关系数表现的 动态影像序列的综合特征图像。
步骤五：分割步骤四中的所得综合特征图像得到分割结果图像。
步骤六：从步骤五得到的分割结果图像选择特定的划分组合构成感兴趣区域。
步骤七：对分割结果图像中的感兴趣区域的每一块划分，以动态序列的第一 帧图像对应区域的灰度信息进行填充。然后，在该划分的范围内分割图像，得到 该划分内的加细分割结果。
步骤八：将步骤七所得的所有划分的加细分割结果按照步骤六中的区域位置 进行拼装，得到最终结果图像。
以下对本发明作进一步的说明，具体内容如下：
所述的步骤四，具体为：
(1)构造研究对象区域中所有像素点的时间强度曲线与步骤三中构造的基准 时间强度曲线之间的线性相关系数矩阵。
(2)由(1)所得的线性相关系数矩阵构造综合特征图像。该图像综合了动 态序列中有效位置的所有时间强度曲线的特征。
所述的步骤五，具体为：
(1)对分割结果图像进行对比度增强操作，获得对应的对比度增强图像。
(2)利用分水岭算法分割(1)中所得的对比度增强图像。
本发明方法主要具有以下特点和作用：(1)提供了综合表现时间强度变化曲 线的手段，其最终结果综合了动态影像序列中隐藏的功能信息。(2)不仅综合了 时间强度变化曲线的特征，而且综合了动态序列中原始采集的图像灰度信息，其 最终结果是全面利用动态影像序列信息所得的结果。(3)是对动态影像序列的全 面的半自动分割方法。本发明适用于对处理对象按时间段进行连续图像采集以反 映处理对象功能的场合，如核磁共振设备采集的动态增强序列和血管灌注序列。
具体实施方式
核磁共振设备采集的血管灌注影像序列是动态影像序列的一种。血管灌注影 像序列中像素点的时间强度曲线(即血管灌注曲线)反映了对应生理位置的血管 灌注功能。实际应用中，医生通过选择生理位置，观察对应位置的血管灌注曲线 的变化情况，确定该位置的生理特性；但目前缺乏全面的自动的反映综合所有血 管灌注曲线特征和血管灌注序列影像灰度特征的方法。本发明可以应用在这种场 合。
血管灌注序列病例影像为某某医院利用马可尼医疗设备公司1.5T核磁共振影 像设备(Eclipse 1.5T，Marconi Medical Systems，Inc，)采集的17岁男孩的腿部 (含骨肉瘤)的血管灌注序列。
运行本发明方法的设备为普通个人计算机，其具体配置为P41.8G，256M内存， 80G硬盘。发明方法的具体实现平台为Windows2000下的Matlab计算平台。
具体实施情况如下：
步骤1、确定骨肉瘤血管灌注影像序列的背景：
(1)利用Sobel算子对血管灌注序列第一幅图像进行边界检测。
(2)利用形态学的膨胀操作完成边界连接为封闭边界。
(3)去掉边界内部的孔隙，以及与边界连接的可能物体。
(4)利用形态学的腐蚀操作进行边界平滑，得到了区分背景与生理位置的 二值图。对应生理位置区域设为ROI(处理对象区域)。
步骤2利用步骤1(4)中的二值图为模板，应用到整个血管灌注序列中，去 除血管灌注序列中所有图像的背景。
步骤3计算基准时间强度曲线：
利用去除背景后的整个生理结构区域ROI的影像中所有像素点的灰度平均值 构成的时间强度曲线作为基准时间强度曲线，计算公式为：
①
其中，
Aijt 为血管灌注影像序列A1、A2......AT中第t幅图像At中的
(i，j)位置对像素的灰度值。
(i，j)∈ROI；t＝1，2，......T；
ROI 为步骤1.4中得到的生理位置区域。
T 为血管灌注影像序列中图像的总帧数；
μAt 为血管灌注影像序列中第t幅图像At的ROI区域中的所有像
素的灰度平均值。
步骤4、构造综合特征图像：
(1)计算影像中所有有效像素点位置的标准化线性相关系数，构造线性相关 系数矩阵。
标准化计算线性相关系数的公式为
②
其中，
NCOR11 为所求标准化线性相关系数。
μA 为序列μA1，μA2，......μAT的平均值。
其余参数同公式①中的说明。
(2)由线性相关系数矩阵转换为每一像素点灰度值的构造综合特征图像。
由于标准化线性相关系数的取值范围为(-1，1)，为获得综合特征图像，需 要将标准化计算线性相关系数映射成为可显示图像灰度范围。映射公式为：
其中，
[k1，k2]-为步骤4(1)中实际计算所得到的标准化线性相关系数的取 值范围。
k-为对应像素点的标准化线性相关系数值。
I(k)-为利用公式③转换所得的对应像素点所得灰度值。
步骤5、利用改进的分水岭算法分割综合特征图像，得到分割结果图像：
(1)对分割结果图像进行对比度增强操作，获得对应的对比度增强图像
①以步骤2所得的中间结果图像为本步骤的原始图像。计算形态学的“头帽” 图像和“底帽”图像。
②计算对比度增强图像，计算公式为：
对比度增强图像＝-(原始图像+头帽图像-底帽图像)
(2)利用分水岭算法分割对比度增强图像。
所得分割结果的每一分块为与基准时间强度曲线的线性相似度相近的区域。 由于血管灌注曲线的功能特性，每一区域反映的血供情况基本相同。步骤5(1) 及步骤5(2)中组成了改进的分水岭算法。
步骤6、由临床影像学医生从步骤5中的分割结果图像选取所有可能包含肿 瘤病变区的分块，共同组成待进一步研究的敏感区域。
步骤7、利用该血管灌注序列的第一幅图像灰度信息填充步骤6步中所得的 整个敏感区域的每一划分。然后，对每一划分所对应影像区域继续使用改进的分 水岭算法进行加细分割。
步骤8、将第7步所得所有加细分割结果按原对应位置拼接起来，就可以得 到对敏感区域内部的加细分割结果，为最终结果图像。
本发明方法取得了较好的实施效果：最终结果图像显示出了对敏感区域内部 的分割结果，该分割结果综合了血管灌注序列反映的生理组织的血供功能信息和 该对应位置的原始采集的影像灰度信息。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 岩石孔隙孔径的测定方法及设备 | 2020-05-08 | 745 |
| 基于面线基元关联模型的光伏面板模板自动生成方法 | 2020-05-12 | 219 |
| 基于动态分配的分水岭图像分割并行方法 | 2020-05-12 | 212 |
| 一种基于滑动窗口算法的纹理特征度量方法 | 2020-05-18 | 222 |
| 一种基于颗粒粗糙度的纹理特征度量方法 | 2020-05-18 | 94 |
| 一种知识驱动的高空间分辨率遥感影像自动变化检测方法 | 2020-05-13 | 231 |
| 一种基于OCT图像处理的枇杷早期瘀伤鉴别方法 | 2020-05-14 | 926 |
| 一种基于方向超像素的快速图像分割方法 | 2020-05-11 | 861 |
| 一种基于连通区域分析的肺部图像分割方法及系统 | 2020-05-08 | 616 |
| 一种基于对抗学习的非感应式餐盘图像数据自动标注方法 | 2020-05-16 | 910 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
