CN110109149B - 一种激光雷达低层消光系数廓线校正方法 - Google Patents

一种激光雷达低层消光系数廓线校正方法 Download PDF

Info

Publication number
CN110109149B
CN110109149B CN201910378942.6A CN201910378942A CN110109149B CN 110109149 B CN110109149 B CN 110109149B CN 201910378942 A CN201910378942 A CN 201910378942A CN 110109149 B CN110109149 B CN 110109149B
Authority
CN
China
Prior art keywords
extinction coefficient
laser radar
profile data
coefficient profile
absorption spectrometer
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
CN201910378942.6A
Other languages
English (en)
Other versions
CN110109149A (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.)
Nanjing University of Information Science and Technology
Original Assignee
Nanjing University of Information Science and Technology
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 Nanjing University of Information Science and Technology filed Critical Nanjing University of Information Science and Technology
Priority to CN201910378942.6A priority Critical patent/CN110109149B/zh
Publication of CN110109149A publication Critical patent/CN110109149A/zh
Application granted granted Critical
Publication of CN110109149B publication Critical patent/CN110109149B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88—Lidar systems specially adapted for specific applications
    • G01S17/95—Lidar systems specially adapted for specific applications for meteorological use
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/497—Means for monitoring or calibrating
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Optical Radar Systems And Details Thereof (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种激光雷达低层消光系数廓线校正方法，先采用同时采用多轴差分吸收光谱仪与激光雷达分别测量出大气低层的气溶胶消光系数数据；然后分别对多轴差分吸收光谱仪与激光雷达所测量的气溶胶消光系数数据进行计算得到多轴差分吸收光谱仪测量的消光系数廓线数据、激光雷达测量的消光系数廓线数据；最后对多轴差分吸收光谱仪测量的消光系数廓线数据进行处理后作为多轴差分吸收光谱仪测量的大气消光系数廓线数据，对之前得到的激光雷达测量的消光系数廓线数据进行校正，得到校正后的激光雷达低层消光系数廓线。本发明具有方法简单、高效、能够减少低层消光廓线数据误差与减少测量成本等优点。

Description

一种激光雷达低层消光系数廓线校正方法
技术领域
本发明属于大气科学探测技术领域，具体涉及一种激光雷达低层消光系数廓线校正方法。
背景技术
气溶胶是指悬浮在大气中的固体或者液体颗粒物，对大气系统中起着重要作用。气溶胶粒子散射并吸收辐射，会对云的光学特性、寿命和形成产生影响。直接排放的烟尘颗粒，以及二次有机气溶胶的形成和凝结和气溶胶参与形成的有害气体会影响空气质量和人体健康。随着人类的生活水平的提高，由于人为因素产生的气溶胶总量增加，对我们生存环境产生了严重的危害。大部分气溶胶来源于地面，对流层中主要的气溶胶含量比例极高，因此底层和对流层的气溶胶含量很高，底层气溶胶含量对于大气中的云和物理化学过程产生很大的影响，以及在大气传输方面有很大影响。因此加强边界层和对流层低空的气溶胶分布和特性研究具有重要意义。目前，国内采用的气溶胶特性探测技术包括地基被动遥感的太阳光度计，以及卫星遥感探测，和主动的地基遥感激光雷达技术。激光雷达采用主动遥感的方式，因为其抗干扰能力强，激光波长短，发散角小，灵敏度和分辨率高的优点被广泛应用于大气气溶胶特性的研究。激光雷达相较太阳光度计可以探测大气气溶胶的垂直廓线信息，是一种有效的大气气溶胶探测方式。但是，激光雷达在探测大气气溶胶特性的过程中会因为系统本身的限制给测量结果造成无法避免的系统误差，例如由于激光雷达激光的发散角与望远镜接受视场部分重合时，发射激光仅有一部分被有效利用，称作激光雷达的几何重叠因子，由于这种系统误差会在探测的近距离处产生测量误差。在大气气溶胶消光特性廓线探测过程中，低层大气探测过程中会出现测量误差区，所以校正激光雷达低层消光廓线的误差对于研究气溶胶分布的特征具有重要意义。
发明内容
本发明的目的是为了解决现有技术中存在的缺陷，提供一种激光雷达低层消光系数廓线校正方法。
本发明为解决上述技术问题采用以下技术方案：
一种激光雷达低层消光系数廓线校正方法，包括以下步骤：
步骤1、同时采用多轴差分吸收光谱仪与激光雷达分别测量出大气低层的气溶胶消光系数数据；
步骤2、分别对多轴差分吸收光谱仪与激光雷达所测量的气溶胶消光系数数据进行计算得到多轴差分吸收光谱仪测量的消光系数廓线数据、激光雷达测量的消光系数廓线数据；
步骤3、对多轴差分吸收光谱仪测量的消光系数廓线数据进行处理后作为多轴差分吸收光谱仪测量的大气消光系数廓线数据，对步骤2中得到的激光雷达测量的消光系数廓线数据进行校正，得到校正后的激光雷达低层消光系数廓线。
作为本发明的优选技术方案：步骤2中的激光雷达测量的消光系数廓线数据的计算方法为基于激光雷达测量出的大气低层的气溶胶消光系数数据，利用Fernald算法反演激光雷达所测得的气溶胶消光系数廓线数据，拟合出激光雷达测量的消光系数廓线数据。
作为本发明的优选技术方案：步骤2中的多轴差分吸收光谱仪测量的消光系数廓线数据的计算方法为基于多轴差分吸收光谱仪测量出的大气低层的气溶胶消光系数数据，利用最优估算方法与在大气辐射传输模式下，插值拟合得到多轴差分吸收光谱仪测量的消光系数廓线数据。
作为本发明的优选技术方案：步骤3的具体步骤如下：
步骤3.1、采集激光雷达的回波功率信号P(R)，当信号P(R)为峰值时对应的激光雷达探测高度为r0，以r0高度作为处理消光系数廓线数据的临界高度；
步骤3.2、对于r0高度以上的多轴差分吸收光谱仪与激光雷达所测得的消光系数廓线数据进行匹配检验，用以下公式计算：
|σ1(ri)-σ2(ri)|≤ε
式中，σ1(ri)为激光雷达在高度ri所测量的气溶胶消光系数廓线数据，σ2(ri)为多轴差分光谱仪在高度ri所测量的气溶胶消光系数廓线数据，ε取值足够满足消光匹配条件，高度ri>高度 r0；
步骤3.3、对于r0高度以下的高度范围，将激光雷达空间分辨率Δr的高度分为多个区间，采用对多轴差分吸收光谱仪测得的相邻的多个区间的消光系数廓线数据取平均的方法，得到多轴差分吸收光谱仪测量处理后的消光系数廓线数据，作为多轴差分吸收光谱仪测量的大气消光系数廓线数据；
步骤3.4、将多轴差分吸收光谱仪测量的大气消光系数廓线数据来近似代替激光雷达空间分辨率Δr内的一个消光系数廓线数据，从而将激光雷达测量的消光系数廓线数据进行校正，得到激光雷达低层消光系数廓线。
本发明所述激光雷达低层消光系数廓线校正方法，采用以上技术方案与现有技术相比，具有以下技术效果：
本发明提出了一种对激光雷达探测的由于几何重叠因子所产生的低层消光廓线数据误差校正方法，利用简单的、小型的多轴差分吸收光谱仪测得的大气低层气溶胶消光廓线数据，通过数据匹配和数据处理方法，用来修正激光雷达由于几何重叠因子误差所造成的低层气溶胶消光廓线误差。通过两种仪器的数据处理可以修正激光雷达低层的气溶胶消光廓线信息，具有方法简单、高效、能够减少低层消光廓线数据误差与减少测量成本等优点。
附图说明
图1为激光雷达大气低层气溶胶消光廓线校正方法流程图；
图2为原始未校正前激光雷达的消光廓线图；
图3为经过校正之后的激光雷达消光廓线图。
具体实施方式
下面结合附图对本发明做进一步的详细说明。
如图1所示，为激光雷达大气低层气溶胶消光廓线校正方法流程图，首先分别通过激光雷达和多轴差分吸收光谱仪进行大气气溶胶消光廓线数据的探测。
基于激光雷达测量数据获取大气气溶胶消光系数数据，气溶胶消光系数反演一般建立在激光雷达方程基础上进行反演，在激光雷达垂直探测气溶胶消光廓线数据的方法中，经过前人的发展和完善，使用Fernald方法反演，其中多用后向积分方法，气溶胶垂直消光系数的后向积分方程可表示为：
Figure GDA0003779022460000031
其中，X(z)＝P(z)z2，为经过距离校正的回波信号，σ1(z)，σ2(z)分别表示气溶胶和大气分子在高度z处的消光系数，β1(z)，β2(z)分别表示气溶胶和大气分子在z高度处的后向散射系数，Zc为计算时用到的标定高度，σ1(zc)，σ2(zc)分别为标定高度处的气溶胶和大气分子的消光系数，S1＝σ1(z)/β1(z)，
Figure GDA0003779022460000032
分别为气溶胶和大气分子的消光散射比，由上式可知，只要β1(z)，Zc，S1，S2已知，就可以由Fernald后向算法逐步迭代求解得到气溶胶的后向散射系数和消光系数。标定高度一般以标准大气后向散射系数或消光系数为标准，通常会引入散射比，即
Figure GDA0003779022460000041
取散射比R为1.01，这样可以通过标准大气分子的后向散射系数计算得到气溶胶的后向散射系数β2(zc)。对于气溶胶消光散射比S1取为 50，标定高度在白天取3-5km。
再利用多轴差分吸收光谱技术求大气气溶胶消光廓线数据，基于最优估算法和大气辐射传输模式的方法求取大气的气溶胶的消光廓线数据，通过接收大气被动太阳辐射，根据大气分子对太阳辐射的吸收能力不同，来进行大气气溶胶消光廓线反演。把多轴差分光谱仪与激光雷达放置在同一片大气区域。通过最优估算法和大气辐射传输模式得到气溶胶消光数据，因为多轴差分吸收光谱以反演波段为477nm，而激光雷达反演气溶胶消光系数所用波段为 532nm，因此需要对差分吸收光谱测量的消光系数进行转化。
利用以下公式求解532nm消光系数：
Figure GDA0003779022460000042
其中A为Angstrom指数，通过477nm波段消光系数得到532nm处消光系数值，表示为σ2(z)。
最后对多轴差分吸收光谱仪的消光数据进行处理，并且将处理好的多轴差分吸收光谱仪测量的消光数据作为大气的消光系数廓线分布，然后将数据应用到激光雷达的消光廓线修正之中。
通过激光雷达的回波功率信号，当信号P(R)为峰值时对应的激光雷达探测高度为r0，以及激光雷达探测时的空间分辨率为Δr。激光雷达峰值信号所对应的高度为r0，激光雷达的大气气溶胶消光随高度分布表示为σ1(z)，多轴差分吸收光谱仪的消光系数用σ2(z)表示，对于 r0高度以上的多轴差分吸收光谱仪和激光雷达测得消光系数数值，首先是对r0以上的消光数值进行匹配检验，用下面公式计算：
|σ1(z)-σ2(z)|＝ε
其中，ε取值足够满足消光匹配条件；ε取值足够小时，确定高度r0以上的气溶胶消光数值匹配良好。
对于r0高度以下的消光数值，对多轴差分吸收光谱仪的数据进行处理，对于r0高度以下的多轴差分吸收光谱仪所测得消光数据，对高度进行分层，激光雷达的空间分辨率为Δr，取高度分层数为：
j＝r0/3Δr
取相邻的三个高度层的消光系数取平均，所得平均值来代替激光雷达的大气消光数据。
计算如下：
Figure GDA0003779022460000051
得到1到
Figure GDA0003779022460000052
个消光数据点，将利用多轴差分吸收光谱仪测得的消光数据作为激光雷达在高度r0以下的气溶胶消光数值，从而可以将激光雷达因为系统误差导致的消光廓线信息校正，得到大气真实气溶胶消光廓线。
具体工作中，选取一次激光雷达数据和MAXDOAS数据进行分析如下：激光雷达校正前所求得的消光数据廓线对应结果如附图2所示，其中黑色虚线为所对应选取的激光雷达回波信号P(R)为峰值时对应的激光雷达探测高度r0，r0对应探测高度为0.652km。以此高度作为分界，按照本发明进行计算处理，对应r0高度下激光雷达消光廓线进行校正。所得校正后激光雷达消光廓线如附图3所示。结果表明，按照本发明可以对激光雷达低层由于系统误差而导致的消光廓线误差进行校正。
以上所述仅为本发明的具体实施方案而已，并非用以限定本发明的范围，任何本领域的技术人员，在不脱离本发明的构思和原则的前提下所做出的等同变化与修改，均应属于本发明保护的范围。

Claims (3)

1.一种激光雷达低层消光系数廓线校正方法，其特征在于，包括以下步骤：
步骤1、同时采用多轴差分吸收光谱仪与激光雷达分别测量出大气低层的气溶胶消光系数数据；
步骤2、分别对多轴差分吸收光谱仪与激光雷达所测量的气溶胶消光系数数据进行计算得到多轴差分吸收光谱仪测量的消光系数廓线数据、激光雷达测量的消光系数廓线数据；
步骤3、对多轴差分吸收光谱仪测量的消光系数廓线数据进行处理后作为多轴差分吸收光谱仪测量的大气消光系数廓线数据，对步骤2中得到的激光雷达测量的消光系数廓线数据进行校正，得到校正后的激光雷达低层消光系数廓线；
步骤3的具体步骤如下：
步骤3.1、采集激光雷达的回波功率信号P(R)，当信号P(R)为峰值时对应的激光雷达探测高度为r0，以r0高度作为处理消光系数廓线数据的临界高度；
步骤3.2、对于r0高度以上的多轴差分吸收光谱仪与激光雷达所测得的消光系数廓线数据进行匹配检验，用以下公式计算：
|σ1(ri)-σ2(ri)|≤ε
式中，σ1(ri)为激光雷达在高度ri所测量的气溶胶消光系数廓线数据，σ2(ri)为多轴差分光谱仪在高度ri所测量的气溶胶消光系数廓线数据，ε取值足够满足消光匹配条件，高度ri>高度r0；
步骤3.3、对于r0高度以下的高度范围，将激光雷达空间分辨率Δr的高度分为多个区间，采用对多轴差分吸收光谱仪测得的相邻的多个区间的消光系数廓线数据取平均的方法，得到多轴差分吸收光谱仪测量处理后的消光系数廓线数据，作为多轴差分吸收光谱仪测量的大气消光系数廓线数据；
步骤3.4、将多轴差分吸收光谱仪测量的大气消光系数廓线数据来近似代替激光雷达空间分辨率Δr内的一个消光系数廓线数据，从而将激光雷达测量的消光系数廓线数据进行校正，得到激光雷达低层消光系数廓线。
2.根据权利要求1所述的校正方法，其特征在于，步骤2中的激光雷达测量的消光系数廓线数据的计算方法为基于激光雷达测量出的大气低层的气溶胶消光系数数据，利用Fernald算法反演激光雷达所测得的气溶胶消光系数廓线数据，拟合出激光雷达测量的消光系数廓线数据。
3.根据权利要求1所述的校正方法，其特征在于，步骤2中的多轴差分吸收光谱仪测量的消光系数廓线数据的计算方法为基于多轴差分吸收光谱仪测量出的大气低层的气溶胶消光系数数据，利用最优估算方法与在大气辐射传输模式下，插值拟合得到多轴差分吸收光谱仪测量的消光系数廓线数据。
CN201910378942.6A 2019-05-08 2019-05-08 一种激光雷达低层消光系数廓线校正方法 Active CN110109149B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910378942.6A CN110109149B (zh) 2019-05-08 2019-05-08 一种激光雷达低层消光系数廓线校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910378942.6A CN110109149B (zh) 2019-05-08 2019-05-08 一种激光雷达低层消光系数廓线校正方法

Publications (2)

Publication Number Publication Date
CN110109149A CN110109149A (zh) 2019-08-09
CN110109149B true CN110109149B (zh) 2022-11-25

Family

ID=67488786

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910378942.6A Active CN110109149B (zh) 2019-05-08 2019-05-08 一种激光雷达低层消光系数廓线校正方法

Country Status (1)

Country Link
CN (1) CN110109149B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110850393B (zh) * 2019-10-31 2023-01-31 无锡中科光电技术有限公司 激光雷达数据的验证方法
CN111220511B (zh) * 2020-01-21 2020-10-27 中国科学院大气物理研究所 基于普适吸湿增长方案的灰霾消光监测方法
CN115032655B (zh) * 2022-08-10 2022-11-29 安徽科创中光科技股份有限公司 基于可见光和紫外光的双波长激光雷达反演臭氧浓度算法
CN116466368B (zh) * 2023-06-16 2023-08-22 成都远望科技有限责任公司 基于激光雷达和卫星资料的沙尘消光系数廓线估算方法
CN118033676A (zh) * 2024-04-11 2024-05-14 长春理工大学 基于bp神经网络的激光雷达消光系数廓线的修正方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106596355A (zh) * 2016-12-30 2017-04-26 无锡中科光电技术有限公司 一种激光雷达反演中低云层下方消光系数的修正方法
KR101751642B1 (ko) * 2016-06-08 2017-06-28 강릉원주대학교산학협력단 대기 관측용 라이다의 소산계수 보정방법
CN107272023A (zh) * 2017-06-14 2017-10-20 中国人民解放军92232部队 一种测量海平面非均匀水平大气气溶胶消光系数的装置
CN108896456A (zh) * 2018-04-28 2018-11-27 南京信息工程大学 基于反馈型rbf神经网络的气溶胶消光系数反演方法
CN109596594A (zh) * 2018-11-27 2019-04-09 南京信息工程大学 基于拉曼-米散射激光雷达的气溶胶消光系数反演方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101751642B1 (ko) * 2016-06-08 2017-06-28 강릉원주대학교산학협력단 대기 관측용 라이다의 소산계수 보정방법
CN106596355A (zh) * 2016-12-30 2017-04-26 无锡中科光电技术有限公司 一种激光雷达反演中低云层下方消光系数的修正方法
CN107272023A (zh) * 2017-06-14 2017-10-20 中国人民解放军92232部队 一种测量海平面非均匀水平大气气溶胶消光系数的装置
CN108896456A (zh) * 2018-04-28 2018-11-27 南京信息工程大学 基于反馈型rbf神经网络的气溶胶消光系数反演方法
CN109596594A (zh) * 2018-11-27 2019-04-09 南京信息工程大学 基于拉曼-米散射激光雷达的气溶胶消光系数反演方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于多轴差分吸收光谱技术的查找表法反演气溶胶消光廓线研究;吴丰成 等;《光学学报》;20130630;第33卷(第6期);第0601002-1至0601002-9页 *

Also Published As

Publication number Publication date
CN110109149A (zh) 2019-08-09

Similar Documents

Publication Publication Date Title
CN110109149B (zh) 一种激光雷达低层消光系数廓线校正方法
CN109596594B (zh) 基于拉曼-米散射激光雷达的气溶胶消光系数反演方法
CN110441777B (zh) 一种基于激光雷达的气溶胶垂直廓线的反演方法
CN103698305B (zh) 一种实时观测大气透射率的方法与系统
KR101751642B1 (ko) 대기 관측용 라이다의 소산계수 보정방법
CN102662162B (zh) X波段多普勒雷达降水估测方法
Liu et al. Transport, vertical structure and radiative properties of dust events in southeast China determined from ground and space sensors
CN112269189B (zh) 一种利用单波长激光雷达探测气溶胶质量浓度廓线的方法
CN110687020B (zh) 基于四聚氧吸收反演气溶胶光学特性的测量方法及其装置
CN112505651B (zh) 一种用于大气探测激光雷达的自动化处理方法
CN113341432B (zh) 基于激光雷达卫星的地基激光雷达气溶胶反演方法及系统
CN105044039B (zh) 一种根据激光雷达数据自动反演水平能见度的方法
CN113777579B (zh) 一种拉曼-米散射激光雷达气溶胶消光系数廓线反演算法
Sayer et al. In-situ and remotely-sensed observations of biomass burning aerosols at Doi Ang Khang, Thailand during 7-SEAS/BASELInE 2015
CN111679293A (zh) 一种激光雷达质量控制方法
CN114624731A (zh) 基于偏振遥感数据的云层上方气溶胶光学厚度反演方法
CN111650128B (zh) 一种基于地表反射率库的高分辨率大气气溶胶反演方法
CN117826635A (zh) 车载激光雷达雾天探测能力仿真方法、系统、终端及介质
Sullivan et al. Optimization of the GSFC TROPOZ DIAL retrieval using synthetic lidar returns and ozonesondes–Part 1: Algorithm validation
CN108646231A (zh) 一种激光雷达信号过渡区的校正方法
CN116136590A (zh) 一种标定多普勒激光雷达焦距的方法、装置及存储介质
CN111366254A (zh) 一种大气透过率检测方法及装置
CN110850392B (zh) 激光雷达退偏增益比的校正方法
Chiang et al. An iterative calculation to derive extinction-to-backscatter ratio based on lidar measurements
CN112649335B (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