CN102662238A - 一种具有在轨自诊断和补偿功能的空间光学相机 - Google Patents
一种具有在轨自诊断和补偿功能的空间光学相机 Download PDFInfo
- Publication number
- CN102662238A CN102662238A CN2012101349684A CN201210134968A CN102662238A CN 102662238 A CN102662238 A CN 102662238A CN 2012101349684 A CN2012101349684 A CN 2012101349684A CN 201210134968 A CN201210134968 A CN 201210134968A CN 102662238 A CN102662238 A CN 102662238A
- Authority
- CN
- China
- Prior art keywords
- wavefront
- focus
- value
- module
- camera
- 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.)
- Granted
Links
Images
Landscapes
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
一种具有在轨自诊断和补偿功能的空间光学相机属于空间光学获取领域,该相机包括:空间相机光学系统与相位差异波前探测。空间相机光学系统是卡塞格林、同轴三反和离轴三反等光学系统形式。相位差异波前探测单元包括:在焦多级拼接探测器、离焦多级拼接探测器、在焦探测器接口、离焦探测器接口、数据采集模块、数据预处理模块、计算模块、存储模块和控制接口;本发明无需额外增加高精度波前测量装置,能够降低高性能空间光学系统对设计、结构、热控等方面的要求,实现空间光学系统的高精度在轨诊断和补偿成像。
Description
技术领域
本发明属于空间光学获取领域,特别涉及一种具有在轨自诊断和补偿功能的空间光学相机。
背景技术
空间光学系统由于卫星发射时的冲击、过载、振动、热环境的变化、微重力环境及大气湍流影响使波面产生畸变,导致光学成像系统成像质量下降。尤其是对未来更高性能的空间相机,随着系统口径及焦距的增加,相应的对轻量化要求也越来越高,各反射镜在轨位置误差也将进一步增大。如果继续采用传统的无诊断和补偿功能的被动光学成像,会使设计分辨率和使用成像分辨率之间差距进一步增大。
现有的空间相机光学系统主要是卡塞格林、同轴三反和离轴三反等光学系统形式。卡塞格林主要包括两块反射镜,反射镜中大的称为主镜,小的称为副镜,主镜和副镜同轴;双通道卡塞格林通常在主镜中央开孔,当地面景物光信息进入主镜后,反射入射进入副镜,经由副镜反射通过主镜中央的开孔,成像于主镜后面。同轴三反光学系统主要包括主镜、次镜和三镜,主镜、次镜和三镜同在一条光轴上,地面景物光信息进入主镜后,依次反射入射进入次镜和三镜后,在光轴上成像。离轴三反光学系统和同轴三反光学系统类似,都是主要包括主镜、次镜和三镜,不同的是主镜、次镜和三镜不在同一条光轴上,地面景物光信息进入主镜后,依次反射入射进入次镜和三镜后,在三镜的光轴上成像。
发明内容
为了解决现有技术中存在的问题,本发明提供了一种具有在轨自诊断和补偿功能的空间光学相机,该光学相机较好的解决了空间相机因受自身及空间环境影响导致成像质量下降,有效的诊断和补偿了光学系统波前误差,提高了系统对地成像的分辨率。
本发明解决技术问题所采用的技术方案如下:
一种具有在轨自诊断和补偿功能的空间光学相机包括空间相机光学系统与相位差异波前探测单元;该相位差异波前探测单元包括:在焦多级拼接探测器、离焦多级拼接探测器、在焦探测器接口、离焦探测器接口、数据采集模块、数据预处理模块、计算模块、存储模块和控制接口;
首先,地面景物经过空间相机光学系统成像在在焦多级拼接探测器和离焦多级拼接探测器上,在数据采集模块的控制下,在焦探测器接口和离焦探测器接口分别从在焦多级拼接探测器和离焦多级拼接探测器同时采集信息,由数据采集模块将两组采集的信息整合成一帧数据发送到数据预处理模块,数据预处理模块将预处理的子图像发送给计算模块,计算模块把得到的目标恢复图像和波前畸变信息发送给存储模块和控制接口,控制接口通过波前畸变信息对光学相机的离焦像差进行补偿。
本发明的有益效果是:本发明在单通道多级线阵拼接探测成像方法的基础上,增加已知离焦量的相位差异成像通道,双通道的使用改善了反演问题的病态性,能对波前信息进行精确的解算,并根据波前信息对图像进行复原、光学系统的在轨自检测和对离焦误差的补偿。并结合空间光学载荷受体积、重量和系统复杂性限制的特点,光路中仅是在焦面位置引入了一个阶梯状的板,在阶梯的两级上分别各有一条线阵拼接探测器,并且阶梯的高度差为已知量,这样便可以对同一地面景物同时采集焦面和离焦面的双通道图像。本发明无需额外增加高精度波前测量装置,能够降低高性能空间光学相机对设计、结构、热控等方面的要求,实现空间光学相机的高精度在轨诊断和补偿成像。
附图说明
图1本发明一种具有在轨自诊断和补偿功能的空间光学相机光学系统的具体实施方式的光路示意图。
图2本发明一种具有在轨自诊断和补偿功能的空间光学相机相位差异波前探测单元的在焦多级拼接探测器和离焦多级拼接探测器。
图3本发明一种具有在轨自诊断和补偿功能的空间光学相机相位差异波前探测单元的工作流程图。
图4本发明一种具有在轨自诊断和补偿功能的空间光学相机相位差异波前探测单元根据已知测得的波前对采集来的图像进行复原处理部分示意图。
图中:1、入瞳,2、离轴主镜,3、离轴次镜,4、视场光阑,5、离轴三镜,6、里奥光阑,7、调焦镜,8、在焦多级拼接探测器,9、离焦多级拼接探测器,10、在焦探测器接口,11、离焦探测器接口,12、数据采集模块,13、数据预处理模块,14、计算模块,15、存储模块,16、控制接口。
具体实施方式
下面结合附图和实施例对本发明做进一步详细说明。
如图1所示,一种具有在轨自诊断和补偿功能的空间光学相机,该相机包括:离轴三反空间光学系统,离轴三反空间光学系统包括:入瞳1、离轴主镜2、离轴次镜3、视场光阑4、离轴三镜5、里奥光阑6和调焦镜7;相机还包括相位差异波前探测单元,该相位差异波前探测单元包括:在焦多级拼接探测器8、离焦多级拼接探测器9、在焦探测器接口10、离焦探测器接口11、数据采集模块12、数据预处理模块13、计算模块14、存储模块15、控制接口16和图像复原模块17;
首先,地面景物光信息经过入瞳1入射到离轴主镜2,经由离轴主镜2反射后出射光入射进入离轴次镜3,经由离轴次镜3反射后出射光通过视场光阑4入射进入离轴三镜5,经由离轴三镜5反射后出射光通过里奥光阑6入射进入调焦镜7,经由调焦镜7反射后出射光分别进入在焦多级拼接探测器8和离焦多级拼接探测器9;如图2所示,双探测器拼接固定在同一个基准面上,双通道分光通过在多级线阵探测器前放置调焦镜7来实现,并且调焦镜7能够引入固定光程位相差,双区域拼接探测器固定在一起,提高双通道成像受失调误差影响的一致性。
如图3所示,在数据采集模块12的控制下,在焦探测器接口10和离焦探测器接口11分别从在焦多级拼接探测器8和离焦多级拼接探测器9同时采集信息,曝光时间与卫星相对地面的运动相关,也与单级线阵探测器的积分时间相关。在焦多级拼接探测器8和离焦多级拼接探测器9采集信息的离焦量为地面景物对应波前PV值的1至1.5λ,其中,λ为对应波长,本系统的样机采用的是1λ。数据采集模块12将两组采集的图像信息整合成一帧数据发送到数据预处理模块13,数据预处理模块13主要工作是平场,在两幅图像上通过图像配准分别抠出细节丰富载波能力强的区域以减少波前解算的工作量,但是解算到波前后对图像的复原是对完整的整幅图像进行的,最后把预处理后的子图像发送给基于相位差异法计算模块14,它以公式(1)为目标函数,以各帧图像所对应的波前为未知变量进行最优化搜索,最终将得到目标的恢复图像以及波前畸变信息。
其中,u表示频域坐标,C表示使用的通道数,如本发明是双通道探测器,那么C=2;N表示单幅图像的像素总数;φc(υ)表示第c通道的波前相位的在光瞳坐标υ处得值,θc表示第c通道已知的固定离焦相位;{α}表示需要求解的泽尼克(zernike)系数,αm∈{α}表示zernike系数向量中的第m项系数;Zm表示第m项泽尼克多项式基底;为通道c的噪声读出方差的倒数;F为傅里叶变换,F-1表示傅里叶逆变换,υ表示光瞳面坐标,P表示光瞳函数;Sc表示通道c的光学传递函数;Dc=F(dc),式中dc表示第c通道采集到的图像;γ表示非负的正则项系数;当搜索算法退出的时候,Fimg为目标理想图像的期望频谱,对其求傅里叶逆变换就得到了目标的恢复图像。把复原后的图像及波前畸变信息发送给复原后图像及测得的波前存储模块15,并根据波前畸变信息中的离焦信息发送给根据波前信息对光学系统离焦进行补偿的控制接口16,对离焦像差进行补偿。
结合在轨光学系统波前随时间变化并不剧烈的特点,在每次使用相位差异法解算出波前之后可以在接下来的一段时间之内用该波前恢复图像,而不需要从新计算波前。如图4所示,在成像质量没有下降并且波前在前一段时间已经测量的前提下,通过数据采集模块12的控制下,在焦探测器接口10和离焦探测器接口11分别从在焦多级拼接探测器8和离焦多级拼接探测器9同时采集两幅图像,并由数据采集模块12将两组采集的图像信息整合成一帧数据发送到数据预处理模块13,数据预处理模块13的主要工作是平场,最后把预处理后的子图像发送给图像复原模块17,图像复原模块17从复原后图像及测得的波前存储模块15中读取波前信息,应用公式(4)直接对图像进行复原,然后把复原后的图像传入15。
一种具有在轨自诊断和补偿功能的空间光学相机,其图像处理的过程如下:
利用在焦多级拼接探测器8和离焦多级拼接探测器9采集焦面图像和已知离焦量的离焦面的图像,并分别通过在焦探测器接口10和离焦探测器接口11由数据采集模块12获得图像数据,并把同时采集的图像数据对合并成一个数据帧发送给数据预处理模块13;
数据预处理模块13对从数据采集模块12发送来的图像数据进行配准及截取共同的子图像、噪声估计及去噪等预处理,并把处理后的结果发送给计算模块14;
计算模块14主要包括数学最优化搜索引擎和代价函数计算两个部分。其中数学最优化搜索引擎使用的是LBFGS(The limited memory variation of theBroyden-Fletcher-Goldfarb-Shanno(BFGS))算法,包含搜索引擎初始化和搜索引擎的主体两部分。求解步骤如下:
步骤一:由数据预处理模块13传来的数据经过计算模块14中的搜索引擎初始化部分转化为频域数据,并把作为未知量的波前值设为0;
步骤二:进入计算模块14中的代价函数计算部分,这部分首先根据波前值和由数据预处理模块13传来的数据根据公式4来估计目标理想图像的期望,并用公式1对这个期望值进行评价得到代价函数值,并求出代价函数对波前αm的偏导数;
步骤三:用求得的代价函数值和偏导数值通过LBFGS算法来调整作为未知量的波前值;
步骤四:根据目标理想图像的期望的代价函数值和目标理想图像的期望的代价函数对公式3中波前αm的偏导数,判断是否满足退出条件,则计算模块14把得到目标恢复图像和波前畸变信息发送给存储模块15,否则重复步骤二。对于本系统,退出条件为,算法迭代次数大于1000次,或目标函数值在10次迭代中的下降总和小于10-15,或各项偏导值小于10-15,满足以上三个条件之一便退出。
把求得的波前值及对目标理想图像期望值的估计发送到存储模块15,并进行显示。
Claims (7)
1.一种具有在轨自诊断和补偿功能的空间光学相机包括空间相机光学系统,其特征在于,该相机还包括相位差异波前探测单元;该相位差异波前探测单元包括:在焦多级拼接探测器(8)、离焦多级拼接探测器(9)、在焦探测器接口(10)、离焦探测器接口(11)、数据采集模块(12)、数据预处理模块(13)、计算模块(14)、存储模块(15)和控制接口(16);
地面景物经过空间相机光学系统成像在在焦多级拼接探测器(8)和离焦多级拼接探测器(9)上,在数据采集模块(12)的控制下,在焦探测器接口(10)和离焦探测器接口(11)分别从在焦多级拼接探测器(8)和离焦多级拼接探测器(9)同时采集地面景物信息,由数据采集模块(12)将两组采集的信息整合成一帧数据发送到数据预处理模块(13),数据预处理模块(13)将预处理的子图像发送给计算模块(14),计算模块(14)把得到的目标恢复图像和波前畸变信息发送给存储模块(15)和控制接口(16),控制接口(16)通过波前畸变信息对光学相机的离焦像差进行补偿。
2.如权利要求1所述的一种具有在轨自诊断和补偿功能的空间光学相机,其特征在于,所述在焦多级拼接探测器(8)和离焦多级拼接探测器(9)固定在一个阶梯面上,阶梯的高度差为定值。
3.如权利要求1所述的一种具有在轨自诊断和补偿功能的空间光学相机,其特征在于,所述在焦多级拼接探测器(8)和离焦多级拼接探测器(9)采集信息的离焦量为地面景物对应波前PV值的1~1.5λ,其中,λ为对应波长。
4.如权利要求1所述的一种具有在轨自诊断和补偿功能的空间光学相机,其特征在于,所述计算模块(14)包括数学最优化搜索引擎和代价函数计算部分。
5.如权利要求4所述的一种具有在轨自诊断和补偿功能的空间光学相机,其特征在于,所述数学最优化搜索引擎包含搜索引擎初始化和搜索引擎的主体两部分。
6.如权利要求1或4或5所述的一种具有在轨自诊断和补偿功能的空间光学相机,其特征在于,所述计算模块(14)工作方法如下:
步骤一:计算模块(14)接收数据预处理模块(13)传送的数据后,一部分数据经过数学最优化搜索引擎初始化转化为频域数据,把未知量的波前值设为0;
步骤二:计算模块(14)接收数据预处理模块(13)传送的另一部分数据进入代价函数计算部分,代价函数计算部分首先根据波前值设定值和进入代价函数计算部分的数据,再结合公式来估计目标理想图像的期望,并用公式 对这个期望进行评价得到代价函数值,并求出代价函数对波前各项泽尼克多项式的值的偏导数,其中C表示使用的通道数;N表示单幅图像的像素总数;为通道c的噪声读出方差的倒数;Dc=F(dc),式中dc表示第c通道采集到的图像;γ表示非负的正则项系数;{α}表示需要求解的泽尼克(zernike)系数;
步骤三:用步骤二中求得的代价函数值和偏导数值通过LBFGS(Thelimited memory variation of the Broyden-Fletcher-Goldfarb-Shanno(BFGS))算法来调整作为未知量的波前值;
步骤四:根据目标理想图像期望的代价函数值和目标理想图像期望的代价函数对波前各项泽尼克多项式值的偏导数,判断是否满足退出条件,如果满足,计算模块(14)把得到目标恢复图像和波前畸变信息发送给存储模块(15),否则重复步骤二。
7.如权利要求6所述的一种具有在轨自诊断和补偿功能的空间光学相机,其特征在于,所述步骤四中的退出条件为:算法迭代次数大于1000次,或目标函数值在10次迭代中的下降总和小于10-15,或各项偏导值小于10-15中的任意一项。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210134968.4A CN102662238B (zh) | 2012-05-03 | 2012-05-03 | 一种具有在轨自诊断和补偿功能的空间光学相机 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210134968.4A CN102662238B (zh) | 2012-05-03 | 2012-05-03 | 一种具有在轨自诊断和补偿功能的空间光学相机 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102662238A true CN102662238A (zh) | 2012-09-12 |
CN102662238B CN102662238B (zh) | 2014-01-15 |
Family
ID=46771758
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210134968.4A Active CN102662238B (zh) | 2012-05-03 | 2012-05-03 | 一种具有在轨自诊断和补偿功能的空间光学相机 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102662238B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102917243A (zh) * | 2012-10-22 | 2013-02-06 | 中国科学院西安光学精密机械研究所 | 空间相机彩色图像质量评价检测系统 |
WO2016102731A1 (es) * | 2014-12-22 | 2016-06-30 | Universidad De La Laguna | Método para determinar la amplitud compleja del campo electromagnético asociado a una escena |
CN106768876A (zh) * | 2016-11-29 | 2017-05-31 | 中国科学院长春光学精密机械与物理研究所 | 基于星点光斑的空间天文望远镜波前探测方法 |
CN107505798A (zh) * | 2017-06-16 | 2017-12-22 | 中国科学院长春光学精密机械与物理研究所 | 基于线性拟合方式的航空遥感器实时匀速斜距调焦方法和装置 |
CN107655507A (zh) * | 2017-09-27 | 2018-02-02 | 北京空间机电研究所 | 一种焦平面弧形拼接方法 |
CN110345921A (zh) * | 2019-06-12 | 2019-10-18 | 中国农业大学 | 立体视场视觉测量及垂轴像差和轴向像差校正方法及系统 |
CN114217434A (zh) * | 2021-11-22 | 2022-03-22 | 北京空间机电研究所 | 一种高分辨率大视场成像方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101266334A (zh) * | 2008-05-07 | 2008-09-17 | 北京理工大学 | 基于随机并行梯度下降算法的无波前探测器自适应光学系统 |
CN101738721A (zh) * | 2009-12-25 | 2010-06-16 | 中国科学院光电技术研究所 | 一种哈特曼波前传感器无系统误差的自适应光学成像系统 |
-
2012
- 2012-05-03 CN CN201210134968.4A patent/CN102662238B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101266334A (zh) * | 2008-05-07 | 2008-09-17 | 北京理工大学 | 基于随机并行梯度下降算法的无波前探测器自适应光学系统 |
CN101738721A (zh) * | 2009-12-25 | 2010-06-16 | 中国科学院光电技术研究所 | 一种哈特曼波前传感器无系统误差的自适应光学成像系统 |
Non-Patent Citations (1)
Title |
---|
陈荣利等: "高分辨率TDICCD相机轻量化技术", 《航天返回与遥感》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102917243A (zh) * | 2012-10-22 | 2013-02-06 | 中国科学院西安光学精密机械研究所 | 空间相机彩色图像质量评价检测系统 |
CN107209061B (zh) * | 2014-12-22 | 2020-04-14 | 拉古纳大学 | 用于确定场景相关电磁场的复振幅的方法 |
US10230940B2 (en) | 2014-12-22 | 2019-03-12 | Universidad De La Laguna | Method for determining the complex amplitude of the electromagnetic field associated with a scene |
KR20170099959A (ko) * | 2014-12-22 | 2017-09-01 | 유니버시다드 데 라 라구나 | 장면과 관련된 전자기장의 복소 진폭을 결정하는 방법 |
CN107209061A (zh) * | 2014-12-22 | 2017-09-26 | 拉古纳大学 | 用于确定场景相关电磁场的复振幅的方法 |
KR102501402B1 (ko) | 2014-12-22 | 2023-02-20 | 유니버시다드 데 라 라구나 | 장면과 관련된 전자기장의 복소 진폭을 결정하는 방법 |
WO2016102731A1 (es) * | 2014-12-22 | 2016-06-30 | Universidad De La Laguna | Método para determinar la amplitud compleja del campo electromagnético asociado a una escena |
TWI687661B (zh) * | 2014-12-22 | 2020-03-11 | 拉古納大學 | 用於確定場景相關電磁場的複振幅的方法及裝置 |
CN106768876A (zh) * | 2016-11-29 | 2017-05-31 | 中国科学院长春光学精密机械与物理研究所 | 基于星点光斑的空间天文望远镜波前探测方法 |
CN107505798A (zh) * | 2017-06-16 | 2017-12-22 | 中国科学院长春光学精密机械与物理研究所 | 基于线性拟合方式的航空遥感器实时匀速斜距调焦方法和装置 |
CN107505798B (zh) * | 2017-06-16 | 2019-09-03 | 中国科学院长春光学精密机械与物理研究所 | 基于线性拟合方式的航空遥感器实时匀速斜距调焦方法和装置 |
CN107655507B (zh) * | 2017-09-27 | 2019-10-22 | 北京空间机电研究所 | 一种焦平面弧形拼接方法 |
CN107655507A (zh) * | 2017-09-27 | 2018-02-02 | 北京空间机电研究所 | 一种焦平面弧形拼接方法 |
CN110345921A (zh) * | 2019-06-12 | 2019-10-18 | 中国农业大学 | 立体视场视觉测量及垂轴像差和轴向像差校正方法及系统 |
CN114217434A (zh) * | 2021-11-22 | 2022-03-22 | 北京空间机电研究所 | 一种高分辨率大视场成像方法 |
CN114217434B (zh) * | 2021-11-22 | 2024-03-29 | 北京空间机电研究所 | 一种高分辨率大视场成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102662238B (zh) | 2014-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102662238B (zh) | 一种具有在轨自诊断和补偿功能的空间光学相机 | |
US11778289B2 (en) | Multi-camera imaging systems | |
Le Bouquin et al. | PIONIER: a 4-telescope visitor instrument at VLTI | |
CN106803892B (zh) | 一种基于光场测量的光场高清晰成像方法 | |
CN105589210B (zh) | 一种基于光瞳调制的数字化合成孔径成像方法 | |
Kendrick et al. | Flat-panel space-based space surveillance sensor | |
US7274442B2 (en) | Closed-loop wavefront sensor using field programmable gate array | |
US9305378B1 (en) | Lenslet, beamwalk and tilt diversity for anisoplanatic imaging by large-aperture telescopes | |
CN100568045C (zh) | 一种适合于白天工作的高分辨力成像自适应光学望远镜 | |
CN105954734B (zh) | 大口径激光雷达光轴监测装置 | |
US20110157600A1 (en) | Optical wave-front recovery for active and adaptive imaging control | |
US5350911A (en) | Wavefront error estimation derived from observation of arbitrary unknown extended scenes | |
CN110146180A (zh) | 一种基于焦面哈特曼波前传感器的大视场图像清晰化装置及方法 | |
CN104019898A (zh) | 一种超灵敏光谱成像天文望远镜及天文光谱成像方法 | |
CN110794576A (zh) | 一种基于相位调制的光学合成孔径成像望远镜阵列偏心误差探测方法 | |
US20150192411A1 (en) | Two axis interferometer tracking device and method | |
Buscher | Getting the most out of COAST | |
Surdej et al. | The 4m international liquid mirror telescope (ILMT) | |
Nain et al. | Improving band to band registration accuracy of SEVIRI level 1.5 products | |
Marchetti et al. | ERIS adaptive optics system design | |
US12080000B2 (en) | Initial orbit determination using angular velocity and angular acceleration measurements | |
Van et al. | Optical design with non-rotationally symmetric field mapping for infrared search and track system | |
Drexler et al. | Digital adaptive optics with interferometric homodyne encoding for mitigating atmospheric turbulence for imaging applications | |
Wilson et al. | Profiling atmospheric turbulence using a non-cooperative target and a camera bank and validating with sonic anemometers | |
Filippov et al. | Long-Range Thermal 3D Perception in Low Contrast Environments |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |