CN107765245B - 基于机载sar立体像对计算地面点三维位置的方法 - Google Patents
基于机载sar立体像对计算地面点三维位置的方法 Download PDFInfo
- Publication number
- CN107765245B CN107765245B CN201711002479.2A CN201711002479A CN107765245B CN 107765245 B CN107765245 B CN 107765245B CN 201711002479 A CN201711002479 A CN 201711002479A CN 107765245 B CN107765245 B CN 107765245B
- Authority
- CN
- China
- Prior art keywords
- image
- dimensional position
- airplane
- dimensional
- 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.)
- 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/904—SAR modes
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本公开提供了一种基于机载SAR立体像对计算地面点三维位置的方法,包括以下步骤:提取地面点在SAR立体像对上的图像坐标及飞机速度,根据所述图像坐标确定获取立体像对时对应的飞机三维位置;根据所述飞机三维位置计算基线参数;根据所述图像坐标、飞机三维位置及飞机速度计算变换参数;以及根据所述基线参数及变换参数计算地面点的三维位置。本公开基于机载SAR立体像对计算地面点三维位置的方法,通过机载SAR立体像对的平台参数和系统参数直接显式地计算地面点的三维位置,无需地面点的位置初始值,计算中也无任何的迭代和近似过程,提高了计算效率,保证了三维位置计算结果的严密性。
Description
技术领域
本公开涉及航空遥感图像处理技术领域,尤其涉及一种基于机载SAR立体像对计算地面点三维位置的方法。
背景技术
机载合成孔径雷达(英文简称SAR,下文都用SAR表示)遥感图像经过几何定位或几何校正后只能获取地面二维图像信息,为了得到地面的三维位置信息,可以采用立体像对的方法来得到地面点的三维信息。立体像对是由覆盖同一区域的左右各1幅共两幅图像构成,并利用空间几何关系来重建地面点的高度信息。
目前国内外对利用机载SAR立体像对直接计算地面点高度的研究极少,一般都是通过间接方法采用迭代求解的技术来实现地面点的高度提取。如邱春平(邱春平、张红敏、靳国旺、秦志远,机载SAR影像距离-斜视角模型立体定位,《测绘科学技术学报》2014年,第31卷第4期)建立了SAR影像的距离-斜视角模型立体定位方法,是根据SAR立体像对的左右影像建立4个方程通过给定地面点初始坐标采用最小二乘迭代解算地面点的三维坐标。高力(高力、赵杰、王仁礼,利用Leberl模型进行机载SAR图像的立体定位,《测绘学院学报》,第21卷第4期2004年12月)为实现SAR影像立体定位的方法,提出基于Leberl的SAR构像模型和解算过程,通过在立体SAR图像重叠区域的同名像点分别可以建立两个方程,再进行迭代计算来得到地面点的三维坐标;张红敏(张红敏、靳国旺、徐青、李向英,利用单个地面控制点的SAR图像高精度立体定位,《雷达学报》,第3卷第1期,2014年2月)则根据SAR建立了斜侧视距离-多普勒模型,得到了同名像点与相应地面点关系的非线性方程组,需要采用线性化迭代解算得到地面点的三维位置。方勇则(方勇、常本义、胡海彦、陈虹,机载SAR图像数字测图技术研究,《测绘通报》,2006年第8期)采用了光学立体摄影测量的方法处理SAR立体像对,通过等效共线方程重建SAR图像外方位元素和立体模型,也需要进行初值计算、误差方程建立和迭代解算等过程。
上述这些方法都是采用相应的模型来计算得到地面点的高度,但是解算过程中都是采用了复杂的非线性方程组,需要给定地面点三维位置的初始近似值,并通过多次迭代计算才能给出高度的最佳值,工程应用中十分不便。
发明内容
(一)要解决的技术问题
鉴于上述技术问题,本公开提供了一种基于机载SAR立体像对计算地面点三维位置的方法,通过机载SAR立体像对的平台参数和系统参数直接显式地计算地面点的三维位置,无需地面点的位置初始值,计算中也无任何的迭代和近似过程,提高了计算效率,保证了三维位置计算结果的严密性。
(二)技术方案
根据本公开的一个方面,提供了一种基于机载SAR立体像对计算地面点三维位置的方法,包括以下步骤:提取地面点在SAR立体像对上的图像坐标及飞机速度,根据所述图像坐标确定获取立体像对时对应的飞机三维位置;根据所述飞机三维位置计算基线参数;根据所述图像坐标、飞机三维位置及飞机速度计算变换参数;以及根据所述基线参数及变换参数计算地面点的三维位置。
在本公开的一些实施例中,所述提取地面点在SAR立体像对上的图像坐标及飞机速度,根据所述图像坐标确定获取立体像对时对应的飞机三维位置的步骤,包括:提取地面点在SAR立体像对的左图像上的二维图像坐标(x1,y1)及右图像上的二维图像坐标(x2,y2);提取SAR立体像对的左图像对应的飞机速度及右图像对应的飞机速度以及根据飞机平台上GPS记录的位置数据和图像坐标拟合得到左图像所对应的飞机三维位置和右图像所对应的飞机三维位置
在本公开的一些实施例中,所述基线参数为:
在本公开的一些实施例中,利用图像匹配软件或人工方法,提取地面点在SAR立体像对的左图像上的二维图像坐标(x1,y1)和右图像上的二维图像坐标(x2,y2)。
在本公开的一些实施例中,根据左图像上二维图像坐标(x1,y1)和右图像上的二维图像坐标(x2,y2)中的列向坐标x1,x2、左图像对应的飞机速度和右图像对应的飞机速度左图像对应的飞机三维位置和右图像对应的飞机三维位置计算变换参数。
在本公开的一些实施例中,所述变换参数为:
式中,R0为近地点斜距,Mx为斜距的分辨率,R0,Mx均为机载SAR的已知参数,x1为提取的左图像的列方向坐标值,x2为提取的右图像的列方向坐标值。
在本公开的一些实施例中,所述地面点的三维位置为:
(三)有益效果
从上述技术方案可以看出,本公开基于机载SAR立体像对计算地面点三维位置的方法至少具有以下有益效果其中之一:
(1)本公开基于机载SAR立体像对计算地面点三维位置的方法,通过机载SAR立体像对的平台参数和系统参数直接显式地计算地面点的三维位置,无需地面点的位置初始值,计算中也没有任何的迭代和近似过程,提高了计算效率。
(2)本公开基于机载SAR立体像对计算地面点三维位置的方法,通过SAR立体像对严密的成像几何关系进行计算,没有进行任何逼近和替换,保证了三维位置计算结果的严密性。
附图说明
通过附图所示,本公开的上述及其它目的、特征和优势将更加清晰。在全部附图中相同的附图标记指示相同的部分,并未刻意按实际尺寸等比例缩放绘制附图,重点在于示出本公开的主旨。
图1为根据本公开实施例基于机载SAR立体像对计算地面点三维位置的方法的流程图。
具体实施方式
为使本公开的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本公开进一步详细说明。
需要说明的是,在附图或说明书描述中,相似或相同的部分都使用相同的图号。附图中未绘示或描述的实现方式,为所属技术领域中普通技术人员所知的形式。另外,虽然本文可提供包含特定值的参数的示范,但应了解,参数无需确切等于相应的值,而是可在可接受的误差容限或设计约束内近似于相应的值。实施例中提到的方向用语,例如“上”、“下”、“前”、“后”、“左”、“右”等,仅是参考附图的方向。因此,使用的方向用语是用来说明并非用来限制本公开的保护范围。
本公开提供了一种基于机载SAR立体像对计算地面点三维位置的方法,通过机载SAR立体像对的平台参数和系统参数直接显式地计算地面点的三维位置,无需地面点的位置初始值,计算中也无任何的迭代和近似过程,提高了计算效率,保证了三维位置计算结果的严密性。
图1为根据本公开实施例基于机载SAR立体像对计算地面点三维位置的方法的流程图。如图1所示,所述基于机载SAR立体像对计算地面点三维位置的方法,包括:
S1,提取地面点在SAR立体像对上的图像坐标及飞机速度,根据所述图像坐标确定获取立体像对时对应的飞机三维位置;
S2,根据所述飞机三维位置计算基线参数;
S3,根据所述图像坐标、飞机三维位置及飞机速度计算变换参数;以及
S4,根据所述基线参数及变换参数计算地面点的三维位置。
以下进一步详细介绍本公开基于机载SAR立体像对计算地面点三维位置的方法的具体步骤:
在步骤S1中,利用图像匹配软件或人工的方法,提取地面点在SAR立体像对的左图像上二维图像坐标(x1,y1)和右图像上的二维图像坐标(x2,y2);左图像对应的飞机速度从飞机上GPS记录的辅助数据中读取,右图像对应的飞机速度也从飞机上GPS记录的辅助数据中读取;并根据飞机平台上GPS记录的位置数据和图像行方向坐标值y1、y2进行拟合得到左图像所对应的飞机三维位置和右图像所对应的飞机三维位置
在步骤S2中,根据所述飞机三维位置计算基线参数,所述立体像对的基线参数为:
在步骤S3中,根据图像坐标、立体图像对应的飞机三维位置和飞机速度计算变换参数,所述变换参数的计算公式如下:
其中,P1,P2,P3为计算的三个变换参数,是左图像对应的飞机速度,可以从通过飞机上GPS记录的辅助数据中读取,是右图像对应的飞机速度,可以从通过飞机上GPS记录的辅助数据中读取,左幅图像对应的飞机三维位置,右幅图像对应的飞机三维位置,R0为近地点斜距,Mx为斜距的分辨率,这两个都是机载SAR的已知参数,x1为提取的左幅SAR图像的列方向坐标值,x2为提取的右幅SAR图像的列方向坐标值。
步骤S4,利用基线参数和变换参数计算地面点的三维位置,计算公式如下,
其中,(X、Y、Z)表示计算的地面点三维位置,ΔXa,ΔYa,ΔZa表示SAR立体像对的三
个基线参数,由步骤S1计算得到,P1,P2,P3为三个变换参数,由步骤S2计算得到,是左图像对应的飞机速度,从飞机上GPS记录的辅助数据中读取,是右图像对应的飞机速度,也从飞机上GPS记录的辅助数据中读取。
至此,已经结合附图对本实施例进行了详细描述。依据以上描述,本领域技术人员应当对本公开利用机载SAR立体像对计算地面点三维位置的方法有了清楚的认识。
此外,上述对各元件和方法的定义并不仅限于实施方式中提到的各种具体结构、形状或方式,本领域的普通技术人员可对其进行简单地熟知地替换。
综上所述,本公开利用机载SAR立体像对计算地面点三维位置的方法应用了获取SAR立体像对的飞机平台参数和系统参数直接显式地计算地面点的三维位置信息,无需给出任何初值,计算方法也是采用了严密的成像几何关系,没有任何迭代和近似,是十分严密的计算过程。
类似地,应当理解,为了精简本公开并帮助理解各个公开方面中的一个或多个,在上面对本公开的示例性实施例的描述中,本公开的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而,并不应将该公开的方法解释成反映如下意图:即所要求保护的本公开要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说,如下面的权利要求书所反映的那样,公开方面在于少于前面公开的单个实施例的所有特征。因此,遵循具体实施方式的权利要求书由此明确地并入该具体实施方式,其中每个权利要求本身都作为本公开的单独实施例。
以上所述的具体实施例,对本公开的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本公开的具体实施例而已,并不用于限制本公开,凡在本公开的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本公开的保护范围之内。
Claims (4)
1.一种基于机载SAR立体像对计算地面点三维位置的方法,包括以下步骤:
提取地面点在SAR立体像对上的图像坐标及飞机速度,根据所述图像坐标确定获取立体像对时对应的飞机三维位置;
根据所述飞机三维位置计算基线参数;
根据所述图像坐标、飞机三维位置及飞机速度计算变换参数;以及
根据所述基线参数及变换参数计算地面点的三维位置;
所述提取地面点在SAR立体像对上的图像坐标及飞机速度,根据所述图像坐标确定获取立体像对时对应的飞机三维位置的步骤,包括:
提取地面点在SAR立体像对的左图像上的二维图像坐标(x1,y1)及右图像上的二维图像坐标(x2,y2);
根据飞机平台上GPS记录的位置数据和图像坐标拟合得到左图像所对应的飞机三维位置和右图像所对应的飞机三维位置
所述基线参数为:
根据左图像上二维图像坐标(x1,y1)和右图像上的二维图像坐标(x2,y2)中的列方向坐标值x1,x2、左图像对应的飞机速度和右图像对应的飞机速度左图像对应的飞机三维位置和右图像对应的飞机三维位置计算变换参数;
所述变换参数为:
式中,R0为近地点斜距,Mx为斜距的分辨率,R0,Mx均为机载SAR的已知参数,x1为提取的左图像的列方向坐标值,x2为提取的右图像的列方向坐标值;
所述地面点的三维位置为:
2.根据权利要求1所述的基于机载SAR立体像对计算地面点三维位置的方法,其中,利用图像匹配软件或人工方法,提取地面点在SAR立体像对的左图像上的二维图像坐标(x1,y1)和右图像上的二维图像坐标(x2,y2)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711002479.2A CN107765245B (zh) | 2017-10-24 | 2017-10-24 | 基于机载sar立体像对计算地面点三维位置的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711002479.2A CN107765245B (zh) | 2017-10-24 | 2017-10-24 | 基于机载sar立体像对计算地面点三维位置的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107765245A CN107765245A (zh) | 2018-03-06 |
CN107765245B true CN107765245B (zh) | 2020-01-21 |
Family
ID=61270115
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711002479.2A Active CN107765245B (zh) | 2017-10-24 | 2017-10-24 | 基于机载sar立体像对计算地面点三维位置的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107765245B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887122A (zh) * | 2010-06-29 | 2010-11-17 | 上海大学 | 可消除地面高程误差的星载sar图像目标定位方法 |
CN104237887A (zh) * | 2014-09-29 | 2014-12-24 | 中国测绘科学研究院 | 一种sar遥感影像匹配方法 |
CN104865573A (zh) * | 2015-06-12 | 2015-08-26 | 西安电子科技大学 | 一种双基前视高机动平台sar扩展场景成像优化方法 |
CN105913074A (zh) * | 2016-04-05 | 2016-08-31 | 西安电子科技大学 | 基于幅度与径向速度联合的sar图像动目标聚类方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8816896B2 (en) * | 2012-05-11 | 2014-08-26 | Raytheon Company | On-board INS quadratic correction method using maximum likelihood motion estimation of ground scatterers from radar data |
-
2017
- 2017-10-24 CN CN201711002479.2A patent/CN107765245B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887122A (zh) * | 2010-06-29 | 2010-11-17 | 上海大学 | 可消除地面高程误差的星载sar图像目标定位方法 |
CN104237887A (zh) * | 2014-09-29 | 2014-12-24 | 中国测绘科学研究院 | 一种sar遥感影像匹配方法 |
CN104865573A (zh) * | 2015-06-12 | 2015-08-26 | 西安电子科技大学 | 一种双基前视高机动平台sar扩展场景成像优化方法 |
CN105913074A (zh) * | 2016-04-05 | 2016-08-31 | 西安电子科技大学 | 基于幅度与径向速度联合的sar图像动目标聚类方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107765245A (zh) | 2018-03-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108648240B (zh) | 基于点云特征地图配准的无重叠视场相机姿态标定方法 | |
CN107316325B (zh) | 一种基于图像配准的机载激光点云与影像配准融合方法 | |
CN107180215B (zh) | 基于库位和二维码的停车场自动建图与高精度定位方法 | |
US8660338B2 (en) | Wide baseline feature matching using collobrative navigation and digital terrain elevation data constraints | |
EP3678095B1 (en) | Determination of position from images and associated camera positions | |
CN105300362B (zh) | 一种应用于rtk接收机的摄影测量方法 | |
US9959625B2 (en) | Method for fast camera pose refinement for wide area motion imagery | |
CN111815765B (zh) | 一种基于异构数据融合的图像三维重建方法 | |
CN107564057B (zh) | 顾及大气折光校正的高轨面阵光学卫星在轨几何标定方法 | |
CN103278138A (zh) | 一种复杂结构薄部件三维位置及姿态的测量方法 | |
JP2009501996A (ja) | 映像の幾何補正方法およびその装置 | |
CN107024687B (zh) | 一种离线快速实现pos/激光雷达安装误差标定的方法 | |
CN108759788B (zh) | 无人机影像定位定姿方法及无人机 | |
WO2022179094A1 (zh) | 车载激光雷达外参数联合标定方法、系统、介质及设备 | |
CN115187798A (zh) | 一种多无人机高精度匹配定位方法 | |
Lee | Assessing the positioning accuracy of high density point clouds produced from rotary wing quadrocopter unmanned aerial system based imagery | |
CN117253029B (zh) | 基于深度学习的图像匹配定位方法及计算机设备 | |
CN115423863B (zh) | 相机位姿估计方法、装置及计算机可读存储介质 | |
Liu et al. | A new approach to fast mosaic UAV images | |
CN116184430A (zh) | 一种激光雷达、可见光相机、惯性测量单元融合的位姿估计算法 | |
KR100520275B1 (ko) | Gps 측량에 의한 수치지도 제작시 촬영기구회전모델을 이용한 선형 스캐닝 영상의 기하 보정방법 | |
CN110030968B (zh) | 一种基于星载立体光学影像的地面遮挡物仰角测量方法 | |
Habib et al. | Linear features in photogrammetric activities | |
CN107765245B (zh) | 基于机载sar立体像对计算地面点三维位置的方法 | |
CN110503604B (zh) | 一种基于高精度pos的航空面阵影像实时正射拼接方法 |
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 |