CN103675790A - 一种基于高精度DEM提高InSAR技术监测地表形变精度的方法 - Google Patents
一种基于高精度DEM提高InSAR技术监测地表形变精度的方法 Download PDFInfo
- Publication number
- CN103675790A CN103675790A CN201310717105.4A CN201310717105A CN103675790A CN 103675790 A CN103675790 A CN 103675790A CN 201310717105 A CN201310717105 A CN 201310717105A CN 103675790 A CN103675790 A CN 103675790A
- Authority
- CN
- China
- Prior art keywords
- error
- phase
- radar
- error phase
- function
- 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
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/93—Radar or analogous systems specially adapted for specific applications for anti-collision purposes
- G01S13/933—Radar or analogous systems specially adapted for specific applications for anti-collision purposes of aircraft or spacecraft
- G01S13/935—Radar or analogous systems specially adapted for specific applications for anti-collision purposes of aircraft or spacecraft for terrain-avoidance
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
Abstract
一种基于高精度DEM提高InSAR技术监测地表形变精度的方法,它有五大步骤:步骤一、由雷达数据生成干涉纹图;步骤二、差分干涉相位图的生成;步骤三、误差相位的构成及特征分析;步骤四、误差相位最优函数校正模型的建立;步骤五、根据步骤二和步骤四的结果恢复监测区的地表形变信息。该方法通过提取研究区不同区域的误差相位和高程值或误差相位、高程以及沿距离/方位向的坐标值,基于最小二乘法,分别建立相应区域误差相位的最优函数校正模型,最后,将模拟的误差相位从差分干涉图中去除,进而恢复监测区沿雷达视线向的形变信息。本发明在星载合成孔径雷达监测地表形变技术应用领域具有实用价值和广阔的应用前景。
Description
技术领域
本发明涉及一种基于高精度数字高程模型(DEM)提高合成孔径雷达干涉测量(InSAR)技术监测地表形变精度的方法,属于星载合成孔径雷达监测地表形变技术领域,它适用于基于高分辨率雷达数据和外部高精度DEM,根据误差相位与研究区高程之间的某种函数关系,建立研究区的区域最优误差相位校正模型,恢复研究区地面目标沿雷达视线向的形变信息,从而进一步提高InSAR技术监测地表形变的能力。
背景技术
InSAR技术监测地表形变的核心是通过获取雷达与地面目标之间的相位信息恢复地面监测目标的几何及形变特征,对地表微小形变具有很强的敏感性。由于雷达两次观测条件的变化,包括雷达轨道偏差、地物后向散射特性变化以及噪声的影响,使得两次获取的地面目标的雷达信号相干性降低,难以完成空间离散分布的雷达目标的相位解缠。此外,雷达卫星的基线误差、引入的外部DEM数据不准确以及雷达卫星两次观测时刻大气波动的影响等,这些因素都会降低InSAR技术监测地表形变结果的精度。针对这一问题,当前常用的解决方法有两种:一种是传统合成孔径雷达差分干涉技术(D-InSAR);一种是基于具有稳定散射特性的相干点目标时间序列分析技术。
传统D-InSAR技术的基本思路是,在保证主辅影像配准误差在1/8个像元以下,并对主辅影像进行带通滤波处理后,通过主辅影像共轭相乘及差分处理去除平地相位和地形相位,最终得到的差分干涉图中的误差相位主要是由雷达数据两次获取期间较大的大气波动引起的。在可获取外部大气数据的情况下,将大气层中的水汽含量转换为大气相位,将其从差分相位图中消除,进而得到监测区地表的形变信息。这种方法假设轨道基线及外部DEM数据完全精确,数据在处理过程中产生的误差忽略不计,单纯依靠外部地面气象参数进行水汽延迟校正,但现有的大气水汽监测网络以及物理模型并不能对大气中的水汽含量沿微波信号的路径进行足够准确的外推;基于具有稳定散射特性的相干点目标时间序列分析技术的基本思路是,基于具有稳定散射特性的相干点目标,根据其干涉相位的构成及各相位分量的时空特征,对相干目标的干涉相位进行时间序列分析,估算DEM误差、大气波动与轨道误差以及噪声等,逐步将这些非形变相位分离,从而准确获取地表的形变信息。这些方法包括永久散射体干涉测量技术(PSI)、短基线干涉测量法(SBAS)、点目标干涉测量分析(IPTA)等,它们对雷达数据量的要求较高,一般需要25景甚至更多的SAR影像,当数据量较少且监测区非线性形变梯度较大时,应用中存在一定的局限。
发明内容
本发明的目的是,针对现有方法没有考虑轨道基线及外部DEM数据不准确引起的误差,使得InSAR技术监测地表形变结果的精度不高,甚至有时难以恢复短时间段内地表形变信息的问题,提供了一种基于高精度DEM提高InSAR技术监测地表形变精度的方法,它是一种通过建立差分干涉图中误差相位与研究区高程之间的最优函数校正模型去除误差相位提高InSAR监测地表形变精度的方法。通过推导、分析发现,差分干涉图中的误差相位主要由两部分组成:一是轨道基线和外部DEM数据不准确引起的误差相位;一是大气波动,即大气湿分量变化引起的延迟相位,两者都与研究区的高程密切相关。另外,在地形起伏较大、降水丰沛的高山峡谷区的局部区域,大气延迟相位也是地理位置的函数。于是根据研究区误差相位的分布特征不同,分别提取相应区域内的误差相位、高程信息以及沿距离向/方位向的坐标,基于最小二乘法,对误差相位与高程、沿距离向/方位向的坐标值进行函数拟合,在保证函数模型简单、残差平方和最小的前提下,建立相应区域误差相位的最优多项式校正模型。最后,将模型模拟的误差相位从差分干涉图中去除,进而较为准确地恢复出形变区域沿雷达视线向的形变信息。
技术解决方案
本发明的技术解决方案见附图1,包括以下五个步骤:
(1)由雷达数据生成干涉纹图;
(2)差分干涉相位图的生成;
(3)误差相位的构成及特征分析;
(4)误差相位最优函数校正模型的建立;
(5)根据步骤二和步骤四的结果恢复监测区的地表形变信息。
本发明一种基于高精度DEM提高InSAR技术监测地表形变精度的方法,该方法具体步骤如下:
步骤一:由雷达数据生成干涉纹图
通过星载SAR传感器获取的雷达数据,基于相干性进行主辅影像选择。将主辅影像进行精确配准,根据配准多项式将辅影像与主影像的配准误差控制在1/8个像元以下,并对两者进行带通滤波处理,然后将它们逐像素共轭相乘生成干涉纹图。
步骤二:差分干涉相位图的生成
基于卫星轨道参数及干涉纹图的条纹频率,分别估算平行基线和垂直基线;根据SAR影像的坐标系统及投影等信息,将外部高精度DEM转换到SAR坐标系下,结合卫星轨道信息模拟研究区的地形相位、平地相位,将其从干涉纹图中去除,得到差分干涉相位图。
步骤三:误差相位的构成及特征分析
经过差分处理并对差分干涉相位图进行相位解缠,解缠后的差分干涉相位图中任一像素x的差分干涉相位△φdiff(x)可由形变相位△φdef(x)和误差相位△φerror(x)两部分组成,其中误差相位△φerror(x)可由如下公式表示。
△φerror(x)=△φε-h(x)+△φb(x)+△φatm(x)+△w(x) (1)
其中,△φε-h(x)为垂直基线及外部DEM数据不准确造成的地形误差相位,△φb(x)为平行基线不准确引起的平地误差相位,△φatm(x)为大气湿延迟引起的误差相位,△w(x)为随机的噪声误差相位。
垂直基线及外部DEM数据不准确造成的地形误差相位的计算公式如下所示。
其中,h为雷达影像覆盖区任一点x处的高程值,△h为由于外部的DEM数据不准确,雷达影像覆盖区相应点x处的高程误差值,B⊥为两次雷达数据获取时刻雷达卫星之间的垂直基线长度,△B⊥为垂直基线误差,R为雷达影像覆盖区域中任一点x与第一副天线之间的距离,θ0为雷达卫星的入射角,λ为雷达波的波长。
平行基线不准确引起的平地误差相位的计算公式如下所示。
其中,B//为雷达数据两次获取时刻雷达卫星之间的平行基线长度,△B//为平行基线误差,H为雷达卫星至地心的距离,R为雷达影像覆盖区域中任一点x与第一副天线之间的距离,R'为雷达影像覆盖范围内地面任一点至地心的距离,α为基线与水平面的夹角,θ0为雷达卫星的入射角,λ为雷达波的波长。
由于雷达卫星在两次获取雷达影像时,雷达天线之间的距离B近似等于雷达影像覆盖区域中任一点x到第一副天线之间的距离R与B和R之间夹角余弦的乘积,其计算公式可表达如下所示
于是,△φb(x)≈0
大气湿延迟引起的误差相位△φatm(x),在空间域呈低频变化,但在空气湿度变化较大的高山峡谷区(如水库库区),大气湿分量在空间域的分布特征不尽相同,不仅与研究区高程h有关,在某些局部区域还与研究区的位置x密切相关。
雷达数据处理过程中产生的随机噪声误差相位△w(x),通过滤波处理可基本消除。
可见,差分干涉相位图中的误差相位主要由轨道基线和外部DEM数据不准确引起的误差相位和大气波动引起的误差相位两部分构成,因此,误差相位可进一步用下列公式表达。
步骤四:误差相位最优函数校正模型的建立
差分干涉相位图中的误差相位与高程密切相关,通过提取相应区域(除形变区)的误差相位△φerror(x)和高程值h,基于最小二乘法,尝试采用如下多种类型函数,分别进行曲线拟合。
多项式函数:△φerror(x)=a0+a1×h+a2×h2+a3×h3+… (6)
指数函数:△φerror(x)=a×eb×h (7)
高斯函数:
幂函数:△φerror(x)=a×hb+c (9)
其中,a0,a1,a2,a3,…,a,b,c是各拟合函数的系数。
通过多次试验比较发现,二次多项式函数刻画的误差相位校正模型,在研究区高程范围内没有明显的转折点,残差相位相对较小,且随着多项式最高次数继续增加,其最高次数的系数小于10-9,这使得拟合多项式的误差平方和减小的程度不明显,甚至误差平方和基本不变。其他类型的函数模型,指数函数在研究区高程范围内有明显的转折点,不符合误差相位随高程递减变化的规律;高斯函数模型不但有明显转折点,而且误差相位的残留误差平方和较大;幂函数模型描述的误差相位随高程的变化趋势以及残差相位的最小平方和与多项式函数模型基本一致,但其表达方式不如多项式简单。因此,误差相位的最优函数校正模型可以用关于高程h的二次多项式函数较为准确地表达。
△φerror(x)=a0+a1×h+a2×h2 (10)
其中,a0,a1,a2是二次多项式函数的系数。
在研究区的某些局部区域,虽然大气湿分量在空间域呈低频变化,但在地形起伏较大的高山峡谷区,空气湿度变化较大,误差相位的变化特征不仅与高程有关,还与研究区的位置密切相关,此局部区域误差相位的最优函数校正模型应用如下关于高程h的二次和关于位置x的一次多项式函数进行表达。
△φerror(x)=a0+a1×h+a2×h2+b1×x+b2×x×h (11)
其中,a0,a1,a2,b1,b2是二次多项式函数的系数。
步骤五:根据步骤二和步骤四的结果恢复监测区的地表形变信息
根据建立的不同区域的最优函数校正模型,模拟出研究区整个区域的误差相位分布图,将其从差分干涉图中去除。最后,由卫星轨道参数将去除误差相位后的相位信息转换为沿雷达视线向的形变信息,进而恢复出地表形变区的范围及形变量。
优点及功效:本发明一种基于高精度DEM提高InSAR技术监测地表形变精度的方法,其优点是:
(1)本方法基于高精度DEM,利用传统合成孔径雷达差分干涉(D-InSAR)技术,根据误差相位与高程之间的最优函数模型消除差分干涉图中的误差相位,从而提高干涉测量的精度。
(2)本方法在无需订购多景雷达数据进行回归分析,只需两景雷达数据,且无需获取研究区外部大气数据的情况下,通过建立研究区误差相位关于高程的最优函数校正模型即可获取研究区误差相位的分布特征。
(3)本方法利用误差相位与高程之间的函数模型消除差分干涉图中的误差相位,这在地形起伏较大、大气湿延迟现象严重的高山地区极为有效,可以准确地获取地表短时间内的非线性形变,从而极大地扩大了传统差分干涉雷达的测量范围。
附图说明
图1为本发明的流程图。
图2为有严重误差相位存在的解缠后的差分干涉相位图。在树坪滑坡体区域,误差相位将滑坡体的形变信息完全掩盖。其中,等值线上数值代表差分干涉图中不同位置相位值的大小,单位为弧度;粗虚线范围内为树坪滑坡体,细实线和细虚线囊括的范围分别代表区域I和区域II。整个研究区中的误差相位与高程密切相关,且随着高程的增加误差相位值不断减小,但树坪滑坡体区域和区域II中的误差相位不仅与高程密切相关,还是研究区沿雷达距离向坐标的函数,区域A为区域I、II的共同区域。Az代表雷达方位向,Rg代表雷达距离向。
图3为由最优函数校正模型模拟的研究区的误差相位分布图,该误差相位的分布特征与差分干涉图(即图2)中系统误差相位的分布规律基本一致。其中,等值线上数值代表研究区不同位置误差相位值的大小,单位为弧度。Az代表雷达方位向,Rg代表雷达距离向。
图4为树坪滑坡体在11天中沿雷达视线向的形变信息图。其中,粗虚线范围内为树坪滑坡体区域,等值线上数值代表获得的树坪滑坡体不同区域沿雷达视线向上的形变值,两个形变中心明显的表示出来,单位为毫米;树坪滑坡体区域外为稳定区域,等值线上数值代表残留的误差值,单位为毫米,此时整个研究区残留的误差不但其数值较小,且只在局部区域集中存在。Az代表雷达方位向,Rg代表雷达距离向。
具体实施方式
见附图,为了更好的说明本发明的方法与步骤,利用三峡水库南岸树坪滑坡区两景分辨率为1m的间隔时间为11天的TerraSAR-X雷达数据与精度皆为1m的Lidar数字高程模型(DEM)数据为例,进行误差相位去除,恢复树坪滑坡体滑动区形变信息的试验。
(1)试验所用的设备为图形工作站,规格型号为Dell Precision T7400。
(2)见图1,本发明一种基于高精度DEM提高InSAR技术监测地表形变精度的方法,该方法具体步骤如下:
步骤一:由雷达数据生成干涉纹图
基于传统差分干涉测量方法(D-InSAR),以2012年1月2日获取的分辨率为1m的TerraSAR-X雷达数据为主影像,将2012年1月13日获取的同样分辨率的TerraSAR-X雷达数据与之进行精确配准,在保证主辅影像的配准误差控制在1/8个像元以下,并对配准后的这两幅影像进行带通滤波处理,最后将两者逐像素共轭相乘生成干涉纹图。
步骤二:差分干涉相位图的生成
基于卫星轨道参数及干涉纹图的条纹频率,分别进行空间基线的估算,由空间基线计算干涉图的平地相位;根据SAR影像的坐标系统及投影等信息,对分辨率为1m的数字高程模型(DEM)数据进行坐标转换,即地理编码,通过进行SAR图像和外部DEM之间坐标系统的转换,使外部DEM处于SAR坐标系下,结合卫星轨道信息模拟出研究区正确的地形相位,最后将平地相位、地形相位从干涉纹图中去除,得到差分干涉相位图。
步骤三:误差相位的构成及特征分析
为准确恢复树坪滑坡体的形变信息,首先以某一稳定点为相位解缠基准点,采用最小费用流算法对差分干涉图进行相位解缠,解缠后的差分干涉相位图见图2。
从图2中可以看出,差分干涉图中存在严重的误差相位,树坪滑坡体滑动区域的形变相位完全被误差相位“淹没”,无法将滑坡体的形变信息和周围稳定区域的误差相位信息区分开来。由于误差相位主要由轨道基线和外部DEM数据不准确引起的误差相位和大气波动引起的误差相位两部分构成,差分干涉相位图中误差相位是研究区高程的函数,且随着高程的增加,误差相位值不断降低。但II区域及树坪滑坡区,它们的误差相位分布特征基本一致,不仅与高程密切相关,还是研究区沿雷达距离向坐标的函数。
于是,为恢复整幅差分干涉相位图中误差相位的分布状态,针对I、II两区域的误差相位分布特征不同,分别建立相应区域误差相位的最优函数校正模型。
步骤四:误差相位最优函数校正模型的建立
在I、II区域中,分别提取高相干点上的误差相位 和相应的高程值hI、hII,并提取II区域中相应点处沿距离向的坐标值x。基于最小二乘法,采用式(10)和式(11)分别建立I、II区域的误差相位校正模型如下所示。
它们的残差相位都介于[-2.0,1.0]之间。由于树坪滑坡体区域的误差相位分布特征与II区域基本一致,则树坪滑坡区的误差相位分布特征也遵循式(13)的表达规律。
由提取的I、II区域及树坪滑坡体区域高相干点上的误差相位 和相应的高程值hI、hII以及II区域和树坪滑坡体区域中相应点处沿距离向的坐标值x,根据建立的I、II区域误差相位的最优函数校正模型,即式(12)、(13),计算整个研究区高相干点上的误差相位值,其中,I、II区域的共同区域A中高相干点上的误差相位值为由式(12)、(13)分别计算得到的误差相位的平均值。最后,通过最邻近插值法,将计算得到的整个区域高相干点上的误差相位值进行插值,进而得到研究区误差相位的分布特征图,见图3。在树坪滑坡体区域,误差相位的分布特征也形象地表达了出来。
步骤五:根据步骤二和步骤四的结果恢复监测区的地表形变信息
将由最优函数校正模型式(12)、(13)模拟出的研究区的误差相位,从解缠后的差分干涉相位图中消除,此时,在树坪滑坡体区域残留的相位只是由滑坡体的滑动形变引起的。最后基于卫星轨道参数将去除误差相位后的相位信息转换为沿雷达视线向的形变信息,树坪滑坡体在2012年1月2日至2012年1月13日间的滑动区域及滑动量就准确、鲜明地表达了出来,见图4。
Claims (1)
1.一种基于高精度DEM提高InSAR技术监测地表形变精度的方法,其特征在于:该方法具体步骤如下:
步骤一:由雷达数据生成干涉纹图
通过星载SAR传感器获取的雷达数据,基于相干性进行主辅影像选择;将主辅影像进行精确配准,根据配准多项式将辅影像与主影像的配准误差控制在1/8个像元以下,并对两者进行带通滤波处理,然后将它们逐像素共轭相乘生成干涉纹图;
步骤二:差分干涉相位图的生成
基于卫星轨道参数及干涉纹图的条纹频率,分别估算平行基线和垂直基线;根据SAR影像的坐标系统及投影信息,将外部高精度DEM转换到SAR坐标系下,结合卫星轨道信息模拟研究区的地形相位、平地相位,将其从干涉纹图中去除,得到差分干涉相位图;
步骤三:误差相位的构成及特征分析
经过差分处理并对差分干涉相位图进行相位解缠,解缠后的差分干涉相位图中任一像素x的差分干涉相位△φdiff(x)由形变相位△φdef(x)和误差相位△φerror(x)两部分组成,其中误差相位△φerror(x)由如下公式表示:
△φerror(x)=△φε-h(x)+△φb(x)+△φatm(x)+△w(x) (1)
其中,△φε-h(x)为垂直基线及外部DEM数据不准确造成的地形误差相位,△φb(x)为平行基线不准确引起的平地误差相位,△φatm(x)为大气湿延迟引起的误差相位,△w(x)为随机的噪声误差相位;
垂直基线及外部DEM数据不准确造成的地形误差相位的计算公式如下所示:
其中,h为雷达影像覆盖区任一点x处的高程值,△h为由于外部的DEM数据不准确,雷达影像覆盖区相应点x处的高程误差值,B⊥为两次雷达数据获取时刻雷达卫星之间的垂直基线长度,△B⊥为垂直基线误差,R为雷达影像覆盖区域中任一点x与第一副天线之间的距离,θ0为雷达卫星的入射角,λ为雷达波的波长;
平行基线不准确引起的平地误差相位的计算公式如下所示:
其中,B//为雷达数据两次获取时刻雷达卫星之间的平行基线长度,△B//为平行基线误差,H为雷达卫星至地心的距离,R为雷达影像覆盖区域中任一点x与第一副天线之间的距离,R'为雷达影像覆盖范围内地面任一点至地心的距离,α为基线与水平面的夹角,θ0为雷达卫星的入射角,λ为雷达波的波长;
由于雷达卫星在两次获取雷达影像时,雷达天线之间的距离B近似等于雷达影像覆盖区域中任一点x到第一副天线之间的距离R与B和R之间夹角余弦的乘积,其计算公式如下所示:
于是,△φb(x)≈0
大气湿延迟引起的误差相位△φatm(x),在空间域呈低频变化,但在空气湿度变化较大的高山峡谷区如水库库区,大气湿分量在空间域的分布特征不尽相同,不仅与研究区高程h有关,在某些局部区域还与研究区的位置x密切相关;
雷达数据处理过程中产生的随机噪声误差相位△w(x),通过滤波处理基本消除;
可见,差分干涉相位图中的误差相位主要由轨道基线和外部DEM数据不准确引起的误差相位和大气波动引起的误差相位两部分构成,因此,误差相位可进一步用下列公式表达:
步骤四:误差相位最优函数校正模型的建立
差分干涉相位图中的误差相位与高程密切相关,通过提取相应区域的误差相位△φerror(x)和高程值h,基于最小二乘法,尝试采用如下多种类型函数,分别进行曲线拟合;
多项式函数:△φerror(x)=a0+a1×h+a2×h2+a3×h3+… (6)
指数函数:△φerror(x)=a×eb×h (7)
高斯函数:
幂函数:△φerror(x)=a×hb+c (9)
其中,a0,a1,a2,a3,…,a,b,c是各拟合函数的系数;
通过多次试验比较发现,二次多项式函数刻画的误差相位校正模型,在研究区高程范围内没有明显的转折点,残差相位相对较小,且随着多项式最高次数继续增加,其最高次数的系数小于10-9,这使得拟合多项式的误差平方和减小的程度不明显,甚至误差平方和基本不变;其他类型的函数模型,指数函数在研究区高程范围内有明显的转折点,不符合误差相位随高程递减变化的规律;高斯函数模型不但有明显转折点,而且误差相位的残留误差平方和较大;幂函数模型描述的误差相位随高程的变化趋势以及残差相位的最小平方和与多项式函数模型基本一致,但其表达方式不如多项式简,因此,误差相位的最优函数校正模型用关于高程h的二次多项式函数较为准确地表达;
△φerror(x)=a0+a1×h+a2×h2 (10)
其中,a0,a1,a2是二次多项式函数的系数;
在研究区的某些局部区域,虽然大气湿分量在空间域呈低频变化,但在地形起伏较大的高山峡谷区,空气湿度变化较大,误差相位的变化特征不仅与高程有关,还与研究区的位置密切相关,此局部区域误差相位的最优函数校正模型应用如下关于高程h的二次和关于位置x的一次多项式函数进行表达;
△φerror(x)=a0+a1×h+a2×h2+b1×x+b2×x×h (11)
其中,a0,a1,a2,b1,b2是二次多项式函数的系数;
步骤五:根据步骤二和步骤四的结果恢复监测区的地表形变信息
根据建立的不同区域的最优函数校正模型,模拟出研究区整个区域的误差相位分布图,将其从差分干涉图中去除,最后,由卫星轨道参数将去除误差相位后的相位信息转换为沿雷达视线向的形变信息,进而恢复出地表形变区的范围及形变量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310717105.4A CN103675790B (zh) | 2013-12-23 | 2013-12-23 | 一种基于高精度DEM提高InSAR技术监测地表形变精度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310717105.4A CN103675790B (zh) | 2013-12-23 | 2013-12-23 | 一种基于高精度DEM提高InSAR技术监测地表形变精度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103675790A true CN103675790A (zh) | 2014-03-26 |
CN103675790B CN103675790B (zh) | 2016-01-20 |
Family
ID=50313930
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310717105.4A Expired - Fee Related CN103675790B (zh) | 2013-12-23 | 2013-12-23 | 一种基于高精度DEM提高InSAR技术监测地表形变精度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103675790B (zh) |
Cited By (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104062660A (zh) * | 2014-07-14 | 2014-09-24 | 中南大学 | 一种基于时域离散InSAR干涉对的矿区地表时序形变监测方法 |
CN104122553A (zh) * | 2014-07-23 | 2014-10-29 | 中国国土资源航空物探遥感中心 | 一种集成多轨道、长条带CTInSAR的区域性地面沉降监测方法 |
CN104123464A (zh) * | 2014-07-23 | 2014-10-29 | 中国国土资源航空物探遥感中心 | 一种高分辨率InSAR时序分析反演地物高程与地面沉降量的方法 |
CN104391296A (zh) * | 2014-10-15 | 2015-03-04 | 淮海工学院 | 基于广义最小二乘平差的雷达三维形变场重构技术 |
CN104730521A (zh) * | 2015-04-01 | 2015-06-24 | 北京航空航天大学 | 一种基于非线性优化策略的SBAS-DInSAR方法 |
CN105158760A (zh) * | 2015-08-10 | 2015-12-16 | 中南大学 | 一种利用InSAR反演地下流体体积变化和三维地表形变的方法 |
CN105204079A (zh) * | 2015-08-31 | 2015-12-30 | 中国科学院测量与地球物理研究所 | 一种利用TanDEM-X双站InSAR提取地震滑坡体积的方法 |
CN105678716A (zh) * | 2016-02-25 | 2016-06-15 | 内蒙古工业大学 | 一种地基sar大气干扰相位校正方法及装置 |
CN106556834A (zh) * | 2016-11-24 | 2017-04-05 | 首都师范大学 | 一种从两平行轨道sar数据集中精确提取地面垂直形变方法 |
CN107037428A (zh) * | 2017-03-27 | 2017-08-11 | 中国科学院遥感与数字地球研究所 | 一种提高星载双站差分InSAR提取形变精度的方法 |
CN107202985A (zh) * | 2017-04-13 | 2017-09-26 | 长安大学 | 一种基于干涉图闭合环的InSAR解缠误差探测方法 |
CN107621636A (zh) * | 2017-11-13 | 2018-01-23 | 河海大学 | 一种基于psi的大型铁路桥梁健康监测方法 |
CN108050964A (zh) * | 2018-01-30 | 2018-05-18 | 长沙深之瞳信息科技有限公司 | 一种基于微波干涉的二维面形变监测方法及系统 |
CN108088358A (zh) * | 2017-12-18 | 2018-05-29 | 电子科技大学 | 一种基于多基线雷达轨道形变检测方法 |
CN108280880A (zh) * | 2018-01-24 | 2018-07-13 | 长春工程学院 | 一种利用遥感影像提高山体的数字高程数据分辨率的方法 |
CN108445458A (zh) * | 2018-03-14 | 2018-08-24 | 中煤航测遥感集团有限公司 | 一种合成孔径雷达轨道误差消除方法和装置 |
CN109031300A (zh) * | 2018-09-03 | 2018-12-18 | 中科卫星应用德清研究院 | 合成孔径雷达监测危岩体变形方法及系统 |
CN109035312A (zh) * | 2018-07-17 | 2018-12-18 | 中国人民解放军国防科技大学 | 一种dem辅助的sar图像高精度配准方法 |
CN109061641A (zh) * | 2018-07-06 | 2018-12-21 | 中南大学 | 一种基于序贯平差的InSAR时序地表形变监测方法 |
CN109541593A (zh) * | 2018-10-30 | 2019-03-29 | 北京航空航天大学 | 一种改进的最小费用流InSAR相位解缠方法 |
CN110018476A (zh) * | 2019-05-20 | 2019-07-16 | 月明星(北京)科技有限公司 | 时间差分基线集时序干涉sar处理方法 |
CN110456345A (zh) * | 2019-06-28 | 2019-11-15 | 深圳市水务规划设计院股份有限公司 | 一种基于InSAR技术的建筑物倾斜监测方法 |
CN110780297A (zh) * | 2019-09-30 | 2020-02-11 | 中国国土勘测规划院 | 一种采用星载合成孔径雷达监测土壤侵蚀的方法与系统 |
CN110806170A (zh) * | 2019-11-26 | 2020-02-18 | 武汉大学 | 监测高堆石坝变形的方法 |
CN110826518A (zh) * | 2019-11-14 | 2020-02-21 | 中国地质科学院矿产资源研究所 | 一种遥感影像隐伏地质构造信息提取方法 |
CN111059998A (zh) * | 2019-12-31 | 2020-04-24 | 中国地质大学(北京) | 一种基于高分辨率的时序InSAR形变监测方法及系统 |
CN111239736A (zh) * | 2020-03-19 | 2020-06-05 | 中南大学 | 基于单基线的地表高程校正方法、装置、设备及存储介质 |
CN111505626A (zh) * | 2020-04-26 | 2020-08-07 | 中国科学院国家空间科学中心 | 一种利用底视差分干涉测量二维地形坡度的方法 |
CN111580101A (zh) * | 2020-04-14 | 2020-08-25 | 中国测绘科学研究院 | 一种基于外部DEM的InSAR基线误差无控改正方法及装置 |
CN111999733A (zh) * | 2020-08-11 | 2020-11-27 | 中国海洋大学 | 海岸带大坝稳定性监测与淹没脆弱性评估方法、系统 |
CN112596055A (zh) * | 2020-12-08 | 2021-04-02 | 中国地质大学(武汉) | 一种改正InSAR DEM残余系统误差的方法 |
CN112649807A (zh) * | 2020-04-09 | 2021-04-13 | 湖南科技学院 | 一种基于小波多尺度相关性分析的机载InSAR轨道误差去除方法 |
CN112986988A (zh) * | 2021-01-25 | 2021-06-18 | 中国电子科技集团公司第二十九研究所 | 一种特殊地形条件下sar干涉测量形变量解析方法 |
CN113340191A (zh) * | 2021-04-08 | 2021-09-03 | 河北省高速公路延崇筹建处 | 时间序列干涉sar的形变量测量方法及sar系统 |
CN113705025A (zh) * | 2021-10-29 | 2021-11-26 | 中科星睿科技(北京)有限公司 | 信息提取与传输方法、装置、电子设备和计算机可读介质 |
CN114487024A (zh) * | 2021-12-31 | 2022-05-13 | 河南省日立信股份有限公司 | 一种基于幂函数的钯合金氢气传感器的校准拟合方法 |
CN114791273A (zh) * | 2022-06-24 | 2022-07-26 | 中国铁道科学研究院集团有限公司铁道建筑研究所 | 一种针对滑坡的InSAR形变监测结果解释方法 |
CN115201825A (zh) * | 2022-09-16 | 2022-10-18 | 眉山环天智慧科技有限公司 | 一种InSAR震间形变监测中的大气延迟校正方法 |
CN115471144A (zh) * | 2022-11-14 | 2022-12-13 | 四川省公路规划勘察设计研究院有限公司 | 基于多源数据融合的泥石流监测预警方法、装置及介质 |
CN115629384A (zh) * | 2022-12-08 | 2023-01-20 | 中南大学 | 一种时序InSAR误差的改正方法及相关设备 |
CN115856889A (zh) * | 2023-02-03 | 2023-03-28 | 北京数慧时空信息技术有限公司 | 一种误差自动校正的InSAR时序变形监测方法 |
CN116736306A (zh) * | 2023-08-15 | 2023-09-12 | 成都理工大学 | 一种基于高分三号的时序雷达干涉监测方法 |
CN117214898A (zh) * | 2023-11-09 | 2023-12-12 | 中国电子科技集团公司第十四研究所 | 一种基于多维电磁信息深度融合的宽尺度地表形变精细化遥感方法 |
CN117274342A (zh) * | 2023-11-21 | 2023-12-22 | 中铁水利水电规划设计集团有限公司 | 一种基于卫星数据的水利工程形变监测方法 |
CN117310706A (zh) * | 2023-11-28 | 2023-12-29 | 中山大学 | 一种地基雷达间断形变监测方法及系统 |
CN117665809A (zh) * | 2023-12-21 | 2024-03-08 | 西南林业大学 | 反演森林冠层高度方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100425283B1 (ko) * | 2003-04-18 | 2004-03-31 | (주)충청측량설계공사 | 간섭측정용 합성 개구 레이다의 기하학적 특성을 이용한지형고도 측정방법 |
KR100441590B1 (ko) * | 2003-04-18 | 2004-07-23 | (주)충청측량설계공사 | 간섭측정용 합성 개구 레이다의 기하학적 특성을 이용하여지형고도를 측정하기 위한 디지털 고도모형 생성방법 |
CN102608584A (zh) * | 2012-03-19 | 2012-07-25 | 中国测绘科学研究院 | 基于多项式反演模型的时间序列InSAR形变监测方法及装置 |
CN102645650A (zh) * | 2012-03-06 | 2012-08-22 | 北京北科安地科技发展有限公司 | 一种基于D-InSAR干涉差分的滑坡动态识别及监测技术 |
CN103091676A (zh) * | 2013-01-22 | 2013-05-08 | 中国矿业大学 | 矿区地表开采沉陷合成孔径雷达干涉测量的监测及解算方法 |
-
2013
- 2013-12-23 CN CN201310717105.4A patent/CN103675790B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100425283B1 (ko) * | 2003-04-18 | 2004-03-31 | (주)충청측량설계공사 | 간섭측정용 합성 개구 레이다의 기하학적 특성을 이용한지형고도 측정방법 |
KR100441590B1 (ko) * | 2003-04-18 | 2004-07-23 | (주)충청측량설계공사 | 간섭측정용 합성 개구 레이다의 기하학적 특성을 이용하여지형고도를 측정하기 위한 디지털 고도모형 생성방법 |
CN102645650A (zh) * | 2012-03-06 | 2012-08-22 | 北京北科安地科技发展有限公司 | 一种基于D-InSAR干涉差分的滑坡动态识别及监测技术 |
CN102608584A (zh) * | 2012-03-19 | 2012-07-25 | 中国测绘科学研究院 | 基于多项式反演模型的时间序列InSAR形变监测方法及装置 |
CN103091676A (zh) * | 2013-01-22 | 2013-05-08 | 中国矿业大学 | 矿区地表开采沉陷合成孔径雷达干涉测量的监测及解算方法 |
Non-Patent Citations (2)
Title |
---|
MAN LI等: "EXPERIMENTAL STUDY OF ATMOSPHERIC CORRECTION TO INTERFEROGRAM WITH HIGH-RESOLUTION RADAR IMAGES AND DEM OVER SHUPING LANDSLIDE", 《IGARSS 2013》, 26 July 2013 (2013-07-26) * |
李曼等: "基于精细DEM的InSAR大气相位改正试验研究", 《国土资源遥感》, 30 June 2013 (2013-06-30) * |
Cited By (71)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104062660A (zh) * | 2014-07-14 | 2014-09-24 | 中南大学 | 一种基于时域离散InSAR干涉对的矿区地表时序形变监测方法 |
CN104123464B (zh) * | 2014-07-23 | 2017-02-22 | 中国国土资源航空物探遥感中心 | 一种高分辨率InSAR时序分析反演地物高程与地面沉降量的方法 |
CN104122553A (zh) * | 2014-07-23 | 2014-10-29 | 中国国土资源航空物探遥感中心 | 一种集成多轨道、长条带CTInSAR的区域性地面沉降监测方法 |
CN104123464A (zh) * | 2014-07-23 | 2014-10-29 | 中国国土资源航空物探遥感中心 | 一种高分辨率InSAR时序分析反演地物高程与地面沉降量的方法 |
CN104391296A (zh) * | 2014-10-15 | 2015-03-04 | 淮海工学院 | 基于广义最小二乘平差的雷达三维形变场重构技术 |
CN104730521A (zh) * | 2015-04-01 | 2015-06-24 | 北京航空航天大学 | 一种基于非线性优化策略的SBAS-DInSAR方法 |
CN104730521B (zh) * | 2015-04-01 | 2017-03-15 | 北京航空航天大学 | 一种基于非线性优化策略的SBAS‑DInSAR方法 |
CN105158760A (zh) * | 2015-08-10 | 2015-12-16 | 中南大学 | 一种利用InSAR反演地下流体体积变化和三维地表形变的方法 |
CN105158760B (zh) * | 2015-08-10 | 2017-04-26 | 中南大学 | 一种利用InSAR反演地下流体体积变化和三维地表形变的方法 |
CN105204079A (zh) * | 2015-08-31 | 2015-12-30 | 中国科学院测量与地球物理研究所 | 一种利用TanDEM-X双站InSAR提取地震滑坡体积的方法 |
CN105678716A (zh) * | 2016-02-25 | 2016-06-15 | 内蒙古工业大学 | 一种地基sar大气干扰相位校正方法及装置 |
CN105678716B (zh) * | 2016-02-25 | 2018-09-18 | 内蒙古工业大学 | 一种地基sar大气干扰相位校正方法及装置 |
CN106556834A (zh) * | 2016-11-24 | 2017-04-05 | 首都师范大学 | 一种从两平行轨道sar数据集中精确提取地面垂直形变方法 |
CN107037428A (zh) * | 2017-03-27 | 2017-08-11 | 中国科学院遥感与数字地球研究所 | 一种提高星载双站差分InSAR提取形变精度的方法 |
CN107037428B (zh) * | 2017-03-27 | 2019-11-12 | 中国科学院遥感与数字地球研究所 | 一种提高星载双站差分InSAR提取形变精度的方法 |
CN107202985A (zh) * | 2017-04-13 | 2017-09-26 | 长安大学 | 一种基于干涉图闭合环的InSAR解缠误差探测方法 |
CN107202985B (zh) * | 2017-04-13 | 2019-10-11 | 长安大学 | 一种基于干涉图闭合环的InSAR解缠误差探测方法 |
CN107621636A (zh) * | 2017-11-13 | 2018-01-23 | 河海大学 | 一种基于psi的大型铁路桥梁健康监测方法 |
CN107621636B (zh) * | 2017-11-13 | 2021-08-06 | 河海大学 | 一种基于psi的大型铁路桥梁健康监测方法 |
CN108088358B (zh) * | 2017-12-18 | 2019-08-20 | 电子科技大学 | 一种基于多基线雷达轨道形变检测方法 |
CN108088358A (zh) * | 2017-12-18 | 2018-05-29 | 电子科技大学 | 一种基于多基线雷达轨道形变检测方法 |
CN108280880A (zh) * | 2018-01-24 | 2018-07-13 | 长春工程学院 | 一种利用遥感影像提高山体的数字高程数据分辨率的方法 |
CN108050964B (zh) * | 2018-01-30 | 2023-04-18 | 长沙深之瞳信息科技有限公司 | 一种基于微波干涉的二维面形变监测方法及系统 |
CN108050964A (zh) * | 2018-01-30 | 2018-05-18 | 长沙深之瞳信息科技有限公司 | 一种基于微波干涉的二维面形变监测方法及系统 |
CN108445458A (zh) * | 2018-03-14 | 2018-08-24 | 中煤航测遥感集团有限公司 | 一种合成孔径雷达轨道误差消除方法和装置 |
CN109061641A (zh) * | 2018-07-06 | 2018-12-21 | 中南大学 | 一种基于序贯平差的InSAR时序地表形变监测方法 |
CN109035312A (zh) * | 2018-07-17 | 2018-12-18 | 中国人民解放军国防科技大学 | 一种dem辅助的sar图像高精度配准方法 |
CN109035312B (zh) * | 2018-07-17 | 2021-11-09 | 中国人民解放军国防科技大学 | 一种dem辅助的sar图像高精度配准方法 |
CN109031300A (zh) * | 2018-09-03 | 2018-12-18 | 中科卫星应用德清研究院 | 合成孔径雷达监测危岩体变形方法及系统 |
CN109541593A (zh) * | 2018-10-30 | 2019-03-29 | 北京航空航天大学 | 一种改进的最小费用流InSAR相位解缠方法 |
CN110018476A (zh) * | 2019-05-20 | 2019-07-16 | 月明星(北京)科技有限公司 | 时间差分基线集时序干涉sar处理方法 |
CN110456345A (zh) * | 2019-06-28 | 2019-11-15 | 深圳市水务规划设计院股份有限公司 | 一种基于InSAR技术的建筑物倾斜监测方法 |
CN110456345B (zh) * | 2019-06-28 | 2020-11-10 | 深圳市水务规划设计院股份有限公司 | 一种基于InSAR技术的建筑物倾斜监测方法 |
CN110780297A (zh) * | 2019-09-30 | 2020-02-11 | 中国国土勘测规划院 | 一种采用星载合成孔径雷达监测土壤侵蚀的方法与系统 |
CN110780297B (zh) * | 2019-09-30 | 2023-05-23 | 中国国土勘测规划院 | 一种采用星载合成孔径雷达监测土壤侵蚀的方法与系统 |
CN110826518A (zh) * | 2019-11-14 | 2020-02-21 | 中国地质科学院矿产资源研究所 | 一种遥感影像隐伏地质构造信息提取方法 |
CN110826518B (zh) * | 2019-11-14 | 2022-06-03 | 中国地质科学院矿产资源研究所 | 一种遥感影像隐伏地质构造信息提取方法 |
CN110806170A (zh) * | 2019-11-26 | 2020-02-18 | 武汉大学 | 监测高堆石坝变形的方法 |
CN111059998A (zh) * | 2019-12-31 | 2020-04-24 | 中国地质大学(北京) | 一种基于高分辨率的时序InSAR形变监测方法及系统 |
CN111239736A (zh) * | 2020-03-19 | 2020-06-05 | 中南大学 | 基于单基线的地表高程校正方法、装置、设备及存储介质 |
CN112649807A (zh) * | 2020-04-09 | 2021-04-13 | 湖南科技学院 | 一种基于小波多尺度相关性分析的机载InSAR轨道误差去除方法 |
CN111580101A (zh) * | 2020-04-14 | 2020-08-25 | 中国测绘科学研究院 | 一种基于外部DEM的InSAR基线误差无控改正方法及装置 |
CN111505626B (zh) * | 2020-04-26 | 2022-07-26 | 中国科学院国家空间科学中心 | 一种利用底视差分干涉测量二维地形坡度的方法 |
CN111505626A (zh) * | 2020-04-26 | 2020-08-07 | 中国科学院国家空间科学中心 | 一种利用底视差分干涉测量二维地形坡度的方法 |
CN111999733A (zh) * | 2020-08-11 | 2020-11-27 | 中国海洋大学 | 海岸带大坝稳定性监测与淹没脆弱性评估方法、系统 |
CN111999733B (zh) * | 2020-08-11 | 2022-06-03 | 中国海洋大学 | 海岸带大坝稳定性监测与淹没脆弱性评估方法、系统 |
CN112596055A (zh) * | 2020-12-08 | 2021-04-02 | 中国地质大学(武汉) | 一种改正InSAR DEM残余系统误差的方法 |
CN112596055B (zh) * | 2020-12-08 | 2023-04-25 | 中国地质大学(武汉) | 一种改正InSAR DEM残余系统误差的方法 |
CN112986988A (zh) * | 2021-01-25 | 2021-06-18 | 中国电子科技集团公司第二十九研究所 | 一种特殊地形条件下sar干涉测量形变量解析方法 |
CN112986988B (zh) * | 2021-01-25 | 2022-09-20 | 中国电子科技集团公司第二十九研究所 | 一种特殊地形条件下sar干涉测量形变量解析方法 |
CN113340191A (zh) * | 2021-04-08 | 2021-09-03 | 河北省高速公路延崇筹建处 | 时间序列干涉sar的形变量测量方法及sar系统 |
CN113705025A (zh) * | 2021-10-29 | 2021-11-26 | 中科星睿科技(北京)有限公司 | 信息提取与传输方法、装置、电子设备和计算机可读介质 |
CN114487024A (zh) * | 2021-12-31 | 2022-05-13 | 河南省日立信股份有限公司 | 一种基于幂函数的钯合金氢气传感器的校准拟合方法 |
CN114487024B (zh) * | 2021-12-31 | 2023-11-03 | 河南省日立信股份有限公司 | 一种基于幂函数的钯合金氢气传感器的校准拟合方法 |
CN114791273A (zh) * | 2022-06-24 | 2022-07-26 | 中国铁道科学研究院集团有限公司铁道建筑研究所 | 一种针对滑坡的InSAR形变监测结果解释方法 |
CN114791273B (zh) * | 2022-06-24 | 2022-09-13 | 中国铁道科学研究院集团有限公司铁道建筑研究所 | 一种针对滑坡的InSAR形变监测结果解释方法 |
CN115201825B (zh) * | 2022-09-16 | 2023-01-17 | 眉山环天智慧科技有限公司 | 一种InSAR震间形变监测中的大气延迟校正方法 |
CN115201825A (zh) * | 2022-09-16 | 2022-10-18 | 眉山环天智慧科技有限公司 | 一种InSAR震间形变监测中的大气延迟校正方法 |
CN115471144A (zh) * | 2022-11-14 | 2022-12-13 | 四川省公路规划勘察设计研究院有限公司 | 基于多源数据融合的泥石流监测预警方法、装置及介质 |
CN115471144B (zh) * | 2022-11-14 | 2023-04-07 | 四川省公路规划勘察设计研究院有限公司 | 基于多源数据融合的泥石流监测预警方法、装置及介质 |
CN115629384A (zh) * | 2022-12-08 | 2023-01-20 | 中南大学 | 一种时序InSAR误差的改正方法及相关设备 |
CN115856889A (zh) * | 2023-02-03 | 2023-03-28 | 北京数慧时空信息技术有限公司 | 一种误差自动校正的InSAR时序变形监测方法 |
CN116736306A (zh) * | 2023-08-15 | 2023-09-12 | 成都理工大学 | 一种基于高分三号的时序雷达干涉监测方法 |
CN116736306B (zh) * | 2023-08-15 | 2023-10-24 | 成都理工大学 | 一种基于高分三号的时序雷达干涉监测方法 |
CN117214898A (zh) * | 2023-11-09 | 2023-12-12 | 中国电子科技集团公司第十四研究所 | 一种基于多维电磁信息深度融合的宽尺度地表形变精细化遥感方法 |
CN117214898B (zh) * | 2023-11-09 | 2024-01-23 | 中国电子科技集团公司第十四研究所 | 一种基于多维电磁信息深度融合的宽尺度地表形变精细化遥感方法 |
CN117274342A (zh) * | 2023-11-21 | 2023-12-22 | 中铁水利水电规划设计集团有限公司 | 一种基于卫星数据的水利工程形变监测方法 |
CN117274342B (zh) * | 2023-11-21 | 2024-02-13 | 中铁水利水电规划设计集团有限公司 | 一种基于卫星数据的水利工程形变监测方法 |
CN117310706A (zh) * | 2023-11-28 | 2023-12-29 | 中山大学 | 一种地基雷达间断形变监测方法及系统 |
CN117310706B (zh) * | 2023-11-28 | 2024-02-02 | 中山大学 | 一种地基雷达间断形变监测方法及系统 |
CN117665809A (zh) * | 2023-12-21 | 2024-03-08 | 西南林业大学 | 反演森林冠层高度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103675790B (zh) | 2016-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103675790B (zh) | 一种基于高精度DEM提高InSAR技术监测地表形变精度的方法 | |
Nuth et al. | Co-registration and bias corrections of satellite elevation data sets for quantifying glacier thickness change | |
CN106772342B (zh) | 一种适用于大梯度地表沉降监测的时序差分雷达干涉方法 | |
Mukherjee et al. | Evaluation of vertical accuracy of open source Digital Elevation Model (DEM) | |
JP6106285B2 (ja) | 衛星レーダー干渉度のイオン歪み補正方法及びその装置 | |
Catalão et al. | Merging GPS and atmospherically corrected InSAR data to map 3-D terrain displacement velocity | |
CN107102333B (zh) | 一种星载InSAR长短基线融合解缠方法 | |
CN104111457B (zh) | 一种升降轨PSInSAR地面沉降监测结果的互检验与时序融合方法 | |
CN111398959B (zh) | 基于地表应力应变模型的InSAR时序地表形变监测方法 | |
CN105204079B (zh) | 一种利用TanDEM-X双站InSAR提取地震滑坡体积的方法 | |
Zhang et al. | Satellite SAR geocoding with refined RPC model | |
Liu et al. | Glacier elevation changes (2012–2016) of the Puruogangri Ice Field on the Tibetan Plateau derived from bi-temporal TanDEM-X InSAR data | |
Liao et al. | Ionospheric correction of InSAR data for accurate ice velocity measurement at polar regions | |
CN108663678B (zh) | 基于混合整数优化模型的多基线InSAR相位解缠算法 | |
Tang et al. | Atmospheric correction in time-series SAR interferometry for land surface deformation mapping–A case study of Taiyuan, China | |
Zhao et al. | A multiple-regression model considering deformation information for atmospheric phase screen compensation in ground-based SAR | |
Mao et al. | Estimation and compensation of ionospheric phase delay for multi-aperture InSAR: An azimuth split-spectrum interferometry approach | |
CN111650570B (zh) | 一种地基干涉雷达三维大气校正方法及系统 | |
Zhou et al. | Comparison and validation of different DEM data derived from InSAR | |
Frolich et al. | Synthetic aperture radar interferometry over Rutford Ice Stream and Carlson Inlet, Antarctica | |
Li et al. | Co-registration and residual correction of digital elevation models: A comparative study | |
Chang et al. | InSAR atmospheric distortions mitigation: GPS observations and NCEP FNL data | |
CN113341410B (zh) | 一种大范围林下地形估计方法、装置、设备及介质 | |
Haque et al. | Time series analysis of subsidence in Dhaka City, Bangladesh using Insar | |
CN111580101A (zh) | 一种基于外部DEM的InSAR基线误差无控改正方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160120 Termination date: 20171223 |
|
CF01 | Termination of patent right due to non-payment of annual fee |