掌桥专利:专业的专利平台
掌桥专利
首页

一种主动脉血管体积测量方法、装置、设备及存储介质

文献发布时间:2024-04-18 19:52:40


一种主动脉血管体积测量方法、装置、设备及存储介质

技术领域

本发明属于医疗影像处理技术领域,具体涉及一种主动脉血管体积测量方法、装置、设备及存储介质。

背景技术

随着医疗影像技术的持续进步,诸如CT(Computed Tomography,即电子计算机断层扫描)图像分析与处理技术等在疾病诊断或治疗中发挥着越来越重要的作用。主动脉血管是人体内最粗大的动脉管,从心脏的左心室发出,向上向右再向下略呈弓状,再沿脊柱向下行,在胸腔和腹腔内分出很多较小的动脉。主动脉血管是向全身各部输送血液的主要导管。

目前,在对主动脉血管进行疾病诊断或治疗的过程中,主动脉血管体积是一个非常重要的参数。现有的主动脉血管体积测量技术主要包括如下步骤:先沿着主动脉血管轴状面对医疗影像进行分割切片,然后分别对经分割切片所得的每一个主动脉血管轴状面切片进行体积计算(即主动脉血管轴状面积与薄型切片的厚度之积),最后将所有主动脉血管轴状面切片的体积累加起来,得到整个主动脉血管的总体积。但是这种技术仅在血管较直时有效,而当血管出现如图1所示的弯曲或变形情况时,其体积测量结果的精度可能会受到影响(即由于沿主动脉血管轴状面的分割切片可能不会完全贴合呈弯曲状或异形状的主动脉血管结构),导致存在一定误差,进而可能误导医生做出错误的疾病诊断或治疗决策。

发明内容

本发明的目的是提供一种主动脉血管体积测量方法、装置、计算机设备及计算机可读存储介质,用以解决现有主动脉血管体积测量技术在面临血管弯曲或变形情况时存在体积测量结果精度低和误差大的问题。

为了实现上述目的,本发明采用以下技术方案:

第一方面,提供了一种主动脉血管体积测量方法,包括:

获取患者的主动脉血管医疗影像;

根据所述主动脉血管医疗影像,确定主动脉血管的且在三维空间中的血管中心线以及血管壁体外表面;

确定待测血管段的两个端点,其中,所述待测血管段为在所述主动脉血管上的任意管段,所述两个端点分别位于所述血管中心线上的不同位置;

针对在所述两个端点中的各个端点,分别作对应的且垂直于所述血管中心线的法平面;

确定由所述血管壁体外表面和与所述两个端点一一对应的两个法平面包围的空间,并将该空间作为所述待测血管段的所占空间;

在所述所占空间内随机选择一个初始体素,并基于所述初始体素,采用区域增长算法确定位于所述所占空间内的所有体素;

根据所述所有体素的体素总数目和每个体素的已知体积,计算得到所述所有体素的总体积,并将该总体积作为所述待测血管段的体积测量结果。

基于上述发明内容,提供了一种基于中心线法平面和区域增长算法实现主动脉血管体积测量的新方案,即先根据主动脉血管医疗影像确定主动脉血管的且在三维空间中的血管中心线以及血管壁体外表面,然后针对在待测血管段的两个端点中的各个端点,分别作对应的且垂直于所述血管中心线的法平面,再然后将由所述血管壁体外表面和与所述两个端点一一对应的两个法平面包围的空间作为所述待测血管段的所占空间,再然后基于初始体素,采用区域增长算法确定位于所述所占空间内的所有体素,最后计算所述所有体素的总体积,得到待测血管段体积测量结果,如此由于不再需要考虑分割切片是否贴合血管结构,因此相比较于现有主动脉血管体积测量技术可具有更高的精确度,特别是对于异形或弯曲的血管,误差可降低10%以上,进而可为医生提供更为精确的主动脉血管体积数据,利于医生做出正确的疾病诊断或治疗决策,便于实际应用和推广。

在一个可能的设计中,获取患者的主动脉血管医疗影像,包括:

通过电子计算机断层扫描机扫描获取患者的主动脉血管医疗三维影像。

在一个可能的设计中,根据所述主动脉血管医疗影像,确定主动脉血管的且在三维空间中的血管中心线以及血管壁体外表面,包括:

对所述主动脉血管医疗影像进行沿着主动脉血管轴状面的分割切片处理,得到连续的多个主动脉血管轴状面切片图像;

针对在所述多个主动脉血管轴状面切片图像中的各个主动脉血管轴状面切片图像,先根据预设像素阈值将对应图像处理成对应的二值化图像,然后对该二值化图像的连通域进行基于Canny算法的边缘轮廓提取处理,得到对应的连通域边缘轮廓,并基于几何中心法确定该连通域边缘轮廓的中心,得到对应的连通域中心点,其中,所述预设像素阈值为与血管颜色对应的像素值;

根据所述多个主动脉血管轴状面切片图像的连续性,依次连接所有的所述连通域中心点,得到主动脉血管在三维空间中的且以所述连通域中心点为离散点的血管中心线,以及还依次连接所有的所述连通域边缘轮廓,得到所述主动脉血管在所述三维空间中的且以所述连通域边缘轮廓为离散管圈的血管壁体外表面。

在一个可能的设计中,在得到主动脉血管在三维空间中的且以所述连通域中心点为离散点的血管中心线之后,所述方法还包括:

对所述血管中心线进行平滑处理或和重采样处理,得到新的血管中心线。

在一个可能的设计中,针对在所述两个端点中的某个端点,作对应的且垂直于所述血管中心线的法平面,包括:

根据所述某个端点在所述三维空间中的三维坐标(x

根据所述法向量(dx

dx

式中,(x,y,z)表示在所述法平面上的任意一点在所述三维空间中的三维坐标。

在一个可能的设计中,基于所述初始体素,采用区域增长算法确定位于所述所占空间内的所有体素,包括有如下步骤S61~S65:

S61.将所述初始体素添加到初始为空的当前区域中,然后将所述初始体素作为当前体素,执行步骤S62;

S62.判断与当前体素相邻的且未添加到所述当前区域中的某个其它体素是否位于所述所占空间内,若是,则将该某个其它体素添加到所述当前区域中,然后执行步骤S63,否则直接执行步骤S63;

S63.判断是否还存在与当前体素相邻的、位于所述所占空间内且未添加到所述当前区域中的其它体素,若存在,则返回执行步骤S62,否则执行步骤S64;

S64.将基于当前体素添加到所述当前区域中的各个其它体素分别作为新的当前体素,然后返回执行步骤S62;

S65.在无任何新体素添加到所述当前区域后,将已添加到所述当前区域中的所有体素确定为位于所述所占空间内的所有体素。

在一个可能的设计中,当所述主动脉血管医疗影像为通过电子计算机断层扫描机扫描获取的主动脉血管医疗三维影像时,根据所述所有体素的体素总数目和每个体素的已知体积,计算得到所述所有体素的总体积,包括:

根据所述电子计算机断层扫描机的分辨率,确定每个体素的体积大小,其中,所述每个体素的体积大小与所述分辨率负相关;

根据所述所有体素的体素总数目和所述每个体素的体积大小,按照如下公式计算得到所述所有体素的总体积V:

V=N×V

式中,N表示所述所有体素的体素总数目,V

第二方面,提供了一种主动脉血管体积测量装置,包括有医疗影像获取模块、血管线面确定模块、管段端点确定模块、法平面生成模块、管段空间确定模块、空间内体素确定模块和管段体积计算模块;

所述医疗影像获取模块,用于获取患者的主动脉血管医疗影像;

所述血管线面确定模块,通信连接所述医疗影像获取模块,用于根据所述主动脉血管医疗影像,确定主动脉血管的且在三维空间中的血管中心线以及血管壁体外表面;

所述管段端点确定模块,通信连接所述血管线面确定模块,用于确定待测血管段的两个端点,其中,所述待测血管段为在所述主动脉血管上的任意管段,所述两个端点分别位于所述血管中心线上的不同位置;

所述法平面生成模块,通信连接所述管段端点确定模块,用于针对在所述两个端点中的各个端点,分别作对应的且垂直于所述血管中心线的法平面;

所述管段空间确定模块,分别通信连接所述血管线面确定模块和所述管段端点确定模块,用于确定由所述血管壁体外表面和与所述两个端点一一对应的两个法平面包围的空间,并将该空间作为所述待测血管段的所占空间;

所述空间内体素确定模块,通信连接所述管段空间确定模块,用于在所述所占空间内随机选择一个初始体素,并基于所述初始体素,采用区域增长算法确定位于所述所占空间内的所有体素;

所述管段体积计算模块,通信连接所述空间内体素确定模块,用于根据所述所有体素的体素总数目和每个体素的已知体积,计算得到所述所有体素的总体积,并将该总体积作为所述待测血管段的体积测量结果。

第三方面,本发明提供了一种计算机设备,包括有依次通信连接的存储器、处理器和收发器,其中,所述存储器用于存储计算机程序,所述收发器用于收发消息,所述处理器用于读取所述计算机程序,执行如第一方面或第一方面中任意可能设计所述的主动脉血管体积测量方法。

第四方面,本发明提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有指令,当所述指令在计算机上运行时,执行如第一方面或第一方面中任意可能设计所述的主动脉血管体积测量方法。

第五方面,本发明提供了一种包含指令的计算机程序产品,当所述指令在计算机上运行时,使所述计算机执行如第一方面或第一方面中任意可能设计所述的主动脉血管体积测量方法。

上述方案的有益效果:

(1)本发明创造性提供了一种基于中心线法平面和区域增长算法实现主动脉血管体积测量的新方案,即先根据主动脉血管医疗影像确定主动脉血管的且在三维空间中的血管中心线以及血管壁体外表面,然后针对在待测血管段的两个端点中的各个端点,分别作对应的且垂直于所述血管中心线的法平面,再然后将由所述血管壁体外表面和与所述两个端点一一对应的两个法平面包围的空间作为所述待测血管段的所占空间,再然后基于初始体素,采用区域增长算法确定位于所述所占空间内的所有体素,最后计算所述所有体素的总体积,得到待测血管段体积测量结果,如此由于不再需要考虑分割切片是否贴合血管结构,因此相比较于现有主动脉血管体积测量技术可具有更高的精确度,特别是对于异形或弯曲的血管,误差可降低10%以上,进而可为医生提供更为精确的主动脉血管体积数据,利于医生做出正确的疾病诊断或治疗决策,便于实际应用和推广。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。

图1为在实际情况中呈弯曲状及异形状的主动脉血管结构的示例图。

图2为本申请实施例提供的主动脉血管体积测量方法的流程示意图。

图3为本申请实施例提供的血管中心线及血管壁体外表面与两端法平面的位置关系示例图。

图4为本申请实施例提供的主动脉血管体积测量装置的结构示意图。

图5为本申请实施例提供的计算机设备的结构示意图。

具体实施方式

为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将结合附图和实施例或现有技术的描述对本发明作简单地介绍,显而易见地,下面关于附图结构的描述仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在此需要说明的是,对于这些实施例方式的说明用于帮助理解本发明,但并不构成对本发明的限定。

应当理解,尽管本文可能使用术语第一和第二等等来描述各种对象,但是这些对象不应当受到这些术语的限制。这些术语仅用于区分一个对象和另一个对象。例如可以将第一对象称作第二对象,并且类似地可以将第二对象称作第一对象,同时不脱离本发明的示例实施例的范围。

应当理解,对于本文中可能出现的术语“和/或”,其仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A、单独存在B或者同时存在A和B等三种情况;又例如,A、B和/或C,可以表示存在A、B和C中的任意一种或他们的任意组合;对于本文中可能出现的术语“或和”,其是描述另一种关联对象关系,表示可以存在两种关系,例如,A或和B,可以表示:单独存在A或者同时存在A和B等两种情况;另外,对于本文中可能出现的字符“/”,一般表示前后关联对象是一种“或”关系。

实施例:

如图2~3所示,本实施例第一方面提供的所述主动脉血管体积测量方法,可以但不限于由具有一定计算资源的计算机设备执行,例如由平台服务器、CT扫描仪的数据处理端、个人计算机(Personal Computer,PC,指一种大小、价格和性能适用于个人使用的多用途计算机;台式机、笔记本电脑到小型笔记本电脑和平板电脑以及超级本等都属于个人计算机)、智能手机、个人数字助理(Personal Digital Assistant,PDA)或可穿戴设备等电子设备执行。如图2所示,所述主动脉血管体积测量方法,可以但不限于包括有如下步骤S1~S7。

S1.获取患者的主动脉血管医疗影像。

在所述步骤S1中,所述主动脉血管医疗影像即为针对主动脉血管探测而得的医疗影像,具体的,获取患者的主动脉血管医疗影像,包括但不限于有如下方式:通过电子计算机断层扫描机扫描获取患者的主动脉血管医疗三维影像。

S2.根据所述主动脉血管医疗影像,确定主动脉血管的且在三维空间中的血管中心线以及血管壁体外表面。

在所述步骤S2中,所述血管中心线即反映了所述主动脉血管的整体路径及方向,所述血管壁体外表面即反映了所述主动脉血管的外周边界。具体的,根据所述主动脉血管医疗影像,确定主动脉血管的且在三维空间中的血管中心线以及血管壁体外表面,包括但不限于有如下步骤S21~S23。

S21.对所述主动脉血管医疗影像进行沿着主动脉血管轴状面的分割切片处理,得到连续的多个主动脉血管轴状面切片图像。

在所述步骤S21中,所述分割切片处理的具体方式可参照现有主动脉血管体积测量技术常规得到,于此不再赘述。

S22.针对在所述多个主动脉血管轴状面切片图像中的各个主动脉血管轴状面切片图像,先根据预设像素阈值将对应图像处理成对应的二值化图像,然后对该二值化图像的连通域进行基于Canny算法的边缘轮廓提取处理,得到对应的连通域边缘轮廓,并基于几何中心法确定该连通域边缘轮廓的中心,得到对应的连通域中心点,其中,所述预设像素阈值为与血管颜色对应的像素值。

在所述步骤S22中,所述预设像素阈值可以具体为一个区间值,并在二值化图像处理时,可具体按照如下方式进行:若在所述主动脉血管轴状面切片图像中的某个像素点的像素值属于该区间值,则将该某个像素点在二值化图像中的像素值设置为255,反之设置为0,如此可得到二值分别为0和255的二值化图像。在所述二值化图像中的连通域即为血管轴状面区域,因此可以通过现有的Canny算法(其是由John F.Canny于1986年开发出来的一个多级边缘检测算法)提取得到作为血管轴状面边界的所述连通域边缘轮廓。此外,所述几何中心法是一种常用的不规则物体重心确定方法,因此可以准确得到作为血管中心点的所述连通域中心点。

S23.根据所述多个主动脉血管轴状面切片图像的连续性,依次连接所有的所述连通域中心点,得到主动脉血管在三维空间中的且以所述连通域中心点为离散点的血管中心线,以及还依次连接所有的所述连通域边缘轮廓,得到所述主动脉血管在所述三维空间中的且以所述连通域边缘轮廓为离散管圈的血管壁体外表面。

在所述步骤S23中,考虑所述血管中心线当前所具有的离散性,可能存在噪声、不规则性以及点分布不均匀等缺陷,会影响后续的测量结果精度,因此为了消除这些缺陷,优选的,在得到主动脉血管在三维空间中的且以所述连通域中心点为离散点的血管中心线之后,所述方法还包括但不限于有如下步骤:对所述血管中心线进行平滑处理或和重采样处理,得到新的血管中心线。前述平滑处理或和重采样处理是指基于现有的平滑算法对所述血管中心线进行平滑处理(主要消除可能的噪声或不规则性),得到平滑的且新的血管中心线,或者是指先基于现有的平滑算法对所述血管中心线进行平滑处理,再基于现有的重采样算法对平滑处理结果进行重采样处理,得到点均匀分布的且新的血管中心线。此外,类似的,还可以基于现有面平滑算法对所述血管壁体外表面进行平滑处理,得到平滑的且新的血管壁体外表面。

S3.确定待测血管段的两个端点,其中,所述待测血管段为在所述主动脉血管上的任意管段,所述两个端点分别位于所述血管中心线上的不同位置。

在所述步骤S3中,所述两个端点用于限定所述待测血管段的两端,它们可以具体但不限于通过常规的人机交互方式来由医生选择确定。所述待测血管段可以是在所述主动脉血管上的部分管段,也可以是整个所述主动脉血管。此外,还可以根据需要将所述主动脉血管拆分成多个所述待测血管段,最后汇总多个所述待测血管段的体积,得到所述主动脉血管的总体积。

S4.针对在所述两个端点中的各个端点,分别作对应的且垂直于所述血管中心线的法平面。

在所述步骤S4中,具体的,针对在所述两个端点中的某个端点,作对应的且垂直于所述血管中心线的法平面,包括但不限于有如下步骤S41~S42。

S41.根据所述某个端点在所述三维空间中的三维坐标(x

在所述步骤S41中,所述法向量(dx

S42.根据所述法向量(dx

dx

式中,(x,y,z)表示在所述法平面上的任意一点在所述三维空间中的三维坐标。

S5.确定由所述血管壁体外表面和与所述两个端点一一对应的两个法平面包围的空间,并将该空间作为所述待测血管段的所占空间。

在所述步骤S5中,举例如图3所示,由血管壁体外表面S0、第一法平面S1和第二法平面S2包围的空间即为所述待测血管段的所占空间。

S6.在所述所占空间内随机选择一个初始体素,并基于所述初始体素,采用区域增长算法确定位于所述所占空间内的所有体素。

在所述步骤S6中,体素是体积元素(Volume Pixel)的简称,具体是数字数据于三维空间分割上的最小单位。所述区域增长算法的基本思想是将具有相似性质的像素点合并到一起,即对每一个区域要先指定一个种子点作为生长的起点,然后将位于种子点周围邻域内的像素点与种子点进行对比,将具有相似性质的点合并起来继续向外生长,直到没有满足条件的像素点被包括进来为止。具体的,基于所述初始体素,采用区域增长算法确定位于所述所占空间内的所有体素,包括有如下步骤S61~S65。

S61.将所述初始体素添加到初始为空的当前区域中,然后将所述初始体素作为当前体素,执行步骤S62。

S62.判断与当前体素相邻的且未添加到所述当前区域中的某个其它体素是否位于所述所占空间内,若是,则将该某个其它体素添加到所述当前区域中,然后执行步骤S63,否则直接执行步骤S63。

S63.判断是否还存在与当前体素相邻的、位于所述所占空间内且未添加到所述当前区域中的其它体素,若存在,则返回执行步骤S62,否则执行步骤S64。

S64.将基于当前体素添加到所述当前区域中的各个其它体素分别作为新的当前体素,然后返回执行步骤S62。

S65.在无任何新体素添加到所述当前区域后,将已添加到所述当前区域中的所有体素确定为位于所述所占空间内的所有体素。

S7.根据所述所有体素的体素总数目和每个体素的已知体积,计算得到所述所有体素的总体积,并将该总体积作为所述待测血管段的体积测量结果。

在所述步骤S7中,所述每个体素即代表了在所述三维空间中的一个小立方体,因此其体积大小可以提前已知,例如若一个体素的长宽高分别为1mm×1mm×1mm,则它的体积大小可以已知为1mm

S71.根据所述电子计算机断层扫描机的扫描分辨率,确定每个体素的体积大小,其中,所述每个体素的体积大小与所述扫描分辨率负相关。

S72.根据所述所有体素的体素总数目和所述每个体素的体积大小,按照如下公式计算得到所述所有体素的总体积V:

V=N×V

式中,N表示所述所有体素的体素总数目,V

由此基于前述步骤S1~S7所描述的主动脉血管体积测量方法,提供了一种基于中心线法平面和区域增长算法实现主动脉血管体积测量的新方案,即先根据主动脉血管医疗影像确定主动脉血管的且在三维空间中的血管中心线以及血管壁体外表面,然后针对在待测血管段的两个端点中的各个端点,分别作对应的且垂直于所述血管中心线的法平面,再然后将由所述血管壁体外表面和与所述两个端点一一对应的两个法平面包围的空间作为所述待测血管段的所占空间,再然后基于初始体素,采用区域增长算法确定位于所述所占空间内的所有体素,最后计算所述所有体素的总体积,得到待测血管段体积测量结果,如此由于不再需要考虑分割切片是否贴合血管结构,因此相比较于现有主动脉血管体积测量技术可具有更高的精确度,特别是对于异形或弯曲的血管,误差可降低10%以上,进而可为医生提供更为精确的主动脉血管体积数据,利于医生做出正确的疾病诊断或治疗决策,便于实际应用和推广。

如图4所示,本实施例第二方面提供了一种实现第一方面所述的主动脉血管体积测量方法的虚拟装置,包括有医疗影像获取模块、血管线面确定模块、管段端点确定模块、法平面生成模块、管段空间确定模块、空间内体素确定模块和管段体积计算模块;

所述医疗影像获取模块,用于获取患者的主动脉血管医疗影像;

所述血管线面确定模块,通信连接所述医疗影像获取模块,用于根据所述主动脉血管医疗影像,确定主动脉血管的且在三维空间中的血管中心线以及血管壁体外表面;

所述管段端点确定模块,通信连接所述血管线面确定模块,用于确定待测血管段的两个端点,其中,所述待测血管段为在所述主动脉血管上的任意管段,所述两个端点分别位于所述血管中心线上的不同位置;

所述法平面生成模块,通信连接所述管段端点确定模块,用于针对在所述两个端点中的各个端点,分别作对应的且垂直于所述血管中心线的法平面;

所述管段空间确定模块,分别通信连接所述血管线面确定模块和所述管段端点确定模块,用于确定由所述血管壁体外表面和与所述两个端点一一对应的两个法平面包围的空间,并将该空间作为所述待测血管段的所占空间;

所述空间内体素确定模块,通信连接所述管段空间确定模块,用于在所述所占空间内随机选择一个初始体素,并基于所述初始体素,采用区域增长算法确定位于所述所占空间内的所有体素;

所述管段体积计算模块,通信连接所述空间内体素确定模块,用于根据所述所有体素的体素总数目和每个体素的已知体积,计算得到所述所有体素的总体积,并将该总体积作为所述待测血管段的体积测量结果。

本实施例第二方面提供的前述装置的工作过程、工作细节和技术效果,可以参见第一方面所述的主动脉血管体积测量方法,于此不再赘述。

如图5所示,本实施例第三方面提供了一种执行如第一方面所述的主动脉血管体积测量方法的计算机设备,包括有依次通信连接的存储器、处理器和收发器,其中,所述存储器用于存储计算机程序,所述收发器用于收发消息,所述处理器用于读取所述计算机程序,执行如第一方面所述的主动脉血管体积测量方法。具体举例的,所述存储器可以但不限于包括随机存取存储器(Random-Access Memory,RAM)、只读存储器(Read-Only Memory,ROM)、闪存(Flash Memory)、先进先出存储器(First Input First Output,FIFO)和/或先进后出存储器(First Input Last Output,FILO)等等;所述处理器可以但不限于采用型号为STM32F105系列的微处理器。此外,所述计算机设备还可以但不限于包括有电源模块、显示屏和其它必要的部件。

本实施例第三方面提供的前述计算机设备的工作过程、工作细节和技术效果,可以参见第一方面所述的主动脉血管体积测量方法,于此不再赘述。

本实施例第四方面提供了一种存储包含如第一方面所述的主动脉血管体积测量方法的指令的计算机可读存储介质,即所述计算机可读存储介质上存储有指令,当所述指令在计算机上运行时,执行如第一方面所述的主动脉血管体积测量方法。其中,所述计算机可读存储介质是指存储数据的载体,可以但不限于包括软盘、光盘、硬盘、闪存、优盘和/或记忆棒(Memory Stick)等计算机可读存储介质,所述计算机可以是通用计算机、专用计算机、计算机网络或者其他可编程装置。

本实施例第四方面提供的前述计算机可读存储介质的工作过程、工作细节和技术效果,可以参见如第一方面所述的主动脉血管体积测量方法,于此不再赘述。

本实施例第五方面提供了一种包含指令的计算机程序产品,当所述指令在计算机上运行时,使所述计算机执行如第一方面所述的主动脉血管体积测量方法。其中,所述计算机可以是通用计算机、专用计算机、计算机网络或者其他可编程装置。

最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

技术分类

06120116331730