CN105954764B - 一种基于椭球的gnss电离层层析投影矩阵获取方法 - Google Patents
一种基于椭球的gnss电离层层析投影矩阵获取方法 Download PDFInfo
- Publication number
- CN105954764B CN105954764B CN201610273936.0A CN201610273936A CN105954764B CN 105954764 B CN105954764 B CN 105954764B CN 201610273936 A CN201610273936 A CN 201610273936A CN 105954764 B CN105954764 B CN 105954764B
- Authority
- CN
- China
- Prior art keywords
- gnss
- inverting
- voxel
- ellipsoid
- waits
- 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
Links
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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种基于椭球的GNSS电离层层析投影矩阵获取方法,包括以下的步骤:S1:确定电离层层析待反演区域及待反演时段;S2:将待反演区域离散化为一系列三维像素;S3:准备电离层层析待反演区域及待反演时段内的GNSS卫星和GNSS接收机的大地空间直角坐标;S4:计算待反演时段内GNSS卫星到GNSS接收机的信号传播路径与三维像素各个面的交点;S5:计算三维像素内GNSS卫星信号的传播路径的截距长度;S6:构成投影矩阵。本发明通过基于椭球的GNSS电离层层析获取投影矩阵,克服了目前基于圆球提取投影矩阵方法时将地球椭球近似圆球的问题,保证了相对于地球椭球的投影矩阵能真实反映出GNSS信号传播路径在反演像素内的真实截距。
Description
技术领域
本发明涉及电离层层析投影矩阵获取方法,特别是涉及一种基于椭球的GNSS电离层层析投影矩阵获取方法。
背景技术
在高度约为60~1000km范围内的大气层中,由于太阳的紫外线、高能粒子和X射线等作用下,形成的存在着大量自由电子和正离子的区域,这个区域就是电离层。当电磁波在电离层中传播时会发生反射、折射、散射和吸收导致信号的路径产生弯曲、传播速度发生变化等现象。为了解释这些现象的产生需要对电离层的结构进行探测,电离层探测主要是为了获取表征电离层特征的参数(主要研究电子密度)。
随着全球卫星导航系统(Global Navigation Satellite System,GNSS)技术的发展,为电离层的研究带来新的技术和方法。利用GNSS观测数据对电离层探测,主要是获得沿卫星信号传播路径上的总电子含量(total electron content,TEC),通常TEC被认为是投影函数与相应垂直总电子含量(Vertical Total Electronic Content,VTEC)的函数。电离层TEC的空间上的分布只能反映出电离层某一层的水平结构,却难以描述其垂直结构。层次技术(Computerized Tomography,CT)在电离层探测中的应用形成了电离层层析成像技术(Computerized Ionospheric Tomography,CIT),该技术对电离层电子密度的时空变化进行分层,克服了单层模型的局限性,且适合于大尺度的电离层电子密度垂直分布及其扰动监测,实现了对电离层垂直结构的研究。目前,利用GNSS的电离层层析成像技术取得了很大的发展,但由于GNSS的电离层层析复杂性,电离层层析还有很多问题需要进一步的研究。
发明内容
发明目的:本发明的目的是提供一种有利于分析高分辨率电离层时空变化的基于椭球的GNSS电离层层析投影矩阵获取方法。
技术方案:为达到此目的,本发明采用以下技术方案:
本发明所述的基于椭球的GNSS电离层层析投影矩阵获取方法,包括以下的步骤:
S1:确定电离层层析待反演区域及待反演时段;
S2:将待反演区域离散化为一系列三维像素;
S3:准备电离层层析待反演区域及待反演时段内的GNSS卫星和GNSS接收机的大地空间直角坐标;
S4:计算待反演时段内GNSS卫星到GNSS接收机的信号传播路径与三维像素各个面的交点:首先确定折半搜索算法在三维像素各个面的目标位置,然后计算GNSS卫星到GNSS接收机的信号传播路径与三维像素各个面的交点;
S5:计算三维像素内GNSS卫星信号的传播路径的截距长度;
S6:构成投影矩阵。
进一步,所述步骤S1包括以下的步骤:
S1.1:确定电离层层析待反演区域;
确定电离层层析待反演区域的经度范围[L1 Ln]、纬度范围[B1 Bm]和椭球高范围[H1 Hq];
S1.2:确定电离层层析待反演时段[T1 Tt],电离层层析待反演的时刻为[T1 T2... Tt]。
进一步,所述步骤S2包括以下的步骤:
S2.1:根据电离层层析待反演区域的经度范围、纬度范围和椭球高范围,将待反演区域离散化为一系列的三维像素,即:
经度被离散化为[L1 L2 ... Ln];
纬度被离散化为[B1 B2 ... Bm];
椭球高被离散化为[H1 H2 ... Hq];
S2.2:对三维像素按照经度从西到东、纬度从南到北和椭球高从下到上的方式进行编号:当三维像素的左下角的经度为Li、纬度为Bj和椭球高为Hk,根据式(1)确定这个三维像素编号f:
f=(k-1)×((m-1)×(n-1))+(j-1)×(n-1)+i (1)
其中i=1,2,…,n、j=1,2,…,m、k=1,2,…,q。
进一步,所述步骤S4中确定折半搜索算法在三维像素各个面的目标位置包括以下的步骤:
S4.1:确定折半搜索算法在三维像素椭球高面的目标位置:
折半搜索算法在三维像素椭球高面的目标位置分别为S2.1中被间隔的椭球高[H1H2 ... Hq];
S4.2:判断待反演时段内GNSS卫星到GNSS接收机的信号传播路径是否属于电离层层析待反演区域;
根据折半搜索算法计算出Ts时刻GNSS信号传播路径与椭球高面H1、Hq的交点的大地空间直角坐标分别为[Xh1 Yh1 Zh1]T、[Xhq Yhq Zhq]T,将大地空间直角坐标分别转换成大地坐标[Lh1 Bh1 Hh1]T、[Lhq Bhq Hhq]T;根据式(2)判断Ts时刻GNSS卫星到GNSS接收机的信号传播路径是否属于电离层层析待反演区域:若满足式(2),则Ts时刻GNSS卫星到GNSS接收机的信号传播路径属于电离层层析待反演区域;若不满足式(2),则Ts时刻GNSS卫星到GNSS接收机的信号传播路径不属于电离层层析待反演区域;
其中s=1,2,…,t;
S4.3:确定折半搜索算法在三维像素经度面和纬度面的目标位置:
根据S4.2,若Ts时刻GNSS卫星到GNSS接收机的信号传播路径属于电离层层析待反演区域,将S2.1中满足式(3)的经度Li和纬度Bj分别作为折半搜索算法在三维像素经度面和纬度面的目标位置;
进一步,所述步骤S3包括以下的步骤:
S3.1:在电离层层析待反演时段[T1 Tt]内,准备待反演区域内以及周边的GNSS接收机的大地空间直角坐标;
S3.2:获取电离层层析待反演时段[T1 Tt]内GNSS卫星星历数据,并根据星历计算待反演时段内GNSS卫星不同时刻的大地空间直角坐标。
进一步,所述步骤S4中计算GNSS卫星到GNSS接收机的信号传播路径与三维像素各个面的交点包括以下的步骤:
S4.4:初始化折半搜索算法的参数:将Ts时刻GNSS卫星的大地直角坐标作为上边界将GNSS接收机的大地直角坐标作为下边界设置目标位置LBH,依据实际电离层层析精度需求设置阀值e,其中s=1,2,…,t;
S4.5:计算探测位置,并判定是否继续折半搜索:根据式(4)计算探测位置的大地空间直角坐标将探测位置的大地空间直角坐标转换成大地坐标若目标位置LBH与探测位置相应大地坐标值P.LBH的差的绝对值|P.LBH-LBH|<e,进行步骤S4.7;若|P.LBH-LBH|≥e,进行步骤S4.6;
S4.6:更新上下边界条件:若P.LBH-LBH<0,则上边界保持不变,将探测位置的大地空间直角坐标作为下边界;若P.LBH-LBH>0,则下边界保持不变,将探测位置的大地空间直角坐标作为上边界,然后回到步骤S4.5;
S4.7:停止折半搜索,则当前的探测位置被认为是Ts时刻GNSS信号传播路径与目标位置的交点。
进一步,所述步骤S5包括以下的步骤:
S5.1:对Ts时刻GNSS卫星到GNSS接收机的信号传播路径与三维像素各个面的交点进行排序:依据交点的坐标大小进行降序或者升序排列;
S5.2:根据S5.1的排序结果,确定相邻交点所穿过的三维像素编号;
将相邻交点的大地空间直角坐标[Xfl Yfl Zfl]T、[Xfu Yfu Zfu]T分别转换成大地坐标[Lfl Bfl Hfl]T、[Lfu Bfu Hfu]T,计算相邻交点的中点的大地坐标[L B H]T,若L、B和H满足式(5),则所穿过的三维像素左下角的经度为Li、纬度为Bj和椭球高为Hk,根据式(1)确定这个三维像素编号f;
其中i=1,2,…,n-1、j=1,2,…,m-1、k=1,2,…,q-1。
S5.3:根据S5.1的排序结果,计算GNSS信号传播路径在三维像素f的截距:
相邻交点的大地空间坐标分别为[Xfl Yfl Zfl]T、[Xfu Yfu Zfu]T,根据式(6)计算GNSS信号传播路径在三维像素f的截距df;
进一步,所述步骤S6包括以下的步骤:
S6.1:构造Ts时刻的投影向量As:
若GNSS信号传播路径穿过三维像素f,则GNSS信号传播路径在三维像素f的截距为df,即As(f)=df;若GNSS信号传播路径没有穿过三维像素f,则As(f)=0,f=1,2,…,(n-1)·(m-1)·(q-1),s=1,2,…,t;
S6.2:构造在电离层层析待反演时段[T1 Tt]内待反演区域的投影矩阵A:
构造电离层层析待反演的时刻[T1 T2 ... Tt]相应的投影向量[A1 A2 ... At],则式(7)就是电离层层析待反演时段[T1 Tt]内待反演区域的投影矩阵A;
有益效果:本发明通过基于椭球的GNSS电离层层析获取投影矩阵,克服了目前基于圆球提取投影矩阵方法时将地球椭球近似圆球的问题,保证了相对于地球椭球的投影矩阵能真实反映出GNSS信号传播路径在反演像素内的真实截距。并且,本发明提出的基于椭球的GNSS电离层层析获取投影矩阵有利于分析高分辨率电离层时空变化。
附图说明
图1为本发明的流程图;
图2为本发明的电离层层析中GNSS信号传播路径穿过三维像素的二维示意图。
具体实施方式
下面结合具体实施方式对本发明的技术方案作进一步的介绍。
本发明公开了一种基于椭球的GNSS电离层层析投影矩阵获取方法,如图1所示,包括以下的步骤:
S1:确定电离层层析待反演区域及待反演时段;
步骤S1包括以下的步骤:
S1.1:确定电离层层析待反演区域;
确定电离层层析待反演区域的经度范围[116.0° 123.0°];纬度范围[30.0°36.0°];椭球高范围[70km 1000km];
S1.2:确定电离层层析待反演时段[05:01:15UT 05:01:45UT],电离层层析待反演的时刻为[05:01:15UT 05:01:30UT 05:01:45UT]。
S2:将待反演区域离散化为一系列三维像素;
步骤S2包括以下步骤:
S2.1:根据电离层层析待反演区域的经度范围、纬度范围和椭球高范围,将待反演区域离散化为一系列的三维像素,即:
经度被离散化为(116.0+(i-1))°E,(i=1,2,…,8);
纬度被离散化为(30.0+(j-1))°N,(j=1,2,…,7);
椭球高被离散化为(70+30·(k-1))km,(k=1,2,…,32);
S2.2:对三维像素按照经度从西到东、纬度从南到北和椭球高从下到上的方式进行编号:当三维像素的左下角的经度为117.0°E、纬度为31.0°N和椭球高为100km,可知i=2,j=2,k=2,根据式(1)确定这个三维像素编号f为:f=(2-1)×((7-1)×(8-1))+(2-1)×(8-1)+2=51。
S3:准备电离层层析待反演区域及待反演时段内的GNSS卫星和GNSS接收机的大地空间直角坐标;
步骤S3包括以下步骤:
S3.1:在电离层层析待反演时段[05:01:15UT 05:01:45UT]内,准备待反演区域GNSS接收机的坐标如表1;
表1 GNSS接收机坐标
S3.2:获取电离层层析待反演时段[05:01:15UT 05:01:45UT]内GNSS卫星不同时刻的大地空间直角坐标如表2。
表2 卫星不同时刻的大地空间直角坐标
S4:计算待反演时段内GNSS卫星到GNSS接收机的信号传播路径与三维像素各个面的交点:根据图2所示的GNSS信号传播路径穿过三维像素的二维示意图,首先确定折半搜索算法在三维像素各个面的目标位置,然后计算GNSS卫星到GNSS接收机的信号传播路径与三维像素各个面的交点。
步骤S4中确定折半搜索算法在三维像素各个面的目标位置包括以下的步骤:
S4.1:确定折半搜索算法在三维像素椭球高面的目标位置:
折半搜索算法在三维像素椭球高面的目标位置分别为S2.1中被间隔的椭球高(70+30·(k-1))km,(k=1,2,…,32);
S4.2:判断待反演时段内GNSS卫星到GNSS接收机的信号传播路径是否属于电离层层析待反演区域;
根据折半搜索算法计算出Ts时刻GNSS信号传播路径与椭球高面70km、1000km的交点的大地空间直角坐标如表3所示;
表3 GNSS信号传播路径与椭球高面70km、1000km的交点大地空间直角坐标和其大地坐标
根据式(2)判断Ts时刻GNSS卫星到GNSS接收机的信号传播路径是否属于电离层层析待反演区域;
S4.3:确定折半搜索算法在三维像素经度面和纬度面的目标位置:
将在S2.1中满足式(3)的被间隔经度和纬度分别作为折半搜索算法在三维像素经度面和纬度面的目标位置。
5:01:15UT时刻:经度面目标位置为空;纬度面目标位置为[33.0°N 34.0°N 35.0°N];
5:01:30UT时刻:经度面目标位置为空;纬度面目标位置为[33.0°N 34.0°N 35.0°N];
5:01:45UT时刻:经度面目标位置为空;纬度面目标位置为[33.0°N 34.0°N 35.0°N];
S4.4:初始化折半搜索算法的参数:将Ts时刻GNSS卫星的大地直角坐标作为上边界,将GNSS接收机的大地直角坐标作为下边界,设置目标位置LBH,依据实际电离层层析精度需求设置阀值e;
其中,经度面和纬度面的阀值e设置为(5.7296×10-8)°;椭球高度面的阀值e设置为0.001m;
S4.5:计算探测位置,并判定是否继续折半搜索:根据式(4)计算探测位置的大地空间直角坐标将探测位置的大地空间直角坐标转换成大地坐标若目标位置LBH与探测位置相应大地坐标值P.LBH的差的绝对值|P.LBH-LBH|<e,进行步骤S4.7;若|P.LBH-LBH|≥e,进行步骤S4.6;
S4.6:更新上下边界条件:若P.LBH-LBH<0,则上边界保持不变,将探测位置的大地空间直角坐标作为下边界;若P.LBH-LBH>0,则下边界保持不变,将探测位置的大地空间直角坐标作为上边界,然后回到步骤S4.5;
S4.7:停止折半搜索,则当前的探测位置被认为是Ts时刻GNSS信号传播路径与目标位置的交点。
S5:计算三维像素内GNSS卫星信号的传播路径的截距长度;
步骤S5包括以下步骤:
S5.1:对GNSS卫星到GNSS接收机的信号传播路径与三维像素各个面交点排序:5:01:15UT时刻依据交点的X轴坐标大小进行升序排列结果如表4。
S5.2:根据S5.1的排序结果,确定相邻交点所穿过的三维像素编号。5:01:15UT时刻的相邻交点所穿过的三维像素编号如表4;
S5.3:根据S5.1的排序结果,计算GNSS信号传播路径在三维像素f的截距。5:01:15UT时刻的GNSS信号传播路径在三维像素f的截距如表4;
表4 (5:01:15)UT时刻交点坐标、相邻交点的像素编号及像素的截距
S6:构成投影矩阵;
步骤S6包括以下步骤:
S6.1:构造Ts时刻的投影向量As:
根据S2.1,Ts时刻的投影向量As为1×1302的行向量。若GNSS信号传播路径穿过在三维像素f,根据S5.3可知在三维像素f的截距为df,则As(f)=df;若GNSS信号传播路径没有穿过在三维像素f,则As(f)=0,(f=1,2,…,1302);
S6.2:构造在电离层层析待反演时段[05:01:15UT 05:01:45UT]内待反演区域的投影矩阵A:
根据S6.1构造电离层层析待反演的时刻[05:01:15UT 05:01:30UT 05:01:45UT]相应的投影向量[A1 A2 A3],则电离层层析待反演时段[05:01:15UT 05:01:45UT]内待反演区域的投影矩阵为:
因此,A就是在待反演区域的经度范围[116.0° 123.0°]、纬度范围[30.0°36.0°]、椭球高范围[70km 1000km]和电离层层析待反演时段[05:01:15UT 05:01:45UT]内的投影矩阵。通过基于椭球的GNSS电离层层析获取投影矩阵A,保证了GNSS信号传播路径在相对地球椭球的反演像素内真实截距,并能够快速定位截距所在像素的位置。本发明提出的基于椭球的GNSS电离层层析获取投影矩阵有利于分析高分辨率电离层时空变化。
Claims (7)
1.一种基于椭球的GNSS电离层层析投影矩阵获取方法,其特征在于:包括以下的步骤:
S1:确定电离层层析待反演区域及待反演时段;所述步骤S1包括以下的步骤:
S1.1:确定电离层层析待反演区域;
确定电离层层析待反演区域的经度范围[L1 Ln]、纬度范围[B1 Bm]和椭球高范围[H1Hq];
S1.2:确定电离层层析待反演时段[T1 Tt],电离层层析待反演的时刻为[T1 T2 ...Tt];
S2:将待反演区域离散化为一系列三维像素;
S3:准备电离层层析待反演区域及待反演时段内的GNSS卫星和GNSS接收机的大地空间直角坐标;
S4:计算待反演时段内GNSS卫星到GNSS接收机的信号传播路径与三维像素各个面的交点:首先确定折半搜索算法在三维像素各个面的目标位置,然后计算GNSS卫星到GNSS接收机的信号传播路径与三维像素各个面的交点;
S5:计算三维像素内GNSS卫星信号的传播路径的截距长度;
S6:构成投影矩阵。
2.根据权利要求1所述的基于椭球的GNSS电离层层析投影矩阵获取方法,其特征在于:所述步骤S2包括以下的步骤:
S2.1:根据电离层层析待反演区域的经度范围、纬度范围和椭球高范围,将待反演区域离散化为一系列的三维像素,即:
经度被离散化为[L1 L2 ... Ln];
纬度被离散化为[B1 B2 ... Bm];
椭球高被离散化为[H1 H2 ... Hq];
S2.2:对三维像素按照经度从西到东、纬度从南到北和椭球高从下到上的方式进行编号:当三维像素的左下角的经度为Li、纬度为Bj和椭球高为Hk,根据式(1)确定这个三维像素编号f:
f=(k-1)×((m-1)×(n-1))+(j-1)×(n-1)+i (1)
其中i=1,2,…,n、j=1,2,…,m、k=1,2,…,q。
3.根据权利要求2所述的基于椭球的GNSS电离层层析投影矩阵获取方法,其特征在于:所述步骤S4中确定折半搜索算法在三维像素各个面的目标位置包括以下的步骤:
S4.1:确定折半搜索算法在三维像素椭球高面的目标位置:
折半搜索算法在三维像素椭球高面的目标位置分别为S2.1中被间隔的椭球高[H1 H2... Hq];
S4.2:判断待反演时段内GNSS卫星到GNSS接收机的信号传播路径是否属于电离层层析待反演区域;
根据折半搜索算法计算出Ts时刻GNSS信号传播路径与椭球高面H1、Hq的交点的大地空间直角坐标分别为[Xh1 Yh1 Zh1]T、[Xhq Yhq Zhq]T,将大地空间直角坐标分别转换成大地坐标[Lh1 Bh1 Hh1]T、[Lhq Bhq Hhq]T;根据式(2)判断Ts时刻GNSS卫星到GNSS接收机的信号传播路径是否属于电离层层析待反演区域:若满足式(2),则Ts时刻GNSS卫星到GNSS接收机的信号传播路径属于电离层层析待反演区域;若不满足式(2),则Ts时刻GNSS卫星到GNSS接收机的信号传播路径不属于电离层层析待反演区域;
其中s=1,2,…,t;
S4.3:确定折半搜索算法在三维像素经度面和纬度面的目标位置:
根据S4.2,若Ts时刻GNSS卫星到GNSS接收机的信号传播路径属于电离层层析待反演区域,将S2.1中满足式(3)的经度Li和纬度Bj分别作为折半搜索算法在三维像素经度面和纬度面的目标位置;
4.根据权利要求1所述的基于椭球的GNSS电离层层析投影矩阵获取方法,其特征在于:所述步骤S3包括以下的步骤:
S3.1:在电离层层析待反演时段[T1 Tt]内,准备待反演区域内以及周边的GNSS接收机的大地空间直角坐标;
S3.2:获取电离层层析待反演时段[T1 Tt]内GNSS卫星星历数据,并根据星历计算待反演时段内GNSS卫星不同时刻的大地空间直角坐标。
5.根据权利要求1所述的基于椭球的GNSS电离层层析投影矩阵获取方法,其特征在于:所述步骤S4中计算GNSS卫星到GNSS接收机的信号传播路径与三维像素各个面的交点包括以下的步骤:
S4.4:初始化折半搜索算法的参数:将Ts时刻GNSS卫星的大地直角坐标作为上边界将GNSS接收机的大地直角坐标作为下边界设置目标位置LBH,依据实际电离层层析精度需求设置阀值e,其中s=1,2,…,t;
S4.5:计算探测位置,并判定是否继续折半搜索:根据式(4)计算探测位置的大地空间直角坐标将探测位置的大地空间直角坐标转换成大地坐标若目标位置LBH与探测位置相应大地坐标值P.LBH的差的绝对值|P.LBH-LBH|<e,进行步骤S4.7;若|P.LBH-LBH|≥e,进行步骤S4.6;
S4.6:更新上下边界条件:若P.LBH-LBH<0,则上边界保持不变,将探测位置的大地空间直角坐标作为下边界;若P.LBH-LBH>0,则下边界保持不变,将探测位置的大地空间直角坐标作为上边界,然后回到步骤S4.5;
S4.7:停止折半搜索,则当前的探测位置被认为是Ts时刻GNSS信号传播路径与三维像素各个面的交点。
6.根据权利要求2所述的基于椭球的GNSS电离层层析投影矩阵获取方法,其特征在于:所述步骤S5包括以下的步骤:
S5.1:对Ts时刻GNSS卫星到GNSS接收机的信号传播路径与三维像素各个面的交点进行排序:依据交点的坐标大小进行降序或者升序排列;
S5.2:根据S5.1的排序结果,确定相邻交点所穿过的三维像素编号;
将相邻交点的大地空间直角坐标[Xfl Yfl Zfl]T、[Xfu Yfu Zfu]T分别转换成大地坐标[LflBfl Hfl]T、[Lfu Bfu Hfu]T,计算相邻交点的中点的大地坐标[L B H]T,若L、B和H满足式(5),则所穿过的三维像素左下角的经度为Li、纬度为Bj和椭球高为Hk,根据式(1)确定这个三维像素编号f;
其中i=1,2,…,n-1、j=1,2,…,m-1、k=1,2,…,q-1;
S5.3:根据S5.1的排序结果,计算GNSS信号传播路径在三维像素f的截距:
相邻交点的大地空间坐标分别为[Xfl Yfl Zfl]T、[Xfu Yfu Zfu]T,根据式(6)计算GNSS信号传播路径在三维像素f的截距df;
7.根据权利要求6所述的基于椭球的GNSS电离层层析投影矩阵获取方法,其特征在于:所述步骤S6包括以下的步骤:
S6.1:构造Ts时刻的投影向量As:
若GNSS信号传播路径穿过三维像素f,则GNSS信号传播路径在三维像素f的截距为df,即As(f)=df;若GNSS信号传播路径没有穿过三维像素f,则As(f)=0,f=1,2,…,(n-1)·(m-1)·(q-1),s=1,2,…,t;
S6.2:构造在电离层层析待反演时段[T1 Tt]内待反演区域的投影矩阵A:
构造电离层层析待反演的时刻[T1 T2 ... Tt]相应的投影向量[A1 A2 ... At],则式(7)就是电离层层析待反演时段[T1 Tt]内待反演区域的投影矩阵A;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610273936.0A CN105954764B (zh) | 2016-04-27 | 2016-04-27 | 一种基于椭球的gnss电离层层析投影矩阵获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610273936.0A CN105954764B (zh) | 2016-04-27 | 2016-04-27 | 一种基于椭球的gnss电离层层析投影矩阵获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105954764A CN105954764A (zh) | 2016-09-21 |
CN105954764B true CN105954764B (zh) | 2018-10-23 |
Family
ID=56916819
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610273936.0A Active CN105954764B (zh) | 2016-04-27 | 2016-04-27 | 一种基于椭球的gnss电离层层析投影矩阵获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105954764B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106646607B (zh) * | 2016-12-22 | 2018-11-27 | 中国矿业大学 | 一种提高ct反演分辨率和效率的自适应不等间距网格划分方法 |
CN111273335B (zh) * | 2019-12-20 | 2021-09-17 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于垂测数据约束的电离层层析成像方法 |
CN111505702B (zh) * | 2020-06-15 | 2023-08-11 | 华东交通大学 | 一种基于垂直边界截断射线的电离层层析方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080232694A1 (en) * | 2007-03-21 | 2008-09-25 | Peter Sulatycke | Fast imaging data classification method and apparatus |
DE102009042327A1 (de) * | 2009-09-21 | 2011-04-07 | Siemens Aktiengesellschaft | Effiziente Visualisierung von Objekteigenschaften mittels Volume Rendering |
CN102930562A (zh) * | 2011-08-10 | 2013-02-13 | 中国科学院电子学研究所 | 电离层的压缩层析成像方法 |
CN103454695B (zh) * | 2013-08-20 | 2015-11-25 | 河海大学 | 一种gps电离层tec层析方法 |
CN103679806B (zh) * | 2013-12-19 | 2016-06-08 | 北京北科光大信息技术股份有限公司 | 自适应可视外壳生成方法及装置 |
CN104007479B (zh) * | 2014-06-13 | 2016-08-31 | 东南大学 | 一种基于多尺度剖分的电离层层析和电离层延迟改正方法 |
CN104933737B (zh) * | 2015-06-03 | 2016-04-27 | 北京航空航天大学 | 一种基于共轭梯度法的电离层层析成像混合反演方法 |
-
2016
- 2016-04-27 CN CN201610273936.0A patent/CN105954764B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105954764A (zh) | 2016-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Voxel-optimized regional water vapor tomography and comparison with radiosonde and numerical weather model | |
Kuang et al. | Quantifying the heat flux regulation of metropolitan land use/land cover components by coupling remote sensing modeling with in situ measurement | |
Champollion et al. | GPS water vapour tomography: preliminary results from the ESCOMPTE field experiment | |
Rohm et al. | Local tomography troposphere model over mountains area | |
CN104007479B (zh) | 一种基于多尺度剖分的电离层层析和电离层延迟改正方法 | |
Colin et al. | Aerodynamic roughness length estimation from very high-resolution imaging LIDAR observations over the Heihe basin in China | |
CN105954764B (zh) | 一种基于椭球的gnss电离层层析投影矩阵获取方法 | |
CN104155635B (zh) | 一种探地雷达单道电磁波谱三维定位方法 | |
CN103454695B (zh) | 一种gps电离层tec层析方法 | |
Wang et al. | Applications of spatial technology in schistosomiasis control programme in the People's Republic of China | |
Liu et al. | High-resolution mapping of mainland China’s urban floor area | |
Peng et al. | Modeling spatiotemporal patterns of understory light intensity using airborne laser scanner (LiDAR) | |
CN107783157A (zh) | 基于导航卫星的外辐射源智能选星方法、系统及雷达平台 | |
Yu et al. | Using a global reference sample set and a cropland map for area estimation in China | |
De Oliveira et al. | Assessment of radargrammetric DSMs from TerraSAR-X Stripmap images in a mountainous relief area of the Amazon region | |
Lee et al. | Creation of river terrain data using region growing method based on point cloud data from UAV photography | |
Zhao et al. | Combining ICESat-2 photons and Google Earth Satellite images for building height extraction | |
Shao et al. | Drivers of global surface urban heat islands: Surface property, climate background, and 2D/3D urban morphologies | |
Vather et al. | Cosmic ray neutrons provide an innovative technique for estimating intermediate scale soil moisture | |
Wang et al. | Comparison of changes in glacier area and thickness on the northern and southern slopes of Mt. Bogda, eastern Tianshan Mountains | |
Zhang et al. | A study of the interrelation between surface water and groundwater using isotopes and chlorofluorocarbons in Sanjiang plain, Northeast China | |
Liu et al. | Effects of land use and topography on spatial variety of soil organic carbon density in a hilly, subtropical catchment of China | |
Saberi et al. | Accuracy assessment and improvement of SRTM, ASTER, FABDEM, and MERIT DEMs by polynomial and optimization algorithm: A case study (Khuzestan Province, Iran) | |
CN108132140A (zh) | 一种星载光学遥感器全视场柱辐射率确定方法 | |
Zreda | Land-surface hydrology with cosmic-ray neutrons: Principles and applications |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |