CN115661293A - 一种电磁散射中目标瞬态感应特征图像生成方法 - Google Patents
一种电磁散射中目标瞬态感应特征图像生成方法 Download PDFInfo
- Publication number
- CN115661293A CN115661293A CN202211339858.1A CN202211339858A CN115661293A CN 115661293 A CN115661293 A CN 115661293A CN 202211339858 A CN202211339858 A CN 202211339858A CN 115661293 A CN115661293 A CN 115661293A
- Authority
- CN
- China
- Prior art keywords
- grid
- curved surface
- point
- transient
- target
- 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
- 230000001052 transient effect Effects 0.000 title claims abstract description 38
- 230000006698 induction Effects 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000013507 mapping Methods 0.000 claims abstract description 12
- 230000003287 optical effect Effects 0.000 claims abstract description 11
- 230000011218 segmentation Effects 0.000 claims abstract description 4
- 230000008569 process Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 2
- 230000010287 polarization Effects 0.000 claims description 2
- 230000005284 excitation Effects 0.000 abstract description 6
- 230000000007 visual effect Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 5
- 238000012512 characterization method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种电磁散射中目标瞬态感应特征图像生成方法,包括以下步骤:S1.给定规则曲面受脉冲平面电磁波辐射的场景信息;S2.将规则曲面在XOY平面进行投影并进行网格生成,得到平面二维离散矩形网格;S3.遍历平面二维离散矩形网格中的每一个网格,得到每一个网格的中心点坐标和在曲面上的映射点,计算平面波波前从O点到达曲面的最小和最大光程,并进行区间分割得到时间序列和幅度序列;S4.对于二维离散矩形网格中的每一个网格,确定网格中心点在规则曲面上的映射点信息列表;S5.生成电磁散射中目标瞬态感应特征图像。本发明通过数字图像表征目标在瞬态脉冲平面波激励条件的感应场特征图像,能够实现目标在瞬态脉冲平面波下感应场特征的可视化展示。
Description
技术领域
本发明涉及电磁散射,特别是涉及一种电磁散射中目标瞬态感应特征图像生成方法。
背景技术
随着短脉冲通信和超宽带雷达系统的广泛应用,采用时域方法来分析目标的瞬态电磁散射已经成为计算电磁学中的热点之一。了能够快速准确分析电大尺寸目标的瞬态散射,近年来发展出与各种频域高频算法相对应的时域方法。其中,时域物理光学是物理光学法的时域版本,其主要思路是通过傅里叶变换将频域物理光学公式转化为时域公式,然后利用把目标表面的时域等效电流密度代入到时域散射电场表达式中,得到散射远场时域分布。然而众多学者从瞬态电磁散射的结果即目标的最终雷达散射截面的快速计算层面,而鲜有对瞬态电磁散射过程中平面波与散射体的时域作用瞬态过程进行分析。
发明内容
本发明的目的在于克服现有技术的不足,提供一种电磁散射中目标瞬态感应特征图像生成方法,采用数字采样的方式将目标表面的时域响应信号离散成时间与空间的二维序列,通过数字图像表征了目标在瞬态脉冲平面波激励条件下的感应场特征图像,能够实现目标在瞬态脉冲平面波下感应场特征的可视化展示。
本发明的目的是通过以下技术方案来实现的:S1.给定规则曲面受脉冲平面电磁波辐射的场景信息;
S2.将规则曲面在XOY平面进行投影,得到投影的矩形区域RXOY,并在矩形区域RXOY内进行网格生成,得到平面二维离散矩形网格;
S3.遍历平面二维离散矩形网格中的每一个网格,得到每一个网格的中心点坐标和在曲面上的映射点,计算平面波波前从O点到达曲面的最小和最大光程,并进行区间分割得到时间序列和幅度序列;
S5.生成电磁散射中目标瞬态感应特征图像。
本发明的有益效果是:
本专利给出了目标在瞬态脉冲平面波激励条件下的感应场表达式,该表达式相比时谐频域表达式的最大不同是采用了时域卷积的形式,而后采用数字采样的方式将目标表面的时域响应信号离散成时间与空间的二维序列,通过数字图像表征了目标在瞬态脉冲平面波激励条件下的感应场特征图像,能够实现目标在瞬态脉冲平面波下感应场特征的可视化展示。
附图说明
图1为本发明的方法流程图;
图2为本发明的辐射场景示意图;
图3为规则曲面投影示意图;
图4为网格剖分示意图;
图5为平面波波前从O点到达曲面的最小和最大光程示意图;
图6为实施例中得到的Iz示意图;
图7为实施例中得到的Ix示意图;
图8为实施例中得到的Ity示意图;
图9为实施例中得到的Itz示意图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
如图1所示,一种电磁散射中目标瞬态感应特征图像生成方法,其特征在于:包括以下步骤:
S1.给定规则曲面受脉冲平面电磁波辐射的场景信息:如图2所示,由脉冲平面电磁波照射到曲面表面,产生电流感应,设所述的场景信息包括:
(2)规则曲面的表征公式为z=f(x,y),参数范围为:x1≤x≤x2和y1≤y≤y2;
(3)网格剖分密度为D,在X轴和Y轴上的网格剖分长度Δx=Δy=λ/D,其中λ=c0/f0为工作频率对应的波长,c0为自由空间波速。
S2.将规则曲面在XOY平面进行投影,得到投影的矩形区域RXOY,如图3所示;然后在矩形区域RXOY内进行网格生成,得到平面二维离散矩形网格,如图4所示;
在矩形区域RXOY内进行网格生成,得到平面二维离散矩形网格的过程包括:
根据设置的网格密度和得到的X轴网格剖分长度Δx和Y轴网格剖分长度Δy,计算X轴剖分数nX和Y轴剖分数nY:
S3.遍历平面二维离散矩形网格中的每一个网格,得到每一个网格的中心点坐标和在曲面上的映射点,计算平面波波前从O点到达曲面的最小和最大光程,并进行区间分割得到时间序列和幅度序列;
所述步骤S3包括:
S301:遍历X轴上的(nX-1)个网格和Y轴上的(nY-1)个网格,得到X轴上的第i个,Y轴上的第j个网格Gij,Gij网格中心点二维坐标Rij(xi,yj)为:
xi=(i+0.5)Δx,yj=(j+0.5)Δy
S302:网格中心点二维坐标(xi,yi)代入规则曲面表征公式得到该网格中心点在曲面上的映射点Rij(xi,yj,zij),其在XOY平面上的投影点与网格中心点重合,其中
zij=f(xi,yj)。
以使得电磁波经过每段的时间长度为Δt,每段的长度Δl=(Te-Tb)/Ns;
所述步骤S4包括:
S401.遍历所有曲面上的网格中心映射点Rij,计算Rij的瞬态感应电流Jij(t);
由时域物理光学近似,得到规则曲面被脉冲平面电磁波照亮区域上上任意一点r′在某一个时刻t的表面感应电流JPO(r′,t)为:
以上步骤完成了规则曲面到XOY平面的平面投影矩形网格模型的映射,得到了矩形网格模型上每个矩形网格对应的规则曲面所受瞬态平面波激励导致的t=iCur*Δt时刻瞬态表面感应场Ju、,Jv、Jtx、Jty和Jtz。
S5.生成电磁散射中目标瞬态感应特征图像:
分为为列表Z、Jtx、Jty和Jtz创建四个分辨率都为nX×nY的数字图像数据ImagZ,ImagJtx,ImagJty,ImagJtz,并将列表Z、Jtx、Jty和Jtz的矩阵数据按照行列的二重循环的顺序分别逐一赋值给ImagZ,ImagJtx,ImagJty,ImagJtz的各个像素点,这四个数字图像采用分辨率为nX×nY,每个像素代表的实际物理尺寸分别是Δx和Δy;设置数字图像的颜色深度和显示的通道数为1;
分别将这四个数字图像数据分别保存至四个BMP文件,取名分别为Iz.bmp、Itx.bmp、Ity.bmp和Itz.bmp。
在本申请的实施例中,按照本申请的方法,进行高斯曲面目标表面感应瞬态电流特征图像表征进行处理,涉及到的参数如下
其中,τ=4/fBW,t0=0.8τ
得到的Iz、Itx、Ity、Itz如图6~9所示;从图8和图9中可以清楚且直观地观察到目标曲面在瞬态平面波映射下的散射中心即为图像中清晰显示的波动中心位置点,结合图6的高度特征图像可以看到,波动中心位置点即与图6中所示的曲面高度的峰值点基本重合,由此可见,曲面的峰值点是构成目标散射中心的一类非常重要的因素。图7可以看到没有散射中心的存在,说明散射中心具有矢量特性,不是在各个方向都存在,而是与入射和散射方向有关,例如在本例的X方向上就不存在散射中心。
以上所述是本发明的优选实施方式,应当理解本发明并非局限于本文所披露的形式,不应该看作是对其他实施例的排除,而可用于其他组合、修改和环境,并能够在本文所述构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
Claims (6)
1.一种电磁散射中目标瞬态感应特征图像生成方法,其特征在于:包括以下步骤:
S1.给定规则曲面受脉冲平面电磁波辐射的场景信息;
S2.将规则曲面在XOY平面进行投影,得到投影的矩形区域RXOY,并在矩形区域RXOY内进行网格生成,得到平面二维离散矩形网格;
S3.遍历平面二维离散矩形网格中的每一个网格,得到每一个网格的中心点坐标和在曲面上的映射点,计算平面波波前从O点到达曲面的最小和最大光程,并进行区间分割得到时间序列和幅度序列;
S4.对于二维离散矩形网格中的每一个网格,确定网格中心点在规则曲面上的映射点信息列表;
S5.生成电磁散射中目标瞬态感应特征图像。
4.根据权利要求1所述的一种电磁散射中目标瞬态感应特征图像生成方法,其特征在于:所述步骤S3包括:
S301:遍历X轴上的(nX-1)个网格和Y轴上的(nY-1)个网格,得到X轴上的第i个,Y轴上的第j个网格Gij,Gij网格中心点二维坐标Pij(xi,yj)为:
xi=(i+0.5)Δx,yj=(j+0.5)Δy
S302:网格中心点二维坐标(xi,yi)代入规则曲面表征公式得到该网格中心点在曲面上的映射点Rij(xi,yj,zij),其在XOY平面上的投影点与网格中心点重合,其中
zij=f(xi,yj)。
以使得电磁波经过每段的时间长度为Δt,每段的长度Δl=(Te-Tb)/Ns;
5.根据权利要求1所述的一种电磁散射中目标瞬态感应特征图像生成方法,其特征在于:所述步骤S4包括:
S401.遍历所有曲面上的网格中心映射点Rij,计算Rij的瞬态感应电流Jij(t);
由时域物理光学近似,得到规则曲面被脉冲平面电磁波照亮区域上上任意一点r′在某一个时刻t的表面感应电流JPO(r′,t)为:
6.根据权利要求1所述的一种电磁散射中目标瞬态感应特征图像生成方法,其特征在于:所述步骤S5包括:
分为为列表Z、Jtx、Jty和Jtz创建四个分辨率都为nX×nY的数字图像数据ImagZ,ImagJtx,ImagJty,ImagJtz,并将列表Z、Jtx、Jty和Jtz的矩阵数据按照行列的二重循环的顺序分别逐一赋值给ImagZ,ImagJtx,ImagJty,ImagJtz的各个像素点,这四个数字图像采用分辨率为nX×nY,每个像素代表的实际物理尺寸分别是Δx和Δy;设置数字图像的颜色深度和显示的通道数为1;
分别将这四个数字图像数据分别保存至四个BMP文件,取名分别为Iz.bmp、Itx.bmp、Ity.bmp和Itz.bmp。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211339858.1A CN115661293B (zh) | 2022-10-27 | 2022-10-27 | 一种电磁散射中目标瞬态感应特征图像生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211339858.1A CN115661293B (zh) | 2022-10-27 | 2022-10-27 | 一种电磁散射中目标瞬态感应特征图像生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115661293A true CN115661293A (zh) | 2023-01-31 |
CN115661293B CN115661293B (zh) | 2023-09-19 |
Family
ID=84993467
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211339858.1A Active CN115661293B (zh) | 2022-10-27 | 2022-10-27 | 一种电磁散射中目标瞬态感应特征图像生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115661293B (zh) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040156558A1 (en) * | 2003-02-04 | 2004-08-12 | Kim Sang Yeon | Image warping method and apparatus thereof |
JP2005149165A (ja) * | 2003-11-17 | 2005-06-09 | Sony Corp | 画像合成方法、プログラム、記録媒体、画像合成装置およびシステム |
US20100220895A1 (en) * | 2007-06-07 | 2010-09-02 | Zvi Koren | Device and method for displaying full azimuth angle domain image data |
US20130272600A1 (en) * | 2010-10-15 | 2013-10-17 | Iee International Electronics & Engineering S.A. | Range image pixel matching method |
CN105930567A (zh) * | 2016-04-14 | 2016-09-07 | 西安电子科技大学 | 一种基于子区域自适应积分的电磁散射特性获取方法 |
US20190128739A1 (en) * | 2017-11-02 | 2019-05-02 | Fluke Corporation | Focus and/or Parallax Adjustment in Acoustic Imaging Using Distance Information |
CN109981920A (zh) * | 2017-12-28 | 2019-07-05 | 夏普株式会社 | 图像处理装置、图像处理程序、图像处理方法及图像形成装置 |
CN111199097A (zh) * | 2019-12-25 | 2020-05-26 | 西安电子科技大学 | 瞬态矢量辐射传输理论箔条云散射处理方法 |
US20200366838A1 (en) * | 2017-08-03 | 2020-11-19 | Hangzhou Hikvision Digital Technology Co., Ltd. | Panoramic image generation method and device |
US20210166353A1 (en) * | 2019-11-29 | 2021-06-03 | Samsung Electronics Co., Ltd. | Generation of bokeh images using adaptive focus range and layered scattering |
US20220044474A1 (en) * | 2019-11-21 | 2022-02-10 | Dalian University Of Technology | Method for constructing grid map by using binocular stereo camera |
US11287390B1 (en) * | 2019-01-30 | 2022-03-29 | Reservoir Labs, Inc. | System and method for electromagnetic beamforming and imaging at low frequency |
CN114450716A (zh) * | 2019-08-29 | 2022-05-06 | 皇家飞利浦有限公司 | 用于卒中表征的图像处理 |
-
2022
- 2022-10-27 CN CN202211339858.1A patent/CN115661293B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040156558A1 (en) * | 2003-02-04 | 2004-08-12 | Kim Sang Yeon | Image warping method and apparatus thereof |
JP2005149165A (ja) * | 2003-11-17 | 2005-06-09 | Sony Corp | 画像合成方法、プログラム、記録媒体、画像合成装置およびシステム |
US20100220895A1 (en) * | 2007-06-07 | 2010-09-02 | Zvi Koren | Device and method for displaying full azimuth angle domain image data |
US20130272600A1 (en) * | 2010-10-15 | 2013-10-17 | Iee International Electronics & Engineering S.A. | Range image pixel matching method |
CN105930567A (zh) * | 2016-04-14 | 2016-09-07 | 西安电子科技大学 | 一种基于子区域自适应积分的电磁散射特性获取方法 |
US20200366838A1 (en) * | 2017-08-03 | 2020-11-19 | Hangzhou Hikvision Digital Technology Co., Ltd. | Panoramic image generation method and device |
CN109756720A (zh) * | 2017-11-02 | 2019-05-14 | 弗兰克公司 | 使用距离信息的声学成像中的焦点和/或视差调整 |
US20190128739A1 (en) * | 2017-11-02 | 2019-05-02 | Fluke Corporation | Focus and/or Parallax Adjustment in Acoustic Imaging Using Distance Information |
CN109981920A (zh) * | 2017-12-28 | 2019-07-05 | 夏普株式会社 | 图像处理装置、图像处理程序、图像处理方法及图像形成装置 |
US11287390B1 (en) * | 2019-01-30 | 2022-03-29 | Reservoir Labs, Inc. | System and method for electromagnetic beamforming and imaging at low frequency |
CN114450716A (zh) * | 2019-08-29 | 2022-05-06 | 皇家飞利浦有限公司 | 用于卒中表征的图像处理 |
US20220044474A1 (en) * | 2019-11-21 | 2022-02-10 | Dalian University Of Technology | Method for constructing grid map by using binocular stereo camera |
US20210166353A1 (en) * | 2019-11-29 | 2021-06-03 | Samsung Electronics Co., Ltd. | Generation of bokeh images using adaptive focus range and layered scattering |
CN111199097A (zh) * | 2019-12-25 | 2020-05-26 | 西安电子科技大学 | 瞬态矢量辐射传输理论箔条云散射处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115661293B (zh) | 2023-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112819962B (zh) | 数字图像相关中非均匀网格划分及局部网格疏密方法 | |
KR100305419B1 (ko) | 영상 위치맞춤 방법, 영상 위치맞춤 장치 및 기록 매체 | |
CN106548462B (zh) | 基于薄板样条插值的非线性sar图像几何校正方法 | |
AU2009323838A2 (en) | Provision of image data | |
CN112612024A (zh) | 微波阵列快速成像方法 | |
CN111161338A (zh) | 一种基于二维图像灰度进行深度预测的点云密度提升方法 | |
CN113390514B (zh) | 一种基于多传感器阵列的三维红外温度测量方法 | |
CN113177356A (zh) | 一种基于深度学习的目标电磁散射特性快速预测方法 | |
CN114996658A (zh) | 一种基于投影的高超声速飞行器气动热预测方法 | |
Tian et al. | Performance evaluation of deception against synthetic aperture radar based on multifeature fusion | |
CN106054190B (zh) | 基于频谱优化建模的双基地前视sar频域成像方法 | |
CN105139432B (zh) | 基于高斯模型的红外弱小目标图像仿真方法 | |
CN115661293A (zh) | 一种电磁散射中目标瞬态感应特征图像生成方法 | |
CN115656219B (zh) | 一种规则曲面表面电流感应特征图像生成方法 | |
CN115908624B (zh) | 一种天线布局感应场特征图像生成方法 | |
CN115698690A (zh) | 图像处理方法、形状检查方法、图像处理系统以及形状检查系统 | |
Flores et al. | Choice of an entropy-like function for range-Doppler processing | |
CN103630895B (zh) | 一种用于毫米波近距离三维成像系统的成像方法 | |
CN113687099B (zh) | 一种基于激光线阵的piv图像标定装置及方法 | |
CN115205354A (zh) | 基于ransac和icp点云配准的相控阵激光雷达成像方法 | |
CN111553876B (zh) | 一种气动光学视线误差图像处理方法及系统 | |
CN114280602A (zh) | 用于毫米波综合孔径辐射计的近场三维成像系统及方法 | |
CN114252876A (zh) | 基于多散射点模型的高分辨一维成像方法 | |
CN110310312B (zh) | 一种图像配准方法、图像配准装置及终端 | |
CN106780581B (zh) | 一种sar图像解译能力的评估方法 |
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 |