CN106405547A - 一种多极化sar简缩极化模式定标方法 - Google Patents

一种多极化sar简缩极化模式定标方法 Download PDF

Info

Publication number
CN106405547A
CN106405547A CN201610696435.3A CN201610696435A CN106405547A CN 106405547 A CN106405547 A CN 106405547A CN 201610696435 A CN201610696435 A CN 201610696435A CN 106405547 A CN106405547 A CN 106405547A
Authority
CN
China
Prior art keywords
polarization
calibration
delta
corner reflectors
receiving channel
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.)
Granted
Application number
CN201610696435.3A
Other languages
English (en)
Other versions
CN106405547B (zh
Inventor
胡继伟
谭小敏
党红杏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN201610696435.3A priority Critical patent/CN106405547B/zh
Publication of CN106405547A publication Critical patent/CN106405547A/zh
Application granted granted Critical
Publication of CN106405547B publication Critical patent/CN106405547B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904SAR modes
    • G01S13/9076Polarimetric features in SAR
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904SAR modes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/024Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using polarisation effects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • G01S7/4004Means for monitoring or calibrating of parts of a radar system

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提出了一种多极化SAR简缩极化模式定标方法,该方法包括:定标场简缩极化SAR数据预处理和极化失真参量定标两部分。定标场简缩极化SAR数据预处理主要完成简缩极化SAR数据的辐射校正和法拉第(Faraday)旋转校正;极化失真参量定标利用常见的四种不同类型角反射器作为参考定标器,基于定标模型推导了简缩极化模式定标方法,即简缩极化模式极化失真参量标定方法。该方法可以精确地对极化失真参量进行标定,且定标算法不依赖于角反射器的绝对RCS值,可以有效避免不同角反射器RCS不一致引起的定标误差,算法稳健性好,实用性强。

Description

一种多极化SAR简缩极化模式定标方法
技术领域
本发明涉及一种多极化SAR简缩极化模式定标方法,属于合成孔径雷达(SAR)数据处理技术领域。
背景技术
星载合成孔径雷达(SAR)作为一种重要的高分辨率微波遥感对地观测手段,近几年来受到越来越多的重视。SAR具有全天时、全天候的工作能力,且可以在多频段、多极化方式下获得目标的高分辨率图像,因此在地形测绘、海洋监测、减灾防灾、农业、林业和军事侦察等领域发挥着重要作用。
简缩极化模式(Compact Polarimetry)是多极化SAR的一种重要工作模式,通过利用双极化SAR系统来构建观测目标的伪全极化信息。本发明针对的简缩极化模式是采用发射圆极化信号,同时接收H和V线极化信号的工作方式。相对于传统双极化SAR系统,简缩极化模式不但可以获得两个正交接收极化通道(H和V)的幅度信息,还可以获得他们的相对相位信息。可以在大幅降低SAR系统复杂度、成本、重量、数据率等条件下获得接近全极化模式的系统性能,同时测绘带宽可以达到全极化模式的两倍,因而作为一种新的极化工作模式具有广阔的发展应用前景。
由于SAR系统的非理想特性,不同极化通道之间存在极化串扰和幅相不平衡误差,这些误差将导致简缩极化模式测量得到的目标极化数据失真,从而无法真实准确地反映目标极化特性和散射机理,不能满足多极化雷达实际应用需求,必须通过简缩极化定标技术来对系统极化失真参量进行标定和校正,从而获取地物目标真实的极化散射矩阵,进而进行极化相关应用模型的反演。
A.Freeman等在2008年EUSAR会议上发表了题为“较长波长简缩极化模式定标”的文章,初步提出了误差影响下简缩极化模式的定标模型和可能的定标方法思路,但未给出具体定标方法。
陈杰等在2011年IEEE期刊上发表了题为“利用混合雷达定标器进行低频段星载简缩极化模式定标”的文章,该文章提出了一种利用角反射器和有源定标器等不同定标器的组合来解算极化失真参量的简缩极化模式定标方法。该方法假设所有定标器的绝对幅相值已知,归一化后等效为理想散射矩阵然后基于简缩极化定标模型来解算极化失真参量。而实际上不同定标器的绝对幅相值是很难精确获得的,且不同定标器间绝对幅相值也会不一致,这将导致利用该方法进行实际定标时会引入误差,实用性较差。
发明内容
本发明的技术解决问题:克服现有简缩极化模式定标方法的不足,提供了一种改进的多极化SAR简缩极化模式定标方法,本方法仅利用角反射器即可完成对简缩极化模式极化失真参量的定标,设备需求简单;定标算法利用角反射器不同极化通道分量的相对比值来解算极化失真参量误差,不需要精确已知每个角反射器的绝对幅相值,且对于不同角反射器的绝对幅相值一致性没有要求,因而对角反射器加工精度要求相对较低,实用性更强。该方法能够实现高精度的多极化SAR系统极化失真参量定标,有效提高目标观测信号的准确性和数据质量,对于简缩极化模式SAR系统具有重要意义。
本发明的一种多极化SAR简缩极化模式定标方法,主要为多极化SAR简缩极化模式极化失真参量定标,极化失真参量定标的步骤如下:
(b1)在数据预处理后定标场简缩极化SAR图像中提取四种不同类型角反射器的极化散射矩阵;
(b2)基于简缩极化模式定标算法进行接收通道幅相不平衡定标;
(b3)对不同类型角反射器的极化散射矩阵进行接收通道不平衡校正;
(b4)利用经过接收通道不平衡校正的不同类型角反射器的极化散射矩阵进行接收通道极化串扰定标和发射通道极化串扰定标;
(b5)由此完成多极化SAR简缩极化模式全部极化失真参量定标,后续可以利用这些极化失真参数对简缩极化SAR图像进行校正。
本发明技术解决方案:
一种多极化SAR简缩极化模式定标方法,该方法的步骤包括:定标设备布设、定标场简缩极化SAR数据预处理和极化失真参量定标三部分;该定标方法基于布设有四种不同类型角反射器的定标场简缩极化SAR图像来进行,因而需要在定标场布设四种不同类型的角反射器,包括:三面角反射器、0°二面角反射器、22.5°二面角反射器和45°二面角反射器;四种角反射器布设在测绘带中心区域;
一种多极化SAR简缩极化模式定标方法,该方法的步骤为:
(1)布设定标设备;定标设备包括四种不同类型的角反射器,分别为三面角反射器、0°二面角反射器、22.5°二面角反射器和45°二面角反射器;四种角反射器布设在覆盖定标场的测绘带中心区域,四种角反射器互相之间的距离不小于100m;
(2)卫星过境获取步骤(1)中布设有四种角反射器的定标场的简缩极化SAR原始数据,得到的原始数据经过成像处理后,生成定标场的简缩极化SAR图像;
(3)对步骤(2)得到的定标场的简缩极化SAR图像进行辐射校正,所述的辐射校正包括天线方向图校正和距离功率衰减校正;
所述的天线方向图校正方法为:步骤(2)得到的SAR图像中的每个像素点乘以该像素点对应的天线方向图的幅度值G(θ)的倒数,即所述的天线方向图是利用SAR辐射外定标测量获得的;
所述的距离功率衰减校正方法为:步骤(2)得到的SAR图像中的每个像素点乘以其中Ri为距离向第i像素位置对应斜距,Rn为近距端斜距,c为光速,Fs为距离向采样频率,i为距离向像素位置,R0为测绘带中心对应斜距;
如果SAR系统工作频段为P频段或L频段,还需要对辐射校正后的定标场简缩极化SAR图像进行法拉第(Faraday)旋转校正,所述的法拉第(Faraday)旋转校正方法为:经过辐射校正后SAR图像的散射矩阵左侧乘以Faraday旋转角校正矩阵;所述的Faraday旋转角校正矩阵为:
其中,其中,K为常数,大小为2.365×104,B表示磁通量密度,表示入射波法线方向与地球磁场方向的夹角,θ表示入射波法线方向与卫星垂直向下方向的夹角,表示磁场相关参数按照400km高度计算,TEC表示总电子浓度;
(4)在步骤(3)经过辐射校正后的定标场简缩极化SAR图像中提取四种不同类型角反射器的极化散射矩阵,其中三面角反射器的极化散射矩阵为0°二面角反射器的极化散射矩阵为22.5°二面角反射器的极化散射矩阵为45°二面角反射器的极化散射矩阵为其中,S1RH表示三面角反射器RH通道图像对应值,S2RH表示0°二面角反射器RH通道图像对应值,S3RH表示22.5°二面角反射器RH通道图像对应值,S4RH表示45°二面角反射器RH通道图像对应值;
S1RV表示三面角反射器RV通道图像对应值,S2RV表示0°二面角反射器RV通道图像对应值,S3RV表示22.5°二面角反射器RV通道图像对应值,S4RV表示45°二面角反射器RV通道图像对应值;
(5)利用步骤(4)提取的四种角反射器的极化散射矩阵进行接收通道幅相不平衡定标;
所述的接收通道幅相不平衡定标方法为:
则所述的接收通道幅相不平衡f为:
其中
A=P-2K3N
B=Q+jNK3(K1+K2)
C=-jM(K1+K2)+(K1-K2)N
M=K2K4+j(K1+K2)(K4-K3)/2-1
N=j+(K1-K2)(K3-K4)/2-jK2K4
P=(K3-K4)[2+K3(K1-K2)]
Q=M[2+K3(K1-K2)]-j(K1+K2)(K3-K4)
(6)利用步骤(5)获得的接收通道幅相不平衡f来对步骤(4)获得的四种不同类型角反射器的极化散射矩阵进行接收通道不平衡校正;校正方法为:对步骤(4)获得的四种不同类型角反射器的极化散射矩阵 的左侧分别乘以矩阵
(7)利用步骤(6)获得的经过接收通道不平衡校正后的四种不同类型角反射器的极化散射矩阵进行接收通道极化串扰定标和发射通道极化串扰定标;
所述的接收通道极化串扰包括接收通道H极化到V极化的串扰δ1和接收通道V极化到H极化的串扰δ2
所述δ1和δ2的定标方法为:
所述的发射通道极化串扰δc的定标方法为:
δc=(1-X)/2
其中,
(8)由此完成多极化SAR简缩极化模式全部极化失真参量定标,后续可以利用这些极化失真参数对简缩极化SAR图像进行校正,从而可以有效提高目标观测信号的准确性和数据质量。
本发明与现有技术相比的有益效果是:
(1)本发明提出的多极化SAR简缩极化模式定标方法仅利用角反射器即可完成对简缩极化模式极化失真参量的定标,设备需求简单,成本低,且角反射器重量小,便于外场布设。
(2)相比于现在定标方法,本发明提出的简缩极化模式定标算法利用角反射器不同极化通道分量的相对比值来解算极化失真参量误差,因而不需要精确已知每个角反射器的绝对幅相值。
(3)本发明提出的简缩极化模式定标方法对于不同角反射器的绝对幅相值一致性没有要求,因而对角反射器加工精度要求相对较低,实用性更强。
(4)本发明提出的简缩极化模式定标方法仅利用角反射器作为定标器,成本低,可以利用多组角反射器在不同观测区域进行定标,从而抑制随意误差的影响,提高极化失真参量定标精度。
(5)本发明提出的简缩极化模式定标方法可以有效克服上述定标方法的不足,且设备需求简单,实用性强,精度高。
(6)本发明提出了一种多极化SAR简缩极化模式定标方法,该方法包括:定标场简缩极化SAR数据预处理和极化失真参量定标两部分。定标场简缩极化SAR数据预处理主要完成简缩极化SAR数据的辐射校正和法拉第(Faraday)旋转校正;极化失真参量定标利用常见的不同类型角反射器作为参考定标器,基于定标模型推导了简缩极化模式定标方法,即简缩极化模式极化失真参量标定方法。该方法可以精确地对极化失真参量进行标定,且定标算法不依赖于角反射器的绝对RCS值,可以有效避免不同角反射器RCS不一致引起的定标误差,算法稳健性好,实用性强。
附图说明
图1是本发明方法数据处理流程图;
图2是本发明实施例中布设的三角形反射器;
图3是本发明实施例中布设的二面角反射器;
图4是本发明实施例中RADARSAT-2定标场SAR图像;
表1是本发明实施例中基于本发明定标方法定标精度。
具体实施方式
如图1所示,一种多极化SAR简缩极化模式定标方法,该方法的步骤为:
(1)布设定标设备;定标设备包括四种不同类型的角反射器,分别为三面角反射器、0°二面角反射器、22.5°二面角反射器和45°二面角反射器;四种角反射器布设在覆盖定标场的测绘带中心区域,四种角反射器互相之间的距离不小于100m;
(2)卫星过境获取步骤(1)中布设有四种角反射器的定标场的简缩极化SAR原始数据,得到的原始数据经过成像处理后,生成定标场的简缩极化SAR图像;
(3)对步骤(2)得到的定标场的简缩极化SAR图像进行辐射校正,所述的辐射校正包括天线方向图校正和距离功率衰减校正;
所述的天线方向图校正方法为:步骤(2)得到的SAR图像中的每个像素点乘以该像素点对应的天线方向图的幅度值G(θ)的倒数,即所述的天线方向图是利用SAR辐射外定标测量获得的;
所述的距离功率衰减校正方法为:步骤(2)得到的SAR图像中的每个像素点乘以其中Ri为距离向第i像素位置对应斜距,Rn为近距端斜距,c为光速,Fs为距离向采样频率,i为距离向像素位置,R0为测绘带中心对应斜距;
如果SAR系统工作频段为P频段或L频段,还需要对辐射校正后的定标场简缩极化SAR图像进行法拉第(Faraday)旋转校正,所述的法拉第(Faraday)旋转校正方法为:经过辐射校正后SAR图像的散射矩阵左侧乘以Faraday旋转角校正矩阵;所述的Faraday旋转角校正矩阵为:
其中,其中,K为常数,大小为2.365×104,B表示磁通量密度,表示入射波法线方向与地球磁场方向的夹角,θ表示入射波法线方向与卫星垂直向下方向的夹角,表示磁场相关参数按照400km高度计算,TEC表示总电子浓度;
(4)在步骤(3)经过辐射校正后的定标场简缩极化SAR图像中提取四种不同类型角反射器的极化散射矩阵,其中三面角反射器的极化散射矩阵为0°二面角反射器的极化散射矩阵为22.5°二面角反射器的极化散射矩阵为45°二面角反射器的极化散射矩阵为其中,S1RH表示三面角反射器RH通道图像对应值,S2RH表示0°二面角反射器RH通道图像对应值,S3RH表示22.5°二面角反射器RH通道图像对应值,S4RH表示45°二面角反射器RH通道图像对应值;
S1RV表示三面角反射器RV通道图像对应值,S2RV表示0°二面角反射器RV通道图像对应值,S3RV表示22.5°二面角反射器RV通道图像对应值,S4RV表示45°二面角反射器RV通道图像对应值;
(5)利用步骤(4)提取的四种角反射器的极化散射矩阵进行接收通道幅相不平衡定标;
所述的接收通道幅相不平衡定标方法为:
则所述的接收通道幅相不平衡f为:
其中
A=P-2K3N
B=Q+jNK3(K1+K2)
C=-jM(K1+K2)+(K1-K2)N
M=K2K4+j(K1+K2)(K4-K3)/2-1
N=j+(K1-K2)(K3-K4)/2-jK2K4
P=(K3-K4)[2+K3(K1-K2)]
Q=M[2+K3(K1-K2)]-j(K1+K2)(K3-K4)
(6)利用步骤(5)获得的接收通道幅相不平衡f来对步骤(4)获得的四种不同类型角反射器的极化散射矩阵进行接收通道不平衡校正;校正方法为:对步骤(4)获得的四种不同类型角反射器的极化散射矩阵 的左侧分别乘以矩阵
(7)利用步骤(6)获得的经过接收通道不平衡校正后的四种不同类型角反射器的极化散射矩阵进行接收通道极化串扰定标和发射通道极化串扰定标;
所述的接收通道极化串扰包括接收通道H极化到V极化的串扰δ1和接收通道V极化到H极化的串扰δ2
所述δ1和δ2的定标方法为:
所述的发射通道极化串扰δc的定标方法为:
δc=(1-X)/2
其中,
(8)由此完成多极化SAR简缩极化模式全部极化失真参量定标,后续可以利用这些极化失真参数对简缩极化SAR图像进行校正,从而可以有效提高目标观测信号的准确性和数据质量。
下面本发明结合附图和以一个具体实施例进一步说明本发明的工作原理、工作过程以及预期可以达到的效果。
极化失真参量定标的步骤如下:
(b1)在数据预处理后定标场简缩极化SAR图像中提取四种不同类型角反射器的极化散射矩阵;
四种不同类型的角反射器分别为:三面角反射器、0°二面角反射器、22.5°二面角反射器和45°二面角反射器。简缩极化模式下四种角反射器的归一化理想极化散射矩阵分别为:
(b2)基于简缩极化模式定标算法进行接收通道幅相不平衡定标;
具体定标算法如下:
存在极化失真参量误差时,简缩极化模式下目标观测信号模型可以表示为
其中,为复系数,表示绝对幅相值;δ12为接收通道极化串扰;f为接收通道幅相不平衡;δc为发射通道极化串扰;NH和NV表示H和V通道噪声分量。
为了便于计算,忽略绝对幅度和相位项和噪声项,同时假设极化串扰δ<<1,则误差影响下四种角反射器对应散射矩阵分别为:
则经过复杂求解过程可以得到:
其中
A=P-2K3N
B=Q+jNK3(K1+K2)
C=-jM(K1+K2)+(K1-K2)N
M=K2K4+j(K1+K2)(K4-K3)/2-1
N=j+(K1-K2)(K3-K4)/2-jK2K4
P=(K3-K4)[2+K3(K1-K2)]
Q=M[2+K3(K1-K2)]-j(K1+K2)(K3-K4)
由此可以获得接收通道幅相不平衡误差f。
(b3)对不同类型角反射器的极化散射矩阵进行接收通道不平衡校正;
校正方法如下:
其中,分别表示经过接收通道幅相不平衡误差校正后不同类型角反射器极化散射矩阵。
(b4)利用经过接收通道不平衡校正的不同类型角反射器的极化散射矩阵进行接收通道极化串扰定标和发射通道极化串扰定标;
定标算法如下:
δ1=[K1+K2+j(K2-K12]/2
δc=(1-X)/2
其中,
(b5)由此完成多极化SAR简缩极化模式全部极化失真参量定标,后续可以利用这些极化失真参数对简缩极化SAR图像进行校正。
实施例
如图2和图3所示,首先将四种不同类型角反射器布设在定标场(本实施例定标场选在机场跑道边的草地上),然后获取定标场的SAR图像,如图4所示。
本实施例采用的定标场简缩极化SAR图像是基于购买的RADARSAT-2全极化SAR图像模拟仿真生成,利用RADARSAT-2全极化图像可以来计算简缩极化模式下RH和RV通道SAR图像,然后加入极化失真参量误差来模拟误差影响下简缩极化SAR图像数据。由于购买的RADARSAT-2全极化SAR图像已经经过辐射校正,同时RADARSAT-2工作在C频段,因而定标方法中的定标场简缩极化SAR数据预处理步骤(3)中的部分内容省略。
极化失真参量定标的步骤如下:
(b1)在数据预处理后定标场简缩极化SAR图像中提取四种不同类型角反射器的极化散射矩阵,分别为:
(b2)基于简缩极化模式定标算法进行接收通道幅相不平衡定标,可以获得接收通道幅相不平衡误差f=1.0094+0.3032j。
(b3)对不同类型角反射器的极化散射矩阵进行接收通道不平衡校正,校正方法如下:
其中,分别表示经过接收通道幅相不平衡误差校正后不同类型角反射器极化散射矩阵。
(b4)利用经过接收通道不平衡校正的不同类型角反射器的极化散射矩阵进行接收通道极化串扰定标和发射通道极化串扰定标,计算结果如下:
δ2=-0.0456+0.0332j
δ1=0.0414+0.0381j
δc=-0.0357-0.0021j
(b5)由此完成多极化SAR简缩极化模式全部极化失真参量定标,然后可以利用这些极化失真参数对简缩极化SAR图像进行校正。通过比较校正后标准角反射器极化散射矩阵与理想标准角反射器极化散射矩阵的误差来对本发明提出的多极化SAR简缩极化模式定标方法精度及有效性进行评估和说明,定标结果如表1所示。
表1
极化通道 理想值 定标后 定标精度
RH 1 1.0799+0.0804j 0.67dB∠1.65°
RV -j 0.4086-1.0657j 0.59dB∠1.54°
由定标结果可以看出,本发明提出的定标方法可以获得很高的定标精度,经过定标后目标散射矩阵幅度误差约0.6dB,相位误差小于1.7°,从而验证了本发明提出的多极化SAR简缩极化模式定标方法的有效性和可行性,可以应用于多极化SAR简缩极化模式定标。
本发明可以用于任意频段星载SAR简缩极化模式定标处理。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。

Claims (6)

1.一种多极化SAR简缩极化模式定标方法,其特征在于该方法的步骤包括:
(1)布设定标设备;定标设备包括四种不同类型的角反射器,分别为三面角反射器、0°二面角反射器、22.5°二面角反射器和45°二面角反射器;四种角反射器布设在覆盖定标场的测绘带中心区域,四种角反射器互相之间的距离不小于100m;
(2)卫星过境获取步骤(1)中布设有四种角反射器的定标场的简缩极化SAR图像;
(3)对步骤(2)得到的定标场的简缩极化SAR图像进行辐射校正,所述的辐射校正包括天线方向图校正和距离功率衰减校正;
所述的天线方向图校正方法为:步骤(2)得到的SAR图像中的每个像素点乘以该像素点对应的天线方向图的幅度值G(θ)的倒数,即所述的天线方向图是利用SAR辐射外定标测量获得的;
所述的距离功率衰减校正方法为:步骤(2)得到的SAR图像中的每个像素点乘以其中Ri为距离向第i像素位置对应斜距,Rn为近距端斜距,c为光速,Fs为距离向采样频率,i为距离向像素位置,R0为测绘带中心对应斜距;
(4)在步骤(3)经过辐射校正后的定标场简缩极化SAR图像中提取四种不同类型角反射器的极化散射矩阵,其中三面角反射器的极化散射矩阵为0°二面角反射器的极化散射矩阵为22.5°二面角反射器的极化散射矩阵为45°二面角反射器的极化散射矩阵为其中,S1RH表示三面角反射器RH通道图像对应值,S2RH表示0°二面角反射器RH通道图像对应值,S3RH表示22.5°二面角反射器RH通道图像对应值,S4RH表示45°二面角反射器RH通道图像对应值;
S1RV表示三面角反射器RV通道图像对应值,S2RV表示0°二面角反射器RV通道图像对应值,S3RV表示22.5°二面角反射器RV通道图像对应值,S4RV表示45°二面角反射器RV通道图像对应值;
(5)利用步骤(4)提取的四种角反射器的极化散射矩阵进行接收通道幅相不平衡定标;
所述的接收通道幅相不平衡定标方法为:
则所述的接收通道幅相不平衡f为:
f = - B + B 2 - 4 A C 2 A
其中
A=P-2K3N
B=Q+jNK3(K1+K2)
C=-jM(K1+K2)+(K1-K2)N
M=K2K4+j(K1+K2)(K4-K3)/2-1
N=j+(K1-K2)(K3-K4)/2-jK2K4
P=(K3-K4)[2+K3(K1-K2)]
Q=M[2+K3(K1-K2)]-j(K1+K2)(K3-K4)
(6)利用步骤(5)获得的接收通道幅相不平衡f来对步骤(4)获得的四种不同类型角反射器的极化散射矩阵进行接收通道不平衡校正;校正方法为:对步骤(4)获得的四种不同类型角反射器的极化散射矩阵 的左侧分别乘以矩阵
(7)利用步骤(6)获得的经过接收通道不平衡校正后的四种不同类型角反射器的极化散射矩阵进行接收通道极化串扰定标和发射通道极化串扰定标;
所述的接收通道极化串扰包括接收通道H极化到V极化的串扰δ1和接收通道V极化到H极化的串扰δ2
所述δ1和δ2的定标方法为:
δ1=[K1+K2+j(K2-K12]/2,
所述的发射通道极化串扰δc的定标方法为:
δc=(1-X)/2
其中,
(8)根据步骤(5)获得的接收通道幅相不平衡f定标结果和步骤(7)获得的接收通道极化串扰δ1,δ2定标结果和发射通道极化串扰δc定标结果,对简缩极化SAR图像进行校正。
2.根据权利要求1所述的一种多极化SAR简缩极化模式定标方法,其特征在于:所述的步骤(3)中,如果SAR系统工作频段为P频段或L频段,还需要对辐射校正后的定标场简缩极化SAR图像进行法拉第(Faraday)旋转校正。
3.根据权利要求2所述的一种多极化SAR简缩极化模式定标方法,其特征在于:所述的法拉第(Faraday)旋转校正方法为:经过辐射校正后SAR图像的散射矩阵左侧乘以Faraday旋转角校正矩阵;所述的Faraday旋转角校正矩阵为:
c o s Ω s i n Ω - s i n Ω c o s Ω - 1
其中,TEC,其中,K为常数,大小为2.365×104,B表示磁通量密度,表示入射波法线方向与地球磁场方向的夹角,θ表示入射波法线方向与卫星垂直向下方向的夹角,表示磁场相关参数按照400km高度计算,TEC表示总电子浓度。
4.根据权利要求1所述的一种多极化SAR简缩极化模式定标方法,其特征在于:所述的步骤(1)中,三面角反射器、0°二面角反射器、22.5°二面角反射器和45°二面角反射器的理想极化散射矩阵分别为:
S 1 R H S 1 R V = 1 - j , S 2 R H S 2 R V = - 1 - j , S 3 R H S 3 R V = - 1 - j 1 - j , S 4 R H S 4 R V = - j 1 .
5.根据权利要求1所述的一种多极化SAR简缩极化模式定标方法,其特征在于:所述的步骤(5)中,接收通道幅相不平衡定标的获得方法如下:
存在极化失真参量误差时,简缩极化模式下目标观测信号模型表示为:
M R H M R V ≅ A ( r , θ ) e j φ α 2 1 δ 2 δ 1 f S h h S h v S v h S v v [ 1 - j + δ c 1 j ] + N H N V
其中,为复系数,表示绝对幅相值;δ12为接收通道极化串扰;f为接收通道幅相不平衡;δc为发射通道极化串扰;NH和NV表示H和V通道噪声分量;
为了便于计算忽略绝对幅度和相位项和噪声项,同时假设串扰δ1<<1,δ2<<1,δc<<1,则误差影响下四种角反射器目标对应散射矩阵分别为:
S 1 R H S 1 R V = ( 1 + δ c ) 1 - j δ 2 δ 1 - j f ( 1 - 2 δ c )
S 2 R H S 2 R V = - ( 1 + δ c ) 1 + j δ 2 δ 1 + j f ( 1 - 2 δ c )
S 3 R H S 3 R V = ( 1 + δ c ) - 1 - j ( 1 - 2 δ c ) + ( 1 - j ) δ 2 ( - 1 - j ) δ 1 + f - j f ( 1 - 2 δ c )
S 4 R H S 4 R V = ( 1 + δ c ) δ 2 - j ( 1 - 2 δ c ) f - jδ 1
X=1-2δc,则经过求解过程得到:
f = - B + B 2 - 4 A C 2 A
其中
A=P-2K3N
B=Q+jNK3(K1+K2)
C=-jM(K1+K2)+(K1-K2)N
M=K2K4+j(K1+K2)(K4-K3)/2-1
N=j+(K1-K2)(K3-K4)/2-jK2K4
P=(K3-K4)[2+K3(K1-K2)]
Q=M[2+K3(K1-K2)]-j(K1+K2)(K3-K4)
由此可以获得接收通道幅相不平衡f。
6.根据权利要求1所述的一种多极化SAR简缩极化模式定标方法,其特征在于:所述的步骤(6)中,不同类型角反射器的极化散射矩阵进行接收通道不平衡校正方法如下:
S i R H f S i R V f = 1 0 0 1 / f S i R H S i R V , i = 1 , 2 , 3 , 4
其中,分别表示经过接收通道幅相不平衡误差校正后不同类型角反射器极化散射矩阵。
CN201610696435.3A 2016-08-19 2016-08-19 一种多极化sar简缩极化模式定标方法 Active CN106405547B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610696435.3A CN106405547B (zh) 2016-08-19 2016-08-19 一种多极化sar简缩极化模式定标方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610696435.3A CN106405547B (zh) 2016-08-19 2016-08-19 一种多极化sar简缩极化模式定标方法

Publications (2)

Publication Number Publication Date
CN106405547A true CN106405547A (zh) 2017-02-15
CN106405547B CN106405547B (zh) 2018-11-23

Family

ID=58004511

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610696435.3A Active CN106405547B (zh) 2016-08-19 2016-08-19 一种多极化sar简缩极化模式定标方法

Country Status (1)

Country Link
CN (1) CN106405547B (zh)

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106950551A (zh) * 2017-03-31 2017-07-14 中国科学院电子学研究所 基于辐射稳定特性挖掘的sar辐射定标方法
CN107229933A (zh) * 2017-05-11 2017-10-03 西安电子科技大学 自适应体散射模型的freeman/特征值分解方法
CN107907880A (zh) * 2017-09-30 2018-04-13 西安空间无线电技术研究所 一种星载简缩极化sar模糊度分析系统、方法及存储器
CN108007476A (zh) * 2017-11-20 2018-05-08 中国科学院空间应用工程与技术中心 一种天基干涉成像雷达高度计的干涉定标方法及系统
CN108037504A (zh) * 2017-11-28 2018-05-15 中国科学院遥感与数字地球研究所 一种基于三个点目标的全极化合成孔径雷达快速极化校正的方法
CN108051790A (zh) * 2017-12-13 2018-05-18 厦门大学 一种用于混合紧缩极化探地雷达系统的极化校准方法
CN108469614A (zh) * 2018-03-27 2018-08-31 武汉大学 基于未定标极化雷达影像的角反射器探测方法
CN108931762A (zh) * 2018-09-04 2018-12-04 中国安全生产科学研究院 边坡位移监测雷达精度标定装置
CN109388897A (zh) * 2018-10-24 2019-02-26 北京环境特性研究所 一种利用旋转二面角修正rcs测量中角误差的方法
CN109459736A (zh) * 2019-01-04 2019-03-12 北京环境特性研究所 一种雷达靶标设计方法和装置
CN109633576A (zh) * 2018-10-30 2019-04-16 上海无线电设备研究所 一种通道噪声下的极化sar高精度定标算法
CN110261853A (zh) * 2019-06-12 2019-09-20 中国测绘科学研究院 一种极化合成孔径雷达的定标处理方法及装置
CN110703251A (zh) * 2019-11-07 2020-01-17 中国科学院电子学研究所 一种矫正通道误差的方法、装置、设备和存储介质
CN111103572A (zh) * 2019-12-25 2020-05-05 中国科学院遥感与数字地球研究所 基于分布目标的星载sar极化定标方法及设备
CN111781568A (zh) * 2020-06-15 2020-10-16 南京航空航天大学 一种基于模拟定标体的极化气象雷达外定标方法
CN111948615A (zh) * 2020-06-30 2020-11-17 中国资源卫星应用中心 一种星载全极化sar数据的极化定标方法及装置
CN112068086A (zh) * 2020-10-17 2020-12-11 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于外定标试验数据的岸基多通道雷达幅相校正方法
CN113945929A (zh) * 2021-12-20 2022-01-18 中国科学院空天信息创新研究院 一种未定标全极化sar无旋性区域选择方法
CN115128564A (zh) * 2022-09-02 2022-09-30 中国科学院空天信息创新研究院 对月球轨道合成孔径雷达进行极化定标的方法及装置
CN117214843A (zh) * 2023-11-07 2023-12-12 中国科学院空天信息创新研究院 一种基于定标器组的全极化sar系统性能评估方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090102704A1 (en) * 2007-09-20 2009-04-23 Takashi Fujimura Synthetic aperture radar, compact polarimetric sar processing method and program
CN102183743A (zh) * 2011-03-04 2011-09-14 北京航空航天大学 一种长波长星载ctlr模式简缩极化sar的定标方法
CN103091666A (zh) * 2011-11-07 2013-05-08 中国科学院电子学研究所 非理想定标器条件下的机载p波段极化sar定标方法
CN103941245A (zh) * 2014-05-09 2014-07-23 中国科学院电子学研究所 简缩极化干涉数据的Freeman-Durden目标分解方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090102704A1 (en) * 2007-09-20 2009-04-23 Takashi Fujimura Synthetic aperture radar, compact polarimetric sar processing method and program
CN102183743A (zh) * 2011-03-04 2011-09-14 北京航空航天大学 一种长波长星载ctlr模式简缩极化sar的定标方法
CN103091666A (zh) * 2011-11-07 2013-05-08 中国科学院电子学研究所 非理想定标器条件下的机载p波段极化sar定标方法
CN103941245A (zh) * 2014-05-09 2014-07-23 中国科学院电子学研究所 简缩极化干涉数据的Freeman-Durden目标分解方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
HONG TAN等: "Correction of Transmit Crosstalk in Reconstruction of Quad-Pol Data From Compact Polarimetry Data", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 *
JIE CHEN等: "Calibration of Spaceborne CTLR Compact Polarimetric Low-Frequency SAR Using Mixed Radar Calibrators", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 *
李亮等: "高分辨率条件下延时有源定标器SAR成像质量分析与校正", 《电子与信息学报》 *
谭洪等: "原始数据压缩对简缩极化SAR极化信息的影响", 《系统工程与电子技术》 *
陈琳等: "单发双收SAR系统通用极化定标方法", 《雷达学报》 *

Cited By (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106950551B (zh) * 2017-03-31 2019-12-31 中国科学院电子学研究所 基于辐射稳定特性挖掘的sar辐射定标方法
CN106950551A (zh) * 2017-03-31 2017-07-14 中国科学院电子学研究所 基于辐射稳定特性挖掘的sar辐射定标方法
CN107229933B (zh) * 2017-05-11 2020-08-04 西安电子科技大学 自适应体散射模型的freeman/特征值分解方法
CN107229933A (zh) * 2017-05-11 2017-10-03 西安电子科技大学 自适应体散射模型的freeman/特征值分解方法
CN107907880A (zh) * 2017-09-30 2018-04-13 西安空间无线电技术研究所 一种星载简缩极化sar模糊度分析系统、方法及存储器
CN107907880B (zh) * 2017-09-30 2019-08-09 西安空间无线电技术研究所 一种星载简缩极化sar模糊度分析系统、方法及存储器
CN108007476A (zh) * 2017-11-20 2018-05-08 中国科学院空间应用工程与技术中心 一种天基干涉成像雷达高度计的干涉定标方法及系统
CN108007476B (zh) * 2017-11-20 2020-05-01 中国科学院空间应用工程与技术中心 一种天基干涉成像雷达高度计的干涉定标方法及系统
CN108037504A (zh) * 2017-11-28 2018-05-15 中国科学院遥感与数字地球研究所 一种基于三个点目标的全极化合成孔径雷达快速极化校正的方法
CN108051790A (zh) * 2017-12-13 2018-05-18 厦门大学 一种用于混合紧缩极化探地雷达系统的极化校准方法
CN108469614A (zh) * 2018-03-27 2018-08-31 武汉大学 基于未定标极化雷达影像的角反射器探测方法
CN108469614B (zh) * 2018-03-27 2021-09-14 武汉大学 基于未定标极化雷达影像的角反射器探测方法
CN108931762A (zh) * 2018-09-04 2018-12-04 中国安全生产科学研究院 边坡位移监测雷达精度标定装置
CN109388897A (zh) * 2018-10-24 2019-02-26 北京环境特性研究所 一种利用旋转二面角修正rcs测量中角误差的方法
CN109633576A (zh) * 2018-10-30 2019-04-16 上海无线电设备研究所 一种通道噪声下的极化sar高精度定标算法
CN109633576B (zh) * 2018-10-30 2020-09-08 上海无线电设备研究所 一种通道噪声下的极化sar高精度定标算法
CN109459736A (zh) * 2019-01-04 2019-03-12 北京环境特性研究所 一种雷达靶标设计方法和装置
CN109459736B (zh) * 2019-01-04 2020-10-13 北京环境特性研究所 一种雷达靶标设计方法和装置
CN110261853B (zh) * 2019-06-12 2021-01-05 中国测绘科学研究院 一种极化合成孔径雷达的定标处理方法及装置
CN110261853A (zh) * 2019-06-12 2019-09-20 中国测绘科学研究院 一种极化合成孔径雷达的定标处理方法及装置
CN110703251B (zh) * 2019-11-07 2022-02-18 中国科学院电子学研究所 一种矫正通道误差的方法、装置、设备和存储介质
CN110703251A (zh) * 2019-11-07 2020-01-17 中国科学院电子学研究所 一种矫正通道误差的方法、装置、设备和存储介质
CN111103572A (zh) * 2019-12-25 2020-05-05 中国科学院遥感与数字地球研究所 基于分布目标的星载sar极化定标方法及设备
CN111103572B (zh) * 2019-12-25 2023-11-07 中国科学院遥感与数字地球研究所 基于分布目标的星载sar极化定标方法及设备
CN111781568A (zh) * 2020-06-15 2020-10-16 南京航空航天大学 一种基于模拟定标体的极化气象雷达外定标方法
CN111948615A (zh) * 2020-06-30 2020-11-17 中国资源卫星应用中心 一种星载全极化sar数据的极化定标方法及装置
CN111948615B (zh) * 2020-06-30 2023-10-31 中国资源卫星应用中心 一种星载全极化sar数据的极化定标方法及装置
CN112068086A (zh) * 2020-10-17 2020-12-11 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于外定标试验数据的岸基多通道雷达幅相校正方法
CN112068086B (zh) * 2020-10-17 2022-03-01 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种基于外定标试验数据的岸基多通道雷达幅相校正方法
CN113945929A (zh) * 2021-12-20 2022-01-18 中国科学院空天信息创新研究院 一种未定标全极化sar无旋性区域选择方法
CN113945929B (zh) * 2021-12-20 2022-03-25 中国科学院空天信息创新研究院 一种未定标全极化sar无旋性区域选择方法
CN115128564A (zh) * 2022-09-02 2022-09-30 中国科学院空天信息创新研究院 对月球轨道合成孔径雷达进行极化定标的方法及装置
CN117214843A (zh) * 2023-11-07 2023-12-12 中国科学院空天信息创新研究院 一种基于定标器组的全极化sar系统性能评估方法

Also Published As

Publication number Publication date
CN106405547B (zh) 2018-11-23

Similar Documents

Publication Publication Date Title
CN106405547A (zh) 一种多极化sar简缩极化模式定标方法
van Zyl et al. The effect of topography on SAR calibration
CN103091666B (zh) 非理想定标器条件下的机载p波段极化sar定标方法
CN104898100B (zh) 基于弱布设角反射器的机载sar辐射外定标处理方法
JP6349937B2 (ja) 変動検出装置、変動検出方法および変動検出用プログラム
Pierdicca et al. SAVERS: A simulator of GNSS reflections from bare and vegetated soils
Kim et al. Correcting distortion of polarimetric SAR data induced by ionospheric scintillation
CN111948615B (zh) 一种星载全极化sar数据的极化定标方法及装置
CN102183743B (zh) 一种长波长星载ctlr模式简缩极化sar的定标方法
CN103235301A (zh) 基于复数域平差理论的POLInSAR植被高度反演方法
CN110412573A (zh) 极化sar图像分解方法及存储介质
CN108375770B (zh) 一种基于旋转对称性的全极化sar数据的极化定标方法
CN102253374B (zh) 一种长波长星载π/4模式简缩极化SAR的定标方法
CN103645466B (zh) 一种基于平台姿态时变性补偿的极化定标方法
Schuler et al. Terrain topography measurement using multipass polarimetric synthetic aperture radar data
CN110133657A (zh) 基于PolInSAR地形效应补偿的植被参数反演方法
CN106526555A (zh) 基于分布目标的全极化sar隔离度评价方法
CN112346027A (zh) 用于确定合成孔径雷达图像的散射特性的方法和系统
CN111983575A (zh) 一种有源和无源融合定标方法及装置
Touzi et al. The SSCM for ship characterization using polarimetric SAR
Wang et al. In-Orbit External Calibration Method for Synthetic Aperture Radiometer
CN106815837A (zh) 基于超复数的极化sar图像多通道联合匹配方法
Fischer et al. Polarization measurements of Saturn Electrostatic Discharges with Cassini/RPWS below a frequency of 2 MHz
Xing et al. Experimental research of unsupervised Cameron/maximum-likelihood classification method for fully polarimetric synthetic aperture radar data
Touzi et al. On the use of symmetric scatterers for calibration and validation of PALSAR polarimetric modes

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