CN106093932A - 一种扫描波束高分辨率的雷达散射计 - Google Patents
一种扫描波束高分辨率的雷达散射计 Download PDFInfo
- Publication number
- CN106093932A CN106093932A CN201610656894.9A CN201610656894A CN106093932A CN 106093932 A CN106093932 A CN 106093932A CN 201610656894 A CN201610656894 A CN 201610656894A CN 106093932 A CN106093932 A CN 106093932A
- Authority
- CN
- China
- Prior art keywords
- resolution
- formula
- distance
- scatterometer
- band
- 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
Links
- 238000000034 method Methods 0.000 claims abstract description 78
- 238000012545 processing Methods 0.000 claims abstract description 45
- 230000008569 process Effects 0.000 claims abstract description 36
- 238000005070 sampling Methods 0.000 claims abstract description 21
- 230000003321 amplification Effects 0.000 claims abstract description 4
- 238000001514 detection method Methods 0.000 claims abstract description 4
- 239000000284 extract Substances 0.000 claims abstract description 4
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 4
- 238000007906 compression Methods 0.000 claims description 10
- 230000006835 compression Effects 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 8
- 230000007246 mechanism Effects 0.000 claims description 6
- 230000005855 radiation Effects 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000005316 response function Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 5
- 210000001367 artery Anatomy 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 210000003462 vein Anatomy 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 238000010257 thawing Methods 0.000 description 2
- 208000004350 Strabismus Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 238000010297 mechanical methods and process Methods 0.000 description 1
- 230000005226 mechanical processes and functions Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
- G01S13/00—Systems 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/42—Simultaneous measurement of distance and other co-ordinates
- G01S13/422—Simultaneous measurement of distance and other co-ordinates sequential lobing, e.g. conical scan
-
- 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
- G01S13/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- 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
- G01S13/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种扫描波束高分辨率的雷达散射计,所述雷达散射计采用笔形波束圆锥扫描体制,其包括:天线单元、发射机单元、接收机单元、信号处理与控制单元、多普勒锐化处理单元和高分辨处理单元;该天线单元采用窄波束圆锥扫描体制,该发射机单元生成线性调频脉冲信号,发射该脉冲时对应一个观测方位角该接收机单元将接收到的该脉冲信号进行低噪声放大、中频接收、混频去斜坡和相位检波处理,提取出回波信号,该数字处理与控制单元记录所述发射信号对应观测方位角根据该观测方位角大小,判断所述雷达散射计的天线指向的波束照射区域;即刈幅两侧或刈幅中间;从而确定采用多普勒波束锐化方法或者高采样重建方法,得到高分辨率的后向散射系数。
Description
技术领域
本发明涉及雷达散射计的技术领域,特别涉及一种扫描波束高分辨率的雷达散射计。
背景技术
雷达散射计是一种用来测量地球表面归一化后向散射系数的雷达***,其最初的设计目标应用于大尺度海洋表面风的测量,即通过地球物理模型方程(GMF)将测得的后向散射数据反演出海面风场。散射计作为一种主动微波遥感器,具有观测刈幅宽、全球覆盖速度快、测量精度高等优点,已经成为陆地定量遥感的理想传感器。
根据天线配置和观测方位角测量实现方式,星载微波散射计分为固定扇形波束散射计、笔形波束旋转扫描散射计和扇形波束旋转扫描散射计三类。由于采用高增益的笔形波束天线,使笔形波束扫描散射计所需功率较小,从而使其逐渐成为星载雷达散射计的主流。由于雷达散射计在地面上形成分辨单元具有两维尺度,即距离向分辨率和方位向分辨率;所以分辨单元需要在两个不同的方向上实现高分辨率,即通过提高距离向和方位向的分辨率来提高空间分辨率。目前,传统的笔形波束散射计只使用了脉冲压缩技术来提高距离向分辨率,而方位向分辨率由于受天线波束宽度的限制,能达到的分辨率相对较低。随着雷达散射计观测技术的发展,兼顾宽刈幅和高分辨率要求成为观测需求和技术发展的重要趋势。
宽刈幅的实现,采用波束扫描方式是一种有效的解决方案,而围绕星下方向的圆锥扫描刈幅宽度决定于入射角外沿,能够实现宽的观测刈幅。
为了提高空间分辨率,采用调频脉冲和脉冲压缩技术,可以有效提高距离方向(俯仰方向)的分辨率,已经成为现代星载雷达提高分辨率的有效手段。对于方位方向分辨率的提高,采用多普勒波束锐化方案的综合孔径信号处理和采用超采样高分辨率重建是提高分辨率的途径。
综合孔径处理技术在距离-方位二维向上可以采用联合距离-多普勒处理算法来提高空间分辨率,这与合成孔径雷达(SAR)获得高分辨率的方法相类似。虽然,在原理上,圆锥扫描微波散射计采用的多普勒滤波技术与SAR的情形相类似;但是,两者仍具有重要差别:(1)SAR***通常采用阵列天线,以固定的方位角(一般是侧视)观测地面,存在观测刈幅窄、方位角单一的问题;而圆锥扫描微波散射计的天线足迹连续旋转扫描目标场景,能够实现宽刈幅、多方位角的观测;(2)对于圆锥扫描微波散射计,测量方位角随着测量幅宽而变化,这相当于每个测量单元的斜视角不同,导致幅宽内分辨率性能的变化。
提高散射计空间分辨率的另一种方法是高采样重建方法。众所周知,随着散射计天线扫描获得的真实孔径条带后向散射系数实际上是视场范围内地表散射强度被天线方向图加权模糊的结果,因此地表条带后向散射系数高分辨重建可等价为对天线方向图的反卷积过程。利用这种反卷积方法的前提是已知天线方向图且满足重建区域充分采样。通常,天线方向图可在发射前测得。对于传统的用于风场测量的散射计而言,重建区域的充分采样可通过覆盖这一区域多轨、多天的观测数据实现。但对于如冻融监测这种昼夜变化较大的陆地观测情况,如果采用多天的观测数据无疑会使反演精度下降。特别是随着星载微波散射计发展,逐步提高其采样频率已成为散射计发展的重要趋势。因此,通过提高散射计采样频率并结合反卷积处理算法实现条带后向散射系数高分辨重建是提高雷达散射计空间分辨率的重要方法。
发明内容
本发明的目的在于,为解决现有的合成孔径雷达观测刈幅窄,以及传统的雷达散射计空间分辨率低的观测缺陷,进而提供一种扫描波束高空间分辨率的雷达散射计;所述雷达散射计采用旋转扫描天线,其包括:天线单元、发射机单元,接收机单元、数字处理与控制单元、多普勒锐化处理单元和高分辨处理单元;所述发射机单元生成并发射线性调频脉冲信号,在发射该线性调频脉冲信号时,对应一个观测方位角所述接收机单元接收该脉冲信号并对其进行低噪声放大、中频接收、混频去斜坡和相位检波的处理,提取出回波信号,所述数字处理与控制单元记录所述发射信号对应的观测方位角当所述观测方位角为刈幅两侧时,波束照射区域为多普勒波束锐化区域,在所述多普勒锐化处理单元,采用多普勒波束锐化方法,将去载频后的回波进行方位向和距离向的两维去调频,再进行距离向和方位向脉冲压缩,得到由多个分辨单元组成的二维的高分辨率的处理数据,从而得到高分辨后向散射系数;当所述观测方位角在刈幅中间时,所述波束照射区域为高采样重建区域,在所述高分辨处理单元中,采用高采样重建方法,测得条带回波功率,并计算条带后向散射系数,再根据反卷积法,得到高分辨后向散射系数。
所述天线单元采用窄波束圆锥扫描体制,用于实现旋转笔形波束扫描;其包括反射面天线和扫描伺服机构;该扫描伺服机构用于控制天线旋转。
当所述波束照射区域为多普勒波束锐化区域时,即刈幅两侧,采用多普勒波束锐化方法,所述多普勒锐化处理单元对经过所述接收机单元和所述数字处理与控制单元的回波信号去载频,将所述去载频后的回波进行方位向和距离向的两维去调频,再进行距离向和方位向脉冲压缩,得到由多个分辨单元组成的二维的高空间分辨率的处理数据,该处理数据需要下传到星下进行处理,进一步得到观测的后向散射系数,所述后向散射系数作为观测目标特性的反演数据。
当所述波束照射区域为高采样重建区域时,即刈幅中间,采用高采样重建方法,由所述高分辨处理单元下传的地面反射的条带回波功率,计算条带后向散射系数,利用反卷积方法,进而得到高分辨后向散射系数。
所述多普勒波束锐化方法进一步包括:
(a)将所述多普勒锐化处理单元对经过所述接收机单元和所述数字处理与控制单元的回波信号进行去载频,得到雷达在一个Burst数据块内获得的基带回波数据sr(t,τ),记为(1)式,
其中,
其中,积分区域为扫描场景范围,c为光速,σ(x,r)表示地面(x,r)处的条带后向散射系数,为天线方向图,即 为天线方位角,G0为天线的最大增益,βaz为天线方位向波束宽度;t和τ分别对应表示方位向的“慢时间”和表示距离向的“快时间”,τp表示发射脉冲宽度,γ表示发射脉冲的调频率,λ表示波长,Tb为脉冲簇发射长度,R(t;r,x)表示斜距;
(b)将步骤(a)中的所述去载频后的回波在距离向和方位向分别进行脉冲压缩处理。该方法根据线性调频脉冲信号时间和频率的对应关系,进行方位向和距离向的两维去调频;即将所述(1)式乘以exp[-j2πf0τr-jπγ(t-τr)2],记为(2)式,
其中f0为中心频率,τr为距离向的延迟时间;
随后,在所述多普勒锐化处理单元中,对该(2)式进行距离向FFT,记为(3)式,
FFT[sr_dechirp(t,τ)] (3)
然后,将(3)式再乘以exp[-jπfrg 2/γ],记为(4)式,
其中,frg为距离向频率;
最后,再进行距离向IFFT,记为(5)式
上述的FFT/IFFT运算有很好的高效性,比在时域作卷积运算有效。
通过上述操作将步骤(a)中的所述去载频后的回波进行方位向和距离向的两维去调频,得到回波基带信号在距离快时间域和方位慢时间域的表述形式;
通过上述的距离向FFT/IFFT运算,消除残余视频相位误差,校正距离向带宽随着距离向延伸产生的带宽展宽;
最后,再对上述结果进行距离向FFT、方位向FFT处理,即
FFTτ{FFTt[srd(t,τ)]}
得到由多个分辨单元组成的二维的高空间分辨率的处理数据;
(c)旋转所述具有扫描体制的雷达散射计的观测几何可知,所述观测方位角与方位向分辨率之间的关系,用(6)式表示为:
其中,δaz为考虑地球斜视影响和观测***的方位向分辨率,f(θaz)为由于斜视变形而影响方位向分辨率的降低因子,θaz天线方位角,λ为载波波长,R0为观测斜距,Vsat为卫星速度,Tdw为有效驻留时间;
因此,根据该方位向分辨率来划分观测单元,实现高空间分辨率的观测;而所述观测单元是由所述步骤(b)中得到的分辨单元的处理数据组成;
(d)在完成所述步骤(c)后,所述步骤(b)中得到的分辨单元的处理数据下传到星下的地面接收站进行处理,进一步得到位于刈幅两侧的观测的后向散射系数,所述后向散射系数作为观测目标特性的反演数据。
所述高采样重建方法进一步包括以下具体处理过程:
(aa)由所述高分辨处理单元下传得到的地面反射的条带回波功率数据来计算条带后向散射系数,记为那么,所述条带后向散射系数用(7)式表示为:
其中,Ps为散射计条带回波功率,X为定标因子;
(bb)所述步骤(aa)中的条带后向散射系数实际为高分辨的后向散射系数的空间加权平均值,因此步骤(aa)中的条带后向散射系数进一步表示为(8)式:
其中,代表地面坐标系,用来表征条带范围内高分辨散射点的位置,为位于处的散射计空间响应函数,在散射计发射前通过实验手段获得,为位于处的高分辨后向散射系数,A表示条带照射区域,n表征散射计测得的条带后向散射系数受到热噪声、仪器测量误差、数据处理误差影响引入的测量误差,该误差服从均值为零方差为的高斯分布,其中,Kp为归一化测量误差,且满足 表示的方差,表示的均值;
(cc)在所述反卷积子单元中,在刈幅中间区域,即高采样重建区域,随着波束足印的运动轨迹,利用相同条带的多次重叠测量,增大采样频率,采用反卷积重建方法实现散射计方位向分辨率的提高;根据所述步骤(bb)给出的(8)式,采用反卷积方法对(8)式进行反卷积重建得到高分辨的σ0,所述过程等价为以下优化问题的求解:
其中,为数值保真项,保证估计的σ0符合(8)式,Ψ(σ0)表示对σ0施加的约束条件,估计的σ0满足已知的或假设的先验信息,α为可调参数;
α值越大,解越接近先验信息项;α值越小,解越接近于数值保真项;α为0,(8)式退化为最小二乘问题。
(dd)在完成所述步骤(cc)后,所得的估计的σ0即为具有高分辨率的位于刈幅中间的后向散射系数。
本发明的优点在于:
采用旋转扫描天线波束,实现宽刈幅观测(>1000km),能够实现全天时持续工作以及快速的全球覆盖,同时结合采用脉冲压缩和高分辨率信号处理能够实现高分辨率(优于2-5km);根据观测方位角的值来分区域采用多普勒波束锐化和高采样重建两种方法,实现散射计全刈幅范围的高分辨处理。
附图说明
图1是本发明的扫描波束高分辨率的雷达散射计的观测几何;
图2是本发明的扫描波束高分辨率的雷达散射计的***处理示意图;
图3是本发明的扫描波束高分辨率的雷达散射计的高分辨处理流程图;
1、地面足印 2、刈幅
3、反馈面 4、斜距
5、观测方位角 6、天线馈源
7、波束足印 8、多普勒锐化区域
9、高采样重建区域
具体实施方式
以下结合附图对本发明作进一步的详细说明。
本发明提供了一种扫描波束高分辨率的雷达散射计,同时,本发明提出一种结合波束圆锥扫描和高分辨率处理的技术方案,实现兼顾宽刈幅和高分辨率的观测要求。所述雷达散射计采用波束旋转扫描天线,如图2所示,其包括:天线单元、发射机单元、接收机单元、数字处理与控制单元、多普勒锐化处理单元和高分辨处理单元;所述天线单元采用窄波束圆锥扫描体制,用于实现旋转笔形波束扫描,所述发射机单元生成并发射线性调频脉冲信号,在发射该线性调频脉冲信号时,对应一个观测方位角所述接收机单元接收该脉冲信号并对其进行低噪声放大、中频接收、混频去斜坡和相位检波的处理,提取出回波信号,所述数字处理与控制单元记录所述发射信号对应的观测方位角当所述观测方位角为刈幅2两侧时,波束照射区域为多普勒波束锐化区域8,在所述多普勒锐化处理单元8中,采用多普勒波束锐化方法,将去载频后的回波进行方位向和距离向的两维去调频,再进行距离向和方位向脉冲压缩,得到由多个分辨单元组成的二维的高分辨率的处理数据,从而得到高分辨后向散射系数;当所述观测方位角在刈幅2中间时,所述波束照射区域为高采样重建区域9,在所述高分辨处理单元9中,采用高采样重建方法,测得条带回波功率,并计算条带后向散射系数,再根据反卷积法,得到高分辨后向散射系数。
所述天线单元采用窄波束圆锥扫描体制,用于实现旋转笔形波束扫描;其包括反射面天线和扫描伺服机构;该扫描伺服机构用于控制天线旋转。
如图1所示,地面足印1反射信号到反馈面3,反射路径称为斜距4,当所述笔形波束天线的观测方向与飞行方向垂直时,即所述观测方位角为90°,则所述波束照射区域为多普勒波束锐化区域8,即刈幅2两侧,在实际操作过程中,如图2所示,在星上进行高分辨处理,采用多普勒波束锐化方法来提高所述雷达散射计的分辨率;
(3)如图1所示,当所述笔形波束天线的观测方向与卫星方向平行时,即所述观测方位角为0°,则所述波束照射区域为高采样重建区域9,即刈幅2的中间,所述高采样重建区域9为非多普勒波束锐化区域,如图2所示,在星下进行高空间分布处理,采用高采样重建方法来提高所述雷达散射计的分辨率;
如图3所示,所述多普勒波束锐化方法进一步包括以下具体处理过程:
(a)将所述多普勒锐化处理单元对经过所述接收机单元和所述数字处理与控制单元的回波信号进行去载频,得到雷达在一个Burst数据块内获得的基带回波数据sr(t,τ),记为(1)式,
其中,
其中积分区域为扫描场景范围,c为光速,σ(x,r)表示地面(x,r)处的条带后向散射系数,为天线方向图,即 为天线方位角,G0为天线的最大增益,βaz为天线方位向波束宽度;t和τ分别对应表示方位向的“慢时间”和表示距离向的“快时间”,τp表示发射脉冲宽度,γ表示发射脉冲的调频率,λ表示波长,Tb为脉冲簇发射长度,R(t;r,x)表示斜距;
(b)将步骤(a)中的所述去载频后的回波在距离向和方位向分别进行脉冲压缩处理。该方法根据线性调频脉冲信号时间和频率的对应关系,进行方位向和距离向的两维去调频;即将所述(1)式乘以exp[-j2πf0τr-jπγ(t-τr)2],记为(2)式,
其中f0为中心频率,τr为距离向的延迟时间;
随后,在所述多普勒锐化处理单元中,对该(2)式进行距离向FFT,记为(3)式,
FFT[sr_dechirp(t,τ)] (3)
然后,将(3)式再乘以exp[-jπfrg 2/γ],记为(4)式,
其中,frg为距离向频率;
最后,再进行距离向IFFT,记为(5)式
上述的FFT/IFFT运算有很好的高效性,比在时域作卷积运算有效。
通过上述操作将步骤(a)中的所述去载频后的回波进行方位向和距离向的两维去调频,得到回波基带信号在距离快时间域和方位慢时间域的表述形式;
通过上述的距离向FFT/IFFT运算,消除残余视频相位误差,校正距离向带宽随着距离向延伸产生的带宽展宽;
最后,再对上述结果进行距离向FFT、方位向FFT处理,即
FFTτ{FFTt[srd(t,τ)]}
得到由多个分辨单元组成的二维的高空间分辨率的处理数据;
(c)旋转所述具有扫描体制的雷达散射计的观测几何可知,如图3所示,所述观测方位角与方位向分辨率之间的关系,用(6)式表示为:
其中,δaz为考虑地球斜视影响和观测***的方位向分辨率,f(θaz)为由于斜视变形而影响方位向分辨率的降低因子,θaz天线方位角,λ为载波波长,R0为观测斜距,Vsat为卫星速度,Tdw为有效驻留时间;
因此,根据该方位向分辨率来划分观测单元,实现高空间分辨率的观测;而所述观测单元是由所述步骤(b)中得到的分辨单元的处理数据组成;
(d)在完成所述步骤(c)后,所述步骤(b)中得到的分辨单元的处理数据下传到星下的地面接收站进行处理,进一步得到位于刈幅两侧的观测的后向散射系数,所述后向散射系数作为观测目标特性的反演数据。
所述多普勒波束锐化方法应用于空对地中提供高分辨率的地图测绘和高分辨率的局部放大测绘,以及空对空敌情判断状态可分辨出密集编队的群目标。
如图3所示,所述高采样重建方法进一步包括以下具体处理方法:
(aa)由所述图2中的高分辨处理单元下传得到的地面反射的条带回波功率数据来计算条带后向散射系数,记为那么,所述条带后向散射系数可以用(7)式表示为:
其中,Ps为散射计条带回波功率,X为定标因子;
(bb)所述步骤(aa)中的条带后向散射系数实际为高分辨的后向散射系数的空间加权平均值,因此步骤(aa)中的条带后向散射系数可进一步表示为(8)式:
其中,代表地面坐标系,用来表征条带范围内高分辨散射点的位置,为位于处的散射计空间响应函数,可在散射计发射前通过实验手段获得,为位于处的高分辨后向散射系数,A表示条带照射区域,n表征散射计测得的条带后向散射系数受到热噪声、仪器测量误差、数据处理误差等影响引入的测量误差,该误差服从均值为零方差为的高斯分布,其中,Kp为归一化测量误差,且满足 表示的方差,表示的均值;
(cc)如图1所示,在所述反卷积子单元中,在刈幅中间区域,即高采样重建区域9,随着波束足印7的运动轨迹,可利用相同条带的多次重叠测量,增大采样频率,采用反卷积重建方法实现散射计方位向分辨率的提高;根据所述步骤(bb)给出的(8)式,可采用反卷积方法对(8)式进行反卷积重建得到高分辨的σ0,所述过程可等价为以下优化问题的求解:
其中,为数值保真项,保证估计的σ0符合(8)式,Ψ(σ0)表示对σ0施加的约束条件,使估计的σ0满足已知的(或假设的)先验信息,α为可调参数;
α值越大,解越接近先验信息项;α值越小,解越接近于数值保真项;当α为0时,(8)式退化为最小二乘问题。
(dd)在完成所述步骤(cc)后,所得的估计的σ0即为具有高分辨率的位于刈幅中间的后向散射系数。
所述高分辨雷达散射计,除了可应用在诸如雪水当量、地表冻融、极地冰盖冰架等陆地观测中还可应用在海面风场、热带气旋等海洋观测。
最后需要说明的是,具体实施方式中所述的实验用图仅用来说明本发明的技术方案软件算法的可行性而非局限于此例,算法已经经过大量实验数据验证,是真实可靠的,搭配硬件便可实现本发明的技术方案。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (4)
1.一种扫描波束高分辨率的雷达散射计,其特征在于,所述雷达散射计采用笔形波束圆锥扫描体制并且实现全球宽刈幅覆盖和高分辨率观测,其包括:天线单元、发射机单元、接收机单元、信号处理与控制单元、多普勒锐化处理单元和高分辨处理单元;所述天线单元采用窄波束圆锥扫描体制,用于实现旋转笔形波束扫描,所述发射机单元生成并发射线性调频脉冲信号,在发射该线性调频脉冲信号时,对应一个观测方位角所述接收机单元接收该脉冲信号并对其进行低噪声放大、中频接收、混频去斜坡和相位检波的处理,提取出回波信号,所述数字处理与控制单元记录所述发射信号对应的观测方位角当所述观测方位角为刈幅两侧时,波束照射区域为多普勒波束锐化区域,在所述多普勒锐化处理单元,采用多普勒波束锐化方法,将去载频后的回波进行方位向和距离向的两维去调频,再进行距离向和方位向脉冲压缩,得到由多个分辨单元组成的二维的高分辨率的处理数据,从而得到高分辨后向散射系数;当所述观测方位角在刈幅中间时,所述波束照射区域为高采样重建区域,在所述高分辨处理单元中,采用高采样重建方法,测得条带回波功率,并计算条带后向散射系数,再根据反卷积法,得到高分辨后向散射系数。
2.根据权利要求1所述的一种扫描波束高分辨率的雷达散射计,其特征在于,所述天线单元包括一个反射面天线和扫描伺服机构;该扫描伺服机构用于控制天线旋转。
3.根据权利要求1所述的一种扫描波束高分辨率的雷达散射计,其特征在于,所述多普勒波束锐化方法进一步包括:
(a)将所述多普勒锐化处理单元对经过所述接收机单元和所述数字处理与控制单元的回波信号进行去载频,得到雷达在一个Burst数据块内获得的基带回波数据sr(t,τ),记为(1)式,
其中,
其中,积分区域为扫描场景范围,c为光速,σ(x,r)表示地面(x,r)处的条带后向散射系数,为天线方向图,即 为天线方位角,G0为天线的最大增益,βaz为天线方位向波束宽度;t和τ分别对应表示方位向的“慢时间”和表示距离向的“快时间”,τp表示发射脉冲宽度,γ表示发射脉冲的调频率,λ表示波长,Tb为脉冲簇发射长度,R(t;r,x)表示斜距;
(b)将步骤(a)中的所述去载频后的回波在距离向和方位向分别进行脉冲压缩处理;该方法根据线性调频脉冲信号时间和频率的对应关系,进行方位向和距离向的两维去调频;即将所述(1)式乘以exp[-j2πf0τr-jπγ(t-τr)2],记为(2)式,
其中f0为中心频率,τr为距离向的延迟时间;
随后,在所述多普勒锐化处理单元中,对该(2)式进行距离向FFT,记为(3)式,
FFT[sr_dechirp(t,τ)] (3)
然后,将(3)式再乘以exp[-jπfrg 2/γ],记为(4)式,
其中,frg为距离向频率;
最后,再进行距离向IFFT,记为(5)式
上述的FFT/IFFT运算有很好的高效性,比在时域作卷积运算有效;
通过上述操作将步骤(a)中的所述去载频后的回波进行方位向和距离向的两维去调频,得到回波基带信号在距离快时间域和方位慢时间域的表述形式;
通过上述的距离向FFT/IFFT运算,消除残余视频相位误差,校正距离向带宽随着距离向延伸产生的带宽展宽;
最后,再对上述结果进行距离向FFT、方位向FFT处理,即
FFTτ{FFTt[srd(t,τ)]}
得到由多个分辨单元组成的二维的高空间分辨率的处理数据;
(c)旋转所述具有扫描体制的雷达散射计的观测几何可知,所述观测方位角与方位向分辨率之间的关系,用(6)式表示为:
其中,δaz为考虑地球斜视影响和观测***的方位向分辨率,f(θaz)为由于斜视变形而影响方位向分辨率的降低因子,θaz天线方位角,λ为载波波长,R0为观测斜距,Vsat为卫星速度,Tdw为有效驻留时间;
因此,根据该方位向分辨率来划分观测单元,实现高空间分辨率的观测;而所述观测单元是由所述步骤(b)中得到的分辨单元的处理数据组成;
(d)在完成所述步骤(c)后,所述步骤(b)中得到的分辨单元的处理数据下传到星下的地面接收站进行处理,进一步得到位于刈幅两侧的观测的后向散射系数,所述后向散射系数作为观测目标特性的反演数据。
4.根据权利要求1所述的一种扫描波束高分辨率的雷达散射计,其特征在于,所述高采样重建方法进一步包括:
(aa)由所述高分辨处理单元下传得到的地面反射的条带回波功率数据来计算条带后向散射系数,记为那么,所述条带后向散射系数用(7)式表示为:
其中,Ps为散射计条带回波功率,X为定标因子;
(bb)所述步骤(aa)中的条带后向散射系数实际为高分辨的后向散射系数的空间加权平均值,因此步骤(aa)中的条带后向散射系数进一步表示为(8)式:
其中,代表地面坐标系,用来表征条带范围内高分辨散射点的位置,为位于处的散射计空间响应函数,在散射计发射前通过实验手段获得,为位于处的高分辨后向散射系数,A表示条带照射区域,n表征散射计测得的条带后向散射系数受到热噪声、仪器测量误差、数据处理误差影响引入的测量误差,该误差服从均值为零方差为的高斯分布,其中,Kp为归一化测量误差,且满足 表示的方差,表示的均值;
(cc)在所述反卷积子单元中,在刈幅中间区域,即高采样重建区域,沿着波束足印的运动轨迹,利用相同条带的多次重叠测量,增大采样频率,采用反卷积重建方法实现散射计方位向分辨率的提高;根据所述步骤(bb)给出的(8)式,采用反卷积方法对(8)式进行反卷积重建得到高分辨的σ0,所述过程等价为以下优化问题的求解:
其中,为数值保真项,保证估计的σ0符合(8)式,Ψ(σ0)表示对σ0施加的约束条件,估计的σ0满足已知的或假设的先验信息,α为可调参数;
α值越大,解越接近先验信息项;α值越小,解越接近于数值保真项;α为0,(8)式退化为最小二乘问题;
(dd)在完成所述步骤(cc)后,所得的估计的σ0即为具有高分辨率的位于刈幅中间的后向散射系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610656894.9A CN106093932B (zh) | 2016-02-29 | 2016-08-11 | 一种扫描波束高分辨率的雷达散射计 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2016101133090 | 2016-02-29 | ||
CN201610113309 | 2016-02-29 | ||
CN201610656894.9A CN106093932B (zh) | 2016-02-29 | 2016-08-11 | 一种扫描波束高分辨率的雷达散射计 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106093932A true CN106093932A (zh) | 2016-11-09 |
CN106093932B CN106093932B (zh) | 2019-01-25 |
Family
ID=57457361
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610656894.9A Active CN106093932B (zh) | 2016-02-29 | 2016-08-11 | 一种扫描波束高分辨率的雷达散射计 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106093932B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106501777A (zh) * | 2016-12-09 | 2017-03-15 | 北京环境特性研究所 | 基于射线追踪的散射中心来源诊断方法 |
CN109709548A (zh) * | 2018-12-29 | 2019-05-03 | 内蒙古工业大学 | 一种全极化旋转微变监测雷达成像方法和雷达*** |
CN110221299A (zh) * | 2019-07-04 | 2019-09-10 | 南京信息工程大学 | 一种星载双通道双波束InSAR测流*** |
CN110609287A (zh) * | 2018-06-14 | 2019-12-24 | 中国科学院国家空间科学中心 | 一种双频雷达散射计及同时测量海面风场和流场的方法 |
CN111096761A (zh) * | 2018-10-29 | 2020-05-05 | 上海西门子医疗器械有限公司 | 修正楔形滤波器散射的方法、装置和相关设备 |
CN111175720A (zh) * | 2020-01-15 | 2020-05-19 | 中国科学院国家空间科学中心 | 一种海面风场星上快速反演方法及*** |
CN114518563A (zh) * | 2022-03-01 | 2022-05-20 | 国家卫星海洋应用中心 | 卫星微波散射计雷达后向散射系数获取方法及相应的装置 |
CN116502476A (zh) * | 2023-06-28 | 2023-07-28 | 中国科学院空天信息创新研究院 | 一种基于非线性频率扫描响应的sar***设计方法 |
CN117518167A (zh) * | 2024-01-04 | 2024-02-06 | 中国科学院空天信息创新研究院 | 一种基于多通道体制的宽幅sar扫描模式***设计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1301970A (zh) * | 1999-12-30 | 2001-07-04 | 中国科学院空间科学与应用研究中心 | 实时变参数滤波笔形波束扫描雷达散射计及其设计方法 |
CN101029928A (zh) * | 2006-02-27 | 2007-09-05 | 中国科学院空间科学与应用研究中心 | 一种收发双波束天线星载扫描雷达散射计 |
CN101672914A (zh) * | 2009-10-26 | 2010-03-17 | 西安空间无线电技术研究所 | 一种圆锥扫描高分辨率微波散射计的信号处理方法 |
CN103675788A (zh) * | 2013-12-05 | 2014-03-26 | 中国科学院空间科学与应用研究中心 | 散射计回波信号中降水回波及海面后向散射的分离方法 |
-
2016
- 2016-08-11 CN CN201610656894.9A patent/CN106093932B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1301970A (zh) * | 1999-12-30 | 2001-07-04 | 中国科学院空间科学与应用研究中心 | 实时变参数滤波笔形波束扫描雷达散射计及其设计方法 |
CN101029928A (zh) * | 2006-02-27 | 2007-09-05 | 中国科学院空间科学与应用研究中心 | 一种收发双波束天线星载扫描雷达散射计 |
CN101672914A (zh) * | 2009-10-26 | 2010-03-17 | 西安空间无线电技术研究所 | 一种圆锥扫描高分辨率微波散射计的信号处理方法 |
CN103675788A (zh) * | 2013-12-05 | 2014-03-26 | 中国科学院空间科学与应用研究中心 | 散射计回波信号中降水回波及海面后向散射的分离方法 |
Non-Patent Citations (3)
Title |
---|
GANG WANG等: "Azimuth High Resolution for a Conically Scanned Pencil-Beam Scatterometer Using Rotating Azimuth Doppler Discrimination", 《IEEE XPLORE》 * |
MICHAEL W.SPENCER 等: "Improved Resolution Backscatter Measurements with the SeaWinds Pencil-Beam Scatterometer", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
刘凡 等: "一种基于最小二乘直线拟合的高分辨率DBS成像算法", 《电子信息学报》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106501777B (zh) * | 2016-12-09 | 2019-03-26 | 北京环境特性研究所 | 基于射线追踪的散射中心来源诊断方法 |
CN106501777A (zh) * | 2016-12-09 | 2017-03-15 | 北京环境特性研究所 | 基于射线追踪的散射中心来源诊断方法 |
CN110609287B (zh) * | 2018-06-14 | 2021-06-08 | 中国科学院国家空间科学中心 | 一种双频雷达散射计及同时测量海面风场和流场的方法 |
CN110609287A (zh) * | 2018-06-14 | 2019-12-24 | 中国科学院国家空间科学中心 | 一种双频雷达散射计及同时测量海面风场和流场的方法 |
CN111096761B (zh) * | 2018-10-29 | 2024-03-08 | 上海西门子医疗器械有限公司 | 修正楔形滤波器散射的方法、装置和相关设备 |
CN111096761A (zh) * | 2018-10-29 | 2020-05-05 | 上海西门子医疗器械有限公司 | 修正楔形滤波器散射的方法、装置和相关设备 |
CN109709548B (zh) * | 2018-12-29 | 2023-02-17 | 内蒙古工业大学 | 一种全极化旋转微变监测雷达成像方法和雷达*** |
CN109709548A (zh) * | 2018-12-29 | 2019-05-03 | 内蒙古工业大学 | 一种全极化旋转微变监测雷达成像方法和雷达*** |
CN110221299B (zh) * | 2019-07-04 | 2023-02-10 | 南京信息工程大学 | 一种星载双通道双波束InSAR测流*** |
CN110221299A (zh) * | 2019-07-04 | 2019-09-10 | 南京信息工程大学 | 一种星载双通道双波束InSAR测流*** |
CN111175720A (zh) * | 2020-01-15 | 2020-05-19 | 中国科学院国家空间科学中心 | 一种海面风场星上快速反演方法及*** |
CN111175720B (zh) * | 2020-01-15 | 2022-03-08 | 中国科学院国家空间科学中心 | 一种海面风场星上快速反演方法及*** |
CN114518563A (zh) * | 2022-03-01 | 2022-05-20 | 国家卫星海洋应用中心 | 卫星微波散射计雷达后向散射系数获取方法及相应的装置 |
CN114518563B (zh) * | 2022-03-01 | 2022-09-09 | 国家卫星海洋应用中心 | 卫星微波散射计雷达后向散射系数获取方法及相应的装置 |
CN116502476A (zh) * | 2023-06-28 | 2023-07-28 | 中国科学院空天信息创新研究院 | 一种基于非线性频率扫描响应的sar***设计方法 |
CN116502476B (zh) * | 2023-06-28 | 2023-09-08 | 中国科学院空天信息创新研究院 | 一种基于非线性频率扫描响应的sar***设计方法 |
CN117518167A (zh) * | 2024-01-04 | 2024-02-06 | 中国科学院空天信息创新研究院 | 一种基于多通道体制的宽幅sar扫描模式***设计方法 |
CN117518167B (zh) * | 2024-01-04 | 2024-03-12 | 中国科学院空天信息创新研究院 | 一种基于多通道体制的宽幅sar扫描模式***设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106093932B (zh) | 2019-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106093932A (zh) | 一种扫描波束高分辨率的雷达散射计 | |
EP2491418B1 (en) | Resolution enhancement system for networked radars | |
CN110109117B (zh) | 基于频控阵的星载合成孔径雷达卷积调制干扰方法 | |
CN104950306B (zh) | 一种海杂波背景下前视海面目标角超分辨成像方法 | |
CN103529437B (zh) | 系留气球载相控阵雷达在多目标下分辨空地目标的方法 | |
CN107976660B (zh) | 弹载多通道雷达超低空目标分析与多径回波建模方法 | |
CN103487803B (zh) | 迭代压缩模式下机载扫描雷达成像方法 | |
CN102313887B (zh) | 一种星机联合双基地合成孔径雷达成像方法 | |
CN107966688A (zh) | 基于相位干涉技术的宽带雷达目标速度解模糊方法 | |
CN102707269B (zh) | 一种机载雷达距离走动校正方法 | |
CN104977582A (zh) | 一种实现扫描雷达方位超分辨成像的解卷积方法 | |
CN109655819B (zh) | 一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法 | |
CN103033812B (zh) | 机载双基地sar的波束追赶空间同步方法 | |
CN105738887A (zh) | 基于多普勒通道划分的机载雷达杂波功率谱的优化方法 | |
CN104280566A (zh) | 基于空时幅相估计的低空风切变风速估计方法 | |
Nohara et al. | SAR-GMTI processing with Canada's Radarsat 2 satellite | |
CN112098999A (zh) | 一种高动态雷达导引头掠海目标电磁信号建模方法 | |
CN107607945A (zh) | 一种基于空间嵌入映射的扫描雷达前视成像方法 | |
CN110879391B (zh) | 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法 | |
US20180074180A1 (en) | Ultrafast target detection based on microwave metamaterials | |
CN112147593B (zh) | 一种高速密集***物破片目标的四维参数估计方法 | |
CN106526544B (zh) | 基于高超声速平台的mimosar杂波抑制方法 | |
CN103616670B (zh) | 基于浮空平台改进的sar散射波干扰方法 | |
CN110308445B (zh) | 一种基于车载数字阵调频连续波雷达的成像方法 | |
CN112213699A (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 |