本构模型与有限元结合的脉络膜新生血管生长预测方法技术

本发明专利技术公开了一种本构模型与有限元结合的脉络膜新生血管生长预测方法，图像预处理；区域提取和划分：将图像分割为CNV区域、外视网膜层、内视网膜层和脉络膜层4个区域；网格化：对4个区域进行四面体网格生成；建模：运用超弹性生物力学模型与反应扩散方程建模，将脉络膜新生血管生长后的质量改变作为源项加入到方程中，使变形梯度张量根据新生血管的生长持续变化；优化模型，算出最佳准确率，进行参数检验；根据每个时间点预测的参数拟合一条参数曲线，预测最后一个时间点的生长参数，得到预测结果。本发明专利技术的方法能够更加灵活和个性化的模式生物机械模型，模型假设的组织是正交各向异性，对非线性、大变形区域也有很好的预测结果，准确度高。

【技术实现步骤摘要】
本构模型与有限元结合的脉络膜新生血管生长预测方法
本专利技术涉及计算机视觉、图像处理与分析
，属于生长预测的建模方法，尤其是应用于SD-OCT(频域光学相干断层成像)的脉络膜新生血管(CNV)生长预测的建模方法。
技术介绍
脉络膜位于视网膜和巩膜之间，由纤维组织、小血管和毛细血管组成。脉络膜的血循环营养视网膜外层、晶状体和玻璃体等，由于脉络膜内血流量大、流速较慢、而其本身软而薄的特性，使其对外界影响或冲击不耐受，导致病原体在此处易滞留，造成脉络膜疾病及各种眼底疾病，因此，研究脉络膜新生血管(CNV)的生长可以得到很多重要的数据，为后续医学研究提供基本的支撑数据。由于新生血管的生长部位在对外界环境极其敏感和脆弱的眼部视网膜内，实时入侵性的对新生血管进行检测显然是不可能的。现有的新生血管生长预测方法通常仅采用单一的增长模型如有限元法进行预测，这种方法存在一定的缺点和局限性：①基于各向同性的假设，②没有考虑到不同组织结构的不同生物学特性，③只能预测线性变化，④只能预测小变形的变化，对大变形变化的预测不准确，等等。
技术实现思路
本专利技术所要解决的技术问题是克服现有技术的缺陷，本专利技术首次提出了一个综合性的方法——基于超弹性材料的生物力学模型和反应扩散方程相结合的方法，进行脉络膜新生血管的生长预测分析，准确度高。为解决上述技术问题，本专利技术提供一种脉络膜新生血管三维生成预测方法。该方法主要包括6个步骤：步骤1.图像预处理：对三维视网膜频域光学相干断层成像SD-OCT图像进行下采样处理；步骤2.区域提取和划分：将图像分割为脉络膜新生血管CNV区域、外视网膜层、内视网膜层和脉络膜层4个区域；步骤3.网格化：对CNV区域、内视网膜层、外视网膜层和脉络膜层进行多尺度、自适应的四面体网格生成；步骤4.运用超弹性生物力学模型与反应扩散方程建模：在超弹性生物力学模型结构中，整合各区域不同的体积模量、剪切模量及其对应的体积弹性响应、等容弹性响应，将脉络膜新生血管生长后的质量改变作为源项加入到方程中，使变形梯度张量根据新生血管的生长持续变化；步骤5.优化模型，算出最佳准确率，进行参数检验；步骤6.根据每个时间点预测的参数拟合一条参数曲线，预测最后一个时间点的生长增值率参数，得到预测结果。步骤1中的图像预处理包括以下步骤：(1a)下采样处理：(1b)使用CAVASS软件，进行OCT图像配准。步骤3中采用一个三维曲面和体积的网格生成器对CNV区域、内视网膜层、外视网膜层和脉络膜层进行多尺度、自适应的四面体网格生成。步骤4包括以下步骤：(4a)反应扩散方程的入侵和扩散逻辑回归型增长的反应扩散方程表示为：其中，t为时间，D是各向异性的扩散张量，包含三个组成部分Dx，Dy，Dz，这三个分量分别为x轴，y轴，z轴的分量；ρ代表生长增殖率，N是每单位体积的细胞数，Kmax代表最大细胞容量，当血管内皮生长因子的各向异性浓度时，有如下公式：(4b)质量效应与超弹性生物力学模型材料的非线性应力-应变行为的应变能密度函数为：其中，表示格林-拉格朗日应变张量，FT代表变形梯度张量F的转置张量，I为单位矩阵；J代表变形梯度张量F的行列式，即J＝det(F)；κ和μ分别代表体积模量和剪切模量；Tr表示矩阵的迹；I1是柯西-格林变形张量的第一主不变量；通过公式4中的梯度力f，将柯西应力张量σ＝J-1FSFT和血管内皮生长因子的各向异性浓度c联系起来：其中，f代表由c产生的梯度力即密度梯度力，ξ代表一个需要被估计的常数，该常数取决于材料的生物学特性；S代表PK2-第二Piola-Kirchhoff应力张量。步骤5中，以一阶马尔科夫序列优化生长增值率ρ参数：前n个时间点的数据，每个时间点均与之前的一个时间点相关，计算ρ的值。从之前所有的时间点进行生长增值率参数的估计，然后使用下面的公式来计算预测的准确性：是基于之前所有时间点Ii,c(i＝1,2,...，10)预测得到的结果，Ii+1,c是当前时间点的金标准；最后，估计参数：νi＝{Dx,Dy,Dz,ρ；c}。定义预测结果与金标准的重合像素点计算反映预测结果准确性的各项指标：(1)召回率(2)错误率(3)相似度DICE系数步骤6中，利用得出的前n个时间点的生长增值率ρ值，画出拟合曲线，从而预测出第n+1个时间点的生长增值率ρ；将预测出的第n+1个时间点的生长增值率ρ代入第n个时间点，得出的预测图像与真实的第n+1个图像重叠，计算预测准确率。本专利技术所达到的有益效果：本专利技术首次提出了一个综合性的生长预测方法——基于超弹性材料的生物力学模型和反应扩散方程相结合的方法，进行脉络膜新生血管的生长预测分析。本专利技术的方法结合超弹性材料和反应反应扩散方程，能够更加灵活和个性化的模式生物机械模型。不同于以往的材料是各向同性假设，该模型假设的组织是正交各向异性。在疾病部位(脉络膜新生血管区域)和周围的组织(包括外层视网膜和脉络膜，视网膜内层等)模拟出不同组织结构的不同生物学特性，且基于正交各向异性的假设，对非线性、大变形区域也有很好的预测结果，准确度高。附图说明图1是脉络膜新生血管预测的整体框图；图2a是OCT的原始图像(512×1024×128)网格化后的结果；图2b是下采样的图像(64×64×64)网格化后的结果；图2c表示原始图像的节点信息和单元信息；图2d表示下采样图像的节点信息和单元信息；图3a是从下采样数据重构的网格划分；图3b是脉络膜新生血管的细胞内体积(ICV)；图3c是ICV的等值面图像；图3d是应力分布的显示：每一个切片上，CNV区域的粗线条圆形和细线条椭圆形分别代表强应力(一般分布在CNV边缘)和弱应力(一般分布在CNV内部)的位置；图4是不同病人的CNV体积差异示意图。具体实施方式下面对本专利技术的具体实施方式作进一步详细的描述：本方法的基本框图如图1所示，主要包括以下步骤：1)图像预处理(a)下采样处理：为了提高处理效率，优化内存使用，本实施例中将每个像素512*1024*128的三维视网膜SD-OCT图像下采样到像素64*64*64。(b)使用CAVASS软件(MIPG开发的软件)，进行OCT图像配准。2)区域提取和划分手动标记金标准，进行视网膜感兴趣区域(CNV区域)的提取、周围组织(内视网膜层、外视网膜层、脉络膜层)的划分；3)网格化，用一个三维曲面和体积的网格生成器对CNV区域、内视网膜层、外视网膜层和脉络膜层进行多尺度、自适应的四面体网格生成，见图2a、图2b、图2c、图2d；4)对CNV区域、内视网膜层、外视网膜层和脉络膜层这4个区域建模，其中各个区域进行不同的超弹性生物力学模型与反应扩散方程模型进行参数训练，使其根据不同的生理学特性进行生长：(a)反应扩散方程的入侵和扩散逻辑回归型(Logistic)增长的反应扩散方程可以表示为：其中，t为时间，D是各向异性的扩散张量，包含三个组成部分Dx，Dy，Dz，这三个分量分别为x轴，y轴，z轴的分量。ρ代表生长增殖率，N是每单位体积的细胞数，Kmax代表最大细胞容量，当血管内皮生长因子(VEGF)的各向异性浓度时，有如下公式：(b)质量效应与超弹性生物力学模型在模型结构中，通过方程(3)比较全面地整合各区域不同的体积模量、剪切模量及其在本构方程中对本文档来自技高网...

