CN106483511B - 一种输电线路在星载合成孔径雷达成像中的电磁仿真方法 - Google Patents
一种输电线路在星载合成孔径雷达成像中的电磁仿真方法 Download PDFInfo
- Publication number
- CN106483511B CN106483511B CN201510543429.XA CN201510543429A CN106483511B CN 106483511 B CN106483511 B CN 106483511B CN 201510543429 A CN201510543429 A CN 201510543429A CN 106483511 B CN106483511 B CN 106483511B
- Authority
- CN
- China
- Prior art keywords
- conducting wire
- satellite
- axis
- vector
- transmission line
- 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
- 230000005540 biological transmission Effects 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000005611 electricity Effects 0.000 title claims abstract description 36
- 238000004088 simulation Methods 0.000 title claims abstract description 22
- 238000003384 imaging method Methods 0.000 title claims abstract description 13
- 230000005855 radiation Effects 0.000 claims abstract description 34
- 239000013598 vector Substances 0.000 claims description 72
- 238000005070 sampling Methods 0.000 claims description 34
- 239000004020 conductor Substances 0.000 claims description 16
- 229910000831 Steel Inorganic materials 0.000 claims description 15
- 239000010959 steel Substances 0.000 claims description 15
- 238000005286 illumination Methods 0.000 claims description 3
- 238000010304 firing Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 5
- 238000012423 maintenance Methods 0.000 abstract description 4
- 230000006872 improvement Effects 0.000 abstract description 3
- 238000007689 inspection Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000012876 topography Methods 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 208000032370 Secondary transmission Diseases 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000003595 mist Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
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
- 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/40—Means for monitoring or calibrating
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种输电线路在星载合成孔径雷达成像中的电磁仿真方法,该方法通过对卫星轨道多项式拟合求解和输电线路导线三维建模解算,高精度的解出任意时刻的卫星轨道位置和其辐射中心对应的导线位置;通过SAR几何构像和散射模型建模,计算导线点目标对应的辐射能量值,最终通过能量辐射叠加确定最终的导线目标中心点位置辐射能量值。本发明的仿真方法,各个步骤采用的模型均为成熟且严谨的模型,有较好的鲁棒性和正确性,其计算过程准确性及可靠性高,同时便于后期的改进和维护,方便应用在后续的具体应用中;且能够识别输电线路导线是否存在松弛、断线等相关损坏情况,保证了输电线路的可行性运行。
Description
技术领域
本发明涉及星载合成孔径雷达遥感图像处理领域,具体涉及一种输电线路在星载合成孔径雷达成像中的电磁仿真方法。
背景技术
为了掌握线路的运行状况和沿线的环境情况,及时发现是否存在缺陷和安全隐患,需要对线路进行定期巡检。尤其是在气候变化(大雾、冰雪、强风、暴雨、雷暴等)、自然灾害(地震、洪水、泥石流、森林火灾等)、外力破坏、异常运行和其他特殊情况时,需要及时进行特殊巡检,发现和解决线路的异常情况。
我国输电线路常绵延数百甚至上千公里,多数位于人口稀少地区,跨越高山峻岭,沿线地形复杂,进行检测及巡检难度大且成本高。人工巡检、机器人巡检、直升机巡检和基于多种传感器的在线检测装置等巡检监测技术,在极端天气条件和大范围自然灾害条件下,难以快速、大范围、低成本和高安全系数进行巡检。
卫星平台的快速发展,使得基于光学遥感、多光谱遥感、合成孔径雷达遥感等遥感技术在电力系统的应用变成可能。基于合成孔径雷达的卫星遥感技术具有全天时、全天候、大范围探测能力,可以在大范围灾害下对输电线路进行广域监测,解决各种复杂气象条件和复杂地形条件的巡检检测难题。
合成孔径雷达观测技术,是主动微波遥感技术,可以不受时间限制,可以穿云透雾,能够在各种极端环境下获得高质量图像,但雷达图像与光学图像的差异,不同地物特征的辐射特性各不相同,铁塔及电线的辐射特性更是难以定量分析,造成了雷达图像解译的困难,如何基于星载合成孔径雷达图像提取有效信息,进行输电线路巡检和检测,分析和发现输电线路的异常区域,具有非常重要的意义。
目前,国内外在相关领域的研究工作主要集中在雷达图像中铁塔位置检测和状态识 别,通过不同时间雷达图像中铁塔相关信息的变化,分析铁塔的工作状态,这在一定程度上可以解决目前输电线路巡检和监测的需求,但是输电线路地形复杂,存在遮挡和背景干扰,有些区域无法全面有效的监测,而且对于输电线路导线是否存在松弛、断线等相关损坏情况,是不能够识别的。
发明内容
有鉴于此,本发明提供的一种输电线路在星载合成孔径雷达成像中的电磁仿真方法,该方法基于输电导线空间位置模型和雷达后向散射模型仿真导线的高分辨率SAR图像成像特性,解决了雷达影像中输电线路导线成像状态分析的问题,为分析输电线路的工作运行状态和损坏情况提供了理论模型和仿真算法。
本发明的目的是通过以下技术方案实现的:
一种输电线路在星载合成孔径雷达成像中的电磁仿真方法,所述方法包括如下步骤:
步骤1.根据卫星轨道离散采样状态矢量数据,计算得到星载合成孔径雷达的卫星轨道任意时刻的轨道状态矢量;
步骤2.构建输电线路中的导线的斜抛物线模型;
步骤3.根据所述导线的斜抛物线模型及卫星轨道的位置坐标,求得所述导线的入射角;
步骤4.根据所述导线的入射角,构建并计算雷达辐射区域导线计算模型,得到对应的导线目标中心点位置辐射能量值;
步骤5.根据计算结果,生成仿真区域背景,并叠加和仿真导线。
优选的,所述步骤1,包括:
1-1.从卫星平台提供的卫星雷达数据头文件中读取卫星轨道离散采样状态矢量数据;所述卫星轨道离散采样状态矢量数据包括离散采样卫星位置矢量RS(ti)、离散采样卫星速度矢量VS(ti)及采样对应的时刻ti(i=1,2,…,n);
RS(ti)=(XS(ti),YS(ti),ZS(ti))T;
VS(ti)=(VSX(ti),VSY(ti),VSZ(ti))T;
式中,XS(ti),YS(ti),ZS(ti)为卫星在x轴、y轴和Z轴上的离散采样位置矢量;
VSX(ti),VSY(ti),VSZ(ti)分别为卫星在x轴、y轴和Z轴上的离散采样速度矢量;
1-2.构建离散采样卫星位置矢量RS(ti)的三次多项式拟合轨道方程:
构建离散采样卫星速度矢量VS(ti)的三次多项式拟合轨道方程:
式中,(a0,a1,a2,a3,b0,b1,b2,b3,c0,c1,c2,c3)为三次多项式拟合轨道方程的拟合参数;
1-3.根据RS(ti)及VS(ti),采用最小估计二乘法求解所述三次多项式拟合轨道方程,得到拟合参数(a0,a1,a2,a3,b0,b1,b2,b3,c0,c1,c2,c3)的值;
1-4.根据拟合参数(a0,a1,a2,a3,b0,b1,b2,b3,c0,c1,c2,c3)的值,求得卫星轨道任意时刻星载合成孔径雷达的卫星轨道任意时刻t的轨道位置矢量RS(t)和速度矢量VS(t),即求得任意时刻的轨道位置坐标(xs,ys,zs)和速度坐标(Vsx,Vsy,Vsz);
其中,RS(t)的拟合轨道方程为:
VS(t)的拟合轨道方程为:
式中,XS(t),YS(t),ZS(t)分别为卫星在x轴、y轴和Z轴上的任意时刻t的位置矢量;VSX(t),VSY(t),VSZ(t)分别为卫星在x轴、y轴和Z轴上的任意时刻t的速度矢量。
优选的,所述步骤2,包括:
2-1.构建所述输电线路中的导线进行斜抛物线模型:
当所述输电线路不平行于X轴时:
当所述输电线路平行于X轴时:
式中,m、n和R为所述导线所在重力平面方程系数,r1,r2,r3为导线在以z轴为中心的旋转抛物面方程系数,其中g和ρ0为导线相关系数参数;x,y,z为导线铁塔空间坐标。
2-2.若所述导线不平行于x坐标轴,则r1,r2,r3的值为:
r2=mr1+n;
r3=zt1-A[(xt1-r1)2+(yt1-r2)2];
其中,n=yt1-mxt1;xt1,yt1,zt1和xt2,yt2,zt2分别为所述导线两端铁塔的空间坐标;
2-3.若所述导线平行于x坐标轴,则r1,r2,r3的值为:
r1=R;
r3=zt1-A(yt1-r2)2。
优选的,所述中的求解方法包括:
根据已知所述导线的两端铁塔的位置坐标(xt1,yt1,zt1)和(xt2,yt2,zt2),求出高差角的余弦值:
优选的,所述步骤3,包括:
3-1.任意时刻的卫星轨道位置坐标(xs,ys,zs)和目标点位置坐标(xp,yp,zp),基于距离公式计算SAR数据中对应的斜距r:
r2=(xs-xp)2+(ys-yp)2+(zs-zp)2;
3-2.对于所述导线上的任一点,对所述导线的斜抛物线模型方程求偏导,获得导线的方向向量
若所述导线平行于X轴,则所述方向向量为(0,1,2A(y-r2));
若所述导线不平行于X轴,
则所述方向向量为(1,m,2A(x-r1)+2mA(mx+n-r2));
3-3.利用入射向量和导线方向向量基于空间几何夹角公式计算所述导线的入射角θ:
式中,雷达电磁波的入射向量
优选的,所述步骤4,包括:
4-1.根据所述导线的入射角θ,构建雷达辐射区域点目标散射模型:
e=σ0(cosθ)2;
式中,e为雷达辐射区域点目标散射模型,σ0为后向散射系数,θ为雷达照射目标区域入射角;
4-2.确定仿真数据星载SAR相关参数;
4-3.构建雷达辐射区域导线计算模型;
对于给定时刻t卫星平台位置(XS(t),YS(t),ZS(t))T及正侧视中心导线位置坐标(xp,yp,zp)T,确定所述雷达辐射带宽范围内的导线坐标为:
[(xp-,yp-,zp-)T,(xp+,yp+,zp+)T];
4-4.计算并叠加雷达辐射带宽区域所有导线上的点目标的辐射能量值,生成对应的导线目标中心点位置辐射能量值。
优选的,所述相关参数包括卫星飞行高度约为、辐射波长、波束宽度方位向、距离向、在方位向辐射宽度及辐射分辨率。
优选的,所述步骤5,包括:
5-1.采用光学灰度图像作为导线区域背景,对图像叠加高斯乘性噪声模拟SAR图像;
5-2.基于斜距确定导线目标在背景图像中的坐标位置,合并背景区域和导线目标仿真能量值,生成仿真图像。
从上述的技术方案可以看出,本发明提供了一种输电线路在星载合成孔径雷达成像中的电磁仿真方法,该方法通过对卫星轨道多项式拟合求解和输电线路导线三维建模解算,高精度的解出任意时刻的卫星轨道位置和其辐射中心对应的导线位置;通过SAR几何构像和散射模型建模,计算导线点目标对应的辐射能量值,最终通过能量辐射叠加确定最终的导线目标中心点位置辐射能量值。本发明的仿真方法,各个步骤采用的模型均为成熟且严谨的模型,有较好的鲁棒性和正确性,其计算过程准确性及可靠性高,同时便于后期的改进和维护,方便应用在后续的具体应用中;且能够识别输电线路导线是否存在松弛、断线等相关损坏情况,保证了输电线路的可行性运行。
与最接近的现有技术比,本发明提供的技术方案具有以下优异效果:
1、本发明所提供的技术方案中,通过对卫星轨道多项式拟合求解和输电线路导线三维建模解算,高精度的解出任意时刻的卫星轨道位置和其辐射中心对应的导线位置;通过SAR几何构像和散射模型建模,计算导线点目标对应的辐射能量值,最终通过能量辐射叠加确定最终的导线目标中心点位置辐射能量值。本发明的仿真方法,各个步骤采用的模型均为成熟且严谨的模型,有较好的鲁棒性和正确性,其计算过程准确性及可靠性高,同时便于后期的改进和维护,方便应用在后续的具体应用中;且能够识别输电线路导线是否存在松弛、断线等相关损坏情况,保证了输电线路的可行性运行。
2、本发明所提供的技术方案,相较于大多数仿真和应用方法都是针对铁塔等输电线路目标进行相关分析的思路,本文方法与以往输电线路仿真应用方法有本质不同,第一次对输电线路导线进行分析、建模和仿真,提出了一整套的比较严密的理论模型和仿真方法,为基于导线对输电线路进行分析和运行状态景象评估提出了一种新的思路和模型
3、本发明所提供的技术方案,输电线路导线区域的辐射模型也是严格基于雷达电磁辐射模型进行解算,因此,保证了计算过程的准确性和可靠性。
4、本发明所提供的技术方案,模块化分解计算,所以便于后期的改进和维护,方便应用在后续的具体应用中。
5、本发明提供的技术方案,应用广泛,具有显著的社会效益和经济效益。
附图说明
图1是本发明的一种输电线路在星载合成孔径雷达成像中的电磁仿真方法的流程图;
图2是本发明的电磁仿真方法中步骤1的流程图;
图3是本发明的电磁仿真方法中步骤2的流程图;
图4是本发明的电磁仿真方法中步骤3的流程图;
图5是本发明的电磁仿真方法中步骤4的流程图;
图6是本发明的电磁仿真方法中步骤5的流程图;
图7是本发明的仿真方法的应用例的整体流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整 地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供一种输电线路在星载合成孔径雷达成像中的电磁仿真方法,包括如下步骤:
步骤1.根据卫星轨道离散采样状态矢量数据,计算得到星载合成孔径雷达的卫星轨道任意时刻的轨道状态矢量;
步骤2.构建输电线路中的导线的斜抛物线模型;
步骤3.根据导线的斜抛物线模型及卫星轨道的位置坐标,求得导线的入射角;
步骤4.根据导线的入射角,构建并计算雷达辐射区域导线计算模型,得到对应的导线目标中心点位置辐射能量值;
步骤5.根据步骤1至4的计算结果,生成仿真区域背景,并叠加和仿真导线。
如图2所示,步骤1,包括:
1-1.从卫星平台提供的卫星雷达数据头文件中读取卫星轨道离散采样状态矢量数据;卫星轨道离散采样状态矢量数据包括离散采样卫星位置矢量RS(ti)、离散采样卫星速度矢量VS(ti)及采样对应的时刻ti(i=1,2,…,n);
RS(ti)=(XS(ti),YS(ti),ZS(ti))T;
VS(ti)=(VSX(ti),VSY(ti),VSZ(ti))T;
式中,XS(ti),YS(ti),ZS(ti)为卫星在x轴、y轴和Z轴上的离散采样位置矢量;
VSX(ti),VSY(ti),VSZ(ti)分别为卫星在x轴、y轴和Z轴上的离散采样速度矢量;
1-2.构建离散采样卫星位置矢量RS(ti)的三次多项式拟合轨道方程:
构建离散采样卫星速度矢量VS(ti)的三次多项式拟合轨道方程:
式中,(a0,a1,a2,a3,b0,b1,b2,b3,c0,c1,c2,c3)为三次多项式拟合轨道方程的拟合参数;
1-3.根据RS(ti)及VS(ti),采用最小估计二乘法求解三次多项式拟合轨道方程,得到拟合参数(a0,a1,a2,a3,b0,b1,b2,b3,c0,c1,c2,c3)的值;
1-4.根据拟合参数(a0,a1,a2,a3,b0,b1,b2,b3,c0,c1,c2,c3)的值,求得卫星轨道任意时刻星载合成孔径雷达的卫星轨道任意时刻t的轨道位置矢量RS(t)和速度矢量VS(t),即求得任意时刻的轨道位置坐标(xs,ys,zs)和速度坐标(Vsx,Vsy,Vsz);
其中,RS(t)的拟合轨道方程为:
VS(t)的拟合轨道方程为:
式中,XS(t),YS(t),ZS(t)分别为卫星在x轴、y轴和Z轴上的任意时刻t的位置矢量;VSX(t),VSY(t),VSZ(t)分别为卫星在x轴、y轴和Z轴上的任意时刻t的速度矢量。
如图3所示,步骤2,包括:
2-1.构建输电线路中的导线进行斜抛物线模型:
当输电线路不平行于X轴时:
当输电线路平行于X轴时:
式中,m、n和R为导线所在重力平面方程系数,r1,r2,r3为导线在以z轴为中心的旋转抛物面方程系数,其中g和ρ0为导线相关系数参数;x,y,z为导线铁塔空间坐标。
2-2.若导线不平行于x坐标轴,则r1,r2,r3的值为:
r2=mr1+n;
r3=zt1-A[(xt1-r1)2+(yt1-r2)2];
其中,n=yt1-mxt1;xt1,yt1,zt1和xt2,yt2,zt2分别为所述导线两端铁塔的空间坐标;
2-3.若导线平行于x坐标轴,则r1,r2,r3的值为:
r1=R;
r3=zt1-A(yt1-r2)2。
中的求解方法包括:
根据已知导线的两端铁塔的位置坐标(xt1,yt1,zt1)和(xt2,yt2,zt2),求出高差角的余弦值:
如图4所示,步骤3,包括:
3-1.任意时刻的卫星轨道位置坐标(xs,ys,zs)和目标点位置坐标(xp,yp,zp),基于距离公式计算SAR数据中对应的斜距r:
r2=(xs-xp)2+(ys-yp)2+(zs-zp)2;
3-2.对于导线上的任一点,对导线的斜抛物线模型方程求偏导,获得导线的方向向量
若导线平行于X轴,则方向向量为(0,1,2A(y-r2));
若导线不平行于X轴,
则方向向量为(1,m,2A(x-r1)+2mA(mx+n-r2));
3-3.利用入射向量和导线方向向量基于空间几何夹角公式计算导线的入射角θ:
式中,雷达电磁波的入射向量
如图5所示,步骤4,包括:
4-1.根据导线的入射角θ,构建雷达辐射区域点目标散射模型:
e=σ0(cosθ)2;
式中,e为雷达辐射区域点目标散射模型,σ0为后向散射系数,θ为雷达照射目标区域入射角;
4-2.确定仿真数据星载SAR相关参数;
4-3.构建雷达辐射区域导线计算模型;
对于给定时刻t卫星平台位置(XS(t),YS(t),ZS(t))T及正侧视中心导线位置坐标(xp,yp,zp)T,确定雷达辐射带宽范围内的导线坐标为:
[(xp-,yp-,zp-)T,(xp+,yp+,zp+)T];
4-4.计算并叠加雷达辐射带宽区域所有导线上的点目标的辐射能量值,生成对应的导线目标中心点位置辐射能量值。
其中,相关参数包括卫星飞行高度约为、辐射波长、波束宽度方位向、距离向、在方位向辐射宽度及辐射分辨率。
如图6所示,步骤5,包括:
5-1.采用光学灰度图像作为导线区域背景,对图像叠加高斯乘性噪声模拟SAR图像;
5-2.基于斜距确定导线目标在背景图像中的坐标位置,合并背景区域和导线目标仿 真能量值,生成仿真图像。
如图7所示,本发明提供一种输电线路在星载合成孔径雷达成像中的电磁仿真方法的应用例,具体为:
(一)卫星轨道状态矢量多项式拟合计算
1)从卫星雷达数据头文件中读取卫星轨道状态矢量;
为了拟合计算目标图像成像时刻的卫星轨道,我们需要卫星平台提供至少5个点的卫星状态矢量,包括位置和速度矢量以及其对应的时刻,为此我们从卫星雷达数据头文件中读取卫星轨道离散采样状态矢量,分别用卫星位置矢量
RS(ti)=(XS(ti),YS(ti),ZS(ti))T和速度矢量VS(ti)=(VSX(ti),VSY(ti),VSZ(ti))T表示,
其中ti(i=1,2,…,n)为离散采样点对应的时间;
2)多项式拟合参数确定
构建三次多项式拟合轨道方程,其表达式为:
其中,(a0,a1,a2,a3,b0,b1,b2,b3,c0,c1,c2,c3)为轨道多项式拟合参数,利用离散卫星轨道离散采样状态矢量RS(ti)和VS(ti),采用最小二乘法估计解算这些待定参数;
3)卫星轨道任意时刻t轨道状态矢量计算
通过步骤2)中计算的(a0,a1,a2,a3,b0,b1,b2,b3,c0,c1,c2,c3),精确拟合t时刻的卫星状态矢量RS(t)和VS(t),其中
(二)导线三维建模及精确计算
1)基于斜抛物线模型的导线位置矢量精确计算
在保证计算精度的前提下,用斜抛物线模型简化计算导线空间位置。其建模方程为,
(如果导线不平行x轴)
特殊情况下,
(如果导线平行x轴)
式中,m,n和R为导线所在重力平面方程系数,r1,r2,r3为导线在以z轴为中心的旋转抛物面方程系数,其中g和ρ0为导线相关系数参数,不同输电线路都有先验值或者可以实时测量。那么,已知导线两端铁塔的位置坐标(xt1,yt1,zt1)和(xt2,yt2,zt2),则可以求出高差角的余弦值。
当导线不平行于x坐标轴的时候;n=yt1-mxt1,
r2=mr1+n,r3=zt1-A[(xt1-r1)2+(yt1-r2)2]
当导线平行于x坐标轴的时候;
(三)SAR几何构像和散射模型建模计算
1)斜距计算
成孔径雷达雷达成像方式与光学成像不一样,SAR成像为等斜距成像,分辨率反应在斜距上,只能区分不同距离目标的回波信号。
其简化距离公式为:
r2=(xs-xp)2+(ys-yp)2+(zs-zp)2
式中(xs,ys,zs)为卫星轨道位置坐标,(xp,yp,zp)为目标点位置坐标。
2)对于导线任一点,基于导线三维建模,对导线建模方程求偏导,获得导线的方向向量:(0,1,2A(y-r2))(导线平行x轴)或者(1,m,2A(x-r1)+2mA(mx+n-r2))(导线不平行x轴)
3)计算导线入射角θ
关于入射角θ的计算,我们引入空间几何计算方法,在三维空间中两条直线的夹角可以用公式:
式中为雷达电磁波入射向量,为对应导线点的方向向量, 或者
(四)导线点目标雷达辐射建模计算
1)雷达辐射区域点目标散射模型为e=σ0(cosθ)2,式中σ0为后向散射系数,θ为雷达照射目标区域入射角;
2)仿真数据星载SAR相关参数确定
以经典星载SAR条带模式为例分析,参考德国TerraSAR-X卫星相关参数,飞行高度约为514km,辐射波长为3.2cm,波束宽度方位向0.33°和距离向2.3°,在方位向辐射宽度为2000-3000m,辐射分辨率1-3m;所以本仿真算法选取与TerraSAR-X相同参数,方位向辐射宽度为2000m,辐射分辨率1m;
3)雷达辐射区域导线计算模型
对于给定时刻t卫星平台位置(XS(t),YS(t),ZS(t))T,及正侧视中心导线位置坐标(xp,yp,zp)T,确定其辐射带宽范围内的导线坐标[(xp-,yp-,zp-)T,(xp+,yp+,zp+)T],分别计算1m间隔采样点集所有点的入射角和后向散射系数,基于入射角大于等于89°的采样点进行辐射能量叠加,生成该时刻对应的导线目标位置辐射能量值;
(五)仿真区域背景生成
1)采用光学灰度图像作为导线区域背景,对图像叠加高斯乘性噪声模拟SAR图像;
2)基于斜距确定导线目标在背景图像中的坐标位置,合并背景区域和导线目标仿真能量值,生成仿真图像;
至此输电线路导线区域电磁辐射能量值计算完毕,获得了较为理想的导线区域仿真SAR图像,完成了本发明技术的所有步骤。
以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员依然可以对本发明的具体实施方式进行修改或者等同替换,而这些未脱离本发明精神和范围的任何修改或者等同替换,其均 在申请待批的本发明的权利要求保护范围之内。
Claims (7)
1.一种输电线路在星载合成孔径雷达成像中的电磁仿真方法,其特征在于,所述方法包括如下步骤:
步骤1.根据卫星轨道离散采样状态矢量数据,计算得到星载合成孔径雷达的卫星轨道任意时刻的轨道状态矢量;
步骤2.构建输电线路中的导线的斜抛物线模型;
步骤3.根据所述导线的斜抛物线模型及卫星轨道的位置坐标,求得所述导线的入射角;
步骤4.根据所述导线的入射角,构建并计算雷达辐射区域导线计算模型,得到对应的导线目标中心点位置辐射能量值;
步骤5.根据计算结果,生成仿真区域背景,并叠加和仿真导线;
所述步骤3,包括:
3-1.任意时刻的卫星轨道位置坐标(xs,ys,zs)和目标点位置坐标(xp,yp,zp),基于距离公式计算SAR数据中对应的斜距r:
r2=(xs-xp)2+(ys-yp)2+(zs-zp)2;
3-2.对于所述导线上的任一点,对所述导线的斜抛物线模型方程求偏导,获得导线的方向向量
若所述导线平行于X轴,则所述方向向量为(0,1,2A(y-r2));
若所述导线不平行于X轴,
则所述方向向量为(1,m,2A(x-r1)+2mA(mx+n-r2));
3-3.利用入射向量和导线方向向量基于空间几何夹角公式计算所述导线的入射角θ:
式中,雷达电磁波的入射向量
2.如权利要求1所述的方法,其特征在于,所述步骤1,包括:
1-1.从卫星平台提供的卫星雷达数据头文件中读取卫星轨道离散采样状态矢量数据;所述卫星轨道离散采样状态矢量数据包括离散采样卫星位置矢量RS(ti)、离散采样卫星速度矢量VS(ti)及采样对应的时刻ti(i=1,2,…,n);
RS(ti)=(XS(ti),YS(ti),ZS(ti))T;
VS(ti)=(VSX(ti),VSY(ti),VSZ(ti))T;
式中,XS(ti),YS(ti),ZS(ti)为卫星在x轴、y轴和Z轴上的离散采样位置矢量;
VSX(ti),VSY(ti),VSZ(ti)分别为卫星在x轴、y轴和Z轴上的离散采样速度矢量;
1-2.构建离散采样卫星位置矢量RS(ti)的三次多项式拟合轨道方程:
构建离散采样卫星速度矢量VS(ti)的三次多项式拟合轨道方程:
式中,(a0,a1,a2,a3,b0,b1,b2,b3,c0,c1,c2,c3)为三次多项式拟合轨道方程的拟合参数;
1-3.根据RS(ti)及VS(ti),采用最小估计二乘法求解所述三次多项式拟合轨道方程,得到拟合参数(a0,a1,a2,a3,b0,b1,b2,b3,c0,c1,c2,c3)的值;
1-4.根据拟合参数(a0,a1,a2,a3,b0,b1,b2,b3,c0,c1,c2,c3)的值,求得卫星轨道任意时刻星载合成孔径雷达的卫星轨道任意时刻t的轨道位置矢量RS(t)和速度矢量VS(t),即求得任意时刻的轨道位置坐标(xs,ys,zs)和速度坐标(Vsx,Vsy,Vsz);
其中,RS(t)的拟合轨道方程为:
VS(t)的拟合轨道方程为:
式中,XS(t),YS(t),ZS(t)分别为卫星在x轴、y轴和Z轴上的任意时刻t的位置矢量;VSX(t),VSY(t),VSZ(t)分别为卫星在x轴、y轴和Z轴上的任意时刻t的速度矢量。
3.如权利要求2所述的方法,其特征在于,所述步骤2,包括:
2-1.构建所述输电线路中的导线进行斜抛物线模型:
当所述输电线路不平行于X轴时:
当所述输电线路平行于X轴时:
式中,m、n和R为所述导线所在重力平面方程系数,r1,r2,r3为导线在以z轴为中心的旋转抛物面方程系数,其中g和ρ0为导线相关系数参数;x,y,z为导线铁塔空间坐标;
2-2.若所述导线不平行于x坐标轴,则r1,r2,r3的值为:
r2=mr1+n;
r3=zt1-A[(xt1-r1)2+(yt1-r2)2];
其中,n=yt1-mxt1;xt1,yt1,zt1和xt2,yt2,zt2分别为所述导线两端铁塔的空间坐标;
2-3.若所述导线平行于x坐标轴,则r1,r2,r3的值为:
r1=R;
r3=zt1-A(yt1-r2)2。
4.如权利要求3所述的方法,其特征在于,所述中的求解方法包括:
根据已知所述导线的两端铁塔的位置坐标(xt1,yt1,zt1)和(xt2,yt2,zt2),求出高差角的余弦值:
5.如权利要求1所述的方法,其特征在于,所述步骤4,包括:
4-1.根据所述导线的入射角θ,构建雷达辐射区域点目标散射模型:
e=σ0(cosθ)2;
式中,e为雷达辐射区域点目标散射模型,σ0为后向散射系数,θ为雷达照射目标区域入射角;
4-2.确定仿真数据星载SAR相关参数;
4-3.构建雷达辐射区域导线计算模型;
对于给定时刻t卫星平台位置(XS(t),YS(t),ZS(t))T及正侧视中心导线位置坐标(xp,yp,zp)T,确定所述雷达辐射带宽范围内的导线坐标为:
[(xp-,yp-,zp-)T,(xp+,yp+,zp+)T];
4-4.计算并叠加雷达辐射带宽区域所有导线上的点目标的辐射能量值,生成对应的导线目标中心点位置辐射能量值。
6.如权利要求5所述的方法,其特征在于,所述相关参数包括卫星飞行高度、辐射波长、波束宽度方位向、距离向、在方位向辐射宽度及辐射分辨率。
7.如权利要求6所述的方法,其特征在于,所述步骤5,包括:
5-1.采用光学灰度图像作为导线区域背景,对图像叠加高斯乘性噪声模拟SAR图像;
5-2.基于斜距确定导线目标在背景图像中的坐标位置,合并背景区域和导线目标仿真能量值,生成仿真图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510543429.XA CN106483511B (zh) | 2015-08-28 | 2015-08-28 | 一种输电线路在星载合成孔径雷达成像中的电磁仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510543429.XA CN106483511B (zh) | 2015-08-28 | 2015-08-28 | 一种输电线路在星载合成孔径雷达成像中的电磁仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106483511A CN106483511A (zh) | 2017-03-08 |
CN106483511B true CN106483511B (zh) | 2018-12-11 |
Family
ID=58235085
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510543429.XA Active CN106483511B (zh) | 2015-08-28 | 2015-08-28 | 一种输电线路在星载合成孔径雷达成像中的电磁仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106483511B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108919261A (zh) * | 2018-03-27 | 2018-11-30 | 中国电力科学研究院有限公司 | 一种基于合成孔径雷达影像确定高压输电线路弧垂的方法及系统 |
CN110471061A (zh) * | 2019-07-16 | 2019-11-19 | 青岛擎鹰信息科技有限责任公司 | 一种实现机载合成孔径雷达成像的仿真方法及其系统 |
CN112990071B (zh) * | 2021-03-31 | 2024-01-05 | 香港中文大学(深圳) | 一种输电线路异常检测方法、装置及计算机可读存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101907704A (zh) * | 2010-06-11 | 2010-12-08 | 西安电子科技大学 | 多模式合成孔径雷达仿真成像评估方法 |
CN103163509A (zh) * | 2013-02-03 | 2013-06-19 | 北京航空航天大学 | 一种基于电磁散射的高频近似法的模拟合成孔径雷达 |
CN103197291A (zh) * | 2013-03-12 | 2013-07-10 | 北京航空航天大学 | 一种基于非停走模型的星载sar回波信号仿真方法 |
EP2762918A2 (en) * | 2014-01-16 | 2014-08-06 | Institute of Electronics, Chinese Academy of Sciences | Method for analyzing effect of sub-band interference on imaging performance in synthetic aperture radar |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101320508B1 (ko) * | 2012-09-13 | 2013-10-23 | 국방과학연구소 | 항공기에 탑재된 영상 레이더의 성능 시험 장치 |
-
2015
- 2015-08-28 CN CN201510543429.XA patent/CN106483511B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101907704A (zh) * | 2010-06-11 | 2010-12-08 | 西安电子科技大学 | 多模式合成孔径雷达仿真成像评估方法 |
CN103163509A (zh) * | 2013-02-03 | 2013-06-19 | 北京航空航天大学 | 一种基于电磁散射的高频近似法的模拟合成孔径雷达 |
CN103197291A (zh) * | 2013-03-12 | 2013-07-10 | 北京航空航天大学 | 一种基于非停走模型的星载sar回波信号仿真方法 |
EP2762918A2 (en) * | 2014-01-16 | 2014-08-06 | Institute of Electronics, Chinese Academy of Sciences | Method for analyzing effect of sub-band interference on imaging performance in synthetic aperture radar |
Also Published As
Publication number | Publication date |
---|---|
CN106483511A (zh) | 2017-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108181635B (zh) | 一种用于输电线路交叉跨越分析的激光点云分类方法 | |
CN103135550B (zh) | 用于电力巡线的无人机多重避障控制方法 | |
US20210055180A1 (en) | Apparatuses and methods for gas flux measurements | |
Teng et al. | Mini-UAV LiDAR for power line inspection | |
CN202046439U (zh) | 电力巡线无人直升机超低空飞行障碍规避子系统 | |
CN104932001A (zh) | 一种实时3d核辐射环境重建监测系统 | |
CN109829908B (zh) | 基于双目影像的电力线下方地物安全距离检测方法及设备 | |
WO2012087387A2 (en) | Thermal powerline rating and clearance analysis using thermal imaging technology | |
CN109521792A (zh) | 一种基于输变电设备三维模型的无人机地面飞行控制系统 | |
CN106483511B (zh) | 一种输电线路在星载合成孔径雷达成像中的电磁仿真方法 | |
CN106683039B (zh) | 一种生成火情态势图的系统 | |
CN104122560A (zh) | 一种输电线路广域冰情监测方法 | |
CN114743021A (zh) | 一种输电线路图像与点云数据的融合方法及系统 | |
CN103697883B (zh) | 一种基于天际线成像的飞行器水平姿态确定方法 | |
Wood et al. | Measurement of three dimensional volcanic plume properties using multiple ground based infrared cameras | |
Jiang et al. | UAV-based oblique photogrammetry for 3D reconstruction of transmission line: Practices and applications | |
CN103743488B (zh) | 遥感卫星地球临边背景特性的红外成像仿真方法 | |
Zhang et al. | Power line simulation for safety distance detection using point clouds | |
CN115828616A (zh) | 一种深空红外微动目标模型构建方法及成像仿真系统 | |
Hao et al. | Extraction and analysis of tree canopy height information in high-voltage transmission-line corridors by using integrated optical remote sensing and LiDAR | |
CN110389353A (zh) | 一种森林火灾火点精准定位方法及系统 | |
US20240028051A1 (en) | Closed surface flight pattern generation for unmanned aerial vehicle (uav) flux plane assessment of large facilities | |
CN209746157U (zh) | 一种小型移动涡动相关湍流水热通量观测系统 | |
Safdarinezhad et al. | A photogrammetric solution for measurement of power lines sag via integration of image and accelerometer data of a smartphone | |
Filippone et al. | GIS Systems and LIDAR technology for the operation of HV lines. Sicilian Transmission network applications |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |