CN110018476B - 时间差分基线集时序干涉sar处理方法 - Google Patents
时间差分基线集时序干涉sar处理方法 Download PDFInfo
- Publication number
- CN110018476B CN110018476B CN201910421195.XA CN201910421195A CN110018476B CN 110018476 B CN110018476 B CN 110018476B CN 201910421195 A CN201910421195 A CN 201910421195A CN 110018476 B CN110018476 B CN 110018476B
- Authority
- CN
- China
- Prior art keywords
- phase
- deformation
- processing
- interference
- elevation
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 51
- 238000000034 method Methods 0.000 claims abstract description 33
- 238000012937 correction Methods 0.000 claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 230000001427 coherent effect Effects 0.000 claims description 9
- 238000010276 construction Methods 0.000 claims description 4
- 238000004804 winding Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000000611 regression analysis Methods 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 description 14
- 230000002123 temporal effect Effects 0.000 description 7
- 230000002596 correlated effect Effects 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000010587 phase diagram Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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/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/9017—SAR image acquisition techniques with time domain processing of the SAR signals in azimuth
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)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种时间差分基线集时序干涉SAR处理方法,其特征在于包括:对(N+1)景已配准的SAR影像进行时间差分干涉处理形成干涉对;利用外部高程模型确定处理区域地形相位和平地相位,用于去除处理区域的地形和平地相位;根据已形成的干涉对确定空间基线;通过一维回归方式确定高程校正;通过一维回归方式确定形变;获取大气相位和非线性形变,经高程校正和形变解算,确定线性形变、非线性形变、大气相位和高程校正信号量。本发明方法中利用临近影像,有效避免了时间失相干的影响,提高了处理精度;通过将二维回归分析分解为两个一维回归分析,极大提高了运算效率,降低了运算复杂度。
Description
技术领域
本发明涉及合成孔径雷达处理技术,进一步涉及一种时间差分基线集时序干涉SAR处理方法。
背景技术
近几十年来,随着时间序列干涉合成孔径雷达技术(Time SeriesInterferometric Synthetic Aperture Radar,TSInSAR)的提出,大大扩展了InSAR的应用,该方法通过识别场景中的高相干散射体,例如房屋、高楼、道路等人工建筑及裸露的岩石等,仅在这些散射特性非常稳定的相干散射点(Coherent Scatterers,CS)上进行形变信号、大气延时、高程误差等信号的提取,在一定程度上克服了时间失相干与空间失相干的影响。但是在非城区该技术可识别到的CS点非常有限,且无法保证足够的形变监测精度,严重影响了该技术应用范围。随后,为了扩展该技术在非城区的应用,发展出了第二代时序干涉SAR技术,通过一套自适应算法,在不影响城区CS点识别的基础上,同时能够识别出分布散射体(Distributed Scatterers,DS),从而极大提升了CS点的空间密度和相干性,最后根据传统的时序干涉SAR技术求解形变、高程误差等信息。图1给出了时序干涉SAR技术的处理流程。但是该时序干涉SAR技术的缺点就是监测点密度低,算法复杂度高,计算量大,效率偏低。为了进一步抑制时间失相干与空间失相干的影响,人们提出了小基线集技术,图2给出了处理流程,该技术将所有时间基线和空间基线在一定阈值内的影像进行干涉组合,得到一个或几个干涉对集合,进行相位解缠后通过奇异值分解(Singular ValueDecomposition,SVD)方法求出时间序列形变,但是该方法需要牺牲分辨率并且需要相位解缠。
发明内容
(一)要解决的技术问题
有鉴于此,本发明的目的在于提供一种时间差分基线集时序干涉SAR处理方法,以解决以上所述的至少部分技术问题。
(二)技术方案
为实现上述目的,本发明提供一种时间差分基线集时序干涉SAR处理方法,包括:
对(N+1)景已配准的SAR影像进行时间差分干涉处理形成干涉对;
利用外部高程模型确定处理区域地形相位和平地相位,用于去除处理区域的地形和平地相位;
根据已形成的干涉对确定空间基线;
通过一维回归方式确定高程校正;
通过一维回归方式确定形变;
获取大气相位和非线性形变,经高程校正和形变解算,确定线性形变、非线性形变、大气相位和高程校正信号量。
在进一步的实施方案中,(N+1)景已配准的SAR影像确定方式包括:在(N+1)景SAR影像中选取一景影像作为主影像,其他N景影像作为从影像与分别与主影像进行配准,最终得到N+1景已经配准的SAR影像;通过对(N+1)景已经配准的SAR影像进行处理,获取处理区域的高相干散射点。
在进一步的实施方案中,所述外部高程模型通过以下至少一种方式获取:引入外部DEM,通过地理编码转换到SAR坐标系;任选SAR影像中两景影像数据,通过InSAR处理获取。
在进一步的实施方案中,对(N+1)景已配准的SAR影像进行时间差分干涉处理形成干涉对,包括:对SAR影像进行二轨差分处理,确定处理区域的形变信息。
在进一步的实施方案中,对(N+1)景已配准的SAR影像进行时间差分干涉处理形成干涉对,包括:将(N+1)幅已经配准的影像按时间顺序排列为:{φ0,φ1,φ2,...,φN};相邻时间的影像组合进行干涉,生成N幅干涉图{δφ1,δφ2,...,δφN};确定相干散射点的干涉相位,表示如下:δφk=φk-φk-1,k=1,2,...,N;确定相位信号,表示如下:
其中,δφk,i表示第k幅干涉图的第i个点的干涉相位,注意它是缠绕的,范围在-π到π之间,{·}2π表示缠绕操作,/>表示高程误差相位,由于不精确的外部数字高程模型(Digital Elevation Model,DEM)引起,/>表示形变信号相位,/>表示轨道误差相位,/>表示大气延迟相位,/>表示随机噪声相位。
在进一步的实施方案中,高程误差相位为:
所述形变相位包括线性形变项与非线性形变项,表示如下:
其中,δTk表示第k幅干涉图的时间基线,表示沿视线方向的线性形变速率,/>表示非线性形变相位。通常,非线性形变速率非常小,又注意到我们是通过邻近时间影像生成的干涉图,故时间基线δTk非常小,那么形变相位/>就很小。
在进一步的实施方案中,通过一维回归方式确定高程校正包括:
对空间上相近的高相干散射点(CS)做二阶差分,得到如下相位:
则有:
建立如下的目标函数:
在进一步的实施方案中,一维回归方式确定形变,包括:获取线性形变速率,经过高程校正后,得到包含形变信息的残余相位积分后如下式:
(三)有益效果
本发明的方法是一项全新的时序干涉SAR处理技术,它完全不同于传统的时序干涉SAR技术。从影像配对方式到参数求解方法与当前已有技术都不同,是一套完整的处理方法,从理论分析与实际数据处理上验证了该技术的可行性和有效性;
本发明方法中利用临近影像,有效避免了时间失相干的影响,提高了处理精度;
本发明中,通过将二维回归分析分解为两个一维回归分析,极大提高了运算效率,降低了运算复杂度;
本发明的方法简化了数据处理流程,提升了处理效率。
附图说明
图1是现有技术的一种时序干涉SAR处理方法流程图。
图2是现有技术的另一种时序干涉SAR处理方法流程图。
图3是本发明实施例的时间差分基线集时序干涉SAR处理方法流程图。
图4A,图4B和图4C是本发明实施例的图1、图2和图3三种方法形成的干涉图生成方式对比图。
图5是实施例1的线性形变速率图。
图6是实施例1的高程校正图。
图7是实施例1的非线性形变结果图。
图8是实施例1的形变序列图。
图9是实施例1的大气相位示意图。
图10是实施例1的任意两个相干散射点的形变时序图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明作进一步的详细说明。
随着星载SAR技术的发展,可以获得越来越多的SAR图像,并且其图像质量也在逐步提高,一方面,其控轨精度不断提高,避免了以往信号处理遭遇的空间失相干问题。另一方面,重访周期逐步缩短(最快4天),所以可以在短时间内获取大量的同一地区的影像进行形变监测,不用担心影像数量不够多而对处理精度产生影响。
为了克服现有技术的缺点,本发明实施例提出的时间差分集时序干涉SAR处理方法,该方法创造性地采用相邻时间的影像进行配对生成干涉图,称为时间差分基线集(TimeDifference Baseline Set,TDBS)。一方面,由于TDBS只采用时间最邻近的影像进行干涉图的生成,这样显然保证了所有干涉图都是保持最短时间基线,有效避免了时间失相干的影响;另一方面,TDBS技术将原本需要在高程与形变上进行的二维空间求解的问题转化为分别在高程上与形变速率上的一维求解问题,大大提高了计算效率。最后通过时空滤波分离出非线性形变与大气延迟。
本发明的总体方案是利用临近时间影像干涉去除时间去相干误差,并将常规的二维回归分析转换为两个一维回归分析降低运算量,提升运算效率,简化处理流程。
如图3和4所示,本发明实施例提供一种时间差分集时序干涉SAR处理方法,其中可以包括如下步骤:
第一步:在N+1幅SAR影像中选取一景影像作为主影像,其他N景影像作为从影像与分别与主影像进行配准,最终得到N+1景已经配准的SAR影像。
第二步:获取处理区域高精度DEM,用于处理区域地形相位和平地相位去除。有两种方式:第一种是引入外部DEM,通过地理编码转换到SAR坐标系。第二种是任选两景影像数据,通过InSAR处理获取。
第三步:对影像进行二轨差分处理,获取处理区域厘米级形变信息。
第四步:通过对N+1景已经配准的SAR影像进行处理,获取该地区的高相干散射点(CS)。
第五步:针对N+1幅已经配准的同一地区不同时间获取的SAR影像,按时间顺序排列为:
{φ0,φ1,φ2,...,φN} (1)
相邻时间的影像组合进行干涉,可以生成N幅干涉图:
{δφ1,δφ2,...,δφN} (2)
相干散射点(CS)的干涉相位表示如下:
δφk=φk-φk-1,k=1,2,...,N (3)
上面这些步骤完成以后,将得到的相位信号表示如下:
其中,δφk,i表示第k幅干涉图的第i个点的干涉相位,注意它是缠绕的,范围在-π到π之间,{·}2π表示缠绕操作,表示高程误差相位,由于不精确的外部数字高程模型(Digital Elevation Model,DEM)引起,/>表示形变信号相位,/>表示轨道误差相位,/>表示大气延迟相位,/>表示随机噪声相位。
下面分析一下每个相位分量。
对于高程误差相位经过推导可以得到:
形变相位是在时间和空间都具有相关性的信号。常将它分解为线性形变项与非线性形变项:
其中,δTk表示第k幅干涉图的时间基线,表示沿视线方向的线性形变速率,/>表示非线性形变相位。通常,非线性形变速率非常小,又注意到是通过邻近时间影像生成的干涉图,故时间基线δTk非常小,那么形变相位/>就很小。
第六步:针对已生成的干涉对估计空间基线。
第七步:对空间上相近的高相干散射点(CS)做二阶差分,得到如下相位:
所以可以认为下式成立:
有:
因此,通过建立如下的目标函数:
对上式优化即可得到任意两点之间的高程误差值Δεi,j,然后在空间上构网进行平差处理得到任意一点的高程误差值Δεi。
第八步:获取线性形变速率,经过高程校正后,得到了包含形变信息的残余相位。为了从中求解出形变速率,对残余相位做如下积分操作:
其中:
由于线性形变在整个时间范围上是相关的,故积分使其加强;而非线性形变、大气相位、噪声在时间上不相关,故积分操作操作不会改变它们的大小。
于是,满足:
与高程类似的,建立如下目标函数:
第九步:获取大气相位和非线性形变,经过高程校正与线性形变求解后,现在残余相位包含非线性形变、大气延迟、噪声:
大气延迟在空间上表现为相关的,同一个点的大气延迟在不同时间上是不相关的,非线性形变信号在时间和空间上显然都是相关的,噪声表现为随机的高斯白噪声。基于这几个特性,通过以下滤波方法分离出非线性形变与大气延时。
其中,HPT{·}表示时间上的高通滤波,LPS{·}表示空间上的低通滤波。
非线性形变:
第十步:通过二次差分获取了全部信号分量,包括毫米级精度的线性形变和非线性形变,大气相位,高程校正。
传统的TSInSAR技术,小基线集技术(SBAS),本发明实施例的(TDBS)的干涉图生成方式说明如图4A-4C所示:从结果看出,本发明实施例的方法即TDBS技术有效地求解出了形变信息。
实施例1
为了验证该方法的有效性,发明人进行了实际数据的处理验证。采用从2012年1月到2016年3月共32景北京通州地区的高分辨率卫星影像,利用TDBS技术形变解算,我们得到图5-10的形变速率结果图。如图所示,通过上述的本发明实施例的方法,获得了线性形变速率图、高程校正图、非线性形变结果图、形变序列图、大气相位示意图和相干散射点形变时序图。图中可以看出,所求出的形变信息数据大幅提高了处理精度。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种时间差分基线集时序干涉SAR处理方法,其特征在于包括:
对N+1景已配准的SAR影像进行时间差分干涉处理形成干涉对;
利用外部高程模型确定处理区域地形相位和平地相位,用于去除处理区域的地形和平地相位;
根据已形成的干涉对确定空间基线;
通过一维回归方式确定高程校正:
对空间上相近的高相干散射点CS做二阶差分,得到如下相位:
则有:
建立如下的目标函数:
通过一维回归方式确定形变;
获取大气相位和非线性形变,经高程校正和形变解算,确定线性形变、非线性形变、大气相位和高程校正信号量。
2.根据权利要求1所述的方法,其特征在于:
N+1景已配准的SAR影像确定方式包括:
在N+1景SAR影像中选取一景影像作为主影像,其他N景影像作为从影像分别与主影像进行配准,最终得到N+1景已经配准的SAR影像;
通过对N+1景已经配准的SAR影像进行处理,获取处理区域的高相干散射点。
3.根据权利要求1所述的方法,其特征在于,所述外部高程模型通过以下至少一种方式获取:
引入外部DEM,通过地理编码转换到SAR坐标系;
任选SAR影像中两景影像数据,通过InSAR处理获取。
4.根据权利要求2所述的方法,其特征在于,对N+1景已配准的SAR影像进行时间差分干涉处理形成干涉对,包括:
对SAR影像进行二轨差分处理,确定处理区域的形变信息。
5.根据权利要求1所述的方法,其特征在于,对N+1景已配准的SAR影像进行时间差分干涉处理形成干涉对,包括:
将N+1幅已经配准的影像按时间顺序排列为:
{φ0,φ1,φ2,...,φN};
相邻时间的影像组合进行干涉,生成N幅干涉图:
{δφ1,δφ2,...,δφN};
确定相干散射点的干涉相位,表示如下:
δφk=φk-φk-1,k=1,2,...,N;
确定相位信号,表示如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910421195.XA CN110018476B (zh) | 2019-05-20 | 2019-05-20 | 时间差分基线集时序干涉sar处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910421195.XA CN110018476B (zh) | 2019-05-20 | 2019-05-20 | 时间差分基线集时序干涉sar处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110018476A CN110018476A (zh) | 2019-07-16 |
CN110018476B true CN110018476B (zh) | 2023-03-24 |
Family
ID=67194064
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910421195.XA Active CN110018476B (zh) | 2019-05-20 | 2019-05-20 | 时间差分基线集时序干涉sar处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110018476B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110673145B (zh) * | 2019-10-24 | 2021-07-27 | 中国地质大学(北京) | 一种基于间断相干的InSAR地表形变监测方法及*** |
CN111538006A (zh) * | 2020-05-13 | 2020-08-14 | 深圳大学 | 基于动态基线的InSAR数字高程模型构建方法及*** |
CN111965645B (zh) * | 2020-08-10 | 2022-04-05 | 中南大学 | 一种顾及几何约束的多基线植被高度反演方法及装置 |
CN112946647A (zh) * | 2021-02-02 | 2021-06-11 | 河海大学 | 大气误差改正InSAR干涉图堆叠地质灾害普查方法和装置 |
CN117473395A (zh) * | 2023-11-21 | 2024-01-30 | 北京市市政工程研究院 | 一种基于时间序列的形变分析方法和*** |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ITMI991154A1 (it) * | 1999-05-25 | 2000-11-25 | Milano Politecnico | Procedimento per misure radar di spostamento di aere urbane e zone franose |
EP2304466B1 (en) * | 2008-07-04 | 2015-10-28 | Telespazio S.p.A. | Identification and analysis of persistent scatterers in series of sar images |
CN103675790B (zh) * | 2013-12-23 | 2016-01-20 | 中国国土资源航空物探遥感中心 | 一种基于高精度DEM提高InSAR技术监测地表形变精度的方法 |
CN104122553B (zh) * | 2014-07-23 | 2017-01-25 | 中国国土资源航空物探遥感中心 | 一种集成多轨道、长条带CTInSAR的区域性地面沉降监测方法 |
CN104123464B (zh) * | 2014-07-23 | 2017-02-22 | 中国国土资源航空物探遥感中心 | 一种高分辨率InSAR时序分析反演地物高程与地面沉降量的方法 |
CN109031300B (zh) * | 2018-09-03 | 2020-06-05 | 中科卫星应用德清研究院 | 合成孔径雷达监测危岩体变形方法及*** |
-
2019
- 2019-05-20 CN CN201910421195.XA patent/CN110018476B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110018476A (zh) | 2019-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110018476B (zh) | 时间差分基线集时序干涉sar处理方法 | |
CN106772342B (zh) | 一种适用于大梯度地表沉降监测的时序差分雷达干涉方法 | |
CN107102333B (zh) | 一种星载InSAR长短基线融合解缠方法 | |
US8482453B2 (en) | Identification and analysis of persistent scatterers in series of SAR images | |
CN101770027B (zh) | 基于InSAR与GPS数据融合的地表三维形变监测方法 | |
CN110412574A (zh) | 一种时空相干性增强的分布式目标InSAR时序处理方法和装置 | |
CN106249236B (zh) | 一种星载InSAR长短基线图像联合配准方法 | |
CN111856459B (zh) | 一种改进的DEM最大似然约束多基线InSAR相位解缠方法 | |
CN109752715B (zh) | 一种sar数据全散射体探测方法及装置 | |
CN111239736B (zh) | 基于单基线的地表高程校正方法、装置、设备及存储介质 | |
CN109696675B (zh) | 基于迪杰斯特拉算法的InSAR时序影像集合配准方法 | |
CN109509219B (zh) | 基于最小生成树的InSAR时序影像集合的配准方法 | |
Fu et al. | A wavelet decomposition and polynomial fitting-based method for the estimation of time-varying residual motion error in airborne interferometric SAR | |
CN104793196A (zh) | 基于距离徙动改进算法的实时sar成像方法 | |
Choi et al. | Large-scale forest height mapping by combining TanDEM-X and GEDI data | |
Fu et al. | Underlying topography estimation over forest areas using single-baseline InSAR data | |
CN112363165B (zh) | 一种林下地形反演方法、装置、设备及介质 | |
CN106646471B (zh) | 基于方位空变误差补偿的机载高分辨sar成像方法 | |
CN109946682B (zh) | 基于ICESat/GLAS的GF3数据基线估计方法 | |
CN110657742A (zh) | 含水层形变信号分离方法、装置、设备及可读存储介质 | |
CN109143188B (zh) | Tops哨兵-1数据电离层校正方法 | |
Buckley et al. | Evaluating ScanSAR interferometry deformation time series using bursted stripmap data | |
Mao et al. | Ionospheric phase delay correction for time series multiple-aperture InSAR constrained by polynomial deformation model | |
Lombardi et al. | Accuracy of high resolution CSK interferometric Digital Elevation Models | |
CN113341410A (zh) | 一种大范围林下地形估计方法、装置、设备及介质 |
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 |