CN108508439B - 双机载sar对目标协同成像立体定位的方法 - Google Patents
双机载sar对目标协同成像立体定位的方法 Download PDFInfo
- Publication number
- CN108508439B CN108508439B CN201810406947.0A CN201810406947A CN108508439B CN 108508439 B CN108508439 B CN 108508439B CN 201810406947 A CN201810406947 A CN 201810406947A CN 108508439 B CN108508439 B CN 108508439B
- Authority
- CN
- China
- Prior art keywords
- target
- imaging
- image
- coordinates
- points
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- 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
- G01S13/904—SAR modes
- G01S13/9058—Bistatic or multistatic SAR
-
- 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
- G01S13/9004—SAR image acquisition techniques
- G01S13/9005—SAR image acquisition techniques with optical processing of the SAR signals
Abstract
本发明公开了一种双机载SAR对目标协同成像立体定位的方法,包括:对场景区域返回的原始回波数据进行处理,经过距离向压缩、距离徙动校正及方位向压缩处理完成对目标场景的成像;分别对两幅高分辨率SAR图像进行双参数CFAR处理、根据其几何特征检测出目标区域,并分别取目标区域几何中心,输出目标灰度区域对应的像点坐标;将像点坐标矩阵代入定位解算模型,结合两SAR的系统参数、定向参数,经过牛顿迭代解算出目标在直角坐标系中的坐标。本发明可以在得到目标在图像中的像点坐标后直接进行解算,从而大大减少计算量,提高实时性。
Description
技术领域
本发明属于SAR定位技术领域,具体涉及一种双机载SAR对目标协同成像立体定位的方法。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)具有全天候、全天时的对地观测能力,SAR在海洋监测方面得到了广泛的应用,并发挥了巨大的社会、经济和军事效益。因成像波段不同,SAR系统可以获取不同于光学系统的地物目标信息。
在一般情况下,军事目标在雷达图像上都是由具有较亮的散射中心组成,因而可以采用阈值分割的技术实现目标的检测。然而在具有随机性和复杂性的杂波背景下,SAR所成的图像中,背景与目标的对比度随之变化,对固定阈值的图像分割技术往往会产生较多的虚假目标。基于滑动窗口的双参数恒虚警目标检测方法,是一种自适应门限检测方法,能够适应背景杂波的变化,是常用的目标检测算法之一。
SAR图像定位与三维信息提取技术可以对较大范围的区域实现全天时、全天候和高精度的目标定位。现有的双基机载SAR技术常用的是两部雷达分别作为发射机和接收机的工作模式,这样的方法无法由目标在SAR图像的像素位置信息,通过直接解析方程组获得目标的实际位置信息。
发明内容
本发明的目的在于提供一种双机载SAR对目标协同成像立体定位的方法。
实现本发明目的的技术方案为:一种双机载SAR对目标协同成像立体定位的方法,该方法包括以下步骤:
采用距离多普勒成像算法,对两个机载平台SAR1、SAR2回波数据进行成像处理;
采用目标检测算法,对SAR1、SAR2图像同时进行目标检测,输出目标分别在两幅图像中的像点坐标;
将像点坐标代入双SAR协同立体定位模型,对目标的实际三维坐标进行牛顿迭代解算。
本发明与现有技术相比,具有以下优点:本发明能够更加有效地对SAR图像中目标进行检测及定位;采用双参数CFAR检测算法可以自适应背景杂波的变化,得出最适合的门限;将背景窗口分为四个部分进行杂波均值分析可以大大减少计算量;采用两机载SAR数据融合定位的三维解算方法,可以不受相干性的限制,得到高精度定位结果。
附图说明
图1是本发明中双机载SAR对目标协同成像立体定位方法原理图。
图2是本发明中双机载SAR对目标协同成像立体定位方法总流程图。
图3是本发明中目标检测算法流程图。
图4是本发明中地面点三维坐标解算算法流程图。
图5是本发明中双参数CFAR检测窗口结构示意图。
图6是本发明中双参数CFAR背景窗口分块示意图。
具体实施方式
一种双机载SAR对目标协同成像立体定位的方法,该方法包括以下步骤:
步骤1,采用距离多普勒成像算法,对两个机载平台SAR1、SAR2回波数据进行成像处理;
步骤2,采用目标检测算法,在复杂背景下对SAR1、SAR2图像同时进行目标检测,输出目标分别在两幅图像中的像点坐标;具体为:
采用双参数CFAR算法分别对SAR1、SAR2所成图像检测指定目标,然后采用形态学滤波算法处理以消除检测结果中的虚假目标,再采用函数regionprops统计被标记的区域的面积分布,显示区域总数,输出检测结果,即目标在图中的像点坐标范围,最后分别取两组像点坐标矩阵的几何中心点作为目标具体像点坐标。
双参数CFAR算法的检测过程具体为:
(1a)设定一个检测单元由三个层层嵌套的窗口构成,三个窗口分别是目标窗口T、保护窗口P和背景窗口B,其中目标窗口内是待检测的点,上述3个窗口均正方形;
设定目标窗口边长为a,保护窗口的边长为b,背景窗口环形宽度为c,则得出用于杂波区域计算的像素数量记为numpix,表达式为:numpix=2c·(2c+2b);检测器边长记为d,表达式为:d=b+2c;
(1c)确定CFAR阈值
(1d)利用CFAR检测器求解局部阈值,执行单个像素判断,具体为:由杂波区域得到均值和标准偏差,计算双参数CFAR检测判别式,遍历所有像素点,判别是否为目标点;
将CFAR背景窗口划分为四个部分;
对每个像素(i,j)所对应的背景窗口区域估计,然后分别累加最后求背景区域的均值ub以及标准偏差σb;将该像素(i,j)对应的灰度值x代入双参数CFAR检测判别公式进行计算,满足则判别该像素点为目标,赋值为255,否则判别该像素点为背景,赋值为0;由此得到双参数CFAR检测后的二值矩阵Ia1、Ia2。
所述的形态学滤波具体步骤为:
(2a)创建一个以半径为r1的圆形结构元素矩阵B1,以结构元素B1对二值矩阵Ia1、Ia2进行闭运算,填补轮廓线断裂,得到结果Ib1、Ib2;
(2b)创建一个以半径为r2的圆形结构元素矩阵B2,以结构元素B2对二值矩阵Ib1、Ib2进行腐蚀运算,得到结果Ic1、Ic2;
(2c)以结构元素B2对二值矩阵Ic1、Ic2进行膨胀运算,得到结果Id1、Id2。
步骤3,将像点坐标代入双SAR协同立体定位模型,对目标的实际三维坐标进行牛顿迭代解算;具体为:
将两组同名点像点坐标代入距离多普勒成像模型,结合SAR1、SAR2的定向参数及控制点坐标,采用牛顿迭代算法最后迭代出最符合实际条件的目标三维坐标;其中同名点即两幅图像中指代同一目标的点,定向参数包括机载SAR的飞行速度和实时坐标。
所述的距离多普勒模型导入牛顿迭代解算三维坐标,具体步骤为:
(4a)由两幅SAR图像中同名点的像点T1(iL,jL)、T2(iR,jR)代入距离公式及多普勒公式,得到同名像点坐标与对应地面点坐标(X,Y,Z)的关系由如下四个方程构成的方程组表示:
即:
其中,分别代表像点T1、像点T2成像瞬间SAR1、SAR2天线相位中心位置,分别代表像点T1、像点T2成像瞬间SAR1、SAR2天线相位中心速度,分别为SAR1、SAR2成像时的近距延迟,分别为SAR1、SAR2的斜距向采样间隔;
(4b)计算同名点成像瞬间两天线相位中心位置、速度
式中,t'为每行间的时间间隔;tL和tR分别为像点T1、像点T2在左右图像上的成像时刻;分别为SAR1、SAR2的天线相位中心加速度初值;分别为SAR1、SAR2的天线相位中心速度初值;分别为SAR1、SAR2的天线相位中心位置初值;
(4c)给定地面点三维坐标初值
假设获取测区n个控制点,每个控制点对应的像点坐标(ik,jk)及实际坐标(Xk,Yk,Zk),分别为:pGC1(i1,j1,X1,Y1,Z1),pGC2(i2,j2,X2,Y2,Z2),...pGCn(in,jn,Xn,Yn,Zn),k=1,2,3...n;
则地面点三维坐标的初值为:
(4d)构建误差方程组
由R-D模型对地面点坐标的线性化形式,可知同名像点T1(iLjL)、T2(iR,jR)与所对应的地面点P(X,Y,Z)的线性化关系为:
C·△G-L=0
其中,C是关于地面点坐标改正量的系数矩阵即:
式中,分别是关于像点T1成像时刻对应的SAR1位置的函数;分别是关于像点T1成像时刻对应的SAR1速度的函数;分别是关于像点T2成像时刻对应的SAR2位置的函数; 分别是关于像点T2成像时刻对应的SAR2速度的函数。
系数阵C中的各元素分别为:
△G是地面点坐标的改正量向量,△G=[△X △Y △Z]T;
(4e)计算三维坐标改正量
C矩阵的法方程式表示为:
CTC△G-CTL=0
解法方程,可求得地面点三维坐标的改正量向量△G:
△G=(CTC)-1CTL
再在上一次迭代的基础上修正三维坐标初值:
(4f)限差判断
判断改正量是否小于给定限差,若改正量大于限差则返回(4d),利用改正后的三维坐标重新组误差方程式计算其改正量;若改正量小于等于限差则停止迭代,输出计算出的地面点三维坐标
下面结合附图和实施例对本发明进行详细说明。
实施例
本实施例以杂波海洋环境下定位舰船目标为例,对本发明的示例性实施方式进行详细描述。
图1为本发明提出的一种双机载SAR对目标协同成像立体定位的结构示意图,主要由两架固定翼无人机作为平台的两部SAR雷达组成,两架无人机以相同方向在目标区域两侧飞行,同时对扫描区域进行成像。
图2为本发明提出的一种双机载SAR对目标协同成像立体定位的总流程图,具体实施步骤为:
第一步,采用距离多普勒成像算法,对两个机载平台SAR1、SAR2接收到的海面扫描区域的回波数据进行高分辨率成像,分辨率达到1米。
所述的距离多普勒成像算法,包括:距离向压缩、距离徙动校正及方位向压缩处理。
距离向压缩的流程为,在频域中通过回波信号与匹配滤波函数相乘完成脉冲压缩,匹配滤波函数为接收频域信号所含的二次相位项的共轭函数。距离徙动校正的解决方法为在频域上乘以一个线性相位。再在方位向上通过与匹配滤波函数相乘完成脉冲压缩。
第二步,采用目标检测算法,流程图如图3所示,对SAR1、SAR2所成像进行杂波背景下舰船目标检测。所述的目标检测算法,具体步骤为:采用双参数CFAR算法分别对SAR1、SAR2所成图像检测舰船目标,然后采用形态学滤波算法处理以消除检测结果中的虚假目标,再采用函数regionprops统计被标记的区域的面积分布,输出检测结果,即舰船目标在图中的像点坐标范围,再分别取两组像点坐标矩阵的几何中心点作为目标具体像点坐标T1(iL,jL)、T2(iR,jR)。
所述的双参数CFAR检测舰船目标的流程为,确定CFAR检测器参数,包括窗口尺寸,保护区宽度,杂波区宽度,遍历扩展后的图像检测得到图像上的高亮像素。
选取合适的滑动窗口,对SAR图像中背景窗内参数的高斯分布进行估计,确定阈值将目标估计参数与之做比较,得出图像上的高亮像素区域。
所述的形态学滤波的流程为,从图像中提取连通区域,即可能的舰船目标,对各个区域进行几何特征提取,再基于特征排除不是舰船的目标。
所述的采用函数regionprops对细节进行进一步计算的流程为,采用函数regionprops统计被标记的区域的面积分布,显示区域总数。
所述的双参数CFAR检测舰船目标,具体步骤为:
1a)设定一个检测单元由三个层层嵌套的窗口构成,如图5所示。三个窗口分别是目标窗口T、保护窗口P和背景窗口B,其中目标窗口内是待检测的点,海面背景信息从背景窗口获取,保护窗口是为了确保舰船目标部分不会被包含在背景窗口中。上述3个窗口均正方形。
设定目标窗口边长为a,保护窗口的边长为b,背景窗口环形宽度为c,则得出用于杂波区域计算的像素数量记为numpix,表达式为:numpix=2c·(2c+2b)。检测器边长记为d,表达式为:d=b+2c。
1c)确定CFAR阈值
1d)利用CFAR检测器求解局部阈值,执行单个像素判断(由杂波区域得到均值和标准偏差,计算双参数CFAR检测判别式,遍历所有像素点,判别是否为目标点)
将CFAR背景窗口划分为四个部分,如图6所示。
对每个像素(i,j)所对应的背景窗口区域估计,然后分别累加最后求背景区域的均值ub以及标准偏差σb。将该像素(i,j)对应的灰度值x代入双参数CFAR检测判别公式进行计算,满足则判别该像素点为目标,赋值为255,否则判别该像素点为背景,赋值为0。由此得到双参数CFAR检测后的二值矩阵Ia1、Ia2。
所述的形态学滤波,具体步骤为:
2a)创建一个以半径为r1的圆形结构元素矩阵B1,以结构元素B1对二值矩阵Ia1、Ia2进行闭运算,填补轮廓线断裂,得到结果Ib1、Ib2。
2b)创建一个以半径为r2的圆形结构元素矩阵B2,以结构元素B2对二值矩阵Ib1、Ib2进行腐蚀运算,得到结果Ic1、Ic2。
2c)以结构元素B2对二值矩阵Ic1、Ic2进行膨胀运算,得到结果Id1、Id2。
所述的采用函数regionprops统计被标记的区域的面积分布,显示区域总数,具体步骤为:采用'Area'作为度量数据对图像各个区域中像素总个数进行计算,并确定每个连通区域的坐标范围。
第三步,采用地面点三维坐标解算算法,对舰船目标的实际坐标进行解算。流程图如图4所示,具体步骤为:将两组同名点像点坐标代入距离多普勒成像模型,结合SAR1、SAR2的定向参数及控制点坐标,采用牛顿迭代算法解算三维坐标的流程为计算同名点成像瞬间天线相位中心位置及速度,代入距离多普勒基本方程并构建误差方程,最后迭代出最符合实际条件的目标三维坐标。
距离多普勒基本方程包括距离公式及多普勒频率方程:
距离公式:Rs 2=(X-Xs)2+(Y-Ys)2+(Z-Zs)2=(R0+Mslant·j)2
可记:F1=(X-Xs)2+(Y-Ys)2+(Z-Zs)2-(R0+Mslant·j)2
其中,(X,Y,Z)表示地面点目标坐标,(Xs,Ys,Zs)为成像瞬间天线相位中心位置,R0为近距延迟,Mslant为斜距向采样间隔,j为像点的距离向坐标。
多普勒频率方程:
其中,(Vx,Vy,Vz)表示成像瞬间天线相位中心速度,λ为雷达发射波波长,Rs为地面点目标到雷达平台的瞬时位置距离,fdc为多普勒频移参数。
牛顿迭代解算三维坐标的流程为计算同名点成像瞬间天线相位中心位置及速度,代入距离多普勒基本方程并构建误差方程。
所述的距离多普勒模型导入牛顿迭代解算三维坐标,具体步骤为:
4a)由两幅SAR图像中同名点的像点T1(iL,jL)、T2(iR,jR)代入距离公式及多普勒公式,得到同名像点坐标与对应地面点坐标(X,Y,Z)的关系由如下四个方程构成的方程组表示:
即:
其中,分别代表像点T1、像点T2成像瞬间SAR1、SAR2天线相位中心位置,分别代表像点T1、像点T2成像瞬间SAR1、SAR2天线相位中心速度,分别为SAR1、SAR2成像时的近距延迟,分别为SAR1、SAR2的斜距向采样间隔。
4b)计算同名点成像瞬间两天线相位中心位置、速度
式中,t'为每行间的时间间隔;tL和tR分别为像点T1、像点T2在左右图像上的成像时刻;分别为SAR1、SAR2的天线相位中心加速度初值;分别为SAR1、SAR2的天线相位中心速度初值;分别为SAR1、SAR2的天线相位中心位置初值。
4c)给定地面点三维坐标初值
假设获取测区n个控制点,每个控制点对应的像点坐标(ik,jk)(k=1,2,3...n)及实际坐标(Xk,Yk,Zk),分别为:pGC1(i1,j1,X1,Y1,Z1),pGC2(i2,j2,X2,Y2,Z2),...pGCn(in,jn,Xn,Yn,Zn)
则地面点三维坐标的初值可取为:
4d)构建误差方程组
由R-D模型对地面点坐标的线性化形式,可知同名像点T1(iLjL)、T2(iR,jR)与所对应的地面点P(X,Y,Z)的线性化关系为:
C·△G-L=0
其中,C是关于地面点坐标改正量的系数矩阵即:
式中,分别是关于像点T1成像时刻对应的SAR1位置的函数;分别是关于像点T1成像时刻对应的SAR1速度的函数;分别是关于像点T2成像时刻对应的SAR2位置的函数; 分别是关于像点T2成像时刻对应的SAR2速度的函数。
系数阵C中的各元素分别为:
△G是地面点坐标的改正量向量,△G=[△X △Y △Z]T。
4e)计算三维坐标改正量
C矩阵的法方程式可表示为:
CTC△G-CTL=0
解法方程,可求得地面点三维坐标的改正量向量△G:
△G=(CTC)-1CTL
再在上一次迭代的基础上修正三维坐标初值:
4f)限差判断
判断改正量是否小于给定限差,若改正量大于限差则返回(4d),利用改正后的三维坐标重新组误差方程式计算其改正量;若改正量小于或等于限差则停止迭代,输出计算出的地面点三维坐标。
Claims (4)
1.一种双机载SAR对目标协同成像立体定位的方法,其特征在于,该方法包括以下步骤:
采用距离多普勒成像算法,对两个机载平台SAR1、SAR2回波数据进行成像处理;
采用目标检测算法,对SAR1、SAR2图像同时进行目标检测,输出目标分别在两幅图像中的像点坐标;
将像点坐标代入双SAR协同立体定位模型,对目标的实际三维坐标进行牛顿迭代解算;
所述的目标检测算法,具体步骤为:
采用双参数CFAR算法分别对SAR1、SAR2所成图像检测指定目标,然后采用形态学滤波算法处理以消除检测结果中的虚假目标,再采用函数regionprops统计被标记的区域的面积分布,显示区域总数,输出检测结果,即目标在图中的像点坐标范围,最后分别取两组像点坐标矩阵的几何中心点作为目标具体像点坐标;
双参数CFAR算法的检测过程具体为:
(1a)设定一个检测单元由三个层层嵌套的窗口构成,三个窗口分别是目标窗口T、保护窗口P和背景窗口B,其中目标窗口内是待检测的点,上述3个窗口均正方形;
设定目标窗口边长为a,保护窗口的边长为b,背景窗口环形宽度为c,则得出用于杂波区域计算的像素数量记为numpix,表达式为:numpix=2c·(2c+2b);检测器边长记为d,表达式为:d=b+2c;
(1c)确定CFAR阈值
(1d)利用CFAR检测器求解局部阈值,执行单个像素判断,具体为:由杂波区域得到均值和标准偏差,计算双参数CFAR检测判别式,遍历所有像素点,判别是否为目标点;
将CFAR背景窗口划分为四个部分;
2.根据权利要求1所述的双机载SAR对目标协同成像立体定位的方法,其特征在于,所述的形态学滤波具体步骤为:
(2a)创建一个以半径为r1的圆形结构元素矩阵B1,以结构元素B1对二值矩阵Ia1、Ia2进行闭运算,填补轮廓线断裂,得到结果Ib1、Ib2;
(2b)创建一个以半径为r2的圆形结构元素矩阵B2,以结构元素B2对二值矩阵Ib1、Ib2进行腐蚀运算,得到结果Ic1、Ic2;
(2c)以结构元素B2对二值矩阵Ic1、Ic2进行膨胀运算,得到结果Id1、Id2。
3.根据权利要求1所述的双机载SAR对目标协同成像立体定位的方法,其特征在于,将像点坐标代入双SAR协同立体定位模型,具体步骤为:
将两组同名点像点坐标代入距离多普勒成像模型,结合SAR1、SAR2的定向参数及控制点坐标,采用牛顿迭代算法最后迭代出最符合实际条件的目标三维坐标;所述同名点即两幅图像中指代同一目标的点,定向参数包括机载SAR的飞行速度和实时坐标。
4.根据权利要求3所述的双机载SAR对目标协同成像立体定位的方法,其特征在于,所述的距离多普勒模型导入牛顿迭代解算三维坐标,具体步骤为:
(4a)由两幅SAR图像中同名点的像点T1(iL,jL)、T2(iR,jR)代入距离公式及多普勒公式,得到同名像点坐标与对应地面点坐标(X,Y,Z)的关系由如下四个方程构成的方程组表示:
即:
其中,分别代表像点T1、像点T2成像瞬间SAR1、SAR2天线相位中心位置,分别代表像点T1、像点T2成像瞬间SAR1、SAR2天线相位中心速度,分别为SAR1、SAR2成像时的近距延迟,分别为SAR1、SAR2的斜距向采样间隔;
(4b)计算同名点成像瞬间两天线相位中心位置、速度
式中,t'为每行间的时间间隔;tL和tR分别为像点T1、像点T2在左右图像上的成像时刻;分别为SAR1、SAR2的天线相位中心加速度初值;分别为SAR1、SAR2的天线相位中心速度初值;分别为SAR1、SAR2的天线相位中心位置初值;
(4c)给定地面点三维坐标初值
假设获取测区n个控制点,每个控制点对应的像点坐标(ik,jk)及实际坐标(Xk,Yk,Zk),分别为:pGC1(i1,j1,X1,Y1,Z1),pGC2(i2,j2,X2,Y2,Z2),...pGCn(in,jn,Xn,Yn,Zn),k=1,2,3...n;
则地面点三维坐标的初值为:
(4d)构建误差方程组
由R-D模型对地面点坐标的线性化形式,可知同名像点T1(iLjL)、T2(iR,jR)与所对应的地面点P(X,Y,Z)的线性化关系为:
C·ΔG-L=0
其中,C是关于地面点坐标改正量的系数矩阵即:
式中,分别是关于像点T1成像时刻对应的SAR1位置的函数;分别是关于像点T1成像时刻对应的SAR1速度的函数;分别是关于像点T2成像时刻对应的SAR2位置的函数; 分别是关于像点T2成像时刻对应的SAR2速度的函数;
系数阵C中的各元素分别为:
ΔG是地面点坐标的改正量向量,ΔG=[ΔX ΔY ΔZ]T;
(4e)计算三维坐标改正量
C矩阵的法方程式表示为:
CTCΔG-CTL=0
解法方程,可求得地面点三维坐标的改正量向量ΔG:
ΔG=(CTC)-1CTL
再在上一次迭代的基础上修正三维坐标初值:
(4f)限差判断
判断改正量是否小于给定限差,若改正量大于限差则返回(4d),利用改正后的三维坐标重新组误差方程式计算其改正量;若改正量小于等于限差则停止迭代,输出计算出的地面点三维坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810406947.0A CN108508439B (zh) | 2018-05-01 | 2018-05-01 | 双机载sar对目标协同成像立体定位的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810406947.0A CN108508439B (zh) | 2018-05-01 | 2018-05-01 | 双机载sar对目标协同成像立体定位的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108508439A CN108508439A (zh) | 2018-09-07 |
CN108508439B true CN108508439B (zh) | 2022-02-18 |
Family
ID=63399880
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810406947.0A Active CN108508439B (zh) | 2018-05-01 | 2018-05-01 | 双机载sar对目标协同成像立体定位的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108508439B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109579843B (zh) * | 2018-11-29 | 2020-10-27 | 浙江工业大学 | 一种空地多视角下的多机器人协同定位及融合建图方法 |
CN110109113B (zh) * | 2019-05-07 | 2021-01-12 | 电子科技大学 | 一种基于级联对消的双基前视sar非平稳杂波抑制方法 |
CN111580105B (zh) * | 2020-06-02 | 2022-05-13 | 电子科技大学 | 一种用于太赫兹雷达高分辨成像的自适应处理方法 |
CN111896954A (zh) * | 2020-08-06 | 2020-11-06 | 华能澜沧江水电股份有限公司 | 一种船载sar影像的角反射器坐标定位方法 |
CN114339993B (zh) * | 2022-03-16 | 2022-06-28 | 北京瑞迪时空信息技术有限公司 | 基于天线距离约束的陆基定位方法、装置、设备和介质 |
CN114463365B (zh) * | 2022-04-12 | 2022-06-24 | 中国空气动力研究与发展中心计算空气动力研究所 | 红外弱小目标分割方法、设备及介质 |
CN115932823B (zh) * | 2023-01-09 | 2023-05-12 | 中国人民解放军国防科技大学 | 基于异源区域特征匹配的飞行器对地面目标定位方法 |
CN116985143B (zh) * | 2023-09-26 | 2024-01-09 | 山东省智能机器人应用技术研究院 | 一种打磨机器人的打磨轨迹生成系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105137431A (zh) * | 2015-08-06 | 2015-12-09 | 中国测绘科学研究院 | 一种sar立体模型构建与量测方法 |
-
2018
- 2018-05-01 CN CN201810406947.0A patent/CN108508439B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105137431A (zh) * | 2015-08-06 | 2015-12-09 | 中国测绘科学研究院 | 一种sar立体模型构建与量测方法 |
Non-Patent Citations (3)
Title |
---|
"SAR图像高精度定位技术研究";张红敏;《中国博士学位论文全文数据库 信息科技辑》;20140115;第34-51页 * |
"机载SAR图像的海上目标检测方法的研究";李昭瑞;《中国优秀硕士学位论文全文数据库 信息科技辑》;20110915;第41-42,51-52页 * |
张红敏."SAR图像高精度定位技术研究".《中国博士学位论文全文数据库 信息科技辑》.2014,第34-51页. * |
Also Published As
Publication number | Publication date |
---|---|
CN108508439A (zh) | 2018-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108508439B (zh) | 双机载sar对目标协同成像立体定位的方法 | |
CN109188433B (zh) | 基于无控制点的双机载sar图像目标定位的方法 | |
CN101738614B (zh) | 基于时空图像序列的isar目标转动估计方法 | |
CA2579898C (en) | Method for the processing and representing of ground images obtained by synthetic aperture radar systems (sar) | |
CN111913158B (zh) | 复杂杂波背景下检测低慢小目标的雷达信号处理方法 | |
CN108761484A (zh) | 一种基于多源卫星遥感数据的海雾监测方法 | |
CN104076360B (zh) | 基于压缩感知的二维sar稀疏目标成像方法 | |
CN107783104A (zh) | 一种基于粒子滤波的多异步传感器单目标检测前跟踪方法 | |
CN113484860B (zh) | 基于多普勒中心异常的sar图像海洋锋检测方法及系统 | |
CN108230375A (zh) | 基于结构相似性快速鲁棒的可见光图像与sar图像配准方法 | |
CN112184749A (zh) | 基于视频sar跨域联合的动目标跟踪方法 | |
JPH0980146A (ja) | レーダ装置 | |
CN116027318A (zh) | 多传感器信号融合的方法、装置、电子设备及存储介质 | |
CN110554377B (zh) | 基于多普勒中心偏移的单通道sar二维流场反演方法及系统 | |
CN101620272A (zh) | 一种逆合成孔径雷达的目标转速估计方法 | |
CN113589287B (zh) | 合成孔径雷达稀疏成像方法、装置、电子设备及存储介质 | |
CN109190647B (zh) | 一种有源无源数据融合方法 | |
CN113177593B (zh) | 一种水上交通环境中雷达点云与影像数据的融合方法 | |
CN107729903A (zh) | 基于区域概率统计和显著性分析的sar图像目标检测方法 | |
CN110780299A (zh) | 散度场的获取方法、装置、计算机设备和存储介质 | |
Lu et al. | Research on rainfall identification based on the echo differential value from X-band navigation radar image | |
CN111289953B (zh) | 一种基于模糊矩阵更新的天基雷达距离/速度解模糊方法 | |
CN110618403B (zh) | 一种基于双波束雷达的着陆飞行器参数测量方法 | |
CN110231603B (zh) | 一种基于gmti的快速解算目标速度的方法 | |
CN113419240A (zh) | 基于双通道sar的动目标检测方法,双通道sar及存储介质 |
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 |