CN115442534B - 空间目标天基成像自动多级曝光控制方法及系统 - Google Patents
空间目标天基成像自动多级曝光控制方法及系统 Download PDFInfo
- Publication number
- CN115442534B CN115442534B CN202210933688.3A CN202210933688A CN115442534B CN 115442534 B CN115442534 B CN 115442534B CN 202210933688 A CN202210933688 A CN 202210933688A CN 115442534 B CN115442534 B CN 115442534B
- Authority
- CN
- China
- Prior art keywords
- target
- star
- space
- reflectivity
- exposure time
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000002310 reflectometry Methods 0.000 claims abstract description 54
- 239000013598 vector Substances 0.000 claims abstract description 44
- 238000012544 monitoring process Methods 0.000 claims abstract description 27
- 239000000463 material Substances 0.000 claims abstract description 26
- 210000001747 pupil Anatomy 0.000 claims abstract description 20
- 230000004927 fusion Effects 0.000 claims abstract description 7
- 238000001228 spectrum Methods 0.000 claims description 18
- 239000013077 target material Substances 0.000 claims description 10
- 230000009466 transformation Effects 0.000 claims description 10
- 238000013507 mapping Methods 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000005286 illumination Methods 0.000 description 6
- 230000003044 adaptive effect Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000011248 coating agent Substances 0.000 description 2
- 238000000576 coating method Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Images
Landscapes
- Studio Devices (AREA)
Abstract
Description
技术领域
本发明涉及太空态势感知领域,具体地,涉及一种空间目标天基成像自动多级曝光控制方法及系统。
背景技术
曝光时间是相机成像的重要参数,曝光时间越长相机焦面获得的能量越多,反之则焦面获得的能量越少。一般来说,在光学成像过程中,目标亮度越大,为了避免相机焦面探测元过饱和,积分时间应当减小;当目标亮度较小时,则可适当增加光学相机的积分时间。
在空间目标成像天基成像过程中,目标的成像条件复杂:太阳光照的照射相位与时间相关、监视星与成像对象的相位位置也是变化的。当前空间目标天基成像主要采用的是固定曝光时间的方法,成像结果往往存在过曝光或者欠曝光的情况,成像效果不佳,而已有应用的自动曝光方法均是基于图像信息的方式,在轨应用对图像处理的要求高,不容易实现。因此,本发明提出了一种空间目标天基成像自动多级曝光控制方法,通过在轨成像光照及位置关系,计算最优曝光时间序列,可以有效地进行空间目标天基成像系统的曝光时间控制,实现空间目标的清晰成像。方法简单有效,在轨实现复杂度低。
针对相机的成像曝光时间控制问题,文献1(.杨海涛,《一种基于亮度直方图的自动曝光控制方法》,光学学报,2007)、文献2(王雷斌,《一种自适应快速曝光控制方法》,测绘科学,2009)、文献3(陈绕青《一种空间面阵CCD成像的曝光时间计算方法》)均是通过拍摄图像信息,获取目标直方图等质量评价参数后进行曝光时间计算;文献4(赵晓琳《一种自适应精确调节CCD曝光时间的方法》,传感器技术学报,2011年)为曝光时间的控制电路设计方法。
专利文献CN104301622A(申请号:201410522278.5)公开了一种天基空间目标成像曝光时间自适应控制方法及装置,方法包括1)计算根据预设的最长、中间以及最短曝光时间所获取的对应图像的信息熵;2)根据所有所述信息熵通过加权平均方式获取理论最优曝光时间,进而在预设曝光时间档中选取最接近的一档曝光时间作为当前最优曝光时间,并记录所述当前最优曝光时间的获取次数;3)判断所述获取次数是否小于等于预设阈值,若是,则执行步骤4),否则执行步骤6);4)计算根据当前最优曝光时间所获取的对应图像的信息熵;5)根据当前最优曝光时间更新最长、中间档以及最短曝光时间,以及更新对应的信息熵,并返回执行步骤2);6)控制成像相机以当前最优曝光时间作为最优曝光时间对空间目标进行持续成像。
专利文献CN104184957B(申请号:201410476250.2)公开了一种适用于空间探测成像的基于自适应预期图像平均亮度的自动曝光控制方法,首先对当前帧拍摄图像采用自适应阈值分割,识别出当前帧拍摄主体图像,由当前帧拍摄主体图像平均亮度得到当前帧拍摄图像平均亮度,然后将识别出的当前帧拍摄主体图像进行直方图均衡化处理,得到直方图均衡化后的预期图像平均亮度,再根据当前帧拍摄图像平均亮度与预期图像平均亮度的差值的正负和大小来选取曝光时间,采用连续帧不断循环迭代使当前帧拍摄图像平均亮度达到直方图均衡化后的自适应预期图像平均亮度。
目前网上可查询的文件均是以获取先验图像后根据图像参数进行曝光时间设置,与本发明提出的基于光照条件及空间观测相位的方法不同,本发明无需先验图像进行即可合理地设置在轨多级曝光时间。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种空间目标天基成像自动多级曝光控制方法及系统。
根据本发明提供的一种空间目标天基成像自动多级曝光控制方法,包括:
步骤S2:预估观测目标面的典型材料的反射率,并将占比面积最大材料的反射率参数作为反射率基准值ρ0;
步骤S4:根据目标在相机入瞳处的辐照度以及设定预存的定标结果设置最优曝光时间T0;
步骤S5:在最优曝光时间T0的基础上,根据不同表面反射率参数确定多级曝光时间,最终通过多幅不同曝光图像融合达到动态范围扩展的目的。
优选地,所述步骤S3采用:
其中,Ec表示目标在相机入瞳处的辐照度;Esun(λ)表示光谱为λ值的太阳辐照强度;R表示目标卫星与监视卫星的距离;λ1、λ2分别表示天基成像平台的传感器的探测光谱范围的上下限;ρ(λ)为光谱为λ值的反射率。
优选地,所述设定预存的定标结果是根据卫星发射前或在轨对已知参数进行标定,得到目标亮度与最佳曝光时间的映射表。
优选地,所述步骤S5采用:根据目标表面材料的反射率估计值设置多级曝光时间;
其中,ρi表示预估的不同目标材料的表面反射率;ρ0表示反射率基准值;T0表示最优曝光时间。
优选地,所述步骤S5采用:采用离散小波变换将每张图像分解为高频子图和低频子图,并将高频子图和低频子图分别进行融合。
根据本发明提供的一种空间目标天基成像自动多级曝光控制系统,包括:
模块M2:预估观测目标面的典型材料的反射率,并将占比面积最大材料的反射率参数作为反射率基准值ρ0;
模块M4:根据目标在相机入瞳处的辐照度以及设定预存的定标结果设置最优曝光时间T0;
模块M5:在最优曝光时间T0的基础上,根据不同表面反射率参数确定多级曝光时间,最终通过多幅不同曝光图像融合达到动态范围扩展的目的。
优选地,在所述模块M3中:
其中,Ec表示目标在相机入瞳处的辐照度;Esun(λ)表示光谱为λ值的太阳辐照强度;R表示目标卫星与监视卫星的距离;λ1、λ2分别表示天基成像平台的传感器的探测光谱范围的上下限;ρ(λ)为光谱为λ值的反射率。
优选地,所述设定预存的定标结果是根据卫星发射前或在轨对已知参数进行标定,得到目标亮度与最佳曝光时间的映射表。
优选地,在所述模块M5中:根据目标表面材料的反射率估计值设置多级曝光时间;
其中,ρi表示预估的不同目标材料的表面反射率;ρ0表示反射率基准值;T0表示最优曝光时间。
优选地,在所述模块M5中:采用离散小波变换将每张图像分解为高频子图和低频子图,并将高频子图和低频子图分别进行融合。
与现有技术相比,本发明具有如下的有益效果:
1、本发明通过计算最优曝光时间序列,可以有效地进行空间目标天基成像系统的曝光时间控制,实现空间目标的清晰成像;
2、相对传统基于图像直方图的自动曝光方法,本发明基于空间几何关系及光照条件,计算方案简单有效,在轨实现复杂度低。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为空间目标天基成像自动多级曝光示意图。
图2为空间目标成像示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
本发明提供了一种空间目标天基成像自动多级曝光控制方法及系统,通过在轨成像光照及位置关系,计算最优曝光时间序列,可以有效地进行空间目标天基成像系统的曝光时间控制,实现空间目标的清晰成像。本发明解决现有天基空间目标成像曝光时间设置不合适导致的目标过曝光或者欠曝光引起的成像质量下降问题。方法简单有效,在轨实现复杂度低。
实施例1
根据本发明提供的一种空间目标天基成像自动多级曝光控制方法,如图1至2所示,包括:
步骤S1:根据目标星轨道、监视星轨道以及成像时刻计算在惯性坐标系下的目标星指向太阳的单位矢量目标星指向监视星的单位矢量/>以及垂直于目标星的被观测面且沿着太阳方向的单位矢量/>其中,单位尺量/>的定义在惯性坐标系中,/>单位矢量来源被观测目标的本体坐标系下单位尺量/>可根据其与惯性坐标系的转换关系进行变换。
步骤S2:预估观测目标面的典型材料的反射率,并将占比面积最大材料的反射率参数作为反射率基准值ρ0;其中,目标的反射率针对合作目标可根据已知的目标表面参数,针对未知非合作目标其参数采用常见的卫星表面热控包覆材料表面发射参数,材料0:ρ0、材料1:ρ1、材料2:ρ2…
具体地,所述步骤S3采用:
其中,Ec表示目标在相机入瞳处的辐照度;Esun(λ)表示光谱为λ值的太阳辐照强度;R表示目标卫星与监视卫星的距离;λ1、λ2分别表示天基成像平台的传感器的探测光谱范围的上下限;ρ(λ)为光谱为λ值的反射率;一般可设置为常数,公式可简化为:
步骤S4:根据目标在相机入瞳处的辐照度以及设定预存的定标结果设置最优曝光时间T0;
具体地,所述设定预存的定标结果是根据卫星发射前或在轨对已知参数进行标定,得到目标亮度与最佳曝光时间的映射表;目标亮度与最佳曝光时间的映射表可在轨期间根据情况进行定期地调整。
步骤S5:在最优曝光时间T0的基础上,根据不同表面反射率参数确定多级曝光时间,最终通过多幅不同曝光图像融合达到动态范围扩展的目的。本发明能够准确设置对空间目标的曝光时间设置,实现目标的清晰成像。
具体地,所述步骤S5采用:根据目标表面材料的反射率估计值设置多级曝光时间;
其中,ρi表示预估的不同目标材料的表面反射率;ρ0表示反射率基准值;T0表示最优曝光时间。
更为具体地,所述步骤S5采用:采用离散小波变换将每张图像分解为高频子图和低频子图,并将高频子图和低频子图分别进行融合。即将多积分图像序列分别进行离散小波变换,对于每张图像将被分解为三个高频子图和一个低频子图,对序列中的各低频子带之间进行融合,各高频子带之间进行融合。
进一步地,针对不同反射率的曝光时间计算结果可进行必要的合并,单个目标某一时刻下的曝光分级数一般不大于5次。
根据本发明提供的一种空间目标天基成像自动多级曝光控制系统,包括:
模块M1:根据目标星轨道、监视星轨道以及成像时刻计算在惯性坐标系下的目标星指向太阳的单位矢量目标星指向监视星的单位矢量/>以及垂直于目标星的被观测面且沿着太阳方向的单位矢量/>其中,单位尺量/>的定义在惯性坐标系中,/>单位矢量来源被观测目标的本体坐标系下单位尺量/>可根据其与惯性坐标系的转换关系进行变换。
模块M2:预估观测目标面的典型材料的反射率,并将占比面积最大材料的反射率参数作为反射率基准值ρ0;其中,目标的反射率针对合作目标可根据已知的目标表面参数,针对未知非合作目标其参数采用常见的卫星表面热控包覆材料表面发射参数,材料0:ρ0、材料1:ρ1、材料2:ρ2…
具体地,在所述模块M3中:
其中,Ec表示目标在相机入瞳处的辐照度;Esun(λ)表示光谱为λ值的太阳辐照强度;R表示目标卫星与监视卫星的距离;λ1、λ2分别表示天基成像平台的传感器的探测光谱范围的上下限;ρ(λ)为光谱为λ值的反射率;一般可设置为常数,公式可简化为:
模块M4:根据目标在相机入瞳处的辐照度以及设定预存的定标结果设置最优曝光时间T0;
具体地,所述设定预存的定标结果是根据卫星发射前或在轨对已知参数进行标定,得到目标亮度与最佳曝光时间的映射表;目标亮度与最佳曝光时间的映射表可在轨期间根据情况进行定期地调整。
模块M5:在最优曝光时间T0的基础上,根据不同表面反射率参数确定多级曝光时间,最终通过多幅不同曝光图像融合达到动态范围扩展的目的。本发明能够准确设置对空间目标的曝光时间设置,实现目标的清晰成像。
具体地,在所述模块M5中:根据目标表面材料的反射率估计值设置多级曝光时间;
其中,ρi表示预估的不同目标材料的表面反射率;ρ0表示反射率基准值;T0表示最优曝光时间。
更为具体地,在所述模块M5中:采用离散小波变换将每张图像分解为高频子图和低频子图,并将高频子图和低频子图分别进行融合。即将多积分图像序列分别进行离散小波变换,对于每张图像将被分解为三个高频子图和一个低频子图,对序列中的各低频子带之间进行融合,各高频子带之间进行融合。
进一步地,针对不同反射率的曝光时间计算结果可进行必要的合并,单个目标某一时刻下的曝光分级数一般不大于5次。
本领域技术人员知道,除了以纯计算机可读程序代码方式实现本发明提供的系统、装置及其各个模块以外,完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统、装置及其各个模块以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同程序。所以,本发明提供的系统、装置及其各个模块可以被认为是一种硬件部件,而对其内包括的用于实现各种程序的模块也可以视为硬件部件内的结构;也可以将用于实现各种功能的模块视为既可以是实现方法的软件程序又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (10)
1.一种空间目标天基成像自动多级曝光控制方法,其特征在于,包括:
步骤S4:根据目标星在相机入瞳处的辐照度以及设定预存的定标结果设置最优曝光时间T0;
步骤S5:在最优曝光时间T0的基础上,根据不同目标材料的表面反射率估计值确定多级曝光时间,最终通过多幅不同曝光图像融合达到动态范围扩展的目的。
3.根据权利要求1所述的空间目标天基成像自动多级曝光控制方法,其特征在于,所述设定预存的定标结果是根据卫星发射前或在轨对已知参数进行标定,得到目标亮度与最佳曝光时间的映射表。
5.根据权利要求1所述的空间目标天基成像自动多级曝光控制方法,其特征在于,所述步骤S5采用:采用离散小波变换将每张图像分解为高频子图和低频子图,并将高频子图和低频子图分别进行融合。
6.一种空间目标天基成像自动多级曝光控制系统,其特征在于,包括:
模块M4:根据目标星在相机入瞳处的辐照度以及设定预存的定标结果设置最优曝光时间T0;
模块M5:在最优曝光时间T0的基础上,根据不同目标材料的表面反射率估计值确定多级曝光时间,最终通过多幅不同曝光图像融合达到动态范围扩展的目的。
8.根据权利要求6所述的空间目标天基成像自动多级曝光控制系统,其特征在于,所述设定预存的定标结果是根据卫星发射前或在轨对已知参数进行标定,得到目标亮度与最佳曝光时间的映射表。
10.根据权利要求6所述的空间目标天基成像自动多级曝光控制系统,其特征在于,在所述模块M5中:采用离散小波变换将每张图像分解为高频子图和低频子图,并将高频子图和低频子图分别进行融合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210933688.3A CN115442534B (zh) | 2022-08-04 | 2022-08-04 | 空间目标天基成像自动多级曝光控制方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210933688.3A CN115442534B (zh) | 2022-08-04 | 2022-08-04 | 空间目标天基成像自动多级曝光控制方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115442534A CN115442534A (zh) | 2022-12-06 |
CN115442534B true CN115442534B (zh) | 2023-07-14 |
Family
ID=84242882
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210933688.3A Active CN115442534B (zh) | 2022-08-04 | 2022-08-04 | 空间目标天基成像自动多级曝光控制方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115442534B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004340784A (ja) * | 2003-05-16 | 2004-12-02 | Mitsubishi Electric Corp | スターセンサ |
CN103675794A (zh) * | 2013-12-04 | 2014-03-26 | 北京空间机电研究所 | 基于时空统一特性的航天光学遥感器成像仿真方法 |
CN104301622A (zh) * | 2014-09-30 | 2015-01-21 | 上海微小卫星工程中心 | 一种天基空间目标成像曝光时间自适应控制方法及装置 |
CN106441237A (zh) * | 2015-08-10 | 2017-02-22 | 北京空间飞行器总体设计部 | 一种光学遥感卫星相机成像参数在轨自主调整方法 |
CN109581371A (zh) * | 2018-11-09 | 2019-04-05 | 中国科学院空间应用工程与技术中心 | 遥感相机成像参数的自动调整方法 |
CN110913150A (zh) * | 2019-11-18 | 2020-03-24 | 西北工业大学 | 基于空间平台的自适应曝光成像方法 |
-
2022
- 2022-08-04 CN CN202210933688.3A patent/CN115442534B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004340784A (ja) * | 2003-05-16 | 2004-12-02 | Mitsubishi Electric Corp | スターセンサ |
CN103675794A (zh) * | 2013-12-04 | 2014-03-26 | 北京空间机电研究所 | 基于时空统一特性的航天光学遥感器成像仿真方法 |
CN104301622A (zh) * | 2014-09-30 | 2015-01-21 | 上海微小卫星工程中心 | 一种天基空间目标成像曝光时间自适应控制方法及装置 |
CN106441237A (zh) * | 2015-08-10 | 2017-02-22 | 北京空间飞行器总体设计部 | 一种光学遥感卫星相机成像参数在轨自主调整方法 |
CN109581371A (zh) * | 2018-11-09 | 2019-04-05 | 中国科学院空间应用工程与技术中心 | 遥感相机成像参数的自动调整方法 |
CN110913150A (zh) * | 2019-11-18 | 2020-03-24 | 西北工业大学 | 基于空间平台的自适应曝光成像方法 |
Non-Patent Citations (2)
Title |
---|
一种空间相机快速自动曝光 技术;解静等;《光电工程》;第45卷(第11期);1-11 * |
空间相机自适应曝光;曹旗磊等;《光子学报》;第45卷(第8期);1-6 * |
Also Published As
Publication number | Publication date |
---|---|
CN115442534A (zh) | 2022-12-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110858286B (zh) | 一种用于目标识别的图像处理方法及装置 | |
CN104917975B (zh) | 一种基于目标特征的自适应自动曝光方法 | |
CN106373105B (zh) | 基于低秩矩阵恢复的多曝光图像去伪影融合方法 | |
CN111062876B (zh) | 用于校正模型训练和图像校正的方法和装置及电子设备 | |
US20150302575A1 (en) | Sun location prediction in image space with astronomical almanac-based calibration using ground based camera | |
CN100568278C (zh) | 结合帧选择和盲解卷积的自适应光学图像高分辨率复原方法 | |
WO2019127059A1 (zh) | 基于引导滤波和高通滤波的红外图像非均匀性校正方法 | |
CN112862728B (zh) | 伪影去除方法、装置、电子设备和存储介质 | |
JP2014241505A (ja) | 撮影装置および撮影方法、テンプレート生成装置およびテンプレート生成方法、並びにプログラム | |
Kinoshita et al. | Automatic exposure compensation for multi-exposure image fusion | |
CN115442534B (zh) | 空间目标天基成像自动多级曝光控制方法及系统 | |
CN102708559B (zh) | 使用多核卷积的模糊差估计 | |
CN108174055B (zh) | 一种智能监控方法、系统、设备及存储介质 | |
CN103868601A (zh) | Irfpa探测器非均匀响应的双边全变分正则化校正方法 | |
JP7167708B2 (ja) | 距離情報生成装置 | |
US6859550B2 (en) | Robust method for image feature estimation | |
CN110673428B (zh) | 一种结构光补偿方法、装置及设备 | |
JPWO2011086594A1 (ja) | 画像処理装置、及びその方法 | |
CN104166960B (zh) | 一种基于场景的自适应非均匀固定噪声去除方法 | |
CN115752760A (zh) | 一种适用于微振动环境的相位恢复算法 | |
CN109146966B (zh) | 视觉slam前端处理方法、系统、存储介质及计算机设备 | |
US9900491B2 (en) | Method and system for automatically focusing a camera | |
JP2019086325A (ja) | 日射量推定システム及び日射量推定方法 | |
CN114859099A (zh) | 一种无暗像元的可见光线阵探测器的暗电流实时估计方法 | |
CN113793276A (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 | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20230116 Address after: 200240 No. 3666, Yuanjiang Road, Minhang District, Shanghai Applicant after: SHANGHAI INSTITUTE OF SATELLITE ENGINEERING Applicant after: 63921 TROOPS OF PLA Address before: 200240 No. 3666, Yuanjiang Road, Minhang District, Shanghai Applicant before: SHANGHAI INSTITUTE OF SATELLITE ENGINEERING |
|
GR01 | Patent grant | ||
GR01 | Patent grant |