申请/专利权人：电子科技大学

申请日：2022-05-19

公开（公告）日：2024-05-17

公开（公告）号：CN114996996B

主分类号：G06F30/23

分类号：G06F30/23;G06T17/20;G06F111/10

优先权：

专利状态码：有效-授权

法律状态：2024.05.17#授权;2022.09.20#实质审查的生效;2022.09.02#公开

摘要：该发明公开了一种基于相场模型的3D毛细血管生长方法，属于数据处理领域，特别是模拟毛细血管生长的方法。本发明提出的3D血管生成方法能够以较为简洁的求解办法实现相场控制方程的求解，在对血管模型的实现流程上，本发明能够清晰精准的反映3D血管生长过程中的形态变化，能够为血管形成的原理和血管生长的实验提供更加准确的3D血管模型。

主权项：1.一种基于相场模型的3D毛细血管生长方法，其特征在于，相场序参量φ与血管结构之间的关系为：φ≥0.9代表血管，φ-0.9代表胞外基质，而中间的过渡区域代表血管壁；该方法包括：毛细血管连续生长方法和毛细血管尖端移动方法；1所述毛细血管连续生长方法包括：相场序参量的控制和血管内皮生长因子的控制；所述相场序参量的控制方程为： 其中，相场序参量φ即为内皮细胞密度，M表示内皮细胞扩散系数，表示梯度运算，ε为内皮细胞密度扩散系数，αp为内皮细胞增殖速率，Θ为海维赛德函数，c代表血管内皮生长因子浓度，t为时间；所述血管内皮生长因子的控制方程为： 其中，Dc为血管内皮生长因子扩散系数，αT为内皮细胞对血管内皮生长因子的消耗速率；采用如下方法相场序参量的控制方程进行计算；步骤A1：对相场序参量的控制方程作变量代换，以降低微分次数，得到微分方程组： 其中，步骤A2：对公式进行离散，得到： 其中，表示在时间步数为q时，空间位置l，m，n处，物理量μ的值，上下标q、l、m、n分别表示时间步、空间位置的x轴坐标、空间位置的y轴坐标、空间位置的z轴坐标，h表示对空间进行网格离散后每个网格的边长；步骤A3：对进行离散，得到： 对进行离散，得到： 对离散化后的公式进行离散计算，在计算过程中相场序参量值超过N的值格点，一律将此格点的值设为N；2所述毛细血管尖端移动方法为：步骤B1：设初始时刻t0尖端细胞集合ΩTECt0为： 式中，NTEC表示尖端细胞个数，i表示第i个尖端细胞，表示第i个尖端细胞在初始时刻t0所占据的球形区域；步骤B2：对满足如下条件的内皮细胞激活为尖端细胞；φ≥φact,c≥cact,G≥Gact, 其中，下标act表示对应符号设定的阈值，G表示内皮细胞所处的VEGF浓度梯度，表示第i个内皮细胞的圆心位置，x*表示区域中已存在的尖端细胞的圆心位置，Rc表示内皮细胞的半径大小；步骤B3：尖端细胞进行移动的方法为：对第i个尖端细胞，其在某一时间段的移动将由新位置 其中，Δt表示间隔时间，表示移动速度矢量；所述移动速度矢量的大小受血管内皮生长因子浓度梯度和趋化系数共同决定： 式中，K为扰动矩阵，χc为趋化常数，χφ为相场序参量负梯度对尖端细胞移动方向的贡献值，χp为尖端细胞原移动方向对尖端细胞移动方向的贡献值；所述移动速度矢量的方向由：前一时刻移动的速度方向、中心所处的相场序参量的负梯度、趋化因子提供的梯度确定；根据单个尖端细胞的移动方法，对所有尖端细胞进行移动；步骤B4：对尖端细胞进行分支的方法为：在每一个时间步内在血管占据的区域中寻找满足激活尖端细胞条件的格点，并且将以此点为原心，RTEC为半径的区域作为新尖端细胞所占据的区域，此区域便形成了一个能够移动的独立单元，在确定此区域成为尖端细胞单元后，此区域的φ值保持不变，但是作为新生尖端细胞，加入到tn时刻的尖端细胞域中；步骤B5：尖端细胞的融合的方法为：对于每个尖端细胞单元的中心坐标经过的位置，实时的更新到血管区域轨迹域中，当尖端细胞的伪足探测域探测到血管轨迹域时，此尖端细胞引导的血管便融入相遇的血管中，此尖端细胞将从尖端细胞域中移除。

全文数据：

权利要求：

百度查询： 电子科技大学 一种基于相场模型的3D毛细血管生长方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD