CN108226982B - 单线阵卫星激光联合高精度定位处理方法 - Google Patents

单线阵卫星激光联合高精度定位处理方法 Download PDF

Info

Publication number
CN108226982B
CN108226982B CN201711419908.6A CN201711419908A CN108226982B CN 108226982 B CN108226982 B CN 108226982B CN 201711419908 A CN201711419908 A CN 201711419908A CN 108226982 B CN108226982 B CN 108226982B
Authority
CN
China
Prior art keywords
image
linear array
laser
encryption
single linear
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
CN201711419908.6A
Other languages
English (en)
Other versions
CN108226982A (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.)
Space Star Technology Co Ltd
Original Assignee
Aerospace Tianhui Technology Co ltd
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 Aerospace Tianhui Technology Co ltd filed Critical Aerospace Tianhui Technology Co ltd
Priority to CN201711419908.6A priority Critical patent/CN108226982B/zh
Publication of CN108226982A publication Critical patent/CN108226982A/zh
Application granted granted Critical
Publication of CN108226982B publication Critical patent/CN108226982B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/48Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/87Combinations of systems using electromagnetic waves other than radio waves

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Image Processing (AREA)

Abstract

本发明提供单线阵卫星激光联合高精度定位处理方法,包括:利用单线阵卫星沿轨道前后摆动,以设定交会角获取地面目标的立体影像,及通过立体影像配准获得同名像点的像方坐标;利用初始定向参数和线阵影像像坐标计算加密点初始物方坐标;将激光测高数据与线阵影像配准得到激光点相应的像方坐标;构建误差方程式和对应的权矩阵;建立法方程,并答解影像定位模型参数和加密点地面坐标的改正数,并修正影像定位模型参数和加密点地面坐标;计算中误差,判断两次中误差之差是否小于给定限差,及根据判断结果进行迭代,并输出计算出的影像定位模型参数和加密点地面坐标。本发明在提高定位精度的同时,也能增强数据获取能力,保证遥感信息实时性。

Description

单线阵卫星激光联合高精度定位处理方法
技术领域
本发明属于单线阵卫星通信的技术领域,涉及一种单线阵卫星激光联合高精度定位处理方法。
背景技术
随着航天遥感信息实时性与精确性需求的提高,敏捷卫星成为近年来国际上发展十分迅速的一类高性能卫星。其所搭载的单线阵相机结构简单、体积小、质量轻,在实现高分辨率成像的同时,借助整星大范围、快速姿态机动能力,可以通过沿轨道方向前后摆动获取无明显时间差、无辐射差异的同轨立体像对。但是,与采用固定交会角的多线阵卫星相比,单线阵卫星为实现立体成像,姿态机动角度大且频繁,卫星姿态稳定度较差,使得卫星姿态往往带有很大的随机误差,也会直接影响定位精度。
另外,我国未来测绘卫星的发展趋势是实现无控条件上1:1万比例尺全球地形图的测制。在无控条件下,对于多线阵立体测绘卫星(如天绘一号和资源三号等),其高程精度也很难达到要求,则对单线阵卫星而言更是难上加难。于是,提出了把激光测高数据引入到线阵影像光束法平差中,但是最近几年对于线阵影像和激光数据联合处理的研究主要是针对的多线阵测绘卫星,对单线阵卫星与激光数据的研究尚缺乏。因此,开展单线阵卫星激光联合高精度定位技术的研究尤为重要和迫切。
发明内容
发明所要解决的课题是,如何将对单线阵卫星与激光数据联合,将激光测高数据引入到线阵影像光束法平差中以提高定位精度。
用于解决课题的技术手段是,提出单线阵卫星激光联合高精度定位处理方法,适用于单线阵卫星激光联合高精度定位的处理技术,该方法既能够提高高程定位精度,实现1:1万无控测图要求;也能够充分发挥单线阵相机的优势,保证高分辨率成像和遥感信息的实时性。
本发明的单线阵卫星激光联合高精度定位处理方法,包括如下步骤:
步骤(1)利用单线阵卫星沿轨道前后摆动,以设定交会角获取地面目标的立体影像,及通过立体影像配准获得同名像坐标,用于步骤(2)中计算加密点初始物方坐标;
步骤(2)根据卫星初始参数得到单线阵对应的初始外方位元素,利用初始外方位元素和单线阵影像同名像坐标计算加密点初始物方坐标;
步骤(3)将激光点物方坐标反投影到单线阵影像得到激光点像方坐标,用于步骤(4)中构建误差方程式;
步骤(4)依据有理多项式构建误差方程式和对应的权矩阵,这里构建误差方程式用到加密点物方坐标和像方坐标、激光点物方坐标和激光点像方坐标;
步骤(5)根据所建立的误差方程式和对应的权矩阵,建立法方程,并答解影像外方位元素和加密点物方坐标改正数,并修正影像外方位元素和加密点物方坐标;
步骤(6)根据所述修正的影像外方位元素和加密点地面坐标计算中误差,判断两次中误差之差是否小于给定限差,及根据判断结果进行迭代,并输出计算出的影像外方位元素和加密点物方坐标。
进一步地,作为本发明的一种优选技术方案,所述步骤(4)构建误差方程式和对应的权矩阵包括:
对于通过单线阵影像立体构像产生的连接点,按照有理多项式建立像点观测方程;
对于落在激光光斑范围内的连接点,建立高程约束观测方程。
进一步地,作为本发明的一种优选技术方案,所述步骤(6)计算中误差采用公式:
Figure GDA0002735980890000021
式中,RMS表示中误差,r表示多余观测;V表示误差方程式。
进一步地,作为本发明的一种优选技术方案,所述步骤(6)根据判断结果进行迭代,包括:
若判断两次中误差之差大于给定限差,则利用修正后的影像外方位元素和加密点物方坐标作为初始值重新组误差方程式计算其改正量;
若判断两次中误差之差小于给定限差,则停止迭代。
发明效果是,本发明与现有技术相比的优点在于:
(1)本发明中立体观测是通过单线阵相机沿轨道方向前后摆动实现的,充分发挥单线阵相机高分辨率成像和大角度姿态快速机动性能,使卫星具有更多的姿态工作模式,能够实现快速观测、拼幅成像、立体成像和推扫成像,保证了遥感信息的实时性,增强了数据获取能力,同时降低了成本。
(2)本发明中把激光测高数据引入到光束法区域网平差中,激光测高仪是一种能够精确获得地面三维信息的主动遥感手段,当激光测距精度足够高时,可以把其作为控制约束条件,参与影像的区域网平差计算,优化影像的外方位元素,提高地面点的高程精度,达到1:1万无控条件下地形图测制的精度要求。
附图说明
图1为本发明单线阵卫星激光联合高精度定位处理方法的流程图。
图2为本发明反投影中搜索窗口确定原理图。
图3(a)为本发明实验结果中激光测高数据对单线阵影像定位精度影像图;图3(b)为本发明实验结果中单线阵卫星激光联合平差结果图。
具体实施方式
以下,基于附图针对本发明进行详细地说明。
如图1所示,本发明设计了一种单线阵卫星激光联合高精度定位处理方法,该方法包括如下步骤:
步骤(1)利用单线阵卫星沿轨道前后摆动,以设定的一定交会角获取地面目标的立体影像,及通过立体影像配准获得同名像方坐标,用于步骤(2)中计算加密点初始物方坐标。
步骤(2)根据卫星初始参数得到单线阵影像对应的初始外方位元素,利用初始外方位元素和单线阵影像像坐标计算加密点初始物方坐标,加密点初始物方坐标的计算过程如下:
有理多项式公式如下:
Figure GDA0002735980890000041
式中,(Y,X)是指加密点像方坐标,(P,L,H)是指加密点物方坐标,NumL、DenL、NumS、DenS是关于(P,L,H)的多项式,计算公式如下:
NumL=a1+a2L+a3P+a4H+a5LP+a6LH+a7PH+a8L2+a9P2+a10H2+a11PLH+a12L3+a13LP2+a14LH2+a15L2P+a16P3+a17PH2+a18L2H+a19P2H+a20H3
DenL=b1+b2L+b3P+b4H+b5LP+b6LH+b7PH+b8L2+b9P2+b10H2+b11PLH+b12L3+b13LP2+b14LH2+b15L2P+b16P3+b17PH2+b18L2H+b19P2H+b20H3
Nums=c1+c2L+c3P+c4H+c5LP+c6LH+c7PH+c8L2+c9P2+c10H2+c11PLH+c12L3+c13LP2+c14LH2+c15L2P+c16P3+c17PH2+c18L2H+c19P2H+c20H3
Dens=d1+d2L+d3P+d4H+d5LP+d6LH+d7PH+d8L2+d9P2+d10H2+d11PLH+d12L3+d13LP2+d14LH2+d15L2P+d16P3+d17PH2+d18L2H+d19P2H+d20H3
其中,a1~a20、b1~b20、c1~c20、d1~d20是有理多项式模型参数,初始值可以从卫星自带数据中获得。
步骤(3)将激光点物方坐标反投影到单线阵影像得到激光点像方坐标,反投影计算激光像坐标的基本过程如下:
如图2所示:设相机焦平面上每条线阵CCD像素数为N,条带影像行数为l,地球平均高程面为H0,在轨道条带上选择a(0,N)、b(0,0)、c(l,0)、d(l,N)、e(l/2,N/2)共5个像点,利用地球平均高程面上对应的地面点A、B、C、D、E坐标,建立变换关系式(2),解算系数c0、c1、c2、c3
lm=c0+c1×lat+c2×lon+c3×lat×lon (2)
①当给定激光点物方坐标M(P,L,H)时,代入公式(2)计算搜索窗口[lm-δ,lm+δ],δ由立体相机基高比确定;
②搜索窗口内取初始扫描行序号i=lm-δ;
③根据扫描行序号i,获取对应行的外方位元素,计算立体相机焦平面坐标(xi,yi);
④计算扫描行增量Δi=(xi-x0)/ps,更新扫描行序号i=i+Δi(x0表示CCD线阵在焦平面上的横坐标,ps表示CCD像素大小);
⑤若Δi<δ(0.2-0.5像素)迭代终止,否则返回步骤③.
步骤(4)依据有理多项式构建误差方程式和对应的权矩阵,这里构建误差方程式用到加密点物方坐标和像方坐标、激光点物方坐标和激光点像方坐标,具体如下:
①对于通过单线阵影像立体构像产生的连接点,按照有理多项式建立像点观测方程:
由有理多项式公式(1)可得到:
Figure GDA0002735980890000051
则像点观测方程为:
Figure GDA0002735980890000052
简写成矩阵形式如下:
V1=A1t+B1X-L1 P1 (5)
式中:
V1=[VX VY]T
Figure GDA0002735980890000054
Figure GDA0002735980890000055
t=[Δai Δbi Δci Δdi]T i=1,…,20
X=[Δdi ΔP ΔL ΔH]T
②影像外方位元素的虚拟观测方程
在计算过程中,由于外方位元素解算参数过多,为了避免各未知数的相关性,每个影像的外方位元素均视为加权“观测值”,建立“虚拟误差方程”进行平差解算:
V2=Et-L2 P2 (6)
式中,E为单位矩阵,权值P2,通常依照卫星初始外方位元素确定。
③对于落在激光光斑范围内的连接点,建立高程约束观测方程
V3=B3Δh-L3 P3 (7)
式中,Δh为高程改正数,权值P3根据激光测距精度进行设定
综上所述:误差方程式为:
Figure GDA0002735980890000061
其中:
第一组方程是基于有理多项式的像点观测方程:A1和B1为像点计算值对未知数(影像外方位元素和加密点物方坐标)的偏导数,t为影像外方位元素的改正数,X为加密点物方坐标的改正数,L1为观测方程的常数项,P1为权阵,V1为残差。
第二组方程是影像外方位元素的虚拟观测方程,权值取决于卫星初始外方位元素的精度。
第三组方程是高程约束方程,对落入激光光斑范围内的连接点,其高程值应与激光测量值保持一致,权值取决于激光测高精度。
步骤(5)根据所建立的误差方程式和对应的权矩阵,建立法方程,并答解影像外方位元素和加密点物方坐标的改正数,并修正影像外方位元素和加密点物方坐标。
本步骤利用法化答解外方位元素和加密点物方坐标,具体如下:
根据最小二乘原理,误差方程式的法方程如下:
Figure GDA0002735980890000062
把法方程中的系数矩阵和常数项用新的符号代替,写成如下形式:
Figure GDA0002735980890000071
用消元法消去加密点物方坐标改正数得到改化法方程式,即:
Figure GDA0002735980890000072
则求解得到定向参数和加密点地面坐标改正数如下所示:
Figure GDA0002735980890000073
步骤(6)首先,求解得到未知数t和X,则可得到总的误差方程式V=At+BX-L.按照公式(13)计算中误差。根据所述修正的影像外方位元素和加密点物方坐标计算中误差,计算中误差采用公式:
Figure GDA0002735980890000074
式中,RMS表示中误差,r表示多余观测;V表示误差方程式。
然后,判断两次中误差之差是否小于给定限差,及根据判断结果进行迭代,并输出计算出的影像外方位元素和加密点物方坐标。若大于限差则返回(4)步,利用修正后的影像外方位元素和加密点物方坐标作为初始值重新组误差方程式计算其改正量,继续迭代使得下次迭代过程中的初始值取上次迭代完成修正后的值;若小于限差,则停止迭代,并输出计算出的影像外方位元素和加密点物方坐标。
为了验证本发明能够实现高程定位精度和保证高分辨率成像和遥感信息的实时性,特列举一验证例进行说明。
为了验证本发明联合定位的精度指标,仿真激光测高数据,并运用到单线阵卫星XM的区域网平差中,表1为仿真激光测高数据所用到的技术指标:
表1
激光脉冲频率 激光光斑直径 激光发射数量 激光脉冲夹角 激光测距精度 侧摆角
40HZ 20m 3 0.7° 0.5m 0~42°
根据表1中技术指标参数仿真得到3束在单线阵卫星影像上分布均匀的激光条带数据,然后将得到的激光数据作为高程约束加入到区域网平差中,其平面和高程中误差的变化如图3(a)和图3(b)所示:
综上可以看出:激光测高数据可以显著提高立体影像平差结果中的高程精度。根据当前设备技术指标,测高数据辅助的无控平差高程精度可以保持在2-3米的平均水平。但测高数据的约束对提升平面精度没有帮助。星载激光测高仪在卫星侧摆角较大的情况下会产生显著的测高误差。在侧摆角小于25度时,联合平差可以使高程精度提升30%以上;但是当侧摆角大于25度,高程精度的提升幅度大大下降,仅有15%-20%的水平。
综上,本发明通过激光测高数据约束线阵影像,优化影像的外方位元素,提高高程定位精度,从而达到无控条件下的精度要求。既能够提高高程定位精度,实现1:1万无控测图要求;也能够充分发挥单线阵相机的优势。保证了遥感信息的实时性,增强了数据获取能力,同时降低了成本。单线阵卫星激光联合高精度定位处理技术,将单线阵卫星影像数据和激光测高数据联合处理,在提高定位精度的同时,也能增强数据获取能力,保证遥感信息实时性。
需要说明的是,以上说明仅是本发明的优选实施方式,应当理解,对于本领域技术人员来说,在不脱离本发明技术构思的前提下还可以做出若干改变和改进,这些都包括在本发明的保护范围内。

Claims (4)

1.单线阵卫星激光联合高精度定位处理方法,其特征在于,包括如下步骤:
步骤(1)利用单线阵卫星沿轨道前后摆动,以设定交会角获取地面目标的立体影像,及通过立体影像配准获得同名像坐标,用于步骤(2)中计算加密点初始物方坐标;
步骤(2)根据卫星初始参数得到单线阵对应的初始外方位元素,利用初始外方位元素和单线阵影像同名像坐标计算加密点初始物方坐标;
步骤(3)将激光点物方坐标反投影到单线阵影像得到激光点像方坐标,用于步骤(4)中构建误差方程式;
步骤(4)依据有理多项式构建误差方程式和对应的权矩阵,这里构建误差方程式用到加密点物方坐标和像方坐标、激光点物方坐标和激光点像方坐标;
步骤(5)根据所建立的误差方程式和对应的权矩阵,建立法方程,并答解影像外方位元素和加密点物方坐标改正数,并修正影像外方位元素和加密点物方坐标;
步骤(6)所述修正的影像外方位元素和加密点地面坐标计算中误差,判断两次中误差之差是否小于给定限差,及根据判断结果进行迭代,并输出计算出的影像外方位元素和加密点物方坐标。
2.根据权利要求1所述单线阵卫星激光联合高精度定位处理方法,其特征在于,所述步骤(4)构建误差方程式和对应的权矩阵包括:
对于通过单线阵影像立体构像产生的连接点,按照有理多项式建立像点观测方程;
对于落在激光光斑范围内的连接点,建立高程约束观测方程。
3.根据权利要求1所述单线阵卫星激光联合高精度定位处理方法,其特征在于,所述步骤(6)计算中误差采用公式:
Figure FDA0001522751830000021
式中,RMS表示中误差,r表示多余观测;V表示误差方程式。
4.根据权利要求1所述单线阵卫星激光联合高精度定位处理方法,其特征在于,所述步骤(6)根据判断结果进行迭代,包括:
若判断两次中误差之差大于给定限差,则利用修正后的影像外方位元素和加密点物方坐标作为初始值重新组误差方程式计算其改正量;
若判断两次中误差之差小于给定限差,则停止迭代。
CN201711419908.6A 2017-12-25 2017-12-25 单线阵卫星激光联合高精度定位处理方法 Active CN108226982B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711419908.6A CN108226982B (zh) 2017-12-25 2017-12-25 单线阵卫星激光联合高精度定位处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711419908.6A CN108226982B (zh) 2017-12-25 2017-12-25 单线阵卫星激光联合高精度定位处理方法

Publications (2)

Publication Number Publication Date
CN108226982A CN108226982A (zh) 2018-06-29
CN108226982B true CN108226982B (zh) 2020-12-25

Family

ID=62648730

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711419908.6A Active CN108226982B (zh) 2017-12-25 2017-12-25 单线阵卫星激光联合高精度定位处理方法

Country Status (1)

Country Link
CN (1) CN108226982B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109239745A (zh) * 2018-09-11 2019-01-18 中铁二院工程集团有限责任公司 一种高分辨率卫星影像及有理多项式参数变换方法
CN110940966B (zh) * 2019-11-25 2021-09-03 同济大学 一种基于激光测高卫星足印影像的激光足印平面定位方法
CN112270714B (zh) * 2020-08-17 2022-10-28 中国人民解放军61540部队 卫星线阵立体相机参数在轨渐进标定方法
CN113324527B (zh) * 2021-05-28 2021-11-26 自然资源部国土卫星遥感应用中心 一种同轨激光测高点与三线阵立体影像联合测绘处理方法
CN116886178B (zh) * 2023-09-06 2024-01-19 北京融为科技有限公司 轨道预报修正方法及装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102519436B (zh) * 2011-12-28 2013-07-24 武汉大学 一种嫦娥一号立体相机与激光高度计数据联合平差方法
CN107014399B (zh) * 2017-04-25 2020-09-29 自然资源部国土卫星遥感应用中心 一种星载光学相机-激光测距仪组合***联合检校方法

Also Published As

Publication number Publication date
CN108226982A (zh) 2018-06-29

Similar Documents

Publication Publication Date Title
CN108226982B (zh) 单线阵卫星激光联合高精度定位处理方法
CN107341778B (zh) 基于卫星控制点库和dem的sar影像正射纠正方法
KR101965965B1 (ko) 위성영상과 제공 rpc로부터 제작된 수치표고모델의 자동 오차보정 방법
CN102506824B (zh) 一种城市低空无人机***生成数字正射影像图的方法
CN103823981B (zh) 一种数字高程模型辅助的卫星影像区域网平差方法
CN103679711B (zh) 一种遥感卫星线阵推扫光学相机在轨外方位参数标定方法
CN110388898B (zh) 构建虚拟控制点约束的多源多重覆盖遥感影像平差方法
CN109709551B (zh) 一种星载合成孔径雷达影像的区域网平面平差方法
CN113538595B (zh) 利用激光测高数据辅助提升遥感立体影像几何精度的方法
CN111145227B (zh) 一种地下隧道空间多视点云的可迭代整体配准方法
CN109696182A (zh) 一种星载推扫式光学传感器内方位元素定标方法
CN103673995A (zh) 一种线阵推扫式相机在轨光学畸变参数标定方法
CN113514829B (zh) 面向InSAR的初始DSM的区域网平差方法
Fraser et al. Sub-metre geopositioning with Ikonos GEO imagery
CN104807477B (zh) 一种基于靶标控制点的卫星ccd阵列影像几何检校方法
CN104820984B (zh) 一种卫星遥感三线阵立体影像处理***及方法
CN111473802A (zh) 一种基于线阵推扫的光学传感器内方位元素定标方法
CN102322863B (zh) 一种遥感卫星多星联合逆向定轨定姿方法
CN107330927A (zh) 机载可见光图像定位方法
CN108447100B (zh) 一种机载三线阵ccd相机的偏心矢量和视轴偏心角标定方法
CN111508028A (zh) 光学立体测绘卫星相机的自主在轨几何定标方法及***
CN110030968B (zh) 一种基于星载立体光学影像的地面遮挡物仰角测量方法
CN109188483B (zh) 一种时序化高精度外方位元素自动定标方法
CN113255740B (zh) 一种多源遥感影像平差定位精度分析方法
CN113532377B (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230425

Address after: 100086 No. 82, Haidian District, Beijing, Zhichun Road

Patentee after: SPACE STAR TECHNOLOGY Co.,Ltd.

Address before: Unit 1406, 14 / F, building 51, 63 Zhichun Road, Haidian District, Beijing 100190

Patentee before: AEROSPACE TIANHUI TECHNOLOGY CO.,LTD.