CN115113258B - 一种星载sar幅度闪烁条纹电离层闪烁参数测量方法 - Google Patents

一种星载sar幅度闪烁条纹电离层闪烁参数测量方法 Download PDF

Info

Publication number
CN115113258B
CN115113258B CN202210731818.5A CN202210731818A CN115113258B CN 115113258 B CN115113258 B CN 115113258B CN 202210731818 A CN202210731818 A CN 202210731818A CN 115113258 B CN115113258 B CN 115113258B
Authority
CN
China
Prior art keywords
amplitude
scintillation
image
fringe
measuring
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
Application number
CN202210731818.5A
Other languages
English (en)
Other versions
CN115113258A (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202210731818.5A priority Critical patent/CN115113258B/zh
Publication of CN115113258A publication Critical patent/CN115113258A/zh
Application granted granted Critical
Publication of CN115113258B publication Critical patent/CN115113258B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16Measuring radiation intensity
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • G01R29/0807Measuring electromagnetic field characteristics characterised by the application
    • G01R29/0814Field measurements related to measuring influence on or from apparatus, components or humans, e.g. in ESD, EMI, EMC, EMP testing, measuring radiation leakage; detecting presence of micro- or radiowave emitters; dosimetry; testing shielding; measurements related to lightning
    • G01R29/0857Dosimetry, i.e. measuring the time integral of radiation intensity; Level warning devices for personal safety use
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/12Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Engineering & Computer Science (AREA)
  • Electromagnetism (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种星载SAR幅度闪烁条纹电离层闪烁参数测量方法,包括:S1、对SAR幅度闪烁条纹幅度图像取对数,并进行镜像对称、二维快速傅里叶变换操作,得到图像频谱域;S2、在条纹能量方向上构造带阻高斯滤波器;S3、利用构造的带阻高斯滤波器对频域数据进行滤波处理,采用傅里叶逆变换将频域数据转换至图像域,得到双程的幅度闪烁误差;S4、利用周期图法进行单程强度功率谱估计;S5、取测量功率谱的低频成分进行拟合,测得闪烁强度CkL和谱指数p。本发明能够基于包含幅度闪烁条纹的星载SAR图像实现电离层闪烁参数的测量。

Description

一种星载SAR幅度闪烁条纹电离层闪烁参数测量方法
技术领域
本发明涉及天基微波遥感与地球空间天气的交叉技术领域,更具体地说,特别涉及一种星载SAR幅度闪烁条纹电离层闪烁参数测量方法。
背景技术
相较于更高波段的星载SAR系统,低波段(L波段和P波段)星载SAR一方面具有更强的穿透能力,能够实现植被或浅层地表下隐蔽目标的探测与成像;另一方面具有对生物量更高的敏感性,对研究全球碳循环与气候变化监测具有重要意义。但由于其工作频率低,低波段星载SAR容易受到电离层效应的影响。在特定的地磁场作用下,低纬地区的各向异性的电离层不规则体可导致SAR幅度闪烁条纹,该现象广泛呈现于L波段SAR卫星(如ALOSPALSAR、ALOS-2 PALSAR-2)图像中。
针对星载L波段SAR图像中的幅度闪烁条纹,目前的校正方法主要是在频域对幅度误差能量与低频地物能量的进行分离,最终实现条纹去除。而幅度闪烁误差估计结果作为条纹去除的衍生物,反映了电离层不规则体幅度闪烁效应对电磁波信号的影响程度,因此可以用于测量电离层闪烁参数。为此,有必要开发一种星载SAR幅度闪烁条纹电离层闪烁参数测量方法,基于幅度闪烁误差估计结果进行强度功率谱估计,然后再利用理论的功率谱与测量功率谱进行拟合,从而测量得到闪烁强度CkL以及谱指数p。
发明内容
本发明的目的在于提供一种星载SAR幅度闪烁条纹电离层闪烁参数测量方法,以克服现有技术所存在的缺陷。
为了达到上述目的,本发明采用的技术方案如下:
一种星载SAR幅度闪烁条纹电离层闪烁参数测量方法,包括以下步骤:
S1、对SAR幅度闪烁条纹幅度图像取对数,并进行镜像对称、二维快速傅里叶变换操作,得到图像频谱域;
S2、在条纹能量方向上构造带阻高斯滤波器;
S3、利用构造的带阻高斯滤波器对频域数据进行滤波处理,采用傅里叶逆变换将频域数据转换至图像域,得到双程的幅度闪烁误差;
S4、利用周期图法进行单程强度功率谱估计;
S5、取测量功率谱的低频成分进行拟合,测得闪烁强度CkL和谱指数p。
进一步地,所述步骤S1中对SAR幅度闪烁条纹幅度图像取对数包括:
首先、采用以下公式对输入图像M取幅度:
|M|=Zea
式中,Z表示图像地物目标信息,a0=ea表示幅度闪烁误差;
然后、对幅度图像取对数,得到信号ln|M|=lnZ+a。
进一步地,所述步骤S2中的带阻高斯滤波器由连续的数个带阻滤波器组成。
进一步地,所述步骤S4中利用周期图法进行单程强度功率谱估计的公式为:
式中,κr为电离层相位屏距离向波数,Δκr=2π/(dri·Nr)为距离波数域采样间隔,m=1,2,...,Nr,Nr为SAR图像距离向点数,dri=dr(Hs-Hi)/Hs为相位屏距离向采样间隔,Hi为电离层高度,FFTr()表示距离向快速傅里叶变换。
进一步地,所述步骤S5中拟合的公式为:
式中,||||表示求向量的范数,为球面波路径修正因子,λ为电磁波波长,θ为入射角,Hs为卫星高度。
与现有技术相比,本发明的优点在于:本发明提供的一种星载SAR幅度闪烁条纹电离层闪烁参数测量方法,能够基于包含幅度闪烁条纹的星载SAR图像实现电离层闪烁参数的测量。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明星载SAR幅度闪烁条纹电离层闪烁参数测量方法的流程图。
图2中(a)为一景没有受到电离层闪烁效应影响的PALSAR图像,(b)为仿真的双程幅度闪烁误差,(c)仿真的SAR幅度闪烁条纹图像。
图3中(a)为仿真的SAR幅度闪烁条纹图像幅度闪烁误差估计结果,(b)为残差。
图4为仿真的SAR幅度闪烁条纹图像强度功率谱拟合及参数测量结果。
图5为实测数据验证实验中所用的PALSAR幅度闪烁条纹图像。
图6为PALSAR幅度闪烁条纹图像幅度闪烁误差估计结果。
图7为PALSAR幅度闪烁条纹图像强度功率谱拟合及参数测量结果。
具体实施方式
下面结合附图对本发明的优选实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
输入受电离层幅度闪烁条纹影响的低波段星载SAR图像数据M。并且已知电磁波波长λ,入射角为θ,卫星高度Hs,星载SAR图像距离向采样间隔dr。
参阅图1所示,本实施例公开了一种星载SAR幅度闪烁条纹电离层闪烁参数测量方法,包括以下步骤:
步骤S1、对SAR幅度闪烁条纹幅度图像取对数,并进行镜像对称、二维快速傅里叶变换操作,得到图像频谱域。
具体的,对输入图像M取幅度:
|M|=Zea
式中,Z表示图像地物目标信息,a0=ea表示幅度闪烁误差。再做对数操作将幅度闪烁误差由乘性噪声转化为加性噪声,得到信号:
ln|M|=ln Z+a
为了避免傅里叶变换的边际效应,以图像中心进行镜像翻折,填充为原图像两倍尺寸。最后对图像做二维快速傅里叶变换,得到图像频谱域。
步骤S2、在条纹能量方向上构造带阻高斯滤波器。
具体的,在条纹能量方向上构造滤波器。该滤波器由连续的数个带阻滤波器组成,起点距离原点一定距离以避开低频的图像信息部分。另外由于图像进行了镜像翻折,产生了对称方向的条纹,因此滤波器也要添加对称的一条带阻滤波器。
优选地,本实施例采用高斯滤波器。
步骤S3、利用构造的带阻高斯滤波器对频域数据进行滤波处理,采用傅里叶逆变换将频域数据转换至图像域,得到双程的幅度闪烁误差。
具体的,利用构造的滤波器对频域图像做滤波处理,然后做傅里叶逆变换将频域数据转换至图像域,滤波后的图像为M′,滤掉的成分为幅度闪烁误差
注意,为双程幅度闪烁误差,单程幅度闪烁误差为/>单程强度闪烁误差为
步骤S4、利用周期图法进行单程强度功率谱估计。
具体的,利用周期图法进行单程强度功率谱估计:
式中,κr为电离层相位屏距离向波数Δκr=2π/(dri·Nr)为距离波数域采样间隔,m=1,2,...,Nr,Nr为SAR图像距离向点数,dri=dr(Hs-Hi)/Hs为相位屏距离向采样间隔,Hi为电离层高度,一般取经验值350公里,FFTr()表示距离向快速傅里叶变换。
步骤S5、取测量功率谱的低频成分进行拟合,测得闪烁强度CkL和谱指数p。
具体的,理论的强度功率谱可表示为:
式中,re为常数,代表经典电子半径,κ0=2π/L0为外尺度L0对应的波数,L0通常取10公里,为球面波路径修正因子,Γ()表示伽马函数。上式中仅CkL和p为未知参数。
由于幅度误差估计结果中仍然包含高频地物能量,因此取功率谱的低频成分进行拟合,拟合过程可描述为:
式中,||||表示求向量的范数。
图2(a)为一景没有受到电离层闪烁效应影响的PALSAR图像,可见PALSAR图像中无明显的幅度闪烁条纹现象。图2(b)为仿真的双程幅度闪烁误差,其均值为1.00,标准差为0.15,仿真参数:CkL=3×1034、log CkL=34.48,p=3,L0=10公里,且其他参数采用PALSAR系统参数。图2(c)为图2(a)与图2(b)的乘积,即仿真的SAR幅度闪烁条纹图像。
图3(a)为幅度闪烁条纹误差估计结果,其均值为1.01,标准差为0.13。图3(b)为幅度闪烁条纹误差估计结果与仿真的双程幅度闪烁误差之比,即为估计残差,其均值为1.01,标准差为0.07,标准差的下降说明了幅度误差估计结果与仿真注入的误差具有较好的一致性。
图4为强度功率谱拟合及参数测量结果。图中蓝色实线为利用周期图法估计得到的单程强度功率谱,绿色虚线为实际功率谱曲线,红色实线为拟合功率谱曲线,可见实际功率谱曲线与拟合功率谱曲线具有很好地一致性。另外,图中标出了电离层闪烁参数测量结果,测得CkL=2.63×1034或者log CkL=34.42,谱指数为p=3.10,与参数真值相近,从而验证了该方法的有效性。
图5为包含幅度闪烁条纹的实际PALSAR图像。该图像获取于2008年3月26日的亚马逊雨林地区,编号为ALPSRP115537040。
图6为PALSAR幅度闪烁条纹图像幅度闪烁误差估计结果,其形态与实际PALSAR图像中条纹形态相似。
图7为PALSAR幅度闪烁条纹图像强度功率谱拟合及参数测量结果,最终测得CkL=4.47×1033或者log CkL=33.65,谱指数为p=2.93。
虽然结合附图描述了本发明的实施方式,但是专利所有者可以在所附权利要求的范围之内做出各种变形或修改,只要不超过本发明的权利要求所描述的保护范围,都应当在本发明的保护范围之内。

Claims (4)

1.一种星载SAR幅度闪烁条纹电离层闪烁参数测量方法,其特征在于,包括以下步骤:
S1、对SAR幅度闪烁条纹幅度图像取对数,并进行镜像对称、二维快速傅里叶变换操作,得到图像频谱域;
S2、在条纹能量方向上构造带阻高斯滤波器;
S3、利用构造的带阻高斯滤波器对频域数据进行滤波处理,采用傅里叶逆变换将频域数据转换至图像域,得到双程的幅度闪烁误差;
S4、利用周期图法进行单程强度功率谱估计;
S5、取测量功率谱的低频成分进行拟合,测得闪烁强度CkL和谱指数p;
所述步骤S4中利用周期图法进行单程强度功率谱估计的公式为:
式中,κr为电离层相位屏距离向波数,Δκr=2π/(dri·Nr)为距离波数域采样间隔,m=1,2,...,Nr,Nr为SAR图像距离向点数,dri=dr(Hs-Hi)/Hs为相位屏距离向采样间隔,Hi为电离层高度,FFTr( )表示距离向快速傅里叶变换,dr为星载SAR图像距离向采样间隔,Hs为卫星高度,为单程强度闪烁误差。
2.根据权利要求1所述的星载SAR幅度闪烁条纹电离层闪烁参数测量方法,其特征在于,所述步骤S1中对SAR幅度闪烁条纹幅度图像取对数包括:
首先、采用以下公式对输入图像M取幅度:
|M|=Zea
式中,Z表示图像地物目标信息,a0=ea表示幅度闪烁误差;
然后、对幅度图像取对数,得到信号ln|M|=lnZ+a。
3.根据权利要求1所述的星载SAR幅度闪烁条纹电离层闪烁参数测量方法,其特征在于,所述步骤S2中的带阻高斯滤波器由连续的数个带阻滤波器组成。
4.根据权利要求1所述的星载SAR幅度闪烁条纹电离层闪烁参数测量方法,其特征在于,所述步骤S5中拟合的公式为:
式中,||·||表示求向量的范数,为球面波路径修正因子,λ为电磁波波长,θ为入射角,Hs为卫星高度。
CN202210731818.5A 2022-06-25 2022-06-25 一种星载sar幅度闪烁条纹电离层闪烁参数测量方法 Active CN115113258B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210731818.5A CN115113258B (zh) 2022-06-25 2022-06-25 一种星载sar幅度闪烁条纹电离层闪烁参数测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210731818.5A CN115113258B (zh) 2022-06-25 2022-06-25 一种星载sar幅度闪烁条纹电离层闪烁参数测量方法

Publications (2)

Publication Number Publication Date
CN115113258A CN115113258A (zh) 2022-09-27
CN115113258B true CN115113258B (zh) 2024-05-17

Family

ID=83329994

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210731818.5A Active CN115113258B (zh) 2022-06-25 2022-06-25 一种星载sar幅度闪烁条纹电离层闪烁参数测量方法

Country Status (1)

Country Link
CN (1) CN115113258B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118278218B (zh) * 2024-05-31 2024-08-27 中国科学院长春光学精密机械与物理研究所 大气扰动成像的数字孪生方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101124492A (zh) * 2005-01-31 2008-02-13 阿克伦大学 地震预测的氡监测系统
CN101806911A (zh) * 2010-01-14 2010-08-18 国家卫星气象中心 一种电离层闪烁监测方法及电离层闪烁监测仪
WO2018192151A1 (zh) * 2017-04-17 2018-10-25 苏州瑞派宁科技有限公司 一种闪烁脉冲数字化信号的拟合方法
CN110533665A (zh) * 2019-09-03 2019-12-03 北京航空航天大学 一种抑制扇贝效应和子带拼接效应的sar图像处理方法
CN114488157A (zh) * 2021-12-30 2022-05-13 中国人民解放军国防科技大学 低波段星载sar电离层闪烁相位误差方位向解耦合方法
CN114488156A (zh) * 2021-12-30 2022-05-13 中国人民解放军国防科技大学 一种低波段星载sar图像电离层闪烁效应校正方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101490981B1 (ko) * 2012-12-28 2015-02-09 서울시립대학교 산학협력단 위성레이더 간섭도의 이온왜곡 보정방법 및 그 장치

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101124492A (zh) * 2005-01-31 2008-02-13 阿克伦大学 地震预测的氡监测系统
CN101806911A (zh) * 2010-01-14 2010-08-18 国家卫星气象中心 一种电离层闪烁监测方法及电离层闪烁监测仪
WO2018192151A1 (zh) * 2017-04-17 2018-10-25 苏州瑞派宁科技有限公司 一种闪烁脉冲数字化信号的拟合方法
CN110533665A (zh) * 2019-09-03 2019-12-03 北京航空航天大学 一种抑制扇贝效应和子带拼接效应的sar图像处理方法
CN114488157A (zh) * 2021-12-30 2022-05-13 中国人民解放军国防科技大学 低波段星载sar电离层闪烁相位误差方位向解耦合方法
CN114488156A (zh) * 2021-12-30 2022-05-13 中国人民解放军国防科技大学 一种低波段星载sar图像电离层闪烁效应校正方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
The Influence of Equatorial Scintillation on L-Band SAR Image Quality and Phase;Franz J.Meyer .et al;IEEE Transactions on Geoscience and Remote Sensing;第54卷(第2期);第869-880页 *
低频段星载SAR闪烁效应误差分析与优化补偿方法;曾虹程 等;第七届高分辨率对地观测学术年会论文集;第147-164页 *

Also Published As

Publication number Publication date
CN115113258A (zh) 2022-09-27

Similar Documents

Publication Publication Date Title
Mertens et al. Improved upper limits on the 21 cm signal power spectrum of neutral hydrogen at z≈ 9.1 from LOFAR
Cloude Dual-baseline coherence tomography
Ferretti et al. Multibaseline InSAR DEM reconstruction: The wavelet approach
Chakraborty et al. Detailed study of ELAIS N1 field with the uGMRT–II. Source properties and spectral variation of foreground power spectrum from 300–500 MHz observations
Misra et al. Analysis of radio frequency interference detection algorithms in the angular domain for SMOS
CN102401892B (zh) 一种极化干涉合成孔径雷达系统性能评估方法
CN115113258B (zh) 一种星载sar幅度闪烁条纹电离层闪烁参数测量方法
CN103091665B (zh) 一种压缩感知合成孔径雷达射频干扰抑制处理方法
Zhai Inversion of organic matter content in wetland soil based on Landsat 8 remote sensing image
CN103258330A (zh) 一种高光谱图像端元丰度的估计方法
Belcher et al. Measurement of the ionospheric scintillation parameter $ C_ {k} L $ from SAR images of clutter
Shen et al. Removal of thin clouds using cirrus and QA bands of Landsat-8
CN106595863A (zh) 一种提高液晶可调谐滤波器光谱重建精度和分辨率的方法
CN103808413A (zh) 基于噪声分离的温度与比辐射率分离的方法及装置
CN111273289A (zh) 一种沙漠参数反演方法、装置、设备及存储介质
CN106646471A (zh) 基于方位空变误差补偿的机载高分辨sar成像方法
CN114677401B (zh) 一种基于极化雷达自身影像特征的水体提取方法及系统
CN111257854B (zh) 一种通用的基于遥感影像分割单元的地形校正优化方法
CN109143188B (zh) Tops哨兵-1数据电离层校正方法
Monti-Guarnieri et al. Passive sensing by Sentinel-1 SAR: Methods and applications
CN106248209A (zh) 一种基于仪器特征矩阵的干涉光谱仪光谱复原方法
CN114998126B (zh) 一种低波段星载sar图像电离层幅度闪烁条纹去除方法
Bovenga et al. Multichromatic analysis of satellite wideband SAR data
Oliva et al. Assessment of SMOS RFI Mitigation by means of a Triple collocation technique
Zhao et al. The synergetic estimation approach of forest above ground biomass based on X-band insar and P-band polsar data

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