CN106093902B - 相似目标的散射特征变换方法 - Google Patents

相似目标的散射特征变换方法 Download PDF

Info

Publication number
CN106093902B
CN106093902B CN201610436570.4A CN201610436570A CN106093902B CN 106093902 B CN106093902 B CN 106093902B CN 201610436570 A CN201610436570 A CN 201610436570A CN 106093902 B CN106093902 B CN 106093902B
Authority
CN
China
Prior art keywords
target
scattering
signatures
scattering center
center
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.)
Expired - Fee Related
Application number
CN201610436570.4A
Other languages
English (en)
Other versions
CN106093902A (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical 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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201610436570.4A priority Critical patent/CN106093902B/zh
Publication of CN106093902A publication Critical patent/CN106093902A/zh
Application granted granted Critical
Publication of CN106093902B publication Critical patent/CN106093902B/zh
Expired - Fee Related 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
    • 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/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

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

Abstract

本发明涉及一种相似目标的散射特征变换方法,对提取的散射中心建立一一映射,将相似目标间的联系简化成散射中心的联系,实现了由一个目标的散射特征变换到另一个目标的散射特征的过程。本发明的目的是为了解决在航天航空领域中一些典型目标的散射特性计算和存储量大的问题,通过本发明提供的方法,可以快速进行散射特性计算,并且存储的数据量远小于原数据量。

Description

相似目标的散射特征变换方法
技术领域
本发明属于雷达目标电磁散射领域,具体涉及一种相似目标的散射特征变换方法。
背景技术
目标的散射特性是必须考察的重要指标之一。当目标外形或结构发生微小变化后,其散射特征一般都会随之有所变化。对所有目标在每个研制阶段都进行仿真计算、静态测试,显然是非常耗时、花费很高的。另一方面,同一目标两个研制阶段的外形或结构的小改变所带来的散射特性的变化可能会有一定规律,或者两个散射特性直接可能存在一定的变换关系。而当前并没有这方面的技术,如果能够通过深入研究,理解和掌握这种变换关系,必将会为工程实践带来极大的好处。
发明内容
要解决的技术问题
为了解决在航天航空领域中一些典型目标的散射特性计算和存储量大的问题,本发明提出一种相似目标的散射特征变换方法。
技术方案
一种相似目标的散射特征变换方法,其特征在于步骤如下:
步骤1:通过对相似目标a和目标b进行雷达成像和散射中心提取,获取目标a和目标b的强散射点数据,包括位置、幅度和相位信息;
步骤2:根据目标a、目标b的各个散射中心幅度之间的变换关系,得到幅度映射参数(aAn,bθn):
其中,A1n、A2n分别为目标a、目标b第n个散射中心的复幅度,θ1n、θ2n分别为目标a、目标b第n个散射中心的复幅度的相位;
根据目标a、目标b各个散射中心位置之间的平移变换关系,得到位置映射参数(bxn,byn);
其中,(x1n,y1n)为目标a第n个散射中心的位置,(x2n,y2n)为目标b第n个散射中心的位置;
存储目标b的散射场数据和幅度映射参数(aAn,bθn)、位置映射参数(bxn,byn);
步骤3:将目标b的散射场数据和幅度映射参数(aAn,bθn)、位置映射参数(bxn,byn)代入散射特征变换关系式反演到目标a的散射特征:
其中:是目标a的总散射场,是目标b的第n个散射中心的散射场,是传播方向与x坐标轴的夹角,N为散射中心的总数量。
有益效果
本发明提出的一种相似目标的散射特征变换方法,对提取的散射中心建立一一映射,将相似目标间的联系简化成散射中心的联系,实现了由一个目标的散射特征变换到另一个目标的散射特征的过程。通过本发明提供的方法,可以快速进行散射特性计算,并且存储的数据量远小于原数据量。
附图说明
图1平移变换
图2极坐标格式下位置变换
图3一组点目标模型:(a)相似目标a;(b)相似目标b
图4目标散射特征对比:(a)中心频率下散射特征;(b)中心角度下散射特征
图5本发明方法流程图
具体实施方式
现结合实施例、附图对本发明作进一步描述:
散射特征变换关系式的推导:
根据散射中心理论,在高频区,目标总的电磁散射可以认为是由某些局部位置上电磁散射的相干合成,故目标的后向散射场可以表示为:
其中:An是第n个散射中心的复幅度,xn,yn是第n个散射中心的位置,是第n个散射中心位置的位移矢量,k是波数,kx、ky分别为k在x、y轴上的分量,相位φ∈(-π,π)。
根据散射中心模型,建立两个目标各散射中心之间的变换关系:
具体包括以下步骤,
1):进行幅度变换,将幅度表示成如下形式:
故可以得到两个散射中心幅度之间的变换关系:
需要确定的映射参数为:(aAn,bθn),该映射参数可通过幅度变换关系及散射场数据得到。
2):进行位置变换,二维雷达像的散射中心位置即为成像平面上一个点的二维坐标,二维平面上的位置变换总能通过下面两种变换中的任一种来实现,一种是平移变换(图1);另一种变换(图2)可用极坐标形式实现,其中只有缩放和旋转变换,没有平移变换,径向坐标通过缩放变换得到,周向坐标通过旋转变换得到。显然这种形式没有平移变换简便,故下面将采用平移变换的形式,将各散射点之间的位置关系表示为
由此,我们可以得到散射特征变换关系式:
其中,n表示强散射点的个数,k表示波数,A1n和A2n分别表示目标a和目标b第n个散射中心的复幅度,是目标a和目标b第n个散射中心位置的位移矢量,(x1n,y1n)、(x2n,y2n)是直角坐标系上相应的位置,是传播方向与x坐标轴的夹角。
利用得到散射特征变换关系式由目标b反演目标a:
步骤1:通过对相似目标a和目标b进行雷达成像和散射中心提取,获取目标a和目标b的强散射点数据,包括位置、幅度和相位信息;
步骤2:根据目标a、目标b的各个散射中心幅度之间的变换关系,得到幅度映射参数(aAn,bθn):
其中,A1n、A2n分别为目标a、目标b第n个散射中心的复幅度,θ1n、θ2n分别为目标a、目标b第n个散射中心的复幅度的相位;
根据目标a、目标b各个散射中心位置之间的平移变换关系,得到位置映射参数(bxn,byn);
其中,(x1n,y1n)为目标a第n个散射中心的位置,(x2n,y2n)为目标b第n个散射中心的位置;
存储目标b的散射场数据和幅度映射参数(aAn,bθn)、位置映射参数(bxn,byn);
步骤3:将目标b的散射场数据和幅度映射参数(aAn,bθn)、位置映射参数(bxn,byn)代入散射特征变换关系式反演到目标a的散射特征:
通过该式对目标a进行目标散射特性反演,从而实现了由一个目标的散射特征变换到另一个目标的散射特征。
仿真过程及仿真验证:
一组点目标模型如图3所示,仿真条件为:起始频率为35GHz,带宽为5GHz,频率采样点数为150,观测方位角为:-3.819719°~3.819719°,角度采样点数为150
首先,分别对目标a和目标b进行雷达成像和散射中心提取,提取5个散射中心的散射场数据。根据提取的散射中心的数据确定幅度映射参数和位置映射参数并存储,然后将目标b的散射中心数据及存储的幅度映射参数和位置映射参数代入到散射特征变换关系式,依据此变换关系式将目标b反演目标a,得到结果如图4所示。
从图4可以看出,由所建立的两个目标之间的变换方法得到的散射特征与目标的原始散射特征非常吻合,证明了该方法的正确性。根据散射中心理论,目标总的电磁散射可以认为是由某些局部位置上电磁散射的相干合成,主要贡献来自于强散射源,该散射特征变换关系式是强散射源的叠加,故可以体现目标的散射特征。另外,从散射特征变换关系式可以看出,与常规的电磁计算相比,该方法的计算量少,可以快速得到目标的散射特征,并且从仿真中可以看出,该方法仅需存储目标b的5个散射中心的数据和幅度映射参数及位置映射参数,而原来则需要存储全部的散射场数据,这大大降低了存储量。

Claims (1)

1.一种相似目标的散射特征变换方法,其特征在于步骤如下:
步骤1:通过对相似目标a和目标b进行雷达成像和散射中心提取,获取目标a和目标b的强散射点数据,包括位置、幅度和相位信息;
步骤2:根据目标a、目标b的各个散射中心幅度之间的变换关系,得到幅度映射参数(aAn,bθn):
其中,A1n、A2n分别为目标a、目标b第n个散射中心的复幅度,θ1n、θ2n分别为目标a、目标b第n个散射中心的复幅度的相位;
根据目标a、目标b各个散射中心位置之间的平移变换关系,得到位置映射参数(bxn,byn);
其中,(x1n,y1n)为目标a第n个散射中心的位置,(x2n,y2n)为目标b第n个散射中心的位置;
存储目标b的散射场数据和幅度映射参数(aAn,bθn)、位置映射参数(bxn,byn);
步骤3:将目标b的散射场数据和幅度映射参数(aAn,bθn)、位置映射参数(bxn,byn)代入散射特征变换关系式反演到目标a的散射特征:
其中:是目标a的总散射场,是目标b的第n个散射中心的散射场,是传播方向与x坐标轴的夹角,N为散射中心的总数量。
CN201610436570.4A 2016-06-17 2016-06-17 相似目标的散射特征变换方法 Expired - Fee Related CN106093902B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610436570.4A CN106093902B (zh) 2016-06-17 2016-06-17 相似目标的散射特征变换方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610436570.4A CN106093902B (zh) 2016-06-17 2016-06-17 相似目标的散射特征变换方法

