基于认知任务测试的近红外脑成像波谱分类方法
【专利摘要】本发明公开了一种基于认知任务测试的近红外脑成像图谱分类方法，该近红外脑成像图谱中至少包括对应于氧合血红蛋白浓度的第一曲线和对应于脱氧血红蛋白浓度的第二曲线，该分类方法具体步骤包括：第一步骤，读取近红外脑成像仪获得的近红外脑成像图谱，提取所述近红外脑成像图谱的第一曲线及第二曲线的波形特征；第二步骤，将第一步骤获得的所述波形特征与预设值进行比较并根据比较结果对所述近红外脑成像图谱进行分类。本发明所提供的基于认知任务测试的近红外脑成像图谱分类方法，可应用于心理学领域的基础研究或临床应用上对近红外脑成像图谱进行自动分类。
【专利说明】
基于认知任务测试的近红外脑成像波谱分类方法
技术领域
[0001]本发明涉及一种图像处理方法，特别是涉及基于认知任务测试的近红外脑成像图谱分类方法。该方法可应用于心理学领域的基础研究或临床应用上对近红外脑成像图谱进行自动分类。
【背景技术】
[0002]目前，心理障碍的心理症状的识别、分析与鉴别过程主要遵循“症状-综合征-病史”模式，主要是依靠症状识别的经验性诊断，在临床上还缺乏便捷的辅助症状识别的技术。随着例如功能近红外光谱脑成像(fNIRS)、功能性磁共振(fMRI)、正电子发射断层显像术(PET)和单电子发射断层显像术(SPECT)等无创检测方法的逐渐应用，深化了研究工作者对心理障碍患者的发病机理和疾病过程的认识。比如，现已经证实抑郁症患者的前额功能影像有异常。但由于fMR1、PET、SPECT等设备使用时或多或少会有诸如体积偏大、价格偏高、噪声大、时间长、不耐受、幽闭(可能影响测试结果)与不自然(可能影响测试结果)等方面原因，因此在目前的心理障碍诊断过程中并没有得到实际的临床应用。
[0003]近红外光谱脑成像(fNIRS)技术，其基本原理是由于近红外光700-1400nm的波长能够透过烦骨2-3cm，其中780nm被脱氧血红蛋白(Deoxy-Hb)吸收和散射，而830nm被氧合血红蛋白(Oxy-Hb))吸收和散射，这能间接反映脑皮层神经元的活动和脑认知功能活动。相比较于以上所述的其他几种脑成像技术，fNIRS有着安静无噪声、使用成本低、容易操作和维护以及对被试的限制少等诸多优点，因此正逐渐的成为该领域重要成像手段。但功能近红外光谱脑成像所得图谱一般难以辨识，工作繁琐，使用者难以有效理解使用。具体来说，以心理障碍临床研究中的基于认知任务测试为例，在现有的传统的图谱分类分析方法中，多是比较平均后的任务期与静息期的oxy-Hb浓度差或使用地形图(Topographic map)，来间接反映区域皮层的激活情况(Noda et al，2012;Liu et al ,2014)，虽然得到一些发现，但是平均后掩盖了一些信息，未能反映oxy-Hb、deoxy_Hb随时间的变化(生态效应)。再者，由于NIRS具有时间分辨率高的特点，如果按照传统的做法沿用fMRI的数据解读方式，则不能充分发挥其实时测量被试者完成心理任务时脑皮层氧合和脱氧血红蛋白变化的特点，以及血红蛋白在任务状态下的转换关系。
[0004]因此，亟需一种图谱解读及分类方法，使得在例如认知任务测试等心理任务中，被测试者的近红外脑成像图谱可以被合理的分类，进而方便相关工作者对心理障碍的进一步深入研究以及促进之后的临床应用。
【发明内容】
[0005]鉴于现有技术存在的上述问题，本发明的目的在于提供一种基于认知任务测试的近红外脑成像图谱分类方法，该方法可以实现在预设的多项认知任务的刺激呈现下，对连续获取的受检者脑部对应的近红外脑成像图谱，进行后期处理并按照波形进行分类，方便进行下一步研究工作。
[0006]为了实现上述目的，本发明提供的一种基于认知任务测试的近红外脑成像图谱分类方法，该近红外脑成像图谱中至少包括对应于氧合血红蛋白浓度的第一曲线和对应于脱氧血红蛋白浓度的第二曲线，该分类方法具体步骤包括:
[0007]第一步骤，读取近红外脑成像仪获得的近红外脑成像图谱，提取所述近红外脑成像图谱的第一曲线及第二曲线的波形特征；
[0008]第二步骤，将第一步骤获得的所述波形特征与预设值进行比较并根据比较结果对所述近红外脑成像图谱进行分类。
[0009]在本发明中，近红外脑成像图谱的波形特征的提取主要是利用已有的数学方法对波形进行基于数字信号的特征分析，例如傅里叶变换、希尔伯特变换、Gabor变换以及小波变换等。
[0010]作为优选，所述波形特征包括位相、波幅、频率或趋势。
[0011]作为优选，所述认知任务包括在特定时间周期内执行的情绪图片任务、语言流畅任务、对指运动任务和症状图片任务(看图描述)。在本发明中，情绪图片任务是作为试验环境的适应性任务。语言流畅任务(VFT)为行业内普遍用于检测词语检索与语言执行功能的测试手段;对指运动任务为检查前额叶和运动区的运动调节功能;症状图片刺激任务为检测前额叶情绪调节功能。通过认知任务刺激，可以检测心理障碍患者前额叶的激活模式。
[0012]作为优选，所述认知任务的时间周期模块包括持续15s的任务前静息期、持续30s的任务期和持续15s的任务后静息期。一般4-5个刺激模块组合成一组等效认知任务的刺激范式。这种模块组合范式保证了两任务期之间有30秒的静息期。语言流畅任务与对指运动任务刺激模块组合各为4次，总任务时长各4分钟。症状图片任务(看图描述)刺激模块组合为5次，总任务时长5分钟。
[0013]作为优选，在所述第二步骤中，具体依照如下步骤进行分类:
[0014]I)判断位相步骤:根据所述第一曲线与所述第二曲线的相对位置关系对所述近红外脑成像图谱进行位相标记；
[0015]2)判断波形步骤:根据所述第一曲线与所述第二曲线在任务期内的波幅、频率和/或趋势对所述近红外脑成像图谱进行波形标记；
[0016]3)依据所述位相标记和波形标记，统计所述近红外脑成像图谱中的正相的双曲线波形、高波幅波形及正相的开口波形，以及负相的低平波形、交叉分离波形及负相的开口状波形、正相高波幅波形以及正相开口波形，进而根据统计结果计算波形构成比并根据所述波形构成比对所述近红外脑成像图谱进行分类。
[0017]具体地，可以依据波形构成比最大的波形作为分类参考。例如，正相双曲线波形的波形构成比最高，则可以将该图谱分类为正相双曲线波形图谱。其余以此类推，不难判断，依照本发明提出的分类方法，基于认知任务测试所产生的近红外脑成像图谱可以通过本发明技术方案被至少分为6类，即正相双曲线波形图谱、负相低平波形图谱、负相交叉分离波形图谱、负相开口状波形图谱、正相高波幅波形图谱和正相开口状波形图谱。
[0018]本发明还提供了另一种针对近红外脑成像图谱分类的实施方式，在这一方法中，不同于前一实施方式，在第二步骤中，将第一步骤获得的所述波形特征与对照图谱库中的对照图谱直接进行比较，进而根据比较结果对所述近红外脑成像图谱进行分类。
[0019]作为优选，所述对照图谱包括正相双曲线波形图谱、负相低平波形图谱、负相交叉分离波形图谱、负相开口状波形图谱、正相高波幅波形图谱和正相开口状波形图谱。
[0020]本发明所提供的基于认知任务测试的近红外脑成像图谱分类方法，可应用于医学领域的基础研究或临床应用上对近红外脑成像图谱进行自动分类。在认知任务的刺激呈现下，获取高参考价值的近红外脑成像图谱经过分类之后，便于理解使用。例如，在一组测试中，通过统计经过分类的各种近红外脑成像图谱所占的比例，可能揭示心理障碍在影像上的具体表现，具有较高的研究价值和商业前景。
【附图说明】
[0021]图1为本发明的基于认知任务测试的近红外脑成像图谱分类方法的流程图。
[0022]图2为本发明的基于认知任务测试的近红外脑成像图谱分类方法的另一实施方式的流程图。
[0023]图3为本发明的基于认知任务测试的近红外脑成像图谱分类方法的再一实施方式的流程图。
[0024]图4为本发明的基于认知任务测试的近红外脑成像图谱分类方法归类的典型的正相双曲线波形图谱。
[0025]图5为本发明的基于认知任务测试的近红外脑成像图谱分类方法归类的负相低平波形图谱。
[0026]图6为本发明的基于认知任务测试的近红外脑成像图谱分类方法的归类的负相交叉分离波形图谱。
[0027]图7为本发明的基于认知任务测试的近红外脑成像图谱分类方法的归类的负相开口状波形图谱。
[0028]图8为本发明的基于认知任务测试的近红外脑成像图谱分类方法的归类的正相高波幅波形图谱。
[0029]图9为本发明的基于认知任务测试的近红外脑成像图谱分类方法的归类的正相开口状波形图谱。
【具体实施方式】
[0030]以下参照附图对本发明所提出的技术方案做进一步详细的说明。
[0031]如图1所示，本发明提供的一种基于认知任务测试的近红外脑成像图谱分类方法，该近红外脑成像图谱中至少包括对应于氧合血红蛋白浓度的第一曲线和对应于脱氧血红蛋白浓度的第二曲线，通常来说，在图谱中第一曲线通常以红色曲线表示，而第二曲线则通常以蓝色表示;本发明提供的该分类方法具体步骤包括:
[0032]第一步骤，读取近红外脑成像仪(本文中未示出)获得的近红外脑成像图谱，提取所述近红外脑成像图谱的第一曲线及第二曲线的波形特征;此处的波形特征通常可以包括所述近红外脑成像图谱中的第一曲线及第二曲线的位相、波幅、频率或趋势等。对于图谱中的第一曲线和第二曲线的波性特征的识别可以采用已有的图谱识别算法，例如傅里叶变换、希尔伯特变换、Gabor变换以及小波变换等。而针对两幅图谱之间的直接比较，也可以采用二值化比较或者直方图形式进行相似性比较。
[0033]在完成上述第一步骤后，在第二步骤中，可以将第一步骤获得的所述波形特征与预设值进行比较并根据比较结果对所述近红外脑成像图谱进行分类。在这一步骤中，预设值的设定是根据人工选定的标准对照图谱的波形特征进行确定。例如，如果以图3所示的图谱作为对照图谱话，可以先人为确定该作为对照图谱的位相、波幅、频率或趋势并各自作为预设值进行存储。同样可以想见，针对如图4至图9的图谱，也需要与之对应的多组位相、波幅、频率或趋势的预设值。通过这种方式，可以将任意图谱较准确的进行图谱分类，便于后续相关工作者做出针对性解读。
[0034]图2和图3给出了本发明的另外两种【具体实施方式】。
[0035]在本发明中，情绪图片任务是作为试验环境的适应性任务。语言流畅任务(VFT)为行业内普遍用于检测词语检索与语言执行功能的测试手段;对指运动任务为检查前额叶与运动区的运动调节功能;症状图片刺激任务为检测前额叶情绪调节功能。可以检测心理障碍患者前额叶的激活模式。所述认知任务包括在特定时间周期内执行的情绪图片任务、语言流畅任务、对指运动任务和症状图片刺激任务。所述认知任务的时间周期包括持续15s的任务前静息期、持续30s的任务期和持续15s的任务后静息期。也包括4-5个刺激模块组合成一组等效认知任务的刺激范式。这种模块组合范式保证了两任务期之间有30秒的静息期。语言流畅任务与对指运动任务刺激模块组合各为4次，总任务时长各4分钟。症状图片任务(看图描述)刺激模块组合为5次，总任务时长5分钟。
[0036]进一步地，在所述第二步骤中，具体依照如下步骤进行分类:
[0037]I)判断位相步骤:根据所述第一曲线与所述第二曲线的相对位置关系对所述近红外脑成像图谱进行位相标记;具体地，当所述第一曲线的位置高于所述第二曲线时，标记为正相;反之则标记为负相；
[0038]2)判断波形步骤:根据所述第一曲线与所述第二曲线在任务期内的波幅、频率和/或趋势对所述近红外脑成像图谱进行波形标记；
[0039]具体地，对应于认知任务测试的任务期，出现明显的波峰，且波幅适度(在1/3以内)，标记为双曲线波形，反之则标记为低平波形;对第一曲线和第二曲线整体趋势开口状的，标记为开口波形，反之标记为交叉分离波形；
[0040]3)依据所述位相标记和波形标记，统计所述近红外脑成像图谱中的正相的双曲线波形、高波幅波形及正相的开口波形，以及负相的低平波形、交叉分离波形及负相的开口状波形、正相高波幅波形以及正相开口波形，进而根据统计结果计算波形构成比并根据所述波形构成比对所述近红外脑成像图谱进行分类。
[0041]具体地，可以依据波形构成比最大的波形作为分类参考。例如，正相双曲线波形的波形构成比最高，则可以将该图谱分类为正相双曲线波形图谱。其余以此类推，不难判断，依照本发明提出的分类方法，基于认知任务测试所产生的近红外脑成像图谱可以通过本发明技术方案被至少分为6类，即正相双曲线波形图谱、负相低平波形图谱、负相交叉分离波形图谱、负相开口状波形图谱、正相高波幅波形图谱和正相开口状波形图谱。
[0042]本发明还提供了另一种针对近红外脑成像图谱分类的实施方式，如图2所示，在这一方法中，不同于前一实施方式，在第二步骤中，将第一步骤获得的所述波形特征与对照图谱库中的对照图谱直接进行比较，进而根据比较结果对所述近红外脑成像图谱进行分类。
[0043]以下提供本发明的基于认知任务测试的近红外脑成像仪具体应用的实验例。
[0044]1、受试对象:
[0045]该研究主要在清华大学玉泉医院进行，受试者的招募工作从2012年5月持续至2015年3月。玉泉医院伦理委员会已经批准通过该项研究，所有受试者均签署了知情同意书。
[0046]入组标准:我们总共招募了45名健康对照志愿者(广告招募58名，有13名不符合入组标准被排除)，按整群抽样的原则选取了 74名患有抑郁障碍(单相抑郁，初步筛检93例，有19例不符合入组标准被排除)及60名患焦虑障碍(初步筛检78例，有18例不符合入组标准被排除)的门诊病人。所招募病人均由2名经验丰富的精神科医生(副高及以上职称)依据DSM-1V。所有受试者的母语均为汉语(普通话)，以保证汉语版本的VFT能够顺利进行。VFT总数，正常对照组与抑郁障碍组及焦虑障碍组差异均有显著性(丨=0.048，？ = 0.001)，但抑郁障碍组及焦虑障碍组之间差异无显著性(P = 0.172)。
[0047]在进行NI RS检查的当天，我们将使用24项汉密尔顿抑郁量表(Hami I tonDepress1n Rating Scale for Depress1n,HAMD)来对受试者的抑郁障碍状进行评定;使用汉密顿焦虑量表(Hamilton Anxiety Scale，HAMA)对其焦虑障碍状进行评定;使用耶鲁强迫量表(Yale-Brown obsessive compulsive scale,Y-B0CS)对强迫症状进行评定。研究期间，病例组中部分患者服用一种或两种小剂量药物(抗抑郁药物SSRI)，但用药时间少于2周，有24例未服药的病人。
[0048]为了控制混杂因素，我们使用单因素方差分析ANOVA对年龄、受教育程度和性别进行了比较，三者差异均无显著性(P>0.05)。结果见表I。
[0049]排除标准:1、既往有双相障碍躁狂发作史或轻躁狂发作史;2、既往为伴精神病性症状的重度抑郁发作患者;3、患有其他心理障碍或神经疾病或躯体疾病(如甲低等)；4、服用抗抑郁或精神病药超过2周以上;5、排除引起认知损伤的脑外伤以及酒精或物质依赖或成瘾;6、受教育程度低于6年，母语为非汉语者。
[0050]2、方法
[0051]2.1检测工具
[0052]近红外光谱成像仪作为近红外脑成像模块采集近红外脑成像图谱，选用日本岛津的LabNIRS_52channel，米用CW(continuous wave)技术模式，测量氧合血红蛋白(Oxy-Hb)浓度、总氧合血红蛋白(t一Hb)浓度和脱氧合血红蛋白(Deoxy—Hb)相对浓度。
[0053]2.2刺激呈现
[0054]言语流畅任务(verbalfluency task，VFTr):任务由E.Prime 软件(2.0)进行呈现。
[0055]静息状态:安装近红外光极后，告知受检者，保持安静，尽量不要思考，持续5分钟。氧合与脱氧血红蛋白相对浓度定标为±0.1 (mmol/L X cm)。观察氧合与脱氧血红蛋白相对浓度双曲线处于较稳定的状态。
[0056]为了保证语言流畅任务的稳定性，在该任务前加用情绪图片任务作为预试验，起环境适应的作用。情绪图片任务:任务负荷轻度，它由4s的图片呈现期和2s的静息期组成，共15幅图(正性5幅+中性5幅+负性5幅)，检测时间为90s。该任务未纳入统计。
[0057]激活任务为语言流畅任务(VFT):依据实验心理任务规范和临床预实验，每个任务block为:15s的静息期、30s的任务期和15s的任务后静息期组成。四组blocks连续完成，时间为240s。进行任务时，受试者被告知需在任务期尽可能地多地说出所提示的汉语词汇(蔬菜、家用电器、四条腿的动物和水果)，所选词汇类别均为中性，四组词汇均为等效。四组b locks能够有效地反应血氧浓度的时间依赖曲线变化。
[0058]3、NIRS 检测
[0059]光极安装参考脑电图10/20系统制作了前额脑区光极帽(4X7光极安放位点)，利用3D坐标定位系统依次记录眉心、颅顶正中、左耳前、右耳前4个参考位置及各个光极，面积14cmX20cm，包括Broca与语言调节相关区45、44、46区，非语言调节相关区:8、9、10等。共获取45标定逻辑通道图谱。
[0060]4、波形分析工具
[0061]软件名称:多项认知任务近红外脑成像图谱分析软件(简称:近红外脑图谱分析软件)V3.0.著作权人:刘破资。权利取得方式:原始取得。权利范围:全部权利。开发完成日期2013年12月16日。首次发表日期:2014年I月3日。中华人民共和国国家版权局，计算机软件著作权登记证书号:软著登字第0959027号。软件登记号:2015SR071941。
[0062]5、NIRS数据分类分析
[0063]我们以4次等效语言流畅任务呈现波形的位相、波幅、频率、趋势为依据对图谱进行分类处理。
[0064]波形分类:基于语言流畅任务下典型病例的4个Blocks任务总时间240s的前额标定脑区的45个通道OXY-HB和deoxy-Hb浓度变化的双曲线连续波谱为参考，进行症状波形分类:以45个通道(logic channel)为原点，形成的45个方框的波形为计量单位，纵坐标为波幅，向上为正相，向下为负相，横坐标为时间(240秒)。静息状态:OXY-HB和deoxy-Hb浓度变化的双曲线重合。任务状态:认知任务状态时，OXY-HB和deoxy-Hb浓度变化的双曲线波动状态。氧合与脱氧血红蛋白相对浓度定标为±0.1 (mmol/L X cm)。
[0065]正相双曲线正常波形范围:双曲线连续图谱，波形红线(Oxy-Hb)在上，蓝线(Deoxy-Hb)在下，标定波幅适度(在1/3以内)，部分通道可见于4个任务相关的波峰，波形比较有规律，类似正弦波。如图3所示。
[0066]负相低平波形:蓝线(Deoxy-Hb)在上，红线(Oxy-Hb)在下，称为负相波形，波幅低平，少见与任务相关的波峰。如图4所示。
[0067]负相交叉分离波形:蓝线(Deoxy-Hb)在上，红线(Oxy-Hb)在下，波形负相、交叉、波动或分离，波形漂移不稳。如图5所示。
[0068]负相开口状波形:蓝线(Deoxy-Hb)在上，波形红线(Oxy-Hb)在下，有些波未见明显回落，有叠加现象，局部呈负相向右开口状，似符号〈。如图6所示。
[0069]正相高波幅波形:红线(Oxy-Hb)在上，蓝线(Deoxy-Hb)在下，标定波幅明显高于正常，可见于4个任务相关的高波峰，有些波幅大于1/3，波形比较有规律。如图7所示。
[0070]正相开口状波形:红线(Oxy-Hb)在上，蓝线(Deoxy-Hb)在下，有些波未见明显回落，有叠加现象，波形呈正相开口状〈。如图8所示。
[0071]自变量:分组，因变量:NIRS变化值
[0072]正常范围波形、负相低平波形、负相开口状波形、负相交叉分离波形、正相高波幅波形、正相开口波形构成比(0-45/45，0-100 % )为自变量，以正常对照组、依据DSM-1V诊断为抑郁障碍(单相抑郁)、焦虑障碍为因变量。
[0073]将上述进行了症状图谱分类的正常对照45例、抑郁障碍组74例、焦虑障碍60例，按照HAMD、HAMA量表总分为症状金标准(HAMD彡7分，HAMA彡7分，YALE彡1分为正常;抑郁障碍组需HAMD多17分;焦虑障碍组需HAMA多14分)，计算前额叶标定脑区NIRS vft症状图谱分类的波形构成比和ROC。
[0074]6、统计分析
[0075]所有的统计分析过程均采用SPSS 19.0软件(SPSS Inc，Bei jing，China)。一般情况及临床量表，通过使用单因素方差分析(ANOVA)来检测各组年龄、受教育情况及临床量表的组间差异。性别差异的检测使用卡方检验，差异显著性设定为P<0.05。正常组与抑郁障碍3个亚组及焦虑障碍2个亚组之间的量表差异使用t检验，P值设定为0.01。正常对照组与抑郁障碍组、焦虑障碍组6种波形构成比的差异使用t检验，P值设定为0.01。
【主权项】
1.基于认知任务测试的近红外脑成像图谱分类方法，该近红外脑成像图谱中至少包括对应于氧合血红蛋白浓度的第一曲线和对应于脱氧血红蛋白浓度的第二曲线，该分类方法具体步骤包括: 第一步骤，读取近红外脑成像仪获得的近红外脑成像图谱，提取所述近红外脑成像图谱的第一曲线及第二曲线的波形特征； 第二步骤，将第一步骤获得的所述波形特征与预设值进行比较并根据比较结果对所述近红外脑成像图谱进行分类。2.如权利要求1所述的基于认知任务测试的近红外脑成像图谱分类方法，所述波形特征包括位相、波幅、频率或趋势。3.如权利要求2所述的基于认知任务测试的近红外脑成像图谱分类方法，所述认知任务包括在特定时间周期内执行的情绪图片任务、语言流畅任务、对指运动任务和症状图片任务单元的组合。4.如权利要求3所述的基于认知任务测试的近红外脑成像图谱分类方法，所述认知任务单元包括情绪图片任务、语言流畅任务、对指运动任务或症状图片任务。5.如权利要求4所述的基于认知任务测试的近红外脑成像图谱分类方法，所述认知任务的刺激模块的时间周期包括持续15s的任务前静息期、持续30s的任务期和持续15s的任务后静息期。6.如权利要求1所述的基于认知任务测试的近红外脑成像图谱分类方法，在所述第二步骤中，具体依照如下步骤进行分类: 1)判断位相步骤:根据所述第一曲线与所述第二曲线的相对位置关系对所述近红外脑成像图谱进行位相标记； 2)判断波形步骤:根据所述第一曲线与所述第二曲线在任务期内的波幅、频率和/或趋势对所述近红外脑成像图谱进行波形标记； 3)依据所述位相标记和波形标记，统计所述近红外脑成像图谱中的正相的双曲线波形、高波幅波形及正相的开口波形，以及负相的低平波形、交叉分离波形及负相的开口状波形、正相高波幅波形以及正相开口波形，进而根据统计结果计算波形构成比并根据所述波形构成比对所述近红外脑成像图谱进行分类。7.如权利要求6所述的基于认知任务测试的近红外脑成像图谱分类方法，依据波形构成比最大的波形作为分类参考。8.如权利要求1所述的基于认知任务测试的近红外脑成像图谱分类方法，在第二步骤中，将第一步骤获得的所述波形特征与对照图谱库中的对照图谱直接进行比较，进而根据比较结果对所述近红外脑成像图谱进行分类。9.如权利要求8所述的基于认知任务测试的近红外脑成像图谱分类方法，所述对照图谱包括正相双曲线波形图谱、负相低平波形图谱、负相交叉分离波形图谱、负相开口状波形图谱、正相高波幅波形图谱和正相开口状波形图谱。
【文档编号】A61B5/00GK106037645SQ201610342318
【公开日】2016年10月26日
【申请日】2016年5月23日
【发明人】刘破资
【申请人】清华大学玉泉医院