CN114488134A - 星载多通道gnss-s雷达视频成像系统及舰船轨迹提取方法 - Google Patents
星载多通道gnss-s雷达视频成像系统及舰船轨迹提取方法 Download PDFInfo
- Publication number
- CN114488134A CN114488134A CN202210314900.8A CN202210314900A CN114488134A CN 114488134 A CN114488134 A CN 114488134A CN 202210314900 A CN202210314900 A CN 202210314900A CN 114488134 A CN114488134 A CN 114488134A
- Authority
- CN
- China
- Prior art keywords
- ship target
- radar
- unit
- signal
- echo signals
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 52
- 238000000605 extraction Methods 0.000 title claims abstract description 26
- 238000012545 processing Methods 0.000 claims abstract description 60
- 238000001514 detection method Methods 0.000 claims abstract description 50
- 230000004927 fusion Effects 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 14
- 238000007499 fusion processing Methods 0.000 claims description 14
- 238000000926 separation method Methods 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 12
- 238000011084 recovery Methods 0.000 claims description 11
- 238000013139 quantization Methods 0.000 claims description 8
- 230000003321 amplification Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000013500 data storage Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000005314 correlation function Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 3
- 230000010365 information processing Effects 0.000 claims description 3
- 230000010287 polarization Effects 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
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/88—Radar or analogous systems specially adapted for specific applications
-
- 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/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
-
- 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/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
-
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/14—Receivers specially adapted for specific applications
-
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
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)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种星载多通道GNSS‑S雷达视频成像系统及舰船轨迹提取方法,系统包括:星载多通道GNSS‑S雷达(10),用于对海面舰船目标的回波信号进行同步接收;方位向智能DBF处理单元(20),用于对回波信号进行处理,形成二维回波信号;多源雷达融合成像单元(30),用于对导航卫星散射信号进行处理,获得高信杂噪比雷达图像;舰船目标检测单元(40),用于对高信杂噪比雷达图像进行舰船目标检测与定位,并提取舰船目标的切片信息;舰船目标视频形成单元(50),用于形成舰船目标的视频图像;舰船目标轨迹提取单元(60),用于形成舰船目标的运动轨迹。本发明具有长时间视频成像、批量舰船目标轨迹提取、自身隐蔽性强、系统功耗低等优势。
Description
技术领域
本发明涉及一种星载多通道GNSS-S雷达视频成像系统及舰船轨迹提取方法。
背景技术
海面舰船目标探测与视频成像一直是科学研究的热点,但是受海上云雨雾等天气影响,光学传感器难以发挥其高分辨率成像与识别优势。而星载SAR系统则能够穿云透雾,可实现全天时全天候对海观测,适合海面舰船目标的高分辨率成像与探测应用,但现有的星载SAR系统难以实现舰船目标的长时间视频成像与运动轨迹提取。而由于利用GNSS-S进行海面舰船目标探测时不需要主动发射信号,仅需接收导航卫星信号的散射信号即可实现目标探测,因此具有低功耗优势,相比于大功耗SAR,可实现长时间探测与跟踪,更适合舰船目标长时间跟踪与视频成像。此外,导航卫星信号还具有全球覆盖优势,且海面任一空间均可同时接收到多颗导航卫星信号,可实现多源信息联合探测。
对此,一些技术提出利用导航卫星信号的反射信号(GNSS-R)实现海面风场探测,并完成了低轨卫星搭载试验。同时,另一些技术提出在地基利用导航卫星信号的散射信号(GNSS-S)实现星发地收双站SAR成像,但受限于导航卫星信号的有效带宽,其对应的成像分辨率一般仅在十米量级。因此,如何利用星载GNSS-S雷达实现海面舰船目标的视频成像与轨迹提取成为领域内亟待解决的问题。
发明内容
本发明的目的在于提供一种星载多通道GNSS-S雷达视频成像系统及舰船轨迹提取方法。
为实现上述发明目的,本发明提供一种星载多通道GNSS-S雷达视频成像系统及舰船轨迹提取方法,系统包括:
星载多通道GNSS-S雷达,用于对海面舰船目标的回波信号进行同步接收;
方位向智能DBF(数字波束形成)处理单元,用于对回波信号进行处理,形成二维回波信号;
多源雷达融合成像单元,用于对导航卫星散射信号进行处理,获得高信杂噪比雷达图像;
舰船目标检测单元,用于对高信杂噪比雷达图像进行舰船目标检测与定位,并提取舰船目标的切片信息;
舰船目标视频形成单元,用于形成舰船目标的视频图像;
舰船目标轨迹提取单元,用于形成舰船目标的运动轨迹。
根据本发明的一个方面,所述星载多通道GNSS-S雷达可获得N个回波信号,记为sn(t),n=1,2,...,N;其中,N为雷达天线沿方位向的接收通道个数;
所述方位向智能DBF处理单元对N个回波信号进行智能DBF处理、多源信号分离与时序恢复,形成M个数字域高增益窄波束,获得M个高信噪比回波信号,记第k帧的M个高信噪比回波信号为rm(χk,t),m=1,2,...,M;并通过多源信号分离处理使每个回波信号分离出P个回波信号,记为:
rm,p(χk,t),m=1,2,...,M,p=1,2,...,P;
再经过时序恢复形成M×P个二维回波信号,记为r’m,p(χk,τk,t);其中,χk为视频帧序号,τk为慢时序列;
所述多源雷达融合成像单元对M个数字波束的P颗导航卫星散射信号进行双站雷达成像与多源雷达图像融合处理,获得第k帧的M幅高信杂噪比雷达图像Im(χk,dAm,dRm),m=1,2,...,M;其中,dAm、dRm分别为第m幅图像的方位向网格单元序列与距离向网格单元序列;
所述舰船目标检测单元对第k帧的M幅高信杂噪比雷达图像进行舰船目标检测与定位,并提取舰船目标的切片信息,记第j个舰船目标的第k帧位置与图像切片分别为Ωm,j(χk,xj,yj)与Jm,j(χk,dAj,dRj);其中,xj为第j个舰船目标在地距平面的方位向位置,yj为第j个舰船目标在地距平面的距离向位置,dAj为第j个舰船目标切片沿方位向的网格单元序列,dRj为第j个舰船目标切片沿距离向的网格单元序列;
根据本发明的一个方面,所述星载多通道GNSS-S雷达包括:
L波段方位向多通道天线,用于对海面舰船目标GNSS-S信号进行接收,输出N路GNSS-S信号;
N个高灵敏度数字接收机,用于对GNSS-S信号进行放大接收、下变频与采样量化,获得N路数字域GNSS-S信号sn(t);
GNSS-S雷达主控单元,用于对GNSS-S信号的放大接收、采样量化、数据存储与传输进行控制;
频率信号发生器,用于产生下变频所需的本振信号fLO、采样量化所需的采样时钟fADC、雷达主控所需的工作时钟fFPGA、处理单元所需的处理时钟fSP;下角标LO的含义为本振频率信号;下角标SP的含义为信号处理的工作时钟信号;(ADC为模数转换器;FPGA为现场可编程门阵列)
分布式数据存储器,用于对多路数字域GNSS-S信号与中间计算变量进行分布式并行存储;
阵列信号与多源信息处理单元,用于对多路GNSS-S信号进行智能DBF、多源雷达融合成像、舰船目标检测与定位。
根据本发明的一个方面,所述方位向智能DBF处理单元包括:
波束指向求解单元,用于根据舰船目标的位置信息计算DBF波束指向,获得第k帧的M个数字窄波束指向,记第k帧的第m个波束指向为每个数字窄波束对应的合成孔径时间为Tm,m=1,2,3,…,M,Tm为1-5秒;其中,为第m个波束的方位角,θm为第m个波束的入射角;
最优权值求解单元,用于对智能DBF的M个数字窄波束对应的M组权值进行计算,获得M组智能DBF复数权值,记第k帧的M组DBF复数权值为wk,m,n;
多源回波信号分离单元,用于对M个回波信号rm(χk,t)进行多源信号分离,获得M×P个更高信噪比回波信号;
多源回波信号时序恢复单元,用于对M×P个回波信号的时序进行恢复,形成M×P个二维回波信号。
根据本发明的一个方面,所述多源回波信号分离单元根据卫星导航接收机捕获跟踪获得的P颗导航卫星信号信息,分别对M个高信噪比回波信号rm(χk,t)进行互相关处理,每个波束区域均可获得P个回波信号,共可获得M×P个更高信噪比回波信号,记第k帧、第m个波束、第p颗导航卫星的回波信号为rm,p(χk,t),m=1,2,...,M,p=1,2,...,P,则有rm,p(χk,t)=R(rm(χk,t),Zp(t)),其中,R(·)为互相关函数,Zp(t)为捕获跟踪获得的第p颗导航卫星信号;
所述多源回波信号时序恢复单元根据卫星导航接收机输出的P颗导航卫星位置与速度信息、低轨GNSS-S雷达天线位置与速度信息,将一维回波信号rm,p(χk,t)转换成二维回波信号r’m,p(χk,τk,t),其中,τk与t分别为方位向慢时变量与距离向快时变量。
根据本发明的一个方面,所述多源雷达融合成像单元包括:
双站后向投影(BP)雷达成像处理单元,用于分别对第k帧的M×P个二维GNSS-S回波信号r’m,p(χk,τk,t)进行双站BP雷达成像,获得第k帧的M×P幅双站雷达图像Im,p(χk,dAm,dRm),其中,dAm与dRm分别为第m波束区域的方位向网格变量与距离向网格变量,χk为视频图像中的图像帧数,共分为M个波束区域,每个波束区域均包含P幅双站雷达图像;
(“上横线”表示归一化处理后的图像,即与处理前的图像进行区分)
多源雷达图像融合处理单元,用于分别对第k帧的M个波束区域的P幅双站雷达图像进行融合处理,获得第k帧的M幅高信杂噪比雷达图像Im(χk,dAm,dRm),m=1,2,...,M,则有:
其中,w'p为融合处理的加权值。
根据本发明的一个方面,所述舰船目标检测单元包括:
杂波功率分布统计单元,用于对第k帧、第m个波束的雷达图像Im(χk,dAm,dRm)的海杂波功率分布进行统计分析,功率分布概率函数为fPD(Γsea),其中,Γsea为海杂波功率变量;下角标PD的含义为功率分布;
检测阈值求取单元,用于对舰船目标检测的阈值进行计算,设定虚警率为fR,取10-2-10-4,则计算检测阈值为Γm,k;下角标R的含义为检测虚警率;
舰船目标功率比较单元,用于对第k帧、第m个波束的雷达图像Im(χk,dAm,dRm)进行舰船目标检测;
舰船目标位置计算单元,用于对舰船目标的位置进行计算,将网格点位置转换成地距平面位置信息,记第k帧、第m个波束的雷达图像中的第j个舰船目标的地距平面位置信息为Ωm,j(χk,xj,yj);
舰船目标图像切片提取单元,用于对舰船目标的雷达图像切片进行提取,记第k帧、第m个波束的雷达图像中的第j个舰船目标的雷达图像切片为Jm,j(χk,dAj,dRj)。
根据本发明的一个方面,所述舰船目标功率比较单元将雷达图像各网格点功率与检测阈值进行功率比较,当Im(χk,dAm,dRm)>Γm,k时,则第k帧、第m个波束的雷达图像的方位向第dAm个网格点、距离向第dRm个网格点处存在舰船目标。
根据本发明的一个方面,星载多通道GNSS-S雷达的轨道高度为H,采用一幅L波段方位向多通道天线实现海面舰船目标回波信号接收,每个天线通道在方位向具有宽波束角特性;沿方位向具有N个天线通道,N取12-36,每个天线通道工作频段为L波段,中心频率为1.575GHz,带宽大于50MHz;每个天线通道在方位向的波束宽度为θA,θA取10°-15°;每个天线通道在距离向的波束宽度为θR,θR取8°-12°;天线采用水平极化对舰船目标回波信号进行接收;天线以正侧视模式对海面舰船目标回波信号进行接收,通过方位向多通道天线的数字波束形成处理,在方位向形成M个数字域高增益窄波束,实现对舰船目标的方位向探测与跟踪;下角标A的含义为方位向;下角标R的含义为距离向;
将星载多通道GNSS-S雷达的星下点轨迹划分成K个时间段,对应K帧雷达图像,帧序号记为χk;在第k个时间段内包含M个波束区域对应的M幅雷达图像,记第k帧、第m个波束区域的雷达图像为Im(χk,dAm,dRm),m=1,2,...,M;
在进行舰船目标视频图像形成和舰船轨迹提取时,根据舰船目标位置Ωm,j(χk,xj,yj),计算获得第k+1帧的天线DBF波束指向为并获得第k+1帧M个高信噪比回波信号为rm(χk+1,t),其中,Xk+1表示下一帧图像的时刻或序号,即由于是视频,因此存在多帧图像;再获得第j个舰船目标的第k+1帧位置信息Ωm,j(χk+1,xj,yj)与雷达图像切片信息Jm,j(χk+1,dAj,dRj);将第j个舰船目标的K帧位置信息Ωm,j(χk,xj,yj),k=1,2,...,K进行关联处理,获得第j个舰船目标的运动轨迹
舰船轨迹提取方法,包括以下步骤:
a、对海面舰船目标的回波信号进行同步接收,并对回波信号进行处理,形成二维回波信号;
b、对导航卫星散射信号进行处理,获得高信杂噪比雷达图像;
c、对高信杂噪比雷达图像进行舰船目标检测与定位,并提取舰船目标的切片信息;
d、形成舰船目标的视频图像及运动轨迹。
根据本发明的构思,提出一种星载多通道GNSS-S雷达视频成像系统及舰船轨迹提取方法。其中,星载GNSS-S雷达采用L波段方位向多通道天线,每个天线通道在方位向具有宽波束角特性,以提升对海面舰船目标的观测时间。并通过对方位向多通道天线的多个回波信号进行智能数字波束形成(DBF)、多源雷达融合成像、舰船目标检测与定位、舰船目标视频图像形成等处理,获得海面舰船目标数量、位置与轨迹,以及舰船目标的视频图像。其中,利用智能DBF处理可以自主形成多个数字窄波束,从而对海面舰船目标进行多帧成像与检测,进而获得海面舰船目标的视频图像与运动轨迹。如此,相比于现有的主动星载雷达SAR成像与探测系统,本发明的星载多通道GNSS-S雷达视频成像系统及舰船轨迹提取方法具有长时间视频成像、批量舰船目标轨迹提取、自身隐蔽性强、系统功耗低等优势,具有较高的海上探测应用价值和广阔的市场应用前景。
附图说明
图1示意性表示本发明的一种实施方式的星载多通道GNSS-S雷达视频成像系统的工作流程图;
图2示意性表示本发明的一种实施方式的星载多通道GNSS-S雷达视频成像系统的工作场景图;
图3示意性表示本发明的一种实施方式的星载多通道GNSS-S雷达系统的组成图;
图4示意性表示本发明的一种实施方式的智能DBF处理流程图;
图5示意性表示本发明的一种实施方式的连续多帧雷达图像的空间位置关系图;
图6示意性表示本发明的一种实施方式的多源雷达融合成像处理流程图;
图7示意性表示本发明的一种实施方式的单个波束区域的海面目标检测与定位流程图。
具体实施方式
为了更清楚地说明本发明实施方式或现有技术中的技术方案,下面将对实施方式中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是本发明的一些实施方式,对于本领域普通技术人员而言,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
下面结合附图和具体实施方式对本发明作详细地描述,实施方式不能在此一一赘述,但本发明的实施方式并不因此限定于以下实施方式。
参见图1,本发明的星载多通道GNSS-S雷达视频成像系统,适用于高分宽幅SAR成像,可应用于天基分布式高分宽幅SAR系统关键技术研究,系统包括:星载多通道GNSS-S雷达10,用于对海面舰船目标的回波信号进行同步接收,可获得N个回波信号,记为sn(t),n=1,2,...,N,其中,N为雷达天线沿方位向的接收通道个数;方位向智能DBF处理单元20,用于对N个回波信号进行智能DBF处理、多源信号分离与时序恢复,形成M个数字域高增益窄波束,获得M个高信噪比回波信号,记第k帧的M个高信噪比回波信号为:rm(χk,t),m=1,2,...,M;并通过多源信号分离处理使每个回波信号分离出P个回波信号,记为:rm,p(χk,t),m=1,2,...,M,p=1,2,...,P;再经过时序恢复形成M×P个二维回波信号,记为r’m,p(χk,τk,t),其中,χk为视频帧序号,τk为慢时序列;多源雷达融合成像单元30,用于对M个数字波束的P颗导航卫星散射信号进行双站雷达成像与多源雷达图像融合处理,获得第k帧的M幅高信杂噪比雷达图像:Im(χk,dAm,dRm),m=1,2,...,M,其中,dAm、dRm分别为第m幅图像的方位向网格单元序列与距离向网格单元序列;舰船目标检测单元40,用于对第k帧的M幅高增益信杂噪比雷达图像进行舰船目标检测与定位,并提取舰船目标的切片信息,记第j个舰船目标的第k帧位置与图像切片分别为Ωm,j(χk,xj,yj)与Jm,j(χk,dAj,dRj),其中,xj为第j个舰船目标在地距平面的方位向位置,yj为第j个舰船目标在地距平面的距离向位置,dAj为第j个舰船目标切片沿方位向的网格单元序列,dRj为第j个舰船目标切片沿距离向的网格单元序列;舰船目标视频形成单元50,用于利用连续K帧图像切片形成舰船目标视频图像,记第j个舰船目标的视频图像为:舰船目标轨迹提取单元60,用于利用连续K帧位置信息形成舰船目标运动轨迹,记第j个舰船目标的运动轨迹为:
参见图2,星载多通道GNSS-S雷达的轨道高度为H,采用一幅L波段方位向多通道天线实现海面舰船目标回波信号接收,每个天线通道在方位向具有宽波束角特性,以提升对海面舰船目标的观测时间;沿方位向具有N个天线通道,N取12-36,每个天线通道工作频段为L波段,中心频率为1.575GHz,带宽大于50MHz;每个天线通道在方位向的波束宽度为θA,θA取10°-15°,以提高对海面目标的观测时间;每个天线通道在距离向的波束宽度为θR,θR取8°-12°,以增大有效成像幅宽;天线采用水平极化对舰船目标回波信号进行接收,以降低海杂波的干扰;天线以正侧视模式对海面舰船目标回波信号进行接收,通过方位向多通道天线的数字波束形成(DBF)处理,在方位向形成M个数字域高增益窄波束,实现对舰船目标的方位向探测与跟踪。
参见图3,星载多通道GNSS-S雷达10包括:L波段方位向多通道天线101,用于对海面舰船目标GNSS-S信号进行接收,输出N路GNSS-S信号;N个高灵敏度数字接收机102,用于对GNSS-S信号进行放大接收、下变频与采样量化,获得N路数字域GNSS-S信号sn(t);GNSS-S雷达主控单元103,用于对GNSS-S信号的放大接收、采样量化、数据存储与传输等进行控制;频率信号发生器104,用于产生下变频所需的本振信号fLO、采样量化所需的采样时钟fADC、雷达主控所需的工作时钟fFPGA、处理单元所需的处理时钟fSP;分布式数据存储器105,用于对多路数字域GNSS-S信号与中间计算变量等进行分布式并行存储;阵列信号与多源信息处理单元106,用于对多路GNSS-S信号进行智能DBF、多源雷达融合成像、舰船目标检测与定位等处理。
参见图4,方位向智能DBF处理单元20用于对海面舰船目标的多个回波信号进行智能DBF处理,包括:波束指向求解单元201,用于根据舰船目标的位置信息优化波束覆盖区域,计算DBF的波束个数与波束指向,获得第k帧的M个数字窄波束指向,记第k帧的第m个波束指向为每个数字窄波束对应的合成孔径时间为Tm,m=1,2,3,…,M,Tm为1-5秒,其中,为第m个波束的方位角,θm为第m个波束的入射角;最优权值求解单元202,用于根据DBF波束个数与波束指向,计算DBF的最优复数权值,即对智能DBF的M个数字窄波束对应的M组权值进行计算,获得M组智能DBF复数权值,记第k帧的M组DBF复数权值为wk,m,n;复数加权求和单元203,用于对多个天线通道输出的N路GNSS-S回波信号sn(t)进行复数加权求和处理,实现智能DBF处理,获得M个高增益窄波束,记第k帧的第m个波束指向为对应的第k帧的M个高信噪比回波信号为rm(χk,t),m=1,2,...,M,则有多源回波信号分离单元204,用于对M个高信噪比回波信号rm(χk,t)进行多源信号分离,获得M×P个更高信噪比回波信号,即根据卫星导航接收机捕获跟踪获得的P颗导航卫星信号信息,分别对M个高信噪比回波信号rm(χk,t)进行互相关处理
(即rm,p(χk,t)=R(rm(χk,t),Zp(t)));
每个波束区域均可以获得P个回波信号,则共可获得M×P个更高信噪比回波信号,记第k帧、第m个波束、第p颗导航卫星的回波信号为rm,p(χk,t),m=1,2,...,M,p=1,2,...,P,则有:
rm,p(χk,t)=R(rm(χk,t),Zp(t));
其中,R(·)为互相关函数,Zp(t)为捕获跟踪获得的第p颗导航卫星信号;多源回波信号时序恢复单元205,用于对M×P个回波信号的时序进行恢复,形成M×P个二维回波信号,即根据卫星导航接收机输出的P颗导航卫星位置与速度信息、低轨GNSS-S雷达天线位置与速度信息等,对一维回波信号rm,p(χk,t)的时序进行恢复,将一维回波信号rm,p(χk,t)转换成二维回波信号r’m,p(χk,τk,t),其中,τk与t分别为方位向慢时变量与距离向快时变量。
参见图5示出的连续多帧雷达图像的空间位置特点,本发明将星载GNSS-S雷达的星下点轨迹划分成K个时间段,对应K帧雷达图像,帧序号记为χk;在第k个时间段内包含M个波束区域对应的M幅雷达图像,记第k帧第m个波束区域的雷达图像为Im(χk,dAm,dRm),m=1,2,...,M。
参见图6,多源雷达融合成像单元30包括:双站BP雷达成像处理单元301,用于分别对第k帧的M×P个二维GNSS-S回波信号r’m,p(χk,τk,t)(即第m个数字窄波束区域的P个二维信号)进行双站BP雷达成像,获得第k帧的M×P幅双站雷达图像Im,p(χk,dAm,dRm)(即第m个波束区域的第p幅双站雷达图像),其中,dAm与dRm分别为第m波束区域的方位向网格变量与距离向网格变量,χk为视频图像中的图像帧数(即表示视频图像的第k帧图像),共分为M个波束区域,每个波束区域均包含P幅双站雷达图像;功率归一化单元302,用于分别对M×P幅双站雷达图像进行功率归一化处理,获得功率范围均小于等于1的M×P幅双站雷达图像则有多源雷达图像融合处理单元303,用于分别对第k帧的M个波束区域的P幅双站雷达图像进行融合处理(即对每个数字窄波束区域的P幅双站雷达图像进行融合处理),获得第k帧的M幅高信杂噪比雷达图像Im(χk,dAm,dRm),m=1,2,...,M,作为第m个波束区域的视频图像的第k帧雷达图像,即w'p为融合处理的加权值。
参见图7,舰船目标检测单元40(单个波束区域的海面目标检测与定位)包括:杂波功率分布统计单元401,用于对第k帧、第m个波束的雷达图像Im(χk,dAm,dRm)的海杂波功率分布进行统计分析(即同时对M个波束区域的第k帧雷达图像Im(χk,dAm,dRm)进行海杂波功率分布统计),记功率分布概率函数为fPD(Γsea),其中,Γsea为海杂波功率变量;检测阈值求取单元402,用于对舰船目标检测的阈值进行计算,设定检测的虚警率为fR,取10-2-10-4,则可以计算舰船目标检测阈值为Γm,k;舰船目标功率比较单元403,用于对第k帧、第m个波束的雷达图像Im(χk,dAm,dRm)进行舰船目标检测,即将M个波束区域的第k帧雷达图像各网格点功率与检测阈值进行功率比较,实现舰船目标检测,当Im(χk,dAm,dRm)>Γm,k时,则第k帧、第m个波束区域的雷达图像的方位向第dAm个网格点、距离向第dRm个网格点处存在舰船目标;舰船目标位置计算单元404,用于对舰船目标的位置进行计算,即将网格点位置转换成地距平面位置信息,记第k帧、第m个波束区域的雷达图像中的第j个舰船目标的地距平面位置信息为Ωm,j(χk,xj,yj);舰船目标图像切片提取单元405,用于对舰船目标的雷达图像切片进行提取,记第k帧、第m个波束的雷达图像中的第j个舰船目标的雷达图像切片为Jm,j(χk,dAj,dRj)。
在进行舰船目标视频图像形成和舰船轨迹提取时,首先根据舰船目标位置Ωm,j(χk,xj,yj),计算获得第k+1帧的天线DBF波束指向为并获得第k+1帧M个高信噪比回波信号为rm(χk+1,t)。经过上述同样处理(即图1中的20、30、40三个处理步骤)后,获得第j个舰船目标的第k+1帧位置信息Ωm,j(χk+1,xj,yj)与雷达图像切片信息Jm,j(χk+1,dAj,dRj);将第j个舰船目标的K帧位置信息Ωm,j(χk,xj,yj),k=1,2,...,K进行关联处理,获得第j个舰船目标的运动轨迹:
本发明的舰船轨迹提取方法,首先对海面舰船目标的回波信号进行同步接收,并对回波信号进行处理,形成二维回波信号;再对导航卫星散射信号进行处理,获得高信杂噪比雷达图像;并对高增益信杂噪比雷达图像进行舰船目标检测与定位,并提取舰船目标的切片信息;最后形成舰船目标的视频图像及运动轨迹。
综上所述,本发明的星载多通道GNSS-S雷达视频成像系统及舰船轨迹提取方法利用方位向多个宽波束天线,增大对海面舰船目标的观测时间,利用智能DBF处理,自主形成多个数字窄波束对海面舰船目标进行多帧成像与检测,进而获得海面舰船目标的视频图像与运动轨迹,具有长时间视频成像、批量舰船目标轨迹提取、自身隐蔽性强、系统功耗低等优势。
以上所述仅为本发明的一个实施方式而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包括在本发明的保护范围之内。
Claims (10)
1.一种星载多通道GNSS-S雷达视频成像系统,其特征在于,包括:
星载多通道GNSS-S雷达(10),用于对海面舰船目标的回波信号进行同步接收;
方位向智能DBF处理单元(20),用于对回波信号进行处理,形成二维回波信号;
多源雷达融合成像单元(30),用于对导航卫星散射信号进行处理,获得高信杂噪比雷达图像;
舰船目标检测单元(40),用于对高信杂噪比雷达图像进行舰船目标检测与定位,并提取舰船目标的切片信息;
舰船目标视频形成单元(50),用于形成舰船目标的视频图像;
舰船目标轨迹提取单元(60),用于形成舰船目标的运动轨迹。
2.根据权利要求1所述的系统,其特征在于,所述星载多通道GNSS-S雷达(10)可获得N个回波信号,记为sn(t),n=1,2,...,N;其中,N为雷达天线沿方位向的接收通道个数;
所述方位向智能DBF处理单元(20)对N个回波信号进行智能DBF处理、多源信号分离与时序恢复,形成M个数字域高增益窄波束,获得M个高信噪比回波信号,记第k帧的M个高信噪比回波信号为rm(χk,t),m=1,2,...,M;并通过多源信号分离处理使每个回波信号分离出P个回波信号,记为:
rm,p(χk,t),m=1,2,...,M,p=1,2,...,P;
再经过时序恢复形成M×P个二维回波信号,记为r′m,p(χk,τk,t);其中,χk为视频帧序号,τk为慢时序列;
所述多源雷达融合成像单元(30)对M个数字波束的P颗导航卫星散射信号进行双站雷达成像与多源雷达图像融合处理,获得第k帧的M幅高信杂噪比雷达图像Im(χk,dAm,dRm),m=1,2,...,M;其中,dAm、dRm分别为第m幅图像的方位向网格单元序列与距离向网格单元序列;
所述舰船目标检测单元(40)对第k帧的M幅高信杂噪比雷达图像进行舰船目标检测与定位,并提取舰船目标的切片信息,记第j个舰船目标的第k帧位置与图像切片分别为Ωm,j(χk,xj,yj)与Jm,j(χk,dAj,dRj);其中,xj为第j个舰船目标在地距平面的方位向位置,yj为第j个舰船目标在地距平面的距离向位置,dAj为第j个舰船目标切片沿方位向的网格单元序列,dRj为第j个舰船目标切片沿距离向的网格单元序列;
3.根据权利要求1所述的系统,其特征在于,所述星载多通道GNSS-S雷达(10)包括:
L波段方位向多通道天线(101),用于对海面舰船目标GNSS-S信号进行接收,输出N路GNSS-S信号;
N个高灵敏度数字接收机(102),用于对GNSS-S信号进行放大接收、下变频与采样量化,获得N路数字域GNSS-S信号sn(t);
GNSS-S雷达主控单元(103),用于对GNSS-S信号的放大接收、采样量化、数据存储与传输进行控制;
频率信号发生器(104),用于产生下变频所需的本振信号fLO、采样量化所需的采样时钟fADC、雷达主控所需的工作时钟fFPGA、处理单元所需的处理时钟fSP;
分布式数据存储器(105),用于对多路数字域GNSS-S信号与中间计算变量进行分布式并行存储;
阵列信号与多源信息处理单元(106),用于对多路GNSS-S信号进行智能DBF、多源雷达融合成像、舰船目标检测与定位。
4.根据权利要求1所述的系统,其特征在于,所述方位向智能DBF处理单元(20)包括:
波束指向求解单元(201),用于根据舰船目标的位置信息计算DBF波束指向,获得第k帧的M个数字窄波束指向,记第k帧的第m个波束指向为每个数字窄波束对应的合成孔径时间为Tm,m=1,2,3,…,M,Tm为1-5秒;其中,为第m个波束的方位角,θm为第m个波束的入射角;
最优权值求解单元(202),用于对智能DBF的M个数字窄波束对应的M组权值进行计算,获得M组智能DBF复数权值,记第k帧的M组DBF复数权值为wk,m,n;
多源回波信号分离单元(204),用于对M个回波信号rm(χk,t)进行多源信号分离,获得M×P个更高信噪比回波信号;
多源回波信号时序恢复单元(205),用于对M×P个回波信号的时序进行恢复,形成M×P个二维回波信号。
5.根据权利要求4所述的系统,其特征在于,所述多源回波信号分离单元(204)根据卫星导航接收机捕获跟踪获得的P颗导航卫星信号信息,分别对M个高信噪比回波信号rm(χk,t)进行互相关处理,每个波束区域均可获得P个回波信号,共可获得M×P个更高信噪比回波信号,记第k帧、第m个波束、第p颗导航卫星的回波信号为rm,p(χk,t),m=1,2,...,M,p=1,2,...,P,则有rm,p(χk,t)=R(rm(χk,t),Zp(t)),其中,R(·)为互相关函数,Zp(t)为捕获跟踪获得的第p颗导航卫星信号;
所述多源回波信号时序恢复单元(205)根据卫星导航接收机输出的P颗导航卫星位置与速度信息、低轨GNSS-S雷达天线位置与速度信息,将一维回波信号rm,p(χk,t)转换成二维回波信号r′m,p(χk,τk,t),其中,τk与t分别为方位向慢时变量与距离向快时变量。
6.根据权利要求1所述的系统,其特征在于,所述多源雷达融合成像单元(30)包括:
双站BP雷达成像处理单元(301),用于分别对第k帧的M×P个二维GNSS-S回波信号r′m,p(χk,τk,t)进行双站BP雷达成像,获得第k帧的M×P幅双站雷达图像Im,p(χk,dAm,dRm),其中,dAm与dRm分别为第m波束区域的方位向网格变量与距离向网格变量,χk为视频图像中的图像帧数,共分为M个波束区域,每个波束区域均包含P幅双站雷达图像;
多源雷达图像融合处理单元(303),用于分别对第k帧的M个波束区域的P幅双站雷达图像进行融合处理,获得第k帧的M幅高信杂噪比雷达图像Im(χk,dAm,dRm),m=1,2,...,M,则有:
其中,w'p为融合处理的加权值。
7.根据权利要求1所述的系统,其特征在于,所述舰船目标检测单元(40)包括:
杂波功率分布统计单元(401),用于对第k帧、第m个波束的雷达图像Im(χk,dAm,dRm)的海杂波功率分布进行统计分析,功率分布概率函数为fPD(Γsea),其中,Γsea为海杂波功率变量;
检测阈值求取单元(402),用于对舰船目标检测的阈值进行计算,设定虚警率为fR,取10-2-10-4,则计算检测阈值为Γm,k;
舰船目标功率比较单元(403),用于对第k帧、第m个波束的雷达图像Im(χk,dAm,dRm)进行舰船目标检测;
舰船目标位置计算单元(404),用于对舰船目标的位置进行计算,将网格点位置转换成地距平面位置信息,记第k帧、第m个波束的雷达图像中的第j个舰船目标的地距平面位置信息为Ωm,j(χk,xj,yj);
舰船目标图像切片提取单元(405),用于对舰船目标的雷达图像切片进行提取,记第k帧、第m个波束的雷达图像中的第j个舰船目标的雷达图像切片为Jm,j(χk,dAj,dRj)。
8.根据权利要求7所述的系统,其特征在于,所述舰船目标功率比较单元(403)将雷达图像各网格点功率与检测阈值进行功率比较,当Im(χk,dAm,dRm)>Γm,k时,则第k帧、第m个波束的雷达图像的方位向第dAm个网格点、距离向第dRm个网格点处存在舰船目标。
9.根据权利要求1所述的系统,其特征在于,星载多通道GNSS-S雷达的轨道高度为H,采用一幅L波段方位向多通道天线实现海面舰船目标回波信号接收,每个天线通道在方位向具有宽波束角特性;沿方位向具有N个天线通道,N取12-36,每个天线通道工作频段为L波段,中心频率为1.575GHz,带宽大于50MHz;每个天线通道在方位向的波束宽度为θA,θA取10°-15°;每个天线通道在距离向的波束宽度为θR,θR取8°-12°;天线采用水平极化对舰船目标回波信号进行接收;天线以正侧视模式对海面舰船目标回波信号进行接收,通过方位向多通道天线的数字波束形成处理,在方位向形成M个数字域高增益窄波束,实现对舰船目标的方位向探测与跟踪;
将星载多通道GNSS-S雷达的星下点轨迹划分成K个时间段,对应K帧雷达图像,帧序号记为χk;在第k个时间段内包含M个波束区域对应的M幅雷达图像,记第k帧、第m个波束区域的雷达图像为Im(χk,dAm,dRm),m=1,2,...,M;
在进行舰船目标视频图像形成和舰船轨迹提取时,根据舰船目标位置Ωm,j(χk,xj,yj),计算获得第k+1帧的天线DBF波束指向为并获得第k+1帧M个高信噪比回波信号为rm(χk+1,t);再获得第j个舰船目标的第k+1帧位置信息Ωm,j(χk+1,xj,yj)与雷达图像切片信息Jm,j(χk+1,dAj,dRj);将第j个舰船目标的K帧位置信息Ωm,j(χk,xj,yj),k=1,2,...,K进行关联处理,获得第j个舰船目标的运动轨迹:
10.一种利用权利要求1-9中任一项所述的星载多通道GNSS-S雷达视频成像系统的舰船轨迹提取方法,包括以下步骤:
a、对海面舰船目标的回波信号进行同步接收,并对回波信号进行处理,形成二维回波信号;
b、对导航卫星散射信号进行处理,获得高信杂噪比雷达图像;
c、对高信杂噪比雷达图像进行舰船目标检测与定位,并提取舰船目标的切片信息;
d、形成舰船目标的视频图像及运动轨迹。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210314900.8A CN114488134B (zh) | 2022-03-28 | 2022-03-28 | 星载多通道gnss-s雷达视频成像系统及舰船轨迹提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210314900.8A CN114488134B (zh) | 2022-03-28 | 2022-03-28 | 星载多通道gnss-s雷达视频成像系统及舰船轨迹提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114488134A true CN114488134A (zh) | 2022-05-13 |
CN114488134B CN114488134B (zh) | 2023-03-24 |
Family
ID=81488262
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210314900.8A Active CN114488134B (zh) | 2022-03-28 | 2022-03-28 | 星载多通道gnss-s雷达视频成像系统及舰船轨迹提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114488134B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115150591A (zh) * | 2022-07-27 | 2022-10-04 | 合肥金人科技有限公司 | 一种基于智能算法的区域视频监控方法 |
CN115728766A (zh) * | 2022-11-18 | 2023-03-03 | 北京卫星信息工程研究所 | 两级智能双星sar系统与海面舰船目标跟踪方法 |
CN115755051A (zh) * | 2022-11-18 | 2023-03-07 | 北京卫星信息工程研究所 | 双星高分宽幅sar信号分布式在轨处理方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101833090A (zh) * | 2010-03-12 | 2010-09-15 | 中国科学院遥感应用研究所 | 利用全球卫星定位系统信号源的机载海洋微波遥感系统 |
CN103278083A (zh) * | 2013-05-08 | 2013-09-04 | 南京信息工程大学 | Gnss-r海冰厚度探测设备及利用该设备探测海冰厚度的方法 |
CN111458711A (zh) * | 2020-04-24 | 2020-07-28 | 北京卫星信息工程研究所 | 星载双波段sar系统和舰船目标的探测方法 |
EP3822655A1 (en) * | 2019-11-13 | 2021-05-19 | Airbus Defence and Space Limited | Maritime surveillance radar |
-
2022
- 2022-03-28 CN CN202210314900.8A patent/CN114488134B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101833090A (zh) * | 2010-03-12 | 2010-09-15 | 中国科学院遥感应用研究所 | 利用全球卫星定位系统信号源的机载海洋微波遥感系统 |
CN103278083A (zh) * | 2013-05-08 | 2013-09-04 | 南京信息工程大学 | Gnss-r海冰厚度探测设备及利用该设备探测海冰厚度的方法 |
EP3822655A1 (en) * | 2019-11-13 | 2021-05-19 | Airbus Defence and Space Limited | Maritime surveillance radar |
CN111458711A (zh) * | 2020-04-24 | 2020-07-28 | 北京卫星信息工程研究所 | 星载双波段sar系统和舰船目标的探测方法 |
Non-Patent Citations (2)
Title |
---|
Z. ZHAO ET AL.: ""A novel method of ship detection by combining space-borne SAR and GNSS-R"", 《IET INTERNATIONAL RADAR CONFERENCE》 * |
郭宇华 等: ""空间目标小型相控阵雷达探测技术"", 《航天器环境工程》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115150591A (zh) * | 2022-07-27 | 2022-10-04 | 合肥金人科技有限公司 | 一种基于智能算法的区域视频监控方法 |
CN115728766A (zh) * | 2022-11-18 | 2023-03-03 | 北京卫星信息工程研究所 | 两级智能双星sar系统与海面舰船目标跟踪方法 |
CN115755051A (zh) * | 2022-11-18 | 2023-03-07 | 北京卫星信息工程研究所 | 双星高分宽幅sar信号分布式在轨处理方法及系统 |
CN115728766B (zh) * | 2022-11-18 | 2023-07-25 | 北京卫星信息工程研究所 | 两级智能双星sar系统与海面舰船目标跟踪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114488134B (zh) | 2023-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114488134B (zh) | 星载多通道gnss-s雷达视频成像系统及舰船轨迹提取方法 | |
US6867731B2 (en) | Noise suppression system and method for phased-array based systems | |
Sun et al. | Recent advances in synthetic aperture radar remote sensing—systems, data processing, and applications | |
CN114488135B (zh) | 低轨小卫星分布式gnss-s雷达系统及在轨处理方法 | |
Li et al. | BeiDou-based passive multistatic radar maritime moving target detection technique via space–time hybrid integration processing | |
CN114518577B (zh) | 星载sar与gnss-s一体化系统及协同探测方法 | |
CN111458711A (zh) | 星载双波段sar系统和舰船目标的探测方法 | |
US20220365171A1 (en) | Maritime surveillance radar | |
Li et al. | Passive radar array with low-power satellite illuminators based on fractional Fourier transform | |
CN109143231B (zh) | 基于循环对消的数字电视无源双基地雷达目标检测方法 | |
CN114545394A (zh) | 星载gnss-s雷达系统及其舰船目标探测跟踪方法 | |
Pisciottano et al. | DVB-S based passive polarimetric ISAR—Methods and experimental validation | |
CN114637004B (zh) | 星载gnss-s多站雷达在轨处理与舰船信息融合系统及方法 | |
CN113608216B (zh) | 一种星载多波段共口径sar及目标联合在轨检测系统及方法 | |
CN114660552B (zh) | 星载gnss-s雷达舰船目标信号接收与直射干扰抑制方法 | |
Frölind | Results of airborne passive SAR ground and sea target imaging using DVB-T signals | |
CN114355346B (zh) | 多星收发组网sar系统及超大幅宽海面舰船目标成像方法 | |
CN114594478B (zh) | 基于星载Ka波段SAR系统的船只目标干涉检测方法 | |
Yang et al. | Study on a novel multiple elevation beam technique for HRWS SAR system | |
Wang et al. | Radial Velocity Estimation Approach of Ship for Spaceborne Multi-Channel HRWS SAR SLC Data | |
CN115840226B (zh) | 一种方位向多通道ScanSAR快速目标检测方法 | |
Suzuki et al. | Development status of PALSAR-2 onboard ALOS-2 | |
Gogineni et al. | Sounding and imaging of fast flowing glaciers and ice-sheet margins | |
CN115728766B (zh) | 两级智能双星sar系统与海面舰船目标跟踪方法 | |
CN115728760B (zh) | 基于张量化散射信息的海面风浪流星载无源探测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |