一种血糖波动评价方法和评价装置与流程

本发明涉及血糖评价
技术领域：
，尤其涉及一种血糖波动评价方法和评价装置。
背景技术：
：糖尿病是由于体内胰岛素分泌缺陷或其生物作用受损所引起的代谢紊乱性疾病。随着人口老龄化和饮食习惯的改变，糖尿病的发病率呈急剧上升趋势。迄今为止，糖尿病尚无切实可行的根治方法，糖尿病患者被建议通过连续监测血糖浓度变化来及时调整口服降糖药物和胰岛素的用量或调整饮食运动等，进而控制血糖的波动水平，预防或减轻并发症。连续血糖监测(cgm)是一种通过葡萄糖传感器连续监测皮下组织间液葡萄糖浓度的技术。一般而言，cgm可连续获取患者24小时的血糖值，这与传统血糖监测方法相比具有明显的优势。现有研究表明，血糖的波动水平与低血糖、糖尿病并发症之间存在相关性。因此，基于cgm数据的血糖波动评估是目前的研究热点。传统的血糖波动评估方法包括计算cgm的血糖标准差(sdbg)，血糖变异系数(cv)，最大血糖波动幅度(lage)，平均血糖波动幅度(mage)等。具体地，现有技术包括以下技术方案：1)、血糖标准差(sdbg)计算方法：其中，g为测量血糖值，为测量血糖值的平均值，n为血糖测量总数2)、血糖变异系数(cv)计算方法：其中，sdbg为血糖标准差，为测量血糖值的平均值3)、最大血糖波动幅度(lage)计算方法：lage＝gmax-gmin其中gmax为测量最大血糖值，gmin为测量最小血糖值。4)、平均血糖波动幅度(mage)计算方法：首先计算患者24小时内的血糖标准差(sdbg)；其次，计算血糖的波动幅度，根据波动幅度确定血糖的波动方向(上升或下降)；再次，根据血糖波动方向确定有效血糖波动幅度；最后，求取有效血糖波动的平均值。然而，上述的计算cgm数据的血糖标准差，血糖变异系数等，只能计算出血糖波动的幅度，而忽略了血糖波动的频率。因而，目前的血糖波动参数很难同时表征血糖波动幅度与血糖波动频率，可能会导致部分信息丢失，无法描述血糖的真正波动情况。技术实现要素：本发明的目的在于克服上述现有技术的缺陷，提供一种血糖波动评价方法和评价装置，是针对cgm数据基于梯度变异的血糖波动计算方法，所提出的梯度变异参数可以同时表征血糖波动幅度与血糖波动频率，为研究血糖波动提供一种新的思路。根据本发明的第一方面，提供血糖波动评价方法。该方法包括以下步骤：将以预设频率采集的多个监测周期的连续血糖监测数据表征为矩阵形式，其中，矩阵的行数对应监测周期数，矩阵的列数对应一个监测周期内采集的连续血糖监测数据的数目，矩阵中各元素是所采集的连续血糖监测数据值；对矩阵进行梯度计算，获得横向梯度和纵向梯度，其中所述横向梯度表征各监测周期内的连续血糖监测数据的波动情况，所述纵向梯度表征各监测周期之间的连续血糖监测数据的波动情况；基于所述横向梯度和所述纵向梯度计算梯度变异值，进而利用该梯度变异值评价血糖波动幅度和波动频率。根据本发明的第二方面，提供一种血糖波动评价装置。该装置包括：数据表征单元，其用于将以预设频率采集的多个监测周期的连续血糖监测数据表征为矩阵形式，其中，矩阵的行数对应监测周期数，矩阵的列数对应一个监测周期内采集的连续血糖监测数据的数目，矩阵中各元素是所采集的连续血糖监测数据值；梯度计算单元，其用于对矩阵进行计算梯度，获得横向梯度和纵向梯度，其中所述横向梯度表征各监测周期内的连续血糖监测数据的波动情况，所述纵向梯度表征每个监测周期之间的连续血糖监测数据的波动情况；血糖波动评价单元，其用于基于所述横向梯度和所述纵向梯度计算梯度变异值，进而利用该梯度变异值评价血糖波动幅度和波动频率。在一个实施例中，所述血糖波动评价单元采用以下公式计算所述梯度变异值：其中，fx(i,j)为横向梯度fx的第(i,j)个元素，fy(i,j)为纵向梯度fy的第(i,j)个元素，δti为相邻两个矩阵元素的横向时间，δtj为相邻两个矩阵元素的纵向时间，origin(i,j)为矩阵origin的第(i,j)个元素，m为矩阵origin的行数，n为矩阵origin的列数。与现有技术相比，本发明的优点在于：提出一种基于梯度变异的血糖波动评价方法和评价装置，能够同时反映出血糖的波动幅度与波动频率，通过实验验证这种基于梯度变异的血糖评价方法在预测低血糖时具有较大优势，为临床上分析血糖的波动提供了一种新的解决方案。附图说明以下附图仅对本发明作示意性的说明和解释，并不用于限定本发明的范围，其中：图1是根据本发明一个实施例的血糖评价方法的流程图；图2是根据本发明一个实施例的周期为20分钟的cgm值；图3是根据本发明一个实施例的周期为40分钟的cgm值；图4是根据本发明一个实施例的周期为60分钟的cgm值；图5是根据本发明一个实施例的不同参数下的roc曲线。具体实施方式为了使本发明的目的、技术方案、设计方法及优点更加清楚明了，以下结合附图通过具体实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅用于解释本发明，并不用于限定本发明。在本文示出和讨论的所有例子中，任何具体值应被解释为仅仅是示例性的，而不是作为限制。因此，示例性实施例的其它例子可以具有不同的值。对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论，但在适当情况下，所述技术、方法和设备应当被视为说明书的一部分。参见图1所示，本发明提出的血糖波动评价方法包括以下步骤：步骤s110，获取并处理cgm数据。具体地，对于连续血糖监测数据(cgm数据)的获取，可采用目前市场上的cgm设备，例如血糖监测频率为5分钟，1小时内可获取12个血糖值，1天内(24小时)则可获取288个血糖值。如果由于cgm设备暂时故障导致某部分数据缺失，对缺失的数据开入采用线性差值法进行拟合。步骤s120，对cgm数据进行矩阵化处理，将以预设频率采集的多个监测周期的cgm数据表征为矩阵形式。具体地，将cgm数据表述为以下矩阵：其中(a0101......a0112)表示第一个小时(即第一个监测周期)的cgm值，共有12个值；(a0201......a0212)表示第二个小时(即第二个监测周期)的cgm值，共有12个值，以此类推。因此，矩阵origin的默认行数m＝24，默认列数n＝12。矩阵的行数m和列数n可根据实际需要进行设置。步骤s130，对矩阵进行梯度计算，获得横向梯度和纵向梯度，横向梯度表征各监测周期内的cgm数据的波动情况，纵向梯度表征各监测周期之间的cgm数据的波动情况。具体地，对矩阵进行梯度计算，包括横向梯度和纵向梯度，即gradient(origin)＝[fx,fy]，其中fx为横向梯度，fy为纵向梯度，分别表示为：步骤s140，根据横向梯度和纵向梯度计算梯度变异值。本发明实施例提出用梯度变异值gv反映血糖波动的新参数，该值计算方法如下：其中fx(i,j)为横向梯度fx的第(i,j)个元素，fy(i,j)为纵向梯度fy的第(i,j)个元素，δti为相邻两个矩阵元素的横向时间，即δti为5分种(即约0.083小时)，δtj为相邻两个矩阵元素的纵向时间，即δti为1小时，origin(i,j)为矩阵origin的第(i,j)个元素，m为矩阵origin的行数，n为矩阵origin的列数。从梯度变异值gv表达式(4)可知，其包含了血糖5分钟内的波动，又能同时反映1小时内的波动，还可以反映血糖的波动幅度。步骤s150，分析梯度变异值与其他血糖波动参数的相关性。进一步地，分析梯度变异值gv与其他血糖波动参数的相关性。具体地，通过获取多个患者的cgm数据，分别计算每个患者的梯度变异值gv，标准差(sdbg)，血糖变异系数(cv)，最大血糖波动幅度(lage)，平均血糖波动幅度(mage)，采用相关性分析方法来分析所提出的梯度变异值gv与其他血糖波动参数的相关性。表1是相关性分析结果。由表1可知，梯度变异值gv与标准差sdbg，血糖变异系数cv，最大血糖波动幅度lage，平均血糖波动幅度mage均具有强的相关性，论证了所提出的梯度变异值gv可以反映出血糖的波动。表1梯度变异值gv与其他血糖波动参数的相关性相关系数p值梯度变异值gv与血糖标准差sdbg0.51380梯度变异值gv与血糖变异系数cv0.70160梯度变异值gv与最大血糖波动幅度lage0.58350梯度变异值gv与平均血糖波动幅度mage0.62710步骤s160，基于梯度变异值分析不同频率的血糖波动。具体地，采用数值模拟方法，模拟了三种不同频率的cgm值，参见图2至图4所示，其中横坐标表示各个周期，纵坐标表示血糖值，图2的cgm值周期为20分钟，图3的cgm值的周期为40分钟，图4的cgm值的周期为60分钟。下表2给出了不同血糖波动参数的计算结果。由表2可知，对于不同周期的cgm值，其血糖标准差sdbg，血糖变异系数cv，最大血糖波动幅度lage，平均血糖波动幅度mage均无任何变化，而所提出的梯度变异值gv随着cgm值的周期变化而变化，从而论证了梯度变异值gv可以同时反映出血糖波动幅度与波动频率。表2不同周期的各种血糖波动参数计算结果参数周期20分钟周期40分钟周期60分钟梯度变异值gv5.50922.26741.1018血糖标准差sdbg3.10533.10533.1053血糖差异系数cv0.44360.44360.4436最大血糖波动幅度lage6.20006.20006.2000平均血糖波动幅度mage5.6725.6725.672步骤s170，分析梯度变异值与低血糖预测的关联性。具体地，通过利用cgm设备分别获取多个患者在发生低血糖事件前与正常情况下的cgm数据，根据cgm数据分别计算其梯度变异值gv，血糖标准差sdbg，血糖变异系数cv，最大血糖波动幅度lage，平均血糖波动幅度mage。通过采用roc曲线(受试者工作特征曲线)分析不同参数在预测低血糖事件时的敏感度与特异度。图5是以1-特异度为横坐标，敏感度为纵坐标的roc曲线。从图5可以看出，梯度变异值gv的工作面积最大，验证了梯度变异值gv在预测低血糖方面具有更大的优势。相应地，本发明还提供一种血糖评价装置，该装置用于实现上述方法的一个方面或多个方面，在此不在具体描述。例如，该装置包括：数据表征单元，其用于将以预设频率采集的多个监测周期的连续血糖监测数据表征为矩阵形式，其中，矩阵的行数对应监测周期数，矩阵的列数对应一个监测周期内采集的连续血糖监测数据的数目，矩阵中各元素是所采集的连续血糖监测数据值；梯度计算单元，其用于对矩阵进行计算梯度，获得横向梯度和纵向梯度，其中所述横向梯度表征各监测周期内的连续血糖监测数据的波动情况，所述纵向梯度表征每个监测周期之间的连续血糖监测数据的波动情况；血糖波动评价单元，其用于基于所述横向梯度和所述纵向梯度计算梯度变异值，进而利用该梯度变异值评价血糖波动幅度和波动频率。本发明装置中的各单元可通过逻辑硬件结构或处理器实现。综上所述，本发明提供一种基于cgm数据评价血糖波动的方法和装置，通过矩阵运算等得出梯度变异值，与其他血糖波动参数相比，该参数能够同时反映出血糖波动的幅度以及波动频率，并且本发明所提出的梯度变异值在预测低血糖事件方面具有更大的优势。需要说明的是，虽然上文按照特定顺序描述了各个步骤，但是并不意味着必须按照上述特定顺序来执行各个步骤，实际上，这些步骤中的一些可以并发执行，甚至改变顺序，只要能够实现所需要的功能即可。本发明可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本发明的各个方面的计算机可读程序指令。计算机可读存储介质可以是保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以包括但不限于电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：便携式计算机盘、硬盘、随机存取存储器(ram)、只读存储器(rom)、可擦式可编程只读存储器(eprom或闪存)、静态随机存取存储器(sram)、便携式压缩盘只读存储器(cd-rom)、数字多功能盘(dvd)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本
技术领域：
的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进，或者使本
技术领域：
的其它普通技术人员能理解本文披露的各实施例。当前第1页12