CN111487622A - 输电铁塔变形监测方法、装置、计算机设备和存储介质 - Google Patents
输电铁塔变形监测方法、装置、计算机设备和存储介质 Download PDFInfo
- Publication number
- CN111487622A CN111487622A CN202010541550.XA CN202010541550A CN111487622A CN 111487622 A CN111487622 A CN 111487622A CN 202010541550 A CN202010541550 A CN 202010541550A CN 111487622 A CN111487622 A CN 111487622A
- Authority
- CN
- China
- Prior art keywords
- power transmission
- scatterer
- sar image
- transmission tower
- coordinates
- 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.)
- Pending
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/9004—SAR image acquisition techniques
- G01S13/9005—SAR image acquisition techniques with optical processing of the SAR signals
-
- 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/9021—SAR image post-processing techniques
- G01S13/9027—Pattern recognition for feature extraction
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)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本申请涉及一种输电铁塔变形监测方法、装置、计算机设备和存储介质。通过获取包括不同拍摄时间的输电铁塔图像的SAR影像,根据SAR影像中的坐标值,确定输电铁塔在影像中的覆盖范围,在覆盖范围内提取输电铁塔表面的散射体对应的坐标,并根据该坐标,得到对应的干涉相位值,还可以根据上述坐标以及输电铁塔底部中心的坐标,得到散射体和底部中心的高度差,在去除了散射体的干涉相位中的高度误差相位和热膨胀相位后,可以根据更新后的散射体的干涉相位值,对输电铁塔进行变形监测。相较于传统的通过人工巡视的方式,本方案根据去除了误差之后的输电铁塔表面散射体的干涉相位值,监测输电铁塔的变形情况,提高了输电铁塔变形监测的可靠性。
Description
技术领域
本申请涉及遥感监测技术领域,特别是涉及一种输电铁塔变形监测方法、装置、计算机设备和存储介质。
背景技术
输电铁塔是一种重要的电力设备,主要结构形式为钢结构,是一种高耸的构筑物,主要负责高压、特高压电力的传输,在输电沿途中,输电铁塔会安装在很多环境条件恶劣、地质地形复杂或者气候多变的地区,同时输电线路周边的施工、采掘作业活动,容易导致输电铁塔及基础的倾斜或沉降,威胁输电线路的安全稳定运行,为了保证输电线路的安全的稳定,对输电铁塔的状态监测显得十分必要。目前,对于输电铁塔的状态变化监测的方法通常是采用巡线员和护线员相结合的人工巡视方式,然而,该方式工作量大、巡视周期长、视野范围小,不能准确地监测输电铁塔的状态变化。
因此,目前的输电铁塔的变形监测方法存在监测可靠性低的缺陷。
发明内容
基于此,有必要针对上述技术问题,提供一种能够提高输电铁塔的变形监测可靠性的输电铁塔变形监测方法、装置、计算机设备和存储介质。
一种输电铁塔变形监测方法,所述方法包括:
获取所述输电铁塔的SAR影像;所述SAR影像中包括至少两个对应不同拍摄时间的输电铁塔图像;
根据所述输电铁塔在SAR影像中的坐标值,得到所述输电铁塔在SAR影像中的覆盖范围;
在所述覆盖范围内获取所述输电铁塔表面的散射体对应的坐标;
根据所述散射体的坐标,从预设相位表中获取所述散射体的坐标对应的干涉相位值;所述预设相位表包括所述坐标与相位的对应关系;
根据在所述SAR影像中的所述散射体的坐标以及所述输电铁塔的底部中心坐标,得到所述散射体与所述输电铁塔底部中心的高度差;
基于所述高度差、所述SAR影像对应的雷达卫星的电磁波波长、所述雷达卫星的垂直基线长度以及所述雷达卫星与所述输电铁塔的距离,得到所述散射体的高度误差以及高度误差相位;
基于两个所述输电铁塔图像对应的拍摄时间的环境温度差、所述输电铁塔的材质对应的热膨胀系数,得到所述散射体的热膨胀相位;
将所述高度误差相位以及所述热膨胀相位从所述散射体的干涉相位值中去除,得到更新后的散射体的干涉相位值;根据更新后的散射体的干涉相位值,对所述输电铁塔进行变形监测。
在其中一个实施例中,所述根据所述输电铁塔在SAR影像中的坐标值,得到所述输电铁塔在SAR影像中的覆盖范围,包括:
获取所述SAR影像中所述输电铁塔的底部中心坐标值;
获取所述SAR影像中所述输电铁塔的距离向长度;
获取所述SAR影像中所述输电铁塔的方位向长度;
基于所述底部中心坐标值、所述距离向长度以及所述方位向长度,得到所述输电铁塔在SAR影像中的覆盖范围。
在其中一个实施例中,所述根据所述输电铁塔在SAR影像中的坐标值,得到所述输电铁塔在SAR影像中的覆盖范围之前,还包括:
基于所述雷达卫星的轨道信息以及所述SAR影像信息,生成所述SAR影像与地面的坐标映射关系,以在所述坐标映射关系中,获取所述输电铁塔在SAR影像中的坐标值。
在其中一个实施例中,所述获取所述SAR影像中所述输电铁塔的距离向长度,包括:
获取所述输电铁塔的实际高度;
获取所述雷达卫星的入射角;
根据所述实际高度、所述入射角以及所述SAR影像的距离向分辨率,获取所述SAR影像中所述输电铁塔的距离向长度。
在其中一个实施例中,所述获取所述SAR影像中所述输电铁塔的方位向长度,包括:
获取所述输电铁塔的实际宽度;
根据所述实际宽度以及所述SAR影像的方位向分辨率,获取所述SAR影像中所述输电铁塔的方位向长度。
在其中一个实施例中,所述在所述覆盖范围内获取所述输电铁塔表面的散射体对应的坐标,包括:
若所述散射体的光强大于预设光强阈值,获取所述散射体对应的坐标,作为散射体的像元坐标;
将多张不同拍摄时间的SAR影像进行非相干叠加,得到所述多张不同拍摄时间的SAR影像对应的平均光强影像;
对所述平均光强影像进行过采样处理,得到更新后的平均光强影像;根据更新后的平均光强影像中每个像素的坐标、每个像素的光强以及二维辛格函数中所述像素对应的坐标的数值,得到所述更新后的平均光强影像与所述二维辛格函数的相关系数;
获取所述相关系数的最大值对应的坐标,作为所述平均光强影像中散射体的子像元坐标;
基于所述散射体的像元坐标、所述子像元坐标以及预设过采样系数,得到所述输电铁塔中的散射体对应的坐标。
在其中一个实施例中,所述根据在所述SAR影像中的所述散射体的坐标以及所述输电铁塔的底部中心坐标,得到所述散射体与所述输电铁塔底部中心的高度差,包括:
获取所述散射体与所述输电铁塔底部中心的距离向坐标的差值;
基于所述差值以及所述SAR影像的距离向分辨率,得到所述散射体与所述输电铁塔底部中心的高度差。
一种输电铁塔变形监测装置,所述装置包括:
坐标获取模块,用于获取所述输电铁塔的SAR影像;所述SAR影像中包括至少两个对应不同拍摄时间的输电铁塔图像;根据所述输电铁塔在SAR影像中的坐标值,得到所述输电铁塔在SAR影像中的覆盖范围;在所述覆盖范围内获取所述输电铁塔表面的散射体对应的坐标;
相位获取模块,用于根据所述散射体的坐标,从预设相位表中获取所述散射体的坐标对应的干涉相位值;所述预设相位表包括所述坐标与相位的对应关系;
误差获取模块,用于根据在所述SAR影像中的所述散射体的坐标以及所述输电铁塔的底部中心坐标,得到所述散射体与所述输电铁塔底部中心的高度差;基于所述高度差、所述SAR影像对应的雷达卫星的电磁波波长、所述雷达卫星的垂直基线长度以及所述雷达卫星与所述输电铁塔的距离,得到所述散射体的高度误差以及高度误差相位;基于两个所述输电铁塔图像对应的拍摄时间的环境温度差、所述输电铁塔的材质对应的热膨胀系数,得到所述散射体的热膨胀相位;
监测模块,用于将所述高度误差相位以及所述热膨胀相位从所述散射体的干涉相位值中去除,得到更新后的散射体的干涉相位值;根据更新后的散射体的干涉相位值,对所述输电铁塔进行变形监测。
一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现上述的方法的步骤。
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现上述的方法的步骤。
上述输电铁塔变形监测方法、装置、计算机设备和存储介质,通过获取包括不同拍摄时间的输电铁塔图像的SAR影像,根据SAR影像中的坐标值,确定输电铁塔在影像中的覆盖范围,在覆盖范围内提取输电铁塔表面的散射体对应的坐标,并根据该坐标,得到对应的干涉相位值,还可以根据上述坐标以及输电铁塔底部中心的坐标,得到散射体和底部中心的高度差,在去除了散射体的干涉相位中的高度误差相位和热膨胀相位后,可以根据更新后的散射体的干涉相位值,对输电铁塔进行变形监测。相较于传统的通过人工巡视的方式,本方案通过在输电铁塔的SAR影像中获取输电铁塔表面的散射体对应的坐标以及干涉相位,并根据去除了误差之后的干涉相位值,监测输电铁塔的变形情况,提高了输电铁塔变形监测的可靠性。
附图说明
图1为一个实施例中输电铁塔变形监测方法的应用环境图;
图2为一个实施例中输电铁塔变形监测方法的流程示意图;
图3为另一个实施例中输电铁塔变形监测方法的流程示意图;
图4为一个实施例中输电铁塔的SAR影像示意图;
图5为一个实施例中散射体干涉相位的示意图;
图6为一个实施例中去除误差后的散射体干涉相位的示意图;
图7为一个实施例中输电铁塔变形监测装置的结构框图;
图8为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
本申请提供的输电铁塔变形监测方法,可以应用于如图1所示的应用环境中。其中,终端102通过网络与服务器104进行通信。终端102可以从服务器104中获取上述输电铁塔的SAR影像,并可以获取SAR影像中输电铁塔的覆盖范围,从而在范围内得到输电铁塔表面的散射体对应的坐标,终端102还可以根据坐标得到散射体的干涉相位值,并根据去除了误差后的干涉相位值,监测输电铁塔的变形情况。其中,终端102可以但不限于是各种个人计算机、笔记本电脑、智能手机和平板电脑,服务器104可以用独立的服务器或者是多个服务器组成的服务器集群来实现。
在一个实施例中,如图2所示,提供了一种输电铁塔变形监测方法,以该方法应用于图1中的终端为例进行说明,包括以下步骤:
步骤S202,获取输电铁塔的SAR影像;SAR影像中包括至少两个对应不同拍摄时间的输电铁塔图像;根据输电铁塔在SAR影像中的坐标值,得到输电铁塔在SAR影像中的覆盖范围;在覆盖范围内获取输电铁塔表面的散射体对应的坐标。
其中,SAR图像可以是SAR(Synthetic Aperture Radar,合成孔径雷达)拍摄的图像,SAR是一种主动式的对地观测系统,可安装在飞机、卫星、宇宙飞船等飞行平台上,全天时、全天候对地实施观测、并具有一定的地表穿透能力。因此,SAR系统在灾害监测、环境监测、海洋监测、资源勘查、农作物估产、测绘和军事等方面的应用上具有独特的优势,可发挥其他遥感手段难以发挥的作用,因此越来越受到世界各国的重视。SAR可以拍摄多组输电铁塔的图像,例如可以拍摄不同角度的输电铁塔的图像,也可以是不同时间拍摄的输电铁塔的图像,SAR图像可以存储在上述服务器104中,其中服务器104中可以至少存储有两个不同拍摄时间的输电铁塔图像。
终端102可以在SAR图像中建立坐标系,具体地,可以基于雷达卫星的轨道信息和上述SAR影像的参数,建立SAR图像与大地坐标之间的映射关系,其中雷达卫星可以是拍摄上述SAR图像的卫星。终端102可以在SAR图像中,根据坐标系,获取上述输电铁塔的坐标值,从而得到输电铁塔在SAR影像中的覆盖范围,具体地,可以根据输电铁塔的底部中心坐标、长宽等数据得到。
终端102在得到SAR影像中输电铁塔的覆盖范围后,可以在范围内获取上述输电铁塔表面的散射体对应的坐标,输电铁塔在SAR影像中可以表现为明显的亮斑,雷达回拨信号较为稳定,可以作为高信杂比的散射体,信杂比可以是在标准照度下摄像机亮度通道输出信号的峰值与视频杂波的有效值之比。终端102可以基于散射体的光强度信息实现铁塔表面散射体的识别以及元坐标提取,然后利用二维sinc函数,即二维辛格函数匹配方法实现铁塔表面散射体子像素坐标提取,最终获得铁塔表面散射体在SAR影像中的精确坐标。二维辛格函数可以是一维辛格函数的拓展,可以是正弦函数和单调递减函数1/x的乘积。
步骤S204,根据散射体的坐标,从预设相位表中获取散射体的坐标对应的干涉相位值;预设相位表包括坐标与相位的对应关系。
其中,散射体的坐标可以通过上述步骤S202得到,散射体的坐标可以包括方位向坐标以及距离向坐标,其中,方位向可以是沿着上述雷达卫星轨道的方向,距离向可以是沿着雷达卫星发射的雷达波的方向。终端102可以存储有一预设相位表,上述SAR影像中的每个坐标可以对应一干涉相位值,即上述每个散射体的坐标可以从预设的相位表中得到对应的相位,终端102还可以将上述每个散射体坐标以及每个坐标对应的干涉相位值以图表方式显示。其中,干涉相位值可以是对于一个波,特定的时刻在它循环中的位置,一种它是否在波峰、波谷或它们之间的某点的标度。相位描述信号波形变化的度量,通常以度(角度)作为单位,也称作相角。
步骤S206,根据在SAR影像中的散射体的坐标以及输电铁塔的底部中心坐标,得到散射体与输电铁塔底部中心的高度差;基于高度差、SAR影像对应的雷达卫星的电磁波波长、雷达卫星的垂直基线长度以及雷达卫星与输电铁塔的距离,得到散射体的高度误差以及高度误差相位;基于两个输电铁塔图像对应的拍摄时间的环境温度差、输电铁塔的材质对应的热膨胀系数,得到散射体的热膨胀相位。
其中,高度差可以是输电铁塔表面的散射体与输电铁塔的底部的实际高度差,终端102可以根据上述SAR影像中散射体的坐标以及输电铁塔的底部中心坐标,再利用上述雷达卫星的入射角以及距离分辨率,得到散射体和输电铁塔底部中心的高度差。其中,距离分辨率可以是在雷达图像中,当两个目标位于同一方位角时,但与雷达的距离不同时,二者被雷达区分出来的最小距离称为距离分辨率。即雷达距离分辨率定义为雷达分辨两个近距离目标的能力。
终端102在得到上述高度差后,可以基于上述高度差、SAR影像对应的雷达卫星的电磁波波长、雷达卫星的垂直基线长度以及雷达卫星至输电铁塔的距离,得到散射体的高度误差以及高度误差相位。具体地,终端102可以利用上述SAR影像的轨道数据和铁塔表面散射体对应的高度差,得到散射体对应的高度误差及其相位,如下式所示:
其中,Δh可以是铁塔表面散射体相对铁塔底部的高度差,λ为雷达卫星电磁波波长,B⊥为干涉图垂直基线长度,R为卫星与上述输电铁塔的距离。SAR系统中,天线的相对位置向量B被叫做空间基线,而垂直基线可以是基线在垂直视线方向的分量。
由于铁塔表面监测点会受热膨胀影响产生相位变化,因此终端102可以基于上述SAR影像中包含的两个不同输电铁塔图像对应的拍摄时间的环境温度差,以及上述输电铁塔的材质对应的热膨胀系数,得到上述散射体的热膨胀相位。具体地,终端102可以将铁塔表面监测点受热膨胀效应影响产生的相位变化量表示为:
其中,temp为热膨胀因子,可通过查询铁塔材质对应的热膨胀系数获取。ΔTemp为上述不同拍摄时间的SAR影像所对应的外界环境温度差值,外界环境温度可通过查询当地气象网站获取。上述两个不同时间拍摄的铁塔图像可以是以重叠方式叠加,具体地,可以是将图像中的输电铁塔进行重合,得到上述SAR影像。
步骤S208,将高度误差相位以及热膨胀相位从散射体的干涉相位值中去除,得到更新后的散射体的干涉相位值;根据更新后的散射体的干涉相位值,对输电铁塔进行变形监测。
其中,高度误差相位可以是上述SAR影像中,输电铁塔表面的散射体存在的高度误差对应的相位,热膨胀相位可以是上述输电铁塔受热膨胀影响产生的误差相位。终端102可以在求得上述高度误差相位以及热膨胀相位后,将高度误差相位和热膨胀相位从散射体的干涉相位值中去除,得到更新后的散射体的干涉相位值,即精确化后的相位,并且可以基于更新后的干涉相位值,对输电铁塔进行变形监测,例如当干涉相位值越低时,输电铁塔的变形越严重。
上述输电铁塔变形监测方法中,通过获取包括不同拍摄时间的输电铁塔图像的SAR影像,根据SAR影像中的坐标值,确定输电铁塔在影像中的覆盖范围,在覆盖范围内提取输电铁塔表面的散射体对应的坐标,并根据该坐标,得到对应的干涉相位值,还可以根据上述坐标以及输电铁塔底部中心的坐标,得到散射体和底部中心的高度差,在去除了散射体的干涉相位中的高度误差相位和热膨胀相位后,可以根据更新后的散射体的干涉相位值,对输电铁塔进行变形监测。相较于传统的通过人工巡视的方式,本方案通过在输电铁塔的SAR影像中获取输电铁塔表面的散射体对应的坐标以及干涉相位,并根据去除了误差之后的干涉相位值,监测输电铁塔的变形情况,提高了输电铁塔变形监测的可靠性。
在一个实施例中,根据输电铁塔在SAR影像中的坐标值,得到输电铁塔在SAR影像中的覆盖范围,包括:获取SAR影像中输电铁塔的底部中心坐标值;获取SAR影像中输电铁塔的距离向长度;获取SAR影像中所述输电铁塔的方位向长度;基于底部中心坐标值、距离向长度以及方位向长度,得到输电铁塔在SAR影像中的覆盖范围。
本实施例中,输电铁塔在SAR影像中的坐标值可以根据预先建立的SAR影像和地面的坐标映射关系中得到,坐标映射关系可以包括距离向以及方位向,终端102可以根据上述输电铁塔在SAR影像中的坐标值,得到输电铁塔在SAR影像中的覆盖范围。具体地,终端102可以获取SAR影像中输电铁塔的底部中心坐标值,该坐标值可以表示为(XRef,YRef),其中XRef为距离向坐标,YRef为方位向坐标,其中距离向坐标可以是上述底部中心在距离向的坐标值,方位向坐标可以是上述底部中心在方位向上的坐标值。终端102还可以获取上述SAR影像中,输电铁塔的距离向长度,即输电铁塔在SAR影像中对应的距离向像素数,例如根据输电铁塔的高度获取。终端102还可以获取上述SAR影像中,输电铁塔的方位向长度,即输电铁塔在SAR影像中对一个内的方位向像素数,例如根据输电铁塔的宽度获取。其中,方位向可以是沿着上述雷达卫星轨道的方向,距离向可以是沿着雷达卫星发射的雷达波的方向。终端102可以基于获取到的输电铁塔在SAR图像中的底部中心坐标值、距离向长度以及方位向长度,得到上述输电铁塔在SAR影像中的覆盖范围,例如可以以底部中心坐标为基准,按照输电铁塔的方位向长度和距离向长度,计算输电铁塔的覆盖范围。
通过本实施例,终端102可以根据输电铁塔的底部中心坐标值、方位向长度和距离向长度,得到输电铁塔在SAR影像中的覆盖范围,从而可以在该范围内进行相应监测,提高了输电铁塔变形监测的可靠性。
在一个实施例中,获取SAR影像中输电铁塔的距离向长度,包括:获取输电铁塔的实际高度;获取雷达卫星的入射角;根据实际高度、入射角以及SAR影像的距离向分辨率,获取SAR影像中输电铁塔的距离向长度。
本实施例中,距离向长度可以是输电铁塔在SAR影像中,距离向的像素数,其中距离向可以是沿着雷达卫星发射的雷达波的方向。终端102可以获取输电铁塔的实际高度,还可以获取上述雷达卫星的入射角,并且根据上述实际高度、入射角以及上述SAR影像的距离向分辨率,得到SAR影像中输电铁塔的距离向长度。具体地,假设输电铁塔高度为H0,雷达卫星入射角为θ,其在SAR影像中对应的距离向像素数XLPix可标识为:
XLPix=H0/(tanθ*RX)
式中RX为SAR影像距离向分辨率,入射角和距离向分辨率参数可从影像参数文件中获取。其中雷达卫星可以是拍摄上述SAR影像对应的卫星。
通过本实施例,终端102可以根据输电铁塔的实际高度、雷达卫星的入射角以及SAR影像的距离向分辨率,得到输电铁塔的距离向长度,从而可以基于距离向长度得到输电铁塔在SAR影像中的覆盖范围,实现在覆盖范围内对输电铁塔进行变形监测,提高了输电铁塔变形监测的可靠性。
在一个实施例中,获取SAR影像中输电铁塔的方位向长度,包括:获取输电铁塔的实际宽度;根据实际宽度以及SAR影像的方位向分辨率,获取SAR影像中输电铁塔的方位向长度。
本实施例中,方位向长度可以是输电铁塔在SAR影像中,方位向的像素数,其中,方位向可以是沿着上述雷达卫星轨道的方向。终端102可以获取输电铁塔的实际宽度,并根据上述实际宽度,以及上述SAR影像的方位向分辨率,得到SAR影像中输电铁塔的方位向长度。具体地,假设输电铁塔的宽度为L,其在SAR影像中对应的方位向像素数YLPix可标识为:
YLPix=L/RY
式中RY为SAR影像方位向分辨率,可从SAR影像参数文件中获取。
通过本实施例,终端102可以根据输电铁塔的实际宽度,以及上述SAR影像的方位向分辨率,得到SAR影像中输电铁塔的方位向长度,从而可以基于方位向长度得到输电铁塔在SAR影像中的覆盖范围,实现在覆盖范围内对输电铁塔进行变形监测,提高了输电铁塔变形监测的可靠性。
在一个实施例中,在覆盖范围内获取输电铁塔表面的散射体对应的坐标,包括:若散射体的光强大于预设光强阈值,获取散射体对应的坐标,作为散射体的像元坐标;将多张不同拍摄时间的SAR影像进行非相干叠加,得到多张不同拍摄时间的SAR影像对应的平均光强影像;对平均光强影像进行过采样处理,得到更新后的平均光强影像;根据更新后的平均光强影像中每个像素的坐标、每个像素的光强以及二维辛格函数中像素对应的坐标的数值,得到更新后的平均光强影像与二维辛格函数的相关系数;获取相关系数的最大值对应的坐标,作为平均光强影像中散射体的子像元坐标;基于散射体的像元坐标、子像元坐标以及预设过采样系数,得到输电铁塔中的散射体对应的坐标。
本实施例中,散射体可以是上述输电铁塔表面的散射体,每个散射体可以有一定的光强度,终端102可以在上述输电铁塔的覆盖范围内获取输电铁塔表面的散射体的坐标。终端102可以获取上述散射体的像元坐标,终端102可以获取大于预设光强阈值的散射体在上述SAR影像中对应的坐标,作为散射体的像元坐标,具体地,输电铁塔在雷达影像中表现为明显的亮斑,雷达回波信号较为稳定,可作为高信杂比的散射体。对SAR强度影像进行分析,采用强度阈值方法实现铁塔表面散射体的识别及像元坐标提取,铁塔表面散射体的像元坐标可标识为(XPix,YPix),其中XPix为距离向像元坐标,YPix为方位向像元坐标。
终端102还可以获取上述散射体的子像元坐标,终端102可以将多张不同拍摄时间的SAR影像进行非相干叠加,例如将多张SAR影像中的输电铁塔进行非相干叠加,从而得到多张不同拍摄时间的SAR影像对应的平均光强影像,其中非相干叠加可以是不减弱叠加的振动强度的叠加方式。终端102可以对上述平均光强影像进行过采样处理,即以远远高于信号带宽两倍或其最高频率对其进行采样的过程,得到更新后的平均光强影像,并根据更新后的平均光强影像中每个像素的坐标、每个像素的光强以及二维辛格函数中该像素对应的坐标的数值,得到更新后平均光强影像和二维辛格函数的相关系数,并且可以将相关系数最大值对应的坐标,作为上述平均光强影像中散射体的子像元坐标。
具体地,终端102可以将收集的多景SAR影像,即上述多张不同拍摄时间的SAR影像,进行非相干叠加,获取平均光强影像,从而减弱噪声的影响,并对平均光强影像做过采样处理,求取其与二维sinc函数,即二维辛格函数的相关系数,以相关系数峰值对应的坐标作为散射体的子像元坐标值,子像元的坐标值可以表示为(XSub,YSub),其中XSub为距离向子像元坐标,YSub为方位向子像元坐标。相关系数计算公式如下所示:
其中,ρ为相关系数,Amn为所匹配强度影像块中坐标为(m,n)像素的强度值,Bmn为二维sinc函数中(m,n)像素对应的数值。通过移动二维sinc函数窗口与平均强度影像进行匹配,并计算相关系数,相关系数最大值对应的坐标即为子像元的坐标。
终端102在得到散射体的像元坐标和子像元坐标后,可以基于像元坐标以及子像元坐标,得到输电铁塔表面散射体在SAR影像中的精确坐标,具体地,该精确坐标可以表示为(X,Y),其计算如下式所示:
X=XPix+XSub/OX
Y=YPix+YSub/OY
其中,OX为设定的距离向过采样系数,OY为设定的方位向过采样系数,过采样系数通常取值为2~4。
通过本实施例,终端102可以利用像元坐标和子像元坐标得到散射体在SAR影像中的精确坐标,从而可以根据散射体对输电铁塔的变形进行监测,提高了输电铁塔变形监测的可靠性。
在一个实施例中,根据在SAR影像中的散射体的坐标以及输电铁塔的底部中心坐标,得到散射体与输电铁塔底部中心的高度差,包括:获取散射体与输电铁塔底部中心的距离向坐标的差值;基于差值以及SAR影像的距离向分辨率,得到散射体与输电铁塔底部中心的高度差。
本实施例中,终端102根据SAR影像中的散射体的坐标,以及输电铁塔的底部中心坐标,得到散射体和输电铁塔底部中心的高度差。终端102可以获取上述散射体和输电铁塔底部中心的距离向坐标的差值,并基于上述差值和SAR影像的距离向分辨率,得到散射体和输电铁塔底部中心的高度差。距离向可以是沿着雷达卫星发射的雷达波的方向。具体地,终端102可以根据铁塔在SAR影像中的投影关系,得到铁塔表面散射体相对于铁塔底部的高度差,终端102可以将散射体的坐标设为(X,Y),铁塔底部中心的坐标可以是(XRef,YRef),其中XRef为距离向坐标,YRef为方位向坐标,铁塔表面散射体到铁塔底部中心像元的距离向坐标差值ΔX可以表示为:
ΔX=X-XRef
铁塔表面散射体相对铁塔底部的高度差Δh可表示为:
Δh=ΔX*tanθ*RX
其中,θ为雷达卫星入射角,RX为SAR影像距离向分辨率。
通过本实施例,终端102可以利用散射体的坐标以及铁塔底部中心坐标得到两者之间的高度差,从而可以基于高度差得到高程误差及其相位,实现了提高输电铁塔变形监测可靠性的效果。
在一个实施例中,如图3所示,图3为另一个实施例中输电铁塔变形监测方法的流程示意图。终端102可以首先计算铁塔在SAR影像中的覆盖范围,具体地,可以利用铁塔底部中心对应的SAR影像坐标值、上述铁塔在SAR影像中的距离向长度以及方位向长度得到;终端102还可以对输电铁塔表面散射体影像坐标进行精确提取,具体地,可以识别输电铁塔表面散射体并提取散射体的像元坐标以及子像元坐标,并根据像元坐标以及子像元坐标得到散射体的精确坐标;终端102还可以获取铁塔表面散射体相对铁塔底部的高度差,具体地,可以额计算散射体到铁塔底部中心像元的距离向坐标差值,基于该差值计算铁塔表面散射体相对铁塔底部的高度差。终端102还可以进行干涉相位精化,即去除误差。具体地,终端102可以将输电铁塔表面散射体的高程误差进行估计并去除,还可以将上述散射体热膨胀效应对应的相位进行估计和去除,从而基于去除误差后的相位,对输电铁塔的变形进行监测。
通过本实施例,终端102可以基于SAR影像中,输电铁塔表面的散射体的坐标及其干涉相位,对输电铁塔的变形情况进行监测,从而可以提高输电铁塔变形监测的可靠性。
为更好理解上述方法,以下详细阐述本发明输电铁塔变形监测方法的应用实例。
在一个应用实施例中,如图4所示,图4为一个实施例中输电铁塔的SAR影像示意图。终端102可以选取某城市500kV输电线路铁塔为试验对象,采用二景高分辨率X波段雷达卫星影像,空间分辨率为距离向3m×方位向3m。影像拍摄时间分别为20150108和20160528,时间间隔为506天,空间垂直基线为-176m。下图为输电线路铁塔的后向散射强度图,其在SAR影像中表现为明显的亮斑,强度值高于周边地物。
终端102可以利用强度阈值法与二维sinc函数,即二维辛格函数,通过匹配实现铁塔表面散射体影像坐标的精确提取,基于散射体坐标提取对应点位的干涉相位值,如图5所示,图5为一个实施例中散射体干涉相位的示意图。横轴可以是距离向坐标轴,数轴可以是方位向坐标轴,图中每个点可以是每个散射体对应的坐标,每个点包含的颜色可以是该点对应的散射体的干涉相位值,终端102可以根据颜色深浅程度确定该散射体的变形程度,具体地,如图中的干涉相位值标准所示,相位越低则颜色越深,即颜色越深则变形程度越严重,然而图5中的干涉相位为高程误差相位、热膨胀效应所引起相位、形变相位、热噪声的叠加,其相位值不能够准确反映铁塔的变形情况。
因此,终端102可以基于铁塔表面散射体与铁塔底部中心所对应像素点的距离向坐标差值,计算铁塔表面散射体的高程误差及对应的高程误差相位,并利用两景影像拍摄时间的外界环境温度差,实现铁塔表面散射体热膨胀所对应相位的去除。最终获取的铁塔表面散射体的形变相位可以如图6所示,图6为一个实施例中去除误差后的散射体干涉相位的示意图。终端102可以根据图中每个散射体对应的点的颜色深浅,即相位大小,确定输电铁塔的变形情况。
通过本实施例,终端102可以利用SAR影像中输电铁塔表面的散射体的坐标及其相位,对输电铁塔的变形情况进行监测,从而提高了输电铁塔变形监测的可靠性。
应该理解的是,虽然图2至图3的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图2至图3中的至少一部分步骤可以包括多个步骤或者多个阶段,这些步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。
在一个实施例中,如图7所示,提供了一种输电铁塔变形监测装置,包括:坐标获取模块500、相位获取模块502、误差获取模块504和监测模块506,其中:
坐标获取模块500,用于获取输电铁塔的SAR影像;SAR影像中包括至少两个对应不同拍摄时间的输电铁塔图像;根据输电铁塔在SAR影像中的坐标值,得到输电铁塔在SAR影像中的覆盖范围;在覆盖范围内获取输电铁塔表面的散射体对应的坐标。
相位获取模块502,用于根据散射体的坐标,从预设相位表中获取散射体的坐标对应的干涉相位值;预设相位表包括坐标与相位的对应关系。
误差获取模块504,用于根据在SAR影像中的散射体的坐标以及输电铁塔的底部中心坐标,得到散射体与输电铁塔底部中心的高度差;基于高度差、SAR影像对应的雷达卫星的电磁波波长、雷达卫星的垂直基线长度以及雷达卫星与输电铁塔的距离,得到散射体的高度误差以及高度误差相位;基于两个输电铁塔图像对应的拍摄时间的环境温度差、输电铁塔的材质对应的热膨胀系数,得到散射体的热膨胀相位。
监测模块506,用于将高度误差相位以及热膨胀相位从散射体的干涉相位值中去除,得到更新后的散射体的干涉相位值;根据更新后的散射体的干涉相位值,对输电铁塔进行变形监测。
在一个实施例中,上述坐标获取模块500,具体用于获取SAR影像中电铁塔的底部中心坐标值;获取SAR影像中输电铁塔的距离向长度;获取SAR影像中输电铁塔的方位向长度;基于底部中心坐标值、距离向长度以及方位向长度,得到输电铁塔在SAR影像中的覆盖范围。
在一个实施例中,上述装置还包括:映射模块,用于基于雷达卫星的轨道信息以及SAR影像信息,生成SAR影像与地面的坐标映射关系,以在坐标映射关系中,获取输电铁塔在SAR影像中的坐标值。
在一个实施例中,上述坐标获取模块500,具体用于获取输电铁塔的实际高度;获取雷达卫星的入射角;根据实际高度、入射角以及SAR影像的距离向分辨率,获取SAR影像中输电铁塔的距离向长度。
在一个实施例中,上述坐标获取模块500,具体用于获取输电铁塔的实际宽度;根据实际宽度以及SAR影像的方位向分辨率,获取SAR影像中输电铁塔的方位向长度。
在一个实施例中,上述坐标获取模块500,具体用于若散射体的光强大于预设光强阈值,获取散射体对应的坐标,作为散射体的像元坐标;将多张不同拍摄时间的SAR影像进行非相干叠加,得到多张不同拍摄时间的SAR影像对应的平均光强影像;对平均光强影像进行过采样处理,得到更新后的平均光强影像;根据更新后的平均光强影像中每个像素的坐标、每个像素的光强以及二维辛格函数中像素对应的坐标的数值,得到更新后的平均光强影像与二维辛格函数的相关系数;获取相关系数的最大值对应的坐标,作为平均光强影像中散射体的子像元坐标;基于散射体的像元坐标、子像元坐标以及预设过采样系数,得到输电铁塔中的散射体对应的坐标。
在一个实施例中,上述误差获取模块504,具体用于获取散射体与输电铁塔底部中心的距离向坐标的差值;基于差值以及SAR影像的距离向分辨率,得到散射体与输电铁塔底部中心的高度差。
关于输电铁塔变形监测装置的具体限定可以参见上文中对于输电铁塔变形监测方法的限定,在此不再赘述。上述输电铁塔变形监测装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,该计算机设备可以是终端,其内部结构图可以如图8所示。该计算机设备包括通过系统总线连接的处理器、存储器、通信接口、显示屏和输入装置。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的通信接口用于与外部的终端进行有线或无线方式的通信,无线方式可通过WIFI、运营商网络、NFC(近场通信)或其他技术实现。该计算机程序被处理器执行时以实现一种输电铁塔变形监测方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏,该计算机设备的输入装置可以是显示屏上覆盖的触摸层,也可以是计算机设备外壳上设置的按键、轨迹球或触控板,还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解,图8中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
在一个实施例中,提供了一种计算机设备,包括存储器和处理器,存储器中存储有计算机程序,该处理器执行计算机程序时实现如上述的输电铁塔变形监测方法。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现如上述的输电铁塔变形监测方法。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(Read-Only Memory,ROM)、磁带、软盘、闪存或光存储器等。易失性存储器可包括随机存取存储器(Random Access Memory,RAM)或外部高速缓冲存储器。作为说明而非局限,RAM可以是多种形式,比如静态随机存取存储器(Static Random Access Memory,SRAM)或动态随机存取存储器(Dynamic Random Access Memory,DRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种输电铁塔变形监测方法,其特征在于,所述方法包括:
获取所述输电铁塔的SAR影像;所述SAR影像中包括至少两个对应不同拍摄时间的输电铁塔图像;
根据所述输电铁塔在SAR影像中的坐标值,得到所述输电铁塔在SAR影像中的覆盖范围;
在所述覆盖范围内获取所述输电铁塔表面的散射体对应的坐标;
根据所述散射体的坐标,从预设相位表中获取所述散射体的坐标对应的干涉相位值;所述预设相位表包括所述坐标与相位的对应关系;
根据在所述SAR影像中的所述散射体的坐标以及所述输电铁塔的底部中心坐标,得到所述散射体与所述输电铁塔底部中心的高度差;
基于所述高度差、所述SAR影像对应的雷达卫星的电磁波波长、所述雷达卫星的垂直基线长度以及所述雷达卫星与所述输电铁塔的距离,得到所述散射体的高度误差以及高度误差相位;
基于两个所述输电铁塔图像对应的拍摄时间的环境温度差、所述输电铁塔的材质对应的热膨胀系数,得到所述散射体的热膨胀相位;
将所述高度误差相位以及所述热膨胀相位从所述散射体的干涉相位值中去除,得到更新后的散射体的干涉相位值;根据更新后的散射体的干涉相位值,对所述输电铁塔进行变形监测。
2.根据权利要求1所述的方法,其特征在于,所述根据所述输电铁塔在SAR影像中的坐标值,得到所述输电铁塔在SAR影像中的覆盖范围,包括:
获取所述SAR影像中所述输电铁塔的底部中心坐标值;
获取所述SAR影像中所述输电铁塔的距离向长度;
获取所述SAR影像中所述输电铁塔的方位向长度;
基于所述底部中心坐标值、所述距离向长度以及所述方位向长度,得到所述输电铁塔在SAR影像中的覆盖范围。
3.根据权利要求1所述的方法,其特征在于,所述根据所述输电铁塔在SAR影像中的坐标值,得到所述输电铁塔在SAR影像中的覆盖范围之前,还包括:
基于所述雷达卫星的轨道信息以及所述SAR影像信息,生成所述SAR影像与地面的坐标映射关系,以在所述坐标映射关系中,获取所述输电铁塔在SAR影像中的坐标值。
4.根据权利要求2所述的方法,其特征在于,所述获取所述SAR影像中所述输电铁塔的距离向长度,包括:
获取所述输电铁塔的实际高度;
获取所述雷达卫星的入射角;
根据所述实际高度、所述入射角以及所述SAR影像的距离向分辨率,获取所述SAR影像中所述输电铁塔的距离向长度。
5.根据权利要求2所述的方法,其特征在于,所述获取所述SAR影像中所述输电铁塔的方位向长度,包括:
获取所述输电铁塔的实际宽度;
根据所述实际宽度以及所述SAR影像的方位向分辨率,获取所述SAR影像中所述输电铁塔的方位向长度。
6.根据权利要求1所述的方法,其特征在于,所述在所述覆盖范围内获取所述输电铁塔表面的散射体对应的坐标,包括:
若所述散射体的光强大于预设光强阈值,获取所述散射体对应的坐标,作为散射体的像元坐标;
将多张不同拍摄时间的SAR影像进行非相干叠加,得到所述多张不同拍摄时间的SAR影像对应的平均光强影像;
对所述平均光强影像进行过采样处理,得到更新后的平均光强影像;根据更新后的平均光强影像中每个像素的坐标、每个像素的光强以及二维辛格函数中所述像素对应的坐标的数值,得到所述更新后的平均光强影像与所述二维辛格函数的相关系数;
获取所述相关系数的最大值对应的坐标,作为所述平均光强影像中散射体的子像元坐标;
基于所述散射体的像元坐标、所述子像元坐标以及预设过采样系数,得到所述输电铁塔中的散射体对应的坐标。
7.根据权利要求1所述的方法,其特征在于,所述根据在所述SAR影像中的所述散射体的坐标以及所述输电铁塔的底部中心坐标,得到所述散射体与所述输电铁塔底部中心的高度差,包括:
获取所述散射体与所述输电铁塔底部中心的距离向坐标的差值;
基于所述差值以及所述SAR影像的距离向分辨率,得到所述散射体与所述输电铁塔底部中心的高度差。
8.一种输电铁塔变形监测装置,其特征在于,所述装置包括:
坐标获取模块,用于获取所述输电铁塔的SAR影像;所述SAR影像中包括至少两个对应不同拍摄时间的输电铁塔图像;根据所述输电铁塔在SAR影像中的坐标值,得到所述输电铁塔在SAR影像中的覆盖范围;在所述覆盖范围内获取所述输电铁塔表面的散射体对应的坐标;
相位获取模块,用于根据所述散射体的坐标,从预设相位表中获取所述散射体的坐标对应的干涉相位值;所述预设相位表包括所述坐标与相位的对应关系;
误差获取模块,用于根据在所述SAR影像中的所述散射体的坐标以及所述输电铁塔的底部中心坐标,得到所述散射体与所述输电铁塔底部中心的高度差;基于所述高度差、所述SAR影像对应的雷达卫星的电磁波波长、所述雷达卫星的垂直基线长度以及所述雷达卫星与所述输电铁塔的距离,得到所述散射体的高度误差以及高度误差相位;基于两个所述输电铁塔图像对应的拍摄时间的环境温度差、所述输电铁塔的材质对应的热膨胀系数,得到所述散射体的热膨胀相位;
监测模块,用于将所述高度误差相位以及所述热膨胀相位从所述散射体的干涉相位值中去除,得到更新后的散射体的干涉相位值;根据更新后的散射体的干涉相位值,对所述输电铁塔进行变形监测。
9.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述的方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至7中任一项所述的方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010541550.XA CN111487622A (zh) | 2020-06-15 | 2020-06-15 | 输电铁塔变形监测方法、装置、计算机设备和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010541550.XA CN111487622A (zh) | 2020-06-15 | 2020-06-15 | 输电铁塔变形监测方法、装置、计算机设备和存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111487622A true CN111487622A (zh) | 2020-08-04 |
Family
ID=71813480
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010541550.XA Pending CN111487622A (zh) | 2020-06-15 | 2020-06-15 | 输电铁塔变形监测方法、装置、计算机设备和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111487622A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112197690A (zh) * | 2020-09-30 | 2021-01-08 | 深圳市水务规划设计院股份有限公司 | 基于角反射器的土石坝变形监测方法、系统及存储介质 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103323848A (zh) * | 2013-06-19 | 2013-09-25 | 中国测绘科学研究院 | 一种提取地面人工建/构筑物高度的方法及装置 |
CN105954747A (zh) * | 2016-06-20 | 2016-09-21 | 中国电力工程顾问集团中南电力设计院有限公司 | 基于电网不良地质体三维形变监测的塔基稳定性分析方法 |
CN105974411A (zh) * | 2016-04-28 | 2016-09-28 | 电子科技大学 | 基于dinsar的高压输电铁塔顶端倾斜位移监测方法 |
CN106772377A (zh) * | 2017-01-18 | 2017-05-31 | 深圳市路桥建设集团有限公司 | 一种基于InSAR的建筑物变形监测方法 |
CN108549080A (zh) * | 2018-02-28 | 2018-09-18 | 中国电力科学研究院有限公司 | 一种输电杆塔位置提取方法及系统 |
CN108983229A (zh) * | 2018-05-03 | 2018-12-11 | 电子科技大学 | 基于sar层析技术的高压输电铁塔高度及形变提取方法 |
CN109583284A (zh) * | 2017-09-29 | 2019-04-05 | 中国科学院遥感与数字地球研究所 | 基于高分辨率sar图像的城市高层建筑物高度提取方法及装置 |
CN110456346A (zh) * | 2019-06-28 | 2019-11-15 | 深圳市水务规划设计院股份有限公司 | 一种基于InSAR技术的输电铁塔倾斜监测方法 |
-
2020
- 2020-06-15 CN CN202010541550.XA patent/CN111487622A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103323848A (zh) * | 2013-06-19 | 2013-09-25 | 中国测绘科学研究院 | 一种提取地面人工建/构筑物高度的方法及装置 |
CN105974411A (zh) * | 2016-04-28 | 2016-09-28 | 电子科技大学 | 基于dinsar的高压输电铁塔顶端倾斜位移监测方法 |
CN105954747A (zh) * | 2016-06-20 | 2016-09-21 | 中国电力工程顾问集团中南电力设计院有限公司 | 基于电网不良地质体三维形变监测的塔基稳定性分析方法 |
CN106772377A (zh) * | 2017-01-18 | 2017-05-31 | 深圳市路桥建设集团有限公司 | 一种基于InSAR的建筑物变形监测方法 |
CN109583284A (zh) * | 2017-09-29 | 2019-04-05 | 中国科学院遥感与数字地球研究所 | 基于高分辨率sar图像的城市高层建筑物高度提取方法及装置 |
CN108549080A (zh) * | 2018-02-28 | 2018-09-18 | 中国电力科学研究院有限公司 | 一种输电杆塔位置提取方法及系统 |
CN108983229A (zh) * | 2018-05-03 | 2018-12-11 | 电子科技大学 | 基于sar层析技术的高压输电铁塔高度及形变提取方法 |
CN110456346A (zh) * | 2019-06-28 | 2019-11-15 | 深圳市水务规划设计院股份有限公司 | 一种基于InSAR技术的输电铁塔倾斜监测方法 |
Non-Patent Citations (2)
Title |
---|
刘晓龙等: "时序PSInSAR研究建筑物高度信息提取", 《武汉大学学报(信息科学版)》 * |
文江平等: "基于高分辨SAR图像区域入射角反演及高度提取方法", 《计算机应用》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112197690A (zh) * | 2020-09-30 | 2021-01-08 | 深圳市水务规划设计院股份有限公司 | 基于角反射器的土石坝变形监测方法、系统及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8154435B2 (en) | Stability monitoring using synthetic aperture radar | |
Gance et al. | Target Detection and Tracking of moving objects for characterizing landslide displacements from time-lapse terrestrial optical images | |
Iervolino et al. | Flooding water depth estimation with high-resolution SAR | |
CN107917880B (zh) | 一种基于地基云图的云底高度反演方法 | |
Morrison et al. | Laboratory measurement of the DInSAR response to spatiotemporal variations in soil moisture | |
JP2021056008A (ja) | 地すべり領域検出装置及びプログラム | |
CN112904337B (zh) | 一种基于Offset Tracking技术的边坡形变时序监测方法 | |
CN112558066A (zh) | 双极化sar图像系统 | |
KR102151362B1 (ko) | 극좌표변환을 이용한 항공기기반 영상복원장치 및 이를 이용한 영상복원방법 | |
Chai et al. | Deformation monitoring using ground-based differential SAR tomography | |
Pieraccini et al. | A GBSAR operating in monostatic and bistatic modalities for retrieving the displacement vector | |
CN111487622A (zh) | 输电铁塔变形监测方法、装置、计算机设备和存储介质 | |
CN111487621A (zh) | 一种基于雷达图像的海表流场反演方法及电子设备 | |
Lingua et al. | Remote monitoring of a landslide using an integration of GB-INSAR and LIDAR techniques | |
Macfarlane et al. | Topographic and thermal mapping of volcanic terrain using the AVTIS ground-based 94-GHz dual-mode radar/radiometric imager | |
CN110244301A (zh) | 一种基于sar图像的输电走廊树高提取方法 | |
Zhang et al. | Application of the method for prediction of the failure location and time based on monitoring of a slope using synthetic aperture radar | |
Toutin et al. | Fusion of Radarsat-2 polarimetric images for improved stereo-radargrammetric DEM | |
CN113030968B (zh) | 基于csar模式提取dem的方法、装置及存储介质 | |
CN115932823A (zh) | 基于异源区域特征匹配的飞行器对地面目标定位方法 | |
CN115980746A (zh) | 一种可行驶区域的检测方法、装置、车辆及存储介质 | |
André et al. | Spatially variant incoherence trimming for improved bistatic SAR CCD | |
Pujianiki et al. | A New Simple Procedure for Extracting Coastline from SAR Image Based on Low Pass Filter and Edge Detection Algorithm | |
Xu et al. | Theoretical modeling of multi-frequency tomography radar observations of snow stratigraphy | |
US20200166626A1 (en) | Information processing device, information processing method, and storage medium having program stored thereon |
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 | ||
AD01 | Patent right deemed abandoned | ||
AD01 | Patent right deemed abandoned |
Effective date of abandoning: 20220802 |