Loading [MathJax]/jax/element/mml/optable/BasicLatin.js
Li Yuqian, Yi Jianxin, Wan Xianrong, Liu Yuqi, Zhan Weijie. Helicopter Rotor Parameter Estimation Method for Passive Radar[J]. Journal of Radars, 2018, 7(3): 313-319. doi: 10.12000/JR17125
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

Imaging Simulation and Damage Assessment Feature Analysis of Ku Band Polarized SAR of Buildings

DOI: 10.12000/JR20061
Funds:  The National Key R&D Program of China (2016YFB0502504), The National Natural Science Foundation of China (41671359), The “Practical Training Plan” Project for Cross Training of High Level Talents in Beijing Colleges and Universities (No. 17)
More Information
  • Corresponding author: ZHANG Fengli, zhangfl@aircas.ac.cn
  • Received Date: 2020-05-13
  • Rev Recd Date: 2020-06-24
  • Available Online: 2020-07-01
  • Publish Date: 2020-06-01
  • Building damage assessment is important in disaster emergency monitoring. In recent years, with the increase of multi-polarization capability of Synthetic Aperture Radar (SAR), Polarimetric Synthetic Aperture Radar (PolSAR) provides more possibilities for building damage assessment, and the polarization-characteristic-based building damage assessment method has gradually become the focus of research. However, because of the limitations of data acquisition in PolSAR, current research mainly focuses on the L, C, X, and other limited bands. To obtain an in depth understanding of the polarization characteristics of damaged buildings in SAR images and develop the application of the polarization characteristics of damaged buildings to other bands, this study conducted a simulation experiment of Ku band polarized SAR of buildings, and performed damage assessment feature analysis using the SAR image polarization decomposition method. In this study, a scale model of real materials was built and the “microwave characteristic measurement and simulation imaging scientific experiment platform” was used to conduct SAR simulation imaging of the target buildings. The Ku band polarized SAR images before and after building damage were obtained. Then, the polarization scattering characteristics of buildings before and after damage were analyzed using various common polarization decomposition methods such as H/A/α decomposition, Yamaguchi decomposition and Touzi decomposition. Results show that the disoriented volume scattering component and the proportion of the disoriented secondary scattering component obtained by the Yamaguchi decomposition and the αs1 component obtained by the Touzi decomposition have good indicative significance for building damage assessment in the Ku band. Compared with the X band measurement results, the Ku band is more sensitive to building damage assessment, which has important implications for future radar remote sensing applications.

     

  • 目标的振动、转动等微动产生的微多普勒效应包含了目标的结构和运动信息,常用于目标的分类和识别[13]。目前,基于外辐射源雷达微多普勒效应目标分类和识别的研究还处于起步状态。外辐射源雷达是一种利用非合作照射源进行目标探测和分类识别的新体制雷达系统,其自身不辐射电磁能量,具有节约频谱资源,隐蔽性好,设备规模小,易于部署和组网等特点[46]。在微多普勒效应目标分类和识别方面,外辐射源雷达表现出得天独厚的优势:(1)收发分置可实现空间分集,有效避免探测盲区。(2)第三方辐射源多为连续波,长时间相干积累可记录多个连续的回波闪烁,同时有利于提高对低雷达散射截面积(Radar Cross-Section, RCS)微动目标的探测与分类识别能力。(3)对微多普勒特征的提取不要求高距离分辨率,参数估计不受第三方辐射源带宽的限制[7,8]

    针对微多普勒效应参数估计问题。文献[9,10]中依据微动目标正弦特征曲线,利用 Hough变换,在参数域中进行多维搜索提取出微动曲线进行参数估计。文献[11,12]通过正交匹配追踪(Orthogonal Matching Pursuit, OMP)算法进行稀疏逼近实现了微动目标参量的估计。上述方法均具有较好的鲁棒性,但由于估计参量维数较高导致计算量巨大。文献[13]利用微动目标在时频域的周期性,采用循环相关系数方法,实现了目标微动周期的估计,但信号周期较长时计算量急剧增加。文献[14]计算了信号的高阶矩函数,通过检测在不同时延下,高阶虚函数部分傅里叶变换累计结果的峰值位置,快速获得目标的旋转速率,相比于图像处理方法和OMP分解方法,计算复杂度较小,但抗噪性能差。而外辐射源雷达所利用的第三方辐射源多为连续波信号,其发射波形不可控,信号能量主要覆盖地面,杂波环境复杂且对空中目标增益低,利用长时间相干积累来提高处理增益会带来数据量巨大的挑战。上述因素决定了外辐射源雷达参数估计方法需要有良好的抗噪性能且计算量要小。

    直升机旋翼旋转时对雷达信号产生周期性调制,当叶片发生镜面反射时,旋翼回波出现峰值,即回波闪烁。闪烁信号在时频图像中表现为一定宽度的频率带,且闪烁时间、闪烁间隔与直升机旋翼微动参数密切相关。针对外辐射源雷达参数估计问题,本文结合上述时频域中闪烁信号的特点,通过时频分析和正交匹配追踪算法实现了直升机旋翼微动参数的估计。本文首先给出了外辐射源雷达直升机旋翼微动信号模型,其次介绍了如何在时频图中提取出闪烁信号参数及正交匹配追踪算法对直升机旋翼微动参数的估计,最后仿真和实测证明了本文方法的有效性。

    直升机旋转叶片与外辐射源雷达的位置关系如图1所示。以直升机旋转叶片的中心点为原点 o ,旋转叶片平面为 xy 面, x 轴平行于发射站与接收站所在直线,建立空间坐标系 (x,y,z) 。直升机相对于发射站和接收站距离为 rT, rR ,方位角为 γ, α ,仰角为 βT, βR (cosβTcosβR=cosβ) 。叶片上某一散射点 p 到原点 o 距离为 lP ,方位角为 φt

    图  1  外辐射源雷达直升机旋翼回波模型
    Figure  1.  Model of helicopter rotors echo for passive radar

    假设直升机平动得到补偿。在 t 时刻,从发射站经散射点 p 到接收站的距离为:

    rP(t)=||RT RP(t)||+||RRRP(t)|| (1)

    其中 RT, RR, RP(t) 分别为发射站、接收站、散射点 p 在坐标系 xyz 中的位置矢量。

    参考文献[1]中单基地直升机建模,将叶片看作线模型,外辐射源雷达直升机旋翼回波可表示为:

    s(t)=Lexp{j2πλ(rR+rT)}Nk=1sinc{ϕk(t)}exp{jϕk(t)} (2)

    其中,

    ϕk(t)=4πλL2cosβcos(αγ2)cos(φk(t)) (3)
    φk(t)=2πfrt+φ0+(k1)2π/Nα+γ2 (4)

    fr 为叶片转速, L 为叶片长度, N 为叶片数量,整数 k (0<kN) 表示第 k 个叶片, φ0 为叶片初相, λ 为照射源信号波长。

    由式(3)得第 k 个叶片引起的瞬时多普勒频移为:

    fk(t)=2πfrLλcosβcos(αγ2)sin(φk(t)) (5)

    由式(2)可知时域信号幅值受 sinc 函数调制,结合式(3)知当 φk(t) 满足式(6)时, ϕk(t)=0 ,时域信号幅值最大,此刻即时域闪烁。

    φk(t)=±π2+2πn (6)

    由式(2)知连续两个闪烁之间的时间间隔为:

    Δt={12Nfr,N1Nfr, N (7)

    直升机旋翼回波的微多普勒呈非线性变化,通过对目标回波信号进行时频分析能够揭示信号频率的时变特性。短时傅里叶变化(Short-Time Fourier Transform, STFT)计算简单,且不产生交叉项。对直升机旋翼回波信号 s(t) 进行STFT到时频域

    TF(t,f)=s(τ)w(τt)ej2πfτdτ (8)

    其中, w(t) 为窗函数。旋翼微多普勒效应特征曲线为正弦曲线,对应时域闪烁出现的时刻出现垂直于时间横轴的频率带,即时频域“闪烁”[15]

    图2为直升机旋翼回波的时频图。当直升机旋翼的叶片数为奇数时(图2(a)),时频域中正负多普勒“闪烁”交替出现;若旋翼叶片数为偶数(图2(b)),则是同时出现。

    图  2  旋翼回波信号时频分析
    Figure  2.  Time-frequency analysis of rotors echo

    设时频域中正频率“闪烁”发生的时间为 t0 ,由式(5)和式(6)知 t0 满足:

    φk(t0)=π2+2πn(n) (9)

    由式(4)和式(9)得第 k 个叶片初相与叶片数量的关系:

    φ0={πt0Δt1N2π(k1)1N+φ1, N2πt0Δt1N2π(k1)1N+φ1,N (10)

    其中

    φ1=α+γ2+2πn+π2(0φ0<2π) (11)

    由于时频图像中闪烁信号频率带垂直于时间横轴,对正频率轴数据幅值进行累加计算,并判断累加后数据局部峰值点,可得到时频域中正频率“闪烁”发生的时间。同样,对负频率轴数据幅值进行累加计算得到时频域中负频率“闪烁”发生的时间。相应的也可得到闪烁间隔。

    由式(7)知,闪烁间隔与旋翼转速、叶片数量密切相关。由式(10)知,闪烁发生的时间与叶片初相、叶片数量、整数 k 密切相关。因此,可根据得到的闪烁间隔,用叶片数量表示出旋翼转速。根据得到的闪烁时间,用叶片数量、整数 k 表示出第 k 个叶片初相。

    由式(2)知时域回波信号可分解为:

    s(t)=Mm=1cmgm(t;Λ)=Dα (12)

    其中, gm 为第 m 个原子, D 为以原子为列张成的字典矩阵 D=[g1 g2 g3···gM] CNt×M, M 为原子个数, Nt 为时间 t 离散后的取值个数, Λ 为要估计的参量, cm 为原子系数, αCM 为系数矢量,是稀疏的。可转化最优 l0 范数问题进行稀疏向量求解。OMP常用于求解此类问题,通过构建字典矩阵,不断选定与信号最匹配的原子进行稀疏逼近[16]。OMP将字典矩阵中原子正交化保证了迭代的最优性。

    由式(2)知直升机旋翼回波信号由参数 (fr,L,φ0,N,k) 确定。利用叶片数量 N 、整数 k (0<kN) 与旋翼转速和叶片初相的关系式(7)和式(10),时域回波可转化为参数 (L,N,k) 来表示。设时间采样点数 Nt ,目标回波为 Nt×1 的矩阵。确定待估参数的取值范围并离散化,叶片长度取值: L(L1,···,Lr,···,LNL) ,叶片数量 N 的可能取值为: N(N1,···,Np,···,NNN) ,整数 k 的取值为 k(1,···,kq,···,kNq) (kNqNNN)

    由OMP算法原理可知,字典中的原子可按照待分解信号的内在特性来构造[16]。根据微动目标的时域回波表达式(2),第 m 个原子可表示为:

    a(m)=sinc(ϕ(Lr,Np,kq))exp{jϕ(Lr,Np,kq)} (13)

    其中

    m=rpq (14)

    并对原子集里的每个原子进行能量归一化:

    a(m)a(m)/a(m)F (15)

    其中, F 表示矩阵的F范数。

    将5参量 (fr,L,φ0,N,k) 的估计转换为3参量 (L,N,k) 估计, NL, NN, Nk 分别为 L, N, k 的取值个数,由于常见直升机主旋翼叶片数量为:3片、5片、7片(奇数),2片、4片、8片(偶数), NN, Nk 较小,降低字典维数为: NL×Nk×NN ,可达到降低计算量的目的。

    直升机旋翼参数估计具体步骤如下:

    步骤1 对直升机旋翼信号进行短时傅里叶变换,得到时频图像 TF(t,f)

    步骤2 对时频图中正频率轴数据幅值进行累加计算,并判断累加后数据局部峰值点,对应时频域正频率“闪烁”发生的时间。同样,对负频率轴数据幅值进行累加计算得到时频域中负频率“闪烁”发生的时间。

    步骤3 根据步骤2中正负频率“闪烁”发生的时间,判别时频域中正负多普勒“闪烁”是否交替出现。若是,则旋翼叶片数为奇数,否则,旋翼叶片数为偶数。

    步骤4 读取某一正频率闪烁发生的时间 t0 及闪烁间隔 Δt 。依据式(7)用叶片数量 N 表示出旋翼转速,依据式(10)和式(11)用叶片数量 N 及整数 k 表示出第 k 个叶片初相。

    步骤5 确定 (L,N,k) 的取值范围并离散化: L(L1,···,Lr,···,LNL) , N(N1,···,Np,···,NNN) , k(1,···,kq,···,kNq) (kNqNNN) 。利用步骤4中表示出的旋翼转速及初相,依据式(13)和式(15)构建字典矩阵。

    步骤6 利用OMP算法寻找叶片数量,叶片长度的最优值,代入式(7)计算出旋翼转速,代入式(10)和式(11)计算出叶片初相。

    结合上述模型对直升机旋翼回波信号进行仿真,仿真参数设置如表1所示。

    表  1  外辐射源雷达直升机旋翼回波模型仿真参数
    Table  1.  Simulation parameters of helicopter rotor echo model for passive radar
    信号载频 叶片数 叶片长度 旋转速率 发射站方位角 接收站方位角 发射站仰角 接收站仰角 SNR
    658 MHz 3 5 m 200 rpm 33° 76° 23° 23° –5 dB
    下载: 导出CSV 
    | 显示表格

    图3(a)显示了信号的联合时频域特征,可看出闪烁信号及噪声严重影响直升机旋翼微多普勒特征曲线的检测,使微多普勒特征曲线提取困难。

    分别对时频图像中正负频率轴数据幅值进行累加计算,得到时频域中正负多普勒“闪烁”时间,如图3(b)所示,图中正负多普勒“闪烁”等间隔交替出现,则旋翼叶片数为奇数。读取闪烁信号时间间隔为0.05 s,根据式(7)表示出旋翼转速为:

    fr=10/N (16)

    读取某一正频率闪烁信号对应时刻为0.066 s(此处选择了图3(b)中的第1个正频率闪烁信号),根据式(10)和式(11)表示出第 k 个叶片初相为:

    φ0=4.14/N6.28×(k1)/N+2.52 (17)

    图3(c)为利用OMP方法对 (L,N,k) 的估计结果,得到叶片数为3片,图中给出了其对应的切面图,3叶片长度分别4.99 m, 5.00 m, 4.98 m,均值4.99 m,与理论基本一致,代入式(17)得3叶片初相分别为1.14 rad, 3.24 rad, 5.34 rad,代入式(16)得旋翼转速为200 rpm,与理论值一致,本文方法准确实现了直升机旋翼参数估计。

    图  3  本文方法参数估计结果
    Figure  3.  Parameter estimation by this article method

    图4为利用常规Hough变换,通过微多普勒曲线 f=fmax 检测对参数 ({f\!_r},{\varphi _0},L) 的估计结果。其中 {f_{\max }} 为最大频移。

    {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受不确定原理的限制,时频图像中时频分辨率受限使参数空间中的局部峰值点扩展范围较大,只能大致估计局部峰值点的位置,估计结果精度较低。

    图  4  常规Hough变换参数估计结果
    Figure  4.  Parameter estimation by traditional Hough transform

    设待处理的时频图像大小为 {N_t} \times {N_f} 像素, {N_f} \approx {N_t} ,利用常规的Hough变换对微多普勒曲线 f = {f_{\max }}\sin(2{{π}} {f\!_r}t{\rm{ + }}{\varphi _0}) 进行检测,参数 ({f\!_r},{\varphi _0},L) 分别量化为 {N_{f_{r}}} , {N_{{\varphi _0}}} , {N_L} 份。乘法次数可近似表示为: 2{N_{f_{r}}}{N_{{\varphi _0}}}{N_L}{N_t}^{\!2}

    直接使用OMP进行参数 ({f\!_r},{\varphi _0},L) 估计时,设迭代次数为K,乘法次数近似表示为: K{N_{f\!_{r}}}{N_{{\varphi _0}}} {N_L}{N_t}^{\!2}

    本文方法计算量集中在OMP阶段,根据提取的时频域中的闪烁时间,依据式(7)和式(10),最终转化为对参数 (L,N,k) 的估计,常见直升机的叶片数只有若干个取值,且由3.1节方法可判断出叶片数量的奇偶性, {N_N}{N_k} 远小于 {N_{f_{r}}}{N_{{\varphi _0}}} 。本文方法乘法次数近似表示为: K{N_N}{N_k}{N_L}{N_t}^{\!2}

    在对直升机旋翼微动参数估计时,一般 {N_N}{N_k} 取值量级为101~102,迭代次数K的取值量级为 {10^0} {\text ~} {10^1} ,当初相 {\varphi _0} 的估计精度为7°时, 2{N_{{\varphi _0}}} 取值量级为102,当转速 {f_r} 的估计精度为10 rpm时, {N_{f_{r}}} 的取值量级为 {10^1} ,在乘法次数上,常规 Hough变换参数估计方法为本文方法 {10^0}{\text ~} {10^2} 倍,当进一步提高 {f_r}, {\varphi _0} 的估计精度时,算法之间的计算量差距将进一步变大。本文在相同的配置环境下,利用Matlab仿真平台,常规Hough变换方法运行时长13006 s,而本文方法运行总时长只有145 s。

    武汉大学电波传播实验室对EC_120B直升机进行了微多普勒效应探究外场实验,EC_120B直升机主旋翼3叶片,叶片长度5 m,额定转速406 rpm,实验中以武汉龟山电视塔数字电视信号为照射源,信号中心频率为658 MHz,带宽8 MHz,接收站位于武汉大学电波传播实验室楼顶,距离发射站7.56 km,实验场景如图5所示。本组实测数据相干积累时间0.8 s,可近似认为目标在这段时间位置不变,直升机旋翼转速为常量。

    图  5  实验场景图
    Figure  5.  Experimental scene map

    图6(a)为去除目标主体影响后,对直升机旋翼回波信号进行短时傅里叶变换后的时频图像。可以观察到闪烁信号,但微多普勒特征曲线已观察不到。分别对时频图像中正负频率轴数据幅值进行累加计算,得到时频域中正负多普勒“闪烁”时间,如图6(b)所示,图中正负多普勒“闪烁”等间隔交替出现,则旋翼叶片数为奇数。

    图  6  本文方法参数估计结果
    Figure  6.  Parameter estimation by this article method

    读取闪烁信号时间间隔为26.2 ms,根据式(7)用 N 表示出旋翼转速。读取某一正频率闪烁时间对应时间为0.25 s,利用式(10)和式(11)表示出第 k 个叶片初相,图6(c)为利用OMP方法对 (L,N,k) 的估计结果,得到叶片数量为3片,3叶片长度分别为4.93 m, 5.00 m, 4.66 m,均值4.86 m,存在较小的误差,与仰角,方位角估计不精确有关,由式(7)知旋翼转速均为382 rpm,符合实际情况。

    本文根据外辐射源雷达直升机旋翼微动信号模型,充分利用时频域中闪烁信号特征和微动信号内在特性进行了参数估计。通过时频分析和正交匹配追踪算法,估计出了旋翼转速、叶片长度、叶片数量和初相。同时开展了外场实验。仿真数据和实测数据处理都表明本文方法对外辐射源雷达直升机旋翼参数估计的可行性。

  • [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
  • 加载中

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(10)  / Tables(2)

    Article views(2884) PDF downloads(242) Cited by()
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    • 表  1  外辐射源雷达直升机旋翼回波模型仿真参数
      Table  1.  Simulation parameters of helicopter rotor echo model for passive radar
      信号载频 叶片数 叶片长度 旋转速率 发射站方位角 接收站方位角 发射站仰角 接收站仰角 SNR
      658 MHz 3 5 m 200 rpm 33° 76° 23° 23° –5 dB
      下载: 导出CSV 
      | 显示表格