
Citation: | PANG Lei, ZHANG Fengli, WANG Guojun, et al. Imaging simulation and damage assessment feature analysis of Ku band polarized SAR of buildings[J]. Journal of Radars, 2020, 9(3): 578–587. doi: 10.12000/JR20061 |
目标的振动、转动等微动产生的微多普勒效应包含了目标的结构和运动信息,常用于目标的分类和识别[1–3]。目前,基于外辐射源雷达微多普勒效应目标分类和识别的研究还处于起步状态。外辐射源雷达是一种利用非合作照射源进行目标探测和分类识别的新体制雷达系统,其自身不辐射电磁能量,具有节约频谱资源,隐蔽性好,设备规模小,易于部署和组网等特点[4–6]。在微多普勒效应目标分类和识别方面,外辐射源雷达表现出得天独厚的优势:(1)收发分置可实现空间分集,有效避免探测盲区。(2)第三方辐射源多为连续波,长时间相干积累可记录多个连续的回波闪烁,同时有利于提高对低雷达散射截面积(Radar Cross-Section, RCS)微动目标的探测与分类识别能力。(3)对微多普勒特征的提取不要求高距离分辨率,参数估计不受第三方辐射源带宽的限制[7,8]。
针对微多普勒效应参数估计问题。文献[9,10]中依据微动目标正弦特征曲线,利用 Hough变换,在参数域中进行多维搜索提取出微动曲线进行参数估计。文献[11,12]通过正交匹配追踪(Orthogonal Matching Pursuit, OMP)算法进行稀疏逼近实现了微动目标参量的估计。上述方法均具有较好的鲁棒性,但由于估计参量维数较高导致计算量巨大。文献[13]利用微动目标在时频域的周期性,采用循环相关系数方法,实现了目标微动周期的估计,但信号周期较长时计算量急剧增加。文献[14]计算了信号的高阶矩函数,通过检测在不同时延下,高阶虚函数部分傅里叶变换累计结果的峰值位置,快速获得目标的旋转速率,相比于图像处理方法和OMP分解方法,计算复杂度较小,但抗噪性能差。而外辐射源雷达所利用的第三方辐射源多为连续波信号,其发射波形不可控,信号能量主要覆盖地面,杂波环境复杂且对空中目标增益低,利用长时间相干积累来提高处理增益会带来数据量巨大的挑战。上述因素决定了外辐射源雷达参数估计方法需要有良好的抗噪性能且计算量要小。
直升机旋翼旋转时对雷达信号产生周期性调制,当叶片发生镜面反射时,旋翼回波出现峰值,即回波闪烁。闪烁信号在时频图像中表现为一定宽度的频率带,且闪烁时间、闪烁间隔与直升机旋翼微动参数密切相关。针对外辐射源雷达参数估计问题,本文结合上述时频域中闪烁信号的特点,通过时频分析和正交匹配追踪算法实现了直升机旋翼微动参数的估计。本文首先给出了外辐射源雷达直升机旋翼微动信号模型,其次介绍了如何在时频图中提取出闪烁信号参数及正交匹配追踪算法对直升机旋翼微动参数的估计,最后仿真和实测证明了本文方法的有效性。
直升机旋转叶片与外辐射源雷达的位置关系如图1所示。以直升机旋转叶片的中心点为原点
假设直升机平动得到补偿。在
rP(t)=||RT −RP(t)||+||RR−RP(t)|| | (1) |
其中
参考文献[1]中单基地直升机建模,将叶片看作线模型,外辐射源雷达直升机旋翼回波可表示为:
s(t)=Lexp{−j2πλ(rR+rT)}N∑k=1sinc{−ϕk(t)}⋅exp{jϕk(t)} | (2) |
其中,
ϕk(t)=−4πλL2cosβcos(α−γ2)cos(φk(t)) | (3) |
φk(t)=2πfrt+φ0+(k−1)2π/N−α+γ2 | (4) |
由式(3)得第
fk(t)=2πfrLλcosβcos(α−γ2)sin(φk(t)) | (5) |
由式(2)可知时域信号幅值受
φk(t)=±π2+2πn | (6) |
由式(2)知连续两个闪烁之间的时间间隔为:
Δt={12N⋅fr,N为奇数1N⋅fr, N为偶数 | (7) |
直升机旋翼回波的微多普勒呈非线性变化,通过对目标回波信号进行时频分析能够揭示信号频率的时变特性。短时傅里叶变化(Short-Time Fourier Transform, STFT)计算简单,且不产生交叉项。对直升机旋翼回波信号
TF(t,f)=∫s(τ)w(τ−t)e−j2πfτdτ | (8) |
其中,
图2为直升机旋翼回波的时频图。当直升机旋翼的叶片数为奇数时(图2(a)),时频域中正负多普勒“闪烁”交替出现;若旋翼叶片数为偶数(图2(b)),则是同时出现。
设时频域中正频率“闪烁”发生的时间为
φk(t0)=π2+2πn(n为整数) | (9) |
由式(4)和式(9)得第
φ0={−πt0Δt1N−2π(k−1)1N+φ1, N为奇数−2πt0Δt1N−2π(k−1)1N+φ1,N为偶数 | (10) |
其中
φ1=α+γ2+2πn+π2(0≤φ0<2π) | (11) |
由于时频图像中闪烁信号频率带垂直于时间横轴,对正频率轴数据幅值进行累加计算,并判断累加后数据局部峰值点,可得到时频域中正频率“闪烁”发生的时间。同样,对负频率轴数据幅值进行累加计算得到时频域中负频率“闪烁”发生的时间。相应的也可得到闪烁间隔。
由式(7)知,闪烁间隔与旋翼转速、叶片数量密切相关。由式(10)知,闪烁发生的时间与叶片初相、叶片数量、整数
由式(2)知时域回波信号可分解为:
s(t)=M∑m=1cmgm(t;Λ)=Dα | (12) |
其中,
由式(2)知直升机旋翼回波信号由参数
由OMP算法原理可知,字典中的原子可按照待分解信号的内在特性来构造[16]。根据微动目标的时域回波表达式(2),第
a(m)=sinc(ϕ(Lr,Np,kq))⋅exp{−jϕ(Lr,Np,kq)} | (13) |
其中
m=rpq | (14) |
并对原子集里的每个原子进行能量归一化:
a(m)←a(m)/‖a(m)‖F | (15) |
其中,
将5参量
直升机旋翼参数估计具体步骤如下:
步骤1 对直升机旋翼信号进行短时傅里叶变换,得到时频图像
步骤2 对时频图中正频率轴数据幅值进行累加计算,并判断累加后数据局部峰值点,对应时频域正频率“闪烁”发生的时间。同样,对负频率轴数据幅值进行累加计算得到时频域中负频率“闪烁”发生的时间。
步骤3 根据步骤2中正负频率“闪烁”发生的时间,判别时频域中正负多普勒“闪烁”是否交替出现。若是,则旋翼叶片数为奇数,否则,旋翼叶片数为偶数。
步骤4 读取某一正频率闪烁发生的时间
步骤5 确定
步骤6 利用OMP算法寻找叶片数量,叶片长度的最优值,代入式(7)计算出旋翼转速,代入式(10)和式(11)计算出叶片初相。
结合上述模型对直升机旋翼回波信号进行仿真,仿真参数设置如表1所示。
信号载频 | 叶片数 | 叶片长度 | 旋转速率 | 发射站方位角 | 接收站方位角 | 发射站仰角 | 接收站仰角 | SNR |
658 MHz | 3 | 5 m | 200 rpm | 33° | 76° | 23° | 23° | –5 dB |
图3(a)显示了信号的联合时频域特征,可看出闪烁信号及噪声严重影响直升机旋翼微多普勒特征曲线的检测,使微多普勒特征曲线提取困难。
分别对时频图像中正负频率轴数据幅值进行累加计算,得到时频域中正负多普勒“闪烁”时间,如图3(b)所示,图中正负多普勒“闪烁”等间隔交替出现,则旋翼叶片数为奇数。读取闪烁信号时间间隔为0.05 s,根据式(7)表示出旋翼转速为:
fr=10/N | (16) |
读取某一正频率闪烁信号对应时刻为0.066 s(此处选择了图3(b)中的第1个正频率闪烁信号),根据式(10)和式(11)表示出第
φ0=−4.14/N−6.28×(k−1)/N+2.52 | (17) |
图3(c)为利用OMP方法对
图4为利用常规Hough变换,通过微多普勒曲线
{f_{\max }} = \frac{{4{{π}} {f\!_r}L}}{\lambda }\cos\beta \cos\left( \frac{{\alpha - \gamma }}{2}\right) | (18) |
图4中给出了参数空间中局部峰值点中心位置。可得到直升机旋翼转速为200 rpm。3叶片最大频移分别为385.6 Hz, 393.9 Hz, 389.8 Hz,平均值为390.0 Hz,由式(18)计算得叶片长度为4.96 m,与理论值基本一致。3叶片初相分别为0.91 rad, 3.16 rad, 5.24 rad,利用式(4)对初相进行修正,得到3叶片初相位为1.86 rad, 4.11 rad, 6.19 rad,存在较大的误差,是由于STFT受不确定原理的限制,时频图像中时频分辨率受限使参数空间中的局部峰值点扩展范围较大,只能大致估计局部峰值点的位置,估计结果精度较低。
设待处理的时频图像大小为
直接使用OMP进行参数
本文方法计算量集中在OMP阶段,根据提取的时频域中的闪烁时间,依据式(7)和式(10),最终转化为对参数
在对直升机旋翼微动参数估计时,一般
武汉大学电波传播实验室对EC_120B直升机进行了微多普勒效应探究外场实验,EC_120B直升机主旋翼3叶片,叶片长度5 m,额定转速406 rpm,实验中以武汉龟山电视塔数字电视信号为照射源,信号中心频率为658 MHz,带宽8 MHz,接收站位于武汉大学电波传播实验室楼顶,距离发射站7.56 km,实验场景如图5所示。本组实测数据相干积累时间0.8 s,可近似认为目标在这段时间位置不变,直升机旋翼转速为常量。
图6(a)为去除目标主体影响后,对直升机旋翼回波信号进行短时傅里叶变换后的时频图像。可以观察到闪烁信号,但微多普勒特征曲线已观察不到。分别对时频图像中正负频率轴数据幅值进行累加计算,得到时频域中正负多普勒“闪烁”时间,如图6(b)所示,图中正负多普勒“闪烁”等间隔交替出现,则旋翼叶片数为奇数。
读取闪烁信号时间间隔为26.2 ms,根据式(7)用
本文根据外辐射源雷达直升机旋翼微动信号模型,充分利用时频域中闪烁信号特征和微动信号内在特性进行了参数估计。通过时频分析和正交匹配追踪算法,估计出了旋翼转速、叶片长度、叶片数量和初相。同时开展了外场实验。仿真数据和实测数据处理都表明本文方法对外辐射源雷达直升机旋翼参数估计的可行性。
[1] |
Tadono T, Ohki M, and Abe T. Summary of natural disaster responses by the Advanced Land Observing Satellite-2(ALOS-2)[C]. The International Archives of the Photogrammetry, Remote Sensing & Spatial Information Sciences, Kyoto, Japan, 2019: 69–72. doi: 10.5194/isprs-archives-XLII-3-W7-69-2019.
|
[2] |
张风丽, 邵芸. 城市目标高分辨率SAR遥感监测技术研究进展[J]. 遥感技术与应用, 2010, 25(3): 415–422. doi: 10.11873/j.issn.1004-0323.2010.3.415
ZHANG Fengli and SHAO Yun. Urban target monitoring using high resolution SAR data[J]. Remote Sensing Technology and Application, 2010, 25(3): 415–422. doi: 10.11873/j.issn.1004-0323.2010.3.415
|
[3] |
CHEN Qihao, LI Linlin, JIANG Ping, et al. Building collapse extraction using modified freeman decomposition from post-disaster polarimetric SAR image[C]. 2016 IEEE International Geoscience And Remote Sensing Symposium (IGARSS), Beijing, China, 2016: 5769–5772. doi: 10.1109/IGARSS.2016.7730507.
|
[4] |
WANG Xiaoqing, DOU Aixia, JIN Dingjian, et al. Study on the comparation of building damage extracted from different RS images acquired after 2010 M = 7.1 Yushu, Qinghai, China earthquake[C]. 2012 IEEE International Geoscience and Remote Sensing Symposium, Munich, Germany, 2012: 935–938. doi: 10.1109/IGARSS.2012.6351401.
|
[5] |
ZHANG Haizhen, WANG Qing, ZENG Qiming, et al. A novel approach to building collapse detection from post-seismic polarimetric SAR imagery by using optimization of polarimetric contrast enhancement[C]. 2015 IEEE International Geoscience And Remote Sensing Symposium (IGARSS), Milan, Italy, 2015: 3270–3273. doi: 10.1109/IGARSS.2015.7326516.
|
[6] |
PARK S E, YAMAGUCHI Y, and KIM D J. Polarimetric SAR remote sensing of the 2011 Tohoku earthquake using ALOS/PALSAR[J]. Remote sensing of Environment, 2013, 132: 212–220. doi: 10.1016/j.rse.2013.01.018
|
[7] |
ZHAI Wei, SHEN Huanfeng, HUANG Chunlin, et al. Building earthquake damage information extraction from a single post-earthquake PolSAR image[J]. Remote Sensing, 2016, 8(3): 171. doi: 10.3390/rs8030171
|
[8] |
SATO M, CHEN Siwei, and SATAKE M. Polarimetric SAR analysis of tsunami damage following the March 11, 2011 East Japan earthquake[J]. Proceedings of the IEEE, 2012, 100(10): 2861–2875. doi: 10.1109/JPROC.2012.2200649
|
[9] |
WATANABE M, MOTOHKA T, MIYAGI Y, et al. 2012. Analysis of urban areas affected by the 2011 off the pacific coast of Tohoku earthquake and tsunami with L-Band SAR full-polarimetric mode[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(3): 472–476. doi: 10.1109/LGRS.2011.2182030
|
[10] |
ZHAI Wei and HUANG Chunlin. Fast building damage mapping using a single post-earthquake PolSAR image: A case study of the 2010 Yushu earthquake[J]. Earth, Planets and Space, 2016, 68(1): 86. doi: 10.1186/s40623-016-0469-2
|
[11] |
ZHAI Wei, HUANG Chunlin, and PEI Wansheng. Two new polarimetric feature parameters for the recognition of the different kinds of buildings in earthquake-stricken areas based on entropy and eigenvalues of PolSAR decomposition[J]. Remote Sensing, 2018, 10(10): 1613. doi: 10.3390/rs10101613
|
[12] |
LI Xinwu, GUO Huadong, ZHANG Lu, et al. A new approach to collapsed building extraction using RADARSAT-2 polarimetric SAR imagery[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(4): 677–681. doi: 10.1109/LGRS.2011.2178392
|
[13] |
SUN Weidong, SHI Lei, YANG Jie, et al. Building collapse assessment in urban areas using texture information from postevent SAR data[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2016, 9(8): 3792–3808. doi: 10.1109/JSTARS.2016.2580610
|
[14] |
ZHAO Lingli, YANG Jie, LI Pingxiang, et al. Damage assessment in urban areas using post-earthquake airborne PolSAR imagery[J]. International Journal of Remote Sensing, 2013, 34(24): 8952–8966. doi: 10.1080/01431161.2013.860566
|
[15] |
GE Pinglan, GOKON H, and MEGURO K. A review on synthetic aperture radar-based building damage assessment in disasters[J]. Remote Sensing of Environment, 2020, 240: 111693. doi: 10.1016/j.rse.2020.111693
|
[16] |
CHEN Siwei, WANG Xuesong, and SATO M. Urban damage level mapping based on scattering mechanism investigation using fully polarimetric SAR data for the 3.11 East Japan earthquake[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(12): 6919–6929. doi: 10.1109/TGRS.2016.2588325
|
[17] |
王雪松, 陈思伟. 合成孔径雷达极化成像解译识别技术的进展与展望[J]. 雷达学报, 2020, 9(2): 259–276. doi: 10.12000/JR19109
WANG Xuesong and CHEN Siwei. Polarimetric synthetic aperture radar interpretation and recognition: Advances and perspectives[J]. Journal of Radars, 2020, 9(2): 259–276. doi: 10.12000/JR19109
|
[18] |
FREEMAN A and DURDEN S L. A three-component scattering model for polarimetric SAR data[J]. IEEE Transactions on Geoscience and Remote Sensing, 1998, 36(3): 963–973. doi: 10.1109/36.673687
|
[19] |
YAMAGUCHI Y, SATO A, BOERNER W M, et al. Four-component scattering power decomposition with rotation of coherency matrix[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(6): 2251–2258. doi: 10.1109/TGRS.2010.2099124
|
[20] |
JI Kefeng and WU Yonghui. Scattering mechanism extraction by a modified cloude-pottier decomposition for dual polarization SAR[J]. Remote Sensing, 2015, 7(6): 7447–7470. doi: 10.3390/rs70607447
|
[21] |
CHIEN L L C, EWE H T, and SAW S H. Understanding the correlation in scattering mechanisms between H-Alpha decomposition and theoretical modelling[C]. 2018 Progress in Electromagnetics Research Symposium (PIERS-Toyama), Toyama, Japan, 2018: 1198–1203. doi: 10.23919/PIERS.2018.8597662.
|
[22] |
TOUZI R. Target scattering decomposition in terms of roll-invariant target parameters[J]. IEEE Transactions on Geoscience and Remote Sensing, 2007, 45(1): 73–84. doi: 10.1109/TGRS.2006.886176
|
[23] |
AN Wentao, XIE Chunhua, YUAN Xinzhe, et al. Four-component decomposition of polarimetric SAR images with deorientation[J]. IEEE Geoscience and Remote Sensing Letters, 2011, 8(6): 1090–1094. doi: 10.1109/LGRS.2011.2157078
|
[24] |
孙翔, 宋红军, 王宇, 等. 基于高分辨率全极化SAR图像的取向角校正方法[J]. 雷达学报, 2018, 7(4): 465–474. doi: 10.12000/JR18026
SUN Xiang, SONG Hongjun, WANG R, et al. POA correction method using high-resolution full-polarization SAR image[J]. Journal of Radars, 2018, 7(4): 465–474. doi: 10.12000/JR18026
|
[25] |
LIU Shan, ZHANG Fengli, WEI Shiying, et al. Building damage mapping based on Touzi decomposition using quad-polarimetric ALOS PALSAR data[J]. Frontiers of Earth Science, 2020: 1–12. doi: 10.1007/s11707-019-0779-3
|
信号载频 | 叶片数 | 叶片长度 | 旋转速率 | 发射站方位角 | 接收站方位角 | 发射站仰角 | 接收站仰角 | SNR |
658 MHz | 3 | 5 m | 200 rpm | 33° | 76° | 23° | 23° | –5 dB |