CN106249238A - 获取地面目标点高度的方法 - Google Patents
获取地面目标点高度的方法 Download PDFInfo
- Publication number
- CN106249238A CN106249238A CN201610866143.XA CN201610866143A CN106249238A CN 106249238 A CN106249238 A CN 106249238A CN 201610866143 A CN201610866143 A CN 201610866143A CN 106249238 A CN106249238 A CN 106249238A
- Authority
- CN
- China
- Prior art keywords
- impact point
- height
- insar
- height value
- 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
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
- 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/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- 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
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)
- Signal Processing (AREA)
- Image Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种获取地面目标点高度的方法,包括:获取两幅星载合成孔径雷达图像数据,采用干涉处理方法得到一目标点的高度值HInSAR;步骤B:根据所述两幅星载合成孔径雷达图像数据,采用立体匹配处理得到同一目标点的高度值HStereoSAR;以及利用干涉处理的目标点高度值和立体匹配得到的目标点高度值进行加权组合计算,获得该目标点的高度H目标。本发明的方法实现了在任意间隔距离和任意时间间隔利用星载SAR图像有效计算地面目标点的高度。
Description
技术领域
本发明涉及星载合成孔径雷达遥感图像处理技术领域,尤其涉及一种获取地面目标点高度的方法。
背景技术
星载合成孔径雷达(英文简称SAR,下文都用SAR表示)遥感是通过装载在卫星上的雷达基于合成孔径的技术来获取地面的高分辨率图像,由于SAR图像具有全天时、全天候的优势,因此在国民经济建设和国防安全都能发挥了不可替代的作用,目前国内外都在大力发展星载SAR遥感技术。
星载干涉合成孔径雷达(英文简称InSAR)通过卫星上的SAR天线获取同一地区的两幅SAR图像数据,基于无线电干涉原理,经配准、干涉图生成与相位滤波、相位解缠等处理就可以获取地面目标点的高度信息,具有精度高的特点,但星载InSAR要求两幅SAR数据具有很好的相干性,因此要求获取两幅SAR数据的天线位置距离比较近,而且两幅SAR数据的获取时间也较短。
星载立体SAR通过星上的SAR天线获取同一地区的两幅SAR图像,基于立体视觉原理,经图像同名点匹配、高度解算而获取地面目标点的高度值,具有应用灵活适应性强、且能够满足一定的精度要求等特点。但是星载立体SAR为了获取较高的精度,常常要求获取两幅SAR数据的天线位置距离比较远。
综合起来,星载干涉SAR和星载立体SAR都是利用了同一地区的两幅图像,而采用了不同的原理来获得地面目标点的高度值。但是实际应用中,为了采用InSAR来获取较高精度的高度值,要求获取两幅SAR数据的天线位置相距比较近、时间相差比较短;而星载立体SAR则要求获取两幅SAR数据的天线位置相距比较远、而对时间没有任何要求。
发明内容
(一)要解决的技术问题
鉴于上述技术问题,本发明提供了一种获得地面目标点高度的方法。该方法依据获取的同一地区的两幅星载SAR图像数据,有效地将InSAR计算的高度和立体SAR计算的高度值进行组合,解决不同条件下地面目标点的高度计算问题。
(二)技术方案
本发明提供一种获取地面目标点高度的方法,包括:
步骤A:获取两幅星载合成孔径雷达图像数据,采用干涉处理方法得到一目标点的高度值HInSAR;
步骤B:根据所述两幅星载合成孔径雷达图像数据,采用立体匹配处理得到同一目标点的高度值HStereoSAR;以及
步骤C:利用干涉处理的目标点高度值和立体匹配得到的目标点高度值进行加权组合计算,获得该目标点的高度H目标。
进一步的,在步骤C中,加权组合计算的具体公式如下:
如果HInSAR等于0或者HStereoSAR等于0,则H目标=max(HInSAR,HStereoSAR);
如果HInSAR不等于0且HStereoSAR不等于0,则
其中,H目标表示计算的目标点的高度值;max()表示从括号内取其中最大的值;HInSAR表示步骤A计算的目标点的高度值;HStereoSAR表示步骤B计算的目标点的高度值;mx表示SAR图像的分辨率;λ表示获取星载合成孔径雷达数据所使用的无线电波长。
进一步的,步骤A中获取两幅星载SAR图像数据时的天线位置之间的距离不受限制。
进一步的,步骤A中获取两幅星载图像数据时的时间间隔为任意时间间隔。
进一步的,在步骤A中干涉处理方法包括:对两幅SAR数据配准、干涉条纹生成与滤波、干涉相位展开和平地效应去除过程。
进一步的,步骤B中立体匹配处理包括:图像配准、同名点提取和立体解算过程。
(三)有益效果
从上述技术方案可以看出,本发明的获取地面目标点高度方法具有以下有益效果:
(1)根据同一地区的两幅SAR图像数据,有效组合地了InSAR和立体SAR的计算目标高度的各自优势,可以很好地适应各种距离条件下的地面目标高度计算;
(2)计算地面目标点的高度中充分考虑InSAR和立体SAR各自的精度特性,进行了加权组合,较好地兼顾了各自高度计算值的精度。
附图说明
图1为根据本发明实施例组合InSAR和立体SAR计算地面目标点高度方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。需要说明的是,在附图或说明书描述中,相似或相同的部分都使用相同的图号。
附图中未绘示或描述的实现方式,为所属技术领域中普通技术人员所知的形式。另外,虽然本文可提供包含特定值的参数的示范,但应了解,参数无需确切等于相应的值,而是可在可接受的误差容限或设计约束内近似于相应的值。
星载合成孔径雷达在实际应用中总能获取同一地区的两幅SAR图像数据,但是很难单独满足InSAR或立体SAR的距离要求,为此实际应用中可以将星载InSAR和立体SAR有效组合,从而能够在各种条件下都可以有效计算目标点的高度。
本发明获取地面目标高度的方法(通过组合InSAR和立体SAR计算),充分考虑了星载InSAR和立体SAR的计算目标高度的各自优势,计算地面目标点的高度中时也兼顾了各自高度计算值的精度,从而很好地适应星载SAR各种距离条件下的地面目标高度计算。
在本发明的一个示例性实施例中,提供了一种组合InSAR和立体SAR计算地面目标点高度的方法。图1为根据本发明实施例组合InSAR和立体SAR计算地面目标点高度的方法的流程图。如图1所示,本实施例获取地面目标点高度的方法包括:
步骤A:利用两幅SAR数据采用成熟的干涉处理,包括SAR数据配准、干涉条纹生成与滤波、干涉相位展开和平地效应去除等过程,得到某个目标点的高度值HInSAR。如果利用两幅SAR数据进行干涉处理而不能计算出该目标点的高度值,则该目标点的高度值设为0,即HInSAR=0;
在该步骤中,获取两幅星载SAR图像数据时的天线位置之间的距离不受限制,具体是指天线可以在运行轨道的任意位置进行图像数据的获取,两幅图像获取时天线可以位于卫星轨道上的任意位置,但两个位置不应该重合。
另外在该步骤中,获取两幅星载图像数据时的时间间隔为任意时间间隔,具体是指天线可以在运行过程中的任意时刻进行图像数据的获取,两幅图像获取可以间隔任意时间段。
步骤B:利用两幅SAR数据采用成熟的立体匹配处理,包括图像配准、同名点提取、立体解算等过程,得到和步骤A相同目标点的高度值HStereoSAR,如果利用两幅SAR图像进行匹配和立体处理而不能计算出该目标点的高度值,则该目标点的高度值也设为0,即HStereoSAR=0;
步骤C:利用干涉处理的目标点高度值和立体匹配得到的目标点高度值进行组合计算,完成该目标点的高度计算。计算的具体公式如下
(1)如果HInSAR等于0或者HStereoSAR等于0,则H目标=max(HInSAR,HStereoSAR)
(2)如果HInSAR不等于0且HStereoSAR不等于0,则
其中,H目标表示计算的目标点的高度值;max()表示从括号内取其中最大的值;HInSAR表示步骤A计算的目标点的高度值;HStereoSAR表示步骤B计算的目标点的高度值;mx表示SAR图像的分辨率,是SAR系统已知的参数,如0.5m;λ表示获取SAR数据所使用的无线电波长,也是已知的参数,比如0.19m。
至此,已经结合附图对本实施例进行了详细描述。依据以上描述,本领域技术人员应当对本发明组合InSAR和立体SAR计算地面目标点高度的方法有了清楚的认识。
此外,上述对各元件和方法的定义并不仅限于实施方式中提到的各种具体结构、形状或方式,本领域的普通技术人员可对其进行简单地熟知地替换。
综上所述,本发明获取地面目标点高度的方法中会根据自动根据同一地区的两幅SAR图像数据,有效组合地InSAR和立体SAR所各自计算的目标高度值,并能兼顾InSAR和立体SAR各自的精度特性进行加权,从而可以适应星载SAR的各种距离条件下进行地面目标点高度的计算。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种获取地面目标点高度的方法,其特征在于包括:
步骤A:获取两幅星载合成孔径雷达图像数据,采用干涉处理方法得到一目标点的高度值HInSAR;
步骤B:根据所述两幅星载合成孔径雷达图像数据,采用立体匹配处理得到同一目标点的高度值HStereoSAR;以及
步骤C:利用干涉处理的目标点高度值和立体匹配得到的目标点高度值进行加权组合计算,获得该目标点的高度H目标。
2.根据权利要求1所述的方法,其特征在于,步骤C中,加权组合计算的具体公式如下:
如果HInSAR等于0或者HStereoSAR等于0,则H目标=max(HInSAR,HStereoSAR);
如果HInSAR不等于0且HStereoSAR不等于0,则
其中,H目标表示计算的目标点的高度值;max()表示从括号内取其中最大的值;HInSAR表示步骤A计算的目标点的高度值;HStereoSAR表示步骤B计算的目标点的高度值;mx表示SAR图像的分辨率;λ表示获取星载合成孔径雷达数据所使用的无线电波长。
3.根据权利要求1所述的方法,其特征在于,步骤A中获取两幅星载SAR图像数据时的天线位置之间的距离不受限制。
4.根据权利要求1所述的方法,其特征在于,步骤A中获取两幅星载图像数据时的时间间隔为任意时间间隔。
5.根据权利要求1所述的方法,其特征在于,步骤A中干涉处理方法包括:对两幅SAR数据配准、干涉条纹生成与滤波、干涉相位展开和平地效应去除过程。
6.根据权利要求1所述的方法,其特征在于,步骤B中立体匹配处理包括:图像配准、同名点提取和立体解算过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610866143.XA CN106249238B (zh) | 2016-09-29 | 2016-09-29 | 获取地面目标点高度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610866143.XA CN106249238B (zh) | 2016-09-29 | 2016-09-29 | 获取地面目标点高度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106249238A true CN106249238A (zh) | 2016-12-21 |
CN106249238B CN106249238B (zh) | 2019-02-15 |
Family
ID=57612049
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610866143.XA Active CN106249238B (zh) | 2016-09-29 | 2016-09-29 | 获取地面目标点高度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106249238B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108761444A (zh) * | 2018-05-24 | 2018-11-06 | 中国科学院电子学研究所 | 联合星载sar和光学图像计算地面点高度的方法 |
CN108983239A (zh) * | 2018-08-08 | 2018-12-11 | 中国科学院电子学研究所 | 星载干涉sar数字高程模型重建方法 |
CN109000770A (zh) * | 2018-06-20 | 2018-12-14 | 威马智慧出行科技(上海)有限公司 | 一种车辆超载检测方法及其装置 |
CN111597496A (zh) * | 2020-05-18 | 2020-08-28 | 中国科学院空天信息创新研究院 | 基于有理多项式系数单幅星载遥感影像目标高度计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101339245A (zh) * | 2008-08-08 | 2009-01-07 | 西安电子科技大学 | 多基线干涉合成孔径雷达干涉相位展开方法 |
CN104515980A (zh) * | 2014-12-08 | 2015-04-15 | 广西大学 | 一种基于InSAR构型的地面运动目标检测方法和装置 |
-
2016
- 2016-09-29 CN CN201610866143.XA patent/CN106249238B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101339245A (zh) * | 2008-08-08 | 2009-01-07 | 西安电子科技大学 | 多基线干涉合成孔径雷达干涉相位展开方法 |
CN104515980A (zh) * | 2014-12-08 | 2015-04-15 | 广西大学 | 一种基于InSAR构型的地面运动目标检测方法和装置 |
Non-Patent Citations (2)
Title |
---|
王峰 等: ""应用于SAR图像配准的自适应SIFT特征均匀分布算法"", 《武汉大学学报 信息科学版》 * |
韩松涛 等: ""一种基于特征点权重的机载InSAR系统区域网干涉参数定标方法"", 《电子与信息学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108761444A (zh) * | 2018-05-24 | 2018-11-06 | 中国科学院电子学研究所 | 联合星载sar和光学图像计算地面点高度的方法 |
CN108761444B (zh) * | 2018-05-24 | 2021-12-21 | 中国科学院电子学研究所 | 联合星载sar和光学图像计算地面点高度的方法 |
CN109000770A (zh) * | 2018-06-20 | 2018-12-14 | 威马智慧出行科技(上海)有限公司 | 一种车辆超载检测方法及其装置 |
CN108983239A (zh) * | 2018-08-08 | 2018-12-11 | 中国科学院电子学研究所 | 星载干涉sar数字高程模型重建方法 |
CN108983239B (zh) * | 2018-08-08 | 2022-01-04 | 中国科学院电子学研究所 | 星载干涉sar数字高程模型重建方法 |
CN111597496A (zh) * | 2020-05-18 | 2020-08-28 | 中国科学院空天信息创新研究院 | 基于有理多项式系数单幅星载遥感影像目标高度计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106249238B (zh) | 2019-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106249238A (zh) | 获取地面目标点高度的方法 | |
JP6106285B2 (ja) | 衛星レーダー干渉度のイオン歪み補正方法及びその装置 | |
CN107102333B (zh) | 一种星载InSAR长短基线融合解缠方法 | |
WO2016086699A1 (zh) | 一种结合局部频率估计的小波域InSAR干涉相位滤波方法 | |
ITTO20110526A1 (it) | Acquisizione di immagini sar per calcolare una quota o un modello digitale di elevazione tramite elaborazioni interferometriche | |
Yague-Martinez et al. | TanDEM-X mission: SAR image coregistration aspects | |
CN102411778A (zh) | 一种机载激光点云与航空影像的自动配准方法 | |
US20170186164A1 (en) | Method for fast camera pose refinement for wide area motion imagery | |
CN104730519B (zh) | 一种采用误差迭代补偿的高精度相位解缠方法 | |
KR101712084B1 (ko) | 멀티간섭기법 기반의 이온왜곡 보정 방법 및 그 장치 | |
CN108132468A (zh) | 一种多基线极化干涉sar建筑物高度提取方法 | |
CN112797886B (zh) | 面向缠绕相位的InSAR时序三维形变监测方法 | |
CN109146940B (zh) | 一种基于权重相位相关峰值拟合的高精度图像配准方法 | |
CN105301588A (zh) | 一种联合StereoSAR和InSAR的DEM提取方法 | |
CN103197316A (zh) | 一种保相的机载干涉sar地形高度自适应运动补偿方法 | |
CN103630898A (zh) | 对多基线干涉sar相位偏置进行估计的方法 | |
Oh et al. | Automatic georeferencing of aerial images using stereo high-resolution satellite images | |
Cefalu et al. | Structureless bundle adjustment with self-calibration using accumulated constraints | |
Mostafa et al. | Optical flow based approach for vision aided inertial navigation using regression trees | |
Choi et al. | Federated‐filter‐based unmanned ground vehicle localization using 3D range registration with digital elevation model in outdoor environments | |
WO2018139297A1 (ja) | カメラ装置 | |
CN115015931A (zh) | 无需外部误差校正的实时差分立体sar几何定位方法及系统 | |
CN111239735B (zh) | 一种基于低频sar影像的方位向形变场电离层校正方法 | |
KR101315243B1 (ko) | 위성레이더 영상의 도플러 중심주파수 변환방법 | |
CN103323844A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |