基于特定的生理参数获取血管压力值的方法及装置与流程

本发明涉及一种基于特定的生理参数获取血管压力值的方法及装置，属于医疗技术领域。

背景技术：

人体血液中的脂类及糖类物质在血管壁上的沉积将在血管壁上形成斑块，继而导致血管狭窄；特别是发生在心脏冠脉附近的血管狭窄将导致心肌供血不足，诱发冠心病、心绞痛等病症，对人类的健康造成严重威胁。据统计，我国现有冠心病患者约1100万人，心血管介入手术治疗患者数量每年增长大于10％。

冠脉造影、ct等常规医用检测手段虽然可以显示心脏冠脉血管狭窄的严重程度，但是并不能准确评价冠脉的缺血情况。为提高冠脉血管功能评价的准确性，1993年pijls提出了通过压力测定推算冠脉血管功能的新指标——血流储备分数(fractionalflowreserve，ffr)，经过长期的基础与临床研究，ffr已成为冠脉狭窄功能性评价的金标准。

血流储备分数(ffr)通常是指心肌血流储备分数，定义为病变冠脉能为心肌提供的最大血流与该冠脉完全正常时最大供血流量之比，研究表明，在冠脉最大充血状态下，血流量的比值可以用压力值来代替。即ffr值的测量可在冠脉最大充血状态下，通过压力传感器对冠脉远端狭窄处的压力和冠脉狭窄近端压力进行测定继而计算得出。近年来，基于压力导丝测量ffr值的方法逐渐进入临床应用，成为冠心病患者获得精准诊断的有效方法；然而，由于压力导丝在介入过程中易对病人的血管造成损伤；同时，通过压力导丝对ffr值进行测定需要注射腺苷/atp等药物保证冠脉达到最大充血状态，部分病人会因药物的注射感到不适，使得基于压力导丝测量ffr值的方法存在较大的局限性。

随着ct与三维造影重建技术的发展及3d冠状动脉几何重建技术在血液力学研究领域的推广应用，同时，为减少ffr值测量过程中对人体带来的伤害及测量成本，基于医疗影像学的ffr计算技术已成为研究重点。

事实上，现有技术虽然从不同角度、不同方法中给出了确定血流储备分数(ffr)的方法，但其实质均是通过感兴趣区域血管段近端终点处的血流压力pa和感兴趣区域血管段近端终点处和远端终点处的血流压力的差值δp来计算ffr。而在血液流动的实际过程中，病变的位置、大小和类型等因素均会对血流压力的差值δp的计算产生影响；特别地，个体的生理参数不同将导致感兴趣区域血管段近端终点处的血流压力pa的获取和感兴趣区域血管段近端终点处和远端终点处的血流压力的差值δp计算产生误差，同时生理参数的不同在图像的获取过程中，也将导致图像获取的质量，从而导致血流压力的差值δp计算过程中的参数获取存在误差，进一步导致现有技术中，通过血流压力pa、血流压力的差值δp计算获得的ffr、血流特征值等偏离实际值，致使对血管系统的功能评价的结果存在误差。

有鉴于此，确有必要提供一种新的获取血管压力值的方法，以解决上述问题。

技术实现要素：

本发明的目的在于提供一种基于特定的生理参数获取血管压力值的方法及装置，以至少解决现有技术中存在的技术问题之一。本发明提供的基于特定的生理参数获取血管压力值的装置，通过引入个体的生理参数，对感兴趣区域血管段两端的压力值计算过程中的各类数据进行修正，明确不同的生理参数等对血管压力值计算的影响，提高血管压力值计算的准确性。

为实现上述发明目的，本发明提供了一种基于特定的生理参数获取血管压力值的方法，所述基于特定的生理参数获取血管压力值的方法包括：

s1、采集至少一部分的血管系统的解剖数据，根据所述解剖数据获取感兴趣区域的几何参数，并建立感兴趣区域的几何模型；

s2、根据感兴趣区域的解剖数据和/或个体特异性数据建立感兴趣区域的血流模型；

s3、基于一个或多个特定的生理参数，对所述几何模型和/或所述血流模型进行修正，获取感兴趣区域血管压力值的计算模型及感兴趣区域的血流速度v；

s4、根据所述血管压力值计算模型、血流速度v并结合血流动力学，获取感兴趣区域近端终点处的血流压力pa，及相应状态下感兴趣区域近端终点处与远端终点处的压力差数值δp。

作为本发明的进一步改进，所述血管压力值的计算模型包括所述几何模型中至少一部分的血管段和/或血管树，以及感兴趣区域近端终点和远端终点之间各个位置处的横截面形态模型。

作为本发明的进一步改进，所述横截面形态模型包括各横截面上斑块的有无、斑块的位置、斑块的大小、斑块形成的角度、斑块的组成及斑块组成的变化、斑块的形状及斑块形状的变化。

作为本发明的进一步改进，所述横截面形态模型的建立包括：

s1、定义感兴趣区域血管段近端终点处的横截面为参考面，通过中心线提取与建立方法获得所述几何模型的中心径线；

s2、以所述参考面的中心点为原点建立坐标系，沿垂直所述中心径线的方向对所述感兴趣区域血管段进行分割，将各横截面内外边缘投影在所述坐标系中，以获取感兴趣区域血管段在各个位置处管腔横截面的平面几何图像，横截面形态模型建立结束。

作为本发明的进一步改进，所述血管压力值的计算模型还包括以感兴趣区域的近端终点为参考点，对不同尺度下的所述横截面形态模型进行拟合，获取的血管管腔的形态差异函数f(x)，所述尺度为计算形态差异函数f(x)时相邻两横截面之间的距离。

作为本发明的进一步改进，所述形态差异函数f(x)的获取包括：

基于横截面形态模型，建立各横截面的形态函数；

对相邻两横截面的形态函数进行拟合，并获取相邻两横截面在不同尺度下的差异变化函数；

以感兴趣区域血管段的近端终点为参考点，根据差异变化函数获取管腔形态随着到参考点的距离x的变化率，对感兴趣区域血管段从近端终点到远端终点范围内的位置参数进行归一化处理，获取形态差异函数f(x)；优选的，所述形态函数包括面积函数、直径函数和边缘位置函数。

作为本发明的进一步改进，所述血流模型包括固定血流模型及个性化血流模型；优选的，所述个性化血流模型包括静息态血流模型和负荷态血流模型。

作为本发明的进一步改进，所述血流模型包括感兴趣区域的血流速度v，且当所述血流模型为静息态血流模型时，所述血流速度v可通过血管中流体充盈的速度计算获得，或者通过血管树的形态计算获得；其中，所述血管树的形态至少包括所述血管树的面积、体积和血管树中血管段的管腔直径中的一种或几种；所述血流速度v通过所述血管树的形态计算获得时，所述几何参数还包括所述血管树中血管段的长度、灌注面积及分支角度中的一种或几种。

作为本发明的进一步改进，所述特定的生理参数包括年龄、性别、血压及身体质量指数等可直接获取的生理信息。

作为本发明的进一步改进，当所述血管系统为冠脉血管系统时，所述特定的生理参数还包括冠状动脉优势类型。

为实现上述发明目的，本发明还提供了一种基于特定的生理参数获取血管压力值的装置，其特征在于，所述基于特定的生理参数获取血管压力值的装置包括：

数据采集器，所述数据采集器用于获取及存储血管装置的解剖模型中感兴趣区域的几何参数；

压力值处理器，所述压力值处理器用于建立感兴趣区域的血流模型和基于所述几何参数，建立对应感兴趣区域的几何模型；

所述压力值处理器还用于，基于所述几何模型和所述血流模型获取血管压力值计算模型；同时，根据所述血管压力值计算模型并结合血流动力学，获取感兴趣区域近端终点处的血流压力pa，及感兴趣区域血管段近端终点与远端终点之间的压力差数值δp。

作为本发明的进一步改进，所述基于特定的生理参数获取血管压力值的装置还包括纠偏处理器，所述纠偏处理器用于接收个体的一个或多个生理参数，处理后生成纠偏参数，并将所述纠偏参数传递至压力值处理器中对所述几何模型和/或所述血流模型进行修正。

作为本发明的进一步改进，所述几何模型为所述压力值处理器通过对所述数据采集器传递的所述解剖模型的几何参数进行测算，并结合所述纠偏处理器传递的纠偏参数经拟合校准获得。

作为本发明的进一步改进，所述血管压力值计算模型包括感兴趣区域近端终点和远端终点之间各个位置处的横截面形态模型及感兴趣区域几何模型中至少一部分的血管段和/或血管树，所述血管树包括至少一段主动脉或者包括至少一段主动脉和由所述主动脉发出的多个冠状动脉。

作为本发明的进一步改进，所述横截面形态模型为所述压力值处理器通过所述几何模型直接/间接获得；所述横截面形态模型包括各横截面上斑块的有无、斑块的位置、斑块的大小、斑块形成的角度、斑块的组成及斑块组成的变化、斑块的形状及斑块形状的变化。

作为本发明的进一步改进，所述基于特定的生理参数获取血管压力值的装置还包括速度采集模块，所述速度采集模块用于获取感兴趣区域的血流速度v，所述血流速度v用以推算所述感兴趣区域近端终点处的第一血流压力p1及感兴趣区域近端终点处的第一血流压力p1与远端终点处的第二血流压力p2之间的压力差数值δp。

作为本发明的进一步改进，所述速度采集器包括速度计算模块及速度提取模块；所述速度提取模块通过数据采集模块直接获取血流速度信息，或通过血流模型直接提取血流速度v。

作为本发明的进一步改进，所述速度计算模块包括速度转换模块及速度测算模块，所述血流速度v通过血管中流体充盈的速度经速度转换模块转换获得，或通过几何模型中血管树的形态经速度测算模块计算获得。

为实现上述发明目的，本发明还提供了一种用于获取患者血管压力差的设备，所述设备具有处理器，其中，所述处理器被设置为使得所述设备执行以下步骤：

收集患者特定生理参数和待检血管的几何参数；

根据所述待检血管的几何参数建立患者的血管模型；

根据患者特定生理参数修正所述血管模型；

提供至少一种血管压力差数值的计算模型；

基于修正后的所述血管模型与所述血管压力差数值的计算模型确定待检血管的血管压力差。

作为本发明的进一步改进，所述患者的特定生理参数包括年龄、性别、体温、身体质量指数等可直接获取的人体生理信息中的一个或多个。

作为本发明的进一步改进，所述血管压力差数值的计算模型是基于多尺度计算方法而建立的。

本发明的有益效果是：本发明的基于特定的生理参数获取血管压力值的方法，通过在血管压力值的计算过程中引入个体特定的生理参数，一方面可保证血管压力值计算初期各类参数的提取和建立准确，另一方面，可对血管压力值计算过程中的血流模型进行修正，保证几何模型和血流模型建立的准确性，进一步保证通过所述几何模型和所述血流模型获取的相关参数准确，使得通过本发明的基于特定的生理参数获取血管压力值的方法计算得到的血管压力值可准确反映感兴趣区域的血流压力，且结果准确可靠。

附图说明

图1是本发明感兴趣区域血管段在一种形态下的几何模型的示意图。

图2是图1中d1位置处横截面形态模型的结构示意图。

图3是图1中d2位置处横截面形态模型的结构示意图。

图4是图2和图3中d1和d2位置处横截面形态模型拟合后的结构示意图。

图5是本发明感兴趣区域血管段在另一种形态下的几何模型的示意图。

图6是图5中d1位置处横截面形态模型的结构示意图。

图7是图5中d2位置处横截面形态模型的结构示意图。

图8是图6和图7中d1和d2位置处横截面形态模型拟合后的结构示意图。

图9是本发明基于特定的生理参数获取血管压力值的装置的结构框图。

具体实施方式

为了使本发明的目的、技术方案和优点更加清楚，下面结合附图和具体实施例对本发明进行详细描述。

本发明提供了一种基于特定的生理参数获取血管压力值的方法，所述基于特定的生理参数获取血管压力值的方法包括一下步骤：

s1、采集至少一部分的血管系统的解剖数据，根据所述解剖数据获取感兴趣区域的几何参数，并建立感兴趣区域的几何模型；

s2、根据感兴趣区域的解剖数据和/或个体特异性数据建立感兴趣区域的血流模型；

s3、基于一个或多个特定的生理参数，对所述几何模型和/或所述血流模型进行修正，获取感兴趣区域血管压力值的计算模型及感兴趣区域的血流速度v；

s4、根据所述血管压力值计算模型、血流速度v并结合血流动力学，获取感兴趣区域近端终点处的血流压力pa，及相应状态下感兴趣区域近端终点处与远端终点处的压力差数值δp。

在本发明中，所述血管压力值包括感兴趣区域血管段近端终点的血流压力pa及相应状态下感兴趣区域近端终点处与远端终点处的压力差数值δp；当然在其他实施例中，所述血管压力值还包括感兴趣区域远端终点血流压力等其它用于表征感兴趣区域血管段压力情况的数值；进一步的，在本发明中所述特定的生理参数包括年龄、性别、血压及身体质量指数(bmi)等可直接获取的个体生理信息。

进一步的，所述几何模型为反映个体至少一部分血管系统的几何形状的三维模型，且所述几何模型为对个体血管系统的解剖数据进行建模后获得；进一步的，所述血流模型为反映个体感兴趣区域血管段中血液流动情况的数据模型，包括感兴趣区域血管段中血液的流动情况、血液的性质等可反映血液特性的数据，且在本发明中所述血流模型为通过感兴趣区域血管段的解剖数据和/或个体的特异性数据建立

需要说明的是，在本发明中所述个体血管系统的解剖数据，可通过ct设备、oct设备、ivus设备和造影设备等常用图像生成设备生成获取，具体来讲，当采用基于x射线的装置获取个体血管系统的解剖数据时，个体通过口服/注射放射剂后，x射线可对含有放射剂的体液进行显影，进一步反映出个体血管系统/脏器的图像。

然而，在使用x射线装置获取个体的解剖数据时，相同的放射剂量作用在肥胖个体上时，x射线装置发出的x射线无法全部穿过个体准确的形成个体的影响，致使肥胖个体经x射线装置获取的解剖模型边缘位置影像噪声大、质量差，进一步致使获取的解剖数据存在误差；因此，在使用x射线装置获取解剖数据时，需根据个体的肥胖程度对解剖数据进行修正，以保证几何模型和/或血流模型建立准确。

具体来讲，个体的肥胖程度可通过身体质量指数(bmi)体现，所述身体质量指数(bmi)为个体体重(kg)与身高(m)的平方的比值，当个体的身体质量指数不同时，调整x射线装置中x射线管的工作参数，可提高x射线的强度，以保证x射线作用在个体上时形成的解剖模型清晰，继而提高解剖数据的准确性。

进一步的，当所述个体的身体质量指数＜18.5时，所述x射线管的管电压为80～120kv；当所述个体的身体质量指数在18.5～23.9之间时，所述x射线管的管电压为120kv；当所述个体的身体质量指数＞24时，所述x射线管的管电压为120～140kv；如此设置，可在保证个体放射剂摄入量不变时，获得的几何参数准确，进一步提高几何模型建立的准确性。

为进一步提高x射线装置获取解剖数据的准确性，当个体的身体质量指数在24～27之间时，所述x射线管的管电压为120～130kv；当个体的身体质量指数在28～32之间时，所述x射线管的管电压为130～135kv；当个体的身体质量指数＞32时，所述x射线管的管电压为135～140kv。

所述血流模型包括固定血流模型及个性化血流模型；且所述血流模型既可为数据计算模型也可为三维流体流动模型；进一步的，所述固定血流模型即为经验值血流模型，为根据临床实际经验，通过大数据采集及模拟的方法直接建立；所述个性化血流模型包括静息态血流模型和负荷态血流模型。

在本发明中，所述血流速度v和所述血流压力pa均通过所述血流模型直接/间接获取。具体来讲，当所述血流模型为固定血流模型或静息态血流模型时，所述血流压力pa可通过获取的个体的收缩压混合舒张压计算获得，且此时所述血流压力pa＝1/3收缩压+2/3舒张压；当所述血流模型为负荷态血流模型，此时所述血流压力pa可通过所述负荷态血流模型直接测量获得。

所述个性化血流模型包括静息态血流模型和负荷态血流模型；当所述血流模型为静息态血流模型时，所述血流速度v可通过流体充盈的速度计算获得；在本发明的一个实施例中，所述静息态血流模型为造影剂血流模型，此时所述血流速度v为利用灰度时间拟合函数获得的目标血管在造影过程中造影剂的平均流动速度；或者利用timi数帧法计算获得的所述目标血管在造影过程中造影剂的平均流动速度。

当所述静息态血流模型为ct血流模型时，所述血流速度v可通过血管树的形态计算获得，所述血管树的形态至少包括所述血管树的面积、体积和血管树中血管段的管腔直径中的一种或几种；且当所述血流速度通过所述血管树的形态计算获得时，所述几何参数还包括所述血管树中血管段的长度、灌注面积及分支角度中的一种或几种。

在本发明的另一实施例中，所述血流模型为负荷态血流模型，此时所述血流速度v为注射腺苷血管充分扩张后的血流速度v，且此时，所述血流速度v为最大血流速度vmax。

特别地，在本发明中所述血流速度v包括目标血管处于最大充血状态下的血流速度vmax和静息状态下的血流速度vqc，当目标血管位于冠脉区域时，所述血流速度v为最大充血状态下的血流速度vmax，进一步的血流速度vmax可直接通过血流模型获取，或通过血流模型计算的血流速度v转换获得；当目标血管位于外周血管系统时，所述血流速度v为静息态下的血流速度vqc。

所述血管压力值的计算模型为经所述几何模型和所述血流模型拟合处理获得，且所述几何模型和/或所述血流模型均为经个体的生理参数修正后几何模型和/或血流模型。

进一步的，所述血管压力值计算模型包括感兴趣区域近端终点和远端终点之间各个位置处的横截面形态模型及感兴趣区域几何模型中至少一部分的血管段和/或血管树，所述血管树包括至少一段主动脉或者包括至少一段主动脉和由所述主动脉发出的多个冠状动脉，且每个所述血管段和/或血管树均包括所述感兴趣区域的形状、直径和面积等解剖数据，进一步的，所述解剖数据还包括血管段的弯曲角度等可以反映感兴趣区域血管段实际形态的参数。

所述横截面形态模型用于表示感兴趣区域血管段横截面的形状、面积和直径的几何形态，且在本发明中所述横截面形态模型为通过所述几何模型直接/间接获得，具体来讲，所述横截面形态模型的建立包括：

s1、定义感兴趣区域血管段近端终点处的横截面为参考面，通过中心线提取与建立方法获得所述几何模型的中心径线；

s2、以所述参考面的中心点为原点建立坐标系，沿垂直所述中心径线的方向对所述感兴趣区域血管段进行分割，将各横截面内外边缘投影在所述坐标系中，以获取目标血管在各个位置处管腔横截面的平面几何图像，横截面形态模型建立结束。

进一步的，所述横截面形态模型还包括各横截面位置处的斑块信息，所述斑块信息即为感兴趣区域的病变信息，大量数据表明：当斑块(即为病变)的长度＞20mm时，将导致感兴趣区域血管段压力值数值δp的升高，偏离实际值；而当同一横截面处斑块的组成复杂或尺寸过大致使感兴趣区域血管段的狭窄率高，则会进一步导致感兴趣区域血管段压力差数值δp的计算；同时当所述斑块处于不同的位置处时，目标血管所供应的心肌体积区域不同，将导致病变位置与非病变位置处的比例发生变化，进一步影响血流速度v，从而导致感兴趣区域血管段压力差数值δp的计算发生偏差。

因此，在建立所述横截面形态模型时，所述斑块信息还需包括斑块的有无、斑块的位置、斑块的大小、斑块形成的角度、斑块的组成及斑块组成的变化、斑块的形状及斑块形状的变化，且在本发明中，各个位置处的管腔横截面的平面几何图像均需以步骤s2中建立的坐标系为参考，明确各横截面上斑块的位置，以方便横截面形态模型的后续拟合。

需要说明的是，在所述横截面形态模型的建立过程中，当所述解剖数据为采用ct、oct、ivus等检测手段获取时，所述横截面形态模型可通过所述几何模型直接获取，只需保证每个所述横截面形态模型的原点及坐标方向一致即可；当所述解剖数据为采用x射线等检测手段获取时，由于所述几何模型为沿血流方向延伸的立体模型，则在通过所述几何模型建立所述横截面形态模型时，需对所述几何模型进行坐标转换，以准确反应各个横截面的截面形态。

为进一步保证所述基于特定的生理参数获取血管压力值的方法获取的压力差数值δp准确，所述基于特定的生理参数获取血管压力值的方法还包括对不同尺度下的所述横截面形态模型进行拟合，计算目标血管管腔的形态差异函数f(x)。其中，所述形态差异函数f(x)用于表示目标血管不同位置处的横截面形态变化随着该位置到参考点的距离x变化的函数；且所述形态差异函数f(x)的获取包括：

基于横截面形态模型，建立各横截面的形态函数；

对相邻两横截面的形态函数进行拟合，并获取相邻两横截面在不同尺度下的差异变化函数；

以目标血管的近端终点为参考点，根据差异变化函数获取管腔形态随着到参考点的距离x的变化率，对目标血管从近端终点到远端终点范围内的位置参数进行归一化处理，以最终获取形态差异函数f(x)。

所述形态函数包括面积函数、直径函数或边缘距离函数，即在本发明中可通过各横截面面积、直径或边缘距离函数之间的拟合，获取相邻两横截面在不同尺度下的差异变化函数；进一步的，通过差异变化函数获取管腔形态随着到参考点的距离x的变化率，获得形态差异函数f(x)。

具体来讲，当所述形态函数为面积函数时，如图1至图4，对d1和d2位置处的两横截面形态模型进行拟合，d1、d2位置处的横截面形态模型拟合后，有血管管腔斑块增加的区域为a1，对应的面积s1；血管管腔减少的区域为a2，对应的面积s2。由于所述d1和d2位置处的血管管腔(斑块)不重叠，因此当血流经d1处流向d2处时，血流压力将随之发生变化；此时，差异变化函数即为血管管腔中非重叠区域(s1、s2)与重叠区域之间面积(s3)的比值，或者为非重叠区域的面积(s1、s2)与总面积(s1、s2、s3)的比值；且此时，所述形态差异函数f(x)＞0，即横截面d1和d2之间存在压力差。

进一步的，当所述d1和d2位置处的血管管腔(斑块)完全重叠时，如图5至图8，所述区域a1与a2完全重叠，即非重叠区域a1与a2的面积s1＝s2＝0，此时，差异变化函数为0，即所述形态差异函数f(x)＝0，此时，横截面d1和d2之间不存在压力差。

进一步的，所述压力差数值δp为通过不同尺度下的形态差异函数f(x)和感兴趣区域血管段的血流速度v计算获得，且所述压力差数值δp在不同尺度下的计算公式为：

δp＝(c1v+c2v²+…+cmv^m)

*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]

其中，c1、c2、…、cm分别代表血流速度v的参数系数，所述参数系数包括血液粘度影响因素、血液湍流影响因素及粘滞系数等多个参数系数；进一步的，m为大于等于1的自然数，以分别代表不同参数系数对血流速度v的影响，以对压力差数值δp进行修正，保证压力差数值δp计算的准确性。优选的，在本发明中m的取值为2，且当m为2时，c1为因血液流动摩擦产生的参数系数，c2为血液湍流产生的参数系数。

α1、α2、…、αn分别为不同尺度下血管管腔的形态差异函数f1(x)、f2(x)、…、fn(x)的加权系数，其中，n为尺度为大于等于1的自然数；进一步的，所述加权系数的增加可进一步对形态差异函数f(x)进行修正，保证两横截面之间形态差异拟合计算的准确性。

具体来讲，所述不同尺度包括第一尺度、第二尺度、……、第n尺度；

所述第一尺度形态差异函数f1(x)用于检测第一种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异；

所述第二尺度形态差异函数f2(x)用于检测第二种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异；

……

所述第n尺度形态差异函数fn(x)用于检测第n种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异。

进一步的，在本发明的另一实施例中，所述血管压力差的计算还可与所述感兴趣区域的血流速度v无关，此时所述δp在不同尺度下的计算公式为：

δp＝k*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]

其中，k为修正参数，且k为常数为基于个体信息直接/间接获取的数值；

α1、α2、…、αn分别为不同尺度下血管管腔的形态差异函数f1(x)、f2(x)、…、fn(x)的加权系数；

优选的，所述不同尺度包括第一尺度、第二尺度、……、第n尺度；

所述第一尺度形态差异函数f1(x)用于检测第一种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异；

所述第二尺度形态差异函数f2(x)用于检测第二种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异；

……

所述第n尺度形态差异函数fn(x)用于检测第n种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异；其中，所述n为大于等于1的自然数。

以下说明书部分将以具体的影响血管压力值的生理参数为例，对本发明的基于特定的生理参数获取血管压力值的方法进行详细阐述。

实施例1：

在本实施例中，所述生理参数为血压，当个体存在血压异常时，异常的血压变化将导致血管中近端终点处的血流压力pa发生变化。

因此，当个体存在血压异常时，应对感兴趣区域的血流压力pa进行修正，此时，所述血流压力pa包括基于所述血流模型直接获取的第一血流压力p1及经修正后获得第二血流压力p2，且所述第一血流压力p1和所述第二血流压力p2之间满足关系式：

p2＝ω×p1；

其中，ω为与个体的血压信息有关的纠偏参数。

实施例2：

在本实施例中，所述生理参数为冠脉优势类型，且所述几何模型为个体冠脉血管系统位置处的模型。具体来讲，冠脉优势类型一般包括左优型、右优型以及均衡型；对左优型而言，由于左冠状动脉灌注的心肌区域大，因此同一心动周期左冠脉的血流量大于右冠脉的血流量；而对于右优型而言，右冠状动脉灌注的心肌区域大，则同一心动周期右冠脉的最大血流量大于左冠脉；而当斑块位于不同的优势型的不同血管内时，由于灌注的体积发生变化因此感兴趣区域的血流速度v也随之发生变化，进一步的，血流速度v的变化将影响血管压力值中压力差数值δp的计算。

因此，当个体感兴趣区域的冠脉优势类型发生变化时，应对感兴趣区域的血流速度v进行修正，具体来讲，在本实施例中，所述血流速度v包括基于所述血流模型直接获取的第一血流速度v0及经修正后获得第二血流速度v1，且所述第一血流速度v0和所述第二血流速度v1之间满足关系式：

v1＝ω*v0；

其中，ω为与个体的冠脉优势类型有关的纠偏参数。

且此时所述压力差数值δp在不同尺度下的计算公式为：

δp＝(c1v1+c2v1²+…+cmv1^m)

*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]

实施例3：

在本实施例中，所述生理参数为性别。具体来讲，男性比女性心肌收缩能力更强，其所能达到的血流速度v也较大，而由于在对目标血管两端的压力差数值δp进行计算时，血流速度v多是通过血流模型直接/间接的获取，而血流模型的建立常常会忽略性别对计算结果造成的影响，因此，在所述压力差数值δp的计算过程中存在误差。

因此，当个体的性别不同时，应对感兴趣区域的血流速度v进行修正，具体来讲，在本实施例中，所述血流速度v包括基于所述血流模型直接获取的第一血流速度v0及经修正后获得第二血流速度v1，且所述第一血流速度v0和所述第二血流速度v1之间满足关系式：

v1＝ω*v0；

其中，ω为与个体的冠脉优势类型有关的纠偏参数。

且此时所述压力差数值δp在不同尺度下的计算公式为：

δp＝(c1v1+c2v1²+…+cmv1^m)

*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]

实施例4：

在本实施例中，所述生理参数为年龄。具体来讲，成年个体随着年龄的增大，心肌微循环功能逐渐降低，单一心动周期流经感兴趣区域的血流量减少，从而感兴趣区域的血流速度v也将随之改变，而常规计算中，血流速度v多是通过血流模型直接/间接的获取，并未考虑年龄因素对血流速度的影响，因此，在所述压力差数值δp的计算过程中存在误差。

因此，当个体的年龄不同时，应对感兴趣区域的血流速度v进行修正，具体来讲，在本实施例中，所述血流速度v包括基于所述血流模型直接获取的第一血流速度v0及经修正后获得第二血流速度v1，且所述第一血流速度v0和所述第二血流速度v1之间满足关系式：

v1＝ω*v0；

其中，ω为与个体的年龄有关的纠偏参数。

且此时所述压力差数值δp在不同尺度下的计算公式为：

δp＝(c1v1+c2v1²+…+cmv1^m)

*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]

需要说明的是，在本发明的实施例中，所述纠偏参数ω是根据临床实际经验，通过大数据采集及模拟的方法获取的经验值。

请参参阅图9所示，本发明还提供了一种基于特定的生理参数获取血管压力值的装置，所述基于特定的生理参数获取血管压力值的装置包括：

所述基于特定的生理参数获取血管压力值的装置包括：

数据采集器，所述数据采集器用于获取及存储血管装置的解剖模型中感兴趣区域的几何参数；

压力值处理器，所述压力值处理器用于建立感兴趣区域的血流模型和基于所述几何参数，建立对应感兴趣区域的几何模型；

所述压力值处理器还用于，基于所述几何模型和所述血流模型获取血管压力值计算模型；同时，根据所述血管压力值计算模型并结合血流动力学，获取感兴趣区域近端终点处的血流压力pa，及感兴趣区域血管段近端终点与远端终点之间的压力差数值δp。

进一步的，所述几何模型为通过对所述解剖模型中感兴趣区域的几何参数进行测算，并拟合校准获得；所述血管压力值的计算模型包括感兴趣区域近端终点和远端终点之间各个位置处的横截面形态模型及感兴趣区域几何模型中至少一部分的血管段和/或血管树，所述血管树包括至少一段主动脉或者包括至少一段主动脉和由所述主动脉发出的多个冠状动脉，且每个所述血管段和/或血管树均包括所述感兴趣区域的形状、直径和面积等解剖数据，进一步的，所述解剖数据还包括血管段的弯曲角度等可以反映感兴趣区域血管段实际形态的参数。

所述横截面形态模型为通过所述几何模型直接/间接获得，所述横截面形态模型包括各横截面上斑块的有无、斑块的位置、斑块的大小斑块形成的角度、斑块的组成及斑块组成的变化、斑块的形状及斑块形状的变化。

进一步的，所述获取血管压力值的装置还包括速度采集器，所述速度采集器用于获取感兴趣区域的血流速度，所述血流速度用以推算所述感兴趣区域近端终点处的第一血流压力pa及感兴趣区域近端终点与远端终点之间的压力差数值δp。

所述速度采集器包括速度计算模块及速度提取模块；所述速度提取模块可通过所述数据采集器直接采集获得血流速度，也可通过所述血流模型直接提取血流速度。

所述速度计算模块包括速度转换模块及速度测算模块，所述血流速度可通过血管中流体充盈的速度经所述速度转换模块转换获得，还可通过几何模型中血管树的形态经速度测算模块计算获得。

优选的，在本发明中所述压力差数值δp通过如下公式计算获得：

δp＝(c1v+c2v²+…+cmv^m)

*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]

其中，c1、c2、…、cm分别代表血流速度v的参数系数，所述参数系数包括血液粘度影响因素、血液湍流影响因素及粘滞系数等多个参数系数；进一步的，m为大于等于1的自然数，以分别代表不同参数系数对血流速度v的影响，以对压力差数值δp进行修正，保证压力差数值δp计算的准确性。优选的，在本发明中m的取值为2，且当m为2时，c1为因血液流动摩擦产生的参数系数，c2为血液湍流产生的参数系数。

所述α1、α2、…、αn分别为不同尺度下血管管腔的形态差异函数f1(x)、f2(x)、…、fn(x)的加权系数，其中，n为尺度为大于等于1的自然数；进一步的，所述加权系数的增加可进一步对形态差异函数f(x)进行修正，保证两横截面之间形态差异拟合计算的准确性。

进一步的，所述基于特定的生理参数获取血管压力值的装置还包括纠偏处理器，所述纠偏处理器用于接收个体的一个或多个生理参数，处理后生成纠偏参数，并将所述纠偏参数传递至压力值处理器中对所述几何模型和/或所述血流模型进行修正。

需要说明的是，所述纠偏参数既可为对所述血流速度v进行修正，继而保证所述血管压力值计算的准确性；还可直接对所述血管压力值(如感兴趣区域近端终点的血流压力pa、压力差数值δp)进行修正。

进一步的，本发明还提供一种用于获取个体血管压力差的设备，所述设备具有处理器，其中，所述处理器被设置为使得所述设备执行以下步骤：

收集个体特定生理参数和待检血管的几何参数；

根据所述待检血管的几何参数建立个体的血管模型；

根据个体特定生理参数修正所述血管模型；

提供至少一种血管压力差的计算模型；

基于修正后的所述血管模型与所述血管压力差的计算模型确定待检血管的血管压力差。

所述“处理器”包括接收和/或生成信号的任意装置，所述处理器处理的数据可以是文本消息、物体/流体运动的指令、应用程序的输入或一些其它信息；所述待检血管的备选术语可以为目标血管或感兴趣血管；且所述待检血管包括冠脉血管、由冠脉血管发出的分支血管、血管树和单支血管段等个体任意位置处的血管组织；所述血管模型至少包括所述第二几何模型和所述第二血流模型中的一种，且所述血管模型的备选术语还可为管腔模型、流体流动模型等可反映个体待检血管形态和血管内流体流动情况的模型，进一步的，所述血管模型包括待检血管的长度、直径、弯曲角度及待检血管中分支血管的存在、分支血管的角度、分支血管的数量等与所述待检血管的几何形貌有关的数据。

在本实施例中，所述管腔形态模型的备选术语还可为横截面形态模型，且所述管腔形态模型包括斑块的有无、斑块的位置、斑块的大小、斑块形成的角度、斑块的组成及斑块组成的变化、斑块的形状及斑块形状的变化；进一步的所述管腔形态模型的建立包括以下步骤：

s1、定义待检近端终点处的横截面为参考面，通过中心线提取方法，建立获取所述血管模型的中心径线；

s2、以所述参考面的中心点为原点建立坐标系，沿垂直所述中心径线的方向对所述待检血管进行分割，将各横截面内外边缘投影在所述坐标系中，以获取待检血管在各个位置处管腔形态的平面几何图像，管腔形态模型建立结束。

在本发明中，各个位置处的管腔形态的平面几何图像均需以步骤s2中建立的坐标系为参考，明确各管腔截面上斑块的位置，以方便管腔形态模型的后续拟合。

需要说明的是，在所述管腔形态模型的建立过程中，当所述解剖数据为采用ct、oct、ivus等检测手段获取时，所述管腔形态模型可通过所述血管模型直接获取，只需保证每个所述管腔形态模型的原点及坐标方向一致即可；当所述解剖数据为采用x射线等检测手段获取时，由于所述血管模型为沿血流方向延伸的立体模型，则在通过所述血管模型建立所述管腔形态模型时，需对所述血管模型进行坐标转换，以准确反映各个截面的截面形态。

所述处理器还用于基于预设的形态差异函数，通过所述管腔形态模型以及所述血管模型确定待检血管任意两位置间的血管压力差。其中，所述形态差异函数通过所述管腔形态模型拟合建立获取，用于表示待检血管不同位置处的管腔形态变化随着该位置到参考点的距离x变化的函数；且所述形态差异函数包括与待检血管的面积、体积、边缘位置和边缘形态有关的可以体现待检血管任意两位置间形态差异的差异函数，且所述差异函数可通过管腔形态模型直接/间接获取。

所述解剖数据在其他实施例中还可定义为解剖数据等可从图像获取装置直接和/或间接获取的可反映管腔形态的参数。

即在另一上下文中，所述处理器、待检血管、解剖数据、管腔形态模型和血管模型可以为具有相同含义的不同名称。

所述尺度为所述尺度为相邻两横截面之间的距离；所述不同尺度包括第一尺度、第二尺度、……、第n尺度；

所述第一尺度下的形态差异函数f1(x)用于检测第一种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异；

所述第二尺度下的形态差异函数f2(x)用于检测第二种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异；

……

所述第n尺度下的形态差异函数fn(x)用于检测第n种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异。

进一步的，在本发明中所述血管模型的建立方式与所述血流模型和所述几何模型的建立方式基本相同，其差别点仅在于所述血管模型可同时包括待检血管段的形态和血流信息，故在本实施方式中，所述血管模型的具体建立方式与此不在赘述。

当然，在本设备中所述特定的生理参数包括年龄、性别、血压及身体质量指数等可直接获取的生理信息中的一个或多个。

进一步的，在本发明中所述处理器还可用于运行如下公式以计算获得所述血管压力差δp：

δp＝(c1v+c2v²+…+cmv^m)

*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]

其中，c1v+c2v²+…+cmv^m可为常数；

v为血流速度，为通过所述第二血流模型直接/间接获取；c1、c2、…、cm分别代表血流速度v的参数系数，所述参数系数包括血液粘度影响因素、血液湍流影响因素及粘滞系数等多个参数系数；进一步的，m为大于等于1的自然数，以分别代表不同参数系数对血流速度v的影响，以对压力差数值δp进行修正，保证血管压力差δp计算的准确性。优选的，在本发明中m的取值为2，且当m为2时，c1为因血液流动摩擦产生的参数系数，c2为血液湍流产生的参数系数。

所述α1、α2、…、αn分别为不同尺度下血管管腔的形态差异函数f1(x)、f2(x)、…、fn(x)的加权系数，其中，n为尺度为大于等于1的自然数；进一步的，所述加权系数的增加可进一步对形态差异函数f(x)进行修正，保证两横截面之间形态差异拟合计算的准确性。

当然，在本装置中所述特定的生理参数包括年龄、性别、血压及身体质量指数等可直接获取的生理信息中的一个或多个。

进一步的，在本发明中所述血管压力差数值的计算模型与所述横截面形态模型的建立方法基本一致，均是是基于多尺度计算方法而建立的。

需要指出的是，上述装置及功能模块仅仅为示例性的给出实现该技术方案的基本结构，而非唯一结构。

综上所述，本发明的基于特定的生理参数获取血管压力值的方法，通过在血管压力值的计算过程中引入个体特定的生理参数，一方面可保证血管压力值计算初期各类参数的提取和建立准确，另一方面，可对血管压力值计算过程中的血流模型进行修正，保证几何模型和血流模型建立的准确性，进一步保证通过所述几何模型和所述血流模型获取的相关参数准确，使得通过本发明的基于特定的生理参数获取血管压力值的方法计算得到的血管压力值可准确反映感兴趣区域的血流压力，且结果准确可靠。

以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围。