Publications (2)

Publication Number Publication Date
CN106093902A CN106093902A (zh) 2016-11-09
CN106093902B true CN106093902B (zh) 2018-06-12

Family

ID=57236191

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610436570.4A Expired - Fee Related CN106093902B (zh) 2016-06-17 2016-06-17 相似目标的散射特征变换方法

Country Status (1)

Country Link
CN (1) CN106093902B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106991215B (zh) * 2017-03-15 2020-04-03 西北工业大学 一种基于gpu的高精度图形电磁散射计算方法
CN107864020B (zh) * 2017-11-07 2021-02-19 哈尔滨工程大学 水下小目标单分量声散射回波的变换域提取方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103809175B (zh) * 2014-02-26 2017-01-18 北京航空航天大学 隐身飞机散射特性现场快速评估系统
CN104515908B (zh) * 2014-12-10 2017-06-16 北京环境特性研究所 一种超宽带目标电磁散射特性测试系统和方法
CN105467370A (zh) * 2015-12-01 2016-04-06 中国人民解放军国防科学技术大学 一种复合双基地雷达进动目标isar图像横向定标方法
CN105572652B (zh) * 2016-03-18 2018-01-05 西北工业大学 一种利用外推获得具有多次散射目标远场rcs的方法

Also Published As

Publication number Publication date
CN106093902A (zh) 2016-11-09

Similar Documents

Publication Publication Date Title
CN106548462B (zh) 基于薄板样条插值的非线性sar图像几何校正方法
CN104730520B (zh) 基于子孔径合成的圆周sar后向投影自聚焦方法
CN103454638B (zh) 一种圆迹合成孔径雷达三维层析成像方法
CN102393183A (zh) 基于控制网的海量点云快速配准方法
CN102944872B (zh) 雷达散射截面近场到近场的变换方法
CN107102324A (zh) 一种近景微波成像方法及系统
CN106093902B (zh) 相似目标的散射特征变换方法
CN102411794A (zh) 一种基于球谐变换的三维模型的二维投影的输出方法
CN103530627B (zh) 基于二维散射中心集网格模型的isar图像恢复方法
Wang et al. An improved measurement model of binocular vision using geometrical approximation
CN109541589A (zh) 空间自旋目标雷达的三维成像方法、装置及存储介质
Ji et al. The mutual inductance calculation between circular and quadrilateral coils at arbitrary attitudes using a rotation matrix for airborne transient electromagnetic systems
CN103631990A (zh) Sar照射区域的仿真场景模型建立方法和系统
CN104062656B (zh) 基于矢量面元模型的微波成像方法
Wang et al. Interference image registration combined by enhanced scale-invariant feature transform characteristics and correlation coefficient
CN107610161B (zh) 一种基于四叉树分割的InSAR快速图像配准方法
CN108629796B (zh) 一种基于波程差序列的InSAR图像配准方法
Borges et al. Inverse scattering reconstruction of a three dimensional sound-soft axis-symmetric impenetrable object
Sun et al. Extended object tracking based on support functions and extended gaussian images
CN104167019B (zh) 一种sar图像精确度计算方法
Feng et al. Multiview ISAR imaging for complex targets based on improved SBR scattering model
CN106291549B (zh) 一种金属多边形的合成孔径雷达图像仿真方法及其应用
CN110515080A (zh) 基于半径角度插值的高超声速机动平台sar成像方法
Luo et al. An Effencient ESM Technique for Three-Dimension Radiation Source Localization
CN104750920B (zh) 一种对地物成像的仿真系统和方法

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

Granted publication date: 20180612

Termination date: 20190617