申请/专利权人：吉林大学;林赐云

申请日：2022-06-12

公开（公告）日：2024-04-26

公开（公告）号：CN115192009B

主分类号：A61B5/1455

分类号：A61B5/1455

优先权：

专利状态码：有效-授权

法律状态：2024.04.26#授权;2022.11.04#实质审查的生效;2022.10.18#公开

摘要：本发明属于COVID‑19病毒的诊断技术领域，具体涉及一种基于动静脉反射光谱和高光谱图像的SaO2检测方法，本发明一种是非接触测量系统，考虑到COVID‑19病毒的快速筛查，非接触的检测方式可以提升检测效率；具有远程识别含氧动脉和脱氧静脉以测量血氧饱和度的能力；本发明是一种区别于传统检测手段的方法，采用高光谱技术进行医疗检测，可靠性高，检测成本低；高光谱检测系统响应时间迅速13至16秒，对于病毒的前期快速筛查具有良好效益。

主权项：1.一种基于动静脉反射光谱和高光谱图像测量SaO2的COVID-19筛查方法，其特征在于：S1、利用手指脉搏血氧仪采集实验对象的SaO2值作为评估高光谱成像系统检测精度的标准；S2、用光谱范围为348-950nm的多色光源照射实验对象的前臂和手腕之间的区域，每人扫描4次，计算平均值；S3、获取的高光谱图像是三维数据立方体，从光谱维度出发，首先利用最小噪声分离变换去除光谱信息噪声，平滑光谱曲线；然后利用特征提取算法分别找出含氧动脉和去氧静脉的特征谱段；S4、利用归一化处理算法对S3中特征谱段下的二维图像进行图像增强以消除伪影和噪声；S5、将S4中增强后的图像采用灰度直方图均衡化，k-means聚类算法进行图像分割，突出含氧动脉和去氧静脉的位置，最后根据二者比率计算SaO2大小；进一步的，所述步骤S1获得检测样本的具体步骤如下：S101、测量血氧前使用酒精清洁血氧探头表面，注意不能将探头全部浸入酒精中；S102、打开脉搏血氧仪，正确设置脉搏血氧仪相关参数；S103、将血氧探头连接线的一端插入脉搏血氧仪顶端的SaO2探头接口，另一端与实验对象的手指相连；S104、为减小误差，每位实验对象进行多次测量，记录数据，最后关闭脉搏血氧仪；所述步骤S2搭建实验环境，获取高光谱图像数据的具体步骤如下：S201、根据氧合血红蛋白的结构特性选取348-950nm的光源，以及光谱分辨率为4nm的高光谱相机进行实验；S202、固定高光谱相机高度为距离测试者手臂20cm处，并调节焦距和多色光光源的空间位置直到感兴趣区域的预览结果清晰；S203、对感兴趣区域进行凝视拍摄，每人连续拍摄四次，得到四个高光谱图像，求出这四个输出的平均值；进一步讲，S3具体步骤如下：S301、利用最小噪声分离变换去除光谱信息噪声，平滑光谱曲线；最小噪声分离变换算法MNF首先对数据噪声进行估计，得到噪声协方差矩阵，用于分离和调节数据中的噪声，去除数据噪声的相关性，并对噪声进行归一化处理，得到方差为1且不相关的噪声数据；然后，分别对噪声数据和信号数据进行PCA变换，并求PCA变换后数据的信噪比；最后，由大到小排列信噪比，得到数据的各主成分信息；MNF根据信噪比大小来确定主成分信息，减少了噪声的影响；通过高通滤波器模板对数据X进行分离，得到噪声XN和信号XS，则有：X＝XN+XS通过最大化变换后数据的信噪比，即最大化信号协方差与噪声协方差的比值，得到变换矩阵V； 式中，YN和YS分别为变换后数据的噪声和信号，CN和CS分别为数据的噪声和信号的协方差，此优化问题可等价为： 式中，C表示数据的总体协方差，C＝CovX，C＝CN+CS；根据拉格朗日乘子法最优解为：CV＝λCNV根据上式从大到小排列特征值，取前d个特征值对应的特征向量，可得到转换矩阵：V＝[v1,v2,···,vd]S302、利用特征选择算法分别寻找含氧动脉和去氧静脉的特征谱段；实验所需高光谱相机波长范围为400-1000nm，光谱分辨率4nm，128个谱段，将这些二维离散点数据用基于决策树的特征选择算法进行处理以选取主特征；决策树的基本思想是，从根节点开始，选取最优特征将实例划分为不同子集，若子集中实例能被正确分类，则构建叶结点，不再继续划分；否则，对子集再次选最优特征进行划分，如此递归划分，直至所有实例能被正确分类，或达到停止条件；特征选择度量是一种决策树的分裂准则，把给定类标记的训练元组的训练集D划分为单独类的启发方式；这里介绍三种常用的特征选择度量：信息增益，增益率，基尼指数；1、信息增益ID3熵是表示随机变量不确定性的度量，熵只依赖于随机变量的分布，而与随机变量的取值无关，并且熵越大，随机变量的不确定性就越大；信息增益定义：特征A对训练集D的信息增益gD,A，定义为集合D的熵HD与特征A给定条件下D的条件熵HD|A的差，即：gD,A＝HD-HD|A信息增益表示得到特征X的信息而使得类Y的信息不确定性减少的程度； 其中|D|表示其样本容量，K表示类的个数,|Ck|表示的属于类Ck的个数，设特征A有n个不同的取值，根据特征A的取值将D划分为n个子集Di,|Di|为Di的样本个数，记子集Di中属于类Ck的样本集合为Dik；2、增益率C4.5信息增益比：特征A对训练集D的信息增益比定义为其信息增益与训练集D关于特征值A的值的熵之比，即 其中3、基尼指数CART基尼指数度量数据区或训练集D的不纯度，定义为： 基尼指数数值越小，数据集的纯度越高；特征A的基尼指数定义为： S4中采取的自定义图像增强算法具体实现步骤如下：S401、为消除伪影和噪声对被调查样本组织的影响，采用白方和暗方两种方式来进行数据归一化处理和高光谱图像的校准； 上式中Rθ是高光谱图像的反射率，ImΦ是获取的高光谱图像，IdΦ是暗立方体图像，IwΦ是白立方体图像；利用高光谱相机采集的图像计算每个波长的光衰减吸收，如下式所示： 其中是光衰减吸收，是被摄物体的反射光强度，Iw是来自参考白板的反射光强度；利用修正的比尔-兰伯特定律测定了生物组织的光衰减吸收，下式所示： 其中x,y是空间像素坐标，ω是波长，是受试者完整的手臂部位的高光谱图像，bx,y,ω是含氧动脉和去氧静脉的高光谱图像，是参考的高光谱图像；S402、接下来采用移动平均滤波进一步滤除噪声；移动平均滤波基于统计规律，将连续的采样数据看成一个长度固定为N的队列，在新的一次测量后，上述队列的首数据去掉，其余N-1个数据依次前移，并将新的采样数据插入，作为新队列的尾，然后对这个队列进行算术运算，并将其结果作为本次测量的结果；假设输入为x，输出为y，则移动平均滤波器的计算公式如下所示： S5将S4中增强后的图像采用灰度直方图均衡化，k-means聚类算法进行图像分割，突出含氧动脉和去氧静脉的位置，最后根据二者比率计算SaO2大小，具体步骤如下：S501、灰度直方图均衡化，在图像识别的过程中，增加灰度对比度可以突出图像重要的特征，直方图均衡化是通过改变每个灰度级上像素点分布，使其都具有相同的象素点数，目的是使图像在整个灰度值动态变化范围内分布均匀化，改善图像的亮度分布状态，增强图像的视觉效果；首先考虑连续函数并且让变量r表示待增强图像的灰度级；假设r被归一化到[0,1]，且r＝0表示黑色，r＝1表示白色；对于连续函数，假设其变换函数为：s＝Tr0≤r≤1在原始图像中，对于每一个r对应着一个灰度值s；其中变换函数要满足以下条件：1、要求Tr在[0,1]中为单值，且单调递增；2、保证0≤r≤1时，0≤Tr≤1；使得输出的灰度级与输入的灰度级限定在同样的范围内；将公式s＝Tr改写：r＝T-1s0≤s≤1令Prr和Pss分别表示随机变量r和s的概率密度函数；由基本概率理论得到一个基本结果： 因此，变换变量s的概率密度函数由输入图像的灰度级概率密度函数和所选择的变换函数所决定；选择变换函数： 进而有： 将上式带入中得到Pss＝1；S502、采用K-means聚类算法对图像进行分割；利用高光谱相机采集到的每一幅图像计算光衰减；最后，在特征谱段的图像上叠加含氧动脉和去氧静脉的轮廓；K-means聚类算法的输入为包含N个数据对象的数据集D＝{x1,x2,…,xn}，输出为k个彼此独立的类簇C＝{C1,C2,…,Cn}，具体步骤如下：step1、从输入数据集D中随机选取出k个数据对象，作为初始的聚类簇中心点；step2、计算数据集D中数据对象与k个簇类中心点之间的相似度，将聚类对象分配到与其相似程度最高的类簇中心点代表的类簇中；step3、统计每个类簇中数据对象信息，取其均值作为新的类簇中心点，更新类簇中心点信息；step4、迭代执行上述的step2和step3直至算法执行完毕，类簇中心点不再发生变化；其中，对于包含m个属性的两个数据对象x＝{x1,x2,…,xm}和y＝{y1,y2,…,ym}，相似度的计算采用Pearson相关系数如下所示： S503、完成图像分割后，根据含氧动脉和去氧静脉的比率计算SaO2大小；根据光吸收和约化散射系数可以计算得出吸收系数： 其中，是光吸收，μaω是吸收系数，μsω是约化散射系数；吸收系数μaω和动脉血氧饱和度SaO2之间的相关性如下式所示： 其中，μaω是吸收系数，是含氧动脉中的血红蛋白成分，εHB是脱氧静脉中的血红蛋白成分。

全文数据：

权利要求：

百度查询： 吉林大学;林赐云 一种基于动静脉反射光谱和高光谱图像测量SaO2的COVID-19诊断方法

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

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