CN104360332A - 基于地基合成孔径雷达干涉的大气相位屏提取方法 - Google Patents
基于地基合成孔径雷达干涉的大气相位屏提取方法 Download PDFInfo
- Publication number
- CN104360332A CN104360332A CN201410633483.9A CN201410633483A CN104360332A CN 104360332 A CN104360332 A CN 104360332A CN 201410633483 A CN201410633483 A CN 201410633483A CN 104360332 A CN104360332 A CN 104360332A
- Authority
- CN
- China
- Prior art keywords
- phase
- image
- radar
- ground
- formula
- 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
- 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/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
-
- 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/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/354—Extracting wanted echo-signals
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于地基合成孔径雷达干涉的大气相位屏提取方法,依据地基雷达差分干涉相位模型,基于最小二乘法拟合线性形变和时空滤波方法提取大气相位,该干涉相位模型充分考虑地基雷达成像机理以及差分干涉原理,从而提取大气相位屏,本发明提出振幅离差指数、相干指数以及相位稳定性三级PS点选取方法,最后经Kriging插值获得每景雷达数据的大气相位屏,准确地提取受大气影响产生的相位对雷达干涉测量精度的提高具有重要作用。通过大坝实测数据验证本方法的正确性,具有一定的实际应用价值和通用性。
Description
技术领域
本发明涉及测量学领域技术领域,具体涉及一种基于地基合成孔径雷达干涉的大气相位屏提取方法。
背景技术
地基合成孔径雷达干涉(Ground Based雷达,GB雷达,简称地基雷达)技术是基于微波探测主动成像方式获取监测区域二维影像,通过合成孔径技术和步进频率技术实现雷达影像方位向和距离向的高空间分辨率,克服了星载SAR影像受时空失相干严重和时空分辨率低的缺点,通过干涉技术可实现优于毫米级微变形监测,采用GB雷达技术能精确测定被测物表面沿雷达视线向(LOS)的微量变形信息,其基本原理是:通过合成孔径雷达技术获取监测区域的二维影像,利用SF-CW技术提高雷达的距离向分辨率,通过比较影像中目标点的电磁波相位信息,采用干涉技术求取监测区域的变形量。
地基雷达进行工程监测具有以下优势:(1)利用干涉技术,地基雷达提取的形变精度高;(2)不受重访周期限制,地基雷达提取的形变具有连续性;(3)地基雷达可根据大型建筑的变形特点,选择合适的时间和位置进行观测,具有灵活性;(4)地基雷达的设站和周边环境易于控制和测量,测量误差具有可控性。但是同时还具有以下缺点,(1)地基雷达数据处理过程相对复杂;(2)地面干涉测量精度虽高,但易受到大气扰动或***噪声影响,而其中大气扰动造成的相位影响最大,导致测量精度低,有必要对其进行提取处理。
发明内容
本发明所解决的技术问题是克服现有的地基雷达数据处理方法中存在数据处理过程相对复杂,易受到大气扰动或***噪声影响的问题。本发明的基于地基合成孔径雷达干涉的大气相位屏提取方法,依据地基雷达差分干涉相位模型,通过最小二乘法拟合线性形变和时空滤波方法提取大气相位,该干涉相位模型充分考虑地基雷达成像机理以及差分干涉原理,从而提取大气相位屏,准确地提取受大气影响产生的相位对雷达干涉测量精度的提高具有重要作用,具有良好的应用前景。
为了解决上述目的,本发明所采用的技术方案是:
一种基于地基合成孔径雷达干涉的大气相位屏提取方法,其特征在于:包括以下步骤,
步骤(1),基于地基雷达差分干涉原理,建立差分相位模型,则第n景影像的第k个像素点的差分干涉相位如公式(1)所示,
其中,为第n景影像的第k个像素点两次雷达观测中像素在视线向形变相位、为第n景影像的第k个像素点在两次观测下大气延迟差的贡献值、为第n景影像的第k个像素点在两次观测像点的反射特性变化产生的相位贡献、为第n景影像的第k个像素点在干涉测量中,由于时间和空间去相干影响而产生的噪声以及SAR影像本身所产生的热噪声;
步骤(2),对雷达数据进行定标后,进行幅度信息比较,将地基雷达数据进行归一化处理,得到归一化像元值rij,如公式(2)所示,
其中,uij为原始影像第i行、第j列像素的幅值,并根据公式(2),得到归一化矩阵R,如公式(3)所示,
其中,m为矩阵R的行数,n为矩阵R的列数;
步骤(3),利用最小二乘法对连续观测的地基雷达数据进行配准,如公式(4)所示,
其中,Mx,y为主影像像元(x,y)强度值、为从影像像元(x2,y2)强度值、a0,a1,a2,b0,b1,b2为几何畸变参数,利用强度差的平方和最小原理进行配准,即Σvv=min,其中,v为主从影像强度值的较差;
步骤(4),将配准后的地基雷达进行差分干涉处理,如公式(5)所示,
其中,为主影像对应像点的相位值、为从影像对应像点的相位值,式中由SAR复数影像获取相位公式为Im(u)和Re(u)分别为雷达观测数据u的虚部和实部;
步骤(5),通过振幅离差指数法,提取永久散射体候选点DA,振幅离差的数学表达式,如公式(6)所示,
其中,μA、σA分别为振幅的均值和标准偏差,若DA的值不大于0.25,作为提取初始的永久散射体点目标的阈值;
步骤(6),将第n景影像的第k个像素点两次雷达观测中像素在视线向形变相位根据公式(7),分为线性视线向形变相位和非线性视线向形变相位的之和
其中,l为雷达波长,t为雷达干涉时间间隔,结合公式(1),则第n景影像的第k个像素点的差分干涉相位表示为公式(8),
其中,为残余相位;
步骤(7),通过时间域低通滤波去除噪声,去除估算线性形变后的残余相位为各差分干涉残余相位均值为作为主大气相位值,由于噪声在时间序列上是一个高频信号,因此,将为主大气相位从残余相位值中去除减掉在时间序列上通过低通滤波的剥离噪声;
步骤(8),将在空间上进行平滑结果,作为时间T上获取的从影像的大气相位APSslave,再将APSslave加上主大气相位值得到对应每个差分干涉图对应的大气相位;
步骤(9),当从差分干涉相位中去除大气相位屏后,根据公式(9)对影像上的每个像素点进行时序分析,
来重新计算时序影像上每个点的整体相干性,通过不大于阈值0.8以选择最终的永久散射体PS点;
步骤(10),利用每一个最终的永久散射体PS点的大气相位,通过Kriging插值法估计出对应于每一景雷达图像中每一个像元的大气相位,插值公式如公式(10)所示,
其中,APS(si)为第i个永久散射体PS点的大气相位、s0为未知像点位置、si为已知像点位置、M为估算大气相位值所需的永久散射体PS的点个数、λi为距离反比加权函数,表达式为di为待估像点与永久散射体PS点的距离值。
前述的基于地基合成孔径雷达干涉的大气相位屏提取方法,其特征在于:步骤(2)将地基雷达数据进行归一化处理,其中地基雷达数据的影像取景数量不少于30景。
前述的基于地基合成孔径雷达干涉的大气相位屏提取方法,其特征在于:步骤(7),通过时间域低通滤波去除噪声,采用均值滤波器。
前述的基于地基合成孔径雷达干涉的大气相位屏提取方法,其特征在于:步骤(5),通过振幅离差指数法,DA等于0.25,作为提取永久散射体候选点目标的阈值。
本发明的有益效果是:本发明的基于地基合成孔径雷达干涉的大气相位屏提取方法,依据地基雷达差分干涉相位模型,基于最小二乘法拟合线性形变和时空滤波方法提取大气相位,该干涉相位模型充分考虑地基雷达成像机理以及差分干涉原理,从而提取大气相位屏,本发明提出振幅离差指数、相干指数以及相位稳定性三级PS点选取方法,最后经Kriging插值获得每景雷达数据的大气相位屏,准确地提取受大气影响产生的相位对雷达干涉测量精度的提高具有重要作用,通过大坝实测数据验证本方法的正确性,具有一定的实际应用价值和通用性。
附图说明
图1是本发明的基于地基合成孔径雷达干涉的大气相位屏提取方法的流程图。
图2是本发明的相位离差的数学模拟结果的示意图。
图3是本发明的一实施例坝顶点P1的大气改正前后的位移时间序列图。
图4是本发明的一实施例坝底点P2的大气改正前后的位移时间序列图。
具体实施方式
下面将结合说明书附图,对本发明作进一步的说明。
本发明的基于地基合成孔径雷达干涉的大气相位屏提取方法,依据地基雷达差分干涉相位模型,基于最小二乘法拟合线性形变和时空滤波方法提取大气相位,该干涉相位模型充分考虑地基雷达成像机理以及差分干涉原理,从而提取大气相位屏,本发明提出振幅离差指数、相干指数以及相位稳定性三级PS点选取方法,最后经Kriging插值获得每景雷达数据的大气相位屏,准确地提取受大气影响产生的相位对雷达干涉测量精度的提高具有重要作用,如图1所示,具体包括以下步骤,
步骤(1),基于地基雷达差分干涉原理,建立差分相位模型,则第n景影像的第k个像素点的差分干涉相位如公式(1)所示,
其中,为第n景影像的第k个像素点两次雷达观测中像素在视线向形变相位、为第n景影像的第k个像素点在两次观测下大气延迟差的贡献值、为第n景影像的第k个像素点在两次观测像点的反射特性变化产生的相位贡献、为第n景影像的第k个像素点在干涉测量中,由于时间和空间去相干影响而产生的噪声以及SAR影像本身所产生的热噪声,实验中地基雷达的参数设置如下表1所示。
表1 地基雷达的参数设置表
步骤(2),对雷达数据进行定标后,进行幅度信息比较,将地基雷达数据进行归一化处理,这里的地基雷达数据,需要连续取景不少于30景,得到归一化像元值rij,如公式(2)所示,
其中,uij为原始影像第i行、第j列像素的幅值,并根据公式(2),得到归一化矩阵R,如公式(3)所示,
其中,m为矩阵R的行数,n为矩阵R的列数;
步骤(3),利用最小二乘法对连续观测的地基雷达数据进行配准,如公式(4)所示,
其中,Mx,y为主影像像元(x,y)强度值、为从影像像元(x2,y2)强度值、a0,a1,a2,b0,b1,b2为几何畸变参数,利用强度差的平方和最小原理进行配准,即Σvv=min,这里的v为主从影像强度值的较差;
步骤(4),将配准后的地基雷达进行差分干涉处理,如公式(5)所示,
其中,为主影像对应像点的相位值、为从影像对应像点的相位值,式中由SAR复数影像获取相位公式为Im(u)和Re(u)分别为雷达观测数据u的虚部和实部;
步骤(5),通过振幅离差指数法,提取永久散射体候选点DA,振幅离差的数学表达式,如公式(6)所示,
其中,μA、σA分别为振幅的均值和标准偏差,若DA的值不大于0.25,作为提取初始的永久散射体点目标的阈值,振幅离差方法是PS技术中的提取相干点目标的方法中的经典方法,它的主要思想是利用振幅离差和相位标准偏差统计特性上的关系,当即在高信噪比g/σn>4的条件下,相位标准偏差近似等于振幅离差DA,如图2所示,DA等于0.25,作为提取永久散射体候选点目标的阈值,
步骤(6),将第n景影像的第k个像素点两次雷达观测中像素在视线向形变相位根据公式(7),分为线性视线向形变相位和非线性视线向形变相位的之和
其中,l为雷达波长,t为雷达干涉时间间隔,结合公式(1),则第n景影像的第k个像素点的差分干涉相位表示为公式(8),
其中,为残余相位;
步骤(7),通过时间域低通滤波去除噪声,去除估算线性形变后的残余相位为各差分干涉残余相位均值为作为主大气相位值,由于噪声在时间序列上是一个高频信号,因此,将(减号,是将主大气相位从残余相位值中去除减掉)在时间序列上通过低通滤波剥离噪声;
步骤(8),将在空间上进行平滑结果作为时间T上获取的从影像的大气相位APSslave,然后再将APSslave加上主大气相位值(这个就是残余相位均值,也被估计为主大气相位值),得到对应每个差分干涉图对应的大气相位;
步骤(9),当从差分干涉相位中去除大气相位屏后,根据公式(9)对影像上的每个像素点进行时序分析,
来重新计算时序影像上每个点的整体相干性,通过不大于阈值0.8以选择最终的永久散射体PS点;
步骤(10),利用每一个最终的永久散射体PS点的大气相位,通过Kriging插值法估计出对应于每一景雷达图像中每一个像元的大气相位,插值公式如公式(10)所示,
其中,APS(si)为第i个永久散射体PS点的大气相位、s0为未知像点位置、si为已知像点位置、M为估算大气相位值所需的永久散射体PS点个数、λi为距离反比加权函数,表达式为di为待估像点与永久散射体PS点的距离值。
本发明的利用matlab程序处理,经过以上步骤对地基雷达进行大气相位提取后,可以对雷达视线向距离进行改正,进而得到更为准确的观测值,下面以具体实施例,经过永久散射体的选取后,取坝顶点P1及坝底点P2两点进行位移时间序列分析如下,如图3所示,坝顶点P1的大气改正前后的位移时间序列图,如图4所示,坝底点P2大气改正前后的位移时间序列图,大气的精确相位无法用实测数据直接获取,但由于测量误差的主值为大气扰动引起的相位值,因此,可以通过大坝实测的垂线数据来验证大气相位提取的精度,地基雷达数据经过大气改正后得到的视线向位移同大坝垂线数据的较差,明显地要小于未经改正的结果,从而也验证本方法所提取大气相位是准确的,具有一定的实用性。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (4)
1.基于地基合成孔径雷达干涉的大气相位屏提取方法,其特征在于:包括以下步骤,
步骤(1),基于地基雷达差分干涉原理,建立差分相位模型,则第n景影像的第k个像素点的差分干涉相位如公式(1)所示,
其中,为第n景影像的第k个像素点两次雷达观测中像素在视线向形变相位、为第n景影像的第k个像素点在两次观测下大气延迟差的贡献值、为第n景影像的第k个像素点在两次观测像点的反射特性变化产生的相位贡献、为第n景影像的第k个像素点在干涉测量中,由于时间和空间去相干影响而产生的噪声以及SAR影像本身所产生的热噪声;
步骤(2),对雷达数据进行定标后,进行幅度信息比较,将地基雷达数据进行归一化处理,得到归一化像元值rij,如公式(2)所示,
其中,uij为原始影像第i行、第j列像素的幅值,并根据公式(2),得到归一化矩阵R,如公式(3)所示,
其中,m为矩阵R的行数,n为矩阵R的列数;
步骤(3),利用最小二乘法对连续观测的地基雷达数据进行配准,如公式(4)所示,
其中,Mx,y为主影像像元(x,y)强度值、为从影像像元(x2,y2)强度值、a0,a1,a2,b0,b1,b2为几何畸变参数,利用强度差的平方和最小原理进行配准,即Σvv=min,其中,v为主从影像强度值的较差;
步骤(4),将配准后的地基雷达进行差分干涉处理,如公式(5)所示,
其中,为主影像对应像点的相位值、为从影像对应像点的相位值,式中由SAR复数影像获取相位公式为Im(u)和Re(u)分别为雷达观测数据u的虚部和实部;
步骤(5),通过振幅离差指数法,提取永久散射体候选点DA,振幅离差的数学表达式,如公式(6)所示,
其中,μA、σA分别为振幅的均值和标准偏差,若DA的值不大于0.25,作为提取初始的永久散射体点目标的阈值;
步骤(6),将第n景影像的第k个像素点两次雷达观测中像素在视线向形变相位根据公式(7),分为线性视线向形变相位和非线性视线向形变相位的之和
其中,l为雷达波长,t为雷达干涉时间间隔,结合公式(1),则第n景影像的第k个像素点的差分干涉相位表示为公式(8),
其中,为残余相位;
步骤(7),通过时间域低通滤波去除噪声,去除估算线性形变后的残余相位为各差分干涉残余相位均值为作为主大气相位值,由于噪声在时间序列上是一个高频信号,因此,将为主大气相位从残余相位值中去除减掉在时间序列上通过低通滤波的剥离噪声;
步骤(8),将在空间上进行平滑结果,作为时间T上获取的从影像的大气相位APSslave,再将APSslave加上主大气相位值得到对应每个差分干涉图对应的大气相位;
步骤(9),当从差分干涉相位中去除大气相位屏后,根据公式(9)对影像上的每个像素点进行时序分析,
来重新计算时序影像上每个点的整体相干性,通过不大于阈值0.8以选择最终的永久散射体PS点;
步骤(10),利用每一个最终的永久散射体PS点的大气相位,通过Kriging插值法估计出对应于每一景雷达图像中每一个像元的大气相位,插值公式如公式(10)所示,
其中,APS(si)为第i个永久散射体PS点的大气相位、s0为未知像点位置、si为已知像点位置、M为估算大气相位值所需的永久散射体PS的点个数、λi为距离反比加权函数,表达式为di为待估像点与永久散射体PS点的距离值。
2.根据权利要求1所述的基于地基合成孔径雷达干涉的大气相位屏提取方法,其特征在于:步骤(2)将地基雷达数据进行归一化处理,其中地基雷达数据的影像取景数量不少于30景。
3.根据权利要求1所述的基于地基合成孔径雷达干涉的大气相位屏提取方法,其特征在于:步骤(7),通过时间域低通滤波去除噪声,采用均值滤波器。
4.根据权利要求1所述的基于地基合成孔径雷达干涉的大气相位屏提取方法,其特征在于:步骤(5),通过振幅离差指数法,DA等于0.25,作为提取永久散射体候选点目标的阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410633483.9A CN104360332B (zh) | 2014-11-11 | 2014-11-11 | 基于地基合成孔径雷达干涉的大气相位屏提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410633483.9A CN104360332B (zh) | 2014-11-11 | 2014-11-11 | 基于地基合成孔径雷达干涉的大气相位屏提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104360332A true CN104360332A (zh) | 2015-02-18 |
CN104360332B CN104360332B (zh) | 2017-02-01 |
Family
ID=52527612
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410633483.9A Expired - Fee Related CN104360332B (zh) | 2014-11-11 | 2014-11-11 | 基于地基合成孔径雷达干涉的大气相位屏提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104360332B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105678716A (zh) * | 2016-02-25 | 2016-06-15 | 内蒙古工业大学 | 一种地基sar大气干扰相位校正方法及装置 |
CN105783754A (zh) * | 2016-03-14 | 2016-07-20 | 河海大学 | 基于三维激光扫描的GBInSAR三维位移场提取方法 |
CN106023157A (zh) * | 2016-05-10 | 2016-10-12 | 电子科技大学 | 一种基于sar图像的山区地表微形变信息提取方法 |
CN106772342A (zh) * | 2017-01-11 | 2017-05-31 | 西南石油大学 | 一种适用于大梯度地表沉降监测的时序差分雷达干涉方法 |
CN107218923A (zh) * | 2017-05-23 | 2017-09-29 | 北京东方至远科技股份有限公司 | 基于PS‑InSAR技术的地铁沿线周边环境历史沉降风险评估方法 |
CN108549080A (zh) * | 2018-02-28 | 2018-09-18 | 中国电力科学研究院有限公司 | 一种输电杆塔位置提取方法及*** |
CN108830242A (zh) * | 2018-06-22 | 2018-11-16 | 北京航空航天大学 | 基于卷积神经网络的sar图像海洋目标分类检测方法 |
CN110334771A (zh) * | 2019-07-09 | 2019-10-15 | 东北大学 | 一种基于四阈值ps点的露天矿gb-sar监测数据误差修正方法 |
CN112986990A (zh) * | 2021-02-04 | 2021-06-18 | 中国地质大学(北京) | 一种大气相位改正方法及*** |
CN113791412A (zh) * | 2021-09-13 | 2021-12-14 | 中国电建集团贵阳勘测设计研究院有限公司 | 一种基于时移的堤防隐患排查方法 |
CN114624708A (zh) * | 2022-05-16 | 2022-06-14 | 中山大学 | 一种复杂环境下的大气校正方法及*** |
CN116482684A (zh) * | 2023-06-21 | 2023-07-25 | 深圳市城市公共安全技术研究院有限公司 | 区域五维成像方法、装置、设备及存储介质 |
CN117849751A (zh) * | 2024-01-10 | 2024-04-09 | 兰州理工大学 | 一种针对星凸形不规则形状扩展目标的跟踪方法及*** |
CN117968589B (zh) * | 2024-03-20 | 2024-08-02 | 苏州科技大学 | 一种地基合成孔径雷达影像质量评价方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104021284A (zh) * | 2014-05-23 | 2014-09-03 | 河海大学 | 一种地面沉降危害性评价方法 |
CN104123470A (zh) * | 2014-07-25 | 2014-10-29 | 首都师范大学 | 一种优化地面沉降监测网的方法 |
-
2014
- 2014-11-11 CN CN201410633483.9A patent/CN104360332B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104021284A (zh) * | 2014-05-23 | 2014-09-03 | 河海大学 | 一种地面沉降危害性评价方法 |
CN104123470A (zh) * | 2014-07-25 | 2014-10-29 | 首都师范大学 | 一种优化地面沉降监测网的方法 |
Non-Patent Citations (6)
Title |
---|
ALESSANDRO FERRETTI,等: ""Nonlinear Subsidence Rate Estimation Using"", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
CARLO COLESANTI,等: ""Monitoring landslides and tectonic motions with the Permanent"", 《ENGINEERING GEOLOGY》 * |
张祥: ""地基SFCW SAR差分干涉测量技术研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
林忠辉,等: ""中国陆地区域气象要素的空间插值"", 《地理学报》 * |
罗小军: ""永久散射体雷达差分干涉理论及在上海地面沉降监测中的应用"", 《中国博士学位论文全文数据库 基础科学辑》 * |
范洪冬: ""InSAR 若干关键算法及其在地表沉降监测中的应用研究"", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105678716B (zh) * | 2016-02-25 | 2018-09-18 | 内蒙古工业大学 | 一种地基sar大气干扰相位校正方法及装置 |
CN105678716A (zh) * | 2016-02-25 | 2016-06-15 | 内蒙古工业大学 | 一种地基sar大气干扰相位校正方法及装置 |
CN105783754A (zh) * | 2016-03-14 | 2016-07-20 | 河海大学 | 基于三维激光扫描的GBInSAR三维位移场提取方法 |
CN105783754B (zh) * | 2016-03-14 | 2018-12-14 | 河海大学 | 基于三维激光扫描的GBInSAR三维位移场提取方法 |
CN106023157A (zh) * | 2016-05-10 | 2016-10-12 | 电子科技大学 | 一种基于sar图像的山区地表微形变信息提取方法 |
CN106023157B (zh) * | 2016-05-10 | 2018-08-07 | 电子科技大学 | 一种基于sar图像的山区地表微形变信息提取方法 |
CN106772342B (zh) * | 2017-01-11 | 2020-02-07 | 西南石油大学 | 一种适用于大梯度地表沉降监测的时序差分雷达干涉方法 |
CN106772342A (zh) * | 2017-01-11 | 2017-05-31 | 西南石油大学 | 一种适用于大梯度地表沉降监测的时序差分雷达干涉方法 |
CN107218923A (zh) * | 2017-05-23 | 2017-09-29 | 北京东方至远科技股份有限公司 | 基于PS‑InSAR技术的地铁沿线周边环境历史沉降风险评估方法 |
CN108549080A (zh) * | 2018-02-28 | 2018-09-18 | 中国电力科学研究院有限公司 | 一种输电杆塔位置提取方法及*** |
CN108830242A (zh) * | 2018-06-22 | 2018-11-16 | 北京航空航天大学 | 基于卷积神经网络的sar图像海洋目标分类检测方法 |
CN110334771A (zh) * | 2019-07-09 | 2019-10-15 | 东北大学 | 一种基于四阈值ps点的露天矿gb-sar监测数据误差修正方法 |
CN110334771B (zh) * | 2019-07-09 | 2023-02-28 | 东北大学 | 一种基于四阈值ps点的露天矿gb-sar监测数据误差修正方法 |
CN112986990A (zh) * | 2021-02-04 | 2021-06-18 | 中国地质大学(北京) | 一种大气相位改正方法及*** |
CN112986990B (zh) * | 2021-02-04 | 2023-02-17 | 中国地质大学(北京) | 一种大气相位改正方法及*** |
CN113791412A (zh) * | 2021-09-13 | 2021-12-14 | 中国电建集团贵阳勘测设计研究院有限公司 | 一种基于时移的堤防隐患排查方法 |
CN114624708A (zh) * | 2022-05-16 | 2022-06-14 | 中山大学 | 一种复杂环境下的大气校正方法及*** |
CN116482684A (zh) * | 2023-06-21 | 2023-07-25 | 深圳市城市公共安全技术研究院有限公司 | 区域五维成像方法、装置、设备及存储介质 |
CN116482684B (zh) * | 2023-06-21 | 2023-08-22 | 深圳市城市公共安全技术研究院有限公司 | 区域五维成像方法、装置、设备及存储介质 |
CN117849751A (zh) * | 2024-01-10 | 2024-04-09 | 兰州理工大学 | 一种针对星凸形不规则形状扩展目标的跟踪方法及*** |
CN117849751B (zh) * | 2024-01-10 | 2024-07-16 | 兰州理工大学 | 一种针对星凸形不规则形状扩展目标的跟踪方法及*** |
CN117968589B (zh) * | 2024-03-20 | 2024-08-02 | 苏州科技大学 | 一种地基合成孔径雷达影像质量评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104360332B (zh) | 2017-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104360332A (zh) | 基于地基合成孔径雷达干涉的大气相位屏提取方法 | |
CN104123464B (zh) | 一种高分辨率InSAR时序分析反演地物高程与地面沉降量的方法 | |
CN106772342B (zh) | 一种适用于大梯度地表沉降监测的时序差分雷达干涉方法 | |
CN103698764B (zh) | 一种稀疏采样条件下的干涉合成孔径雷达成像方法 | |
CN116047519B (zh) | 一种基于合成孔径雷达干涉测量技术的选点方法 | |
CN104698460A (zh) | 一种双频共面合成孔径雷达海面风场反演方法 | |
CN104020456A (zh) | 一种基于多特显点线阵成像雷达***幅相误差校正方法 | |
CN103714546B (zh) | 一种成像光谱仪的数据处理方法 | |
CN103822598A (zh) | 地基sar在时间去相关严重区域的形变监测方法 | |
CN105467390A (zh) | 一种基于地基InSAR的桥梁形变近距离监测方法 | |
CN102288943A (zh) | 基于两视实图像处理的单通道sar-gmti方法 | |
CN111398959A (zh) | 基于地表应力应变模型的InSAR时序地表形变监测方法 | |
CN104678392A (zh) | 一种基于稳定像元加权的地基sar大气延迟改正方法 | |
CN104200471A (zh) | 基于自适应权值图像融合的sar图像变化检测方法 | |
CN107064933B (zh) | 基于循环谱估计的sar层析建筑物高度的方法 | |
CN105866776A (zh) | 一种地基sar的高质量动态ps点选择方法 | |
CN112444188B (zh) | 一种多视角InSAR海堤高精度三维形变测量方法 | |
CN113281742A (zh) | 一种基于滑坡形变信息和气象数据的sar滑坡预警方法 | |
CN113281749B (zh) | 一种顾及同质性的时序InSAR高相干点选取方法 | |
CN116338607B (zh) | 时间域和空间域两步式InSAR对流层延迟矫正方法 | |
CN111650570A (zh) | 一种地基干涉雷达三维大气校正方法及*** | |
CN103630903A (zh) | 基于顺轨干涉sar测量海面流场径向速度的方法 | |
CN105678716B (zh) | 一种地基sar大气干扰相位校正方法及装置 | |
CN108872989A (zh) | 一种基于最大周期图的PS-InSAR精确搜索方法 | |
CN114594479B (zh) | 一种全散射体FS-InSAR方法及*** |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170201 Termination date: 20191111 |
|
CF01 | Termination of patent right due to non-payment of annual fee |