CN1508562A - 机载合成孔径雷达测量区域定位系统 - Google Patents
机载合成孔径雷达测量区域定位系统 Download PDFInfo
- Publication number
- CN1508562A CN1508562A CNA02155885XA CN02155885A CN1508562A CN 1508562 A CN1508562 A CN 1508562A CN A02155885X A CNA02155885X A CN A02155885XA CN 02155885 A CN02155885 A CN 02155885A CN 1508562 A CN1508562 A CN 1508562A
- Authority
- CN
- China
- Prior art keywords
- measured zone
- aircraft
- positioning system
- point
- reference target
- 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
Images
Landscapes
- Navigation (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
机载合成孔径雷达(SAR)测量区域定位系统,涉及雷达技术,由全球定位系统,地理信息定位平台和测量区域定位系统组成,其全球定位系统与地理信息定位平台无线连接,地理信息定位平台与测量区域定位系统电连接;测量区域定位系统包括测量区域定位模块、测量带内参考目标参数计算模块、地球模型模块和坐标转换模块。本发明能够对机载合成孔径雷达的地面测量区域进行精确定位,并完成区域内参考目标的位置参数测量,可用于最优飞行航迹的预先确立,合成孔径雷达图像到地面坐标系的投影变换,以及地面标准参考目标在雷达图像辐射和几何定标中的应用等。
Description
技术领域
本发明涉及雷达技术,特别是涉及孔径雷达测量区域的确定和区域内参考目标的定位技术。
背景技术
机载合成孔径雷达在成像雷达中占有绝对重要的地位,有着广泛的应用前景和发展潜力,作为有源系统,机载合成孔径雷达的全天候、全天时工作能力以及它在不同频段、不同极化下可以得到目标的高分辨率图像,为人们提供各种非常有用的信息。这些信息广泛用于地质、水文、农业、城市、海洋、测绘以及军事侦察等领域,特别对于大面积测量,机载合成孔径雷达具有其它种类传感器无法比拟的优势。机载合成孔径雷达对地观测存在一个测量区域,该测量区域由机载合成孔径雷达的性能和飞行航迹来确定。在进行对地观测前,如何确定测量区域变得很重要,根据需要的观测区域,我们可以设计最佳的飞行路径,调整雷达系统的工作状态;如果在观测区域内设置参考目标用于定量雷达遥感,则需要确定参考目标的准确位置,因此,研究机载合成孔径雷达(SAR)测量区域定位系统就具有重要意义,它提供了理论上的测量参考,可以指导整个测量工作的进行。
目前,确定地面参考标准目标的指北角度时,是通过在地面经纬度图上,于目标参考点同一经度处为度量变化获得指北参考矢量,这样作存在一定的缺点,它没有考虑地球曲率和磁偏角的存在,因此它的精度不是最准确的。
发明内容
本发明的目的在于建立一种机载合成孔径雷达(SAR)测量区域定位系统,用于测量区域的预先确定和区域内参考目标的精确定位。
为达到上述目的,本发明的技术解决方案为提供一种机载合成孔径雷达测量区域定位系统,由全球定位系统,地理信息定位平台和测量区域定位系统组成,其全球定位系统与地理信息定位平台无线连接,地理信息定位平台与测量区域定位系统电连接;测量区域定位系统包括测量区域定位模块、测量带内参考目标参数计算模块、地球模型模块和坐标转换模块;测量区域定位模块的输入端直接与地理信息定位平台电连接,其输出端与测量带内参考目标参数计算模块电连接,地球模型模块和坐标转换模块的输出端也与测量区域定位模块电连接。
所述的测量区域定位系统,其对测量区域的定位由以下述步骤完成:
开始后,机载合成孔径雷达由地理信息定位平台的接收接口读入全球定位系统发送的飞机起始经纬度、飞机高度、近地点斜距、远地点斜距等数据,并输入测量区域定位模块;同时,测量区域定位模块将飞机起始经纬度通过与坐标转换模块交换数据,而转换成惯性坐标系坐标,并计算出飞机起始点的机下点,再以飞机的机下点坐标为初始值,求飞机起始点地面测绘带的近距点、远距点;测量区域定位模块将飞机起始点地面测绘带的近距点、远距点惯性坐标系坐标数据通过与坐标转换模块交换,而转换成经纬度坐标,再根据天线姿态,确定飞机与天线倾斜方向相同的一侧的测绘带各端点,并将测绘带各端点连续排列在地图上,在显示屏中显示出来,结束。
所述的测量区域定位系统,其在对测量带内参考目标测量时,参考目标的定位由以下述步骤完成:
开始后,由地理信息定位平台向测量区域定位模块输入角反射器的点位置,测量区域定位模块将角反射器点位置的数据通过与坐标转换模块交换,而转换成惯性坐标系坐标,将此数据和地球模型模块在惯性坐标系中的数据传入测量带内参考目标参数计算模块;测量带内参考目标参数计算模块先计算角反射器对准的飞机点,再计算此飞机点的机下点,然后计算高度角,进而计算北极点在过角反射器点切面上的投影点,计算出指北角,最后,将测量带内参考目标的位置输入到地图上,在显示屏中显示出来,结束。
所述的测量区域定位系统,其在求地面投影点时,采用了自适应迭代算法。
所述的测量区域定位系统,其在地面参考目标位置指向的确定中,校正了磁偏角和地球曲率的影响,由此确定的指北角,准确度较高。
本发明涉及到地球模型,一般的测量定位系统认为飞机和它在地面投影点的经纬度是一样的,实际使用地球模型时,高度误差会导致飞机和它的地面投影点经纬度不一致。为了更精确测量,需要考虑飞机高度带来的影响,本发明考虑了飞机不同高度时经纬度数据的变化。
计算目标的指北角度,可以通过在地面经纬度图,于目标参考点同一经度处为度量变化获得指北参考矢量;这样做存在一定的缺点,它没有考虑地球曲率和磁偏角的存在,因此它的精度不是最准确的。本发明考虑到地球曲率和磁偏角的影响,获得更为准确的测量结果。
附图说明
图1本发明机载合成孔径雷达测量区域定位系统;
图2本发明机载合成孔径雷达测量区域定位系统在实际航测时拍摄的地面测量区域定位图;
图3惯性地心坐标系;
图4发明测量区域定位流程;
图5发明预先测定的测量区域;
图6本发明测量带内参考目标的定位流程;
图7本发明测量区域内参考目标分布。
具体实施方式
本发明机载合成孔径雷达测量区域定位系统的主要特点是利用全球定位系统(GPS)和地理信息系统(GIS),结合飞机飞行航迹,建立空地几何关系,再利用软件算法计算地面投影点坐标。该系统的组成参考图1,其由全球定位系统(GPS)1,地理信息定位平台(GIS)2和测量区域定位系统3组成(图1虚框内的部分),其中,全球定位系统(GPS)1与地理信息定位平台(GIS)2无线连接,地理信息定位平台(GIS)2与测量区域定位系统3电连接;测量区域定位系统3包括测量区域定位模块4、测量带内参考目标参数计算模块5、地球模型模块6和坐标转换模块7;测量区域定位模块4的输入端直接与地理信息定位平台(GIS)2电连接,其输出端与测量带内参考目标参数计算模块5电连接,地球模型模块6和坐标转换模块7的输出端也与测量区域定位模块4电连接。
本发明涉及到地球模型,一般的测量定位系统认为飞机和它在地面投影点的经纬度是一样的。实际使用地球模型时,高度误差会导致飞机和它的地面投影点经纬度不一致。为了更精确测量,需要考虑飞机高度带来的影响,本发明考虑了飞机不同高度时经纬度数据的变化。
本发明在求解地面投影点时,采用了自适应迭代算法,对三维二次方程,采用非线性方程组的牛顿-拉斐森方法,其中初始值的选择是解法收敛和稳定的关键。在本发明中,该初始值的选择使用了雷达系统本身的先验知识,使该初始值的选择合理,求解投影点的结果收敛性和稳定性更好。
对于测量区域内的地面参考目标在确定其指北位置时,不是在地图经纬度平面上来计算其指北角,而是考虑了地球曲率以及磁偏角的影响,由此确定的指北角准确度更高。
首先建立惯性地心坐标系(见图3),所有计算统一在惯性地心坐标系中完成。
本发明机载合成孔径雷达测量区域定位系统在测量时,对测量区域的定位由以下述步骤完成:
如图4所示,开始后,机载合成孔径雷达由地理信息定位平台(GIS)2的接收接口读入全球定位系统(GPS)1发送的飞机起始经纬度、飞机高度、近地点斜距、远地点斜距等数据,并输入测量区域定位模块4;同时,测量区域定位模块4将飞机起始经纬度通过与坐标转换模块7交换数据,而转换成惯性坐标系坐标,并计算出飞机起始点的机下点,再以飞机的机下点坐标为初始值,用自适应迭代算法求飞机起始点地面测绘带的近距点、远距点;测量区域定位模块4将飞机起始点地面测绘带的近距点、远距点惯性坐标系坐标数据通过与坐标转换模块7交换,而转换成经纬度坐标,再根据天线姿态,确定飞机与天线倾斜方向相同的一侧的测绘带各端点,并将测绘带各端点连续排列在地图上,在显示屏中显示出来,见图5,结束。
所述测量区域定位模块4的运算是统一在惯性地心坐标系统中,并利用地理信息定位平台(GIS)2来完成的,其具体算法如下:
考虑地球模型和飞机的高度选择飞行起点坐标(x1,y1,z1)和终点坐标(x2,y2,z2),然后确定它们各自对应的地面波束点坐标(xg11,yg11,zg11)和(xg21,yg21,zg21)。各波束点坐标满足以下三个条件:
(1)正侧视工作模式时,矢量垂直,即
(xg11-x1)(x2-x1)+(yg11-y1)(y2-y1)+(zg11-z1)(z2-z1)=0 (1)
(2)波束点位于地球面上
其中,α=6378.145km,b=6356.76km。
(3)雷达系统测量的近距点或远距点距离由雷达系统设定,满足
(xg11-x1)2+(yg11-y1)2+(zg11-z1)2=Rnear (3)
(xg11-x1)2+(yg11-y1)2+(zg11-z1)2=Rfar (4)
其中,Rnear为近距点距离,Rfar为远距点距离
根据上述规则,可以得到一个三维二次方程组。解该非线性方程组采用牛顿-拉斐森方法。对于非线性方程组
fi(x1,x2...,xn)=0, i=1,2,...,N (5)
牛顿-拉斐森方法如下:
令x=(x1,x1,...,xn)T,i=1,2,...,N,在x的邻域做泰勒展开,略去二次和二次以上的项,得
α≡[αij]N×N,β≡(β1,...,βn)T
如果detα≠0,那么迭代公式为
其中(δx1,δx2,...,δxN)T≡δx为线性方程组α·δx=β的解。综上所述,使用牛顿一拉斐森方法解非线性方程组的步骤为:
(1)给定根x的初始近似x(0)(靠近x),允许误差为ε1,ε2,假设已经得到第k次近似x(k)。
(2)计算
可得
及
(3)计算
S1=|f1(x(k))|+|f2(x(k))|+...|fN(x(k))| (10)
如果S1<ε1,则计算结束,x(k)作为满足精度要求的近似解;否则,执行(4)。
(4)用LU分解法求线性代数方程组
α(k)·δx(k)=β(k) (11)
得
(5)计算x(k+1)=x(k)+δx(k)及
如果S2<ε2,则计算结束,x(k+1)作为满足精度要求的近似解;否则, 转向(2)继续计算,直到满足精度要求或者迭代次数已经达到给定的迭代次数为止。
经过上面计算处理,完成坐标转换和投影点的确定,相应地完成了测量区域的确定(如图5所示)。
在测量区域确定后,对于其中设置的地面标准参考目标需要确定其方位方向,本发明机载合成孔径雷达测量区域定位系统,在对测量带内参考目标测量时,参考目标的定位由以下述步骤完成:
如图6所示,开始后,由地理信息定位平台(GIS)2向测量区域定位模块4输入角反射器的点位置,测量区域定位模块4将角反射器点位置的数据通过与坐标转换模块7交换,而转换成惯性坐标系坐标,将此数据和地球模型模块6在惯性坐标系中的数据传入测量带内参考目标参数计算模块5;测量带内参考目标参数计算模块5先计算角反射器对准的飞机点,再计算此飞机点的机下点,然后计算高度角,进而计算北极点在过角反射器点切面上的投影点,在此基础上,校正了磁偏角和地球曲率后,计算出指北角,最后,将测量带内参考目标的位置输入到地图上,在显示屏中显示出来,见图7,结束。
实施例1:
本发明已多次成功用于机载合成孔径雷达的飞行测量工作,预先确定的测量区域(如图5所示)与实际飞行结果吻合。工作中曾发现测量区域不满足要求,经过该系统的分析发现飞机的飞行航迹没有满足要求,找到了问题的原因,及时进行了修正,挽回了不必要的损失。实验中使用的系统用于计算地面参考目标的位置参数,再用于实际实验中,在雷达图像上获得了地面参考目标准确的图像数值结果,如图7所示。
实施例2:
图2所示为本发明机载合成孔径雷达测量区域定位系统,在实际航测时拍摄的地面测量区域定位图。
Claims (5)
1.一种机载合成孔径雷达测量区域定位系统,由全球定位系统,地理信息定位平台和测量区域定位系统组成,其特征在于,全球定位系统与地理信息定位平台无线连接,地理信息定位平台与测量区域定位系统电连接;测量区域定位系统包括测量区域定位模块、测量带内参考目标参数计算模块、地球模型模块和坐标转换模块;测量区域定位模块的输入端直接与地理信息定位平台电连接,其输出端与测量带内参考目标参数计算模块电连接,地球模型模块和坐标转换模块的输出端也与测量区域定位模块电连接。
2.如权利要求1所述的测量区域定位系统,其特征在于,对测量区域的定位以下述步骤完成:
开始后,机载合成孔径雷达由地理信息定位平台的接收接口读入全球定位系统发送的飞机起始经纬度、飞机高度、近地点斜距、远地点斜距等数据,并输入测量区域定位模块;同时,测量区域定位模块将飞机起始经纬度通过与坐标转换模块交换数据,而转换成惯性坐标系坐标,并计算出飞机起始点的机下点,再以飞机的机下点坐标为初始值,求飞机起始点地面测绘带的近距点、远距点;测量区域定位模块将飞机起始点地面测绘带的近距点、远距点惯性坐标系坐标数据通过与坐标转换模块交换,而转换成经纬度坐标,再根据天线姿态,确定飞机与天线倾斜方向相同的一侧的测绘带各端点,并将测绘带各端点连续排列在地图上,在显示屏中显示出来,结束。
3.如权利要求1所述的测量区域定位系统,其特征在于,在对测量带内参考目标测量时,参考目标的定位以下述步骤完成:
开始后,由地理信息定位平台向测量区域定位模块输入角反射器的点位置,测量区域定位模块将角反射器点位置的数据通过与坐标转换模块交换,而转换成惯性坐标系坐标,将此数据和地球模型模块在惯性坐标系中的数据传入测量带内参考目标参数计算模块;测量带内参考目标参数计算模块先计算角反射器对准的飞机点,再计算此飞机点的机下点,然后计算高度角,进而计算北极点在过角反射器点切面上的投影点,计算出指北角,最后,将测量带内参考目标的位置输入到地图上,在显示屏中显示出来,结束。
4.如权利要求1所述的测量区域定位系统,其特征在于,在求地面投影点时,采用了自适应迭代算法。
5.如权利要求1所述的测量区域定位系统,其特征在于,在地面参考目标位置指向的确定中,校正了磁偏角和地球曲率的影响,由此确定的指北角,准确度较高。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB02155885XA CN1303431C (zh) | 2002-12-13 | 2002-12-13 | 机载合成孔径雷达测量区域定位系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB02155885XA CN1303431C (zh) | 2002-12-13 | 2002-12-13 | 机载合成孔径雷达测量区域定位系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1508562A true CN1508562A (zh) | 2004-06-30 |
CN1303431C CN1303431C (zh) | 2007-03-07 |
Family
ID=34236073
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB02155885XA Expired - Fee Related CN1303431C (zh) | 2002-12-13 | 2002-12-13 | 机载合成孔径雷达测量区域定位系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1303431C (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1329743C (zh) * | 2005-05-10 | 2007-08-01 | 中国科学院电子学研究所 | 一种机载干涉合成孔径雷达基线的组合测量装置 |
CN100384013C (zh) * | 2005-05-20 | 2008-04-23 | 中国科学院电子学研究所 | 一种机载干涉合成孔径雷达的收发天线分置结构 |
CN100526912C (zh) * | 2006-06-02 | 2009-08-12 | 中国科学院电子学研究所 | 一种宽带合成孔径雷达的有源外定标器及其定标方法 |
CN101082666B (zh) * | 2006-06-02 | 2011-07-27 | 中国科学院电子学研究所 | 基于自动测试系统实现对雷达脉冲信号高精度调制的方法 |
CN102135612A (zh) * | 2010-12-29 | 2011-07-27 | 电子科技大学 | 一种双基地前视合成孔径雷达测绘带范围计算方法 |
CN102736125A (zh) * | 2010-12-02 | 2012-10-17 | 霍尼韦尔国际公司 | 用于改善相关气象确定的系统和方法 |
CN101571593B (zh) * | 2008-04-30 | 2012-11-21 | 北京航空航天大学 | 一种星载sar图像的严密共线方程模型 |
CN103697885A (zh) * | 2013-12-17 | 2014-04-02 | 中国电子科技集团公司第十一研究所 | 自动补偿磁偏角的远程定位方法 |
CN106483499A (zh) * | 2016-09-21 | 2017-03-08 | 深圳智航无人机有限公司 | 无人机定位系统及无人机起飞、降落方法 |
CN113848549A (zh) * | 2021-09-15 | 2021-12-28 | 电子科技大学 | 一种基于合成孔径技术的辐射源定位方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5432520A (en) * | 1993-10-18 | 1995-07-11 | Hughes Aircraft Company | SAR/GPS inertial method of range measurement |
DE4419359A1 (de) * | 1994-06-03 | 1995-12-07 | Wolfram Dipl Ing Kirchner | Verfahren zur Erfassung, Auswertung, Ausmessung und Speicherung von Geo-Informationen |
-
2002
- 2002-12-13 CN CNB02155885XA patent/CN1303431C/zh not_active Expired - Fee Related
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1329743C (zh) * | 2005-05-10 | 2007-08-01 | 中国科学院电子学研究所 | 一种机载干涉合成孔径雷达基线的组合测量装置 |
CN100384013C (zh) * | 2005-05-20 | 2008-04-23 | 中国科学院电子学研究所 | 一种机载干涉合成孔径雷达的收发天线分置结构 |
CN100526912C (zh) * | 2006-06-02 | 2009-08-12 | 中国科学院电子学研究所 | 一种宽带合成孔径雷达的有源外定标器及其定标方法 |
CN101082666B (zh) * | 2006-06-02 | 2011-07-27 | 中国科学院电子学研究所 | 基于自动测试系统实现对雷达脉冲信号高精度调制的方法 |
CN101571593B (zh) * | 2008-04-30 | 2012-11-21 | 北京航空航天大学 | 一种星载sar图像的严密共线方程模型 |
CN102736125A (zh) * | 2010-12-02 | 2012-10-17 | 霍尼韦尔国际公司 | 用于改善相关气象确定的系统和方法 |
CN102135612B (zh) * | 2010-12-29 | 2012-08-29 | 电子科技大学 | 一种双基地前视合成孔径雷达测绘带范围计算方法 |
CN102135612A (zh) * | 2010-12-29 | 2011-07-27 | 电子科技大学 | 一种双基地前视合成孔径雷达测绘带范围计算方法 |
CN103697885A (zh) * | 2013-12-17 | 2014-04-02 | 中国电子科技集团公司第十一研究所 | 自动补偿磁偏角的远程定位方法 |
CN106483499A (zh) * | 2016-09-21 | 2017-03-08 | 深圳智航无人机有限公司 | 无人机定位系统及无人机起飞、降落方法 |
CN106483499B (zh) * | 2016-09-21 | 2023-09-15 | 深圳智航无人机有限公司 | 无人机定位系统及无人机起飞、降落方法 |
CN113848549A (zh) * | 2021-09-15 | 2021-12-28 | 电子科技大学 | 一种基于合成孔径技术的辐射源定位方法 |
CN113848549B (zh) * | 2021-09-15 | 2023-06-23 | 电子科技大学 | 一种基于合成孔径技术的辐射源定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1303431C (zh) | 2007-03-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108921947B (zh) | 生成电子地图的方法、装置、设备、存储介质以及采集实体 | |
JP4685313B2 (ja) | 任意の局面の受動的な体積画像の処理方法 | |
CN102967305B (zh) | 基于大小回字标志物的多旋翼无人机位姿获取方法 | |
AU2021439678A1 (en) | Method for extracting three-dimensional surface deformation by combining unmanned aerial vehicle doms and satellite-borne sar images | |
US20120133550A1 (en) | Method for Determining the Geographic Coordinates of Pixels in SAR Images | |
EP2442275A2 (en) | Method and apparatus for three-dimensional image reconstruction | |
US8817093B2 (en) | Photogrammetric networks for positional accuracy | |
CN112444834A (zh) | 一种定位方法及电子设备 | |
CN104111445A (zh) | 用于室内导航的超声波阵列辅助定位方法与系统 | |
JP6910511B2 (ja) | レーザ計測方法、レーザ計測用標識、及び座標算出プログラム | |
CN105203023A (zh) | 一种车载三维激光扫描系统安置参数的一站式标定方法 | |
CN1303431C (zh) | 机载合成孔径雷达测量区域定位系统 | |
CN101545776A (zh) | 基于数字地图的数码像片方位元素获取方法 | |
CN113341384A (zh) | 雷达校准 | |
CN102346033A (zh) | 基于卫星观测角误差估计的直接定位方法及系统 | |
CN109146958A (zh) | 一种基于二维图像的交通标志空间位置测量方法 | |
Le VAN et al. | Experimental investigation on the performance of DJI phantom 4 RTK in the PPK mode for 3D mapping open-pit mines | |
Liu et al. | Accurate mapping method for UAV photogrammetry without ground control points in the map projection frame | |
CN114137560A (zh) | 基于改进激光描述符的车辆重定位方法、装置及电子设备 | |
Yu et al. | Automatic extrinsic self-calibration of mobile LiDAR systems based on planar and spherical features | |
CN113959437A (zh) | 一种用于移动测量设备的测量方法及系统 | |
CN113960614A (zh) | 一种基于帧-地图匹配的高程图构建方法 | |
Lee et al. | Georegistration of airborne hyperspectral image data | |
CN111397602A (zh) | 一种宽频电磁指纹与组合导航融合的高精度定位方法与设备 | |
CN112985398A (zh) | 目标定位方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |