CN105445708B - 一种极化合成孔径雷达的定标方法 - Google Patents
一种极化合成孔径雷达的定标方法 Download PDFInfo
- Publication number
- CN105445708B CN105445708B CN201510765328.7A CN201510765328A CN105445708B CN 105445708 B CN105445708 B CN 105445708B CN 201510765328 A CN201510765328 A CN 201510765328A CN 105445708 B CN105445708 B CN 105445708B
- Authority
- CN
- China
- Prior art keywords
- mtd
- mrow
- msub
- mtr
- polarization
- 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
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000010287 polarization Effects 0.000 claims abstract description 118
- 238000005388 cross polarization Methods 0.000 claims abstract description 47
- 238000012937 correction Methods 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims description 90
- 239000000654 additive Substances 0.000 claims description 7
- 230000000996 additive effect Effects 0.000 claims description 7
- 230000014509 gene expression Effects 0.000 claims description 5
- 238000004891 communication Methods 0.000 abstract description 2
- 238000000711 polarimetry Methods 0.000 abstract 2
- 238000004422 calculation algorithm Methods 0.000 description 61
- 238000004088 simulation Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002955 isolation Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012360 testing method 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)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于雷达通信技术领域,公开了一种极化合成孔径雷达的定标方法。该方法包括:获取全极化散射回波数据;建立所述极化合成孔径雷达的极化定标模型;计算所述极化定标模型中交叉极化通道的不平衡参数和所述极化通道的串扰参数;根据所述交叉极化通道的不平衡参数和极化通道的串扰参数对所述全极化散射回波数据进行初步校正,得到初步校正后的回波数据;计算所述接收通道的不平衡参数,根据所述接收通道的不平衡参数对所述初步校正后的回波数据进行再次校正,得到全极化合成孔径雷达的定标数据。本发明能够在交叉通道噪声功率不一致以及信噪比较低的情况下仍能保持良好的参数估计性能。
Description
技术领域
本发明涉及雷达通信技术领域,具体涉及一种极化合成孔径雷达的定标方法。
背景技术
极化合成孔径雷达遥感已经在地物分类、表面粗糙度和土壤湿度估计以及农业应用等方面得到广泛的研究和应用。近年来,国内外已成功研制了众多的机载和星载极化SAR***。目前仍在执行飞行试验的机载极化***主要有德国F-SAR***、美国UAVSAR***及日本的PiSAR-2***;星载极化***主要有德国的TanDEM-X***以及日本的ALOS-2***等。利用这些***,国内外已获取了大量的极化数据。在实际应用时,所有的极化分析与应用均是建立在各极化通道间的相对幅度和相对相位已进行精确标定的前提下,因而在利用极化SAR数据进行各种应用时,首要需要校正由***性能不理想引起的极化数据的畸变。
目前,国内外众多学者已提出了大量的极化定标算法及其改进算法,主要包括使用点目标、使用分布目标及混合使用点目标和分布目标等三类定标算法,其中混合使用分布目标和点目标的定标算法对极化***和分布目标进行了一定的假设,从而简化了定标过程,因而该类算法在实际中得到了广泛的应用。
由于目前的极化SAR***可以获得很高的极化隔离度,并能通过精确的内定标技术进一步提高极化数据的精度,因而在2006年,Ainsworth提出了一种假设条件更加宽泛的定标算法,其仅要求目标满足散射互易性,从而大大扩展了可用定标区域的选取范围,其已被成功应用于F-SAR***和UAVSAR***极化数据的标定。
但在实际应用中,我们发现Ainsworth算法在估计不平衡参数的过程中并没有考虑交叉极化通道噪声的影响,因而在交叉极化通道噪声不可忽略时,参数估计精度将会下降,从而导致极化数据的畸变无法被正确校正。
发明内容
针对上述现有技术的缺点,本发明的实施例提供一种极化合成孔径雷达的定标方法,利用分布目标估计交叉极化通道的不平衡参数,进而对极化SAR数据进行标定。本发明充分考虑了交叉极化通道的噪声对交叉极化通道的不平衡参数估计的影响,其在交叉通道噪声功率不一致以及信噪比较低的情况下仍能保持良好的参数估计性能。
为达到上述目的,本发明的实施例采用如下技术方案:
一种极化合成孔径雷达的定标方法,所述极化合成孔径雷达的各极化通道的噪声与信号不相关且各极化通道的噪声之间也不相关,所述方法包括如下步骤:
步骤1,极化合成孔径雷达获取全极化散射回波数据;
步骤2,建立所述极化合成孔径雷达的极化定标模型,所述极化定标模型中至少包含交叉极化通道的不平衡参数、极化通道的串扰参数和接收通道的不平衡参数;
步骤3,计算所述极化定标模型中交叉极化通道的不平衡参数和所述极化通道的串扰参数;
步骤4,根据所述交叉极化通道的不平衡参数和极化通道的串扰参数对所述全极化散射回波数据进行初步校正,得到初步校正后的回波数据;
步骤5,计算所述接收通道的不平衡参数,根据所述接收通道的不平衡参数对所述初步校正后的回波数据进行再次校正,得到全极化合成孔径雷达的定标数据。
本发明的特点和进一步的改进为:
(1)步骤1中,所述全极化散射回波数据,具体包括:
所述合成孔径雷达的HH通道回波数据,所述合成孔径雷达的HV通道回波数据,所述合成孔径雷达的VV通道回波数据,所述合成孔径雷达的VH通道回波数据。
(2)步骤2中,所述建立所述极化合成孔径雷达的极化定标模型,具体包括如下子步骤:
(2a)所述极化合成孔径雷达观测的极化散射矩阵O与目标的真实散射矩阵S间的关系采用极化合成孔径雷达***误差模型进行描述,其中,所述观测的极化散射矩阵O由所述全极化散射回波数据构造;
其中,极化合成孔径雷达***误差模型为:
O=RST+N
R为接收失真矩阵,T为发射失真矩阵,N为***加性噪声矩阵,且所述接收失真矩阵R、所述发射失真矩阵T和所述***加性噪声矩阵N都为2×2矩阵;
(2b)将所述极化合成孔径雷达***误差模型O=RST+N写为如下矩阵形式:
其中,矩阵元素的下标h表示水平极化状态,矩阵元素的下标v表示垂直极化状态,若将矩阵元素用xij表示,则j表示发射波的极化方式,i表示接收波的极化方式;
(2c)将所述极化合成孔径雷达***误差模型的矩阵形式矢量化后,得到如下所述极化合成孔径雷达***误差模型的矢量形式:
其中,k为接收通道的不平衡参数,α为交叉极化通道的不平衡参数,u,v,w,z为极化通道的串扰参数,k=rhh/rvv,α=rvvthh/rhhtvv,u=rvh/rhh,v=tvh/tvv,w=rhv/rvv,z=thv/thh;
(2d)定义以下变量:
o=[Ohh Ohv Ovh Ovv]T,s=[Shh Shv Svh Svv]T
n=[nhh nhv nvh nvv]T,A=diag(α,1,α,1)
K=diag(k2,k,k,1),Y=rvvtvv
则所述极化合成孔径雷达***误差模型的矢量形式表示为:o=YMAKs+n,其中[·]T表示转置操作;
(2e)根据上述子步骤(2a)~(2d),得到所述极化合成孔径雷达的极化定标模型为:
Co=|Y|2MAKCsKHAHMH+Cn
其中,Co为观测协方差矩阵,Cs为真实协方差矩阵,Co=<ooH>,Cs=<ssH>,Cn=<nnH>,[·]H表示共轭转置操作,<·>表示求集合平均,σ11=ShhShh *,σ41=SvvShh *,σ44=SvvSvv *,A=ShvShh *,B=ShvSvv *,β=ShvShv *。
(3)步骤3中,所述计算所述极化定标模型中交叉极化通道的不平衡参数、所述极化通道的串扰参数,具体包括如下子步骤:
(3a)定义C=Co-Cn=MADAHMH;
其中,τ11=|Y|2|k|4σ11,τ41=|Y|2k*σ41,τ44=|Y|2σ44,γ=|Yk|2β,G=k*|Yk|2A,H=|Y|2kB;
(3b)根据下式估计交叉极化通道的不平衡参数α的值;
其中,Cxy表示矩阵C的第x行第y列的元素,m为VH通道和HV通道的噪声功率比;
(3c)计算矩阵A=diag(α,1,α,1),得到矩阵L=A-1CA-H,并表示矩阵G=(L21+L31)/2,H=(L24+L34)/2;
(3d)采用下式估计极化通道的串扰参数u,v,w,z;
其中,κ=(u,v,z,w)T,且有X=(L21-G,L31-G,L24-H,L34-H)T,
(3e)利用交叉极化通道的不平衡参数α的值和极化通道的串扰参数u,v,w,z的值对矩阵D进行估计,将矩阵D作为矩阵C的估计值;
(3f)令m1=m/|α|2,m=m1,重复执行步骤(3b)~(3e),直到迭代次数达到设定的次数,或者所有串扰参数的变化量低于设定阀值,所述串扰参数的变化量是指连续两次迭代得到的串扰参数的差值;
(3g)根据每次迭代的中间值(ui,vi,wi,zi,αi)确定交叉极化通道的不平衡参数和极化通道的串扰参数:
u=u1+u2+u3+···
v=v1+α1·v2+α1·α2·v3+···
w=w1+w2+w3+···
z=z1+α1 -1z2+(α1·α2)-1z3+…
α=α1·α2·α3···
(4)步骤4中,所述初步校正后的回波数据具体为:
初步校正后的回波数据具体为:
初步校正后的回波数据表示为:
其中,A=diag(α,1,α,1),o=[Ohh Ohv Ovh Ovv]T,
(5)步骤5,计算所述接收通道的不平衡参数,根据所述接收通道的不平衡参数对所述初步校正后的回波数据进行再次校正,得到全极化合成孔径雷达的定标数据,具体包括:
采用下式计算所述接收通道的不平衡参数k:
根据所述接收通道的不平衡参数k对所述初步校正后的回波数据进行再次校正,得到全极化合成孔径雷达的定标数据s0具体为:
其中,表示矩阵的第x行第y列的元素。
本发明与现有的技术相比具有以下优点:由于实际中各极化通道的噪声功率可能并不相等,尤其对于H和V极化接收通道分置的***,如PALSAR***等,因而本发明中的改进定标算法在处理中并不要求交叉极化通道的噪声功率相等。本发明假定VH通道和HV通道的噪声功率比可利用符合一阶Bragg散射模型的光滑表面进行估计,一般平静的水面是最佳目标。为了消除交叉极化通道噪声对参数估计精度的影响,本发明修改了Anisworth算法基本步骤中参数的估计公式。不仅考虑了交叉极化通道存在噪声时的情况,而且还考虑了交叉极化通道噪声不一致的情况。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的极化SAR的定标方法流程示意图;
图2为原始Ainsworth算法α幅度估计值的曲线示意图;
图3为原始Ainsworth算法α相位估计值的曲线示意图;
图4为原始Ainsworth算法α幅度误差估计值的曲线示意图;
图5为改进Ainsworth算法α幅度估计值的曲线示意图;
图6为改进Ainsworth算法α相位估计值的曲线示意图;
图7为改进Ainsworth算法α幅度误差估计值的曲线示意图;
图8为α幅度估计值随信噪比的变化的曲线示意图;
图9为α幅度估计误差随信噪比的变化的曲线示意图;
图10为基于德国Oberpfaffenhofen地区全极化数据的α幅度估计值随距离的变化曲线示意图;
图11为基于德国Oberpfaffenhofen地区全极化数据的α相位估计值随距离的变化曲线示意图;
图12为德国Oberpfaffenhofen地区全极化数据合成示意图;
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供一种极化合成孔径雷达的定标方法,如图1所示,为该定标方法的流程示意图,包括:
步骤1,极化合成孔径雷达获取全极化散射回波数据。
所述全极化散射回波数据,具体包括:
所述合成孔径雷达的HH通道回波数据,所述合成孔径雷达的HV通道回波数据,所述合成孔径雷达的VV通道回波数据,所述合成孔径雷达的VH通道回波数据。
具体的,包括如下步骤:
(1a)将极化合成孔径雷达的HH通道数据输入到***;
(1b)将极化合成孔径雷达的HV通道数据输入到***;
(1c)将极化合成孔径雷达的VV通道数据输入到***;
(1d)将极化合成孔径雷达的VH通道数据输入到***。
步骤2,建立所述极化合成孔径雷达的极化定标模型,所述极化定标模型中至少包含交叉极化通道的不平衡参数、极化通道的串扰参数和接收通道的不平衡参数。
极化合成孔径雷达定标算法的有效性与极化合成孔径雷达***模型的有效性紧密相关。通常认为极化合成孔径雷达***模型是建立在***性能稳定这个基础上的,尤其要求发射机、接收机和天线性能稳定。
所述建立所述极化合成孔径雷达的极化定标模型,具体包括如下子步骤:
(2a)所述极化合成孔径雷达观测的极化散射矩阵O与目标的真实散射矩阵S间的关系采用极化合成孔径雷达***误差模型进行描述,其中,所述观测的极化散射矩阵O由所述全极化散射回波数据构造;
其中,极化合成孔径雷达***误差模型为:
O=RST+N (1)
R为接收失真矩阵,T为发射失真矩阵,N为***加性噪声矩阵,且所述接收失真矩阵R、所述发射失真矩阵T和所述***加性噪声矩阵N都为2×2矩阵;
(2b)将所述极化合成孔径雷达***误差模型O=RST+N写为如下矩阵形式:
其中,矩阵元素的下标h表示水平极化状态,矩阵元素的下标v表示垂直极化状态,若将矩阵元素用xij表示,则j表示发射波的极化方式,i表示接收波的极化方式;
(2c)将所述极化合成孔径雷达***误差模型的矩阵形式矢量化后,得到如下所述极化合成孔径雷达***误差模型的矢量形式:
其中,k为接收通道的不平衡参数,α为交叉极化通道的不平衡参数,u,v,w,z为极化通道的串扰参数。
各参数的定义如下:
k=rhh/rvv α=rvvthh/rhhtvv
u=rvh/rhh v=tvh/tvv (4)
w=rhv/rvv z=thv/thh
(2d)定义以下变量:
o=[Ohh Ohv Ovh Ovv]T,s=[Shh Shv Svh Svv]T
n=[nhh nhv nvh nvv]T,A=diag(α,1,α,1)
K=diag(k2,k,k,1),Y=rvvtvv
则所述极化合成孔径雷达***误差模型的矢量形式表示为:
o=YMAKs+n (6)
其中[·]T表示转置操作;
假设各极化通道噪声与信号不相关且各极化通道噪声间也不相关,则观测协方差矩阵和真实协方差矩阵间的关系可用下式表示
Co=|Y|2MAKCsKHAHMH+Cn (7)
(2e)根据上述子步骤(2a)~(2d),得到所述极化合成孔径雷达的极化定标模型为:
Co=|Y|2MAKCsKHAHMH+Cn
其中,Co为观测协方差矩阵,Cs为真实协方差矩阵,Co=<ooH>,Cs=<ssH>,Cn=<nnH>,[·]H表示共轭转置操作,<·>表示求集合平均。
一般假设场景中的目标满足互易性,则目标的真实协方差矩阵Cs具有以下的形式:
σ11=ShhShh *,σ41=SvvShh *,σ44=SvvSvv *,A=ShvShh *,B=ShvSvv *,β=ShvShv *。通常情况下,参数Y,k的求取和参数u,v,w,z,α的求取可以独立进行。Y和k参数的反演需要散射矩阵精确已知的先验目标。
步骤3,计算所述极化定标模型中交叉极化通道的不平衡参数和所述极化通道的串扰参数。
Ainsworth极化定标方法仅仅要求场景中地物满足电磁波传输的互易性,并不要求场景中地物满足同极化通道和交叉极化通道间不相关的假设,从而扩大了定标区域的选取范围。
所述计算所述极化定标模型中交叉极化通道的不平衡参数、所述极化通道的串扰参数,具体包括如下子步骤:
(3a)由(7)式定义得:
C=Co-Cn=MADAHMH (9)
其中,
τ11=|Y|2|k|4σ11,τ41=|Y|2k*σ41,τ44=|Y|2σ44,γ=|Yk|2β,G=k*|Yk|2A,H=|Y|2kB;可以看出D和Cs具有相同的形式。
如若定义:
L=A-1CA-H (11)
则有
其中利用(9)式~(10)式,Ainsworth算法可以通过迭代实现串扰参数u,v,w,z和通道不平衡参数α的估计。原始Ainsworth算法的基本步骤如下所示:
①根据(13)式估计参数α的值:
其中,Cxy表示矩阵C的xy项,下文中所有类似符号的含义与此处相同。
②计算矩阵A=diag(α,1,α,1),并利用(11)式计算矩阵L。
③利用(14)式估计G和H
④利用(15)~(18)式估计串扰参数u,v,w,z
其中,κ=(u,v,z,w)T,且有
X=(L21-G,L31-G,L24-H,L34-H)T (16)
⑤利用步骤①中估计得到的α值、步骤④中估计得到的串扰参数及(12)式完成矩阵D的估计。
利用上述步骤可进行循环迭代,即将每次迭代得到的矩阵D作为矩阵C的一个估计,循环步骤①~⑤。如此循环迭代,直到迭代次数达到设定的次数或串扰参数变化量低于设定阀值。最终的参数估计值可由每次迭代的中间值(ui,vi,wi,zi,αi)给出
在实际应用中,Anisworth算法通常忽略噪声项Cn,而没有像Quegan算法那样考虑交叉通道噪声对参数α估计的影响。当交叉极化通道噪声不可忽略,但在实际处理中不予考虑时,参数α的估计精度将受到影响。本发明中加性噪声被建模为零均值高斯白噪声,即
由于实际中各极化通道的噪声功率可能并不相等,尤其对于H和V极化接收通道分置的***,如PALSAR***等,因而本发明中的改进定标算法在处理中并不要求交叉极化通道的噪声功率相等。本发明假定VH通道和HV通道的噪声功率比为m=N3/N2。m值可利用符合一阶Bragg散射模型的光滑表面进行估计,一般平静的水面是最佳目标。
为了消除交叉极化通道噪声对参数α估计精度的影响,本发明将Anisworth算法基本步骤中参数α的估计公式(13)进行了修改,具体如下所述。
(3b)根据下式估计交叉极化通道的不平衡参数α的值;
其中,Cxy表示矩阵C的第x行第y列的元素,m为VH通道和HV通道的噪声功率比;
(3c)计算矩阵A=diag(α,1,α,1),得到矩阵L=A-1CA-H,并表示矩阵G=(L21+L31)/2,H=(L24+L34)/2;
(3d)采用下式估计极化通道的串扰参数u,v,w,z;
其中,κ=(u,v,z,w)T,且有X=(L21-G,L31-G,L24-H,L34-H)T,
(3e)利用交叉极化通道的不平衡参数α的值和极化通道的串扰参数u,v,w,z的值对矩阵D进行估计,将矩阵D作为矩阵C的估计值;
(3f)令m1=m/|α|2,m=m1,重复执行步骤(3b)~(3e),直到迭代次数达到设定的次数,或者所有串扰参数的变化量低于设定阀值,所述串扰参数的变化量是指连续两次迭代得到的串扰参数的差值;
式(21)在估计α时,不仅考虑了交叉极化通道存在噪声时的情况,而且还考虑了交叉极化通道噪声不一致的情况。当按照上述步骤完成第一次迭代时,由下式
可得转换后的噪声矩阵P,在忽略参数u,v,w,z的高阶项(二阶及以上)的情况下,有
P22≈N2,P23≈0,P32≈0,P33≈N3/|α|2 (23)
因而每次迭代完成后需要按照下式更新m值
m=m/|α|2 (24)
(3g)根据每次迭代的中间值(ui,vi,wi,zi,αi)确定交叉极化通道的不平衡参数和极化通道的串扰参数:
u=u1+u2+u3+···
v=v1+α1·v2+α1·α2·v3+···
w=w1+w2+w3+···
α=α1·α2·α3···
本发明改进的Anisworth定标算法的处理流程为:首先根据式(21)计算参数α的值,然后利用得到α值,按照Anisworth的步骤②,③,④和⑤进行处理,最后利用(24)式更新估计m值。通过迭代此过程,则可利用(19)式得到各参数的最终估计值。
步骤4,根据所述交叉极化通道的不平衡参数和极化通道的串扰参数对所述全极化散射回波数据进行初步校正,得到初步校正后的回波数据。
所述初步校正后的回波数据具体为:
初步校正后的回波数据表示为:
其中,A=diag(α,1,α,1),o=[Ohh Ohv Ovh Ovv]T,
如若忽略串扰参数的二阶及二阶以上的所有高阶项,则有
步骤5,计算所述接收通道的不平衡参数,根据所述接收通道的不平衡参数对所述初步校正后的回波数据进行再次校正,得到全极化合成孔径雷达的定标数据。
经过初步校正后,则仅有发射与接收通道不平衡比值参数k造成极化SAR数据的畸变。参数k的求解需要散射特性精确已知的外部散射体。此外,也可根据理论模型或经验模型来求取,但一般情况下并不可靠。若图像中存在散射矩阵精确已知的三面角目标,则可采用下式计算所述接收通道的不平衡参数k:
其中,表示矩阵的第x行第y列的元素。
忽略绝对参数Y,利用k参数完成数据校正。一般来说,根据上述方法求取的k值有π的误差,即存在模糊,但这个模糊度在实际中影响很小,它可能会把45°线极化解译成135°线极化,不过对于大多数自然地表来说,它们都是各项同性的散射体,因此不会有任何差异。
根据所述接收通道的不平衡参数k对所述初步校正后的回波数据进行再次校正,得到全极化合成孔径雷达的定标数据s0具体为:
其中,表示矩阵的第x行第y列的元素。
下面结合仿真实验对本发明的效果作进一步的描述。
为了验证本发明中改进定标算法在精确估计参数α方面的有效性,我们分别利用仿真数据和德国Oberpfaffenhofen地区全极化数据对改进算法进行了试验分析。
由于实际数据的失真矩阵不可能事先知道,因此本发明利用已标定好的全极化数据,通过设定失真参数值来产生仿真数据。
真实协方差矩阵Cs中的元素被设置为代表植被的典型值,即σ11=1,σ41=0.01ej0.7,σ44=1,β=0.25,A=B=0.1。在本发明中,我们主要分析交叉极化通道噪声对参数α的影响,以及改进算法对于α估计精度的提升性能。
由于目前的极化雷达***大都具有较高的极化隔离度,这意味着串扰参数均较小(幅度值小于0.1),所以串扰参数对参数α估计精度的影响可以忽略。本发明将串扰参数u,v,w,z设置为固定值,如表1所示。参数α的幅度和相位分别在0.9~1.1和0~2π范围内选取。构造噪声矩阵Cn,使其具有diag(N,N/m,N,N)的形式,其中,N的取值范围为0~0.05,m的取值范围为0.5~2。
表1 仿真时串扰参数设置值
为了分析改进算法估计参数α的性能随m值的变化情况,我们设定三个具有代表性的α值,其幅度值分别为0.9,1.0和1.1,相位均为-70.08°,并设定N=0.05。图2~图7表示的是m取不同值时参数α的估计值及估计误差,其中圆圈代表设置的α幅度值为0.9;正方形代表设置的α幅度值为1.0;三角代表设置的α幅度值为1.1。图2为利用原始的Ainsworth算法估计得到的α幅度值,图5为利用改进算法估计得到的α幅度值。从图2和图5可以看出,原始Ainsworth算法的α幅度值的估计精度受交叉极化通道噪声影响明显,而改进后的算法对噪声表现了良好的鲁棒性,α幅度值的估计精度大大提高。从图4可以看出,当交叉通道噪声功率不一致时,原始算法对于α值的估计将会产生较大的误差,达到10-1数量级,而由图7可以看出,利用改进算法可以大大提高α的估计精度,误差下降到10-3数量级。另外,从图3和图6可以看出,原始Ainsworth算法和改进算法对于α相位值的估计精度差异不大,误差均在±1°之内。由于原始Ainsworth算法和改进算法对于α相位值的估计准确度均很高,因此下面将主要在α的幅度值估计精度方面对改进定标算法的估计性能展开讨论。
为了分析交叉极化通道信噪比对参数α估计精度的影响,本发明分别将α和m的值设置为固定值1.1和0.5,而N的值由0到0.05均匀变化。本发明分别利用原始Ainsworth算法和改进算法对参数α值进行估计。由于在上述设置下,HV通道的信噪比要低于VH通道的信噪比(参考噪声矩阵的形式),因而我们仅画出α估计值以及其估计误差随HV通道信噪比的变化趋势,如图8和图9所示。图8为利用原始Ainsworth算法和改进算法估计得到的α幅度值随信噪比的变化,图9为其估计误差随信噪比的变化。从图8和图9可以看出,当交叉极化通道信噪比较高时(此时对应低噪声水平),原始Ainsworth算法可以获得良好的α估计性能,但当信噪比降低时,估计精度明显降低。而改进算法的估计性能对信噪比的变化表现了良好的鲁棒性,低信噪比下估计误差仍然很小,因而改进算法对于参数α的估计性能要明显优于原始算法。
为了进一步验证改进定标算法的性能,本发明利用德国Oberpfaffenhofen地区全极化数据,通过设定失真参数值和加性噪声值产生了一组仿真数据。该组数据中,串扰参数设置值仍然采用表1中所示的值,α值的幅度和相位值分别为1.1和-70.08°,VH通道和HV通道的所加噪声的功率之比m取值为0.14。图10和图11为基于该数据利用原始Ainsworth算法和改进算法估计得到的参数α的幅度和相位随距离的变化。从图10和图11可以看出,基于原始Ainsworth算法和改进算法的α的幅度和相位估计值随距离向均有所波动,这是由于实测数据中噪声沿距离向不均匀造成的。从图11可以看出原始算法和改进算法在α相位值估计方面的性能基本相同,这与上面的分析结果相同。但值得注意的是,基于改进算法的α幅度估计值随距离向的波动明显比基于原始Ainsworth算法的小的多,如图10所示,这说明改进算法的稳健性要明显优于原始算法。
为了更直观地说明改进算法的性能提升,我们给出了德国Oberpfaffenhofen地区全极化数据无失真时的合成图以及分别基于原始算法和改进算法失真参数估计值校正后的合成图,如图12所示。从图中可以看出,相比于原始算法,基于改进算法估计得到的失真参数校正后的合成图更接近于无失真时的合成图,尤其在椭圆标注的区域。
实验分析结果表明,本发明中改进定标算法在交叉通道噪声功率不一致以及信噪比较低的情况下仍能保持良好的估计性能,从而可以有效保证参数的估计精度,实现极化SAR数据的精确标定,利于后续的极化分类等应用的研究。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (5)
1.一种极化合成孔径雷达的定标方法,所述极化合成孔径雷达的各极化通道的噪声与信号不相关且各极化通道的噪声之间也不相关,其特征在于,所述方法包括如下步骤:
步骤1,极化合成孔径雷达获取全极化散射回波数据;
步骤2,建立所述极化合成孔径雷达的极化定标模型,所述极化定标模型中至少包含交叉极化通道的不平衡参数、极化通道的串扰参数和接收通道的不平衡参数;
步骤3,计算所述极化定标模型中交叉极化通道的不平衡参数和所述极化通道的串扰参数;
步骤3中,所述计算所述极化定标模型中交叉极化通道的不平衡参数、所述极化通道的串扰参数,具体包括如下子步骤:
(3a)定义C=Co-Cn=MADAHMH;
其中,τ11=|Y|2|k|4σ11,τ41=|Y|2k*σ41,τ44=|Y|2σ44,γ=|Yk|2β,G=k*|Yk|2A,H=|Y|2kB;Co为观测协方差矩阵,Cs为真实协方差矩阵,Co=<ooH>,Cs=<ssH>,Cn=<nnH>,[·]H表示共轭转置操作,<·>表示求集合平均,
(3b)根据下式估计交叉极化通道的不平衡参数α的值;
<mrow>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mn>33</mn>
</msub>
<mo>/</mo>
<msub>
<mi>C</mi>
<mn>23</mn>
</msub>
<mo>,</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>C</mi>
<mn>23</mn>
<mo>*</mo>
</msubsup>
<mo>/</mo>
<msub>
<mi>C</mi>
<mn>22</mn>
</msub>
</mrow>
<mrow>
<mi>&alpha;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>|</mo>
<mo>-</mo>
<mi>m</mi>
<mo>+</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>|</mo>
<mo>-</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>4</mn>
<mi>m</mi>
<mo>|</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<mn>2</mn>
<mo>|</mo>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
<mfrac>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>|</mo>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
其中,Cxy表示矩阵C的第x行第y列的元素,m为VH通道和HV通道的噪声功率比;
(3c)计算矩阵A=diag(α,1,α,1),得到矩阵L=A-1CA-H,并表示矩阵G=(L21+L31)/2,H=(L24+L34)/2;α为交叉极化通道的不平衡参数;
(3d)采用下式估计极化通道的串扰参数u,v,w,z;
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>Re</mi>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>Im</mi>
<mo>(</mo>
<mi>X</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>Re</mi>
<mrow>
<mo>(</mo>
<mi>&xi;</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>Im</mi>
<mrow>
<mo>(</mo>
<mi>&xi;</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>Im</mi>
<mrow>
<mo>(</mo>
<mi>&xi;</mi>
<mo>+</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>Re</mi>
<mrow>
<mo>(</mo>
<mi>&xi;</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>Re</mi>
<mo>(</mo>
<mi>&kappa;</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>Im</mi>
<mo>(</mo>
<mi>&kappa;</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,κ=(u,v,z,w)T,且有X=(L21-G,L31-G,L24-H,L34-H)T,
<mrow>
<mi>&xi;</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msub>
<mi>L</mi>
<mn>41</mn>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&alpha;L</mi>
<mn>11</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mn>41</mn>
</msub>
<mo>/</mo>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msub>
<mi>L</mi>
<mn>44</mn>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&alpha;L</mi>
<mn>14</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mn>14</mn>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mn>44</mn>
</msub>
<mo>/</mo>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>&tau;</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mn>22</mn>
</msub>
<mo>/</mo>
<msup>
<mi>&alpha;</mi>
<mo>*</mo>
</msup>
</mrow>
</mtd>
<mtd>
<msub>
<mi>L</mi>
<mn>23</mn>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mn>32</mn>
</msub>
<mo>/</mo>
<msup>
<mi>&alpha;</mi>
<mo>*</mo>
</msup>
</mrow>
</mtd>
<mtd>
<msub>
<mi>L</mi>
<mn>33</mn>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&alpha;L</mi>
<mn>23</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>L</mi>
<mn>32</mn>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msup>
<mi>&alpha;</mi>
<mo>*</mo>
</msup>
<msub>
<mi>L</mi>
<mn>33</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
k为接收通道的不平衡参数,u,v,w,z为极化通道的串扰参数;
(3e)利用交叉极化通道的不平衡参数α的值和极化通道的串扰参数u,v,w,z的值对矩阵D进行估计,将矩阵D作为矩阵C的估计值;
(3f)令m1=m/|α2,m=m1,重复执行步骤(3b)~(3e),直到迭代次数达到设定的次数,或者所有串扰参数的变化量低于设定阀值,所述串扰参数的变化量是指连续两次迭代得到的串扰参数的差值;
(3g)根据每次迭代的中间值(ui,vi,wi,zi,αi)确定交叉极化通道的不平衡参数和极化通道的串扰参数:
u=u1+u2+u3+…
v=v1+α1·v2+α1·α2·v3+…
w=w1+w2+w3+…
z=z1+α1 -1z2+(α1·α2)-1z3+…
α=α1·α2·α3…;
步骤4,根据所述交叉极化通道的不平衡参数和极化通道的串扰参数对所述全极化散射回波数据进行初步校正,得到初步校正后的回波数据;
步骤5,计算所述接收通道的不平衡参数,根据所述接收通道的不平衡参数对所述初步校正后的回波数据进行再次校正,得到全极化合成孔径雷达的定标数据。
2.根据权利要求1所述的一种极化合成孔径雷达的定标方法,其特征在于,步骤1中,所述全极化散射回波数据,具体包括:
所述合成孔径雷达的HH通道回波数据,所述合成孔径雷达的HV通道回波数据,所述合成孔径雷达的VV通道回波数据,所述合成孔径雷达的VH通道回波数据。
3.根据权利要求1所述的一种极化合成孔径雷达的定标方法,其特征在于,步骤2中,所述建立所述极化合成孔径雷达的极化定标模型,具体包括如下子步骤:
(2a)所述极化合成孔径雷达观测的极化散射矩阵O与目标的真实散射矩阵S间的关系采用极化合成孔径雷达***误差模型进行描述,其中,所述观测的极化散射矩阵O由所述全极化散射回波数据构造;
其中,极化合成孔径雷达***误差模型为:
O=RST+N
R为接收失真矩阵,T为发射失真矩阵,N为***加性噪声矩阵,且所述接收失真矩阵R、所述发射失真矩阵T和所述***加性噪声矩阵N都为2×2矩阵;
(2b)将所述极化合成孔径雷达***误差模型O=RST+N写为如下矩阵形式:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>O</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>O</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>O</mi>
<mrow>
<mi>v</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>O</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>r</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>r</mi>
<mrow>
<mi>v</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>r</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>t</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>t</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>t</mi>
<mrow>
<mi>v</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>t</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>n</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>n</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>n</mi>
<mrow>
<mi>v</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>n</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,矩阵元素的下标h表示水平极化状态,矩阵元素的下标v表示垂直极化状态,若将矩阵元素用xij表示,则j表示发射波的极化方式,i表示接收波的极化方式;
(2c)将所述极化合成孔径雷达***误差模型的矩阵形式矢量化后,得到如下所述极化合成孔径雷达***误差模型的矢量形式:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>O</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>O</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>O</mi>
<mrow>
<mi>v</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>O</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
<msub>
<mi>t</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>k</mi>
<mn>2</mn>
</msup>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>v</mi>
<mi>k</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>w</mi>
<mi>k</mi>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>v</mi>
<mi>w</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>zk</mi>
<mn>2</mn>
</msup>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mi>k</mi>
</mtd>
<mtd>
<mrow>
<mi>w</mi>
<mi>z</mi>
<mi>k</mi>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mi>w</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>uk</mi>
<mn>2</mn>
</msup>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>u</mi>
<mi>v</mi>
<mi>k</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>k</mi>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mi>v</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<mi>uzk</mi>
<mn>2</mn>
</msup>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>u</mi>
<mi>k</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>z</mi>
<mi>k</mi>
<mi>&alpha;</mi>
</mrow>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>n</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>n</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>n</mi>
<mrow>
<mi>v</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>n</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,k为接收通道的不平衡参数,α为交叉极化通道的不平衡参数,u,v,w,z为极化通道的串扰参数,k=rhh/rvv,α=rvvthh/rhhtvv,u=rvh/rhh,v=tvh/tvv,w=rhv/rvv,z=thv/thh;
(2d)定义以下变量:
o=[Ohh Ohv Ovh Ovv]T,s=[Shh Shv Svh Svv]T
n=[nhh nhv nvh nvv]T,A=diag(α,1,α,1)
K=diag(k2,k,k,1),Y=rvvtvv
<mrow>
<mi>M</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>v</mi>
</mtd>
<mtd>
<mi>w</mi>
</mtd>
<mtd>
<mrow>
<mi>v</mi>
<mi>w</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>z</mi>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mi>w</mi>
<mi>z</mi>
</mrow>
</mtd>
<mtd>
<mi>w</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>u</mi>
</mtd>
<mtd>
<mrow>
<mi>u</mi>
<mi>v</mi>
</mrow>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>v</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>u</mi>
<mi>z</mi>
</mrow>
</mtd>
<mtd>
<mi>u</mi>
</mtd>
<mtd>
<mi>z</mi>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
则所述极化合成孔径雷达***误差模型的矢量形式表示为:o=YMAKs+n,其中[·]T表示转置操作;
(2e)根据上述子步骤(2a)~(2d),得到所述极化合成孔径雷达的极化定标模型为:
Co=|Y|2MAKCsKHAHMH+Cn
其中,Co为观测协方差矩阵,Cs为真实协方差矩阵,Co=<ooH>,Cs=<ssH>,Cn=<nnH>,[·]H表示共轭转置操作,<·>表示求集合平均,σ11=ShhShh *,σ41=SvvShh *,σ44=SvvSvv *,A=ShvShh *,B=ShvSvv *,β=ShvShv *。
4.根据权利要求1所述的一种极化合成孔径雷达的定标方法,其特征在于,步骤4中,所述初步校正后的回波数据具体为:
<mrow>
<mover>
<mi>s</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mi>Y</mi>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>k</mi>
<mn>2</mn>
</msup>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>kS</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>kS</mi>
<mrow>
<mi>v</mi>
<mi>h</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
</mrow>
初步校正后的回波数据表示为:
<mrow>
<mover>
<mi>s</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<msup>
<mi>A</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>o</mi>
</mrow>
其中,A=diag(α,1,α,1),o=[Ohh Ohv Ovh Ovv]T,
<mrow>
<msup>
<mi>M</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>v</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>w</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>v</mi>
<mi>w</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>z</mi>
</mrow>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mi>w</mi>
<mi>z</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>w</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>u</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>u</mi>
<mi>v</mi>
</mrow>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>v</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>u</mi>
<mi>z</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>u</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>z</mi>
</mrow>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&times;</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>u</mi>
<mi>w</mi>
</mrow>
<mo>)</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>v</mi>
<mi>z</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
k为接收通道的不平衡参数,α为交叉极化通道的不平衡参数,u,v,w,z为极化通道的串扰参数。
5.根据权利要求1所述的一种极化合成孔径雷达的定标方法,其特征在于,步骤5,计算所述接收通道的不平衡参数,根据所述接收通道的不平衡参数对所述初步校正后的回波数据进行再次校正,得到全极化合成孔径雷达的定标数据,具体包括:
采用下式计算所述接收通道的不平衡参数k:
<mrow>
<mi>k</mi>
<mo>=</mo>
<msqrt>
<mfrac>
<msub>
<mover>
<mi>s</mi>
<mo>&OverBar;</mo>
</mover>
<mn>11</mn>
</msub>
<msub>
<mover>
<mi>s</mi>
<mo>&OverBar;</mo>
</mover>
<mn>22</mn>
</msub>
</mfrac>
</msqrt>
</mrow>
3
根据所述接收通道的不平衡参数k对所述初步校正后的回波数据进行再次校正,得到全极化合成孔径雷达的定标数据s0具体为:
<mrow>
<msub>
<mi>s</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mi>Y</mi>
<mo>&CenterDot;</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<msub>
<mover>
<mi>s</mi>
<mo>&OverBar;</mo>
</mover>
<mn>11</mn>
</msub>
<msup>
<mi>k</mi>
<mn>2</mn>
</msup>
</mfrac>
</mtd>
<mtd>
<mfrac>
<msub>
<mover>
<mi>s</mi>
<mo>&OverBar;</mo>
</mover>
<mn>21</mn>
</msub>
<mi>k</mi>
</mfrac>
</mtd>
<mtd>
<mfrac>
<msub>
<mover>
<mi>s</mi>
<mo>&OverBar;</mo>
</mover>
<mn>31</mn>
</msub>
<mi>k</mi>
</mfrac>
</mtd>
<mtd>
<msub>
<mover>
<mi>s</mi>
<mo>&OverBar;</mo>
</mover>
<mn>41</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
</mrow>
其中,表示矩阵的第x行第y列的元素。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510765328.7A CN105445708B (zh) | 2015-11-11 | 2015-11-11 | 一种极化合成孔径雷达的定标方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510765328.7A CN105445708B (zh) | 2015-11-11 | 2015-11-11 | 一种极化合成孔径雷达的定标方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105445708A CN105445708A (zh) | 2016-03-30 |
CN105445708B true CN105445708B (zh) | 2017-10-13 |
Family
ID=55556127
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510765328.7A Active CN105445708B (zh) | 2015-11-11 | 2015-11-11 | 一种极化合成孔径雷达的定标方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105445708B (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106501804A (zh) * | 2016-11-25 | 2017-03-15 | 中国石油大学(华东) | 一种利用全极化sar回波数据解析海面风浪谱的方法 |
CN107179534B (zh) * | 2017-06-29 | 2020-05-01 | 北京北科天绘科技有限公司 | 一种激光雷达参数自动标定的方法、装置及激光雷达 |
CN107422324A (zh) * | 2017-08-08 | 2017-12-01 | 四川测绘地理信息局测绘技术服务中心 | 反射非对称信息的合成孔径雷达极化定标方法 |
CN108037504A (zh) * | 2017-11-28 | 2018-05-15 | 中国科学院遥感与数字地球研究所 | 一种基于三个点目标的全极化合成孔径雷达快速极化校正的方法 |
CN108051790A (zh) * | 2017-12-13 | 2018-05-18 | 厦门大学 | 一种用于混合紧缩极化探地雷达***的极化校准方法 |
CN108375770B (zh) * | 2018-03-15 | 2020-07-17 | 中国科学院遥感与数字地球研究所 | 一种基于旋转对称性的全极化sar数据的极化定标方法 |
CN108646226B (zh) * | 2018-05-10 | 2022-03-15 | 北京航空航天大学 | 用于极化散射矩阵测量的无源极化校准体 |
CN109388897A (zh) * | 2018-10-24 | 2019-02-26 | 北京环境特性研究所 | 一种利用旋转二面角修正rcs测量中角误差的方法 |
CN109633576B (zh) * | 2018-10-30 | 2020-09-08 | 上海无线电设备研究所 | 一种通道噪声下的极化sar高精度定标算法 |
CN110261853B (zh) * | 2019-06-12 | 2021-01-05 | 中国测绘科学研究院 | 一种极化合成孔径雷达的定标处理方法及装置 |
CN111103572B (zh) * | 2019-12-25 | 2023-11-07 | 中国科学院遥感与数字地球研究所 | 基于分布目标的星载sar极化定标方法及设备 |
CN111596271B (zh) * | 2020-06-01 | 2022-09-16 | 中国科学院空天信息创新研究院 | 基于有源定标器参考矩阵的合成孔径雷达极化定标方法 |
CN112330560B (zh) * | 2020-11-05 | 2024-02-20 | 中国科学院国家空间科学中心 | 一种合成孔径雷达数据图像可视化增强方法及*** |
CN112698285B (zh) * | 2020-12-24 | 2022-06-24 | 南京航空航天大学 | 一种基于Ainsworth算法的气象雷达外定标方法 |
CN113945928B (zh) * | 2021-12-20 | 2022-04-08 | 中国科学院空天信息创新研究院 | 一种基于无旋性区域的全极化sar定标方法 |
CN115047416B (zh) * | 2022-08-15 | 2022-10-25 | 中国科学院空天信息创新研究院 | 一种基于对称区域表面扭曲成分的全极化sar定标方法 |
CN115542276B (zh) * | 2022-11-28 | 2023-07-04 | 中国科学院空天信息创新研究院 | 一种圆极化***下分布式目标的选取及定标方法 |
CN117214843B (zh) * | 2023-11-07 | 2024-01-12 | 中国科学院空天信息创新研究院 | 一种基于定标器组的全极化sar***性能评估方法 |
CN117253132B (zh) * | 2023-11-14 | 2024-01-26 | 中国科学院空天信息创新研究院 | 基于分布式目标的dfsar全极化数据可用性评估方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183743A (zh) * | 2011-03-04 | 2011-09-14 | 北京航空航天大学 | 一种长波长星载ctlr模式简缩极化sar的定标方法 |
CN103645466A (zh) * | 2013-12-16 | 2014-03-19 | 中国科学院电子学研究所 | 一种基于平台姿态时变性补偿的极化定标方法 |
CN103869299A (zh) * | 2014-04-01 | 2014-06-18 | 武汉大学 | 基于自然裸土的极化合成孔径雷达定标方法 |
-
2015
- 2015-11-11 CN CN201510765328.7A patent/CN105445708B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183743A (zh) * | 2011-03-04 | 2011-09-14 | 北京航空航天大学 | 一种长波长星载ctlr模式简缩极化sar的定标方法 |
CN103645466A (zh) * | 2013-12-16 | 2014-03-19 | 中国科学院电子学研究所 | 一种基于平台姿态时变性补偿的极化定标方法 |
CN103869299A (zh) * | 2014-04-01 | 2014-06-18 | 武汉大学 | 基于自然裸土的极化合成孔径雷达定标方法 |
Non-Patent Citations (5)
Title |
---|
A Unified Algorithm for Phase and Cross-Talk Calibration of Polarimetric Data-Theory and Observations;Shaun Quegan;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;19940131;第32卷(第1期);第89-99页 * |
Calibration of Polarimetric Radar Images Using Only Image Parameters and Trihedral Corner Reflector Responses;JAKOB J.VAN ZYL;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;19900331;第28卷(第3期);第337-348页 * |
Co-polarization channel imbalance determination by the use of bare soil;Lei Shi.etc;《ISPRS Journal of Photogrammetry and Remote Sensing》;20141231;第53-67页 * |
Improvement of Polarimetric SAR Calibration Based on the Ainsworth Algorithm for Chinese Airborne PolSAR Data;Hong Zhang.etc;《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》;20130731;第10卷(第4期);第898-900页 * |
Orientation Angle Preserving A Posteriori Polarimetric SAR Calibration;Thomas L. Ainsworth.etc;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;20060430;第44卷(第4期);第994-1003页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105445708A (zh) | 2016-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105445708B (zh) | 一种极化合成孔径雷达的定标方法 | |
CN103323846B (zh) | 一种基于极化干涉合成孔径雷达的反演方法及装置 | |
CN104391279B (zh) | 基于电离层传播特性的相径扰动抑制方法 | |
CN105158759B (zh) | 基于杂波相消的hrws sar通道相位偏差校正方法 | |
CN103439708B (zh) | 基于广义散射矢量的极化InSAR干涉图估计方法 | |
CN107942331B (zh) | 基于谱分析的多通道sar***通道偏差估计方法 | |
CN103605117B (zh) | 基于干扰抑制的相控阵信号失真实时校正方法 | |
CN103323822A (zh) | 一种估计通道误差的方法及装置 | |
CN108880586B (zh) | 一种宽带弱信号增强方法与装置 | |
CN104539340B (zh) | 一种基于稀疏表示和协方差拟合的稳健波达角估计方法 | |
Zhang et al. | Improvement of polarimetric SAR calibration based on the Ainsworth algorithm for Chinese airborne PolSAR data | |
CN109633576A (zh) | 一种通道噪声下的极化sar高精度定标算法 | |
CN104181513B (zh) | 一种雷达天线阵元位置的校正方法 | |
CN104020456A (zh) | 一种基于多特显点线阵成像雷达***幅相误差校正方法 | |
CN105403871A (zh) | 一种双基地mimo雷达阵列目标角度估计和互耦误差校准方法 | |
CN113945928B (zh) | 一种基于无旋性区域的全极化sar定标方法 | |
CN106707250A (zh) | 基于互耦校正的雷达阵列自适应波束形成方法 | |
CN103760534A (zh) | 一种星载sar数据的电离层色散效应校正方法 | |
CN108375770B (zh) | 一种基于旋转对称性的全极化sar数据的极化定标方法 | |
CN104535987A (zh) | 适用于均匀圆阵列声纳***的幅相误差自校正方法 | |
CN107607915A (zh) | 基于固定地物回波的有源相控阵雷达接收通道校正方法 | |
CN102778672B (zh) | 一种应用于多极化sar的相位误差估计方法 | |
CN109725298A (zh) | 一种基于秩损求根的阵列流型误差校准和波达方向估计方法 | |
CN113030888A (zh) | 一种极化相控阵雷达测量误差的轴向校正方法 | |
CN104155653A (zh) | 一种基于特征距离子空间的sar后向投影成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |