CN107102303A - 机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法 - Google Patents
机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法 Download PDFInfo
- Publication number
- CN107102303A CN107102303A CN201610097343.3A CN201610097343A CN107102303A CN 107102303 A CN107102303 A CN 107102303A CN 201610097343 A CN201610097343 A CN 201610097343A CN 107102303 A CN107102303 A CN 107102303A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- baseline
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000035945 sensitivity Effects 0.000 claims description 8
- 230000003068 static effect Effects 0.000 claims description 4
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 239000011159 matrix material Substances 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 101100001673 Emericella variicolor andH gene Proteins 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供了一种机载混合顺轨‑交轨干涉合成孔径雷达定标方法,该方法针对机载混合顺轨‑交轨干涉合成孔径雷达,建立了新的基线矢量模型,根据此模型和载机姿态角信息能够描述混合干涉SAR系统任意时刻的基线矢量;在新的基线矢量模型的基础上,对混合干涉SAR系统的待定标参数进行了修改,并提出了适用于机载混合顺轨‑交轨干涉合成孔径雷达的定标方法,该方法通过对某一时刻基线分量的标定来实现对基线长度和两个基线角的定标,避免了复杂的公式计算,更加容易实现。
Description
技术领域
本发明属于合成孔径雷达(Synthetic Aperture Radar,SAR)技术领域,具体涉及一种机载混合顺轨-交轨干涉合成孔径雷达定标方法。
背景技术
合成孔径雷达是一种先进的主动式微波对地观测手段,因具有全天时、全天候、高分辨率以及对地面植被具有一定的穿透性等工作能力而得到广泛应用。干涉合成孔径雷达(Interferometric Synthetic Aperture Radar,干涉SAR或InSAR)是无线电干涉测量技术与合成孔径雷达技术的结合,其利用两个天线获取的雷达数据进行干涉处理获得地形高程、目标径向速度、地形形变等信息。干涉SAR系统中两个天线相位中心之间的相对位置称为干涉基线,按照基线类型干涉SAR可以分为两个天线相位中心沿航迹方向放置、对运动目标径向速度敏感的顺轨干涉SAR(Along-Track Interferometric SAR,ATI-SAR)以及两个天线相位中心沿垂直航迹方向放置、对地形高程敏感的交轨干涉SAR(Cross-TrackInterferometric SAR,XTI-SAR)。混合顺轨-交轨干涉SAR(Hybrid Along-track andCross-track Interferometric Synthetic Aperture Radar,简称AT-XI InSAR)系统为两个天线相位中心之间的基线矢量既包含了沿航迹方向的顺轨基线分量、又包含了垂直航迹方向的交轨基线分量的干涉SAR系统,能够同时用于地面动目标检测及地形高程反演。
干涉定标是干涉SAR获得高精度测量结果的关键技术之一。目前关于干涉SAR定标的研究主要是针对机载双天线XTI-SAR系统,如2009年电子所张薇博士的学位论文《机载双天线干涉SAR定标方法研究》,XTI-SAR系统的基线矢量模型为[B,α],其中B为基线长度,α为基线角,待定标参数为[Δφ,B,α],Δφ为干涉相位偏差,定标方法为经典的基于敏感度方程的定标方法,利用定标场控制点的高程信息直接对待定标参数进行定标。
上述干涉SAR定标技术不适用于混合干涉SAR系统的定标,主要体现在:
第一、基线矢量模型对AT-XT InSAR系统而言不充分,具体来说,混合AT-XT InSAR系统由于同时具有顺轨基线分量和交轨基线分量,相比于单纯的ATI-SAR或者XTI-SAR系统,其基线矢量模型更加复杂,上述基线矢量模型不能全面、正确地描述AT-XT InSAR系统的基线特点,因此需要针对混合AT-XT InSAR系统对基线进行重新建模。
第二,同上,由于上述基线矢量模型不适用于机载AT-XT InSAR系统,因此上述待定标的参数也不能涵盖全部的AT-XT InSAR系统待定标参数,待定标参数需要进行相应的修改。
第三、定标方法只是针对XTI-SAR系统进行定标,只属于AT-XT InSAR系统顺轨基线为零时的一种特殊情况,需要研究适用于混合AT-XT InSAR系统的定标方法。
发明内容
(一)要解决的技术问题
本发明的目的在于,提供一种机载混合AT-XT InSAR系统的定标方法,能够正确、全面地对机载混合AT-XT InSAR系统的干涉参数进行定标。
(二)技术方案
本发明提供一种机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法,用于对机载混合顺轨-交轨干涉合成孔径雷达系统中基线矢量模型的干涉参数进行定标,其中,机载混合顺轨-交轨干涉合成孔径雷达系统包括前置天线和后置天线,将前置天线作为参考天线,后置天线与前置天线之间的相对位置关系用基线矢量模型表示,干涉参数包括前置天线和后置天线之间的基线长度、顺轨基线角和交轨基线角,其中,基线长度是指前置天线相位中心与后置天线相位中心之间的长度,顺轨基线角是指基线矢量与yoz平面的夹角,交轨基线角是指基线矢量在yoz平面的投影与y轴之间的夹角,方法包括:
S1,标定顺轨基线分量;
S2,标定干涉相位偏置、交轨基线长度和交轨基线角;
S3,根据交轨基线长度和交轨基线角,计算交轨基线分量;
S4,根据顺轨基线分量和交轨基线分量,计算载机坐标系下的基线矢量;
S5,根据载机坐标系下基线矢量,计算基线长度和顺轨基线角的定标结果。
(三)有益效果
1、本发明针对机载AT-XT InSAR系统建立了包含了基线长度、顺轨基线角和交轨基线角的新基线矢量模型,该模型能够全面地描述机载AT-XT InSAR系统的基线特点,且能够根据载机姿态信息获得任意时刻的基线矢量表达式。
2、本发明的基线矢量模型中的参数能够正确、全面地描述机载AT-XT InSAR系统的定标需求,从而为定标方法研究提供基础。
3、本发明通过标定某一时刻的三个基线分量间接实现了基线长度和两个基线角的标定,避免了复杂的公式计算,实现了机载AT-XT InSAR系统干涉参数的定标。
附图说明
图1是本发明载机坐标系中混合基线InSAR基线矢量模型。
图2A是本发明载机混合AT-XT InSAR系统的工作几何图。
图2B是本发明方位校准后载机混合AT-XT InSAR系统的工作几何图。
图3是本发明提供的机载混合AT-XT InSAR系统的定标方法的流程图。
图4是本发明中基于敏感度方程的交轨基线分量定标算法流程图。
具体实施方式
本发明提供一种机载混合顺轨-交轨干涉合成孔径雷达定标方法,该方法针对机载混合顺轨-交轨干涉合成孔径雷达,建立了新的基线矢量模型,根据此模型和载机姿态角信息能够描述混合干涉SAR系统任意时刻的基线矢量;在新的基线矢量模型的基础上,对混合干涉SAR系统的待定标参数进行了修改,并提出了适用于机载混合顺轨-交轨干涉合成孔径雷达的定标方法,该方法通过对某一时刻基线分量的标定来实现对基线长度和两个基线角的定标,避免了复杂的公式计算,更加容易实现。
图1是本发明载机坐标系中混合基线InSAR基线矢量模型,如图1所示,x方向为机头方向,z为垂直机身向上的方向,y与二者成右手准则。图中A1、A2分别为两个天线相位中心的位置,B代表两个天线相位中心之间的物理基线长度,基线角α1为基线矢量与yoz平面的夹角,基线角α2为基线矢量在yoz平面的投影与y轴之间的夹角。如图1所示,在载机坐标系下,基线矢量可以表示为:
其中,Bx1、By1、Bz1分别为载机坐标系下基线矢量在三个坐标轴方向的分量大小。上式写成矩阵形式为:
其中,[Γ1]、[Γ2]为与基线角相关的旋转矩阵:
由公式(1)可以看出,由于基线角α2只影响交轨基线分量的大小,因此称之为交轨基线角,其与XTI-SAR中的基线角α相对应。与XTI-SAR系统相比,此模型中多了一个基线角α1,此基线角非零时会导致顺轨干涉基线分量Bx的出现,实际上α1也影响交轨基线分量的大小,这里为了与交轨基线角对应(区分),称α1为顺轨基线角。
实际飞行中由于气流等影响,载机会发生偏航、俯仰、滚转的姿态变化,此时瞬时基线分量在大地坐标系中也会发生相应变化。受姿态影响后大地坐标系中瞬态基线矢量变为:
式中为大地坐标系中的瞬态基线矢量,矩阵[Y]、[P]、[R]分别表示与偏航角、俯仰角和滚转角相关的旋转矩阵:
其中θy为偏航角,θp为俯仰角,θr为滚转角。在ATI-SAR系统中,主要考虑偏航角和俯仰角影响,而在XTI-SAR系统中,则主要考虑滚转角的影响。在本发明的混合AT-XT InSAR系统中,偏航角、俯仰角和滚转角对基线分量的影响均要加以考虑。
由上可知,在本发明所建立的包含了基线长度B、顺轨基线角α1和交轨基线角α2的基线矢量模型下,结合载机运动过程中的姿态角信启、,可以根据公式(2)和公式(5)得到任意时刻的顺轨基线分量、跨航向基线分量和竖直方向基线分量的大小。
图2A是本发明载机混合AT-XT InSAR系统的工作原理图,如图2A所示,两个天线先后对同一块区域进行成像,与传统的顺轨干涉SAR系统不同,此处后置天线2并不按照前置天线1的轨迹运动,而是存在跨航向和竖直方向的轨迹差。假设两个天线在飞行方向的基线分量为Bx,载机速度vp设为常数,则经历时间间隔Δt=Bx/vp后,两个天线之间形成传统的对地形高程敏感的交轨干涉系统,其中交轨基线长度记为B⊥,交轨基线角仍为α2。图2B是本发明方位校准后载机混合AT-XT InSAR系统的工作原理图,如图2B所示,A1、A2分别代表两个天线的相位中心,P代表地面目标,H为参考天线(天线1)相对参考地形的高度,h为目标所在高程,R为参考天线与目标之间的最短斜距,θ为下视角,ΔR为两个天线到目标点的斜距之差。
根据AT-XT InSAR系统工作几何,可以得到在时间间隔Δt内,两个天线相位中心到同一点目标的距离差为:
ΔR=vrΔt-Bysinθ+Bzcosθ (6)
由此引入的相位差为:
其中λ为波长,Δt=Bx/vp为两个天线相位中心运动到同一方位位置所需的时间间隔,vr为运动目标的径向速度,假设其在Δt时间间隔内为常数。θ为参考天线(天线1)对应的入射角,且H为平台相对参考高程的高度,h为目标所处的高度,R为参考天线与目标之间的最短斜距。φATI为顺轨干涉相位,φXTI为交轨干涉相位,且:
其中B⊥为交轨基线长度,B⊥=Bcosα1,相当于XTI-SAR系统中的基线长度。
公式(7)写成向量的形式为:
对于机载AT-XT InSAR系统,需要定标的参数主要包括[Δφ,B,α1,α2,θy,θp,θr],即干涉相位偏置Δφ、物理基线长度B、顺轨基线角α1、交轨基线角α2以及姿态角[θy,θp,θr]。其中载机的姿态角信息可以根据载机上的姿态测量系统直接获得,如果测量系统的测量精度满足要求可以不对其进行定标。为便于说明,在本发明中假设姿态测量系统的测量误差满足系统精度要求,此时只需要对参数[Δφ,B,α1,α2]进行定标处理。
由于运动的地面控制点速度不易控制,为减少目标运动速度控制引入的误差,我们采用静止的地面控制点对上述系统参数进行定标,如静止的角反射器。此时干涉相位的表达式为:
即干涉相位对目标高程敏感,目标高程为h=H-Rcosθ,其中θ可根据(11)式通过干涉相位得出。将(11)式展开可以发现,包含了姿态角信息的目标高程的解析表达式很难得出,因此直接根据(11)式标定基线长度B和两个基线角α1、α2不易实现。因此,本发明提出如图3所示的定标方法对[Δφ,B,α1,α2]进行间接定标。
如图3所示,本发明所提定标方法的实施步骤具体如下:
第一步,选择定标场,并沿跨航向放置N(N≥3)个静止的角反射器;
第二步,对定标场获取的数据进行方位向时间校准,根据方位向时间校准量和载机平台速度标定顺轨基线分量Bx,Bx=vp·Δt,其中vp为平台的方位向速度,Δt为后置天线运动到前置天线同样的方位向位置所需的时间间隔;
第三步,利用敏感度方程标定干涉相位偏置、交轨基线长度和交轨基线角[Δφ,B⊥,α2]。将(11)式用[Bx,By,Bz]表示,则公式(11)变为公式(9),进而可推得以及目标高程的表达式:h=H-Rcosθ。
敏感度方程为:
其中,分别为控制点高程关于相位偏置Δφ、交轨基线长度B⊥和交轨基线角α2的偏导数,其构成的矩阵即式(12)中等式左边的第一个矩阵称为敏感度矩阵,用F表示。Δhi为根据干涉相位得到的目标高度与目标实际高度之间的误差。这里i表示定标点目标的编号。
通过迭代地求解敏感度方程(12)的解,即可得到定标后的[Δφ,B⊥,α2]。定标过程如图4所示,具体步骤如下:
步骤①:根据得到的干涉相位和系统参数估计目标高度hest,同时,根据定标点高度信息得到定标点的真实高度hreal,二者相减得到高程差Δh:
步骤②:利用干涉参数值得到敏感度矩阵F,利用方程组(12)计算参数偏差[Δφ,ΔB⊥,Δα2]的最小二乘解;
步骤③:利用得到的[Δφ,ΔB⊥,Δα2]更新系统参数[Δφ,B⊥,α2],重新得到定标点高度的估计值hest;
步骤④:重新计算估计的定标点高度与定标点真实高度之间的误差Δh,如果Δh满足精度要求则退出定标过程,否则进行步骤②继续进行迭代定标。
第四步,根据敏感度方程的算法流程得到的[B⊥,α2]的定标结果,计算交轨基线分量Bv,Bz的值:
第五步,根据第二步和第四步的定标结果可以得到瞬态基线矢量的表达式;根据载机上安装的测量系统记录的姿态角信息(θy,θp,θr)可以计算得到欧拉矩阵E。利用瞬态基线与载机坐标系下的基线矢量之间的关系(如公式(4)所示),计算得到载机坐标系下的基线矢量:
第六步,利用载机坐标系下基线矢量的表达式与基线长度和两个基线角的关系(如公式(1)所示),得到基线长度B和基线角α1的定标结果:
通过上述定标步骤,即可得到定标后的干涉参数[Δφ,B,α1,α2],完成机载混合AT-XT InSAR系统的定标。
与目前常用的机载XTI-SAR系统定标方法相比较,本发明通过对机载混合AT-XTInSAR系统的基线进行重新建模,并修改相应的待定标参数,在此基础上提出了一种适用于机载AT-XT InSAR系统的定标方法。所取得的有益效果为:
1、本发明针对机载AT-XT InSAR系统建立了包含了基线长度和顺轨基线角、交轨基线角的新基线矢量模型,该模型能够全面地描述机载AT-XT InSAR系统的基线特点,且能够根据载机姿态信息获得任意时刻的基线矢量表达式。
2、按照本发明中的基线矢量模型,将待定标参数进行了相应的修改,修改后的待定标参数能够正确、全面地描述机载AT-XT InSAR系统的定标需求,从而为定标方法研究提供基础。
3、本发明提出了一种针对机载AT-XT InSAR系统的定标方法,该方法通过标定某一时刻的三个基线分量间接实现了基线长度和两个基线角的标定,避免了复杂的公式计算,实现了机载AT-XT InSAR系统干涉参数的定标。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法,用于对机载混合顺轨-交轨干涉合成孔径雷达系统中基线矢量模型的干涉参数进行定标,其中,所述机载混合顺轨-交轨干涉合成孔径雷达系统包括前置天线和后置天线,其特征在于,将前置天线作为参考天线,后置天线与前置天线之间的相对位置关系用基线矢量模型表示,所述干涉参数包括前置天线和后置天线之间的基线长度、顺轨基线角和交轨基线角,其中,基线长度是指前置天线相位中心与后置天线相位中心之间的长度,顺轨基线角是指基线矢量与yoz平面的夹角,交轨基线角是指基线矢量在yoz平面的投影与y轴之间的夹角,方法包括:
S1,标定顺轨基线分量;
S2,标定干涉相位偏置、交轨基线长度和交轨基线角;
S3,根据交轨基线长度和交轨基线角,计算交轨基线分量;
S4,根据顺轨基线分量和交轨基线分量,计算载机坐标系下的基线矢量;
S5,根据载机坐标系下基线矢量,计算基线长度和顺轨基线角的定标结果。
2.根据权利要求1所述的机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法,其特征在于,所述步骤S1包括,从定标场获取定标数据,并对该定标数据进行方位向时间校准,得到方位向时间校准量Δt,根据方位向时间校准量Δt和载机平台速度vp标定顺轨基线分量Bx:
Bx=vp·Δt
其中,所述定标场中放置有至少3个静止的角反射器,用于获取所述定标数据。
3.根据权利要求2所述的机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法,其特征在于,所述步骤S2包括,利用敏感度方程标定干涉相位偏置Δφ、交轨基线长度B⊥和交轨基线角α2,其中,所述敏感度方程为:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>&phi;</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>B</mi>
<mo>&perp;</mo>
</msub>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>&phi;</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>B</mi>
<mo>&perp;</mo>
</msub>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>&phi;</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mi>N</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>B</mi>
<mo>&perp;</mo>
</msub>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mi>N</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mi>N</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>&phi;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;B</mi>
<mo>&perp;</mo>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>h</mi>
<msub>
<mo>|</mo>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>h</mi>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>h</mi>
<msub>
<mo>|</mo>
<mi>N</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中,h表示目标所处的高度,其表达式为:
h=H-Rcosθ
H为平台相对参考高程的高度,R为系统中参考天线与目标之间的最短斜距,θ为参考天线对应的入射角,其表达式为:
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mi>arcsin</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&lambda;</mi>
<mi>&phi;</mi>
</mrow>
<mrow>
<mn>4</mn>
<msub>
<mi>&pi;B</mi>
<mo>&perp;</mo>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
其中,λ为波长,φ为两个天线之间的相位差,其表达式为:
<mrow>
<mi>&phi;</mi>
<mo>=</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>A</mi>
<mi>T</mi>
<mi>I</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>X</mi>
<mi>T</mi>
<mi>I</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mi>r</mi>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>B</mi>
<mi>y</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>+</mo>
<msub>
<mi>B</mi>
<mi>z</mi>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
φATI为顺轨干涉相位,φXTI为交轨干涉相位;
分别为控制点高程关于相位偏置Δφ、交轨基线长度B⊥和交轨基线角α2的偏导数,Δhi为根据干涉相位得到的目标高度与目标实际高度之间的误差,i表示定标点目标的编号。
4.根据权利要求3所述的机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法,其特征在于,所述步骤S3中,交轨基线分量By,Bz的表达式为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<msub>
<mi>B</mi>
<mo>&perp;</mo>
</msub>
<msub>
<mi>cos&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mi>z</mi>
</msub>
<mo>=</mo>
<msub>
<mi>B</mi>
<mo>&perp;</mo>
</msub>
<msub>
<mi>sin&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
5.根据权利要求4所述的机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法,其特征在于,所述步骤S4中,载机坐标系下的基线矢量的表达式为:
<mrow>
<msub>
<mover>
<mi>B</mi>
<mo>&RightArrow;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>x</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>y</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>z</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>&rsqb;</mo>
<mo>=</mo>
<msup>
<mi>E</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mover>
<mi>B</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>,</mo>
</mrow>
其中,E为欧拉矩阵,其由姿态角信息计算得到。
6.根据权利要求5所述的机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法,其特征在于,所述步骤S5中,根据载机坐标系下基线矢量计算基线长度B和顺轨基线角α1的定标结果:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>arctan</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>B</mi>
<mrow>
<mi>x</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>B</mi>
<mrow>
<mi>y</mi>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>cos&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>B</mi>
<mo>=</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>x</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>sin&alpha;</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
2
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610097343.3A CN107102303B (zh) | 2016-02-23 | 2016-02-23 | 机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610097343.3A CN107102303B (zh) | 2016-02-23 | 2016-02-23 | 机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107102303A true CN107102303A (zh) | 2017-08-29 |
CN107102303B CN107102303B (zh) | 2024-03-19 |
Family
ID=59659036
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610097343.3A Active CN107102303B (zh) | 2016-02-23 | 2016-02-23 | 机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107102303B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108007476A (zh) * | 2017-11-20 | 2018-05-08 | 中国科学院空间应用工程与技术中心 | 一种天基干涉成像雷达高度计的干涉定标方法及系统 |
CN108061891A (zh) * | 2017-12-04 | 2018-05-22 | 上海无线电设备研究所 | 一种无控制点的干涉sar基线矢量估计方法 |
CN110007302A (zh) * | 2019-03-27 | 2019-07-12 | 北京空间飞行器总体设计部 | 一种星载双天线斜视顺轨干涉sar洋流速度测量方法 |
CN110823191A (zh) * | 2019-10-08 | 2020-02-21 | 北京空间飞行器总体设计部 | 混合基线双天线斜视干涉sar洋流测量性能确定方法及系统 |
CN110865346A (zh) * | 2018-08-28 | 2020-03-06 | 中国人民解放军61646部队 | 一种基于直接定位算法的星载sar时间参数标定方法 |
CN110907932A (zh) * | 2019-11-26 | 2020-03-24 | 上海卫星工程研究所 | 分布式InSAR卫星测高精度影响因素分析方法及系统 |
CN111289957A (zh) * | 2020-03-10 | 2020-06-16 | 上海高仙自动化科技发展有限公司 | 外参标定方法及装置、智能机器人及计算机可读存储介质 |
CN112327261A (zh) * | 2020-10-22 | 2021-02-05 | 上海卫星工程研究所 | 分布式InSAR卫星时间同步在轨测试方法及系统 |
CN115267780A (zh) * | 2022-09-26 | 2022-11-01 | 中国电子科技集团公司第十四研究所 | 一种同时实现顺轨和切轨干涉的多通道极化sar系统及其方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1301969A (zh) * | 1999-12-30 | 2001-07-04 | 中国科学院空间科学与应用研究中心 | 顺轨方向多基线干涉式综合孔径微波辐射计及其设计方法 |
CN103064083A (zh) * | 2011-10-24 | 2013-04-24 | 中国科学院电子学研究所 | 毫米波交轨三孔径稀疏阵sar系统的侧视三维成像法 |
CN103197316A (zh) * | 2013-03-26 | 2013-07-10 | 中国科学院电子学研究所 | 一种保相的机载干涉sar地形高度自适应运动补偿方法 |
CN103728621A (zh) * | 2014-01-24 | 2014-04-16 | 中国科学院电子学研究所 | 一种采用顺轨干涉处理抑制平台振动的机载sal成像方法 |
US20150094944A1 (en) * | 2013-10-01 | 2015-04-02 | Jaybridge Robotics, Inc. | Computer-implemented method and system for dynamically positioning a vehicle relative to another vehicle in motion for on-the-fly offloading operations |
CN105301570A (zh) * | 2015-10-20 | 2016-02-03 | 中国科学院电子学研究所 | 一种机载顺轨干涉sar系统的外场定标方法 |
-
2016
- 2016-02-23 CN CN201610097343.3A patent/CN107102303B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1301969A (zh) * | 1999-12-30 | 2001-07-04 | 中国科学院空间科学与应用研究中心 | 顺轨方向多基线干涉式综合孔径微波辐射计及其设计方法 |
CN103064083A (zh) * | 2011-10-24 | 2013-04-24 | 中国科学院电子学研究所 | 毫米波交轨三孔径稀疏阵sar系统的侧视三维成像法 |
CN103197316A (zh) * | 2013-03-26 | 2013-07-10 | 中国科学院电子学研究所 | 一种保相的机载干涉sar地形高度自适应运动补偿方法 |
US20150094944A1 (en) * | 2013-10-01 | 2015-04-02 | Jaybridge Robotics, Inc. | Computer-implemented method and system for dynamically positioning a vehicle relative to another vehicle in motion for on-the-fly offloading operations |
CN103728621A (zh) * | 2014-01-24 | 2014-04-16 | 中国科学院电子学研究所 | 一种采用顺轨干涉处理抑制平台振动的机载sal成像方法 |
CN105301570A (zh) * | 2015-10-20 | 2016-02-03 | 中国科学院电子学研究所 | 一种机载顺轨干涉sar系统的外场定标方法 |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108007476A (zh) * | 2017-11-20 | 2018-05-08 | 中国科学院空间应用工程与技术中心 | 一种天基干涉成像雷达高度计的干涉定标方法及系统 |
CN108007476B (zh) * | 2017-11-20 | 2020-05-01 | 中国科学院空间应用工程与技术中心 | 一种天基干涉成像雷达高度计的干涉定标方法及系统 |
CN108061891A (zh) * | 2017-12-04 | 2018-05-22 | 上海无线电设备研究所 | 一种无控制点的干涉sar基线矢量估计方法 |
CN108061891B (zh) * | 2017-12-04 | 2019-10-18 | 上海无线电设备研究所 | 一种无控制点的干涉sar基线矢量估计方法 |
CN110865346B (zh) * | 2018-08-28 | 2021-07-20 | 中国人民解放军61646部队 | 一种基于直接定位算法的星载sar时间参数标定方法 |
CN110865346A (zh) * | 2018-08-28 | 2020-03-06 | 中国人民解放军61646部队 | 一种基于直接定位算法的星载sar时间参数标定方法 |
CN110007302A (zh) * | 2019-03-27 | 2019-07-12 | 北京空间飞行器总体设计部 | 一种星载双天线斜视顺轨干涉sar洋流速度测量方法 |
CN110823191A (zh) * | 2019-10-08 | 2020-02-21 | 北京空间飞行器总体设计部 | 混合基线双天线斜视干涉sar洋流测量性能确定方法及系统 |
CN110823191B (zh) * | 2019-10-08 | 2021-12-07 | 北京空间飞行器总体设计部 | 混合基线双天线斜视干涉sar洋流测量性能确定方法及系统 |
CN110907932A (zh) * | 2019-11-26 | 2020-03-24 | 上海卫星工程研究所 | 分布式InSAR卫星测高精度影响因素分析方法及系统 |
CN111289957A (zh) * | 2020-03-10 | 2020-06-16 | 上海高仙自动化科技发展有限公司 | 外参标定方法及装置、智能机器人及计算机可读存储介质 |
CN111289957B (zh) * | 2020-03-10 | 2022-08-16 | 上海高仙自动化科技发展有限公司 | 外参标定方法及装置、智能机器人及计算机可读存储介质 |
CN112327261A (zh) * | 2020-10-22 | 2021-02-05 | 上海卫星工程研究所 | 分布式InSAR卫星时间同步在轨测试方法及系统 |
CN112327261B (zh) * | 2020-10-22 | 2022-10-25 | 上海卫星工程研究所 | 分布式InSAR卫星时间同步在轨测试方法及系统 |
CN115267780A (zh) * | 2022-09-26 | 2022-11-01 | 中国电子科技集团公司第十四研究所 | 一种同时实现顺轨和切轨干涉的多通道极化sar系统及其方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107102303B (zh) | 2024-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107102303A (zh) | 机载混合顺轨-交轨干涉合成孔径雷达系统的定标方法 | |
CN106289246B (zh) | 一种基于位置和姿态测量系统的柔性杆臂测量方法 | |
Huang et al. | Accurate 3-D position and orientation method for indoor mobile robot navigation based on photoelectric scanning | |
EP2828675B1 (en) | Method and device for pose tracking using vector magnetometers | |
CN103323855B (zh) | 一种基线动态测量系统的精度获取方法 | |
CN110907932B (zh) | 分布式InSAR卫星测高精度影响因素分析方法及系统 | |
CN103389099B (zh) | 基于x射线脉冲星的航天器姿态和位置测量系统及方法 | |
CN107270891B (zh) | 基于抗差估计的惯性地磁匹配定位方法 | |
Zhu et al. | Research on gravity vertical deflection on attitude of position and orientation system and compensation method | |
CN108106637A (zh) | 一种分布式pos的精度检校方法和装置 | |
CN106054185B (zh) | 一种基于分布式POS的机载双天线InSAR基线计算方法 | |
CN105203103B (zh) | 地面航天器相对地理方位关系的实时测量方法 | |
CN107991676A (zh) | 星载单航过InSAR系统对流层误差校正方法 | |
CN116540285B (zh) | 惯性辅助的gnss双天线定向方法、装置与电子设备 | |
Li et al. | Fast fine initial self-alignment of INS in erecting process on stationary base | |
Damiani et al. | A high-resolution dual-Doppler technique for fixed multiantenna airborne radar | |
CN103344958A (zh) | 基于星历数据的星载sar高阶多普勒参数估算方法 | |
CN106323271A (zh) | 基于特征奇异值的航天器相对姿态测量矢量选取方法 | |
CN104075710B (zh) | 一种基于航迹预测的机动扩展目标轴向姿态实时估计方法 | |
CN105260610B (zh) | 一种多探测器坐标系转化及误差纠正方法 | |
CN104330078B (zh) | 一种基于三点后方交会模型的联合测量方法 | |
Herman et al. | Nonlinear least-squares estimation for sensor and navigation biases | |
CN103886192B (zh) | 基于二次规划的运动目标单站定位方法 | |
Wu et al. | A sequential converted measurement Kalman filter with Doppler measurements in ECEF coordinate system | |
CN105242262B (zh) | 一种基于天线周期扫描时间差无源定位方法 |
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 |