申请/专利权人：北京航空航天大学

申请日：2017-08-25

公开（公告）日：2020-08-04

公开（公告）号：CN107561534B

主分类号：G01S13/90(20060101)

分类号：G01S13/90(20060101)

优先权：

专利状态码：有效-授权

法律状态：2020.08.04#授权;2018.02.02#实质审查的生效;2018.01.09#公开

摘要：本发明公开了一种基于全极化高轨SAR的电离层时变TEC测量方法，包括步骤一：计算方位向合成孔径点数；步骤二：方位向数据补零；步骤三：方位向傅里叶变换；步骤四：方位向信号解压缩，获取等效方位向回波频谱信号；步骤五：方位向逆傅里叶变换；步骤六：估计法拉第旋转角；步骤七：获取合成孔径时间内的时变TEC；经过以上七个步骤，利用全极化高轨SAR信号完成了对合成孔径时间内的时变TEC精确测量。本发明具有时变TEC测量精度高的特点，由于采用法拉第旋转角模型与TEC反演模型，相对于传统测量方法，具有更高的测量精度。

主权项：1.一种基于全极化高轨SAR的电离层时变TEC测量方法，其特征在于，包括以下几个步骤：步骤一：计算方位向合成孔径点数；根据雷达系统参数，计算方位向合成孔径点数Num_a；步骤二：方位向数据补零；根据步骤一计算得到的方位向合成孔径点数Num_a，分别在每个极化通道图像数据每一列的头部和尾部补零，补零数目均为Num_a2，得到新的四个极化通道图像数据分别为MHH-ins、MHV-ins、MVH-ins、MVV-ins；步骤三：方位向傅里叶变换；将步骤二得到的补零后的四个极化通道图像数据MHH-ins、MHV-ins、MVH-ins、MVV-ins沿每个距离门按列进行快速傅里叶变换FFT，得到方位向频谱数据MHH-FFT、MHV-FFT、MVH-FFT、MVV-FFT；步骤四：方位向信号解压缩，获取等效方位向回波频谱信号；根据雷达系统参数，计算每个距离门对应的参考斜距Rref，进而计算出方位向信号解压缩因子Φdecom，利用该解压缩因子乘以步骤三所得的方位向频谱数据MHH-FFT、MHV-FFT、MVH-FFT、MVV-FFT得到等效方位向回波频谱信号MHH-decom-FFT、MHV-decom-FFT、MVH-decom-FFT、MVV-decom-FFT；步骤五：方位向逆傅里叶变换；将步骤四得到的等效方位向回波频谱信号MHH-decom-FFT、MHV-decom-FFT、MVH-decom-FFT、MVV-decom-FFT沿每个距离门进行快速逆傅里叶变换，得到方位向回波信号MHH-echo、MHV-echo、MVH-echo、MVV-echo；步骤六：估计法拉第旋转角；根据步骤五得到的四个极化通道方位向回波信号MHH-echo、MHV-echo、MVH-echo、MVV-echo，通过线性组合得到一组正交圆极化波信号Z12和Z21，通过该组正交圆极化波信号共轭相乘取幅角FRtemp，对FRtemp进行平滑处理后得到法拉第旋转角估计值FR；步骤七：获取合成孔径时间内的时变TEC；根据步骤六得到的估计法拉第旋转角值FR，利用时变TEC与法拉第旋转角间的线性转换关系，获得合成孔径时间内每个方位时刻的TEC，即合成孔径时间内的时变TEC。

全文数据：一种基于全极化高轨SAR的电离层时变TEC测量方法技术领域[0001]本发明涉及一种基于全极化高轨合成孔径雷达（SAR的电离层时变总电子量TEC测量方法，属于信号处理技术领域。背景技术[0002]近年来，随着科技发展与社会进步，卫星通信、卫星导航和天基雷达系统已广泛应用于军事和民用的各个方面，并成为人类生活必不可少的工具，这使得包括电离层在内的空间环境监测和技术保障愈显重要和迫切。同时，BI0MASS、NISAR等卫星计划将对全球包括电离层在内的大气环境进行探测，该探测数据将为地球空间气象科学以及空间信息系统气象保障等方面提供重要帮助，具有极为重大科学意义和应用价值。[0003]当前国内外针对电离层探测仍主要依靠传统电离层探测手段，即地面站垂测、基于GPS信号掩星反演等。但是，由于存在探测区域受地基站布设网络限制，需要搭载专门的电离层探测载荷等问题，导致无法实现全球电离层高精度测量。针对该问题，近年来国外学者提出了利用低频段SAR信号在电离层传播过程中会引入电离层效应误差，通过误差估计进而对电离层特征参数进行反演，实现电离层高精度探测。现有研究中，利用全极化SAR系统，可以通过估计电离层所引入的法拉第旋转角，根据法拉第旋转角与电离层TEC的映射关系，实现电离层空变TEC单位为TECU反演。然而，由于现有研究中所涉及的均为中低轨星载全极化SAR系统，卫星飞行速度较快且合成孔径时间较短，仅能针对不同照射地区的电离层TEC分布即空变TEC进行探测，而无法实现固定区域的时变TEC精确探测。[0004]同时，目前针对高轨SAR的研究仅仅停留在电离层所引起的成像误差补偿上，多数均采用相位梯度自聚焦的处理方法估计电离层时变TEC所引起的误差相位。但由于该方法处理过程中不采用固定的误差模型，仅基于像素间相位梯度变化估计误差相位，从而完全丢失掉电离层时变TEC所引起的绝对相位误差，因此无法实现时变TEC精确探测。发明内容[0005]本发明的目的是为了解决利用星载SAR系统实现电离层时变TEC测量的瓶颈技术问题，提出一种基于全极化高轨SAR的电离层时变TEC测量方法，充分利用高轨SAR飞行速度较慢且具有较长合成孔径时间，对于同一地区可实现较长时间持续观测，以及采用全极化信号可精确估计法拉第旋转角的特点，实现基于全极化高轨SAR的电离层时变TEC测量新体制。该方法能够实现全球实时电离层时变TEC精确测量，获取毫秒级时变TEC测量结果，扩展了SAR在遥感科学与地球空间气象科学中的应用。[0006]—种基于全极化高轨SAR的电离层时变TEC测量方法，包括以下几个步骤：[0007]步骤一:计算方位向合成孔径点数；[0008]根据雷达系统参数，计算方位向合成孔径点数Num_a;[0009]步骤二:方位向数据补零；[0010]根据步骤一所计算得到的方位向合成孔径点数Num_a，分别在每个极化通道图像数据每一列的头部和尾部补零，补零数目均为Num_a2,得到新的四个极化通道图像数据分力Ij为MhH-ins、MhV-ins、MvH-ins、Mw-ins;[0011]步骤三:方位向傅里叶变换；[0012]将步骤二得到的补零后的四个极化通道图像数据MHH-ins、MHV-ins、MvH-ins、MVV-ins沿每个距离门（按列）进行快速傅里叶变换FFT，得到方位向频谱数据Mhh-fft、Mhv-fft、Mvh-fft、Mvv-FFT；[0013]步骤四:方位向信号解压缩，获取方位向回波信号；[0014]根据雷达系统参数，计算每个距离门海列对应的参考斜距Rref，以此为参考计算出方位向信号解压缩因子，利用该解压缩因子乘以步骤三所得的方位向频谱数据M.FFT、MHV-FFT、MvH-FFT、Mvv-FFT得至丨J等效方位向回波步页谱信号Μ.decomFFT、MHV-decom-FFT、MvH-decom-FFT、Mw-decom-FFT;[0015]步骤五:方位向逆傅里叶变换；[0016]将步骤四得到的等效方位向回波频谱信号Mhh-decom—FFT、Mhv-decom—FFT、Mvh-decom—FFT、Mvv—de_-FFT沿每个距离门（按列）进行快速逆傅里叶变换（IFFT，得5IJ方位向回波信号MhH-echo、MhV-echo、MvH-echo、Mw-echo;[0017]步骤六:估计法拉第旋转角；[0018]根据步骤五得到的四个极化通道方位向回波信号Mffl-ech。、MHV-ech。、MvH-ech。、Mw-ech。，通过线性组合得到一组正交圆极化波信号Z12和Z21，通过该组正交圆极化波信号共辄相乘取幅角FRtemp，对FRtemp进行平滑处理后得到法拉第旋转角估计值FR;[0019]步骤七:获取合成孔径时间内时变TEC;[0020]根据步骤六得到的法拉第旋转角估计值FR，利用时变TEC与法拉第旋转角间的线性转换关系，获得合成孔径时间内每个方位时刻的TEC，即合成孔径时间内的时变TEC。[0021]本发明的优点在于：[0022]⑴本发明提出一种基于全极化高轨SAR的电离层时变TEC测量方法，该方法具有时变TEC测量精度高的特点。由于采用法拉第旋转角模型与TEC反演模型，相对于传统测量方法，具有更高的测量精度。[0023]2本发明提出一种基于全极化高轨SAR的电离层时变TEC测量方法，该方法具有可持续测量固定区域TEC时变特性的特点。由于采用高轨SAR体制，相比较传统测量方法，利用雷达合成孔径时间长的特性，可在分钟量级的合成孔径时间段内对固定区域进行时变TEC持续测量。[0024]⑶本发明提出一种基于全极化高轨SAR的电离层时变TEC测量方法，该方法具有可测量时变TEC短时间间隔时变特性的特点。由于采用高轨SAR体制，相比较传统测量方法，利用雷达脉冲重复间隔短的特性，可实现变化时间间隔为毫秒量级的时变TEC测量。[0025]⑷本发明提出一种基于全极化高轨SAR的电离层时变TEC测量方法，该方法具有实时观测的特点。由于采用全极化高轨SAR体制，相比较传统测量体制，可通过星上实时处理SAR系统获取的全极化数据得到时变TEC实时测量结果。[0026]5本发明提出一种基于全极化高轨SAR的电离层时变TEC测量方法，该方法具有全球覆盖观测的特点。由于采用全极化高轨SAR体制，相比较传统测量体制，不受地面站布设网络的限制，可利用高轨SAR全球覆盖观测的特点实现全球时变TEC测量。[0027]6本发明提出一种基于全极化高轨SAR的电离层时变TEC测量方法，该方法具有观测设备简单的特点。由于采用全极化高轨SAR体制，相比较传统测量体制，无需搭载专门用于电离层探测的探测设备如:无线电等离子成像仪、电离层顶层探测器等），可直接利用全极化高轨SAR获取的全极化数据进行时变TEC反演，极大地节省探测成本且便于卫星搭载。附图说明[0028]图1是本发明提出的一种基于全极化高轨SAR的电离层时变TEC测量方法的方法流程图；[0029]图2是本发明实施例中在仿真全极化数据中加入的合成孔径时间内的时变TEC;[0030]图3是本发明实施例中采用本发明提出方法估计得到的合成孔径时间内的时变TEC0具体实施方式[0031]下面将结合附图对本发明作进一步详细说明。[0032]本发明的一种基于全极化高轨SAR的电离层时变TEC测量方法，该方法基于全极化高轨SAR的四个极化通道数据二维复数矩阵分别表示为施1]、1^、1^、1^，设其大小均为1XNr，一维是方位向，有Na个采样点，PRF表示雷达的脉冲重复频率，1PRF称为雷达的脉冲重复周期，不同的方位向采样点对应不同的方位时刻，相邻两个方位向采样点相差时间间隔为1PRF;另一维是距离向，有Nr个采样点，表示每个方位时刻开启一次回波接收窗，对回波信号进行采样，一次连续距离向采样有Nr个采样点，采样率为fs。[0033]本发明是一种基于全极化高轨SAR的电离层时变TEC测量方法，方法流程图如图1所示，具体包括以下几个步骤：[0034]步骤一:计算方位向合成孔径点数；[0035]根据雷达系统参数，计算方位向合成孔径点数Num_a:[0037]式中，λ表示雷达系统工作波长，R。表示场景中心点的参考斜距，1^表示方位向天线长度，Vrrf表示以场景中心点为参考时的等效速度，PRF表示雷达脉冲重复频率。floorX表示取不大于X的最大偶数运算。[0038]步骤二:方位向数据补零；[0039]根据步骤一计算得到的方位向合成孔径点数Num_a，分别在每个极化通道图像数据每一列的头部和尾部补零，补零数目均为Num_a2,得到新的四个极化通道图像数据分别^MHH-insNMHV-insNMvH-insNMw-ins,：[0040]设原始四个极化通道数据!^、!^、!^、!^分别为：[0045]式中，表示HH通道图像中第（I，1个像素的数值，M☆表示HH通道图像中第I，Nr个像素的数值，Mg表示HH通道图像中第Na，1个像素的数值，表示HH通道图像中第Na，Nr个像素的数值;表示HV通道图像中第（1，1个像素的数值，表示HV通道图像中第（I，Nr个像素的数值，Μ。1表示HV通道图像中第Na，1个像素的数值，My。表示HV通道图像中第Na，Nr个像素的数值;表示VH通道图像中第1，1个像素的数值，M巧表示VH通道图像中第I，Nr个像素的数值，Mf表示VH通道图像中第Na，1个像素的数值，表示VH通道图像中第Na，Nr个像素的数值;表示VV通道图像中第（1，1个像素的数值，Μ”表示VV通道图像中第（I，Nr个像素的数值，表示VV通道图像中第Na，1个像素的数值，表示VV通道图像中第Na，Nr个像素的数值。[0046]结合步骤一计算得到的方位向合成孔径点数Num_a，分别在四个极化通道数据每一列的头部和尾部分别补零Num_a2,得到补零后新的四个极化通道图像数据MHH-ins、Mhv-ins、MvH-ins、Mw-ins，其均扩大为大*Na-insXNr-ins的二维复数矩阵，且Na-ins=Na+Num_a，Nr-ins=Nr，四个通道的信号分别为：[0051]式中，I表示补零后HH通道图像中第（Num_a2+l，1个像素的数值，表示补零后HH通道图像中第Num_a2+l，Nr个像素的数值表示补零后HH通道图像中第Num_a2+Na，1个像素的数值，.'表示补零后HH通道图像中第Num_a2+Na，Nr个像素的数值；1表示补零后HV通道图像中第Num_a2+1，1个像素的数值，表示补零后HV通道图像中第Num_a2+l，Nr个像素的数值，_表示补零后HV通道图像中第Num_a2+Na，l个像素的数值，表示补零后HV通道图像中第Num_a2+Na，Nr个像素的数值；表示补零后VH通道图像中第Num_a2+l，1个像素的数值，〃表示补零后VH通道图像中第Num_a2+1，Nr个像素的数值：表示补零后VH通道图像中第Num_a2+Na，1个像素的数值，表示补零后VH通道图像中第Num_a2+Na，Nr个像素的数值;α_表示补零后VV通道图像中第Num_a2+l，1个像素的数值，表示补零后VV通道图像中第Num_a2+l，Nr个像素的数值，11表示补零后VV通道图像中第Num_a2+Na，l个像素的数值，7v表示补零后VV通道图像中第Num_a2+Na，Nr个像素的数值。[0052]步骤三:方位向傅里叶变换；[0053]将步骤二得到的补零后的四个极化通道图像数据MHH-ins、MHV-ins、MvH-ins、Mw-insVS*·^个距离门（按列）进行快速傅里叶变换FFT，得到方位向频谱数据Mhh-fft、Mhv-fft、Mvh-fft、Mvv-FFT，具体表不为：[0058]式中，.表示HH通道方位向频谱数据中第1，1个像素的数值，表示HH通道方位向频谱数据中第（I，Nr—ins个像素的数值，1表示HH通道方位向频谱数据中第Na-ins，1个像素的数值，羡示HH通道方位向频谱数据中第Na-ins，Nr-ins个像素的数值;表示HV通道方位向频谱数据中第1，1个像素的数值，表示HV通道方位向频谱数据中第（I，Nr-ins个像素的数值，表示HV通道方位向频谱数据中第Na-ins，1个像素的数值，_表示HV通道方位向频谱数据中第Na-ins，Nr-ins个像素的数值;表示VH通道方位向频谱数据中第（1，1个像素的数值，表示VH通道方位向频谱数据中第（l，Nr-ins个像素的数值，表示VH通道方位向频谱数据中第Na-ins，1个像素的数值，表示VH通道方位向频谱数据中第Na-ins，Nr-ins个像素的数值;表示VV通道方位向频谱数据中第（I，1个像素的数值，i.表示VV通道方位向频谱数据中第（l，Nr-ins个像素的数值，，表示VV通道方位向频谱数据中第Na-ins，1个像素的数值：％表示VV通道方位向频谱数据中第Na-ins，Nr-ins个像素的数值。FFT·表示对二维矩阵的每列进行快速傅里叶变换。[0059]在本发明中，快速傅里叶变换FFT和快速逆傅里叶逆变换（IFFT的具体内容，请参考2012年6月电子工业出版社、IanGCumming等著、洪文等译的《合成孔径雷达成像-算法与实现》一书，第18页至第19页。[0060]步骤四:方位向信号解压缩，获取等效方位向回波频谱信号；[0061]根据雷达系统参数，计算每个距离门（每列对应的参考斜距Rref，以此为参考计算出方位向信号解压缩因子Φ^。。》，利用该解压缩因子乘以步骤三所得的方位向频谱数据M.FFT、MhV-FFT、MvH-FFT、Mw-FFT得至丨J等效方位向回波步页谱信号M.decomFFT、MhV-decom-FFT、MvH-decom-FFT、Mw-decom-FFT，本:[0062]I根据雷达系统参数，计算每个距离门（每列对应的参考斜距Rref:[0063][0064]式中，c表示光速。需要注意的是，式（14为每个距离门^^列对应的参考斜距Rref二维矩阵表示方式，式中处于同一距离门（每列不同方位向（每行）的参考斜距数值是相同的，即每个距离门（每列采用相同的参考斜距。[0065]2雷达的脉冲重复频率为PRF，计算方位向频谱数据每个方位向（每行的方位频率匕为：[0067]⑶根据上述计算所得每个距离门海列对应的参考斜距Rref以及每个方位向海行的方位频率fa，可计算得到方位向信号解压缩因子Φd_m:[0069]式中，Λη表示方位频率fa—维数组中的第1个频率值，表示方位频率fa—维数组中的第Na-ins个频率值，〇表示参考斜距Rref二维矩阵中第（1，1个斜距值，表示参考斜距Rref二维矩阵中第（I，Nr-ins个斜距值，乾^ίΒίΙ表示参考斜距Rref二维矩阵中第Na-ins，1个斜距值，」表示参考斜距Rref二维矩阵中第Na-ins，Nr-ins个斜距值。[0070]4将计算得到的方位向信号解压缩因子乘以方位向频谱数据Mhh-fft、MhV-FFT、MvH-FFT、Mw-FFT得至丨J等效方位向回波步页谱信号Mffll-decom-FFT、MhV-decom-FFT、MvH-decom-FFT、Mw-decom-FFT：[0075]式中，j表示虚数单位，SP表示方位向信号解压缩因子二维矩阵中第（1，1个因子数值，：表示方位向信号解压缩因子φ^。。》二维矩阵中第（1，Nr-ins个因子数值表示方位向信号解压缩因子Φ。™二维矩阵中第Na-ins，1个因子数值，表示方位向信号解压缩因子二维矩阵中第Na-ins，Nr-ins个因子数值。[0076]步骤五:方位向逆傅里叶变换；[0077]将步骤四得到的等效方位向回波频谱信号Mhh-decom-FFT、Mhv-decom-FFT、Mvh-decom-FFT、Mvv—de_-FFT沿每个距离门（按列）进行快速逆傅里叶变换（IFFT，得5IJ方位向回波信号具体表示为：[0082]式中，，表示HH通道方位向回波信号中第（1，1个像素的数值：g示HH通道方位向回波信号中第（I，Nr-ins个像素的数值表示HH通道方位向回波信号中第Na-ins，l个像素的数值，_____u表示HH通道方位向回波信号中第Na-ins，U个像素的数值：表示HV通道方位向回波信号中第（1，1个像素的数值，表示HV通道方位向回波信号中第（I，Nr-ins个像素的数值，表示HV通道方位向回波信号中第Na-ins，1个像素的数值表示HV通道方位向回波信号中第Na-ins，Nr-ins个像素的数值；，表示VH通道方位向回波信号中第（1，1个像素的数值：表示VH通道方位向回波信号中第（l，Nr-ins个像素的数值，表示VH通道方位向回波信号中第Na-ins，1个像素的数值，'表示VH通道方位向回波信号中第Na-ins，Nr-ins个像素的数值；表示VV通道方位向回波信号中第（1，1个像素的数值，t表示VV通道方位向回波信号中第（l，Nr-ins个像素的数值，表示VV通道方位向回波信号中第Na-ins，1个像素的数值，^表示VV通道方位向回波信号中第Na-ins，Nr-ins个像素的数值。IFFT·表示对二维矩阵的每列进行快速逆傅里叶变换。[0083]步骤六:估计法拉第旋转角；[0084]根据步骤五得到的四个极化通道方位向回波信号Mffl-ech。、MHV-ech。、MvH-ech。、Mw-ech。，通过线性组合得到一组正交圆极化波信号Z12和Z21，通过该组正交圆极化波信号共辄相乘取幅角FRtemp，对FRtemp进行平滑处理后得到法拉第旋转角估计值FR。具体为：[0085]1对四个极化通道方位向回波信号进行线性组合处理，计算一组正交圆极化波信号Zl2和Z21:[0088]式中，Z12和Z21分别代表计算得到的一组正交圆极化波信号，其均为大小SNa-insXNr-ins的二维复数矩阵。[0089]2对该组正交圆极化波信号共辄相乘取幅角FRtemp:[0091]式中，FRtemp表示利用正交圆极化波信号共辄相乘计算得到的幅角值，其为大小为Na-insXNr-ins的二维实数矩阵，*表示共辄，angleX为取复数X的幅角运算。[0092]⑶对上述计算得到的幅角值FRtemp二维实数矩阵），在仅考虑时变TEC的条件下首先将FRtemp按列叠加求平均值，再对得到的结果一维实数数组进行平滑处理得到法拉第旋转角估计值FR，其为长度为1-的一维实数数组：[0094]式中，FRtemp:，i表示幅角值矩阵FRtemt^i列所有元素，表示对X采用窗口大小为2M+1的窗函数进行平滑处理，即以X所在坐标p为中心，求前后第p-M到p+M个内的所有元素均值作为X值，需要注意的是，对于靠近数组边缘的元素，会存在平滑处理窗口超出数组边界的情况，即P-M彡0，p+MNa-ins，此时，该元素不做平滑处理直接以计算出X值作为该处结果。[0095]步骤七:获取合成孔径时间内的时变TEC;[0096]根据步骤六得到的估计法拉第旋转角值FR，利用时变TEC与法拉第旋转角间的线性转换关系，获得合成孔径时间内每个方位时刻的TEC，即合成孔径时间内的时变TEC:[0098]式中，FRη表示估计法拉第旋转角值FR中第η个元素值，Κω为常数且Κω=2.365XIO4AXm2kg，B表示平行于雷达波束方向的地球磁场强度，单位为Wbm2。[0099]经过以上七个步骤，利用全极化高轨SAR信号完成了对合成孔径时间内的时变TEC精确测量。[0100]实施例[0101]下面将结合附图和实施例对本发明作进一步的详细说明。[0102]本实施例提供一种基于全极化高轨SAR的电离层时变TEC测量方法。由于现有星载SAR系统没有直接获取的全极化高轨SAR数据，因此实施例中采用仿真得到全极化高轨SAR数据，并在信号方位向合成孔径时间内加入电离层时变TEC如图2所示所引起的法拉第旋转角，处理过程中所涉及的参数如表1所示：[0103]表1实施例参数[0104][0106]本实施例具体包含以下步骤：[0107]步骤一:计算方位向合成孔径点数；[0108]根据雷达系统参数，按式⑴计算方位向合成孔径点数为7768。[0109]步骤二:方位向数据补零；[0110]根据步骤一所计算得到的方位向合成孔径点数7768,按式⑹⑺⑻⑼分别在每个极化通道图像数据Mhh、Mhv、Mvh、Mvv每一列的头部和尾部补零，补零数目均为3884,得到新的四个极化通道图像数据分别为]V[HH_ins、MHV-ins、MvH-ins、MvV-ins。原始每个极化通道图像数据1^、1^、1^、1^均为大小为10241024的二维复数矩阵，经过方位向数据补零处理后，得到新的四个极化通道图像数据购11-、1取-、1^-、1^-^均为大小为8792\1024的二维复数矩阵。[0111]步骤三:方位向傅里叶变换；[0112]将步骤二得到的补零后的四个极化通道图像数据MHH-ins、Mm-ins、MvH-ins、Mw-ins沿每个距尚门（按列）进行快速傅里叶变换FFT，得到方位向频谱数据Mhh-fft、Mhv-fft、Mvh-fft、Mvv-FFTO[0113]步骤四:方位向信号解压缩，获取等效方位向回波频谱信号；[0114]根据表1给出的雷达系统参数，计算每个距离门（每列）对应的参考斜距Rref，以此为参考计算出方位向信号解压缩因子Φ。™，利用该解压缩因子乘以步骤三所得的方位向频谱数据Mhh-FFT、Mhv-FFT、Mvh-FFT、Mw-FFT得到等效方位向回波频谱信号MHH-decom-FFT、MHV-de_-FFT、MvH-decom-FFT、Mw-decom-FFT，本:[0115]I根据雷达系统参数，计算每个距离门海列对应的参考斜距Rref:[0116]根据式（14计算得到的每个距离门列对应的参考斜距Rref，式（14中处于同一距离门（每列）不同方位向（每行）的参考斜距数值是相同的，即每个距离门（每列采用相同的参考斜距。[0117]2雷达的脉冲重复频率PRF为120Hz，按式（15计算方位向频谱数据每个方位向每行的方位频率fa。[0118]⑶根据上述计算所得每个距离门海列对应的参考斜距Rref以及每个方位向海行的方位频率fa，按式16计算得到方位向信号解压缩因子Φ。™。[0119]⑷按式（17181920将计算得到的方位向信号解压缩因子乘以方位向步页谱数据MhH-FFT、MhV-FFT、MvH-FFT、Mvv-FFT得至丨J等效方位向回波步页谱信号Μ·-decom-FFT、MhV-decom-FFT、MvH-decom-FFT、MvV-decom-FFTo[0120]步骤五:方位向逆傅里叶变换；[0121]将步骤四得到的等效方位向回波频谱信号MHH—decom—FFT、MhV—decom—FFT、MvH—decom—FFT、Mvv—de_-FFT沿每个距离门（按列）进行快速逆傅里叶变换（IFFT，得5IJ方位向回波信号MhH-echo、MhV-echo、MvH-echo、Mw-echo〇[0122]步骤六:估计法拉第旋转角；[0123]根据步骤五得到的四个极化通道方位向回波信号Mffl-ech。、MHV-ech。、MvH-ech。、Mw-ech。，通过线性组合得到一组正交圆极化波信号Z12和Z21，对该组正交圆极化波信号共辄相乘取幅角FRtemp，进行平滑处理后得到所估计法拉第旋转角FR。具体为：[0124]1对四个极化通道方位向回波信号进行线性组合处理，按式25和26计算得到一组正交圆极化波信号Z12和Z21，其均为大小为8792X1024的二维复数矩阵。[0125]2按式27对该组正交圆极化波信号Z12和Z21进行共辄相乘并取幅角FRtemp。[0126]3对上述计算得到的幅角值FRtemp大小为8792X1024的二维实数矩阵），在仅考虑时变TEC的条件下，按式（28首先将FU安列叠加求平均值，再对得到的结果（长度为8792的一维实数数组采用长度为2M+1=9的平滑窗进行平滑处理，即以FRtemp按列叠加求平均值中以第P项为中心，求前后第P-4到p+4个内的所有元素均值作为第p项值，需要注意的是，对于该数组第1项到第4项，第8789项到8792项，由于平滑处理窗口超出数组边界，因此这些项不做平滑处理，直接以平滑前数值作为该项结果。由此，计算得到法拉第旋转角估计值FR，其为长度为8792的一维实数数组。[0127]步骤七:获取合成孔径时间内的时变TEC;[0128]根据步骤六得到的法拉第旋转角估计值FR，利用时变TEC与法拉第旋转角间的线性转换关系，按式29计算得到合成孔径时间内每个方位时刻的TEC，即合成孔径时间内的时变TEC，如图3所示。[0129]对比图1所示仿真中原始加入的合成孔径时间内时变TEC与根据本发明计算得到的图2所示估计得到的合成孔径时间内时变TEC，可以看出两者极为接近。统计采用本发明方法计算得到的时变TEC估计误差标准差为0.0287TECU，即可实现高精度时变测量。同时，由于所实施例中时变TEC在合成孔径时间内以雷达的脉冲重复周期0.0083s为间隔变化，因此可见本发明提出的测量方法可实现毫秒级的时变TEC测量。[0130]综上，上述实施例说明本发明提出的一种基于全极化高轨SAR的电离层时变TEC测量方法，可以实现时变TEC高精度测量。

权利要求：1.一种基于全极化高轨SAR的电离层时变TEC测量方法，其特征在于，包括以下几个步骤：步骤一:计算方位向合成孔径点数；根据雷达系统参数，计算方位向合成孔径点数Num_a;步骤二:方位向数据补零；根据步骤一计算得到的方位向合成孔径点数Num_a，分别在每个极化通道图像数据每一列的头部和尾部补零，补零数目均为Num_a2,得到新的四个极化通道图像数据分别为MHH-insMhV-ins'«MvH-ins'«MvV-ins；步骤三:方位向傅里叶变换；将步骤二得到的补零后的四个极化通道图像数据MHH-ins、MHV-ins、MVH-ins、MVV-in|每个距离门（按列进行快速傅里叶变换FFT，得到方位向频谱数据Mhh-fft、Mhv-fft、Mvh-fft、Mvv-fft;步骤四:方位向信号解压缩，获取等效方位向回波频谱信号；根据雷达系统参数，计算每个距离门对应的参考斜距Rref，进而计算出方位向信号解压缩因子Φdec™，利用该解压缩因子乘以步骤三所得的方位向频谱数据Mhh-fft、Mhv-fft、Mvh-fft、Mvv-FFT得到等效方位向回波频谱信号MHH-de_-FFT、MHV-de_-FFT、MvH-de_-FFT、Mw-de_-FFT;步骤五:方位向逆傅里叶变换；将步骤四得到的等效方位向回波频谱信号MHH-decom-FFT、MfiV-decom-FFT、MvH-decom-FFT、Mw-decom-FFT沿每个距离门进行快速逆傅里叶变换，得到方位向回波信号MHH-eCho、MHV-ech。、MvH-echo、Mw-echo;步骤六:估计法拉第旋转角；根据步骤五得到的四个极化通道方位向回波信号MhH-ech。、MhV-ech。、MVH-ech。、Mw-ech。，通过线性组合得到一组正交圆极化波信号Z12和Z21，通过该组正交圆极化波信号共辄相乘取幅角FRtemp，对FRtemp进行平滑处理后得到法拉第旋转角估计值FR;步骤七:获取合成孔径时间内的时变TEC;根据步骤六得到的估计法拉第旋转角值FR，利用时变TEC与法拉第旋转角间的线性转换关系，获得合成孔径时间内每个方位时刻的TEC，即合成孔径时间内的时变TEC。2.根据权利要求1所述的一种基于全极化高轨SAR的电离层时变TEC测量方法，其特征在于，所述的步骤一具体包括：根据雷达系统参数，计算方位向合成孔径点数Num_a:式中，λ表示雷达系统工作波长，R。表示场景中心点的参考斜距，La表示方位向天线长度，Vrrf表示以场景中心点为参考时的等效速度，PRF表示雷达脉冲重复频率;fl〇〇rx表示取不大于X的最大偶数运算。3.根据权利要求1所述的一种基于全极化高轨SAR的电离层时变TEC测量方法，其特征在于，所述的步骤二具体为：设原始四个极化通道数据Mhh、Mhv、Mvh、Mvv分别为：式中，表示HH通道图像中第（I，1个像素的数值，表示HH通道图像中第（I，Nr个像素的数值^表示HH通道图像中第Na，l个像素的数值，表示HH通道图像中第Na，Nr个像素的数值；表示HV通道图像中第（1，1个像素的数值，表示HV通道图像中第（I，Nr个像素的数值，表示HV通道图像中第Na，1个像素的数值，表示HV通道图像中第Na，Nr个像素的数值；表示VH通道图像中第1，1个像素的数值，表示VH通道图像中第（I，Nr个像素的数值，表示VH通道图像中第Na，1个像素的数值，表示VH通道图像中第Na，Nr个像素的数值：表示VV通道图像中第（1，1个像素的数值，表示VV通道图像中第（I，Nr个像素的数值，表示VV通道图像中第Na，1个像素的数值，表示VV通道图像中第Na，Nr个像素的数值；结合步骤一计算得到的方位向合成孔径点数Num_a，分别在四个极化通道数据每一列的头部和尾部分别补零Num_a2,得到补零后新的四个极化通道图像数据MHH-ins、MHV-ins、MvH-ins、Mw-ins，其均扩大为大*Na-insXNr-ins的二维复数矩阵，且Na-ins=Na+Num_a,?Nr-ins一Nr，四个通道的信号分别为：式中，表示补零后HH通道图像中第（Num_a2+l，1个像素的数值，表示补零后HH通道图像中第Num_a2+l，Nr个像素的数值，表示补零后HH通道图像中第Num_a2+Na，1个像素的数值，表示补零后HH通道图像中第Num_a2+Na，Nr个像素的数值；表示补零后HV通道图像中第Num_a2+1，1个像素的数值，表示补零后HV通道图像中第Num_a2+l，Nr个像素的数值，表示补零后HV通道图像中第Num_a2+Na，1个像素的数值，表示补零后HV通道图像中第Num_a2+Na，Nr个像素的数值；表示补零后VH通道图像中第Num_a2+l，1个像素的数值，表示补零后VH通道图像中第Num_a2+1，Nr个像素的数值，表示补零后VH通道图像中第Num_a2+Na，1个像素的数值，表示补零后VH通道图像中第Num_a2+Na，Nr个像素的数值；表示补零后VV通道图像中第Num_a2+l，1个像素的数值，表示补零后VV通道图像中第Num_a2+l，Nr个像素的数值，表示补零后VV通道图像中第Num_a2+Na，1个像素的数值，表示补零后VV通道图像中第Num_a2+Na，Nr个像素的数值。4.根据权利要求1所述的一种基于全极化高轨SAR的电离层时变TEC测量方法，其特征在于，所述的步骤三具体为：式中，表示HH通道方位向频谱数据中第（1，1个像素的数值，表示HH通道方位向频谱数据中第（I，Nr-ins个像素的数值，表示HH通道方位向频谱数据中第Na-ins，1个像素的数值，表示HH通道方位向频谱数据中第Na-ins，Nr-ins个像素的数值；表示HV通道方位向频谱数据中第1，1个像素的数值，表示HV通道方位向频谱数据中第（l，Nr-ins个像素的数值，表示HV通道方位向频谱数据中第Na-ins，1个像素的数值，表示HV通道方位向频谱数据中第Na-ins，Nr-ins个像素的数值I表示VH通道方位向频谱数据中第（1，1个像素的数值，.1表示VH通道方位向频谱数据中第（l，Nr-ins个像素的数值，表示VH通道方位向频谱数据中第Na-ins，1个像素的数值，表示VH通道方位向频谱数据中第Na-ins，Nr-ins个像素的数值；表示VV通道方位向频谱数据中第（1，1个像素的数值，.表示VV通道方位向频谱数据中第（l，Nr-ins个像素的数值，表示VV通道方位向频谱数据中第Na-ins，1个像素的数值：表示VV通道方位向频谱数据中第Na-ins，Nr-ins个像素的数值;FFT·表示对二维矩阵的每列进行快速傅里叶变换。5.根据权利要求1所述的一种基于全极化高轨SAR的电离层时变TEC测量方法，其特征在于，所述的步骤四具体包括：1根据雷达系统参数，计算每个距离门对应的参考斜距Rref:式中，c表不光速；⑵雷达的脉冲重复频率为PRF，计算方位向频谱数据每个方位向的方位频率匕为：3根据上述计算所得每个距离门对应的参考斜距Rref以及每个方位向的方位频率fa，得到方位向信号解压缩因子Φd_m:式中，表示方位频率匕一维数组中的第1个频率值，表示方位频率匕一维数组中的第Na-ins个频率值，表示参考斜距Rref二维矩阵中第（1，1个斜距值，表示参考斜距Rref二维矩阵中第（I，Nr-ins个斜距值，表示参考斜距Rref二维矩阵中第Na-ins，1个斜距值，表示参考斜距Rref二维矩阵中第Na-ins，Nr-ins个斜距值；⑷将计算得到的方位向信号解压缩因子Φde_乘以方位向频谱数据Mhh-fft、Mhv-fft、Mvh-FFT、Mw-FFT得到等效方位向回波频谱信号MHH-de_-FFT、MHV-de_-FFT、MvH-de_-FFT、Mw-de_-FFT:式中，j表示虚数单位，即表示方位向信号解压缩因子Φd_m二维矩阵中第（1，1个因子数值，表示方位向信号解压缩因子二维矩阵中第（I，Nr-ins个因子数值，表示方位向信号解压缩因子二维矩阵中第Na-ins，l个因子数值，表示方位向信号解压缩因子Φι™二维矩阵中第Na-ins，Nr-ins个因子数值。6.根据权利要求1所述的一种基于全极化高轨SAR的电离层时变TEC测量方法，其特征在于，所述的步骤五具体为：式中，，表示HH通道方位向回波信号中第（1，1个像素的数值，表示HH通道方位向回波信号中第（l，Nr-ins个像素的数值，表示HH通道方位向回波信号中第Na-ins，1个像素的数值，表示HH通道方位向回波信号中第Na-ins，Nr-ins个像素的数值；表示HV通道方位向回波信号中第（1，1个像素的数值，，表示HV通道方位向回波信号中第（l，Nr-ins个像素的数值，，表示HV通道方位向回波信号中第Na-ins，1个像素的数值，表示HV通道方位向回波信号中第Na-ins，Nr-ins个像素的数值；表示VH通道方位向回波信号中第（1，1个像素的数值，表示VH通道方位向回波信号中第（l，Nr-ins个像素的数值，表示VH通道方位向回波信号中第Na-ins，1个像素的数值，表示VH通道方位向回波信号中第Na-ins，Nr-ins个像素的数值；，表示VV通道方位向回波信号中第（1，1个像素的数值，表示VV通道方位向回波信号中第（l，Nr-ins个像素的数值，表示VV通道方位向回波信号中第Na-ins，1个像素的数值，表示VV通道方位向回波信号中第Na-ins，Nr-ins个像素的数值;IFFT·表示对二维矩阵的每列进行快速逆傅里叶变换。7.根据权利要求1所述的一种基于全极化高轨SAR的电离层时变TEC测量方法，其特征在于，所述的步骤六具体为：1对四个极化通道方位向回波信号进行线性组合处理，计算一组正交圆极化波信号Zl2和Z21:Zl2—MvH-echo_MHV-ech。+jXMhH-echo+Mw-echo25Z21—MhV-echo_MvH-ech。+jXMhH-echo+Mw-echo26式中，Z12和Z21分别代表计算得到的一组正交圆极化波信号，其均为大小SNa-insXNr-ins的二维复数矩阵；⑵对该组正交圆极化波信号共辄相乘取幅角FRtemp:式中，FRtemp表示利用正交圆极化波信号共辄相乘计算得到的幅角值，其为大小为Na-insXNr-ins的二维实数矩阵，*表示共辄，angleX为取复数X的幅角运算；3首先将FRte3mp按列叠加求平均值，再对得到的结果进行平滑处理得到法拉第旋转角估计值FR，其为长度为Na-ins的一维实数数组：式中，FRte3mp:，i表示幅角值矩阵?心_第1列所有元素，表示对X采用窗口大小为2M+1的窗函数进行平滑处理，即以X所在坐标p为中心，求前后第p-M到p+M个内的所有元素均值作为X值。8.根据权利要求1所述的一种基于全极化高轨SAR的电离层时变TEC测量方法，其特征在于，所述的步骤七中，获得合成孔径时间内每个方位时刻的TEC为：式中，FRη表示估计法拉第旋转角值FR中第η个元素值，Κω为常数且Κω=2.365XIO4AXm2kg，B表示平行于雷达波束方向的地球磁场强度，单位为Wbm2。

百度查询： 北京航空航天大学 一种基于全极化高轨SAR的电离层时变TEC测量方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD