(19)国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202111665365.2 (22)申请日 2021.12.31 (65)同一申请的已公布的文献号 申请公布号 CN 114511670 A (43)申请公布日 2022.05.17 (73)专利权人 深圳市铱硙医疗科技有限公司 地址 518000 广东省深圳市南 山区粤海街 道高新区社区高新南六道6号迈科龙 大厦1001 (72)发明人 王思伦　肖焕辉　周竞宇　 (74)专利代理 机构 深圳市诺正鑫泽知识产权代 理有限公司 4 4689 专利代理师 颜燕红 (51)Int.Cl. G06T 17/00(2006.01)G06T 7/10(2017.01) G06K 9/62(2022.01) G06V 10/762(2022.01) (56)对比文件 CN 1094315 32 A,2019.0 3.08 CN 113850755 A,2021.12.28 CN 1094315 31 A,2019.0 3.08 CN 113034441 A,2021.0 6.25 CN 106485706 A,2017.0 3.08 US 2009016587 A1,20 09.01.15 US 201915974 4 A1,2019.0 5.30 审查员 王垚 (54)发明名称 基于动态灌注影像的血管重建方法、 装置、 设备和介质 (57)摘要 本申请涉及一种基于动态灌注影像的血管 重建方法、 装置、 设备和介质。 该方法包括： 从动 态灌注影像的各血管体素中， 筛选属于输入动脉 的第一血管体素和属于输出静脉的第二血管体 素； 分别基于各第一血管体素和各第二血管体素 的造影剂浓度变化曲线生成动脉输入曲线和静 脉输出曲线； 基于动脉输入曲线确定造影剂开始 到达输入动脉的起始时间， 以及基于动脉输入曲 线和静脉输出曲线的交点确定截止时间， 以确定 动脉时间区间； 确定各血管体素在动脉时间 区间 内的平均造影剂浓度， 以生成浓度分布三维图； 对浓度分布 三维图进行分割， 得到候选血管体素 分割图， 从中剔除达峰时间在动脉时间区间外的 血管体素后重建三维动脉血管图像。 本申请方案 能够避免局限性。 权利要求书2页 说明书14页 附图4页 CN 114511670 B 2022.08.30 CN 114511670 B 1.一种基于动态灌注影 像的血管重建方法， 其特 征在于， 所述方法包括： 从动态灌注影像的各血管体素中， 筛选属于输入动脉的第 一血管体素和属于输出静脉 的第二血 管体素； 基于各所述第一血管体素的造影剂浓度变化曲线生成动脉输入曲线， 以及， 基于各所 述第二血 管体素的造影剂浓度变化曲线生成静脉输出曲线； 基于所述动脉输入曲线确定造影剂开始到达所述输入动脉的起始时间， 以及基于所述 动脉输入曲线和所述静脉输出曲线的交点确定截止时间， 并根据所述起始时间和所述截止 时间确定动脉时间区间； 确定各所述血管体素在所述动脉时间区间内分别对应的平均造影剂浓度， 并基于各所 述血管体素的所述平均造影剂浓度生成浓度分布三维图； 对所述浓度分布三维图进行体素分割， 得到三维的候选血 管体素分割图； 从候选血管体素分割图中， 剔除造影剂浓度达峰时间在所述动脉时间区间外的血管体 素； 基于剔除血 管体素后的血 管体素分割图重建三维动脉 血管图像。 2.根据权利要求1所述的方法， 其特征在于， 所述从动态灌注影像的血管体素中， 筛选 属于输入动脉的第一血 管体素和 属于输出静脉的第二血 管体素包括： 获取动态灌注影 像的各血管体素分别对应的造影剂浓度变化曲线； 基于每个所述血管体素的造影剂浓度变化曲线， 生成每个所述血管体素的体素特征向 量； 基于所述体素 特征向量对各 所述血管体素进行聚类， 得到多个聚类簇； 从所述多个聚类簇 中， 筛选满足输入动脉特征条件的第 一聚类簇和满足输出静脉特征 条件的第二聚类簇； 从所述第一聚类簇 中获取所述第 一血管体素， 以及从所述第 二聚类簇中获取所述第 二 血管体素。 3.根据权利要求1所述的方法， 其特征在于， 所述动态灌注影像是在三维血管模型影像 中按照时间维度动态展示造影剂灌注情况的四维动态影 像； 所述确定各所述血管体素在所述动脉时间区间内分别对应的平均造影剂浓度， 并基于 各所述血管体素的所述平均造影剂浓度生成浓度分布三维图包括： 针对所述三维血管模型影像的任意一个所述血管体素， 基于所述血管体素对应的造影 剂浓度变化曲线， 确定所述血 管体素在动脉时间区间内各时间点对应的造影剂浓度； 将每个所述血管体素在动脉时间区间内各时间点对应的造影剂浓度在时间维度上求 平均值， 得到每 个所述血 管体素在所述动脉时间区间内对应的平均造影剂浓度； 在所述三维血管模型影像上标注各所述血管体素的所述平均造影剂浓度， 生成浓度分 布三维图。 4.根据权利要求1所述的方法， 其特征在于， 在所述从动态灌注影像的各血管体素中， 筛选属于输入动脉的第一血管体素和属于输出静脉的第二血管体素之前， 所述方法还包 括： 获取待处理动态灌注影像； 所述待处理动态灌注影像是连续多次地进行层面扫描得到 的动态影 像； 将对所述待处理动态灌注影像中同一位置上多次扫描得到的各影像帧分别作为目标权　利　要　求　书 1/2 页 2 CN 114511670 B 2帧， 针对每个所述 目标帧， 确定所述 目标帧与对应于同一位置的时间相邻帧之间的局部扫 描时间差； 所述时间相邻帧， 是在所述 目标帧的扫描时间的相邻扫描时间上对所述同一位 置扫描得到的影 像帧； 根据所述局部扫描时间差， 自适应地确定所述目标帧对应的时间滤波卷积核参数； 基于所述 时间滤波卷积核参数对所述目标帧进行平滑滤波处理， 得到所述动态灌注影 像。 5.根据权利要求 4所述的方法， 其特 征在于， 所述获取待处 理动态灌注影 像包括： 获取目标对象的原 始动脉灌注影 像； 对原始动脉灌注影像进行对象运动情况分析， 以识别所述原始动脉灌注影像中的偏位 帧； 所述偏位帧， 是表征目标对象发生移动而存在偏差的影 像帧； 对所述偏位帧进行 校正， 得到待处 理动态灌注影 像。 6.根据权利要求1所述的方法， 其特征在于， 在所述从候选血管体素分割图中， 剔除造 影剂浓度达峰时间在所述动脉时间区间外的血 管体素之前， 所述方法还 包括： 针对所述候选血管体素分割图中各候选血管体素， 基于各候选血管体素的造影剂浓度 变化曲线， 计算各候选血 管体素对应的造影剂浓度达峰时间； 各候选血管体素对应的造影剂浓度达峰时间分别与所述动脉时间区间进行比对， 以从 候选血管体素分割图中确定 造影剂浓度达峰时间在所述动脉时间区间外的血 管体素。 7.根据权利要求1至 6中任一项所述的方法， 其特 征在于， 所述方法还 包括： 从重建后的三维动脉 血管图像中提取 各动脉血管的中轴线； 基于预设的标准动脉血管配准到所述中轴 线上， 以将所述重建后的三维动脉血管图像 进行血管分段， 得到属于不同动脉类别的各 段子图像； 从不同动脉类别的各 段子图像中， 选取目标动脉类别的子图像； 基于选取的子图像中各 血管体素的造影剂浓度变化曲线， 生成目标动脉输入曲线。 8.一种基于动态灌注影 像的血管重建装置， 其特 征在于， 所述装置包括： 生成模块， 用于从动态灌注影像的各血管体素中， 筛选属于输入动脉的第一血管体素 和属于输出静脉的第二血管体素； 基于各所述第一血管体素的造影剂浓度变化曲线生成动 脉输入曲线， 以及， 基于各 所述第二血 管体素的造影剂浓度变化曲线生成静脉输出曲线； 动脉期确定模块， 用于基于所述动脉输入曲线确定造影剂开始到达所述输入动脉的起 始时间， 以及基于所述动脉输入曲线和所述静脉输出曲线的交点确定截止时间， 并根据所 述起始时间和所述截止时间确定动脉时间区间； 血管重建模块， 用于确定各所述血管体素在所述动脉时间区间内分别对应的平均造影 剂浓度， 并基于各所述血管体素 的所述平均造影剂浓度生成浓度分布三维 图； 对所述浓度 分布三维图进行体素分割， 得到三 维的候选血管体素分割图； 从候选血管体素分割图中， 剔 除造影剂浓度达峰时间在所述动脉时间区间外的血管体素； 基于剔除血管体素后的血管体 素分割图重建三维动脉 血管图像。 9.一种计算机设备， 包括存储器和处理器， 所述存储器存储有计算机程序， 其特征在 于， 所述处 理器执行所述计算机程序时实现权利要求1至7中任一项所述的方法的步骤。 10.一种计算机可读存储介质， 其上存储有计算机程序， 其特征在于， 所述计算机程序 被处理器执行时实现权利要求1至7中任一项所述的方法的步骤。权　利　要　求　书 2/2 页 3 CN 114511670 B 3