一种无测距经纬仪双站交会的大气折射修正方法技术

本发明专利技术提供一种无测距经纬仪双站交会的大气折射修正方法，包括：根据飞行目标相对于测站中心在北天东坐标系中三个方向上的距离，确定飞行目标与测站站心距离；根据测站的坐标，获得测站与地心距离；通过交会解算获得飞行目标的坐标位置，并根据目标位置获取飞行目标与地心距离；计算地心角将地心角划分为多个且内的折射率不变；从测站位置转动至地心角根据大气分层折射率参数计算每转动的高度变化dh，并累加得到修正后的飞行目标与地心距离；根据修正后的飞行目标与地心距离，计算获得修正后的飞行目标到测站站心的距离以及俯仰角。本发明专利技术采用了更加接近于实际的“地心角”作为积分边界，有效地避免了“伪视在距离”引入的误差，精度更高。精度更高。精度更高。

【技术实现步骤摘要】
一种无测距经纬仪双站交会的大气折射修正方法


[0001]本专利技术涉及飞行目标探测
，具体涉及一种无测距经纬仪双站交会的大气折射修正方法。

技术介绍

[0002]大气为非均匀介质，光线在大气中传播会发生折射，从而使得光电经纬仪等飞行目标光学测量系统测角数据产生误差。大气折射误差是一项较大的系统误差；当飞行目标距测量站址较远、测角线较低时，大气折射对俯仰角的测量影响更大。
[0003]既有研究成果中对雷达研究较多，对电波折射修正的公式误差已很小，常用的以电波射线描迹法为代表的电波折射修正算法已具有很高的精度，采用经验参数和分层大气折射率模型结合地面气象参数的“盲”修正，也能满足大部分测试精度要求。在结合实时分层气象数据的折射修正中，其原理大多从测距信息出发，通过视在距离，结合测角数据进行运算的；这和雷达特别是相控阵雷达在测角方面的精度较低这一因素是分不开的。
[0004]《GJB2234A光电经纬仪事后处理》是目前光电经纬仪事后数据处理方面的标准；其修正方法即先获得视在距离，再通过反向递推方法获得俯仰角的修正。
[0005]专利“一种大气折射修正量计算快速迭代方法及装置”(公开号CN112098953 A)其核心思想是基于雷达测量的视在距离为核心，确定目标的近似地心距并据此确定电磁波传播距离，从而通过迭代确定距离修正量和仰角修正量。
[0006]专利“一种基于参数优化的大气折射误差实时修正方法”(公开号CN111025243 A)提供了一种基于参数优化的大气折射误差实时修正方法，包括：根据测站附近历史测量数据，选择一种合适的大气折射率高度分布模式，并确定其历史经验参数；外测设备跟踪恒星，输出并记录经纬仪测量系的测角数据，同时获取恒星地固系的精确星历数据，并转换为外测设备测量系的理论测元数据；建立模式参数优化方程组；以模式参数的历史经验值为参考值，对优化方程组进行数值求解，得到优化后的模式参数；利用优化后的模式参数对外测设备的测角数据进行实时修正。
[0007]专利“一种面向实时的雷达电波折射积分迭代误差修正方法”(公开号CN113552544 A)设计了一种可用于实时任务的雷达电波折射的误差修正方法，旨在提高实时雷达设备测元精度，通过地面采集的温湿压数据，利用Hopfield模型计算不同高度的对流层折射率，并拟合目标相对测站高度的初始值计算系数，再通过一次积分即可获取火箭视在距离及火箭真实高度，最后得到雷达测距、俯仰的大气折射修正值。为进一步提升算法计算效率，特别是火箭飞行高度较高时的积分计算速度，采用更易于计算机运行实现的高斯积分法，计算中使用高斯积分常数，达到提升积分效率的目的，本方法修正精度远高于现有实时经验修正模型。
[0008]专利“一种高精度大气折射率误差计算方法”(公开号CN109212516 A)公开了一种高精度大气折射率误差计算方法，在地面到离地1km的范围内，大气折射率呈线性衰减，当目标高度在此范围内时，利用线性插值法得到该目标的大气折射率；在离地面1km以上的范
围内，大气折射率成指数衰减，当目标高度在此范围内时，利用指数插值法得到该目标的大气折射率，并利用该大气折射率计算出目标的误差。
[0009]现有的大气折射修正方法均以测试设备能够获得“视在距离”为条件，即通过测距功能，能够测得站址与目标位置之间的距离。但是当采用无测距功能的经纬仪，并不能获得“视在距离”。若采用经纬仪测角原始数据进行交会解算得到目标位置，然后计算得到“伪视在距离”，其与实际距离相差较大，导致大气修正误差较大。

技术实现思路

[0010]为解决上述问题，本专利技术提供了如下的技术方案。
[0011]一种无测距经纬仪双站交会的大气折射修正方法，包括以下步骤：
[0012]确定飞行目标与测站站心距离；
[0013]根据测站的坐标，获得测站与地心距离；通过交会解算获得飞行目标的坐标位置，并根据飞行目标的坐标位置获取飞行目标与地心距离；
[0014]根据测站与地心距离、飞行目标与地心距离以及飞行目标与测站站心距离，计算地心角
[0015]将地心角划分为多个且每一个内的折射率不变；从测站位置转动，根据大气分层折射率参数计算每转动的高度变化dh，并累加dh得到修正后的飞行目标与地心距离；
[0016]根据修正后的飞行目标与地心距离，计算获得修正后的飞行目标到测站站心的距离和俯仰角。
[0017]优选地，所述测站与地心距离、飞行目标与地心距离以及飞行目标与测站站心距离的计算，包括以下步骤：
[0018]假设O(x,y,h)为发射点的大地坐标，P(x,y,h)为测站的大地坐标，M(x,y,h)为飞行目标在发射坐标系下的坐标；
[0019]确定飞行目标相对于测站中心在北天东坐标系中三个方向上的距离，X
pm
、Y
pm
、H
pm
：
[0020][0021]计算测站与地心距离R
op
、飞行目标与地心距离R
om
、飞行目标与测站站心距离R
pm
：
[0022][0023]R
op
＝R0+P(h)
[0024]R
om
＝R0+P(h)+M(h)
[0025]其中，R0表示地球半径。
[0026]优选地，所述地心角的计算，包括以下步骤：
[0027]根据余弦定理：
[0028][0029]优选地，所述飞行目标与地心距离的计算，包括以下步骤：
[0030]根据正弦定理a/sina＝b/sinb＝c/sinc得：
[0031][0032]其中，I
i+1
为下一薄层的入射角，表示为上式表示为：
[0033][0034]根据折射定律得：n
r
(i)
×
sin(I
i+1
)＝n
r
(i+1)
×
sin(π/2
‑
E
i+1
)，即：E
i+1
＝π/2
‑
arcsin[(n
r
(i)
×
sin(I
i+1
))/n
r
(i+1)]；式中，n
r
为折射率；
[0035]由于是从测站点开始转动，上式中R
om(i)
即为测站与地心距离R
op
，根据上式求取修正后的飞行目标与地心距离R
omn
。
[0036]优选地，所述修正后的飞行目标到测站站心的距离以及俯仰角的计算，包括以下步骤：
[0037]计算修正后的飞行目标到测站站心的距离R
pmn
：
[0038]根据余弦定理c2＝a2+b2‑2×
a
×
b
×
cosC得：
[0本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种无测距经纬仪双站交会的大气折射修正方法，其特征在于，包括以下步骤：确定飞行目标与测站站心距离；根据测站的坐标，获得测站与地心距离；通过交会解算获得飞行目标的坐标位置，并根据飞行目标的坐标位置获取飞行目标与地心距离；根据测站与地心距离、飞行目标与地心距离以及飞行目标与测站站心距离，计算地心角将地心角划分为多个且每一个内的折射率不变；从测站位置转动，根据大气分层折射率参数计算每转动的高度变化dh，并累加dh得到修正后的飞行目标与地心距离；根据修正后的飞行目标与地心距离，计算获得修正后的飞行目标到测站站心的距离和俯仰角。2.根据权利要求1所述的无测距经纬仪双站交会的大气折射修正方法，其特征在于，所述测站与地心距离、飞行目标与地心距离以及飞行目标与测站站心距离的计算，包括以下步骤：假设O(x,y,h)为发射点的大地坐标，P(x,y,h)为测站的大地坐标，M(x,y,h)为飞行目标在发射坐标系下的坐标；确定飞行目标相对于测站中心在北天东坐标系中三个方向上的距离，X
pm
、Y
pm
、H
pm
：计算测站与地心距离R
op
、飞行目标与地心距离R
om
、飞行目标与测站站心距离R
pm
：R
op
＝R0+P(h)R
om
＝R0+P(h)+M(h)其中，R0表示地球半径。3.根据权利要求2所述的无测距经纬仪双站交会的大气折射修正方法，其特征在于，所述地心角的计算，包括以下步骤：根据余弦定理：4.根据权利要求3所述的无测距经纬仪双站交会的大气折射修正...

【专利技术属性】
技术研发人员：胡杰，雷斌，吴建敏，李娜，王若晗，陈娟，唐舒，王珊珊，
申请(专利权)人：西安工业大学，
类型：发明
国别省市：