CN107367716A - 一种高精度星载sar几何定标方法 - Google Patents

一种高精度星载sar几何定标方法 Download PDF

Info

Publication number
CN107367716A
CN107367716A CN201710539627.8A CN201710539627A CN107367716A CN 107367716 A CN107367716 A CN 107367716A CN 201710539627 A CN201710539627 A CN 201710539627A CN 107367716 A CN107367716 A CN 107367716A
Authority
CN
China
Prior art keywords
mrow
msub
point
mfrac
sar
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
Application number
CN201710539627.8A
Other languages
English (en)
Other versions
CN107367716B (zh
Inventor
张过
赵瑞山
邓明军
李少宁
徐凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201710539627.8A priority Critical patent/CN107367716B/zh
Publication of CN107367716A publication Critical patent/CN107367716A/zh
Application granted granted Critical
Publication of CN107367716B publication Critical patent/CN107367716B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • 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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及一种高精度星载SAR几何定标方法，将影响斜距测量精度的主要影响因素分为SAR系统时延误差和大气传播延迟误差。SAR系统时延误差主要是雷达信号在SAR载荷内部经过各个部件时产生的固有的时延误差，即从振荡器处产生到SAR天线发射端、再从SAR天线接收端到回波数据的接收采样，主要受雷达信号的带宽和脉冲宽度影响；大气传播延迟误差主要是雷达信号在传播路径上受到大气环境的折射等影响产生的时变的时延误差。由此，根据雷达信号的带宽和脉冲宽度组合，分组进行星载SAR的几何定标，无需考虑成像模式、波位、左右侧视等因素的影响，大大降低了星载SAR几何定标的工作量；考虑了时变的大气环境对斜距测量精度的影响，提高了定标精度。

Description

一种高精度星载SAR几何定标方法
技术领域
本发明属于微波遥感数据处理技术领域，尤其涉及一种高精度星载SAR几何定标方法。
背景技术
随着日本7m分辨率的ALOS‐PALSAR、加拿大lm分辨率的Radarsat‐2、意大利lm分辨率的Cosmo‐SkyMed、德国lm分辨率的军民两用卫星TerraSAR‐X和TanDEM‐X、日本1m分辨率的ALOS‐2、中国1m分辨率的高分三号等SAR卫星的相继成功发射，标志着高分辨率星载SAR时代的到来。随着SAR遥感影像的分辨率提高，对其几何定位能力的需求也相应提高。由于SAR卫星入轨后，环境变化、设备老化等因素严重影响其定位性能，因此，需要通过高精度的几何定标方法，提升SAR系统的无控定位能力。
星载SAR几何定标方法的核心是，利用控制点在SAR影像上的像素坐标与利用严密成像几何模型反算的像素坐标存在偏移量，根据偏移量再反演星载SAR系统的几何定标参数，进而实现提升星载SAR系统的无控定位精度。SAR信号的发射和接收都是在时间尺度上完成的，主要包括快时间(距离向)和慢时间(方位向)。这二维的时间误差主要影响SAR影像在距离向和方位向的几何定位误差，是星载SAR几何定位的主要误差源。经典的星载SAR几何定标主要是针对方位向起始时间和距离向起始斜距两个参数进行标定。由于系统设备时间控制单元的误差引起，系统设备时间控制单元用于记录所有事件的系统日期，它有一定的时间精度及一定的记录频率，它们会对方位向的开始时间精度带来一定影响。由于雷达信号经过信号通道的各个器件时产生的系统时延，不同类型(时宽和带宽组合)雷达信号的系统时延也不同；另外，雷达信号穿越大气到达地面再从地面返回大气时的大气路径双向延迟。因此，雷达信号的时宽和带宽以及大气传播延迟均影响着几何定标精度，而经典的几何定标方法均为充分考虑这两个影响因素。
发明内容
针对现有技术存在的不足，本发明顾及雷达信号的时宽和带宽特点和大气传播延迟影响，提出了一种高精度星载SAR几何定标方法，包括如下步骤：
步骤1，根据SAR系统设计的波位表，查找出与某一种时宽和带宽组合模式相同的波位号；
步骤2，通过对SAR卫星系统的成像任务规划，以查找出的任意波位对几何定标场区域进行成像；
步骤3，远程调整几何定标场区域的所有自动角反射器指向，使角反射器的法线方向与入射波的方向保持一致，根据自动角反射器调整的方位角和俯仰角，精确计算成像时刻各个角反射器顶点的地理坐标，同时，对各个角反射器顶点的地理坐标进行固体潮改正补偿；
步骤4，对几何定标场区域的回波数据进行成像处理，获得SAR影像和辅助参数文件；
步骤5，利用基于质心法的亚像元点提取方法，对每景SAR影像中的角反射器点进行精确提取，获得角反射器点的像方坐标；
步骤6，针对每景SAR影像，计算雷达信号在天线相位中心到各个角反射器点之间的大气传播延迟改正值，其中大气传播延迟包括中性大气延迟和电离层延迟；
步骤7，分别计算每景SAR影像的几何定标参数，所述几何定标参数包括距离向起始时间改正值、方位向起始时间改正值、采样频率改正值、脉冲重复频率改正值；
步骤8，根据几何定标场区域多景SAR影像的几何定标参数，取均值，获取该成像模式的几何定标参数。
进一步，还包括以下步骤，对几何定标后的几何精度验证场区域SAR影像进行精度评价，
步骤9，针对步骤1中的时宽和带宽组合模式，对几何精度验证场区域进行成像任务规划、拍摄；
步骤10，根据步骤8计算出的几何定标参数，对几何精度验证场区域的SAR回波数据进行补偿，并进行成像处理，获取SAR影像和辅助参数文件；
步骤11，获取几何精度验证场的地面控制点的像方坐标，并重复进行步骤6的操作，获得每个地面控制点的大气传播延迟改正值；
步骤12，根据SAR影像中地面控制点的像方坐标和大气传播延迟改正值，利用RD模型的直接定位方法，计算SAR影像中每个控制点的物方坐标，并与每个地面控制点的真实坐标进行对比，所得的差值即为该点的定位误差；统计所有地面控制点定位误差的标准差，即得到SAR系统在该成像模式下经几何定标后无控制点情况的几何定位精度。
进一步的，所述步骤3中根据自动角反射器调整的方位角和俯仰角，精确计算成像时刻各个角反射器顶点的地理坐标，同时，对各个角反射器顶点的地理坐标进行固体潮改正补偿，其实现方式如下，
假设以云台垂直方向的转动轴中心为坐标系原点，以真北方向为坐标系的X轴，建立左手直角坐标系O-XYZ，设P为角反射器顶点，P'为P在XOY平面的投影，P'A和P'B分别为P'到Y轴和X轴的距离，ω为点P的方位角，为点P的俯仰角；
根据高斯投影的方法，将地面控制点WGS‐84系下的大地坐标(Lat0,Lon0,H0)转换成高斯平面坐标(x0,y0,H0)，由此可知，坐标系原点坐标O为(x0,y0,H0+h)，点P的高斯平面坐标为其中，h为云台垂直方向的转动中心相对于地面控制点的高度，PO是角反射器顶点与云台垂直方向的转动中心之间的距离；根据高斯投影关系，反算出点P在WGS‐84系下的大地坐标(Latp,Lonp,Hp)，即获得高精度角反射器顶点的地理坐标；
采用以下公式对角反射器点的重力固体潮进行校正，
δg(t)＝δthG(t)-δfc (1)
式中，δth为潮汐因子，δfc是地心纬度的函数，即
式中，为地理纬度的函数，即Cm为地心至月心的平均距离，rm为地心至月心的距离，Cs为地心至日心的平均距离，rs为地心至日心的距离，Zm为月球对地面点的地心天顶距，Zs为太阳对地面点的地心天顶距。
进一步，所述步骤5的实现方式如下，
1)通过SAR图像间接定位模型，由地面大地坐标计算角反射器对应的SAR图像像点坐标，再由影像地面点坐标和轨道参数获得该地面点对应的方位向时间和角反射器点的星地距离后，解算该点对应的行列号；其中地面点坐标，由外业GPS测量获得，轨道参数由星上GPS获得，下传给地面使用；
2)在步骤1)中计算的行列号附近区域，寻找在此区域像素灰度值最大的亮点，再通过人工目视判读，比较周围地物，确定该点是否为角反射器在SAR图像上的粗略位置，初步确定角反射器的粗位置；
3)以步骤2)中的粗略位置所对应像素点为中心，选取其邻域为目标质心寻找区，将该领域内所有像点对应的坐标和灰度强度信息按顺序提取出来，并把该灰度赋予每个像素的几何中心；
4)对该目标区内每一像素进行的细分，细分尺寸根据不同的精度要求划分，然后采用双线性插值方法或Sinc函数插值方法求得每一细分点的灰度；
5)计算图像回波强度质心。
在计算质心时，采用逐行分析计算，取平均值法计算质心，计算方式如下：
计算每一行质心的坐标
其中，n为被插值点的个数，m为细分后点的个数，xk为被插值点的坐标，xu为细分后点的行坐标，Ik为被插值点的灰度值，f(u)为细分后点的灰度值；
然后，以每一行质心坐标和对应灰度值为参数，求整个目标区的质心，
其中，为每一行质心坐标，yu为细分后点的列坐标，为细分后的点对应平均灰度值；计算得到的(Xc,Yc)即为目标区域的质心，也即角反射器点对应的像方坐标。
进一步的，步骤6中所述的中性大气延迟改正值的计算方式如下，
上式中地表压强Psruf＝Pd+e，PW为大气可降水量；Md、MW分别为干、湿大气分子量，Md＝28.9644kg/kmol，MW＝18.0152kg/kmol；Rm为摩尔气体常量，Rm＝8.31451J/(mol·K)；
k1、k2值与卫星载荷发射电磁波信号的波长λ相关，求解公式：
其中，地表压强和大气可降水量从全球大气模型中获得。
进一步的，步骤6中所述的电离层延迟改正值的计算方式如下，
设目标点T垂直方向上的全部电子都压缩于一点T'处，雷达信号在传播路径上的电离层延迟为，
式中K＝40.28m3/s2，θ'是压缩点T'处的入射角，TEC的单位是1016个电子/米2，计算公式如下：
第一，在IONEX文件中取得全球的TEC含量图，然后以穿刺点的经纬度(lat,lon)为参数，从它附近的网格点做内插，具体的公式为，
式中，I0,0，I1,0，I1,1，I0,1为穿刺点相邻四角网格点的TEC含量，lat0为左下角网格点的纬度，lon0为左下角网格点的精度，Δlat为电离层网格点的纬度间隔，Δlon为电离层网格点的经度间隔；
第二，得到空间内插后进行相应时间点t的时间双线性内插，具体的公式为，
式中，Ti和Ti+1为时间点t两个临近时间格网的时间；
第三，取得穿刺点对应时刻垂直方向上的TEC后，根据公式(7)获得雷达信号在传播路径上的电离层延迟值。
进一步的，步骤7中每景SAR影像的几何定标参数计算方式如下，
构建星载SAR几何定标模型，
式中，tf、ts分别为距离向的快时间和方位向的慢时间，tf0、ts0分别为距离向起始时间的测量值和方位向起始时间测量值，tdelay为大气传播延迟时间，为大气传播延迟值除以光速，Δtf、Δts为系统时延误差，x、y为像素坐标，fs为采样频率，fp为脉冲重复频率；
对ts0进行补偿，
式中，ts'0为星上记录的方位向起始时间，N为雷达信号从发射到接收经历整周期的个数，tsample_delay为星上记录的采样时间延迟；
星载SAR几何定标模型(10)可以表示成如下形式：
式(12)的误差方程为
V＝Bx-l (13)
其中，
设初始值Δtf＝0、Δts＝0，fs和fp由星上下传的辅助参数文件获得，采用迭代运算的方法获取几何定标参数x＝[Δtf Δts Δfs Δfp]T，实现方式如下，
a)根据SAR影像间接定位算法，通过角反射器点的地面坐标和轨道参数，可以计算出该地面角反射器点对应的SAR影像方位向时间ts和该点的星地距离Rst，星地距离Rst即距离向时间tf；
b)根据公式(10)，利用几何定标场区域地面布设的n个角反射器点，组建n个方程组；通过基于谱修正法迭代法的最小二乘平差，精确计算方位向起始时间延迟Δts、距离向起始时间延迟Δtf、脉冲重复频率改正值Δfp＝fp'-fp和采样频率改正值Δfs＝fs'-fs；
c)更新四个几何定标参数，重新执行步骤a)，再次计算四个几何定标参数，判断两次计算的几何定标参数是否小于一个预先设定的极小值，若小于则迭代终止，否则转向a)继续迭代运算；
d)迭代运算结束后的结果为该景SAR影像的几何定标参数。
与现有技术相比，本发明的有益效果：
现有的星载SAR几何定标技术发展仍不成熟，尚无系统的定标方法和解决方案。本发明首先将影响斜距测量精度的主要影响因素进行了分类，主要分为SAR系统时延误差和大气传播延迟误差。SAR系统时延误差主要是雷达信号在SAR载荷内部经过各个部件时产生的固有的时延误差，即从振荡器处产生到SAR天线发射端、再从SAR天线接收端到回波数据的接收采样，主要受雷达信号的带宽和脉冲宽度影响；大气传播延迟误差主要是雷达信号在传播路径上受到大气环境的折射等影响产生的时变的时延误差。由此，根据雷达信号的带宽和脉冲宽度组合，分组进行星载SAR的几何定标，无需考虑成像模式、波位、左右侧视等因素的影响，大大降低了星载SAR几何定标的工作量；考虑了时变的大气环境对斜距测量精度的影响，提高了定标精度。
附图说明
图1为高精度星载SAR几何定标流程图；
图2为自动角反射器顶点坐标补偿示意图；
图3为单层电离层模型内插示意图。
具体实施方式
本发明提供了一种高精度星载SAR几何定标方法，其流程图如图1所示。以一种时宽和带宽组合的成像模式为例作为一种实施方式，具体实现步骤如下：
步骤1，根据SAR系统设计的波位表，查找出与该时宽和带宽组合模式相同的波位号。
步骤2，通过对SAR卫星系统的成像任务规划，以查找出的任意波位对几何定标场区域进行成像，至少成像3次；其中，几何定标场区域的SAR影像用于计算SAR卫星系统的几何定标参数。
步骤3，远程调整几何定标场的所有自动角反射器指向，使角反射器的法线方向与入射波的方向保持一致。根据自动角反射器调整的方位角和俯仰角，精确计算成像时刻各个角反射器顶点的地理坐标。同时，对各个角反射器顶点的地理坐标进行固体潮改正补偿。
由于自动角反射器设备在水平和垂直方向上转动，其顶点位置会随着方位角和俯仰角的变化而变化。为了获取角反射器顶点的精确地面坐标，需要对角反射器顶点的坐标进行补偿修正。假设以云台垂直方向的转动轴中心为坐标系原点，以真北方向为坐标系的X轴，建立左手直角坐标系O-XYZ，如图2所示。其中，P为角反射器顶点，P'为P在XOY平面的投影，P'A和P'B分别为P'到Y轴和X轴的距离，ω为点P的方位角，为点P的俯仰角。
根据高斯投影的方法，将地面控制点WGS‐84系下的大地坐标(Lat0,Lon0,H0)转换成高斯平面坐标(x0,y0,H0)。由此可知，坐标系原点坐标O为(x0,y0,H0+h)，点P的高斯平面坐标为其中，h为云台垂直方向的转动中心相对于地面控制点的高度，PO是角反射器顶点与云台垂直方向的转动中心之间的距离。再根据高斯投影关系，反算出点P在WGS‐84系下的大地坐标(Latp,Lonp,Hp)，即获得高精度角反射器顶点的地理坐标。
由于地球自传，地球上任意点相对于月球和太阳的位置发生周期地改变，故重力固体潮随时间而周期变化。在进行角反射器点的重力固体潮校正时，一般采用下面公式计算：
δg(t)＝δthG(t)-δfc (1)
式中，δth为潮汐因子，一般采用1.15或1.16；δfc是地心纬度的函数，即
式中，为地理纬度的函数，即Cm为地心至月心的平均距离，rm为地心至月心的距离，Cs为地心至日心的平均距离，rs为地心至日心的距离，Zm为月球对地面点的地心天顶距，Zs为太阳对地面点的地心天顶距。
步骤4，对几何定标场区域的回波数据进行成像处理，获得SAR影像和辅助参数文件。
步骤5，利用基于质心法的亚像元点提取方法，对每景SAR影像中的角反射器点进行精确提取，获得角反射器点的像方坐标。
亚像元方法是利用目标像点的灰度分布特性通过内插细分算法确定出目标像点位置，从而使测量精度可以达到亚像元级。它是通过构造一个尽可能精确反映目标区域像素灰度值位置和目标像点的质心位置之间关系的数学模型，从而实现对像点位置的精确估计。为了提高求得的质心位置的精度，需要将读入的图像进行细分，通过对图像各点的灰度值的插值，得到插值曲线，求得细分后各点的灰度值。根据不同的精度要求，图像细分比例可取0.1、0.01、0.001等。
其算法具体步骤如下：
1)通过SAR图像间接定位模型，由地面大地坐标计算角反射器对应的SAR图像像点坐标，再由影像地面点坐标和轨道参数获得该地面点对应的方位向时间和角反射器点的星地距离后，解算该点对应的行列号；其中地面点坐标，由外业GPS测量获得；轨道参数由星上GPS获得，下传给地面使用。
2)在步骤1)中计算的行列号附近区域(如10*10像素)，通过计算机寻找在此区域像素灰度值最大的亮点，再通过人工目视判读，比较周围地物，确定该点是否为角反射器在SAR图像上的粗略位置，初步确定角反射器的粗位置。
3)以步骤2)中的粗略位置所对应像素点为中心，选取其邻域(例如选取5×5领域)为目标质心寻找区。将该领域内所有像点对应的坐标和灰度强度信息按顺序提取出来，并把该灰度赋予每个像素的几何中心，例如对像素坐标点(8000，5000)，其中心(7999.5，4999.5)的灰度为该像素的灰度大小。
4)对该目标区内每一像素进行的细分，细分尺寸根据不同的精度要求可以是10×10或100×100或1000×1000，当细分尺寸为10×10时，插值点步长为0.1，它的精度即可达到0.1个像素。采用一定的数学函数方法(如双线性插值方法、Sinc函数插值方法等)求得每一细分点的灰度。
5)计算图像回波强度质心。
在计算质心时，采用逐行分析计算，取平均值法计算质心，其实现原理为：
对目标区域，逐行对像素图像(行列都以细分后的单位为基准，如细分尺寸为100×100，则行列基本单位为原图像的0.01，即0.01行/列)进行处理(包括插值、图像分割、质心计算等)，并求出每一行的质心坐标和平均灰度值，质心坐标也就是X轴方向坐标，最后将每一行的(取决于图像的高度)质心坐标求平均值，即得原图像质心位移。
计算每一行质心的坐标x_u公式为：
其中，n为被插值点的个数，m为细分后点的个数，xk为被插值点的坐标，xu为细分后点的行坐标，Ik为被插值点的灰度值，f(u)为细分后点的灰度值。
然后，以每一行质心坐标和对应灰度值为参数，求整个目标区的质心公式为：
其中，x_u为每一行质心坐标，yu为细分后点的列坐标，f__(_u_)为细分后的点对应平均灰度值。
计算得到的(Xc,Yc)即为目标区域的质心，也即角反射器点对应的像方坐标。计算质心时，根据插值点步长(如有效步长为0.001)的选取和像素的计算精度要求，精度取到0.0001即可。
步骤6，针对每景SAR影像，计算雷达信号在天线相位中心到各个角反射器点之间的大气传播延迟改正值，大气传播延迟主要分为中性大气延迟和电离层延迟。
James C.Owens结合流体静力学方程和非理想气体公式推导出中性大气中干分量和湿分量的距离改正模型，如下式：
上式中地表压强Psruf＝Pd+e，PW为大气可降水量；Md、MW分别为干、湿大气分子量，Md＝28.9644kg/kmol，MW＝18.0152kg/kmol；Rm为摩尔气体常量，Rm＝8.31451J/(mol·K)。
k1、k2值与卫星载荷发射电磁波信号的波长λ相关，求解公式：
由式(5)可知，干大气和湿大气的延迟量求解还需要获取大气的地表压强和大气可降水量，所述的地表压强和大气可降水量通常从全球大气模型中获得。全球大气模型(Global Atmospheric Models,GAM)可以提供由卫星和气象基站记录的地表气象数据，包括温度、气压、水汽含量、风速风向等。本发明以美国国家环境预报中心(NCEP)的大气分析模型作为外部数据，模型提供了每隔6小时1°×1°的经纬网格点存储的等压面数据。
图3为单层电离层模型示意图，可看作把目标点T垂直方向上的全部电子都压缩于一点T'处，电离层延迟是压缩点T'处的TEC含量和视向与压缩点T'处单层电离层法线方向的夹角有很大关系。由此，可以得到雷达信号在传播路径上的电离层延迟：
式中K＝40.28m3/s2，TEC的单位是1016个电子/米2，θ'是压缩点T'处的入射角。由上式可知，电离层延迟与TEC和电磁波频率fc有关(已知)。
欧洲定轨中心(CODE)以IONEX(Ionosphere Exchange Format)文件格式每天给出，一天中从UTC零时到24时，每过2个小时会产生一幅全球电离层TEC图，每天给出13幅(从2014年10月19日开始，以1小时为间隔每天给出25幅)，IGS计算电离层的格网划分按照全球经度方向间隔5°，纬度方向间隔2.5°，经度范围为‐180°W至180°E，纬度范围为87.5°S至‐87.5°N，共有5183个格网点。如果要获得某一特定时间穿刺点的延迟，要经过以下步骤：首先对历元时刻的TEC含量采用双线性内插原理进行四站网格空间内插；接着按实际的时间对历元之间做内插，具体过程如下：
第一，在IONEX文件中取得全球的TEC含量图，然后以穿刺点的经纬度(lat,lon)为参数，从它附近的网格点做内插，具体的公式为：
式中，I0,0，I1,0，I1,1，I0,1为穿刺点相邻四角网格点的TEC含量，lat0为左下角网格点的纬度，lon0为左下角网格点的精度，Δlat为电离层网格点的纬度间隔，Δlon为电离层网格点的经度间隔。
第二，得到空间内插后进行相应时间点t的时间双线性内插，具体的公式是：
式中，Ti和Ti+1为时间点t两个临近时间格网的时间。
第三，取得穿刺点对应时刻垂直方向上的TEC后，根据公式(7)就能获得雷达信号在传播路径上的电离层延迟值。
步骤7，分别计算每景SAR影像的几何定标参数，即距离向起始时间改正值Δtf、方位向起始时间改正值Δts、采样频率改正值Δf、脉冲重复频率改正值Δfp。
SAR信号的发射和接收都是在时间尺度上完成的，主要包括快时间(距离向)和慢时间(方位向)。这二维的时间误差主要影响SAR影像在距离向和方位向的几何定位误差，是星载SAR几何定位的主要误差源。由此，构建星载SAR几何定标模型：
式中，tf、ts分别为距离向的快时间和方位向的慢时间，tf0、ts0分别为距离向起始时间的测量值和方位向起始时间测量值，tdelay为大气传播延迟时间，为大气传播延迟值除以光速，Δtf、Δts为系统时延误差，x、y为像素坐标，fs为采样频率，fp为脉冲重复频率(PRF)。
星上记录的方位向起始时间ts'0是雷达信号的接收时刻，然而SAR卫星是在持续运动中发射和接收雷达信号，等效的SAR成像时刻ts0为雷达信号的发射时刻与接收时刻的中间时刻，故需进行补偿：
式中，N为雷达信号从发射到接收经历整周期的个数，tsample_delay为星上记录的采样时间延迟。
星载SAR几何定标模型(10)可以表示成如下形式：
式(12)的误差方程为
V＝Bx-l (13)
其中，
设初始值Δtf＝0、Δts＝0，fs和fp由星上下传的辅助参数文件获得，采用迭代运算的方法获取几何定标参数x＝[Δtf Δts Δfs Δfp]T。几何定标参数解算的具体步骤如下：
1)根据SAR影像间接定位算法，通过角反射器点的地面坐标和轨道参数，可以计算出该地面角反射器点对应的SAR影像方位向时间ts和该点的星地距离Rst(即距离向时间tf)。
2)根据公式(10)，利用定标场地面布设的n个角反射器点，可以组建n个方程组。通过基于谱修正法迭代法的最小二乘平差，可精确计算方位向起始时间改正值Δts、距离向起始时间改正值Δtf、PRF改正值Δfp＝fp'-fp和采样频率改正值Δfs＝fs'-fs。
3)更新四个几何定标参数，重新执行步骤1)，再次计算四个几何定标参数，判断两次计算的几何定标参数是否小于一个预先设定的极小值，若小于则迭代终止，否则转向1)继续迭代运算。
4)迭代运算结束后的结果，就是该景SAR影像的几何定标参数。
步骤8，根据几何定标场区域多景SAR影像的几何定标参数，取均值，获取该成像模式的几何定标参数。
步骤9，针对该时宽和带宽组合模式，对几何精度验证场区域进行成像任务规划、拍摄，至少2个精度验证区域；其中，几何精度验证场区域的SAR影像用于对上述步骤获取的几何定标参数的正确性进行验证。
步骤10，根据步骤8计算出的几何定标参数，对几何精度验证场区域的SAR回波数据进行补偿，并进行成像处理，获取SAR影像和辅助参数文件。
步骤11，获取几何精度验证场的地面控制点(角反射器点、GNSS测量点等)的像方坐标，并重复进行步骤6的操作，获得每个地面控制点的大气传播延迟改正值。
步骤12，对几何定标后的几何精度验证场区域SAR影像进行精度评价。根据SAR影像中地面控制点的像方坐标和大气传播延迟改正值，利用RD模型的直接定位方法，计算SAR影像中每个控制点的物方坐标，并与每个地面控制点的真实坐标进行对比，所得的差值即为该点的定位误差。统计所有地面控制点定位误差的标准差，即得到SAR系统在该成像模式下经几何定标后无控制点情况的几何定位精度。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

Claims (7)

1.一种高精度星载SAR几何定标方法，其特征在于，包括如下步骤：
步骤1，根据SAR系统设计的波位表，查找出与某一种时宽和带宽组合模式相同的波位号；
步骤2，通过对SAR卫星系统的成像任务规划，以查找出的任意波位对几何定标场区域进行成像；
步骤3，远程调整几何定标场区域的所有自动角反射器指向，使角反射器的法线方向与入射波的方向保持一致，根据自动角反射器调整的方位角和俯仰角，精确计算成像时刻各个角反射器顶点的地理坐标，同时，对各个角反射器顶点的地理坐标进行固体潮改正补偿；
步骤4，对几何定标场区域的回波数据进行成像处理，获得SAR影像和辅助参数文件；
步骤5，利用基于质心法的亚像元点提取方法，对每景SAR影像中的角反射器点进行精确提取，获得角反射器点的像方坐标；
步骤6，针对每景SAR影像，计算雷达信号在天线相位中心到各个角反射器点之间的大气传播延迟改正值，其中大气传播延迟包括中性大气延迟和电离层延迟；
步骤7，分别计算每景SAR影像的几何定标参数，所述几何定标参数包括距离向起始时间改正值、方位向起始时间改正值、采样频率改正值、脉冲重复频率改正值；
步骤8，根据几何定标场区域多景SAR影像的几何定标参数，取均值，获取该成像模式的几何定标参数。
2.如权利要求1所述的一种高精度星载SAR几何定标方法，其特征在于：还包括以下步骤，对几何定标后的几何精度验证场区域SAR影像进行精度评价，
步骤9，针对步骤1中的时宽和带宽组合模式，对几何精度验证场区域进行成像任务规划、拍摄；
步骤10，根据步骤8计算出的几何定标参数，对几何精度验证场区域的SAR回波数据进行补偿，并进行成像处理，获取SAR影像和辅助参数文件；
步骤11，获取几何精度验证场的地面控制点的像方坐标，并重复进行步骤6的操作，获得每个地面控制点的大气传播延迟改正值；
步骤12，根据SAR影像中地面控制点的像方坐标和大气传播延迟改正值，利用RD模型的直接定位方法，计算SAR影像中每个控制点的物方坐标，并与每个地面控制点的真实坐标进行对比，所得的差值即为该点的定位误差；统计所有地面控制点定位误差的标准差，即得到SAR系统在该成像模式下经几何定标后无控制点情况的几何定位精度。
3.如权利要求1或2所述的一种高精度星载SAR几何定标方法，其特征在于：所述步骤3中根据自动角反射器调整的方位角和俯仰角，精确计算成像时刻各个角反射器顶点的地理坐标，同时，对各个角反射器顶点的地理坐标进行固体潮改正补偿，其实现方式如下，
假设以云台垂直方向的转动轴中心为坐标系原点，以真北方向为坐标系的X轴，建立左手直角坐标系O-XYZ，设P为角反射器顶点，P′为P在XOY平面的投影，P′A和P′B分别为P′到Y轴和X轴的距离，ω为点P的方位角，为点P的俯仰角；
根据高斯投影的方法，将地面控制点WGS-84系下的大地坐标(Lat0，Lon0，H0)转换成高斯平面坐标(x0，y0，H0)，由此可知，坐标系原点坐标O为(x0，y0，H0+h)，点P的高斯平面坐标为其中，h为云台垂直方向的转动中心相对于地面控制点的高度，PO是角反射器顶点与云台垂直方向的转动中心之间的距离；根据高斯投影关系，反算出点P在WGS-84系下的大地坐标(Latp，Lonp，Hp)，即获得高精度角反射器顶点的地理坐标；
采用以下公式对角反射器点的重力固体潮进行校正，
δg(t)＝δthG(t)-δfc (1)
式中，δth为潮汐因子，δfc是地心纬度的函数，即
式中，为地理纬度的函数，即Cm为地心至月心的平均距离，rm为地心至月心的距离，Cs为地心至日心的平均距离，rs为地心至日心的距离，Zm为月球对地面点的地心天顶距，Zs为太阳对地面点的地心天顶距。
4.如权利要求3所述的一种高精度星载SAR几何定标方法，其特征在于：所述步骤5的实现方式如下，
1)通过SAR图像间接定位模型，由地面大地坐标计算角反射器对应的SAR图像像点坐标，再由影像地面点坐标和轨道参数获得该地面点对应的方位向时间和角反射器点的星地距离后，解算该点对应的行列号；其中地面点坐标，由外业GPS测量获得，轨道参数由星上GPS获得，下传给地面使用；
2)在步骤1)中计算的行列号附近区域，寻找在此区域像素灰度值最大的亮点，再通过人工目视判读，比较周围地物，确定该点是否为角反射器在SAR图像上的粗略位置，初步确定角反射器的粗位置；
3)以步骤2)中的粗略位置所对应像素点为中心，选取其邻域为目标质心寻找区，将该领域内所有像点对应的坐标和灰度强度信息按顺序提取出来，并把该灰度赋予每个像素的几何中心；
4)对该目标区内每一像素进行的细分，细分尺寸根据不同的精度要求划分，然后采用双线性插值方法或Sinc函数插值方法求得每一细分点的灰度；
5)计算图像回波强度质心。
在计算质心时，采用逐行分析计算，取平均值法计算质心，计算方式如下：
计算每一行质心的坐标
<mrow> <mover> <msub> <mi>x</mi> <mi>u</mi> </msub> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>x</mi> <mi>k</mi> </msub> <msub> <mi>I</mi> <mi>k</mi> </msub> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>u</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>x</mi> <mi>u</mi> </msub> <mi>f</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>I</mi> <mi>k</mi> </msub> <mo>+</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>u</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <mi>f</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
其中，n为被插值点的个数，m为细分后点的个数，xk为被插值点的坐标，xu为细分后点的行坐标，Ik为被插值点的灰度值，f(u)为细分后点的灰度值；
然后，以每一行质心坐标和对应灰度值为参数，求整个目标区的质心，
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>X</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>u</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <mover> <msub> <mi>x</mi> <mi>u</mi> </msub> <mo>&amp;OverBar;</mo> </mover> <mover> <mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </mrow> <mo>&amp;OverBar;</mo> </mover> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>u</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <mover> <mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </mrow> <mo>&amp;OverBar;</mo> </mover> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>Y</mi> <mi>c</mi> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>u</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>y</mi> <mi>u</mi> </msub> <mover> <mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </mrow> <mo>&amp;OverBar;</mo> </mover> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>u</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <mover> <mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> </mrow> <mo>&amp;OverBar;</mo> </mover> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
其中，为每一行质心坐标，yu为细分后点的列坐标，为细分后的点对应平均灰度值；计算得到的(Xc,Yc)即为目标区域的质心，也即角反射器点对应的像方坐标。
5.如权利要求4所述的一种高精度星载SAR几何定标方法，其特征在于：步骤6中所述的中性大气延迟改正值的计算方式如下，
上式中地表压强Psruf＝Pd+e，PW为大气可降水量；Md、MW分别为干、湿大气分子量，Md＝28.9644kg/kmol，MW＝18.0152kg/kmol；Rm为摩尔气体常量，Rm＝8.31451J/(mol·K)；
k1、k2值与卫星载荷发射电磁波信号的波长λ相关，求解公式：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>0.237134</mn> <mo>+</mo> <mn>68.39397</mn> <mfrac> <mrow> <mo>(</mo> <mn>130</mn> <mo>+</mo> <msup> <mi>&amp;lambda;</mi> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </msup> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mn>130</mn> <mo>-</mo> <msup> <mi>&amp;lambda;</mi> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </msup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mfrac> <mo>+</mo> <mn>0.45473</mn> <mfrac> <mrow> <mo>(</mo> <mn>38.9</mn> <mo>+</mo> <msup> <mi>&amp;lambda;</mi> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </msup> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mn>38.9</mn> <mo>-</mo> <msup> <mi>&amp;lambda;</mi> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </msup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>k</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>0.648731</mn> <mo>+</mo> <mn>0.0174174</mn> <msup> <mi>&amp;lambda;</mi> <mrow> <mo>-</mo> <mn>2</mn> </mrow> </msup> <mo>+</mo> <mn>3.55750</mn> <mo>&amp;times;</mo> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mn>4</mn> </mrow> </msup> <msup> <mi>&amp;lambda;</mi> <mrow> <mo>-</mo> <mn>4</mn> </mrow> </msup> <mo>+</mo> <mn>6.1957</mn> <mo>&amp;times;</mo> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mn>5</mn> </mrow> </msup> <msup> <mi>&amp;lambda;</mi> <mrow> <mo>-</mo> <mn>6</mn> </mrow> </msup> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
其中，地表压强和大气可降水量从全球大气模型中获得。
6.如权利要求5所述的一种高精度星载SAR几何定标方法，其特征在于：步骤6中所述的电离层延迟改正值的计算方式如下，
设目标点T垂直方向上的全部电子都压缩于一点T'处，雷达信号在传播路径上的电离层延迟为，
<mrow> <msub> <mi>&amp;Delta;L</mi> <mi>I</mi> </msub> <mo>=</mo> <mi>K</mi> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <mi>T</mi> <mi>E</mi> <mi>C</mi> </mrow> <msubsup> <mi>f</mi> <mi>c</mi> <mn>2</mn> </msubsup> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mn>1</mn> <mrow> <msup> <mi>cos&amp;theta;</mi> <mo>&amp;prime;</mo> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
式中K＝40.28m3/s2，θ'是压缩点T'处的入射角，TEC的单位是1016个电子/米2，计算公式如下：
第一，在IONEX文件中取得全球的TEC含量图，然后以穿刺点的经纬度(lat,lon)为参数，从它附近的网格点做内插，具体的公式为，
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>I</mi> <mrow> <mo>(</mo> <mrow> <mi>l</mi> <mi>a</mi> <mi>t</mi> <mo>,</mo> <mi>l</mi> <mi>o</mi> <mi>n</mi> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>I</mi> <mrow> <mn>0</mn> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>&amp;times;</mo> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mi>l</mi> <mi>a</mi> <mi>t</mi> <mo>-</mo> <msub> <mi>lat</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>l</mi> <mi>a</mi> <mi>t</mi> </mrow> </mfrac> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mi>l</mi> <mi>o</mi> <mi>n</mi> <mo>-</mo> <msub> <mi>lon</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>l</mi> <mi>o</mi> <mi>n</mi> </mrow> </mfrac> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>I</mi> <mrow> <mn>1</mn> <mo>,</mo> <mn>0</mn> </mrow> </msub> <mo>&amp;times;</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>l</mi> <mi>a</mi> <mi>t</mi> <mo>-</mo> <msub> <mi>lat</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>l</mi> <mi>a</mi> <mi>t</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mi>l</mi> <mi>o</mi> <mi>n</mi> <mo>-</mo> <msub> <mi>lon</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>l</mi> <mi>o</mi> <mi>n</mi> </mrow> </mfrac> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <msub> <mi>I</mi> <mrow> <mn>1</mn> <mo>,</mo> <mn>1</mn> </mrow> </msub> <mo>&amp;times;</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>l</mi> <mi>a</mi> <mi>t</mi> <mo>-</mo> <msub> <mi>lat</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>l</mi> <mi>a</mi> <mi>t</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>l</mi> <mi>o</mi> <mi>n</mi> <mo>-</mo> <msub> <mi>lon</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>l</mi> <mi>o</mi> <mi>n</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>I</mi> <mrow> <mn>0</mn> <mo>,</mo> <mn>1</mn> </mrow> </msub> <mo>&amp;times;</mo> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mi>l</mi> <mi>a</mi> <mi>t</mi> <mo>-</mo> <msub> <mi>lat</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>l</mi> <mi>a</mi> <mi>t</mi> </mrow> </mfrac> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>l</mi> <mi>o</mi> <mi>n</mi> <mo>-</mo> <msub> <mi>lon</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>l</mi> <mi>o</mi> <mi>n</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
式中，I0,0，I1,0，I1,1，I0,1为穿刺点相邻四角网格点的TEC含量，lat0为左下角网格点的纬度，lon0为左下角网格点的精度，Δlat为电离层网格点的纬度间隔，Δlon为电离层网格点的经度间隔；
第二，得到空间内插后进行相应时间点t的时间双线性内插，具体的公式为，
<mrow> <mi>I</mi> <mrow> <mo>(</mo> <mi>l</mi> <mi>a</mi> <mi>t</mi> <mo>,</mo> <mi>l</mi> <mi>o</mi> <mi>n</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>I</mi> <msub> <mrow> <mo>(</mo> <mi>l</mi> <mi>a</mi> <mi>t</mi> <mo>,</mo> <mi>l</mi> <mi>o</mi> <mi>n</mi> <mo>)</mo> </mrow> <mi>i</mi> </msub> <mo>&amp;times;</mo> <mfrac> <mrow> <msub> <mi>T</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <mi>t</mi> </mrow> <mrow> <msub> <mi>T</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>+</mo> <mi>I</mi> <msub> <mrow> <mo>(</mo> <mi>l</mi> <mi>a</mi> <mi>t</mi> <mo>,</mo> <mi>l</mi> <mi>o</mi> <mi>n</mi> <mo>)</mo> </mrow> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>&amp;times;</mo> <mfrac> <mrow> <mi>t</mi> <mo>-</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> <mrow> <msub> <mi>T</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
式中，Ti和Ti+1为时间点t两个临近时间格网的时间；
第三，取得穿刺点对应时刻垂直方向上的TEC后，根据公式(7)获得雷达信号在传播路径上的电离层延迟值。
7.如权利要求6所述的一种高精度星载SAR几何定标方法，其特征在于：步骤7中每景SAR影像的几何定标参数计算方式如下，
构建星载SAR几何定标模型，
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>t</mi> <mi>f</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mrow> <mi>f</mi> <mn>0</mn> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mrow> <mi>d</mi> <mi>e</mi> <mi>l</mi> <mi>a</mi> <mi>y</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;Delta;t</mi> <mi>f</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mi>x</mi> <mo>-</mo> <mn>1</mn> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>t</mi> <mi>s</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>t</mi> <mrow> <mi>s</mi> <mn>0</mn> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;Delta;t</mi> <mi>s</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mi>y</mi> <mo>-</mo> <mn>1</mn> </mrow> <msub> <mi>f</mi> <mi>p</mi> </msub> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
式中，tf、ts分别为距离向的快时间和方位向的慢时间，tf0、ts0分别为距离向起始时间的测量值和方位向起始时间测量值，tdelay为大气传播延迟时间，为大气传播延迟值除以光速，Δtf、Δts为系统时延误差，x、y为像素坐标，fs为采样频率，fp为脉冲重复频率；
对ts0进行补偿，
<mrow> <msub> <mi>t</mi> <mrow> <mi>s</mi> <mn>0</mn> </mrow> </msub> <mo>=</mo> <msubsup> <mi>t</mi> <mrow> <mi>s</mi> <mn>0</mn> </mrow> <mo>&amp;prime;</mo> </msubsup> <mo>-</mo> <mfrac> <mrow> <mi>N</mi> <mo>/</mo> <msub> <mi>f</mi> <mi>p</mi> </msub> <mo>+</mo> <msub> <mi>t</mi> <mrow> <mi>s</mi> <mi>a</mi> <mi>m</mi> <mi>p</mi> <mi>l</mi> <mi>e</mi> <mo>_</mo> <mi>d</mi> <mi>e</mi> <mi>l</mi> <mi>a</mi> <mi>y</mi> </mrow> </msub> </mrow> <mn>2</mn> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
式中，t′s0为星上记录的方位向起始时间，N为雷达信号从发射到接收经历整周期的个数，tsample_delay为星上记录的采样时间延迟；
星载SAR几何定标模型(10)可以表示成如下形式：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>F</mi> <mi>x</mi> </msub> <mo>=</mo> <msub> <mi>t</mi> <mi>f</mi> </msub> <mo>-</mo> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>t</mi> <mrow> <mi>f</mi> <mn>0</mn> </mrow> </msub> <mo>+</mo> <msub> <mi>t</mi> <mrow> <mi>d</mi> <mi>e</mi> <mi>l</mi> <mi>a</mi> <mi>y</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;Delta;t</mi> <mi>f</mi> </msub> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mi>x</mi> <mo>-</mo> <mn>1</mn> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>F</mi> <mi>y</mi> </msub> <mo>=</mo> <msub> <mi>t</mi> <mi>s</mi> </msub> <mo>-</mo> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>t</mi> <mrow> <mi>s</mi> <mn>0</mn> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;Delta;t</mi> <mi>s</mi> </msub> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mrow> <mi>y</mi> <mo>-</mo> <mn>1</mn> </mrow> <msub> <mi>f</mi> <mi>p</mi> </msub> </mfrac> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
式(12)的误差方程为
V＝Bx-l (13)
其中，x＝[Δtf Δts fs' f′p]T，
设初始值Δtf＝0、Δts＝0，fs和fp由星上下传的辅助参数文件获得，采用迭代运算的方法获取几何定标参数x＝[Δtf Δts Δfs Δfp]T，实现方式如下，
a)根据SAR影像间接定位算法，通过角反射器点的地面坐标和轨道参数，可以计算出该地面角反射器点对应的SAR影像方位向时间ts和该点的星地距离Rst，星地距离Rst即距离向时间tf；
b)根据公式(10)，利用几何定标场区域地面布设的n个角反射器点，组建n个方程组；通过基于谱修正法迭代法的最小二乘平差，精确计算方位向起始时间延迟Δts、距离向起始时间延迟Δtf、脉冲重复频率改正值Δfp＝f′p-fp和采样频率改正值Δfs＝f′s-fs；
c)更新四个几何定标参数，重新执行步骤a)，再次计算四个几何定标参数，判断两次计算的几何定标参数是否小于一个预先设定的极小值，若小于则迭代终止，否则转向a)继续迭代运算；
d)迭代运算结束后的结果为该景SAR影像的几何定标参数。
CN201710539627.8A 2017-07-04 2017-07-04 一种高精度星载sar几何定标方法 Active CN107367716B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710539627.8A CN107367716B (zh) 2017-07-04 2017-07-04 一种高精度星载sar几何定标方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710539627.8A CN107367716B (zh) 2017-07-04 2017-07-04 一种高精度星载sar几何定标方法

Publications (2)

Publication Number Publication Date
CN107367716A true CN107367716A (zh) 2017-11-21
CN107367716B CN107367716B (zh) 2019-09-10

Family

ID=60306189

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710539627.8A Active CN107367716B (zh) 2017-07-04 2017-07-04 一种高精度星载sar几何定标方法

Country Status (1)

Country Link
CN (1) CN107367716B (zh)

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108469614A (zh) * 2018-03-27 2018-08-31 武汉大学 基于未定标极化雷达影像的角反射器探测方法
CN108562882A (zh) * 2018-06-21 2018-09-21 武汉大学 一种星载sar影像几何交叉定标方法和系统
CN109581367A (zh) * 2018-12-19 2019-04-05 西安空间无线电技术研究所 一种天基预警雷达重频组优化设计方法
CN109696182A (zh) * 2019-01-23 2019-04-30 张过 一种星载推扫式光学传感器内方位元素定标方法
CN109782280A (zh) * 2019-03-06 2019-05-21 中国科学院电子学研究所 基于定标卫星的中高轨sar天线距离向方向图测量方法
CN109870685A (zh) * 2017-12-05 2019-06-11 南京理工大学 基于改进的rd算法的室内距离向移动sar成像方法
CN110031841A (zh) * 2019-04-01 2019-07-19 中国科学院遥感与数字地球研究所 基于ECMWF的InSAR大气延迟改正的方法和系统
CN110531329A (zh) * 2019-09-25 2019-12-03 中南大学 一种用于冲激脉冲式穿墙雷达的多通道误差自校正方法
CN110703247A (zh) * 2019-10-30 2020-01-17 中国科学院电子学研究所 一种波位信息检测方法、波位信息检测装置和存储介质
CN110865346A (zh) * 2018-08-28 2020-03-06 中国人民解放军61646部队 一种基于直接定位算法的星载sar时间参数标定方法
CN111596272A (zh) * 2020-06-01 2020-08-28 中国科学院空天信息创新研究院 基于角反射器的机载多极化sar系统外定标方法
CN111896954A (zh) * 2020-08-06 2020-11-06 华能澜沧江水电股份有限公司 一种船载sar影像的角反射器坐标定位方法
CN111915570A (zh) * 2020-07-10 2020-11-10 云南电网有限责任公司带电作业分公司 基于反向传播神经网络的大气延迟估计方法
CN112083387A (zh) * 2020-09-02 2020-12-15 南京慧尔视智能科技有限公司 一种雷达标定方法及装置
CN112179850A (zh) * 2020-11-09 2021-01-05 南京信息工程大学 基于地面实测的gf-4 aod产品高时间分辨率检验
CN112764052A (zh) * 2020-12-25 2021-05-07 中国人民解放军32181部队 防空导弹飞行监测系统
CN113093190A (zh) * 2021-04-08 2021-07-09 中国电子科技集团公司第三十八研究所 基于高精度组合惯导系统的机载条带sar图像定位方法
TWI735390B (zh) * 2020-12-01 2021-08-01 財團法人金屬工業研究發展中心 影像定位系統之即時定位補償方法及可即時定位補償之影像定位系統
CN113446902A (zh) * 2021-07-07 2021-09-28 中国人民解放军32382部队 一种自动对消测角误差的防空导弹的外弹道监测系统
CN113671505A (zh) * 2021-09-08 2021-11-19 辽宁工程技术大学 一种基于系统几何误差补偿的合成孔径雷达立体定位方法
CN113742803A (zh) * 2021-09-07 2021-12-03 辽宁工程技术大学 一种中高轨道sar卫星带控几何定位精度仿真分析方法
CN113960545A (zh) * 2021-10-19 2022-01-21 安徽大学 基于对称几何构型约束的星载sar无场几何定标方法及其系统
CN114065531A (zh) * 2021-11-20 2022-02-18 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种实测数据驱动的电离层法拉第旋转效应估计方法
CN114594435A (zh) * 2022-01-18 2022-06-07 中国资源卫星应用中心 一种国产民用sar卫星几何检校与定位精度改进方法
CN115629363A (zh) * 2022-12-16 2023-01-20 西安空间无线电技术研究所 星载探测雷达实时内定标方法及装置
CN115856807A (zh) * 2023-02-22 2023-03-28 中国科学院空天信息创新研究院 高分辨率sar卫星定标器雷达波相位中心精确定位方法
CN117665818A (zh) * 2024-02-02 2024-03-08 北京东方至远科技股份有限公司 一种针对合成孔径雷达卫星的平面位置修正方法及系统
CN117761684A (zh) * 2023-11-15 2024-03-26 中国科学院空间应用工程与技术中心 一种天基干涉成像雷达高度计方位时延确定方法及系统
CN118091567A (zh) * 2024-04-28 2024-05-28 中山大学 基于正射影像能量的sar快时间慢时间参数虚拟定标方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103247029A (zh) * 2013-03-26 2013-08-14 中国科学院上海技术物理研究所 一种用于拼接式探测器生成的高光谱图像几何配准方法
CN105403886A (zh) * 2015-12-02 2016-03-16 中国科学院电子学研究所 一种机载sar定标器图像位置自动提取方法
CN106249766A (zh) * 2016-09-09 2016-12-21 中国科学院电子学研究所 远程监控角反射器系统
CN106526593A (zh) * 2016-12-19 2017-03-22 国家测绘地理信息局卫星测绘应用中心 基于sar严密成像模型的子像素级角反射器自动定位方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103247029A (zh) * 2013-03-26 2013-08-14 中国科学院上海技术物理研究所 一种用于拼接式探测器生成的高光谱图像几何配准方法
CN105403886A (zh) * 2015-12-02 2016-03-16 中国科学院电子学研究所 一种机载sar定标器图像位置自动提取方法
CN106249766A (zh) * 2016-09-09 2016-12-21 中国科学院电子学研究所 远程监控角反射器系统
CN106526593A (zh) * 2016-12-19 2017-03-22 国家测绘地理信息局卫星测绘应用中心 基于sar严密成像模型的子像素级角反射器自动定位方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
RUISHAN ZHANG ETAL: "Multimode Hybrid Geometric Calibration of Spaceborne SAR Considering Atmospheric Propagation Delay", 《REMOTE SENSING》 *
赵瑞山 等: "利用光学遥感影像进行星载SAR影像正射纠正", 《测绘通报》 *

Cited By (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109870685B (zh) * 2017-12-05 2023-01-03 南京理工大学 基于改进的rd算法的室内距离向移动sar成像方法
CN109870685A (zh) * 2017-12-05 2019-06-11 南京理工大学 基于改进的rd算法的室内距离向移动sar成像方法
CN108469614A (zh) * 2018-03-27 2018-08-31 武汉大学 基于未定标极化雷达影像的角反射器探测方法
CN108469614B (zh) * 2018-03-27 2021-09-14 武汉大学 基于未定标极化雷达影像的角反射器探测方法
CN108562882A (zh) * 2018-06-21 2018-09-21 武汉大学 一种星载sar影像几何交叉定标方法和系统
CN110865346A (zh) * 2018-08-28 2020-03-06 中国人民解放军61646部队 一种基于直接定位算法的星载sar时间参数标定方法
CN110865346B (zh) * 2018-08-28 2021-07-20 中国人民解放军61646部队 一种基于直接定位算法的星载sar时间参数标定方法
CN109581367A (zh) * 2018-12-19 2019-04-05 西安空间无线电技术研究所 一种天基预警雷达重频组优化设计方法
CN109581367B (zh) * 2018-12-19 2020-10-23 西安空间无线电技术研究所 一种天基预警雷达重频组优化设计方法
CN109696182A (zh) * 2019-01-23 2019-04-30 张过 一种星载推扫式光学传感器内方位元素定标方法
CN109782280A (zh) * 2019-03-06 2019-05-21 中国科学院电子学研究所 基于定标卫星的中高轨sar天线距离向方向图测量方法
CN110031841B (zh) * 2019-04-01 2021-07-23 中国科学院遥感与数字地球研究所 基于ECMWF的InSAR大气延迟改正的方法和系统
CN110031841A (zh) * 2019-04-01 2019-07-19 中国科学院遥感与数字地球研究所 基于ECMWF的InSAR大气延迟改正的方法和系统
CN110531329B (zh) * 2019-09-25 2023-06-27 中南大学 一种用于冲激脉冲式穿墙雷达的多通道误差自校正方法
CN110531329A (zh) * 2019-09-25 2019-12-03 中南大学 一种用于冲激脉冲式穿墙雷达的多通道误差自校正方法
CN110703247A (zh) * 2019-10-30 2020-01-17 中国科学院电子学研究所 一种波位信息检测方法、波位信息检测装置和存储介质
CN111596272A (zh) * 2020-06-01 2020-08-28 中国科学院空天信息创新研究院 基于角反射器的机载多极化sar系统外定标方法
CN111596272B (zh) * 2020-06-01 2022-11-29 中国科学院空天信息创新研究院 基于角反射器的机载多极化sar系统外定标方法
CN111915570B (zh) * 2020-07-10 2024-03-08 云南电网有限责任公司带电作业分公司 基于反向传播神经网络的大气延迟估计方法
CN111915570A (zh) * 2020-07-10 2020-11-10 云南电网有限责任公司带电作业分公司 基于反向传播神经网络的大气延迟估计方法
CN111896954A (zh) * 2020-08-06 2020-11-06 华能澜沧江水电股份有限公司 一种船载sar影像的角反射器坐标定位方法
CN112083387A (zh) * 2020-09-02 2020-12-15 南京慧尔视智能科技有限公司 一种雷达标定方法及装置
CN112083387B (zh) * 2020-09-02 2022-03-18 南京慧尔视智能科技有限公司 一种雷达标定方法及装置
CN112179850A (zh) * 2020-11-09 2021-01-05 南京信息工程大学 基于地面实测的gf-4 aod产品高时间分辨率检验
TWI735390B (zh) * 2020-12-01 2021-08-01 財團法人金屬工業研究發展中心 影像定位系統之即時定位補償方法及可即時定位補償之影像定位系統
CN112764052A (zh) * 2020-12-25 2021-05-07 中国人民解放军32181部队 防空导弹飞行监测系统
CN113093190A (zh) * 2021-04-08 2021-07-09 中国电子科技集团公司第三十八研究所 基于高精度组合惯导系统的机载条带sar图像定位方法
CN113446902A (zh) * 2021-07-07 2021-09-28 中国人民解放军32382部队 一种自动对消测角误差的防空导弹的外弹道监测系统
CN113742803A (zh) * 2021-09-07 2021-12-03 辽宁工程技术大学 一种中高轨道sar卫星带控几何定位精度仿真分析方法
CN113671505A (zh) * 2021-09-08 2021-11-19 辽宁工程技术大学 一种基于系统几何误差补偿的合成孔径雷达立体定位方法
CN113671505B (zh) * 2021-09-08 2024-01-23 辽宁工程技术大学 一种基于系统几何误差补偿的合成孔径雷达立体定位方法
CN113960545A (zh) * 2021-10-19 2022-01-21 安徽大学 基于对称几何构型约束的星载sar无场几何定标方法及其系统
CN114065531A (zh) * 2021-11-20 2022-02-18 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种实测数据驱动的电离层法拉第旋转效应估计方法
CN114594435A (zh) * 2022-01-18 2022-06-07 中国资源卫星应用中心 一种国产民用sar卫星几何检校与定位精度改进方法
CN115629363A (zh) * 2022-12-16 2023-01-20 西安空间无线电技术研究所 星载探测雷达实时内定标方法及装置
CN115629363B (zh) * 2022-12-16 2023-04-28 西安空间无线电技术研究所 星载探测雷达实时内定标方法及装置
CN115856807A (zh) * 2023-02-22 2023-03-28 中国科学院空天信息创新研究院 高分辨率sar卫星定标器雷达波相位中心精确定位方法
CN117761684A (zh) * 2023-11-15 2024-03-26 中国科学院空间应用工程与技术中心 一种天基干涉成像雷达高度计方位时延确定方法及系统
CN117665818A (zh) * 2024-02-02 2024-03-08 北京东方至远科技股份有限公司 一种针对合成孔径雷达卫星的平面位置修正方法及系统
CN117665818B (zh) * 2024-02-02 2024-04-12 北京东方至远科技股份有限公司 一种针对合成孔径雷达卫星的平面位置修正方法及系统
CN118091567A (zh) * 2024-04-28 2024-05-28 中山大学 基于正射影像能量的sar快时间慢时间参数虚拟定标方法

Also Published As

Publication number Publication date
CN107367716B (zh) 2019-09-10

Similar Documents

Publication Publication Date Title
CN107367716B (zh) 一种高精度星载sar几何定标方法
RU2531802C2 (ru) Способ определения географических координат точек изображения на sar изображениях
CN112305510B (zh) 一种基于dem匹配的合成孔径雷达影像几何定标方法
CN102654576B (zh) 基于sar图像和dem数据的图像配准方法
CN107238824B (zh) 基于先验dem数据的星载sar图像几何精校正方法
CN101551450B (zh) 一种星载极化sar法拉第旋转效应校正平台的构建方法
Brenot et al. Preliminary signs of the initiation of deep convection by GNSS
Schubert et al. Mitigation of atmospheric perturbations and solid Earth movements in a TerraSAR-X time-series
CN102866393B (zh) 一种基于pos与dem数据的sar多普勒参数估计方法
KR102151362B1 (ko) 극좌표변환을 이용한 항공기기반 영상복원장치 및 이를 이용한 영상복원방법
CN111721241A (zh) 一种GNSS-InBSAR和GB-InSAR跨系统融合三维形变测量方法
CN105929400B (zh) 基于精确测量的恒定prf下变航速机载sar数据补偿方法
CN111123345B (zh) 一种基于gnss测量的经验电离层模型数据驱动方法
CN113945955A (zh) 基于大气延迟误差修正提高海面测高精度的方法和系统
CN114594435B (zh) 一种国产民用sar卫星几何检校与定位精度改进方法
Nonaka et al. Evaluation of the geometric accuracy of TerraSAR-X
JP2009509125A (ja) 画像に関連する位置を決定するための方法および装置
RU2406071C1 (ru) Способ навигации движущихся объектов
Frey et al. Geometric error budget analysis for TerraSAR-X
CN111856464B (zh) 一种基于单控制点信息的车载sar的dem提取方法
Capaldo et al. Centimeter range measurement using amplitude data of TerraSAR-X imagery
Hu et al. Scan planning optimization for 2-D beam scanning using a future geostationary microwave radiometer
CN117347971B (zh) 一种激光雷达建筑物边界估测方法及计算机可读介质
CN113093190B (zh) 基于高精度组合惯导系统的机载条带sar图像定位方法
Li et al. High-Resolution Spaceborne SAR Geolocation Accuracy Analysis and Error Correction

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