CN104049243A - 一种合成孔径雷达图像仿真的方法及系统 - Google Patents

一种合成孔径雷达图像仿真的方法及系统 Download PDF

Info

Publication number
CN104049243A
CN104049243A CN201410261722.2A CN201410261722A CN104049243A CN 104049243 A CN104049243 A CN 104049243A CN 201410261722 A CN201410261722 A CN 201410261722A CN 104049243 A CN104049243 A CN 104049243A
Authority
CN
China
Prior art keywords
image
noise energy
system sensitivity
measuring image
under
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
CN201410261722.2A
Other languages
English (en)
Other versions
CN104049243B (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.)
Institute of Electronics of CAS
Original Assignee
Institute of Electronics of CAS
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 Institute of Electronics of CAS filed Critical Institute of Electronics of CAS
Priority to CN201410261722.2A priority Critical patent/CN104049243B/zh
Publication of CN104049243A publication Critical patent/CN104049243A/zh
Application granted granted Critical
Publication of CN104049243B publication Critical patent/CN104049243B/zh
Active 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/40Means for monitoring or calibrating
    • G01S7/4052Means for monitoring or calibrating by simulation of echoes
    • 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
    • G01S13/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques

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)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种合成孔径雷达图像仿真的方法及系统,所述方法包括:根据实测图像获取指定系统灵敏度下的图像噪声能量差;当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像;或者,根据仿真图像获取指定系统灵敏度下的图像噪声能量差;当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像。

Description

一种合成孔径雷达图像仿真的方法及系统
技术领域
本发明涉及合成孔径雷达(SAR,Synthetic Aperture Radar)成像技术领域,尤其涉及一种SAR图像仿真的方法及系统。
背景技术
本申请发明人在实现本申请实施例技术方案的过程中,至少发现相关技术中存在如下技术问题:
星载极化SAR系统的设计包括频率选择、轨道参数、几何参数、雷达参数等技术,而雷达参数是星载SAR系统设计中的关键参数之一。在雷达参数中系统灵敏度(NEσ0)又称为噪声等效后向散射系数,是衡量雷达对弱目标成像能力的指标。NEσ0与SAR成像质量密切相关。不同NEσ0下的雷达图像效果及NEσ0的设计是否满足应用需求是SAR及极化SAR系统设计工作者十分关心的。
在现有的SAR图像仿真技术中,往往给定一个指定的系统灵敏度,在系统仿真时,所输出的SAR图像的仿真结果的实际系统灵敏度与设计人员给定的系统灵敏度具有一定的出入,从而造成SAR的成像能力受一定的影响,使得设计人员得到SAR图像仿真结果常常不是所想要的指定系统灵敏度下的仿真图像。
发明内容
为解决现有存在的上述问题,本发明实施例期望提供一种SAR图像仿真的方法及系统,使得SAR图像仿真结果为指定系统灵敏度下的仿真图像。
为达到上述目的,本发明实施例的技术方案是这样实现的:
本发明实施例提供一种合成孔径雷达图像仿真的方法,所述方法包括:根据实测图像获取指定系统灵敏度下的图像噪声能量差;当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像;或者,根据仿真图像获取指定系统灵敏度下的图像噪声能量差;当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像。
上述方案中,所述根据实测图像获取指定系统灵敏度下的图像噪声能量差包括:选择与待仿真合成孔径雷达系统的系统参数一致的实测图像;根据所述系统参数确定所述实测图像的系统灵敏度;根据所述实测图像统计所述实测图像的噪声能量;根据所述实测图像的系统灵敏度与所述实测图像的噪声能量确定所述图像噪声能量差。
上述方案中,所述根据所述实测图像的系统灵敏度与所述实测图像的噪声能量确定所述图像噪声能量差包括:根据所述实测图像的系统灵敏度、所述实测图像的噪声能量得到目标功率与目标后向散射系数之比;通过所述指定系统灵敏度、所述目标功率与目标后向散射系数之比确定得到指定系统灵敏度下的噪声能量;根据所述指定系统灵敏度下噪声能量、所述实测图像的噪声能量确定所述图像噪声能量差。
上述方案中,所述根据待仿真图像获取指定系统灵敏度下的图像噪声能量包括:获取所述仿真图像的目标后向散射系数和目标功率;根据所述目标后向散射系数、所述目标功率、所述指定系统灵敏度得到指定系统灵敏度下的图像噪声能量差。
上述方案中,所述利用所述图像噪声能量差得到系统灵敏度下的输出图像包括:将所述图像噪声能量差以正太分布噪声的形式加入到实测或仿真图像单视复数据的实部、虚部,得到指定系统灵敏度下的输出图像。
本发明实施例还提供一种合成孔径雷达图像仿真的系统,所述系统包括:噪声能量获取单元、图像输出单元;其中,所述噪声能量获取单元,用于根据实测图像获取指定系统灵敏度下的图像噪声能量差;或者,根据仿真图像获取指定系统灵敏度下的图像噪声能量差;所述图像输出单元,用于当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像。
上述方案中,所述噪声能量获取单元,包括:选择子单元、确定子单元、统计子单元、噪声能量差确定子单元;其中,所述选择子单元,用于选择与待仿真合成孔径雷达系统的系统参数一致的实测图像;所述确定子单元,用于根据所述系统参数确定所述实测图像的系统灵敏度;所述统计子单元,用于根据所述实测图像统计所述实测图像的噪声能量;所述噪声能量差确定子单元,用于根据所述实测图像的系统灵敏度与所述实测图像的噪声能量确定需加入到实测图像的噪声能量差。
上述方案中,所述噪声能量差确定单元,具体用于根据所述实测图像的系统灵敏度、所述实测图像的噪声能量得到目标功率与目标后向散射系数之比;通过所述指定系统灵敏度、所述目标功率与目标后向散射系数之比确定得到指定系统灵敏度下的噪声能量;根据所述指定系统灵敏度下的噪声能量、所述实测图像的噪声能量确定待加入到实测图像的噪声能量差。
上述方案中,所述噪声能量获取单元包括:获取子单元、噪声能量确定子单元;其中,所述获取子单元,用于获取所述仿真图像的目标后向散射系数和目标功率;所述噪声能量确定子单元,用于根据所述目标后向散射系数、所述目标功率、所述指定系统灵敏度得到指定系统灵敏度下的图像噪声能量差。上述方案中,所述图像输出单元,具体用于当输出图像时,将所述图像噪声能量差以正太分布噪声的形式加入到实测或仿真单视复图像数据的实部、虚部,得到指定系统灵敏度下的输出图像。
本发明实施例提供的一种SAR图像仿真的方法,根据实测图像获取指定系统灵敏度下的图像噪声能量差,当输出图像时,利用所述图像噪声能量差得到系统灵敏度的输出图像;或者,根据仿真图像获取指定系统灵敏度下的图像噪声能量差;当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像;由此,使得SAR图像仿真结果为指定系统灵敏度下的仿真图像,解决了SAR图像仿真结果的实际系统灵敏度与要求的系统灵敏度不一致的问题。
附图说明
图1为本发明实施例提供的一种SAR图像仿真的方法的流程示意图;
图2为本发明实施例提供的另一种SAR图像仿真的方法的流程示意图;
图3为本发明实施例提供的基于实测SAR图像仿真指定系统灵敏度对应的SAR图像的流程示意图;
图4为本发明实施例提供的基于仿真图像仿真指定系统灵敏度对应的SAR图像的流程示意图;
图5为一具体的实测图像;
图6为基于图5所示的实测图像仿真的指定系统灵敏度为-30dB下的输出图像;
图7为本发明实施例提供的一种SAR图像仿真的系统。
具体实施方式
在本发明实施例中,根据实测图像获取指定系统灵敏度下的图像噪声能量差,当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像;或者,根据仿真图像获取指定系统灵敏度下的图像噪声能量差,当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像。
在本发明实施例中,得到指定系统灵敏度下的图像噪声能量差分为两种情况:
情况1,基于实测图像输出待仿真系统的图像;
情况2,基于仿真图像输出待仿真系统的图像。
在对本发明实施例的方案进行进一步说明之前,对系统灵敏度系数进行说明。
在雷达参数中,系统灵敏度为噪声等效后向散射系数,与SAR成像质量密切相关,是衡量雷达对弱目标成像能力的指标。
系统灵敏度的定义如下: NEσ 0 = ( 4 π ) 3 R 3 ( kT 0 ) · Loss · 2 v s G 2 · λ 3 · P ‾ · ρ gr - - - ( 1 )
其中,公式(1)中各参数定义选下:
NEσ0,系统灵敏度;
发射机平均功率;
发射和接收天线功率增益;Gτ,发射天线功率增益;接收天线功率增益;通常,收发天线是共用的。
λ,雷达波长;
ρgr,地距分辨率;
R,目标到雷达的距离;
vs,SAR飞行速度;
Loss,系统损耗因子(>1);
k,波尔兹曼常数,等于1.38054×10-23J/K;
T0,常温(290K)。
系统灵敏度与图像信噪比和目标后向散射系数关系如下:
S N = σ 0 NEσ 0 - - - ( 2 )
其中,S为目标功率,N为目标噪声功率,为图像信噪比,σ0为目标后向散射系数。
NEσ0的物理意义为雷达输出信噪比为0dB(,即)时对应的后向散射系数。
下面通过附图及具体实施例对本发明实施例再做进一步的详细说明。
情况1,基于实测图像输出待仿真系统的图像;如图1所示,包括以下步骤:
步骤101,根据实测图像获取图指定系统灵敏度下图像噪声能量差;
这里,指定系统灵敏度为NEσ0,图像噪声能量为NA,根据实测图像获取图像噪声能量误差NA包括以下步骤:
步骤101a,选择与待仿真SAR系统的系统参数一致的实测图像;
这里,系统参数包括发射机平均功率、发射和接收天线功率增益、雷达波长、地距分辨率、目标到雷达的距离、SAR飞行速度、系统损耗因子等系统参数。这里,通过选择的实测图像的系统参数与待仿真SAR系统的系统参数一致的实测图像,可获得待仿真系统的图像噪声能量差NA
步骤101b,根据所述系统参数确定所述实测图像的系统灵敏度;
在选择好实测图像之后,根据所选择的实测图像的系统参数计算实测图像的系统灵敏度其中,计算实测图像的系统灵敏度采用的公式为公式(1)。这里,将实测图像的系统参数通过公式(1)的计算,得到的系统灵敏度为实测图像的系统灵敏度
步骤101c,根据所述实测图像统计所述实测图像的噪声能量差;
具体的,利用实测图像统计实测图像的噪声能量NR;在统计实测图像的噪声能量时,选择实测图像中的弱目标作为统计的对象。
在实际应用中,统计阴影区域所有像素的方差,将所统计的方差作为实测图像的噪声能量;阴影区域可选择背对雷达电磁波照射方向的山体或高建筑物阴影区,且实测图像的阴影区域的像素个数至少为100个,以满足大样本数的统计需求。
在实际应用中,步骤101b与步骤101c的执行不分先后顺序,可先执行步骤101b,再执行步骤101c,也可先执行步骤101c,再执行步骤101b,同时,步骤101b与步骤101c也可同时执行。
步骤101d,根据所述实测图像的系统灵敏度与所述实测图像的噪声能量NR确定所述图像噪声能量差,即:待加入到实测图像中的噪声能量;
这里,通过步骤101b与步骤101c的执行,获得实测图像的系统灵敏度与噪声能量NR后,进一步确定图像噪声量差NA;其中,图像噪声量差NA的确定包括以下步骤:
步骤101d1:根据所述实测图像的系统灵敏度、所述实测图像的噪声能量得到目标功率与目标后向散射系数之比;
具体的,根据公式(2)可确定实测图像的系统灵敏度与噪声能量NR存在关系:其中,KR为目标功率与目标后向散射系数之比。
步骤101d2,通过所述指定系统灵敏度、所述目标功率与目标后向散射系数之比确定得到指定系统灵敏度下的噪声能量;
具体的,由公式(2)得得到指定灵敏度系数NEσ0下的噪声能量N。
步骤101d3:根据所述指定系统灵敏度下噪声能量、所述实测图像的噪声能量确定所述图像噪声能量差;
具体的,根据公式NA=N-NR计算所述图像噪声能量差NA;其中,N为指定系统灵敏度NEσ0下的噪声能量,NR为实测图像的噪声能量。
步骤102,当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像;
具体的,将所述图像噪声能量差NA以正太分布噪声的形式加入到待仿真SAR系统的实测图像单视复数据的实部、虚部,即分别在待仿真SAR系统的实测图像单视复数据的实部、虚部加入均值为0方差为NA/2的正态分布噪声,从而得到指定系统灵敏度NEσ0下的输出图像;其中,待仿真SAR系统的实测图像为目标图像,在目标单视复图像数据的实部、虚部分别加入图像噪声能量差包括以下步骤:
步骤121,产生与目标图像像素个数相等的噪声图像1,也就是噪声图像实部,噪声图像1满足方差为NA/2均值为0的正态分布。
步骤122,将噪声图像1与目标图像的实部按像素值相加,得到加入噪声的仿真图像实部。
步骤123,产生与目标图像像素个数相等的噪声图像2,也就是噪声图像虚部,噪声图像2同样满足方差为NA/2均值为0的正态分布。
步骤124,将噪声图像2与目标图像的虚部按像素值相加,得到加入噪声的仿真图像虚部。
情况2,基于仿真图像输出待仿真系统的图像;如图2所示,包括以下步骤:
步骤201,根据仿真图像获取指定系统灵敏度下的图像噪声能量差;
具体的,获取所述仿真图像的目标后向散射系数和目标功率;根据所述目标后向散射系数、所述目标功率、所述指定系统灵敏度得到指定系统灵敏度下的图像噪声能量差。
根据仿真系统的设置,获取仿真系统本身已输出的仿真图像的目标后向散射系数σ0 s和目标功率Ss;在获取仿真图像的目标后向散射系数σ0 s和目标功率Ss之后,将σ0 s和目标功率Ss通过公式(2)得到指定系统灵敏度NEσ0下的图像噪声能量N:
N s = NEσ 0 S s σ 0 s
对于仿真图像而言,仿真图像中没有系统噪声,因此,仿真图像的系统灵敏度NEσ0=-∞,因此,图像噪声能量N即为指定系统灵敏度NEσ0下的待加入仿真图像的噪声能量差NA,即,NA=N。
在实际应用中,仿真图像中的每区域目标的σ0值及其对应信号强度可知,因此,由(2)式可直接求得NEσ0对应的图像噪声目标功率N。在仿真SAR图像中加入能量为N的噪声,得到的SAR图像对应的系统灵敏度即为NEσ0
步骤202,当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像;
具体的,将所述图像噪声能量差NA以正太分布噪声的形式加入到实测图像单视复数据的实部、虚部,即分别实测图像单视复数据的实部、虚部加入均值为0方差为N/2的正态分布噪声,从而得到指定系统灵敏度NEσ0下的输出图像,具体的过程同步骤102。
在SAR图像仿真的方法中,情况二与情况一的区别主要在于情况一采用实测图像计算图像噪声能量差,而情况二采用仿真图像计算图像噪声能量,而仿真图像与实测图像的不同点在于仿真图像中没有系统噪声,仿真图像实际的系统灵敏度为 NEσ R 0 = - ∞ .
通过本发明实施例的SAR图像的仿真方法,特别是对于弱目标的图像仿真结果来说,能够根据实际的要求,获取到不同系统灵敏度下与实际的图像更接近的SAR仿真图像,并且本发明实施例的方法对于SAR系统灵敏度指标的设计具有一定指导的作用,对SAR图像仿真和系统设计具有重要的实用价值。
图3为基于实测SAR图像仿真指定系统灵敏度对应的SAR图像的流程示意图,如图所示,包括以下步骤:
步骤301,根据系统参数计算实测图像的
这里,系统参数为已知的,选择与待仿真SAR系统的系统参数一致的实测图像,根据公式(1)结合系统参数计算实测图像的系统灵敏度
步骤302,统计实测图像的噪声能量NR
这里,选择实测图像的弱目标作为阴影区域,统计阴影区域的像素的方差,得到实测图像的噪声能量NR
步骤303,根据计算KR
这里,KR为目标功率与目标后向散射系数之比;根据公式由实测图像的系统灵敏度与噪声能量得到KR
步骤304,根据N=KR*NEσ0计算NEσ0下的噪声能量N;
计算指定系统灵敏度下的图像噪声能量;
这里,根据当前仿真需求给定一指定系统灵敏度NEσ0,在给定该指定系统灵敏度后,结合目标功率与目标后向散射系数之比KR由公式N=KR*NEσ0得出该指定系统灵敏度NEσ0下的噪声能量N。
步骤305,计算图像噪声能量差;
具体的,根据公式NA=N-NR计算图像噪声能量差NA,计算出的图像噪声能量差为需要在仿真系统中加入的噪声能量。
步骤306,在目标图像中加入噪声能量;
这里,目标图像为当前待仿真SAR系统的实测图像,根据步骤304中的产生图像噪声能量差生成方差为NA/2均值为0的正态分布噪声1与噪声分布噪声2,从而产生与目标图像像素个数相等的噪声图像1和噪声图像2,其中,噪声图像1满足方差为NA/2均值为0的正态分布,且其图像像素个数与目标图像实部像素个数相等,噪声图像2满足方差为NA/2均值为0的正态分布,且其图像像素个数与目标图像虚部像素个数相等;将产生的噪声图像1和噪声图像2分别加入目标的单视复图像数据中,得到目标仿真图像的实部与虚部,并将得到的目标仿真图像输出,得到输出图像,也就是指定系统灵敏度下的仿真图像。
图4为基于仿真SAR图像仿真指定系统灵敏度对应的SAR图像的流程示意图,如图所示,包括以下步骤:
步骤401,根据计算NEσ0下的噪声能量N。
这里,首先选择一仿真图像,该仿真图像为系统仿真的无噪仿真图像,获取该仿真图像的目标后向散射系数σ0 s和目标功率Ss。选定仿真图像之后,根据当前仿真需求给定一指定系统灵敏度NEσ0,在给定该指定系统灵敏度后,结合目标后向散射系数σ0 S和目标功率SS由公式得出该指定系统灵敏度NEσ0下的噪声能量N,由于仿真图像本身为无噪图像,噪声能量N即为图像噪声能量差。
步骤402,在目标图像中加入噪声能量;
这里,同步骤305过程一致,产生均值为0方差为N/2的正态分布噪声产生噪声图像,分别加入目标仿真图像的实部与虚部,并输出得到输出图像。
下面以指定系统灵敏度为-30dB下的情况一为例对本发明实施例进行进一步说明。这里,采用的实测图像如图5所示,具体过程如下:
步骤501,计算实测图像的系统灵敏度与噪声能量;
这里,该实测图像的SAR系统参数具体为:
G2=3.55×104;λ=0.33m;ρgr=0.5m;R=11km;vs=150m/S;Loss=2.5;k=1.38054×10-23J/K;T0=290K;根据公式(1)结合上述系统参数得到实测图像的系统灵敏度为-40dB。
选取图5中右上角的山体阴影区作为统计噪声能量的阴影区域,统计该区域的方差得到噪声能量NR=0.68。
步骤502,根据实测图像的灵敏度与噪声能量确定目标功率与目标后向散射系数之比;
具体的,根据 K R = N R NEσ R 0 = S σ 0 有:
K R = N R NEσ R 0 = 0.68 0.0001 .
步骤503,计算指定系统灵敏度下的图像噪声能量;
这里,在指定系统灵敏度为-30dB下,得到-30dB下的图像噪声能量 N = K R * NEσ 0 = 0.68 0.0001 * 0.001 = 6.8 .
步骤504,计算图像噪声能量差;
这里,求得需要加入到实测图像噪声能量,也就是图像噪声能量差NA为NA=N-NR=6.12。
步骤505,在目标图像中加入噪声能量;
具体的,在实测图像单视复图像数据的实部、虚部分别加入均值为0方差为3.06的正态分布噪声,则得到NEσ0为-30dB的输出图像,所得输出图像如图4所示。
为实现上述方法,本发明实施例还提供一种SAR图像仿真的系统,如图7所示,所述系统包括:噪声能量获取单元701、图像输出单元702;其中,
噪声能量获取单元701,用于根据实测图像获取指定系统灵敏度下的图像噪声能量差;或者,根据仿真图像获取指定系统灵敏度下的图像噪声能量差;
图像输出单元702,用于当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像。
图像输出单元702,具体用于当输出图像时,将所述图像噪声能量差以正太分布噪声的形式加入到实测图像单视复数据的实部、虚部,得到指定系统灵敏度下的输出图像。
噪声能量获取单元701包括:选择子单元711、确定子单元712、统计子单元713、噪声能量差确定子单元714;其中,
选择子单元711,用于选择与待仿真合成孔径雷达系统的系统参数一致的实测图像;
确定子单元712,用于根据所述系统参数确定所述实测图像的系统灵敏度;
统计子单元713,用于根据所述实测图像统计所述实测图像的噪声能量;
噪声能量差确定子单元714,用于根据所述实测图像的系统灵敏度与所述实测图像的噪声能量确定所述图像噪声能量差;
噪声能量差确定子单元714,具体用于根据所述实测图像的系统灵敏度、所述实测图像的噪声能量得到目标功率与目标后向散射系数之比;通过所述指定系统灵敏度、所述目标功率与目标后向散射系数之比确定得到指定系统灵敏度下的噪声能量;根据所述指定系统灵敏度下的噪声能量、所述实测图像的噪声能量确定所述图像噪声能量差。
如图7所示,噪声能量获取单元701包括:获取子单元715、噪声能量确定子单元716;其中,
获取子单元715,用于获取所述仿真图像的目标后向散射系数和目标功率;
噪声能量确定子单元716,用于根据所述目标后向散射系数、所述目标功率、所述指定系统灵敏度得到指定系统灵敏度下的图像噪声能量差。
在实际应用中,噪声能量获取单元701可只包括选择子单元711、确定子单元712、统计子单元713、噪声能量差确定子单元714;或只包括获取子单元715以及噪声能量确定子单元716,也可同时包括选择子单元711、确定子单元712、统计子单元713、噪声能量差确定子单元714、获取子单元715以及噪声能量确定子单元716。
在实际应用中,本发明提供的装置或系统装置可作为一个单独的系统,还可以是在现有的终端如个人计算机(PC,Personal Computer)中增加完成不同功能的逻辑单元。
当在PC中增加逻辑单元时,噪声能量获取单元701、图像输出单元702,选择子单元711、确定子单元712、统计子单元713、噪声能量差确定子单元714、获取子单元715以及噪声能量确定子单元716可由位于PC中的中央处理器(CPU,Central Processing Unit)、数字信号处理器(DSP,Digital SignalProcessor)、或可编程门阵列(FPGA,Field Programmable Gate Array)实现。
以上所述,仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。

Claims (10)

1.一种合成孔径雷达图像仿真的方法,其特征在于,所述方法包括:
根据实测图像获取指定系统灵敏度下的图像噪声能量差;
当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像;
或者,
根据仿真图像获取指定系统灵敏度下的图像噪声能量差;
当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像。
2.根据权利要求1所述的方法,其特征在于,所述根据实测图像获取指定系统灵敏度下的图像噪声能量差包括:
选择与待仿真合成孔径雷达系统的系统参数一致的实测图像;
根据所述系统参数确定所述实测图像的系统灵敏度;
根据所述实测图像统计所述实测图像的噪声能量;
根据所述实测图像的系统灵敏度与所述实测图像的噪声能量确定所述图像噪声能量差。
3.根据权利要求2所述的方法,其特征在于,所述根据所述实测图像的系统灵敏度与所述实测图像的噪声能量确定所述图像噪声能量差包括:
根据所述实测图像的系统灵敏度、所述实测图像的噪声能量得到目标功率与目标后向散射系数之比;
通过所述指定系统灵敏度、所述目标功率与目标后向散射系数之比确定得到指定系统灵敏度下的噪声能量;
根据所述指定系统灵敏度下噪声能量、所述实测图像的噪声能量确定所述图像噪声能量差。
4.根据权利要求1所述的方法,其特征在于,所述根据待仿真图像获取指定系统灵敏度下的图像噪声能量包括:
获取所述仿真图像的目标后向散射系数和目标功率;
根据所述目标后向散射系数、所述目标功率、所述指定系统灵敏度得到指定系统灵敏度下的图像噪声能量差。
5.根据权利要求1至4中的任一项所述的方法,其特征在于,所述利用所述图像噪声能量差得到系统灵敏度下的输出图像包括:将所述图像噪声能量差以正太分布噪声的形式加入到实测或仿真图像单视复数据的实部、虚部,得到指定系统灵敏度下的输出图像。
6.一种合成孔径雷达图像仿真的系统,其特征在于,所述系统包括:噪声能量获取单元、图像输出单元;其中,
所述噪声能量获取单元,用于根据实测图像获取指定系统灵敏度下的图像噪声能量差;或者,根据仿真图像获取指定系统灵敏度下的图像噪声能量差;
所述图像输出单元,用于当输出图像时,利用所述图像噪声能量差得到指定系统灵敏度下的输出图像。
7.根据权利要求6所述的系统,其特征在于,所述噪声能量获取单元,包括:选择子单元、确定子单元、统计子单元、噪声能量差确定子单元;其中,
所述选择子单元,用于选择与待仿真合成孔径雷达系统的系统参数一致的实测图像;
所述确定子单元,用于根据所述系统参数确定所述实测图像的系统灵敏度;
所述统计子单元,用于根据所述实测图像统计所述实测图像的噪声能量;
所述噪声能量差确定子单元,用于根据所述实测图像的系统灵敏度与所述实测图像的噪声能量确定需加入到实测图像的噪声能量差。
8.根据权利要求7所述的系统,其特征在于,所述噪声能量差确定单元,具体用于根据所述实测图像的系统灵敏度、所述实测图像的噪声能量得到目标功率与目标后向散射系数之比;
通过所述指定系统灵敏度、所述目标功率与目标后向散射系数之比确定得到指定系统灵敏度下的噪声能量;
根据所述指定系统灵敏度下的噪声能量、所述实测图像的噪声能量确定待加入到实测图像的噪声能量差。
9.根据权利要求6所述的系统,其特征在于,所述噪声能量获取单元包括:获取子单元、噪声能量确定子单元;其中,
所述获取子单元,用于获取所述仿真图像的目标后向散射系数和目标功率;
所述噪声能量确定子单元,用于根据所述目标后向散射系数、所述目标功率、所述指定系统灵敏度得到指定系统灵敏度下的图像噪声能量差。
10.根据权利要求6至9中任一项所述的系统,其特征在于,所述图像输出单元,具体用于当输出图像时,将所述图像噪声能量差以正太分布噪声的形式加入到实测或仿真单视复图像数据的实部、虚部,得到指定系统灵敏度下的输出图像。
CN201410261722.2A 2014-06-12 2014-06-12 一种合成孔径雷达图像仿真的方法及系统 Active CN104049243B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410261722.2A CN104049243B (zh) 2014-06-12 2014-06-12 一种合成孔径雷达图像仿真的方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410261722.2A CN104049243B (zh) 2014-06-12 2014-06-12 一种合成孔径雷达图像仿真的方法及系统

Publications (2)

Publication Number Publication Date
CN104049243A true CN104049243A (zh) 2014-09-17
CN104049243B CN104049243B (zh) 2016-08-17

Family

ID=51502329

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410261722.2A Active CN104049243B (zh) 2014-06-12 2014-06-12 一种合成孔径雷达图像仿真的方法及系统

Country Status (1)

Country Link
CN (1) CN104049243B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106646467A (zh) * 2016-11-09 2017-05-10 西南科技大学 一种合成孔径雷达图像噪声功率估计新方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EA014167B1 (ru) * 2005-09-08 2010-10-29 Бьорн А.Й. Ангельсен Способ построения изображения акустических свойств путем нелинейного низкочастотного управления свойствами высокочастотного рассеяния и распространения (варианты) и акустическое устройство для осуществления этого способа
CN103336272A (zh) * 2013-03-26 2013-10-02 中国科学院电子学研究所 一种基于几何结构的复杂目标sar图像仿真方法
CN103630897A (zh) * 2012-08-28 2014-03-12 中国科学院电子学研究所 一种多通道合成孔径雷达成像的方法
WO2014074241A1 (en) * 2012-11-07 2014-05-15 Axsun Technologies, Inc. Oct system with phase sensitive interference signal sampling

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EA014167B1 (ru) * 2005-09-08 2010-10-29 Бьорн А.Й. Ангельсен Способ построения изображения акустических свойств путем нелинейного низкочастотного управления свойствами высокочастотного рассеяния и распространения (варианты) и акустическое устройство для осуществления этого способа
CN103630897A (zh) * 2012-08-28 2014-03-12 中国科学院电子学研究所 一种多通道合成孔径雷达成像的方法
WO2014074241A1 (en) * 2012-11-07 2014-05-15 Axsun Technologies, Inc. Oct system with phase sensitive interference signal sampling
CN103336272A (zh) * 2013-03-26 2013-10-02 中国科学院电子学研究所 一种基于几何结构的复杂目标sar图像仿真方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘向阳等: "一种二维频域SAR回波仿真的改进方法", 《西安电子科技大学学报(自然科学版)》, vol. 40, no. 3, 30 June 2013 (2013-06-30) *
段秋萍: "合成孔径雷达系统参数设计、回波仿真及应用研究", 《中国优秀硕士学位论文全文数据库信息信息科技辑》, no. 10, 15 October 2008 (2008-10-15), pages 1 - 82 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106646467A (zh) * 2016-11-09 2017-05-10 西南科技大学 一种合成孔径雷达图像噪声功率估计新方法
CN106646467B (zh) * 2016-11-09 2019-05-14 西南科技大学 一种合成孔径雷达图像噪声功率估计新方法

Also Published As

Publication number Publication date
CN104049243B (zh) 2016-08-17

Similar Documents

Publication Publication Date Title
CN102279393B (zh) 一种基于多光谱传感器对高光谱传感器交叉辐射定标方法
CN101090301B (zh) 一种无线电波路径损耗仿真测量方法
CN104730503B (zh) 确定高分辨率sar参考目标rcs对定标影响的方法及补偿方法
CN109752696B (zh) 一种高分辨率合成孔径雷达卫星图像中角反射器rcs校正方法
EP3729857A1 (en) Radio coverage map generation
CN104199026A (zh) 基于线迹扫描二维近场成像的反向散射截面测量方法
CN107607945B (zh) 一种基于空间嵌入映射的扫描雷达前视成像方法
CN105929388A (zh) 一种基于Wi-Fi网络的新型室内定位方法
Ren et al. 3D Imaging Algorithm for Down‐Looking MIMO Array SAR Based on Bayesian Compressive Sensing
Geok et al. 3D RT adaptive path sensing Method: RSSI modelling validation at 4.5 GHz, 28 GHz, and 38 GHz
CN106772510B (zh) 一种基于载波相位测量的跳频测距方法
CN106021766B (zh) 基于机电集成优化的圆极化网状天线高指向精度设计方法
CN104407336B (zh) 一种方位敏感的目标电磁回波模拟方法
CN104049243A (zh) 一种合成孔径雷达图像仿真的方法及系统
CN103983948B (zh) 基于稀疏表示的目标角度距离联合估计方法
CN102608591A (zh) 一种基于相位修正的rcs外推方法
CN108061920A (zh) 探地雷达建模的方法
CN116184463A (zh) 一种宽刈幅海面风场反演方法、装置、设备和介质
CN110471030A (zh) 基于改进共轭梯度的射频层析成像被动式定位求解方法
CN104181530B (zh) 一种极化隔离度指标确定方法及装置
Hany et al. Local parametric approach of wireless capsule endoscope localization using randomly scattered path loss based WCL
Zhang et al. A Spaceborne SAR Calibration Simulator Based on Gaofen-3 Data
Zhang et al. A Cluster-Based Statistical Channel Model for Integrated Sensing and Communication Channels
CN113938365B (zh) THz通信探测一体化系统的能限检测方法以及装置
WO2023093825A1 (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
C14 Grant of patent or utility model
GR01 Patent grant