CN104050716A - 一种海上多目标sar图像可视化建模方法 - Google Patents

一种海上多目标sar图像可视化建模方法 Download PDF

Info

Publication number
CN104050716A
CN104050716A CN201410293057.5A CN201410293057A CN104050716A CN 104050716 A CN104050716 A CN 104050716A CN 201410293057 A CN201410293057 A CN 201410293057A CN 104050716 A CN104050716 A CN 104050716A
Authority
CN
China
Prior art keywords
scene
opengl
radar
sea
coordinate system
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
CN201410293057.5A
Other languages
English (en)
Other versions
CN104050716B (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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN201410293057.5A priority Critical patent/CN104050716B/zh
Publication of CN104050716A publication Critical patent/CN104050716A/zh
Application granted granted Critical
Publication of CN104050716B publication Critical patent/CN104050716B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明一种海上多目标SAR图像可视化建模方法,在图形电磁计算(GRECO)方法的基础上,引入由海谱生成的随机海面模型,同时加载多个海上目标并设置其在整个场景中的相对位置和姿态,组成海面多目标场景,计算整个场景中各个子目标及海面对雷达回波的贡献,最终得到整个场景的雷达回波数据,并以此计算场景的合成孔径雷达(SAR)图像。其中,在计算大型海上目标(如航母等)时,采用了分区计算的方法。本发明的优点为:而可以提高电磁计算的精度;得到场景SAR图像可信度更高,具有很高的分辨率;具有较大的灵活性;快捷,成本低;具有更高的通用性。

Description

一种海上多目标SAR图像可视化建模方法
技术领域
本发明涉及一种建模方法,具体来说,是一种海上多目标SAR图像可视化建模方法。
背景技术
合成孔径雷达(SAR)是一种全天候、全天时的现代高分辨率侧视成像雷达,可以不受日照、天气等因素的影响,并且具有一定的穿透力,它是目前成像探测和遥感观测的重要手段,在很多领域得到广泛的应用。SAR成像模拟与系统仿真技术是揭示SAR成像机理、验证和设计新的成像算法的有效途径。
现今SAR成像仿真主要针对地形场景、点目标或者单一复杂目标进行仿真计算。由于海上目标所处的特定环境,海杂波和海上目标复合目标的SAR图像更为实际所需,同时,海上舰船目标的尺寸往往较大,导致获取高精度SAR成像的难度加大。
发明内容
针对上述问题,本发明提出一种海上多目标SAR图像的可视化建模方法,可更快、更便捷的精确获得海上目标的高分辨率SAR图像。
1、一种海上多目标SAR图像可视化建模方法,其特征在于:通过下述步骤实现:
步骤1:获取海面目标的三角面元模型数据文件;
步骤2:生成随机海面模型的三角面元描述的模型数据文件;
步骤3:生成海上多目标的三维场景模型数据文件;
步骤4:显示三维场景;
a、将OpenGL屏幕坐标系O2(x2、y2、z2)的x2轴方向设为屏幕上的水平向右方向,y2轴方向设为屏幕上竖直向上,z2方向设为垂直屏幕向外;
b、读入海面上多目标的三维场景模型数据文件;
c、使场景坐标系的原点为OpenGL屏幕坐标系的原点,且使场景中初始的海平面向上方向为z2方向;
步骤5:设置机载雷达运动参数,对场景进行显示;
在雷达视线的各个姿态下,以机载雷达波束照射点为中心旋转场景,使得雷达视线方向保持与OpenGL屏幕坐标系中的-Z2方向一致;根据机载雷达波束照射的范围,通过OpenGL中的窗口函数设置OpenGL视场中的目标场景;
步骤6:计算在机载雷达运动期间,雷达视线的各个姿态下OpenGL视场中的目标场景相对雷达视线的姿态角序列;进而得到目标场景显示序列,并依次通过步骤7~9,得到雷达视线的各个姿态对应的雷达回波,进而得到机载雷达运动期间的雷达回波序列;
步骤7:对目标场景进行分区显示;
a、获取目标场景在OpenGL屏幕坐标系x2、y2轴方向的矩形投影尺寸,分别为
b、令给定的像素分辨尺寸为l0,对矩形投影在OpenGL屏幕坐标系x2、y2轴方向分别进行划分,得到的划分区间数量分别为:
M x 2 = ROUND ( ( D max x 2 + D min x 2 ) / l 0 ) Q y 2 = ROUND ( ( L max y 2 + L min y 2 ) / l 0 ) ; 其中,ROUND表示四舍五入;则投影矩形总共的分区数为对各个分区进行单独的投影显示;
步骤8:计算各个分区的散射场;
步骤9:将各分区的散射场矢量叠加,得到总雷达回波数据。
步骤10:获得场景的SAR图像;
本发明采用经典的雷达成像算法:距离多普勒成像算法,对机载雷达运动期间的雷达回波序列进行处理,生成海上多目标的高分辨率雷达图像。
本发明的优点在于:
1、本发明海上多目标SAR图像可视化建模方法,利用分区计算大型雷达目标散射特性功能,计算精度更高;
2、本发明海上多目标SAR图像可视化建模方法,利用海谱生成随机海面模型,得到场景SAR图像可信度(逼真度)更高;
3、本发明海上多目标SAR图像可视化建模方法,通过改变每个子目标相对位置及姿态,就可生成不同场景布置的雷达图像,具有较大的灵活性;
4、本发明海上多目标SAR图像可视化建模方法,生成的场景雷达图像具有很高的分辨率;
5、本发明海上多目标SAR图像可视化建模方法,利用计算机场景仿真比实测得到的海面多目标场景的SAR图像数据更为快捷,成本更低;
6、本发明海上多目标SAR图像可视化建模方法,采用三角面元格式文件作为目标的模型数据文件,具有更高的通用性。
附图说明
图1为本发明海上多目标SAR图像可视化建模方法流程图;
图2为本发明整体场景与机载雷达相对关系示意图;
图3为实施例中最终生成的海上多目标场景SAR图像。
具体实施方式
下面结合附图对本发明作进一步的详细说明。
本发明海上多目标SAR图像可视化建模方法,如图1所示,通过下述步骤实现:
步骤1:获取海面目标的三角面元模型数据文件;
通过建模软件(如:AutoCAD、3Dmax等)分别对各个海上目标进行三维几何建模,生成三角面元描述的模型数据文件(如:DXF、MSH等)。
步骤2:生成随机海面模型的三角面元描述的模型数据文件;
海面上海浪的运动是一种复杂的随机过程;在海洋学中,利用谱以随机过程来描述海浪是进行海浪研究的主要途径之一,可将无限个随机的余弦波进行叠加来描述一个定点的波面η(t),则有:
η ( t ) = Σ n = 1 ∞ a n cos ( ω n t + ϵ n ) - - - ( 1 )
式(1)中,an与ωn分别为第n个余弦波的振幅与角频率,n∈[1~∞];εn为第n个余弦波在0~2π范围内均布的随机初位相;t为时刻。
但式(1)仅能描述固定点的波面随时间的变化,而实际的海面是三维的,其能量不仅分布在一定的频率范围内,还分布在相当广的方向范围内,因此通过式(2),可得到随时间以及空间同时变化的固定点波面。
η ( x , y , t ) = Σ n = 1 ∞ a n cos [ ( ω n t - k n ( x cos θ n + y sin θ n ) + ϵ n ) ] - - - ( 2 )
式(2)中,x、y分别为二维海面坐标系O(x、y)的横纵坐标;第n个余弦波传播方向与二维海面坐标系中x轴夹角为θn,-π≤θn≤π;kn为波数。对深水波来说,根据线性波浪理论,
式(2)中,各个波的振幅an,可以由如式(3)中给出的海浪谱求出:
S ( ω , θ ) dωdθ = Σ Δω Σ Δθ 1 2 a n 2 - - - ( 3 )
式(3)中,海浪谱S(ω,θ)的定义如下:
S(ω,θ)=S(ω)G(ω,θ)   (4)
式(4)中,S(ω)为频谱,G(ω,θ)为方向分布函数,简称方向函数,Δω为角频率间隔;Δθ余弦波传播方向与二维海面坐标系中x轴夹角间隔;
式(4)中,海浪频谱通过Pierson-Moscowitz谱,简称P-M谱来模拟:
S ( ω ) = a g 2 ω 5 exp [ - β ( g Uω ) 4 ] - - - ( 5 )
式(5)中,a=8.1×10-3,为无因次常数,β=0.74,g为重力加速度,U为海面上19.5m高处的风速。
式(4)中,海浪方向分布函数通过波浪立体观测计划(简称为SWOP—Stereo Wave Observation Project)得到:
G ( ω , θ ) = 1 π ( 1 + p cos 2 θ + q cos 4 θ ) - - - ( 6 )
式(6)中, p = 0.5 + 0.82 exp [ - 1 2 ( ω ω m ) 4 ] , q = 0.32 exp [ - 1 2 ( ω ω m ) 4 ] , | θ | ≤ π 2 ωm处的谱峰频率。
将通过上述方法生成的海面模型η(x,y,t)存储为三角面元数据文件。
步骤3:生成真实海上多目标的三维场景模型数据文件;
由建模计算机读入步骤1中得到的海上目标的三角面元描述的模型数据文件,与步骤2中生成的海面模型的三角面元描述的模型数据文件,构成整个三维场景。将海面看成场景中的一个特殊目标,并通过参数输入的形式设定海面在场景坐标系O1(x1、y1、z1)中的空间位置和法矢。每个海上目标均具有一个初始姿态;此时,在场景坐标系中,以海面的几何中心为原点,改变各个海上目标的初始姿态至所需初始姿态。随后在场景坐标系中设置各个海上目标在场景坐标系中的位置。由于海上目标基本属于舰船类,按照海上目标实际状态,将海上目标的位置和姿态调整至漂浮在海面上的状态,使得场景符合实际雷达照射场景状态;最后得到模拟真实海面上多目标的三维场景模型数据文件。
步骤4:显示三维场景;
a、设置OpenGL屏幕窗口大小为1024*768像素;同时,将OpenGL屏幕坐标系O2(x2、y2、z2)的x2轴方向设为屏幕上的水平向右方向,y2轴方向设为屏幕上竖直向上,z2方向设为垂直屏幕向外;
b、通过OpenGL的API函数读入步骤3中得到的模拟真实海面上多目标的三维场景模型数据文件;
c、设置三维场景在OpenGL屏幕窗口中的姿态,使场景坐标系的原点为OpenGL屏幕坐标系的原点,且使场景中初始的海平面向上方向为z2方向。
步骤5:设置真实的机载雷达运动参数(包括起点、终点与速度),对场景进行显示;
如图2所示,机载雷达从P1点运动至P2点,随着机载雷达的不断运动,场景相对机载雷达视线的姿态也不断变化,在机载雷达视线的各个姿态下,以机载雷达波束照射点为中心旋转场景,使得机载雷达视线方向保持与OpenGL屏幕坐标系中的-Z2方向一致,即旋转场景坐标系至屏幕坐标系位置,进而保证机载雷达视线的可视目标作为OpenGL视场中的目标场景。则根据机载雷达波束照射的范围,通过OpenGL中的窗口函数设置OpenGL视场中的目标场景。同时,计算在机载雷达运动期间,雷达视线的各个姿态下OpenGL视场中的目标场景相对雷达视线的姿态角序列;进而得到目标场景显示序列,并依次通过步骤7~9,得到雷达视线的各个姿态对应的雷达回波,进而得到机载雷达运动期间的雷达回波序列。
步骤6:计算在机载雷达运动期间,雷达视线的各个姿态下OpenGL视场中的目标场景相对雷达视线的姿态角序列;进而得到目标场景显示序列,并依次通过步骤7~9,得到雷达视线的各个姿态对应的雷达回波,进而得到机载雷达运动期间的雷达回波序列。
步骤7:对OpenGL视场中的目标场景进行分区显示,用于扩充后续图形电磁计算(GRECO)的计算范围;
由于海面上多目标的电尺寸较大,为提高雷达回波的计算精度,通过分区显示法将OpenGL视场中的目标场景进行分区,分别利用OpenGL函数单独显示,具体方式为:
a、获取目标场景在OpenGL屏幕坐标系x2、y2轴方向的矩形投影尺寸,分别为
b、令给定的像素分辨尺寸为l0,对矩形投影在OpenGL屏幕坐标系x2、y2轴方向分别进行划分,得到的划分区间数量分别为:
M x 2 = ROUND ( ( D max x 2 + D min x 2 ) / l 0 ) Q y 2 = ROUND ( ( L max y 2 + L min y 2 ) / l 0 ) ; 其中,ROUND表示四舍五入;则投影矩形总共的分区数为对各个分区进行单独的投影显示。
步骤8:按图形电磁计算(GRECO)的方法来计算各个分区的散射场;
A、计算投影矩形的像素尺寸;
OpenGL屏幕坐标系x2、y2轴方向的像素尺寸分别为:
Δ x 2 = D max x 2 + D min x 2 M x 2 P x 2 - - - ( 7 )
Δ y 2 = L max y 2 + L min y 2 Q y 2 P y 2 - - - ( 8 )
式(7)与式(8)中,分别是每个分区在OpenGL屏幕坐标系x2、y2轴方向的像素数。
B、利用OpenGL的API函数按各分区的投影比例(分区实际大小与屏幕分辨率的比值)分别进行显示,使每个分区显示时屏幕像素代表的电尺寸具有合适的大小;
C、计算各分区像各素位置信息与法矢信息,具体为:
通过OpenGL的窗口大小和OpenGL的API函数获得屏幕上各像素的位置信息和颜色信息;
各像素的位置信息:
令屏幕像素为个,OpenGL屏幕窗口的横向大小为纵向大小为则第m,n个像素的x2,y2大小为 x 2 m , n = - L min x 2 + L max x 2 + L min x 2 M x 2 m , y 2 m , n = - L min y 2 + L max y 2 + L min y 2 N y 2 n , 第m,n个像素的z2值由OpenGL的z-buffer提供。
各像素颜色信息:
利用OpenGL的光照模型,使红、绿、蓝光源分别位于OpenGL屏幕坐标系的x2,y2,z2轴,由OpenGL的颜色buffer可获得每个像素颜色的(R,G,B),分别对应于红绿蓝三种颜色。
通过每个像素位置信息和颜色信息,可分别获取每个像素位置所代表的面元三维坐标和法矢朝向。其中,颜色信息与法矢朝向间的对应关系为:
令物体表面是漫反射并且每种颜色的漫反射系数为1,光源的强度为1,光传播在媒质中的衰减为1,则光照模型为:
L d = n ^ · r ^ i , d = R G , B , - - - ( 9 )
式中,Ld为每种颜色的漫反射亮度,是像素代表的面元法矢,是指向光源的方向矢量。则像素代表的面元法矢的子分量与颜色的对应关系为:
RED : r ^ i = x ^ 2 ⇒ R = n ^ · x ^ 2 = n x 2
GREEN : r ^ i = y ^ 2 ⇒ G = n ^ · y ^ 2 = n y 2
BLUE : r ^ i = z ^ 2 ⇒ B = n ^ · z ^ 2 = n z 2
分别为屏幕坐标系中x2、y2、z2方向的矢量。
由于颜色分量总为正值,而法矢量存在正负问题。在OpenGL中,当面元的法矢量指向光源时,面元会显示相应的颜色;反之,当其法矢量背离光源时,面元将不显示。这样,为确定法矢的正负,必须从正负坐标轴用红绿蓝三种光源分别照射目标分区,显示颜色的一面即为法矢方向。
D、根据得到的分l0区各像素位置信息与法矢信息,并对每个子区域利用高频近似理论中的物理光学计算理论及物理绕射计算理论分区计算其散射场计算分区内所有有效像素的散射场。
则散射场SPO为:
S PO = ∫ s ( n ^ · z ^ 2 ) e j 2 k z 2 ds ≈ Σ PIXELS e j 2 k z 2 Δ s ′ - - - ( 10 )
其中,k为雷达波数,Δs′为每个像素代表的实际投影面积。
步骤9:将各分区的散射场矢量叠加,得到总雷达回波数据。
步骤10:获得场景的SAR图像。
本发明采用经典的雷达成像算法:距离多普勒成像算法,对机载雷达运动期间的雷达回波序列进行处理,生成海上多目标的高分辨率雷达图像。
实施例:
本实例中,海上目标以两艘舰船目标为例,海面波浪以10m/s风速为例,利用高分辨率聚束式机载SAR回波仿真模型对场景进行回波仿真,工作方式为正侧视,雷达发射波采用X波段的线性调频波脉冲,雷达工作参数如下:
中心频率:10GHz
带宽:0.4GHz
飞行距离:200米
采样点数:1024点
生成分辨率分别为0.2米、0.2米的海面及2艘舰船目标的SAR图像。如图3所示,从图中可以看出两艘舰船的基本外形及海面分布。实验结果证明了本发明对海面多目标场景进行SAR成像建模的可行性。

Claims (5)

1.一种海上多目标SAR图像可视化建模方法,其特征在于:通过下述步骤实现:
步骤1:获取海面目标的三角面元模型数据文件;
步骤2:生成随机海面模型的三角面元描述的模型数据文件;
步骤3:生成海上多目标的三维场景模型数据文件;
步骤4:显示三维场景;
a、将OpenGL屏幕坐标系O2(x2、y2、z2)的x2轴方向设为屏幕上的水平向右方向,y2轴方向设为屏幕上竖直向上,z2方向设为垂直屏幕向外;
b、读入海面上多目标的三维场景模型数据文件;
c、使场景坐标系的原点为OpenGL屏幕坐标系的原点,且使场景中初始的海平面向上方向为z2方向;
步骤5:设置机载雷达运动参数,对场景进行显示;
在雷达视线的各个姿态下,以机载雷达波束照射点为中心旋转场景,使得雷达视线方向保持与OpenGL屏幕坐标系中的-Z方向一致;根据机载雷达波束照射的范围,通过OpenGL中的窗口函数设置OpenGL视场中的目标场景;
步骤6:计算在机载雷达运动期间,雷达视线的各个姿态下OpenGL视场中的目标场景相对雷达视线的姿态角序列;进而得到目标场景显示序列,并依次通过步骤7~9,得到雷达视线的各个姿态对应的雷达回波,进而得到机载雷达运动期间的雷达回波序列;
步骤7:对目标场景进行分区显示;
a、获取目标场景在OpenGL屏幕坐标系x2、y2轴方向的矩形投影尺寸,分别为
b、令给定的像素分辨尺寸为l0,对矩形投影在OpenGL屏幕坐标系x2、y2轴方向分别进行划分,得到的划分区间数量分别为:
M x 2 = ROUND ( ( D max x 2 + D min x 2 ) / l 0 ) Q y 2 = ROUND ( ( L max y 2 + L min y 2 ) / l 0 ) ; 其中,ROUND表示四舍五入;则投影矩形总共的分区数为对各个分区进行单独的投影显示;
步骤8:计算各个分区的散射场;
步骤9:将各分区的散射场矢量叠加,得到总雷达回波数据。
步骤10:获得场景的SAR图像。
本发明采用经典的雷达成像算法:距离多普勒成像算法,对机载雷达运动期间的雷达回波序列进行处理,生成海上多目标的高分辨率雷达图像。
2.如权利要求1所述一种海上多目标SAR图像可视化建模方法,其特征在于:所述步骤2具体实现方法为:
将无限个随机的余弦波进行叠加来描述一个定点的波面η(t),则有:
η ( t ) = Σ n = 1 ∞ a n cos ( ω n t + ϵ n ) - - - ( 1 )
式(1)中,an与ωn分别为第n个余弦波的振幅与角频率,n∈[1~∞];εn为第n个余弦波在0~2π范围内均布的随机初位相;t为时刻;
通过式(2),可得到随时间以及空间同时变化的固定点波面;
η ( x , y , t ) = Σ n = 1 ∞ a n cos [ ( ω n t - k n ( x cos θ n + y sin θ n ) + ϵ n ) ] - - - ( 2 )
式(2)中,x、y分别为二维海面坐标系O(x、y)的横纵坐标;第n个余弦波传播方向与二维海面坐标系中x轴夹角为θn,-π≤θn≤π;kn为波数。对深水波来说,
式(2)中,各个波的振幅an,可以由如式(3)中给出的海浪谱求出:
S ( ω , θ ) dωdθ = Σ Δω Σ Δθ 1 2 a n 2 - - - ( 3 )
式(3)中,海浪谱S(ω,θ)的定义如下:
S(ω,θ)=S(ω)G(ω,θ)   (4)
式(4)中,S(ω)为频谱,G(ω,θ)为方向分布函数,简称方向函数,Δω为角频率间隔;Δθ余弦波传播方向与二维海面坐标系中x轴夹角间隔;
式(4)中,海浪频谱通过Pierson-Moscowitz谱来模拟:
S ( ω ) = a g 2 ω 5 exp [ - β ( g Uω ) 4 ] - - - ( 5 )
式(5)中,a=8.1×10-3,为无因次常数,β=0.74,g为重力加速度,U为海面上19.5m高处的风速;
式(4)中,海浪方向分布函数通过波浪立体观测计划得到:
G ( ω , θ ) = 1 π ( 1 + p cos 2 θ + q cos 4 θ ) - - - ( 6 )
式(6)中, p = 0.5 + 0.82 exp [ - 1 2 ( ω ω m ) 4 ] , q = 0.32 exp [ - 1 2 ( ω ω m ) 4 ] , | θ | ≤ π 2 ωm处的谱峰频率;
将通过上述方法生成的海面模型η(x,y,t)存储为三角面元数据文件。
3.如权利要求1所述一种海上多目标SAR图像可视化建模方法,其特征在于:所述步骤3的实现方法为:由建模计算机读入海上目标的三角面元描述的模型数据文件,与海面模型的三角面元描述的模型数据文件,构成整个三维场景;在场景坐标系O1(x1、y1、z1)中,以海面的几何中心为原点,改变各个海上目标的初始姿态至所需初始姿态;随后按照海上目标实际状态,将海上目标的位置和姿态调整至漂浮在海面上的状态,得到模拟真实海面上多目标的三维场景模型数据文件。
4.如权利要求1所述一种海上多目标SAR图像可视化建模方法,其特征在于:步骤8实现方法为:
A、计算投影矩形的像素尺寸;
OpenGL屏幕坐标系x2、y2轴方向的像素尺寸分别为:
Δ x 2 = D max x 2 + D min x 2 M x 2 P x 2 - - - ( 7 )
Δ y 2 = L max y 2 + L min y 2 Q y 2 P y 2 - - - ( 8 )
式(7)与式(8)中,分别是每个分区在OpenGL屏幕坐标系x2、y2轴方向的像素数;
B、按各分区的投影比例分别进行显示;
C、计算各分区像各素位置信息与法矢信息,具体为:
获得屏幕上各像素的位置信息和颜色信息;
各像素的位置信息:
令屏幕像素为个,OpenGL屏幕窗口的横向大小为纵向大小为则第m,n个像素的x2,y2大小为 x 2 m , n = - L min x 2 + L max x 2 + L min x 2 M x 2 m , y 2 m , n = - L min y 2 + L max y 2 + L min y 2 N y 2 n , 第m,n个像素的z2值由OpenGL的z-buffer提供;
各像素颜色信息:
利用OpenGL的光照模型,使红、绿、蓝光源分别位于OpenGL屏幕坐标系的x2,y2,z2轴,由OpenGL的颜色buffer可获得每个像素颜色的(R,G,B),分别对应于红绿蓝三种颜色;
通过每个像素位置信息和颜色信息,可分别获取每个像素位置所代表的面元三维坐标和法矢朝向。其中,颜色信息与法矢朝向间的对应关系为:
令物体表面是漫反射并且每种颜色的漫反射系数为1,光源的强度为1,光传播在媒质中的衰减为1,则光照模型为:
L d = n ^ · r ^ i , d = R G , B , - - - ( 9 )
式中,Ld为每种颜色的漫反射亮度,是像素代表的面元法矢,是指向光源的方向矢量。则像素代表的面元法矢的子分量与颜色的对应关系为:
RED : r ^ i = x ^ 2 ⇒ R = n ^ · x ^ 2 = n x 2
GREEN : r ^ i = y ^ 2 ⇒ G = n ^ · y ^ 2 = n y 2
BLUE : r ^ i = z ^ 2 ⇒ B = n ^ · z ^ 2 = n z 2
分别为屏幕坐标系中x2、y2、z2方向的矢量;
D、根据得到的分区各像素位置信息与法矢信息,并对每个子区域利用高频近似理论中的物理光学计算理论及物理绕射计算理论分区计算其散射场计算分区内所有有效像素的散射场。
则散射场SPO为:
S PO = ∫ s ( n ^ · z ^ 2 ) e j 2 k z 2 ds ≈ Σ PIXELS e j 2 k z 2 Δ s ′ - - - ( 10 )
其中,k为雷达波数,Δs′为每个像素代表的实际投影面积。
5.如权利要求4所述一种海上多目标SAR图像可视化建模方法,其特征在于:步骤C中,为确定法矢的正负,须从正负坐标轴用红绿蓝三种光源分别照射目标分区,显示颜色的一面为法矢方向。
CN201410293057.5A 2014-06-25 2014-06-25 一种海上多目标sar图像可视化建模方法 Expired - Fee Related CN104050716B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410293057.5A CN104050716B (zh) 2014-06-25 2014-06-25 一种海上多目标sar图像可视化建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410293057.5A CN104050716B (zh) 2014-06-25 2014-06-25 一种海上多目标sar图像可视化建模方法

Publications (2)

Publication Number Publication Date
CN104050716A true CN104050716A (zh) 2014-09-17
CN104050716B CN104050716B (zh) 2017-06-16

Family

ID=51503484

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410293057.5A Expired - Fee Related CN104050716B (zh) 2014-06-25 2014-06-25 一种海上多目标sar图像可视化建模方法

Country Status (1)

Country Link
CN (1) CN104050716B (zh)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105356946A (zh) * 2015-09-29 2016-02-24 浙江大学 一种基于空中平台的空中-水下通信系统
CN106682391A (zh) * 2016-11-23 2017-05-17 大连理工大学 基于p‑m谱和jonswap谱的实测风浪情况筛选和理论频谱拟合方法
CN106680812A (zh) * 2016-10-27 2017-05-17 西安空间无线电技术研究所 一种基于解析面元的微波关联成像仿真方法
CN108519075A (zh) * 2018-03-22 2018-09-11 北京控制工程研究所 一种空间多目标位姿测量方法
CN109064544A (zh) * 2018-08-09 2018-12-21 太平洋未来科技(深圳)有限公司 全景视频中虚拟对象的光影渲染方法、装置及电子设备
CN109061571A (zh) * 2018-06-25 2018-12-21 北京无线电测量研究所 一种基于OpenGL的雷达信息显示方法和系统
CN110223335A (zh) * 2019-06-06 2019-09-10 西北工业大学 一种基于视图划分的高精度greco方法
US20190369231A1 (en) * 2018-06-01 2019-12-05 Thales Method and system for determining a characteristic dimension of a ship
CN111025246A (zh) * 2019-11-28 2020-04-17 北京遥测技术研究所 一种静止轨道sar对海面及船只复合场景成像仿真系统及方法
CN112285711A (zh) * 2020-09-28 2021-01-29 西安航空学院 一种基于深度置信网络的三维sar图像目标识别方法
CN112684469A (zh) * 2021-01-14 2021-04-20 江苏恒澄交科信息科技股份有限公司 基于海事雷达图像的航道特征方向识别方法及系统
CN112764030A (zh) * 2020-12-25 2021-05-07 中国人民解放军空军工程大学 基于电磁与统计混合模型的海面sar图像杂波生成方法
CN112816794A (zh) * 2021-01-04 2021-05-18 西安电子科技大学 天线波束照射下近场动态海面电磁波强度分布可视化方法
CN113111758A (zh) * 2021-04-06 2021-07-13 中山大学 一种基于脉冲神经网络的sar图像舰船目标识别方法
CN113238201A (zh) * 2021-05-26 2021-08-10 深圳成谷科技有限公司 一种超分辨雷达定位方法、系统、设备及存储介质

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5192208A (en) * 1989-08-21 1993-03-09 General Electric Company Radar simulation for use with a visual simulator
CN101281249A (zh) * 2008-05-20 2008-10-08 北京航空航天大学 一种高分辨率合成孔径雷达的目标散射特性建模方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LIN YIMING等: "A novel 3D visualization method of SAR data", 《RADAR CONFERENCE 2013, IET INTERNATIONAL》 *
刘鹏 等: "基于谱的三维随机海浪数值模拟研究", 《计算机与数字工程》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105356946B (zh) * 2015-09-29 2018-01-30 浙江大学 一种基于空中平台的空中‑水下通信系统
CN105356946A (zh) * 2015-09-29 2016-02-24 浙江大学 一种基于空中平台的空中-水下通信系统
CN106680812A (zh) * 2016-10-27 2017-05-17 西安空间无线电技术研究所 一种基于解析面元的微波关联成像仿真方法
CN106680812B (zh) * 2016-10-27 2019-04-09 西安空间无线电技术研究所 一种基于解析面元的微波关联成像仿真方法
CN106682391A (zh) * 2016-11-23 2017-05-17 大连理工大学 基于p‑m谱和jonswap谱的实测风浪情况筛选和理论频谱拟合方法
CN106682391B (zh) * 2016-11-23 2019-04-16 大连理工大学 基于p-m谱和jonswap谱的实测风浪情况筛选和理论频谱拟合方法
CN108519075A (zh) * 2018-03-22 2018-09-11 北京控制工程研究所 一种空间多目标位姿测量方法
US20190369231A1 (en) * 2018-06-01 2019-12-05 Thales Method and system for determining a characteristic dimension of a ship
US11635511B2 (en) * 2018-06-01 2023-04-25 Thales Method and system for determining a characteristic dimension of a ship
CN109061571A (zh) * 2018-06-25 2018-12-21 北京无线电测量研究所 一种基于OpenGL的雷达信息显示方法和系统
CN109064544A (zh) * 2018-08-09 2018-12-21 太平洋未来科技(深圳)有限公司 全景视频中虚拟对象的光影渲染方法、装置及电子设备
CN110223335A (zh) * 2019-06-06 2019-09-10 西北工业大学 一种基于视图划分的高精度greco方法
CN111025246A (zh) * 2019-11-28 2020-04-17 北京遥测技术研究所 一种静止轨道sar对海面及船只复合场景成像仿真系统及方法
CN112285711A (zh) * 2020-09-28 2021-01-29 西安航空学院 一种基于深度置信网络的三维sar图像目标识别方法
CN112764030A (zh) * 2020-12-25 2021-05-07 中国人民解放军空军工程大学 基于电磁与统计混合模型的海面sar图像杂波生成方法
CN112764030B (zh) * 2020-12-25 2023-09-08 中国人民解放军空军工程大学 基于电磁与统计混合模型的海面sar图像杂波生成方法
CN112816794A (zh) * 2021-01-04 2021-05-18 西安电子科技大学 天线波束照射下近场动态海面电磁波强度分布可视化方法
CN112684469A (zh) * 2021-01-14 2021-04-20 江苏恒澄交科信息科技股份有限公司 基于海事雷达图像的航道特征方向识别方法及系统
CN113111758A (zh) * 2021-04-06 2021-07-13 中山大学 一种基于脉冲神经网络的sar图像舰船目标识别方法
CN113111758B (zh) * 2021-04-06 2024-01-12 中山大学 一种基于脉冲神经网络的sar图像舰船目标识别方法
CN113238201A (zh) * 2021-05-26 2021-08-10 深圳成谷科技有限公司 一种超分辨雷达定位方法、系统、设备及存储介质
CN113238201B (zh) * 2021-05-26 2024-06-28 深圳成谷科技有限公司 一种超分辨雷达定位方法、系统、设备及存储介质

Also Published As

Publication number Publication date
CN104050716B (zh) 2017-06-16

Similar Documents

Publication Publication Date Title
CN104050716A (zh) 一种海上多目标sar图像可视化建模方法
CN106772290B (zh) 一种全极化双基地雷达目标动态回波模拟方法
CN103207387B (zh) 一种机载相控阵pd雷达杂波的快速模拟方法
CN101281249A (zh) 一种高分辨率合成孔径雷达的目标散射特性建模方法
CN106324571A (zh) 基于正向法的仿三维场景sar雷达回波的快速实现方法
CN102788979B (zh) 一种基于后向投影InSAR成像配准的GPU实现方法
CN114047511B (zh) 一种基于csa算法的时变海面机载sar成像仿真方法
WO2015160289A1 (en) Method and system for rendering a synthetic aperture radar image
CN105260524B (zh) 一种扫描状态下船舶导航雷达二维回波序列像仿真方法
CN103018740B (zh) 一种基于曲面投影的InSAR成像方法
CN105445711A (zh) 一种基于逆Omega-K算法的海面要素SAR原始数据仿真方法
CN109541591B (zh) 一种基于线性滤波法的sar海浪成像仿真方法
CN107607951A (zh) 一种sar图像二次散射特征模拟方法
CN109188384B (zh) 空间目标回波动态观测的电磁仿真方法
CN103616682A (zh) 一种基于曲面投影的多基线InSAR处理方法
El Natour et al. Radar and vision sensors calibration for outdoor 3D reconstruction
Willis et al. Hardware-accelerated SAR simulation with NVIDIA-RTX technology
CN113447896B (zh) 一种基于动态遮挡判断的起伏地形sar回波仿真方法
CN104101869A (zh) 一种极坐标下的地波雷达运动目标仿真建模方法
Dogan et al. Time domain SAR raw data simulation of distributed targets
CN108562899B (zh) 高分辨极化sar目标图像快速仿真方法
Dumont et al. DIONISOS. An end-to-end SAR simulator
CN103996214B (zh) Bp-sar图像重建误差分析方法及bp-sar图像重建方法
Perälä et al. ESPRESS—On efficient bistatic characterization of radar targets
Lei et al. SAR-NeRF: Neural Radiance Fields for Synthetic Aperture Radar Multi-View Representation

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170616

Termination date: 20200625

CF01 Termination of patent right due to non-payment of annual fee