CN106646467B - 一种合成孔径雷达图像噪声功率估计新方法 - Google Patents
一种合成孔径雷达图像噪声功率估计新方法 Download PDFInfo
- Publication number
- CN106646467B CN106646467B CN201610983515.7A CN201610983515A CN106646467B CN 106646467 B CN106646467 B CN 106646467B CN 201610983515 A CN201610983515 A CN 201610983515A CN 106646467 B CN106646467 B CN 106646467B
- Authority
- CN
- China
- Prior art keywords
- noise power
- estimation
- image
- parasitism
- estimates
- 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.)
- Active
Links
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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种估计SAR图像噪声功率的新方法。在SAR图像超分辨率重建中,噪声功率估计是一个非常重要的问题。本发明提供了一种SAR图像噪声功率估计新方法。本发明首先定义了自洽度与无寄生测度两个新的参量,然后利用这两个新参量的乘积定义了一种联合参量,利用该联合参量即可实现对雷达图像噪声功率的准确估计。与现有技术相比,本发明的有益效果是:(1)利用本发明方法对SAR图像噪声功率进行估计,同时考虑了热噪声和非热噪声,估计结果更加全面准确;(2)本发明方法实现了对雷达图像噪声功率的自动估计,无需人工根据雷达接收机温度进行计算确定;(3)利用本发明方法求得的噪声功率对SAR图像进行超分辨率重建,能够得到更好的重建图像。
Description
技术领域
本发明涉及一种估计雷达图像噪声功率的新方法。
背景技术
雷达图像超分辨率重建中,通常需要使用雷达成像系统点扩散函数(PSF)和平均噪声功率等先验知识。雷达PSF可以利用信号带宽及雷达接收机匹配滤波器中降低副瓣的窗函数等知识进行估计,也可通过记录强点散射体(如三面角反射器)的散射数据进行测定。雷达系统噪声功率通常建模为雷达接收机热噪声,噪声的大小取决于接收机的温度,所以理论上如果温度已知,则可计算出噪声功率。但实际系统中的噪声还应包含发射机/接收机非线性、气相干扰以及模型未曾考虑的其它寄生杂散效应等非热噪声的影响。对于宽带雷达,这是尤其需要考虑的情况。
在雷达图像超分辨率重建中,如果假定的噪声太小,则超分辨率重建算法将过拟合数据,这将导致估计场景图像中寄生能量很高但相位相反的交替象素。PSF矩阵的平滑作用可以消除此类现象并产生近乎实测数据的图像,但估计场景将不再是真实场景的准确重建。如果假定的噪声太大,则不管场景包含的是点目标、展布式目标还是杂波,估计场景都将由非常少的点目标构成,这是因为实测数据中任何对正确PSF的偏离,都被认为是噪声而非目标结构导致的。由于我们不知道(或不能测出)非热噪声的影响,因此需要根据实测图像数据本身估计总体系统噪声。
假定雷达目标具有高分辨率二维后向散射系数,当从某一角度对目标进行成像时,将产生一幅特定的实测图像。雷达成像过程可等效为雷达系统PSF与高分辨率目标场景的卷积,在数学上可表示为:
g=Tf+n
这里f为高分辨率目标场景;T是反映雷达成像系统PSF作用的矩阵(即成像算子);g为低分辨率实测图像;n为系统噪声。
设潜在场景的一个估计为fest,则对于每个距离单元,其平均噪声功率的一种估计为:
这里N是距离单元的数目。
说明书附图1为利用上述噪声功率估计算法,对某多目标一维场景噪声功率(相对于最强目标,加入了-30dB的噪声功率;假定目标多普勒频率已知)进行估计所得的结果与输入给上述算法的假定噪声功率的对比曲线(实线),曲线通过多次运行上述噪声功率估计算法得到,且假定噪声功率在-140~40dB范围内。图中同时给出了一条一致线(自洽线),在该线上估计噪声功率等于假定噪声功率。从说明书附图1可以看出,当假定噪声功率处于-30dB~-45dB之间时,由算法获得的噪声功率与假定噪声功率近似一致,该噪声功率区间包含了-30dB的真实噪声功率值。
发明内容
上述噪声功率估计算法只能确定图像真实噪声功率的大致范围,不能确定图像真实噪声功率的具体准确值。为了克服上述噪声功率估计算法的不足,本发明提供了一种新颖的雷达图像噪声功率估计方法,以实现对雷达图像噪声功率的准确估计。
本发明所采用的具体技术方案即噪声功率估计方法如下:
假定噪声的标准偏差是σin,我们首先定义一个自洽度参量:
其中σest为噪声功率估计值的标准偏差。显然,0≤SC≤1。当输入的假定噪声功率与估计噪声功率相等时,SC达到最大值。对一系列输入噪声功率进行测量,获得的自洽度曲线如说明书附图2中的实线所示。可以看出,如果利用自洽度作为估计真实噪声功率大小的唯一依据,将存在不确定性(可能存在多个峰值),还需进一步确定真实噪声功率为-30dB还是-45dB。确定估计场景数据是否过拟合的一种方法是比较估计场景与实测图像的功率大小。当数据被过拟合时,估计场景中将寄生能量很大但相位相反的交替像素,这些像素通常具有比实测数据更大的幅值。这里我们定义一个无寄生测度参量进行度量(这里,单个像素的无寄生测度参量用SFpixel(i)表示,整幅场景图像的无寄生测度用SF表示):当估计场景功率小于实测图像功率时,SFpixel(i)等于1,否则SFpixel(i)等于实测图像功率与估计场景功率的比值,其数学表达式如下:
其中|fest(i)|为估计场景fest中第i个像素的功率;|g(i)|为实测图像g中第i个像素的功率;SFpixel(i)为第i个像素的无寄生测度,显然,该值处于0~1之间。对所有像素的无寄生测度SFpixel(i)求平均即得整个场景图像的无寄生测度SF:
其中N为估计场景图像和实测图像的像素个数。当SF接近1时,数据不会过拟合,估计场景中将不会寄生能量很大但相位相反的交替像素。SF与输入给算法的假定噪声功率的比较曲线如说明书附图2中的点划线所示。当假定的噪声功率很大时,SF近似等于1;随着假定噪声功率的减小,当减小到低于实际噪声功率-30dB时,SF迅速下降到0.1;对于所有低于该值的假定噪声功率,SF基本保持为0.1不变。这说明SF值对导致估计场景寄生交替像素的低假定噪声否决效果良好。
本发明定义一种联合参量对噪声功率进行估计,这种联合参量定义为上述自洽度与无寄生测度两个参量的乘积:
CQ=SC×SF
联合参量曲线如说明书附图2中的虚线所示,从该曲线可以看出,在真实噪声功率-30dB处具有唯一明确的峰值。利用这种噪声功率估计方法,对具有不同目标数目和不同实际噪声功率的图像的噪声功率进行估计,发现噪声功率估计结果与实际噪声功率均基本一致,误差很小,通常在几个dB之内。
现有技术将雷达成像系统噪声功率建模为雷达接收机热噪声,根据接收机的温度人工确定图像噪声功率的大小,但实际雷达成像系统特别是宽带雷达成像系统中的噪声还应包含发射机/接收机非线性、气相干扰以及其它寄生杂散效应等非热噪声的影响。因此,与现有技术相比,本发明的有益效果是:(1)利用本文方法对雷达图像噪声功率进行估计,同时考虑了热噪声和非热噪声,估计结果更全面准确;(2)本文方法实现了对雷达图像噪声功率的自动估计,无需人工根据雷达接收机温度进行计算确定;(3)利用本文方法求得的噪声功率,能得到更好的超分辨率重建图像。
附图说明
图1为估计噪声功率与假定噪声功率对比曲线;图2为噪声功率估计的优化图。
具体实施方式
为了使本发明的技术手段、创作特征、工作流程、使用方法达成目的与功效易于明白了解,下面结合说明书附图2对本发明进一步说明。
本发明所提出的噪声功率估计方法的实现步骤如下:
(1)首先定义一个自洽度参量如下:
其中σin为假定噪声的标准偏差,σest为噪声功率估计值的标准偏差,0≤SC≤1。当输入的假定噪声功率与估计噪声功率相等时,SC达到最大值。对一系列输入噪声功率进行测量,获得的自洽度曲线如说明书附图2中的实线所示。可以看出,利用自洽度作为估计真实噪声功率大小的唯一依据,将存在不确定性,还需进一步确定真实噪声功率为-30dB还是-45dB。
(2)然后定义一个无寄生测度参量(这里,单个像素的无寄生测度参量用SFpixel(i)表示,整幅场景图像的无寄生测度用SF表示):当估计场景功率小于实测图像功率时,SFpixel(i)等于1,否则SFpixel(i)等于实测图像功率与估计场景功率的比值,其数学表达式如下:
其中|fest(i)|为估计场景fest中第i个像素的功率;|g(i)|为实测图像g中第i个像素的功率;SFpixel(i)为第i个像素的无寄生测度,显然,该值处于0~1之间。对所有像素的SFpixel(i)求平均,即可得到整个场景的无寄生测度SF:
其中N为估计场景图像和实测图像的像素个数,无寄生测度SF与输入给算法的假定噪声功率的关系曲线如说明书附图2中的点划线所示。当假定的噪声功率很大时,SF近似等于1;随着假定噪声功率的减小,当减小到低于实际噪声功率-30dB时,SF迅速下降到0.1;对于所有低于该值的假定噪声功率,SF基本保持为0.1不变。这说明无寄生测度SF值对导致估计场景寄生交替像素的低假定噪声否决效果良好。
(3)最后定义一种联合参量(自洽度与无寄生测度的乘积)对噪声功率进行估计:
CQ=SC×SF
联合参量曲线如说明书附图2中的虚线所示,从该曲线可以看出,在真实噪声功率-30dB处具有唯一明确的峰值,从而成功实现了对图像真实噪声功率的估计。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (1)
1.一种合成孔径雷达图像噪声功率估计新方法,与现有的利用估计场景图像确定合成孔径雷达图像平均噪声功率方法不同,其特征在于,本方法利用自洽度和无寄生测度的乘积构成的联合参量对图像噪声功率进行估计,获得的联合参量曲线在真实噪声功率处具有唯一明确的峰值,从而实现对图像真实噪声功率的估计;首先定义一个自洽度参量 其中σest为噪声功率估计值的标准偏差,σin为假定噪声的标准偏差,对一系列输入噪声功率进行测量,即可获得自洽度曲线;利用自洽度作为估计真实噪声功率大小的依据,当自洽度曲线存在多个峰值时,将存在不确定性,为了解决不确定性的问题,接着定义单个像素无寄生测度和整个场景图像无寄生测度两个参量,单个像素的无寄生测度其中|fest(i)|为估计场景fest中第i个像素的功率,|g(i)|为实测图像g中第i个像素的功率,SFpixel(i)为第i个像素的无寄生测度,对所有像素的无寄生测度SFpixel(i)求平均,得到整个场景图像的无寄生测度其中N为估计场景图像和实测图像的像素个数,SF为整个场景图像的无寄生测度;最后利用自洽度和无寄生测度的乘积构成的联合参量CQ=SC×SF对图像噪声功率进行估计,获得的联合参量曲线仅在真实噪声功率处具有唯一明确的峰值,从而成功实现了对图像真实噪声功率的估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610983515.7A CN106646467B (zh) | 2016-11-09 | 2016-11-09 | 一种合成孔径雷达图像噪声功率估计新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610983515.7A CN106646467B (zh) | 2016-11-09 | 2016-11-09 | 一种合成孔径雷达图像噪声功率估计新方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106646467A CN106646467A (zh) | 2017-05-10 |
CN106646467B true CN106646467B (zh) | 2019-05-14 |
Family
ID=58805701
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610983515.7A Active CN106646467B (zh) | 2016-11-09 | 2016-11-09 | 一种合成孔径雷达图像噪声功率估计新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106646467B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108693528B (zh) * | 2018-04-13 | 2021-04-20 | 北博(厦门)智能科技有限公司 | 雷达实时信号处理器 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1213076A (zh) * | 1998-09-30 | 1999-04-07 | 复旦大学 | 一种红外双重富利埃变换的椭圆偏振光谱仪及其设计方法 |
CN101598783A (zh) * | 2009-07-08 | 2009-12-09 | 西安电子科技大学 | 基于ppca模型的强噪声背景下雷达距离像统计识别方法 |
CN102158637A (zh) * | 2011-04-01 | 2011-08-17 | 西安电子科技大学 | 基于Surfacelet变换域的空间自适应阈值视频去噪方法 |
JP2014006082A (ja) * | 2012-06-21 | 2014-01-16 | Furuno Electric Co Ltd | ノイズ抑圧装置、ノイズ抑圧方法およびノイズ抑圧プログラム |
CN104049243A (zh) * | 2014-06-12 | 2014-09-17 | 中国科学院电子学研究所 | 一种合成孔径雷达图像仿真的方法及系统 |
-
2016
- 2016-11-09 CN CN201610983515.7A patent/CN106646467B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1213076A (zh) * | 1998-09-30 | 1999-04-07 | 复旦大学 | 一种红外双重富利埃变换的椭圆偏振光谱仪及其设计方法 |
CN101598783A (zh) * | 2009-07-08 | 2009-12-09 | 西安电子科技大学 | 基于ppca模型的强噪声背景下雷达距离像统计识别方法 |
CN102158637A (zh) * | 2011-04-01 | 2011-08-17 | 西安电子科技大学 | 基于Surfacelet变换域的空间自适应阈值视频去噪方法 |
JP2014006082A (ja) * | 2012-06-21 | 2014-01-16 | Furuno Electric Co Ltd | ノイズ抑圧装置、ノイズ抑圧方法およびノイズ抑圧プログラム |
CN104049243A (zh) * | 2014-06-12 | 2014-09-17 | 中国科学院电子学研究所 | 一种合成孔径雷达图像仿真的方法及系统 |
Non-Patent Citations (2)
Title |
---|
"Super-resolution reconstruction of synthetic aperture radar image using adaptive threshold singular value decomposition technique";Zhu Z W 等;《Journal of Central South University》;20111231;第18卷(第3期);第809-815页 |
"基于联合聚焦/超分辨贝叶斯模型的雷达目标超分辨重建";朱正为 等;《系统工程与电子技术》;20110630;第33卷(第6期);第1261-1264页 |
Also Published As
Publication number | Publication date |
---|---|
CN106646467A (zh) | 2017-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110488393B (zh) | X波段双偏振天气雷达定量测量降水方法及系统 | |
US7132974B1 (en) | Methods and systems for estimating three dimensional distribution of turbulence intensity using radar measurements | |
CN107167781B (zh) | 海杂波幅度对数正态分布参数的分位点估计方法 | |
JP5932643B2 (ja) | 同一領域上で取得されたsar画像から得られるインターフェログラムのフィルタリング処理方法 | |
CN102567996A (zh) | 一种遥感卫星图像质量评价的方法 | |
CN105093925A (zh) | 一种基于被测地形特点的机载激光雷达参数实时自适应调整方法与装置 | |
Schwarz et al. | Exponential decomposition with implicit deconvolution of lidar backscatter from the water column | |
JP7308869B2 (ja) | 航空機搭載レーダを用いて波高を測定するための方法 | |
CN108089166A (zh) | 一种基于毫米波测云雷达的自适应频域检测方法 | |
US20230081137A1 (en) | Spatial metrics for denoising depth image data | |
US11789142B2 (en) | Graph-based array signal denoising for perturbed synthetic aperture radar | |
US8854250B2 (en) | Off-diagonal element echo power estimator for polarization weather radar | |
EP2817655B1 (en) | Systems and methods for image sharpening | |
CN108535719A (zh) | 基于多普勒频谱重心校正的调频连续波着陆雷达测速方法 | |
US7339519B2 (en) | Methods and apparatus for target radial extent determination using deconvolution | |
CN106646467B (zh) | 一种合成孔径雷达图像噪声功率估计新方法 | |
Treece et al. | Ultrasound attenuation measurement in the presence of scatterer variation for reduction of shadowing and enhancement | |
CN109917383A (zh) | 基于回波降采样积累的低信噪比isar成像方法 | |
CN110261837A (zh) | 一种基于航迹信息的复杂目标rcs计算方法 | |
JP2005043333A (ja) | レーダ画像からの波浪方向スペクトル逆推定方法及びシステム | |
JP3626089B2 (ja) | ウィンドプロファイラにおける信号処理装置および信号処理方法 | |
CN108845301A (zh) | 一种用于双基地isar的目标等效旋转中心估计方法 | |
Jones et al. | A method to quantify bedform height and asymmetry from a low-mounted sidescan sonar | |
CN112698330B (zh) | 基于比值对数积的极化sar图像非局部均值方法 | |
CN114252878A (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 |