【技术保护点】
一种本构模型与有限元结合的脉络膜新生血管生长预测方法，其特征是，包括以下步骤：步骤1.图像预处理：对三维视网膜频域光学相干断层成像SD‑OCT图像进行下采样处理；步骤2.区域提取和划分：将图像分割为脉络膜新生血管CNV区域、外视网膜层、内视网膜层和脉络膜层4个区域；步骤3.网格化：对CNV区域、内视网膜层、外视网膜层和脉络膜层进行多尺度、自适应的四面体网格生成；步骤4.运用超弹性生物力学模型与反应扩散方程建模：在超弹性生物力学模型结构中，整合各区域不同的体积模量、剪切模量及其对应的体积弹性响应、等容弹性响应，将脉络膜新生血管生长后的质量改变作为源项加入到方程中，使变形梯度张量根据新生血管的生长持续变化；步骤5.优化模型，算出最佳准确率，进行参数检验；步骤6.根据每个时间点预测的参数拟合一条参数曲线，预测最后一个时间点的生长增值率参数，得到预测结果。

【技术特征摘要】
1.一种本构模型与有限元结合的脉络膜新生血管生长预测方法，其特征是，包括以下步骤：步骤1.图像预处理：对三维视网膜频域光学相干断层成像SD-OCT图像进行下采样处理；步骤2.区域提取和划分：将图像分割为脉络膜新生血管CNV区域、外视网膜层、内视网膜层和脉络膜层4个区域；步骤3.网格化：对CNV区域、内视网膜层、外视网膜层和脉络膜层进行多尺度、自适应的四面体网格生成；步骤4.运用超弹性生物力学模型与反应扩散方程建模：在超弹性生物力学模型结构中，整合各区域不同的体积模量、剪切模量及其对应的体积弹性响应、等容弹性响应，将脉络膜新生血管生长后的质量改变作为源项加入到方程中，使变形梯度张量根据新生血管的生长持续变化；步骤5.优化模型，算出最佳准确率，进行参数检验；步骤6.根据每个时间点预测的参数拟合一条参数曲线，预测最后一个时间点的生长增值率参数，得到预测结果。2.根据权利要求1所述的本构模型与有限元结合的脉络膜新生血管生长预测方法，其特征是，步骤1中的图像预处理包括以下步骤：(1a)下采样处理：(1b)使用CAVASS软件，进行OCT图像配准。3.根据权利要求1所述的本构模型与有限元结合的脉络膜新生血管生长预测方法，其特征是，步骤3中采用一个三维曲面和体积的网格生成器对CNV区域、内视网膜层、外视网膜层和脉络膜层进行多尺度、自适应的四面体网格生成。4.根据权利要求1所述的本构模型与有限元结合的脉络膜新生血管生长预测方法，其特征是，步骤4包括以下步骤：(4a)反应扩散方程的入侵和扩散逻辑回归型增长的反应扩散方程表示为：其中，t为时间，D是各向异性的扩散张量，包含三个组成部分Dx，Dy，Dz，这三个分量分别为x轴，y轴，z轴的分量；ρ代表生长增殖率，N是每单位体积的细胞数，Kmax代表最大细胞容量，当血管内皮生长因子的各向异性浓度时，有如下公式：(4b)质量效应与超弹性生物力学模型材料的非线性应力-应变行为的应变能密度函数为：

【专利技术属性】
技术研发人员：陈新建，左畅，朱伟芳，
申请(专利权)人：苏州大学，
类型：发明
国别省市：江苏,32