CN114298115B - 一种传感器交互运动意图获取方法及系统 - Google Patents

一种传感器交互运动意图获取方法及系统 Download PDF

Info

Publication number
CN114298115B
CN114298115B CN202210213808.2A CN202210213808A CN114298115B CN 114298115 B CN114298115 B CN 114298115B CN 202210213808 A CN202210213808 A CN 202210213808A CN 114298115 B CN114298115 B CN 114298115B
Authority
CN
China
Prior art keywords
sensor
data
output data
movement intention
acquiring
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210213808.2A
Other languages
English (en)
Other versions
CN114298115A (zh
Inventor
张雪波
郭宪
牛慧敏
安博
吕国学
方勇纯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nankai University
Original Assignee
Nankai University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nankai University filed Critical Nankai University
Priority to CN202210213808.2A priority Critical patent/CN114298115B/zh
Publication of CN114298115A publication Critical patent/CN114298115A/zh
Application granted granted Critical
Publication of CN114298115B publication Critical patent/CN114298115B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • User Interface Of Digital Computer (AREA)

Abstract

本发明涉及一种传感器交互运动意图获取方法及系统,包括:采集传感器数据进行初始化处理得到传感器初始数据;利用所述传感器初始数据获取传感器输出数据;根据所述传感器输出数据得到传感器交互运动意图,避免了数据传输过程中产生的噪声干扰问题,在实际应用中人体运动意图需要高精度高准确率识别提供了基础,同时对三维空间内不同方向的受力数据进行综合处理,与运动意图分析相结合,针对各种实际情况均具有可操作性。

Description

一种传感器交互运动意图获取方法及系统
技术领域
本发明涉及交通运输领域,具体涉及一种传感器交互运动意图获取方法及系统。
背景技术
普通传感器存在零点漂移问题,并且漂移量是变化的,人会在(对于系统来说)毫无征兆的突然启动,采样率设置不合理会导致这部分信息缺失,并且人在走路时本身会产生波动,这种波动可以理解为一种低频信息,这部分的信息会被其他频段的信号混杂干扰,同时电路中的感性元件和容性元件会产生一定的高频噪声,人在前进中会产生许多并非与行进方向或行进意图重合的扭动,这些扭动会干扰到对人真是意图的判断,因此急需一种合理有效的方案去分析获取运动意图。
发明内容
针对现有技术的不足,本发明提供了一种传感器交互运动意图获取方法,其特征在于,包括:
采集传感器数据进行初始化处理得到传感器初始数据;
利用所述传感器初始数据获取传感器输出数据;
根据所述传感器输出数据得到传感器交互运动意图;
对所述传感器初始数据基于滤波函数进行滤波处理得到传感器滤波数据;
利用所述传感器滤波数据得到传感器输出数据;
Figure 785412DEST_PATH_IMAGE001
其中,
Figure 109077DEST_PATH_IMAGE002
为传感器初始数据,
Figure 154393DEST_PATH_IMAGE003
为传感器滤波数据,f为传感器初始数据采样频率,
Figure 815182DEST_PATH_IMAGE004
c为传感器初始数据采样截止频率,
Figure 36079DEST_PATH_IMAGE005
为传感器初始数据相邻上一周期的传感器初始数据;
采集传感器滤波数据f(t),则此次采集传感器滤波数据的相邻下一时刻的传感器滤波数据为f(t+1),采集传感器滤波数据后y时刻的传感器滤波数据为f(t+y);
当f(t+1)与f(t)的差值delta大于波动标准阈值fth时,所述传感器滤波数据的相邻下一时刻的传感器输出数据为(1+n%)*f(t);
第一次采集传感器滤波数据后的第y时刻,当f(t+1)至f(t+y)的连续y个数据与f(t)的差值delta全部大于波动标准阈值fth,将f(t+y)作为传感器t+y时刻的输出数据;
对于t时刻与t+y之间的任意采样点t+x时刻(x<y),所述t+y时刻对应的f(t+x)减去f(x)大于delta,则其输出为(1+n%)^x * f(t);如果其值f(t+x)减去f(x)小于delta,则输出f(t+x),并将计数器清零;
其中,n、fth、delta、y为调整参数,x为触发波动阈值报警后至消除阈值报警或阈值报警计数计满间的某个采样点;
当传感器数据在无意识波动阈值内时,所述传感器输出数据对应的运动意图为无动作;
当传感器输出数据的纵向输出数据超出无意识波动阈值时,所述传感器输出数据对应的运动意图为前后运动;
当传感器输出数据的垂直输出数据超出无意识波动阈值时,所述传感器输出数据对应的运动意图为存在动作,传感器输出数据大于0时,交互运动意图相对速度加快,传感器输出数据小于0时,交互运动意图相对速度减慢;
传感器输出数据的水平输出数据超出无意识波动阈值时,所述传感器输出数据对应的运动意图为存在动作,传感器输出数据的左侧输出数据与传感器输出数据的水平右侧输出数据的差值大于0时,得到交互运动意图为向左运动;
传感器输出数据的水平输出数据超出无意识波动阈值时,所述传感器输出数据对应的运动意图为存在动作,传感器输出数据的左侧输出数据与传感器输出数据的水平右侧输出数据的差值小于0时,得到交互运动意图为向右运动;
其中,无意识波动阈值为发生晃动时且不存在运动意图时的设定阈值,传感器输出数据的左侧输出数据与传感器输出数据的水平右侧输出数据的差值的绝对值与交互运动意图左右运动的半径成反比。
优选的,所述采集传感器数据进行初始化处理得到传感器初始数据包括:
采集传感器数据后,利用所述传感器数据计算传感器数据补偿值后,利用所述传感器数据补偿值进行偏移补偿得到传感器补偿数据;
将所述传感器补偿数据作为传感器初始数据。
进一步的,所述利用传感器数据计算传感器数据补偿值进行偏移补偿的计算式如下:
Figure 124121DEST_PATH_IMAGE006
其中,m为传感器数据补偿值,n1、n2
Figure 656733DEST_PATH_IMAGE007
nx为传感器数据,x为传感器数据的数量。
优选的,所述传感器输出数据包括:
传感器输出数据的方向包括纵向输出数据、垂直输出数据与水平输出数据。
基于同一发明构思,本发明还提供了一种传感器交互运动意图获取系统,其特征在于,包括:
处理模块,用于采集传感器数据进行初始化处理得到传感器初始数据;
获取模块,用于利用所述传感器初始数据获取传感器输出数据;
输出模块,用于根据所述传感器输出数据得到传感器交互运动意图。
与最接近的现有技术相比,本发明具有的有益效果:
通过对数据初期的噪声处理,将干扰因素排除在外,后续通过阈值的设定,与初期滤波形成双重过滤,对于所得到的运动意图以及传感数据的准确性均有较大幅度提升,在实际应用中复杂情况的应对有良好的帮助,同时在运动意图分析中,引入多角度对运动意图进行定义,可以充分体现人在运动时所表现出的状态。
附图说明
图1是本发明提供的一种传感器交互运动意图获取方法流程图;
图2是本发明提供的一种传感器交互运动意图获取系统流程图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
实施例1:
本发明提供了一种传感器交互运动意图获取方法,如图1所示,包括:
采集传感器数据进行初始化处理得到传感器初始数据;
利用所述传感器初始数据获取传感器输出数据;
根据所述传感器输出数据得到传感器交互运动意图。
所述采集传感器数据进行初始化处理得到传感器初始数据包括:
采集传感器数据后,利用所述传感器数据计算传感器数据补偿值后,利用所述传感器数据补偿值进行偏移补偿得到传感器补偿数据;
将所述传感器补偿数据作为传感器初始数据。
所述利用传感器数据计算传感器数据补偿值进行偏移补偿的计算式如下:
Figure 121212DEST_PATH_IMAGE006
其中,m为传感器数据补偿值,n1、n2
Figure 695151DEST_PATH_IMAGE007
nx为传感器数据,x为传感器数据的数量。
所述利用传感器初始数据获取传感器输出数据包括:
对所述传感器初始数据基于滤波函数进行滤波处理得到传感器滤波数据;
利用所述传感器滤波数据得到传感器输出数据。
所述滤波函数包括:
Figure 954094DEST_PATH_IMAGE008
其中,
Figure 974003DEST_PATH_IMAGE009
为传感器初始数据,
Figure 117539DEST_PATH_IMAGE010
为传感器滤波数据,f为传感器初始数据采样频率,
Figure 906504DEST_PATH_IMAGE011
为传感器初始数据采样截止频率,
Figure 601927DEST_PATH_IMAGE012
为传感器初始数据相邻上一周期的传感器初始数据。
所述利用传感器滤波数据得到传感器输出数据包括:
采集传感器滤波数据f(t),则此次采集传感器滤波数据的相邻下一时刻的传感器滤波数据为f(t+1),采集传感器滤波数据后y时刻的传感器滤波数据为f(t+y);
当f(t+1)与f(t)的差值delta大于波动标准阈值fth时,所述传感器滤波数据的相邻下一时刻的传感器输出数据为(1+n%)*f(t);
第一次采集传感器滤波数据后的第y时刻,当f(t+1)至f(t+y)的连续y个数据与f(t)的差值delta全部大于波动标准阈值fth,将f(t+y)作为传感器t+y时刻的输出数据;
对于t时刻与t+y之间的任意采样点t+x时刻(x<y),所述t+y时刻对应的f(t+x)减去f(x)大于delta,则其输出为(1+n%)^x * f(t);如果其值f(t+x)减去f(x)小于delta,则输出f(t+x),并将计数器清零;
其中,n、fth、delta、y为调整参数,x为触发波动阈值报警后至消除阈值报警或阈值报警计数计满间的某个采样点。
所述根据传感器输出数据得到传感器交互运动意图包括:
当传感器数据在无意识波动阈值内时,所述传感器输出数据对应的运动意图为无动作;
当传感器输出数据的纵向输出数据超出无意识波动阈值时,所述传感器输出数据对应的运动意图为前后运动;
当传感器输出数据的垂直输出数据超出无意识波动阈值时,所述传感器输出数据对应的运动意图为存在动作,传感器输出数据大于0时,交互运动意图相对速度加快,传感器输出数据小于0时,交互运动意图相对速度减慢;
传感器输出数据的水平输出数据超出无意识波动阈值时,所述传感器输出数据对应的运动意图为存在动作,传感器输出数据的左侧输出数据与传感器输出数据的水平右侧输出数据的差值大于0时,得到交互运动意图为向左运动;
传感器输出数据的水平输出数据超出无意识波动阈值时,所述传感器输出数据对应的运动意图为存在动作,传感器输出数据的左侧输出数据与传感器输出数据的水平右侧输出数据的差值小于0时,得到交互运动意图为向右运动;
其中,无意识波动阈值为发生晃动时且不存在运动意图时的设定阈值,传感器输出数据的左侧输出数据与传感器输出数据的水平右侧输出数据的差值的绝对值与交互运动意图左右运动的半径成反比。
所述传感器输出数据包括:
传感器输出数据的方向包括纵向输出数据、垂直输出数据与水平输出数据。
实施例2:
本发明提供了一种传感器交互运动意图获取系统,如图2所示,包括:
处理模块,用于采集传感器数据进行初始化处理得到传感器初始数据;
获取模块,用于利用所述传感器初始数据获取传感器输出数据;
输出模块,用于根据所述传感器输出数据得到传感器交互运动意图。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。

Claims (5)

1.一种传感器交互运动意图获取方法,其特征在于,包括:
采集传感器数据进行初始化处理得到传感器初始数据;
利用所述传感器初始数据获取传感器输出数据;
对所述传感器初始数据基于滤波函数进行滤波处理得到传感器滤波数据;
利用所述传感器滤波数据得到传感器输出数据;
Figure 291256DEST_PATH_IMAGE001
其中,
Figure 797324DEST_PATH_IMAGE002
为传感器初始数据,
Figure 341438DEST_PATH_IMAGE003
为传感器滤波数据,f为传感器初始数据采样频率,
Figure 223943DEST_PATH_IMAGE004
c为传感器初始数据采样截止频率,
Figure 310848DEST_PATH_IMAGE005
为传感器初始数据相邻上一周期的传感器初始数据;
采集传感器滤波数据f(t),则此次采集传感器滤波数据的相邻下一时刻的传感器滤波数据为f(t+1),采集传感器滤波数据后y时刻的传感器滤波数据为f(t+y);
当f(t+1)与f(t)的差值delta大于波动标准阈值fth时,所述传感器滤波数据的相邻下一时刻的传感器输出数据为(1+n%)*f(t);
第一次采集传感器滤波数据后的第y时刻,当f(t+1)至f(t+y)的连续y个数据与f(t)的差值delta全部大于波动标准阈值fth,将f(t+y)作为传感器t+y时刻的输出数据;
对于t时刻与t+y之间的任意采样点t+x时刻(x<y),所述t+y时刻对应的f(t+x)减去f(x)大于delta,则其输出为(1+n%)^x * f(t);如果其值f(t+x)减去f(x)小于delta,则输出f(t+x),并将计数器清零;
其中,n、fth、delta、y为调整参数,x为触发波动阈值报警后至消除阈值报警或阈值报警计数计满间的某个采样点;
当传感器数据在无意识波动阈值内时,所述传感器输出数据对应的运动意图为无动作;
当传感器输出数据的纵向输出数据超出无意识波动阈值时,所述传感器输出数据对应的运动意图为前后运动;
当传感器输出数据的垂直输出数据超出无意识波动阈值时,所述传感器输出数据对应的运动意图为存在动作,传感器输出数据大于0时,交互运动意图相对速度加快,传感器输出数据小于0时,交互运动意图相对速度减慢;
传感器输出数据的水平输出数据超出无意识波动阈值时,所述传感器输出数据对应的运动意图为存在动作,传感器输出数据的左侧输出数据与传感器输出数据的水平右侧输出数据的差值大于0时,得到交互运动意图为向左运动;
传感器输出数据的水平输出数据超出无意识波动阈值时,所述传感器输出数据对应的运动意图为存在动作,传感器输出数据的左侧输出数据与传感器输出数据的水平右侧输出数据的差值小于0时,得到交互运动意图为向右运动;
其中,无意识波动阈值为发生晃动时且不存在运动意图时的设定阈值,传感器输出数据的左侧输出数据与传感器输出数据的水平右侧输出数据的差值的绝对值与交互运动意图左右运动的半径成反比。
2.如权利要求1所述的方法,其特征在于,所述采集传感器数据进行初始化处理得到传感器初始数据包括:
采集传感器数据后,利用所述传感器数据计算传感器数据补偿值后,利用所述传感器数据补偿值进行偏移补偿得到传感器补偿数据;
将所述传感器补偿数据作为传感器初始数据。
3.如权利要求2所述的方法,其特征在于,所述利用传感器数据计算传感器数据补偿值进行偏移补偿的计算式如下:
Figure 987817DEST_PATH_IMAGE006
其中,m为传感器数据补偿值,n1、n2
Figure 160172DEST_PATH_IMAGE007
nx为传感器数据,x为传感器数据的数量。
4.如权利要求1所述方法,其特征在于,所述传感器输出数据包括:
传感器输出数据的方向包括纵向输出数据、垂直输出数据与水平输出数据。
5.一种传感器交互运动意图获取系统,为实现上述权利要求1-4任一项所述的传感器交互运动意图获取方法,其特征在于,包括:
处理模块,用于采集传感器数据进行初始化处理得到传感器初始数据;
获取模块,用于利用所述传感器初始数据获取传感器输出数据;
输出模块,用于根据所述传感器输出数据得到传感器交互运动意图。
CN202210213808.2A 2022-03-07 2022-03-07 一种传感器交互运动意图获取方法及系统 Active CN114298115B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210213808.2A CN114298115B (zh) 2022-03-07 2022-03-07 一种传感器交互运动意图获取方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210213808.2A CN114298115B (zh) 2022-03-07 2022-03-07 一种传感器交互运动意图获取方法及系统

Publications (2)

Publication Number Publication Date
CN114298115A CN114298115A (zh) 2022-04-08
CN114298115B true CN114298115B (zh) 2022-05-24

Family

ID=80978474

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210213808.2A Active CN114298115B (zh) 2022-03-07 2022-03-07 一种传感器交互运动意图获取方法及系统

Country Status (1)

Country Link
CN (1) CN114298115B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101872259A (zh) * 2010-06-03 2010-10-27 中国人民解放军第二炮兵工程学院 自动录入书写内容的自然交互笔及笔迹检测方法
CN102247151A (zh) * 2011-04-25 2011-11-23 中国科学院合肥物质科学研究院 一种肌肉张力传感器及肌肉张力检测方法
CN104523403A (zh) * 2014-11-05 2015-04-22 陶宇虹 一种判断外骨骼助行机器人穿戴者下肢行动意图的方法
CN109324324A (zh) * 2017-08-01 2019-02-12 广州极飞科技有限公司 数据处理的方法及装置
CN112770286A (zh) * 2021-01-13 2021-05-07 深圳市木千科技有限公司 一种传感器的数据处理方法、装置和计算机设备

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5591945A (en) * 1995-04-19 1997-01-07 Elo Touchsystems, Inc. Acoustic touch position sensor using higher order horizontally polarized shear wave propagation
CN101301250A (zh) * 2008-07-08 2008-11-12 哈尔滨工业大学 五自由度外骨骼式上肢康复机器人交互康复训练控制策略
CN101872260B (zh) * 2010-06-03 2013-07-31 张通达 一种远程交互笔及书写笔迹检测方法
CN103876756A (zh) * 2014-04-18 2014-06-25 南京工程学院 下肢助力外骨骼机器人步态模式识别方法及系统
CN103984416B (zh) * 2014-06-10 2017-02-08 北京邮电大学 一种基于加速度传感器的手势识别方法
US10151604B2 (en) * 2014-08-13 2018-12-11 The Timken Company High speed sensing system
ES2649265T3 (es) * 2014-12-18 2018-01-11 Iveco Magirus Ag Método para controlar un aparato elevador, y aparato elevador con controlador que implementa dicho método
CN107757795A (zh) * 2017-11-21 2018-03-06 南阳师范学院 一种基于myRIO平台的自平衡小车控制系统及方法
CN110812122B (zh) * 2019-12-04 2022-03-29 上海金矢机器人科技有限公司 一种下肢康复机器人的坐站训练方法及系统
CN113946225B (zh) * 2021-12-20 2022-04-26 深圳市心流科技有限公司 一种手势锁定方法、智能仿生手、终端以及存储介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101872259A (zh) * 2010-06-03 2010-10-27 中国人民解放军第二炮兵工程学院 自动录入书写内容的自然交互笔及笔迹检测方法
CN102247151A (zh) * 2011-04-25 2011-11-23 中国科学院合肥物质科学研究院 一种肌肉张力传感器及肌肉张力检测方法
CN104523403A (zh) * 2014-11-05 2015-04-22 陶宇虹 一种判断外骨骼助行机器人穿戴者下肢行动意图的方法
CN109324324A (zh) * 2017-08-01 2019-02-12 广州极飞科技有限公司 数据处理的方法及装置
CN112770286A (zh) * 2021-01-13 2021-05-07 深圳市木千科技有限公司 一种传感器的数据处理方法、装置和计算机设备

Also Published As

Publication number Publication date
CN114298115A (zh) 2022-04-08

Similar Documents

Publication Publication Date Title
CN103674062B (zh) 基于Allan方差与ARMA模型分析提高陀螺仪测量精度的方法
CN103678299B (zh) 一种监控视频摘要的方法及装置
JP2011070623A5 (zh)
US9292123B2 (en) System and method for sample rate adaption
CN105105737A (zh) 基于光电容积描记和谱分析的运动状态心率监测方法
CN103027667A (zh) 脉搏波的特征参数提取
CN104546007B (zh) 胎动检测抗干扰处理方法及装置
EP1973101A1 (en) Pitch extraction with inhibition of harmonics and sub-harmonics of the fundamental frequency
CN107704192B (zh) 一种滤波算法及其探测器
CN105654428A (zh) 图像降噪方法及系统
CN112057080A (zh) 一种基于分阶段特征提取的冻结步态检测方法和系统
CN114298115B (zh) 一种传感器交互运动意图获取方法及系统
CN109034100A (zh) 人脸图样检测方法、装置、设备及存储介质
EP1686561A1 (en) Determination of a common fundamental frequency of harmonic signals
CN111510109B (zh) 一种信号滤波的方法、装置、设备和介质
CN111679091A (zh) 高实时性分段自适应转速测量方法、装置、介质及设备
CN114757855B (zh) 动作数据修正方法、装置、设备及存储介质
CN108492308B (zh) 一种基于相互结构引导滤波的变分光流的确定方法及系统
CN103106192A (zh) 文学作品作者识别方法及装置
CN106197466B (zh) 一种数据处理方法、装置及终端
CN108022429B (zh) 一种车辆检测的方法及装置
CN115048954A (zh) 一种仿视网膜的目标检测方法、装置、存储介质及终端
CN112494008A (zh) 基于ppg信号的呼吸率测量方法及装置
CN114795168B (zh) 生命体征参数心率计算方法及系统
CN112116917A (zh) 基于相位跃变度的电抗器本体与风机声信号分离方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant