CN115047503A - 陆面星载gnss反射信号镜面反射点校准方法及*** - Google Patents
陆面星载gnss反射信号镜面反射点校准方法及*** Download PDFInfo
- Publication number
- CN115047503A CN115047503A CN202210673099.6A CN202210673099A CN115047503A CN 115047503 A CN115047503 A CN 115047503A CN 202210673099 A CN202210673099 A CN 202210673099A CN 115047503 A CN115047503 A CN 115047503A
- Authority
- CN
- China
- Prior art keywords
- reflection point
- current
- specular reflection
- determining
- initial specular
- 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.)
- Pending
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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining 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/40—Correcting position, velocity or attitude
-
- 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)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Aerials With Secondary Devices (AREA)
Abstract
本发明提供了一种陆面星载GNSS反射信号镜面反射点校准方法及***,涉及镜面反射点校准技术领域,包括获取初始镜面反射点集合;初始镜面反射点集合是全球导航卫星***将待测陆地区域等效为圆球面模型或椭球面模型确定的;获取任一初始镜面反射点为当前初始镜面反射点;根据当前初始镜面反射点的高程数据确定当前近似反射面;根据当前近似反射面确定当前初始镜面反射点的判决变量;根据判决变量校准当前初始镜面反射点。本发明通过反射点的高程数据构建当前近似反射面,并基于当前近似反射面校正将待测陆地区域等效为圆球面模型或椭球面模型确定镜面反射点,校正后的镜面反射点适用于陆面星载GNSS,提高了陆面星载GNSS的测量精度。
Description
技术领域
本发明涉及镜面反射点校准技术领域,特别是涉及一种陆面星载GNSS反射信号镜面反射点校准方法及***。
背景技术
星载GNSS-R(全球导航卫星***反射测量)几何关系是GNSS-R分析时空和信号特性的基础,主要描述了GNSS-R的空间体系分别与GNSS卫星、镜面反射点以及低轨卫星三者在空间体系中的绝对关系和相对关系。星载GNSS-R是一种镜向探测技术,不管是利用窄波束天线指向镜面反射点接收反射信号,还是选择在天线覆盖范围内对镜面反射信号进行处理,均需基于对GNSS-R镜面反射点的估计进行指向角计算、信号时延选取、以及多普勒范围内信号捕获处理,因此镜面反射点是星载GNSS-R几何关系的重要参考点。
在星载场景下必须考虑地球曲率的影响:反射面为一个曲面,而现有星载GNSS-R任务的设计初始目标均为海洋探测,而海洋可建模为符合某种概率分布的随机粗糙表面,所以常用坐标系(如WGS-84坐标系)将地球模型等同为椭球模型能够很好的适用海洋遥感。但是对于陆地表面而言,地形地貌复杂且不符合随机分布特性,使用椭球模型计算出的镜面反射点与实际镜面反射点的偏差较大:目前星载GNSS-R任务根据椭球模型计算得到的镜面反射路径、根据椭球模型计算得到的原始镜面反射点、实际真实镜面反射路径和实际镜面反射点如图1所示。由图1可看出:真实镜面反射点和根据椭球模型计算的原始镜面反射点有较大出入。
现有的镜面反射点搜索方法包括:黄金分割法、二分法、Gleason算法、Wu算法、Wagner算法、椭球面算法以及Adagrad自适应算法等,这些算法均能够实现对圆球面或椭球面上镜面反射点的位置估计。但陆地的实际反射面不是光滑的圆球面和椭球面,陆面的地势起伏在计算中需要利用高程模型对真实的地形进行模拟,而地面高程数据产品均具有有限的分辨率,是不连续的函数,无法实现镜面反射点在实际陆面反射面上的位置估计。
发明内容
本发明的目的是提供一种陆面星载GNSS反射信号镜面反射点校准方法及***,能够对利用高程数据对基于圆球面模型或椭球面模型确定的镜面反射点进行校正,提高陆面星载GNSS的测量精度。
为实现上述目的,本发明提供了如下方案:
一种陆面星载GNSS反射信号镜面反射点校准方法,包括:
获取初始镜面反射点集合;所述初始镜面反射点集合是由全球导航卫星***将待测陆地区域等效为圆球面模型或椭球面模型确定的;
获取任一初始镜面反射点为当前初始镜面反射点;
根据当前初始镜面反射点的高程数据确定当前近似反射面;
根据所述当前近似反射面确定所述当前初始镜面反射点的判决变量;
根据所述判决变量校准所述当前初始镜面反射点。
可选的,所述根据当前初始镜面反射点的高程数据确定当前近似反射面,包括:
以所述当前初始镜面反射点为中心,以预设距离为边长构建当前九宫格;
确定当前九宫格中所有边缘格对应区域的平均高程;所述边缘格为九宫格中除当前初始镜面反射点所在格之外的格子;
确定所述平均高程所在平面为当前近似反射面。
可选的,所述根据所述当前近似反射面确定所述当前初始镜面反射点的判决变量,包括:
获取入射信号的入射角θi;
获取出射信号的出射角θr;
确定入射信号在当前近似反射面上的投影与正东方向的第一夹角φi;
确定出射信号在当前近似反射面上的投影与正东方向的第二夹角φr;
根据所述入射角、所述出射角、所述第一夹角和所述第二夹角,根据公式ψ=|θi-θr|+r|φi-φr|确定当前初始镜面反射点的判决变量;其中,ψ为判决变量;r为惩罚系数;0<r<1。
可选的,所述根据所述判决变量校准所述当前初始镜面反射点,包括:
判断所述判决变量是否小于判决阈值,得到第一判断结果;
若所述第一判断结果为是,则确定当前初始镜面反射点为陆面镜面反射点;
若所述第一判断结果为否,则根据当前九宫格中所有边缘格校准所述当前初始镜面反射点。
可选的,所述根据当前九宫格中所有边缘格校准所述当前初始镜面反射点,包括:
确定当前九宫格中所有边缘格的中点为伪镜面反射点;
根据所述当前近似反射面,确定每个伪镜面反射点的判决变量为伪判决变量;
判断最小伪判决变量是否小于判决阈值,得到第二判断结果;
若所述第二判断结果为是,则确定最小伪判决变量对应的伪镜面反射点为陆面镜面反射点;
若所述第二判断结果为否,则确定当前初始镜面反射点为类漫反射点。
可选的,在所述根据所述判决变量校准所述当前初始镜面反射点之后,还包括:
当前初始镜面反射点并返回步骤“根据当前初始镜面反射点的高程数据确定当前近似反射面”直至历遍初始镜面反射点集合,得到陆面镜面反射点集合和类漫反射点集合。
一种陆面星载GNSS反射信号镜面反射点校准***,包括:
初始镜面反射点集获取模块,用于获取初始镜面反射点集合;所述初始镜面反射点集合是全球导航卫星***将待测陆地区域等效为圆球面模型或椭球面模型确定的;
当前初始镜面反射点获取模块,用于获取任一初始镜面反射点为当前初始镜面反射点;
当前近似反射面确定模块,用于根据当前初始镜面反射点的高程数据确定当前近似反射面;
判决变量确定模块,用于根据所述当前近似反射面确定所述当前初始镜面反射点的判决变量;
校准模块,用于根据所述判决变量校准所述当前初始镜面反射点。
可选的,所述当前近似反射面确定模块包括:
当前九宫格构建单元,用于以所述当前初始镜面反射点为中心,以预设距离为边长构建当前九宫格;
平均高程确定单元,用于确定当前九宫格中所有边缘格对应区域的平均高程;所述边缘格为九宫格中除当前初始镜面反射点所在格之外的格子;
当前近似反射面确定单元,用于确定所述平均高程所在平面为当前近似反射面。
可选的,所述判决变量确定模块包括:
入射角获取单元,用于获取入射信号的入射角θi;
出射角获取单元,用于获取出射信号的出射角θr;
第一夹角确定单元,用于确定入射信号在当前近似反射面上的投影与正东方向的第一夹角φi;
第二夹角确定单元,用于确定出射信号在当前近似反射面上的投影与正东方向的第二夹角φr;
判决变量确定单元,用于根据所述入射角、所述出射角、所述第一夹角和所述第二夹角,根据公式ψ=|θi-θr|+r|φi-φr|确定当前初始镜面反射点的判决变量;其中,ψ为判决变量;r为惩罚系数;0<r<1。
可选的,所述校准模块包括:
第一判断单元,用于判断所述判决变量是否小于判决阈值,得到第一判断结果;若所述第一判断结果为是,则调用陆面镜面反射点确定单元;若所述第一判断结果为否,则调用校准单元;
陆面镜面反射点确定单元,用于确定当前初始镜面反射点为陆面镜面反射点;
校准单元,用于根据当前九宫格中所有边缘格校准所述当前初始镜面反射点。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种陆面星载GNSS反射信号镜面反射点校准方法及***,包括获取初始镜面反射点集合;初始镜面反射点集合是全球导航卫星***将待测陆地区域等效为圆球面模型或椭球面模型确定的;获取任一初始镜面反射点为当前初始镜面反射点;根据当前初始镜面反射点的高程数据确定当前近似反射面;根据当前近似反射面确定当前初始镜面反射点的判决变量;根据判决变量校准当前初始镜面反射点。本发明通过反射点的高程数据构建当前近似反射面,并基于当前近似反射面校正将待测陆地区域等效为圆球面模型或椭球面模型确定镜面反射点,校正后的镜面反射点适用于陆面星载GNSS,提高了陆面星载GNSS的测量精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明现有技术中真实镜面反射点和根据椭球模型的原始镜面反射点的关系
图2为本发明实施例1中陆面星载GNSS反射信号镜面反射点校准方法流程图;
图3为本发明实施例2中地表区域待反演信息和观测量的关联图;
图4为本发明实施例3中第一校准原理图;
图5为本发明实施例3中第二校准原理图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种陆面星载GNSS反射信号镜面反射点校准方法及***,能够对利用高程数据对基于圆球面模型或椭球面模型确定的镜面反射点进行校正,提高陆面星载GNSS的测量精度。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1
如图2,本实施例提供了一种陆面星载GNSS反射信号镜面反射点校准方法,包括:
步骤101:获取初始镜面反射点集合;初始镜面反射点集合是全球导航卫星***将待测陆地区域等效为圆球面模型或椭球面模型确定的;
步骤102:获取任一初始镜面反射点为当前初始镜面反射点;
步骤103:根据当前初始镜面反射点的高程数据确定当前近似反射面;
步骤103包括:以当前初始镜面反射点为中心,以预设距离为边长构建当前九宫格;确定当前九宫格中所有边缘格对应区域的平均高程;边缘格为九宫格中除当前初始镜面反射点所在格之外的格子;确定平均高程所在平面为当前近似反射面。
步骤104:根据当前近似反射面确定当前初始镜面反射点的判决变量;
步骤104包括:获取入射信号的入射角θi;获取出射信号的出射角θr;确定入射信号在当前近似反射面上的投影与正东方向的第一夹角φi;确定出射信号在当前近似反射面上的投影与正东方向的第二夹角φr;根据入射角、出射角、第一夹角和第二夹角,根据公式ψ=|θi-θr|+r|φi-φr|确定当前初始镜面反射点的判决变量;其中,ψ为判决变量;r为惩罚系数;0<r<1。
步骤105:根据判决变量校准当前初始镜面反射点。
步骤105包括:判断判决变量是否小于判决阈值,得到第一判断结果。若第一判断结果为是,则确定当前初始镜面反射点为陆面镜面反射点。
若第一判断结果为否,则根据当前九宫格中所有边缘格校准当前初始镜面反射点。
根据当前九宫格中所有边缘格校准当前初始镜面反射点包括:确定当前九宫格中所有边缘格的中点为伪镜面反射点;根据当前近似反射面,确定每个伪镜面反射点的判决变量为伪判决变量;判断最小伪判决变量是否小于判决阈值,得到第二判断结果;若第二判断结果为是,则确定最小伪判决变量对应的伪镜面反射点为陆面镜面反射点;若第二判断结果为否,则确定当前初始镜面反射点为类漫反射点。
在步骤105之后,还包括:当前初始镜面反射点并返回步骤103直至历遍初始镜面反射点集合,得到陆面镜面反射点集合和类漫反射点集合。
实施例2
本实施例提供一种算法效率高,结果偏差小的校准方法,分别基于夹角和最短路径,利用GNSS-R几何关系和非线性规划方法,创建镜点搜索的迭代算法,从而实现对真实镜面反射点位置的估计。
步骤1:获取星载GNSS-R最重要的观测量:双基散射相关功率的时延-多普勒二维变量,即时延多普勒图(Delay Doppler Map,DDM),如公式所示:
此处,<|Y(τ,f)|2>表示GNSS多径反射信号中时延大小为τ、多普勒频率为f的反射信号分量与接收机本地伪随机码之间的相关功率,Ti为接收机的相干积分时间;λ为GNSS信号的载波波长;Pt是GNSS卫星发射功率;Gt为GNSS卫星发射天线增益;Gr(ρ)为GNSS反射信号接收天线增益;R0(ρ)是发射机到地面的范围R1(ρ)是地面到接收机的范围;为直射信号为p极化时反射信号采用q极化接收时的双基雷达散射系数;A表示信号照射区域。
步骤2:DDM包含了反射面的待反演信息。结合节中等延时区与等多普勒区,建立地表区域待反演信息和DDM特征参数(观测量)之间的联系。映射关系如图3。
步骤3:设星载GNSS-R DDM最大功率点的时延为τmax,多普勒频移为Dmax。则定义反射点S的时延为:
τs=τdir-δP (2)
其中τs为GNSS信号到达点s的码延时,τdir为直射信号的码延时,δP为直射反射信号的码延时:
由此可以定义初始镜面反射点S到功率极大值点的时延为:
δτ=τmax-τs (5)
步骤4:定义初始镜面反射点S的多普勒频移Ds:
Ds=DR+DT+Dclk (6)
式中,DT为GNSS卫星运动引起的多普勒频移,DR为GNSS-R卫星运动引起的频移,Dclk为卫星式中偏移引起的多普勒频移。DT和DR可以分别通过公式(7)和(8)计算:
式中,为GNSS-R接收机的速度,为GNSS卫星的速度,f为GNSS信号的载波频率,c为光速。GNSS-R卫星在当前近似反射面RST运动的单位向量,为GNSS卫星在当前近似反射面RST运动的单位向量,其中:
由此,反射点S相对于DDM功率极大值点的多普勒差值可以表示为:
δD=Dmax-DS (11)
步骤5:考虑地形影响,对地面做网格化划分,S1,S2,S3,S4分别为反射点S在北、南、东西四个方向的临接点(伪镜面反射点),如图4-5所示。
图4-5中,P1为GNSS信号的入射向量;P2为GNSS信号出射向量,P′1,P′2为二者在真实地表的投影。由于初始镜面反射点S的位置是根据椭圆模型推算,不满足斯涅尔定律(即不满足θi≠θr)。
将入射向量P1和出射向量P2投影到ENU坐标系中的反射平面:
P′1,E、P′1,N和P′1,U分别表示入射向量P1投影到E坐标轴、N坐标轴和坐标系轴后的入射向量。
P′2,E、P′2,N和P′2,U分别表示入射向量P2投影到E坐标轴、N坐标轴和坐标系轴后的入射向量。
进而求得θi,θr,φi和φr:
步骤6:基于θi,θr,φi和φr建立迭代算法,进行陆面镜面反射点的修正,描述如下:
算法第一步:
1.1)如图4-5,取数据集中第一个初始伪镜面反射点S点,利用DEM计算八个相邻单元平均高程,并取近似反射面。基于近似反射面计算法线、各种角度(包括θi,θr,φi和φr)。
1.2)定义判决变量ψ=|θi-θr|+r|φi-φr|,0<r<1为惩罚系数。设T0为判决阈值,若ψs>T0,计算S点八个邻域单元中心点的ψ值。
否则终止循环,确定S点为镜面反射点的最佳近似,将该点收入镜面反射集合。
1.3)取八个邻域的最小判决变量ψ′,若ψ′<ψs,将该点作为新的S点。否则终止循环,将该点数据收入非镜面反射集合。
1.4)取数据集中下一个伪镜面反射点S点,回到步骤1.1)。若数据集中所有伪镜面反射点S点均已运算完毕,结束算法,得到镜面反射集合和类漫反射集合,后续在反演土壤湿度或植被参数时,可对两类集合采取不同的算法,并确定不同的精度级别。
算法第二步:
2.2)设定能量阈值E0=ratio_k×E,(0<ratio_k<1),将DDM中的元素值Mi,j进行排序,计算出能量之和大于阈值的前m个元素值,记为[M0,M1,...,Mm],即其中M0为DDM元素的最大值,m为满足条件最小的值。其中参数ratio_k是一个大于0小于1的分数。
2.3)统计[M0,M1,...,Mm]出现的位置分布概率,M0位置所在的极大概率区间即为反射能量主区域。
最终得到镜面反射集合的校正后反射点以及类漫反射集合的反射能量主区域。反射能量主区域用于确定探测的主区域,进而建立反演模型。
实施例3
本实施例提供了一种陆面星载GNSS反射信号镜面反射点校准***,包括:
初始镜面反射点集获取模块,用于获取初始镜面反射点集合;初始镜面反射点集合是全球导航卫星***将待测陆地区域等效为圆球面模型或椭球面模型确定的。
当前初始镜面反射点获取模块,用于获取任一初始镜面反射点为当前初始镜面反射点。
当前近似反射面确定模块,用于根据当前初始镜面反射点的高程数据确定当前近似反射面。
判决变量确定模块,用于根据当前近似反射面确定当前初始镜面反射点的判决变量。
校准模块,用于根据判决变量校准当前初始镜面反射点。
其中,当前近似反射面确定模块包括:当前九宫格构建单元,用于以当前初始镜面反射点为中心,以预设距离为边长构建当前九宫格;平均高程确定单元,用于确定当前九宫格中所有边缘格对应区域的平均高程;边缘格为九宫格中除当前初始镜面反射点所在格之外的格子;当前近似反射面确定单元,用于确定平均高程所在平面为当前近似反射面。
具体的,判决变量确定模块包括:入射角获取单元,用于获取入射信号的入射角θi;出射角获取单元,用于获取出射信号的出射角θr;第一夹角确定单元,用于确定入射信号在当前近似反射面上的投影与正东方向的第一夹角φi;第二夹角确定单元,用于确定出射信号在当前近似反射面上的投影与正东方向的第二夹角φr;判决变量确定单元,用于根据入射角、出射角、第一夹角和第二夹角,根据公式ψ=|θi-θr|+r|φi-φr|确定当前初始镜面反射点的判决变量;其中,ψ为判决变量;r为惩罚系数;0<r<1。
此外,校准模块包括:第一判断单元,用于判断判决变量是否小于判决阈值,得到第一判断结果;若第一判断结果为是,则调用陆面镜面反射点确定单元;若第一判断结果为否,则调用校准单元;陆面镜面反射点确定单元,用于确定当前初始镜面反射点为陆面镜面反射点;校准单元,用于根据当前九宫格中所有边缘格校准当前初始镜面反射点。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的***而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种陆面星载GNSS反射信号镜面反射点校准方法,其特征在于,包括:
获取初始镜面反射点集合;所述初始镜面反射点集合是全球导航卫星***将待测陆地区域等效为圆球面模型或椭球面模型确定的;
获取任一初始镜面反射点为当前初始镜面反射点;
根据当前初始镜面反射点的高程数据确定当前近似反射面;
根据所述当前近似反射面确定所述当前初始镜面反射点的判决变量;
根据所述判决变量校准所述当前初始镜面反射点。
2.根据权利要求1所述的一种陆面星载GNSS反射信号镜面反射点校准方法,其特征在于,所述根据当前初始镜面反射点的高程数据确定当前近似反射面,包括:
以所述当前初始镜面反射点为中心,以预设距离为边长构建当前九宫格;
确定当前九宫格中所有边缘格对应区域的平均高程;所述边缘格为九宫格中除当前初始镜面反射点所在格之外的格子;
确定所述平均高程所在平面为当前近似反射面。
3.根据权利要求1所述的一种陆面星载GNSS反射信号镜面反射点校准方法,其特征在于,所述根据所述当前近似反射面确定所述当前初始镜面反射点的判决变量,包括:
获取入射信号的入射角θi;
获取出射信号的出射角θr;
确定入射信号在当前近似反射面上的投影与正东方向的第一夹角φi;
确定出射信号在当前近似反射面上的投影与正东方向的第二夹角φr;
根据所述入射角、所述出射角、所述第一夹角和所述第二夹角,根据公式ψ=|θi-θr|+r|φi-φr|确定当前初始镜面反射点的判决变量;其中,ψ为判决变量;r为惩罚系数;0<r<1。
4.根据权利要求2所述的一种陆面星载GNSS反射信号镜面反射点校准方法,其特征在于,所述根据所述判决变量校准所述当前初始镜面反射点,包括:
判断所述判决变量是否小于判决阈值,得到第一判断结果;
若所述第一判断结果为是,则确定当前初始镜面反射点为陆面镜面反射点;
若所述第一判断结果为否,则根据当前九宫格中所有边缘格校准所述当前初始镜面反射点。
5.根据权利要求4所述的一种陆面星载GNSS反射信号镜面反射点校准方法,其特征在于,所述根据当前九宫格中所有边缘格校准所述当前初始镜面反射点,包括:
确定当前九宫格中所有边缘格的中点为伪镜面反射点;
根据所述当前近似反射面,确定每个伪镜面反射点的判决变量为伪判决变量;
判断最小伪判决变量是否小于判决阈值,得到第二判断结果;
若所述第二判断结果为是,则确定最小伪判决变量对应的伪镜面反射点为陆面镜面反射点;
若所述第二判断结果为否,则确定当前初始镜面反射点为类漫反射点。
6.根据权利要求5所述的一种陆面星载GNSS反射信号镜面反射点校准方法,其特征在于,在所述根据所述判决变量校准所述当前初始镜面反射点之后,还包括:
当前初始镜面反射点并返回步骤“根据当前初始镜面反射点的高程数据确定当前近似反射面”直至历遍初始镜面反射点集合,得到陆面镜面反射点集合和类漫反射点集合。
7.一种陆面星载GNSS反射信号镜面反射点校准***,其特征在于,包括:
初始镜面反射点集获取模块,用于获取初始镜面反射点集合;所述初始镜面反射点集合是全球导航卫星***将待测陆地区域等效为圆球面模型或椭球面模型确定的;
当前初始镜面反射点获取模块,用于获取任一初始镜面反射点为当前初始镜面反射点;
当前近似反射面确定模块,用于根据当前初始镜面反射点的高程数据确定当前近似反射面;
判决变量确定模块,用于根据所述当前近似反射面确定所述当前初始镜面反射点的判决变量;
校准模块,用于根据所述判决变量校准所述当前初始镜面反射点。
8.根据权利要求7所述的一种陆面星载GNSS反射信号镜面反射点校准***,其特征在于,所述当前近似反射面确定模块包括:
当前九宫格构建单元,用于以所述当前初始镜面反射点为中心,以预设距离为边长构建当前九宫格;
平均高程确定单元,用于确定当前九宫格中所有边缘格对应区域的平均高程;所述边缘格为九宫格中除当前初始镜面反射点所在格之外的格子;
当前近似反射面确定单元,用于确定所述平均高程所在平面为当前近似反射面。
9.根据权利要求7所述的一种陆面星载GNSS反射信号镜面反射点校准***,其特征在于,所述判决变量确定模块包括:
入射角获取单元,用于获取入射信号的入射角θi;
出射角获取单元,用于获取出射信号的出射角θr;
第一夹角确定单元,用于确定入射信号在当前近似反射面上的投影与正东方向的第一夹角φi;
第二夹角确定单元,用于确定出射信号在当前近似反射面上的投影与正东方向的第二夹角φr;
判决变量确定单元,用于根据所述入射角、所述出射角、所述第一夹角和所述第二夹角,根据公式ψ=|θi-θr|+r|φi-φr|确定当前初始镜面反射点的判决变量;其中,ψ为判决变量;r为惩罚系数;0<r<1。
10.根据权利要求8所述的一种陆面星载GNSS反射信号镜面反射点校准***,其特征在于,所述校准模块包括:
第一判断单元,用于判断所述判决变量是否小于判决阈值,得到第一判断结果;若所述第一判断结果为是,则调用陆面镜面反射点确定单元;若所述第一判断结果为否,则调用校准单元;
陆面镜面反射点确定单元,用于确定当前初始镜面反射点为陆面镜面反射点;
校准单元,用于根据当前九宫格中所有边缘格校准所述当前初始镜面反射点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210673099.6A CN115047503A (zh) | 2022-06-14 | 2022-06-14 | 陆面星载gnss反射信号镜面反射点校准方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210673099.6A CN115047503A (zh) | 2022-06-14 | 2022-06-14 | 陆面星载gnss反射信号镜面反射点校准方法及*** |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115047503A true CN115047503A (zh) | 2022-09-13 |
Family
ID=83160658
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210673099.6A Pending CN115047503A (zh) | 2022-06-14 | 2022-06-14 | 陆面星载gnss反射信号镜面反射点校准方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115047503A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116540273A (zh) * | 2023-07-05 | 2023-08-04 | 天津大学 | Gnss-r镜面反射点初值的位置确定方法、装置及电子设备 |
CN117826197A (zh) * | 2023-12-29 | 2024-04-05 | 航天天目(重庆)卫星科技有限公司 | 一种偏振gnss-r接收机的极化校准方法 |
-
2022
- 2022-06-14 CN CN202210673099.6A patent/CN115047503A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116540273A (zh) * | 2023-07-05 | 2023-08-04 | 天津大学 | Gnss-r镜面反射点初值的位置确定方法、装置及电子设备 |
CN116540273B (zh) * | 2023-07-05 | 2023-09-19 | 天津大学 | Gnss-r镜面反射点初值的位置确定方法、装置及电子设备 |
CN117826197A (zh) * | 2023-12-29 | 2024-04-05 | 航天天目(重庆)卫星科技有限公司 | 一种偏振gnss-r接收机的极化校准方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110824510B (zh) | 一种提高gnss-r测高卫星接收海面反射信号数量的方法 | |
EP1836507B1 (fr) | Procede et dispositif de positionnement | |
Groves et al. | Intelligent urban positioning using multi-constellation GNSS with 3D mapping and NLOS signal detection | |
Gini et al. | Multibaseline cross-track SAR interferometry: A signal processing perspective | |
US6421003B1 (en) | Attitude determination using multiple baselines in a navigational positioning system | |
CN115047503A (zh) | 陆面星载gnss反射信号镜面反射点校准方法及*** | |
US9229111B2 (en) | Method for estimating the direction of arrival of navigation signals at a receiver after reflection by walls in a satellite positioning system | |
US20060132357A1 (en) | System and technique for calibrating radar arrays | |
WO2003001686A2 (en) | Method and system for calibrating location system | |
CN107783157B (zh) | 基于导航卫星的外辐射源智能选星方法、***及雷达平台 | |
RU2379709C1 (ru) | Способ определения скорости распространения и направления прихода ионосферного возмущения | |
CN112904384B (zh) | 基于机器学习的星载gnss-r高度要素探测装置及方法 | |
CN114879197A (zh) | 一种星上ddm实时定标的方法 | |
CN115184916A (zh) | 一种海面风速联合反演方法、装置、介质及计算设备 | |
Seube et al. | Multibeam echo sounders-IMU automatic boresight calibration on natural surfaces | |
CN114637004A (zh) | 星载gnss-s多站雷达在轨处理与舰船信息融合***及方法 | |
US11428802B2 (en) | Localization using particle filtering and image registration of radar against elevation datasets | |
CN114966690B (zh) | 一种gnss-r sar双星融合成像方法及*** | |
Zhang et al. | Research on sea surface wind speed FM based on CYGNSS and HY-2B microwave scatterometer | |
King | The GPS contribution to the error budget of surface elevations derived from airborne LIDAR | |
Pettersson et al. | Detection and imaging of moving targets in wide band SAS using fast time backprojection combined with space time processing | |
US7142149B2 (en) | Mensuration for the conformal range migration algorithm | |
RU2287169C2 (ru) | Способ радиоконтроля | |
RU2285934C2 (ru) | Способ однопозиционного определения местоположения дкмв передатчиков | |
Lin et al. | Computing specular points over complex land surfaces for airborne GNSS-r applications |
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 |