CN111044489B - 一种基于多波长测量获得大气折射率高度分布剖面的方法 - Google Patents

一种基于多波长测量获得大气折射率高度分布剖面的方法 Download PDF

Info

Publication number
CN111044489B
CN111044489B CN201911132808.4A CN201911132808A CN111044489B CN 111044489 B CN111044489 B CN 111044489B CN 201911132808 A CN201911132808 A CN 201911132808A CN 111044489 B CN111044489 B CN 111044489B
Authority
CN
China
Prior art keywords
refractive index
measurement
height distribution
atmospheric
estimated
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
CN201911132808.4A
Other languages
English (en)
Other versions
CN111044489A (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.)
No63686 Troops Pla
Original Assignee
No63686 Troops Pla
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 No63686 Troops Pla filed Critical No63686 Troops Pla
Priority to CN201911132808.4A priority Critical patent/CN111044489B/zh
Publication of CN111044489A publication Critical patent/CN111044489A/zh
Application granted granted Critical
Publication of CN111044489B publication Critical patent/CN111044489B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00—Features of devices classified in G01N21/00
    • G01N2201/12—Circuits of general importance; Signal processing
    • G01N2201/127—Calibration; base line adjustment; drift compensation
    • 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)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Biochemistry (AREA)
  • Mathematical Optimization (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Mathematical Analysis (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明涉及航天器测控技术领域，提供了一种基于多波长测量获得大气折射率高度分布剖面的方法，包括：使用光电外测设备在不同波长下跟踪预定恒星，输出并记录测量系的测角数据；获取相同时间段恒星地固系精确星历数据；将恒星地固系精确星历数据转换为外测设备测量系的理论测元数据；以大气折射率高度分布的Hopfield模型的模式参数为待估未知数，建立模式参数待估方程组；对模式参数待估方程组进行求解，获得模式参数的数值解；将求解的模式参数代入模型获得大气折射率高度分布剖面。本发明具有成本低廉、精度高、适应性好等优势，能够获得不同情况下的大气折射率高度分布剖面，可以作为一种获得大气折射率高度分布剖面常规方法。

Description

一种基于多波长测量获得大气折射率高度分布剖面的方法
技术领域
本发明涉及航天测控技术领域，具体涉及一种基于多波长测量获得大气折射率高度分布剖面的方法。
背景技术
相比与其他测量手段，光学测量技术具有一些突出的优势：测量精度较高、可测量目标的光辐射特性参数、直观性较强并且可以进行事后复现以及不受“黑障”和地面杂波的干扰等。因而光学测量成为空间目标测量的重要手段之一，值得进行深入研究。由于地球大气折射率在高度上的分布不均匀，光波在大气中的传输路线将发生弯曲，产生大气折射误差。因此，大气折射误差是空间目标观测过程中的重要误差，需进行修正以获得高精度的探测数据。通常的，要对大气折射误差进行精确修正，首先需要获得大气折射率高度分布剖面，然后通过高斯分层积分计算进行误差修正。获得大气折射率高度分布剖面的方法有多种，一类是使用探空气球或者GPS信号反演进行实地测量，该方法可以比较准确地获得大气折射率的高度分布，误差修正精度较高，但是成本较高、计算量大；另一类是根据历史数据总结获得的大气折射率高度分布模型，该方法的模型适用性广、操作简单、成本较低，但精度有限，且难以适应随时变化的大气环境。因此，需要研究一种成本低、精度高、实时性好的获得大气折射率高度分布剖面方法。引用文献：
[1]任席闯等，红外与微波波段大气折射率比较分析，舰船科学技术,2011，33(1)：第65-68页。
[2]Stephen D.H.,The geometrical optics approach to atmosphericpropagation models,Proc.of SPIE,2003,4976:60-69。
[3]王旭良等，一种基于海上探空气象数据的大气折射率模型，飞行器测控学报，2015，34(1)：第36-41页。
萧安家，光波大气折射率计算公式，测绘学报，1988，17(2)：第129-135页。
发明内容
本发明的目的在于，克服现有技术中存在的缺陷，提供成本低、实时性好、精度高的获得大气折射率高度剖面的方法。
为实现上述目的，本发明的技术方案是设计一种基于多波长测量获得大气折射率高度分布剖面的方法，包括以下步骤：
S1、使用光电外测设备在不同波长下跟踪一颗预定恒星，输出并记录外测设备测量系的测角数据Eci(i＝1,...,n)，n为波长数目，同时获取相同时间段恒星地固系的精确星历数据(xei,yei,zei)；
S2、将恒星地固系的精确星历数据(xei,yei,zei)转换为外测设备测量系的理论测元数据Egi；
S3、以大气折射率高度分布的(Hopfield)霍普菲尔德模型的模式参数为待估未知数，利用高斯分层积分法计算外测设备的测角误差，从而建立模式参数待估方程组；
S4、对大气折射高度分布模式参数待估方程组进行求解，获得模式参数的数值解；
S5、利用数值求解的模式参数代入大气折射率高度分布的Hopfield模型，获得大气折射率高度分布剖面。
所述步骤S2具体操作如下：
S21、将恒星地固系的精确星历数据(xei,yei,zei)转换至外测设备测量系，得到精确星历直角坐标的数据
Figure GDA0003702251960000021
其中
Figure GDA0003702251960000031
Figure GDA0003702251960000032
为测站的天文经纬度。
S22、将恒星的精确星历直角坐标的数据(xci,yci,zci)从外测设备测量系的直角坐标系转换至极坐标系，得到Egi
Figure GDA0003702251960000033
其中，Egi为目标恒星的真实高度角；
所述步骤S3具体操作如下：
S31、大气折射率高度分布的Hopfield模型为：
N(h)＝Nd(h)+Nw(h)
Figure GDA0003702251960000034
j＝d，w，h≤Hj
Nj(h)＝0，j＝d，w，h＞Hj (3)
其中，h为高度，hs为折射率分布起始高度，Nd(h)和Nw(h)分别为高度h处干项和湿项折射率，且对于世界各地平均剖面而言，Hd＝40136+148.72t0，单位为米，Hw＝11000米，t0为地表温度；
S32、在可见光与近红外波段的大气折射率可以表示为：
Figure GDA0003702251960000035
其中，P(h)为气压，T(h)为绝对温度，e(h)为水汽压，m1(λ)和m2(λ)分别为波长λ(μm)的函数，
由下式确定：
Figure GDA0003702251960000036
S33、对比式(4)和(5)，可以建立关系式，对于波长为λi(i＝1,...,n)的光波，其对应的折射率分布有：
Figure GDA0003702251960000041
选择i＝1的波长λ1对应的折射率参数(Nd01,Nw01,Hd1,Hw1)作为待求参数，则其他波长对应的折射率参数(Nd0i,Nw0i,Hdi,Hwi)(i＝2,...,n)可以表示为
Figure GDA0003702251960000042
即
Figure GDA0003702251960000043
这样就可以将不同波长对应的折射率表示为折射率参数
Figure GDA0003702251960000047
的函数：
Figure GDA0003702251960000044
S34、根据折射率高度分布剖面，计算不同波长下目标恒星的等效地心角
Figure GDA0003702251960000045
Figure GDA0003702251960000046
其中，n0i测站地面大气折射率,h0为地球半径，Eci为目标恒星的视在高度角，ni(h)＝1+Ni(h)为h处大气折射指数，hm为目标恒星的等效地心距；
S35、由等效地心角
Figure GDA0003702251960000051
计算目标恒星的测量高度角Emi
Figure GDA0003702251960000052
S36、建立以模式参数为未知量(Nd01,Nw01,Hd1,Hw1)的参数待估方程组：
Egi-fi(Nd01,Nw01,Hd1,Hw1)＝0,(i＝1,...,n) (12)
其中，由于待估参数数目为4，所以n≥4，fi(Nd01,Nw01,Hd1,Hw1)是Emi关于模式参数(Nd01,Nw01,Hd1,Hw1)的函数表达形式，其计算公式表示为：
Figure GDA0003702251960000053
其中，(Nd01,Nw01,Hd1,Hw1)即为待求的优化模式参数。
所述步骤S4具体操作如下：
以模式参数的世界各地平均值(Nd0r,Nw0r,Hdr,Hwr)作为参考值，对式(12)的方程组进行数值求解，在求解误差小于设定值ε时，即认为获得了有效的模式参数(Nd01,Nw01,Hd1,Hw1)。
所述步骤S5具体操作如下：
根据模式参数的计算值(Nd01,Nw01,Hd1,Hw1)，代入式(3)则可以得到波长λ1对应的大气折射率高度分布的Hopfield模型，其它波长对应的折射率则可以根据式(9)得出。
本发明的优点和有益效果在于：该发明实现了基于多波长测量获得大气折射率高度分布剖面的方法。在每次探测工作开始前，利用恒星精确实测数据对大气折射率高度分布模式参数进行修正，既避免了使用探空气球、GPS信号反演等常规高成本测量手段，又可以获得特定区域、特定时间段大气折射率高度分布剖面，精度较高、成本低廉、操作简单、适应性好。因此，本发明能够因地制宜地获得不同情况下大气折射率高度分布剖面，可以作为一种获得大气折射率高度分布剖面常规方法。
本发明提出了基于多波长测量获得大气折射率高度分布剖面的方法。发明思路如下：目前应用广泛的Hopfield(霍普菲尔德)模型可以较好地表示大气折射率高度分布剖面，该模型与世界各地平均折射率剖面符合较好，但是与局部地区的折射率剖面存在一定差异。考虑利用实地测量的方法对上述模型的参数进行优化修正，以适应局部地区的折射率剖面。对于光学外测设备，通常可以在多个波长对目标进行观测。尽管不同光波波长对应的折射率是不同的，但是根据折射率的计算公式可以发现在其他条件相同的情况下不同光波波长折射率与波长直接相关。因而可以利用对同一目标恒星在不同波长的测量结果建立方程组，对Hopfield模型的参数进行优化修正，从而获得适应当地当时局部大气环境的大气折射率高度分布剖面。由于该分布剖面是基于实地观测后计算得到的，因而相比于标准的Hopfield模型精度更高。而且由于仅仅需要采用不同波长对一颗预定恒星进行观测，不需要额外设备，因而该方法操作简单、成本低廉，是一种较好地获得高精度大气折射率高度分布剖面的方法。
附图说明
图1是本发明基于多波长测量获得大气折射率高度分布剖面方法的实施流程图。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
如图1所示，本发明是一种基于多波长测量获得大气折射率高度分布剖面的方法包括以下步骤：
S1、使用光电外测设备在不同波长下跟踪一颗预定恒星，输出并记录外测设备测量系的测角数据Eci(i＝1,...,n)，n为波长数目。同时获取相同时间段恒星地固系的精确星历数据(xei,yei,zei)；
S2、将恒星地固系的精确星历数据(xei,yei,zei)转换为外测设备测量系的理论测元数据Egi；具体操作如下：
S21、将恒星地固系的精确星历数据(xei,yei,zei)转换至外测设备测量系，得到精确星历直角坐标的数据
Figure GDA0003702251960000071
其中
Figure GDA0003702251960000072
Figure GDA0003702251960000073
为测站的天文经纬度，
Figure GDA0003702251960000074
为测站原点在地固系中的位置坐标。
S22、将恒星的精确星历直角坐标的数据(xci,yci,zci)从外测设备测量系的直角坐标系转换至极坐标系，得到Egi
Figure GDA0003702251960000075
其中，Egi为目标恒星的真实高度角。
S3、以大气折射率高度分布的Hopfield模型的模式参数为待估未知数，利用高斯分层积分法计算外测设备的测角误差，从而建立模式参数待估方程组；具体操作如下：
S31、大气折射率高度分布的Hopfield模型为
N(h)＝Nd(h)+Nw(h)
Figure GDA0003702251960000081
j＝d，w，h≤Hj
Nj(h)＝0，j＝d，w，h＞Hj (3)
其中，h为高度，hs为折射率分布起始高度，Nd(h)和Nw(h)分别为高度h处干项和湿项折射率，Nd0和Nw0分别为地表的湿项和干项折射率，Hd和Hw分别为湿项和干项折射率的等效高度，且对于世界各地平均剖面而言，Hd＝40136+148.72t0，单位为米，Hw＝11000米，t0为地表温度。
S32、在可见光与近红外波段的大气折射率可以表示为：
Figure GDA0003702251960000082
其中，P(h)为气压，T(h)为绝对温度，e(h)为水汽压，m1(λ)和m2(λ)分别为波长λ(μm)的函数，
由下式确定：
Figure GDA0003702251960000083
S33、对比式(4)和(5)，可以建立关系式，对于波长为λi(i＝1,...,n)的光波，其对应的折射率分布有
Figure GDA0003702251960000084
选择i＝1的波长λ1对应的折射率参数(Nd01,Nw01,Hd1,Hw1)作为待求参数，则其他波长对应的折射率参数(Nd0i,Nw0i,Hdi,Hwi)(i＝2,...,n)可以表示为
Figure GDA0003702251960000091
即
Figure GDA0003702251960000092
这样就可以将不同波长对应的折射率表示为波长为λ1对应的折射率分布：
Figure GDA0003702251960000093
S34、根据折射率高度分布剖面，计算不同波长下目标恒星的等效地心角
Figure GDA0003702251960000094
Figure GDA0003702251960000095
其中,n0i测站地面大气折射率,h0为地球半径，Eci为目标恒星的视在高度角，ni(h)＝1+Ni(h)为h处大气折射指数，hm为目标恒星的等效地心距。
S35、由等效地心角
Figure GDA0003702251960000096
计算目标恒星的测量高度角Emi
Figure GDA0003702251960000097
S36、建立以模式参数为未知量(Nd01,Nw01,Hd1,Hw1)的参数待估方程组：
Egi-fi(Nd01,Nw01,Hd1,Hw1)＝0,(i＝1,...,n) (12)
其中，由于待估参数数目为4，所以n≥4，fi(Nd01,Nw01,Hd1,Hw1)是Emi关于模式参数(Nd01,Nw01,Hd1,Hw1)的函数表达形式，其计算公式表示为：
Figure GDA0003702251960000101
其中，(Nd01,Nw01,Hd1,Hw1)即为待求的优化模式参数。
S4、以模式参数的世界各地平均值(Nd0r,Nw0r,Hdr,Hwr)作为参考值，对大气折射高度分布模式参数待估方程组进行求解，在求解误差小于设定值ε时，即认为获得了有效的模式参数(Nd01,Nw01,Hd1,Hw1)。
S5、将根据模式参数的计算值(Nd01,Nw01,Hd1,Hw1)，代入式(3)得到波长λ1对应的大气折射率高度分布的Hopfield模型，其它波长对应的折射率则可以根据式(9)得出。
综上可知，本发明具有成本低、精度高、适应性好等优势，能够获得不同地域、气候、季节的大气折射率高度分布剖面，可以作为一种获得大气折射率高度分布剖面常规方法。
以上实施方式仅用于说明本发明，而并非对本发明的限制，有关技术领域的普通技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型，因此所有等同的技术方案也属于本发明的范畴，本发明的专利保护范围应由权利要求限定。

Claims (1)

1.一种基于多波长测量获得大气折射率高度分布剖面的方法，其特征在于，所述方法包括如下步骤：
S1、使用光电外测设备在不同波长下跟踪一颗预定恒星，输出并记录外测设备测量系的测角数据，同时获取相同时间段恒星地固系的精确星历数据；
S2、将恒星地固系的精确星历数据转换为外测设备测量系的理论测元数据；
S3、以大气折射率高度分布的霍普菲尔德模型的模式参数为待估未知数，利用高斯分层积分法计算外测设备的测角误差，从而建立模式参数待估方程组；
S4、对大气折射高度分布模式参数待估方程组进行求解，获得模式参数的数值解；
S5、利用数值求解的模式参数代入大气折射率高度分布的霍普菲尔德模型，获得大气折射率高度分布剖面数学模型；
所述步骤S2具体操作如下：
S21、将恒星地固系的精确星历数据(xei,yei,zei)转换至外测设备测量系，得到精确星历直角坐标的数据其计算公式为：
Figure FDA0003702251950000011
其中
Figure FDA0003702251950000012
Figure FDA0003702251950000013
为测站的天文经纬度；
S22、将恒星的精确星历直角坐标的数据(xci,yci,zci)从外测设备测量系的直角坐标系转换至极坐标系，得到Egi
Figure FDA0003702251950000025
其中，Egi为目标恒星的真实高度角；
所述步骤S3具体操作如下：
S31、大气折射率高度分布的霍普菲尔德模型为：
N(h)＝Nd(h)+Nw(h)
Figure FDA0003702251950000021
Nj(h)＝0，j＝d，w，h＞Hj (3)
其中，h为高度，hs为折射率分布起始高度，Nd(h)和Nw(h)分别为高度h处干项和湿项折射率，且对于世界各地平均剖面而言，Hd＝40136+148.72t0，单位为米，Hw＝11000米，t0为地表温度；
S32、在可见光与近红外波段的大气折射率表示为：
Figure FDA0003702251950000022
其中，P(h)为气压，T(h)为绝对温度，e(h)为水汽压，m1(λ)和m2(λ)分别为波长λ(μm)的函数，
由下式确定：
Figure FDA0003702251950000023
S33、对比式(4)和(5)，建立关系式，对于波长为λi(i＝1,...,n)的光波，其对应的折射率分布有：
Figure FDA0003702251950000024
选择i＝1的波长λ1对应的折射率参数(Nd01,Nw01,Hd1,Hw1)作为待求参数，则其他波长对应的折射率参数(Nd0i,Nw0i,Hdi,Hwi)(i＝2,...,n)表示为：
Figure FDA0003702251950000031
即
Figure FDA0003702251950000032
将不同波长对应的折射率表示为(Nd01,Nw01,Hd1,Hw1)的函数：
Figure FDA0003702251950000033
S34、根据折射率高度分布剖面，计算不同波长下目标恒星的等效地心角
Figure FDA0003702251950000034
Figure FDA0003702251950000035
其中，n0i测站地面大气折射率,h0为地球半径，Eci为目标恒星的视在高度角，ni(h)＝1+Ni(h)为h处大气折射指数，hm为目标恒星的等效地心距；
S35、由等效地心角
Figure FDA0003702251950000036
计算目标恒星的测量高度角Emi
Figure FDA0003702251950000037
S36、建立以模式参数为未知量(Nd01,Nw01,Hd1,Hw1)的参数待估方程组：
Egi-fi(Nd01,Nw01,Hd1,Hw1)＝0,(i＝1,...,n) (12)
其中，由于待估参数数目为4，所以n≥4，fi(Nd01,Nw01,Hd1,Hw1)是Emi关于模式参数(Nd01,Nw01,Hd1,Hw1)的函数表达形式，其计算公式表示为：
Figure DEST_PATH_FDA0003658403350000042
其中，(Nd01,Nw01,Hd1,Hw1)即为待求的优化模式参数；
所述步骤S4具体操作如下：
以模式参数的世界各地平均值(Nd0r,Nw0r,Hdr,Hwr)作为参考值，对式(12)的方程组进行数值求解，在求解误差小于设定值ε时，即认为获得了有效的模式参数(Nd01,Nw01,Hd1,Hw1)；
所述步骤S5具体操作如下：
根据模式参数的计算值(Nd01,Nw01,Hd1,Hw1)，代入式(3)则得到波长λ1对应的大气折射率高度分布的霍普菲尔德模型，其它波长对应的折射率则根据式(9)得出。
CN201911132808.4A 2019-11-19 2019-11-19 一种基于多波长测量获得大气折射率高度分布剖面的方法 Active CN111044489B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911132808.4A CN111044489B (zh) 2019-11-19 2019-11-19 一种基于多波长测量获得大气折射率高度分布剖面的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911132808.4A CN111044489B (zh) 2019-11-19 2019-11-19 一种基于多波长测量获得大气折射率高度分布剖面的方法

Publications (2)

Publication Number Publication Date
CN111044489A CN111044489A (zh) 2020-04-21
CN111044489B true CN111044489B (zh) 2022-09-23

Family

ID=70232611

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911132808.4A Active CN111044489B (zh) 2019-11-19 2019-11-19 一种基于多波长测量获得大气折射率高度分布剖面的方法

Country Status (1)

Country Link
CN (1) CN111044489B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113552544B (zh) * 2021-07-06 2023-09-26 中国人民解放军63796部队 一种面向实时的雷达电波折射积分迭代误差修正方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107219193A (zh) * 2017-05-08 2017-09-29 西安电子科技大学 大气折射率剖面的反演方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7265820B1 (en) * 2005-08-23 2007-09-04 G & A Technical Software, Inc. Method of determining atmospheric refraction profile using two spatially separated light sources
US7623065B2 (en) * 2006-05-26 2009-11-24 The Boeing Company Refractivity retrieval via direct measurement of GNSS bending angle
CN108898252A (zh) * 2018-07-02 2018-11-27 河南师范大学 一种全国对流层大气折射率剖面的预测方法
CN109900658B (zh) * 2019-03-06 2019-11-01 中国人民解放军国防科技大学 一种利用恒星观测数据校正光学成像卫星大气折射的方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107219193A (zh) * 2017-05-08 2017-09-29 西安电子科技大学 大气折射率剖面的反演方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
一种新的大气折射率剖面模型构建方法;李振锋等;《强激光与粒子束》;20151015(第10期);全文 *

Also Published As

Publication number Publication date
CN111044489A (zh) 2020-04-21

Similar Documents

Publication Publication Date Title
CN110031877B (zh) 一种基于grnn模型的区域nwp对流层延迟改正方法
CN103698305B (zh) 一种实时观测大气透射率的方法与系统
CN105844090B (zh) 多参数分步海面盐度反演方法以及装置
CN104181515A (zh) 一种基于蓝-黄波段高光谱数据的浅海水深反演方法
CN102609940A (zh) 利用地面激光扫描技术进行测量对象表面重建时点云配准误差处理方法
CN113408111B (zh) 大气可降水量反演方法及系统、电子设备和存储介质
CN114219994A (zh) 一种基于空海协同观测的海洋光学卫星辐射定标方法
CN116519913B (zh) 基于星载和地基平台融合的gnss-r数据土壤水分监测方法
CN115308386B (zh) 一种基于cygnss卫星数据的土壤盐分反演方法和系统
Wang et al. Atmospheric refraction calibration of geometric positioning for optical remote sensing satellite
CN108009392A (zh) 一种浓密植被地表的遥感反射率模型构建及标定应用方法
CN207456381U (zh) 提高激光跟踪仪测量精度的装置
CN111191380B (zh) 一种基于地基光谱仪测量数据的大气气溶胶光学厚度估算方法和装置
CN111044489B (zh) 一种基于多波长测量获得大气折射率高度分布剖面的方法
CN108318458B (zh) 一种适用于不同天气条件下的室外典型地物pBRDF的测量方法
CN115015842B (zh) 一种无测距经纬仪双站交会的大气折射修正方法
CN109764892B (zh) 近红外星等的校准方法
CN116295362A (zh) 一种基于全天时大气模式的恒星光线偏折计算方法
CN114252875B (zh) 一种成像高度计数据的高精度网格化方法
CN115855272A (zh) 一种生成近实时逐小时全天候地表温度的方法
CN111473868B (zh) 一种远距离红外测温方法
CN114663613A (zh) 一种基于双ftir遥测扫描成像系统的泄漏气云3-d重建方法
CN113049006A (zh) 一种基于模型参数估计的星光矢量蒙气差修正方法
CN111123406A (zh) 手持气象仪温度数据拟合方法
CN112161943B (zh) TanSat卫星XCO2反演数据的偏差校正方法和系统

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