Processing math: 0%
高级检索

基于深度图像先验的高光谱图像去噪方法

马飞, 王梓璇, 刘思雨

马飞, 王梓璇, 刘思雨. 基于深度图像先验的高光谱图像去噪方法[J]. 激光技术, 2024, 48(3): 379-386. DOI: 10.7510/jgjs.issn.1001-3806.2024.03.013
引用本文: 马飞, 王梓璇, 刘思雨. 基于深度图像先验的高光谱图像去噪方法[J]. 激光技术, 2024, 48(3): 379-386. DOI: 10.7510/jgjs.issn.1001-3806.2024.03.013
MA Fei, WANG Zixuan, LIU Siyu. Hyperspectral image denoising method based on depth image prior[J]. LASER TECHNOLOGY, 2024, 48(3): 379-386. DOI: 10.7510/jgjs.issn.1001-3806.2024.03.013
Citation: MA Fei, WANG Zixuan, LIU Siyu. Hyperspectral image denoising method based on depth image prior[J]. LASER TECHNOLOGY, 2024, 48(3): 379-386. DOI: 10.7510/jgjs.issn.1001-3806.2024.03.013

基于深度图像先验的高光谱图像去噪方法

基金项目: 

辽宁省教育厅科学研究面上项目 LJKZ0357

辽宁省自然科学基金资助项目 2023-MS-314

详细信息
    通讯作者:

    马飞, femircom@gmail.com

  • 中图分类号: TP751.1

Hyperspectral image denoising method based on depth image prior

  • 摘要: 为了避免现有的高光谱图像去噪优化模型仅考虑有限的高光谱内在结构特点、并未实现图像特征的精确表征的问题, 采用了一种基于空谱深度图像先验与平滑的高光谱图像去噪方法, 将紧框架变换与具有高表达与强学习能力的深度学习模型进行结合, 构建基于深度学习的噪声去除模型。首先在低秩矩阵分解的基础上, 利用特定的深度图像先验学习潜在的空谱特征; 然后分别构建端元与丰度矩阵的紧框架稀疏正则探究空谱局部平滑, 并解决深度图像先验的半拟合问题; 最后设计高效迭代算法实现模型求解。结果表明, 基于空谱深度图像先验的方法在各种复杂的噪声干扰下均表现出较好的视觉恢复性能, 峰值信噪比至少有1 dB以上的提升, 得到了高质量的恢复图像。该方法为高光谱图像去噪提供了参考。
    Abstract: In order to avoid that the problems of the existing hyperspectral image denoising optimization model only considers the limited intrinsic structure characteristics, and does not realize the accurate representation of image features, a denoising method based on spatial spectral depth image prior and smoothing is proposed. The model combines tight frame transform with a deep learning model with high expression and strong learning ability. Firstly, on the basis of low-rank matrix decomposition, the potential-spatial spectral features were learned by using specific depth images prior. Secondly, a tight frame of end and abundance matrix was constructed respectively to explore the local smoothing of the empty spectrum and solve the semi-fitting behavior of the depth image prior. Finally, an efficient iterative algorithm was designed to solve the model. The results show that the method based on space spectrum depth image prior has better performance under various complex noise interference, and the peak signal-to-noise ratio(PSNR) is improved by at least 1 dB, and high quality restored images are obtained. The method provides a reference for hyperspectral image denoising.
  • 激光半自动制导是一种典型的激光制导模式[1],其基本原理是通过位于载机或地面的激光目标指示器发射激光束照射目标,由位于弹头的探测器接收目标漫反射的激光信号, 探测器通过对接收波门内多个脉冲回波的处理以实现对目标的捕获, 实现目标捕获后,通过对目标持续确认以实现对目标的跟踪。

    其它激光目标指示器的散射信号和激光干扰机发射的信号都可能影响探测器的性能。前者是非恶意干扰,具有低重频干扰效果;后者通常工作在高重频工作模式,以保证干扰效果。干扰信号对探测器的捕获和跟踪环节都可能产生影响[2-3],具体的干扰效果与干扰机布设方位[4-5]、干扰信号重频[6-7]、干扰信号编码方式[8-9]有关,也与激光目标指示器信号波形[10]、探测器波门设计[11]、脉冲锁定技术[12]有关。

    激光制导抗干扰的技术途径主要可以分为两种: 一是对指示器发射信号进行编码,提升目标回波信号信息熵,降低干扰信号与目标回波的相似性[13];二是在目标捕获和跟踪处理环节充分挖掘目标回波信号的特征信息,剔除不满足目标特征的干扰信号[14]。研究者深入分析了编码波形抗干扰技术,并分别提出了性能良好的随机编码波形和抗干扰算法[15-17]。WU等人[18]提出了基于接收波门临近回波的高重频干扰排查方法,改善了抗干扰效果。CHEN等人[19]提出了基于干扰重频估计的高重频干扰剔除技术,可有效对抗固定重频的高重频干扰。CAO等人[20]提出了一种自适应扩展实时波门技术,有效避免了干扰脉冲进入波门导致目标脉冲漏检的情况。

    从原理上看,在满足激光目标指示器、探测器工程现实基础上,充分考虑感兴趣目标的行为特征和干扰机工作特点,开展发射波形精细化编码设计、脉间接收波门设计和抗干扰信号处理算法研究,可有效提升抗干扰性能。基于此,本文作者设计了一种随机重频发射信号,在收发准同步条件下,结合激光目标指示器、目标和探测器的运动特性,提出了端到端的基于自适应交叠波门的高重频干扰鉴别、目标探测方法, 仿真分析了该方法抗固定重频和随机重频干扰的性能。

    探测器通过对特定距离上目标脉冲多个回波的连续探测实现目标捕获,并通过连续接收目标脉冲回波的连续观测实现目标跟踪。受激光目标指示器、目标和探测器运动影响,不同脉冲回波通常在时延上存在差异。为了有效实现脉间回波关联,通常设定与激光目标指示器、目标和探测器先验信息相关的具有一定宽度的距离波门。

    当特定波门出现连续脉冲时,探测器进入脉冲锁定状态,并通过对锁定波门内的回波连续探测的关联实现目标捕获。干扰机通过发射高重频干扰信号,在探测器上产生大量密集虚假目标信号充斥于探测器接收波门之内,以较高的概率打断激光雷达目标捕获跟踪链路,严重影响激光制导效能。图 1是高重频干扰的原理图。图 1a中给出了无干扰时探测器接收目标回波信号的情况,此时目标回波信号被设置的波门较好地接收;图 1b中给出了高重频干扰信号的特征,即重复频率越高,在时间轴上分布越密集;图 1c中给出了在高重频干扰场景下探测器设置波门录取回波信号时,部分干扰信号进入了接收波门内。

    图 1 高重频干扰示意图
    图  1  高重频干扰示意图
    Figure  1.  High-repetition-frequency jamming schematic

    在特定漏射率、检测概率条件下,探测器连续观测到特定距离上N个发射信号的回波脉冲时,可认为成功捕获到目标,其中N与指示器漏射率、探测器目标检测概率等有关,在当前工艺下N为3或4,本文中取干扰方案中,通常将脉冲编码技术、时间波门技术、脉冲锁定技术这3种技术相结合来得到更好的抗干扰效果,如图 2所示。假设探测器接收到的首信号是激光目标指示器发射的第l个发射脉冲的回波,分析第l个~第l+3个发射脉冲的回波,目标回波序列为[1, 3, 6, 9]T,其中上标T表示向量的转置。由于高重频干扰信号超前于目标回波进入接收波门,接收波门内的首脉冲被干扰信号占据,接收下一个脉冲回波信号的实时波门向前偏移,导致在接收第l+3个脉冲时无法捕获目标回波7。根据首脉冲锁定技术,此时系统判定的目标回波序列[1, 2, 4, 6]T,而实际目标回波序列为[1, 3, 5, 7]T,此时表明高重频干扰信号已经成功干扰了制导系统。

    图 2 传统抗干扰方案
    图  2  传统抗干扰方案
    Figure  2.  Traditional anti-jamming approach

    在高重频情况下,干扰信号可能在目标回波信号之前就进入探测器接收波门,导致首个接收信号的可靠性下降。此时,以首脉冲为代表的传统脉冲锁定方法存在误判干扰信号的概率高以及易被牵引的问题。因此,在高重频干扰环境中,需要在波形设计、时间波门和信号处理技术方面开展端到端的优化。

    激光目标指示器以随机重频方式发射信号,即脉冲周期为\boldsymbol{T}=\left[T_1, T_2, \cdots, T_{k-1}\right]^{\mathrm{T}},如图 3所示。激光目标指示器在工作时间内发射k个信号,可划分k-1个重复周期,定义Tk为第k-1个和第k个脉冲的重复周期,k取值范围为1≤kk-1,应有Tmin < Tk < TmaxTminTmax分别为激光目标指示器发射信号的最小周期和最大周期,则激光指示器总的工作时间T_{\text {life }}=\sum\limits_{k=1}^{k-1} T_k。受传播损耗、目标散射特征、发射功率、探测器灵敏度等工程限制,通常情况下可探测回波最大时延τmax=2R/cTmin,其中c为光速,R为有效探测距离。发射信号、脉冲重复周期、最大探测距离等参数如图 3所示。

    图 3 发射信号与探测距离关系图
    图  3  发射信号与探测距离关系图
    Figure  3.  Transmission signal and detection range relationship graph

    假设发射脉冲信号波形为x(t),则发射波形为:

    s_0(t)=\sum\limits_{k=0}^{k-1} x\left(t-t_k\right) (1)

    式中: t_k=\sum\limits_{m=0}^k T_m,表示第k个脉冲相对第0个脉冲的发射迟延。第0个脉冲的发射时间为起始时间,T0=0。

    当干扰机侦察到照射源发射信号时,干扰机发射高重频干扰信号以影响探测器效能,在高重频干扰场景下,探测器接收到的信号可以表示为:

    s(t)=s_0{ }^{\prime}(t)+J(t)+n(t) (2)

    式中: s0(t)是发射脉冲信号对应的目标回波; J(t)为探测器接收到的干扰信号; n(t)为加性高斯白噪声。J(t)的表达式为:

    J(t)=\sum\limits_{i=1}^{\infty} A_i x\left(t-{\rm i} t_{\mathrm{jam}}{ }^{\prime}\right) (3)

    式中:tjam是高重频干扰信号重复周期; Ai是干扰机发射的第i个固定高重频干扰信号的幅度。此时探测器接收到的目标回波以及干扰信号如图 4所示。

    图 4 探测器接收到的信号示意图
    图  4  探测器接收到的信号示意图
    Figure  4.  Signal schematic diagram captured by detector

    在随机重频干扰场景下,接收机接收到的干扰信号的重复周期相较于固定重频干扰信号有随机抖动,此时J(t)的表达式为:

    J(t)=\sum\limits_{i=1}^{\infty} A_i x\left(t-\mathrm{i} t_{\mathrm{jam}}{ }^{\prime}-\Delta t_i\right) (4)

    式中: Δti为第i个随机重频干扰脉冲的随机时间抖动,传统的反向重频脉冲信号电路[19]无法有效对抗此类干扰。

    式(2)中的n(t)为加性高斯白噪声,可看作是出现概率较低、且呈现时空随机分布的背景回波,会使探测器出现虚警,在影响机理上看,与距离随机的低重频干扰相似。由于本文中考虑了随机脉冲重复周期情况下高重频J(t)干扰,可以认为噪声n(t)为特定J(t)的一部分,其影响机理可以在大样本统计实验中包含。在后续的分析中,将略去n(t)项的影响。

    传统目标捕获方法常采用固定波门和单脉冲锁定技术[11],即每个脉冲回波的波门宽度相等,而且每一个发射脉冲接收波门按前一个接收波门检测到的某一信号为时间基准设置。该技术易出现波门被高重频干扰引偏而丢失目标问题。与此不同,自适应交叠波门则通过首脉冲确定下一脉冲接收波门的开启时刻,通过末脉冲确定下一脉冲接收波门的关闭时刻[20]。本文作者在参考文献[20]的基础上,采用自适应交叠波门技术,提出基于先验信息的信号判决方法,实现目标有效捕获和抗干扰技术。

    对脉冲重复周期为T=[T1, T2, …, Tk-1]T的发射激光序列,假设第l个脉冲的接收波门内接收到的首信号时间为ts,末信号的时间为te,则可以实时对第l+1个脉冲的波门进行设置:

    \left\{\begin{array}{l} T_{\mathrm{g}, 1}=t_{\mathrm{s}}+\left(T_{l+1}-T_l\right)-\Delta t / 2 \\ T_{\mathrm{g}, 2}=t_{\mathrm{e}}+\left(T_{l+1}-T_l\right)+\Delta t / 2 \end{array}\right. (5)

    式中: Tg, 1Tg, 2是第l+1个脉冲实时波门的开启时刻和关闭时刻; TlTl+1分别是第l个脉冲和第l+1个脉冲距离其下一个脉冲的脉冲重复周期; Δt是波门宽度,如图 5所示。

    图 5 自适应交叠波门
    图  5  自适应交叠波门
    Figure  5.  Adaptive overlapped wavegate

    假设第l个脉冲发射之后,在第l个脉冲有效探测范围内接收到1号回波,以1号回波为基准时间,实时设置第l+1个脉冲的接收波门,并在此波门内接收到2、3号回波;以2号回波和3号回波为基准分别确定第l+2个脉冲接收波门的开启时刻和关闭时刻,并在此波门内接收到4、5、6号回波;以此递推下去,在第l+3个脉冲发射之后,探测器可以接收到7、8、9、10号回波,故将1号回波认为是疑似目标回波之后,通过自适应交叠波门得到的脉冲组合一共有4!=24种。接收机以接收4个发射脉冲的回波作为目标捕获的必要条件,在4个脉冲发射之后,在发射脉冲相对应的波门内收集到多组脉冲组合进入信号判决环节。

    确定自适应交叠波门后,可以得到多个包含目标、干扰和噪声(如上文所述,噪声被看作干扰信号的一部分)的多个回波组合。由于发射脉冲时序T=[T1, T2, …, Tk-1]T已知, 通过编码先验信息和探测器运动特性,可以有效地排除固定重频干扰信号和随机重频干扰信号。

    为简化问题,如第1.1节中所述,假定无干扰情况下,在4个接收波门中接收到至少3个满足时延特性的信号时, 可确认目标(具体探测脉冲数指标与漏设率、背景参数等决定,不是本文重点)。为不失一般性,假设打击目标静止不动,且激光目标指示器和探测器分别以2Ma的速度匀速飞向打击目标。假定某测试回波序列如图 5所示,此时的24种组合方式中,[1, 3, 5, 8]T组合是真实目标回波的组合,满足激光目标指示器和探测器运动特性,显然[1, 3, 6, 9]T不满足运动特性,可以直接鉴别为干扰组合。当有脉冲组合满足条件时,即可认为成功捕获到目标,否则回到目标捕获阶段,以第l+1个脉冲为起始判断脉冲,分析第l+1个脉冲至第l+4个脉冲的回波信号。算法流程图如图 6所示。

    图 6 目标判决流程图
    图  6  目标判决流程图
    Figure  6.  Target discrimination flowchart

    当接收到回波信号后,该算法依据激光探测器和激光目标指示器之间的时间同步关系,确定从第l个脉冲开始进行分析; 再根据发射脉冲先验信息及接收回波信号自适应生成第l个~第l+3个脉冲的交叠波门,以获得脉冲回波组合; 然后基于目标的运动特性对脉冲回波组合进行判决, 若判决结果符合运动规律,可以确定为已成功捕获目标;若未能捕获,则会将起始分析脉冲调整为第l+1个脉冲,并继续进行下一轮目标捕获。

    考虑到工程实际情况,设置激光目标指示器总工作时间Tlife=70 s,假定其从第l个脉冲开始检测回波信号后,判断第l个~第l+3个脉冲接收波门的回波信号是否为真实目标回波信号,为不失一般性,设置干扰信号与目标回波信号为01序列。

    在仿真实验中,激光目标指示器和探测器的速率以680 m/s向目标移动,激光目标指示器的最大有效照射距离为25 km,激光漏射率为1‰,发射脉冲脉宽为10 ns,单个脉冲目标检测概率为98%,基础波门宽度为4 μs,各脉冲的自适应交叠波门由第2.1节中的算法得到。本文中随机生成了1200个重复频率均匀分布在20 Hz~30 Hz之间的发射脉冲信号,激光目标指示器按照生成的随机序列发射激光。全时段发射高重频干扰信号的干扰机位于目标的100 m处。

    仿真生成与干扰信号混杂的目标回波信号,使用本文中算法对高重频干扰信号进行分析筛选,并计算出捕获到的目标的位置。当连续4个脉冲的回波信号符合运动规律时,即可认为目标捕获成功。为了验证本文中方法的性能,采用目标捕获时间(target acquisition time, TAT)指标进行性能评估,即:

    T_{\mathrm{a}}=\sum\limits_{i=T_{\mathrm{s}}}^{T_{\mathrm{e}}} T_i (6)

    式中: Ts为开始接收到目标回波的第1个脉冲间隔; Te为成功捕获到目标的最后一个脉冲间隔; Ta为捕获目标总消耗时间,Ta越小说明抗重频干扰算法性能越佳。

    设激光目标指示器漏射率为1‰、检测概率为98%,干扰参数设置如表 1所示。共设置12种重频干扰仿真场景,其中随机重频干扰在固定重频干扰的基础上加了随机时间抖动, 仿真结果如图 7~图 9所示。图 7表示某一脉冲发射后受到200 kHz高重频干扰后,探测器接收到的回波信号,其中标注的信号为真实目标回波,对应目标距离探测器10 km;图 8表示利用本文中算法成功去除掉高重频干扰后,得到的真实目标回波脉冲组合,其起始时间为第45个脉冲的发射时间;图 9为在不同频率的固定重频干扰和随机重频干扰下各进行104次仿真实验的平均目标捕获时间。

    表  1  干扰场景参数设置
    Table  1.  Jamming scenario parameters configuration
    jamming types transmitted signal jamming frequency/kHz experiment
    fixed frequency random frequency distributed between 20 Hz~30 Hz 0.02, 0.5, 10, 100, 200, 250 10000 trials per frequency
    下载: 导出CSV 
    | 显示表格
    图 7 受到200 kHz高重频干扰的目标回波
    图  7  受到200 kHz高重频干扰的目标回波
    Figure  7.  Target echo under 200 kHz high-repetition-frequency jamming
    图 8 从200 kHz高重频干扰中提取的目标回波
    图  8  从200 kHz高重频干扰中提取的目标回波
    Figure  8.  Target echo extracted from 200 kHz high-repetition-frequency jamming
    图 9 不同固定/随机频率干扰下目标平均捕获时间
    图  9  不同固定/随机频率干扰下目标平均捕获时间
    Figure  9.  Average target acquisition time under different fixed/random frequency jamming

    图 7图 8可以看出,基于自适应交叠波门的全脉冲回波探测及抗干扰方法能够有效提取面对200 kHz的高重频干扰时目标的制导回波,并通过迟延计算确定当前目标与探测器的距离。图 9展示了不同频率干扰下目标平均捕获时间。实验结果表明,多种干扰场景下平均目标捕获时间均在0.146 s左右; 从低频到高频,目标平均捕获时间有轻微的提升; 多种干扰场景下平均目标捕获时间均小于0.2 s,这说明该算法能够在开始检测的前4个脉冲中捕获到真实目标回波的回波序列。

    在多种干扰场景下的12×104次实验中,仿真结果未出现任何虚警情况,表明该算法具有良好的抗重频干扰性能。综上所述,本文中提出的抗干扰算法能够有效适应不同频率的固定重频干扰和随机重频干扰。

    根据发射脉冲的编码规则,本文中设计了接收脉冲的自适应交叠波门,并利用探测器运动的先验信息对自适应交叠波门的接收信号进行了合理性判决。本文中的算法充分利用已知的随机脉冲规则与运动先验信息,分析连续N个脉冲接收波门的信号。由于N个脉冲对应的自适应交叠波门接收到的固定重频干扰信号以及随机重频干扰持续满足真实目标的回波规律的可能性极小,故此方法可以有效筛选去除干扰信号,对固定/随机重频干扰均具有较强的抗干扰性能。

    在不同频率的固定/随机重频干扰场景下,利用仿真实验验证了本文中算法的有效性和优越性。仿真结果表明,在干扰场景下,随机重频激光制导雷达抗高重频干扰技术不仅可以成功捕获目标,而且目标捕获时间与无干扰场景下的捕获时间相当,保证了捕获目标的时效性,并且在12×104次统计试验中未出现虚警问题。此算法的优越性在于: 在高重频干扰环境下仍然具备高探测概率和低虚警概率,可用于提高激光制导雷达系统的抗干扰能力,提高其工作效率和可靠性。

  • 图  1   深度图像先验整体网络结构

    Figure  1.   Overall network flowchart of deep image prior

    图  2   Washington DC Mall数据集第100个波段在方案Ⅰ下的去噪图像

    Figure  2.   The denoised results of the 100th band in Washington DC Mall dataset under case Ⅰ

    图  3   Salinas数据集第100个波段在方案Ⅳ下的去噪图像

    Figure  3.   The denoised results of the 100th band in Salinas dataset under case Ⅳ

    图  4   在方案Ⅰ与Ⅲ下Washington DC Mall数据集所有波段的PSNR与ERGAS值

    Figure  4.   PSNR and ERGAS values of all bands in Washington DC Mall dataset under case Ⅰ and Ⅲ

    图  5   在方案Ⅱ与Ⅳ下Salinas数据集所有波段的PSNR与ERGAS值

    Figure  5.   PSNR and ERGAS values of all bands in Salinas dataset under case Ⅱ and Ⅳ

    图  6   WHU HongHu数据集第2个波段的去噪结果

    Figure  6.   The denoised results of the 2nd band in WHU HongHu dataset

    图  7   WHU HongHu数据集第2个波段的水平均值剖面图

    Figure  7.   Horizontal mean profile of the 2nd band in WHU HongHu dataset

    图  8   MPSNR与各参数值间的关系

    Figure  8.   Relationship between parameter values and MPSNR

    图  9   不同方案下MPSNR与迭代次数间的关系

    Figure  9.   Relationship between iteration number and MPSNR in different cases

    表  1   Washington DC Mall的去噪图像性能对比

    Table  1   Performance comparison of denoising images on Washington DC Mall

    case index LRMR LRTV HySuDeep DIP2D proposed
    MPSNR/dB 32.63 27.63 29.43 26.97 33.31
    MSSIM 0.926 0.859 0.907 0.849 0.938
    SAM 0.0837 0.127 0.153 0.121 0.069
    ERGAS 75.86 152.31 127.10 145.63 70.14
    MPSNR/dB 32.93 32.12 29.72 27.59 33.81
    MSSIM 0.934 0.928 0.913 0.844 0.948
    SAM 0.082 0.094 0.150 0.129 0.065
    ERGAS 73.46 90.12 125.68 140.33 66.39
    MPSNR/dB 31.43 31.96 29.54 32.31 34.15
    MSSIM 0.925 0.934 0.913 0.934 0.949
    SAM 0.101 0.083 0.151 0.071 0.065
    ERGAS 92.39 89.38 126.43 78.87 63.85
    MPSNR/dB 30.24 26.42 28.92 25.04 32.73
    MSSIM 0.906 0.829 0.901 0.789 0.931
    SAM 0.114 0.166 0.156 0.154 0.075
    ERGAS 103.71 181.25 136.37 187.67 74.83
    下载: 导出CSV

    表  2   Salinas的去噪图像性能对比

    Table  2   Performance comparison of denoising images on Salinas

    scase index LRMR LRTV HySuDeep DIP2D proposed
    MPSNR/dB 34.69 29.26 31.98 29.62 36.03
    MSSIM 0.848 0.778 0.907 0.763 0.895
    SAM 0.040 0.089 0.069 0.068 0.034
    ERGAS 89.86 628.62 107.14 338.15 66.70
    MPSNR/dB 34.35 33.95 32.22 32.01 35.86
    MSSIM 0.857 0.880 0.917 0.823 0.894
    SAM 0.043 0.068 0.069 0.053 0.035
    ERGAS 186.45 583.89 112.05 271.69 65.24
    MPSNR/dB 30.68 31.42 32.16 35.39 36.69
    MSSIM 0.815 0.855 0.926 0.884 0.905
    SAM 0.083 0.088 0.071 0.035 0.033
    ERGAS 127.43 109.72 86.83 78.07 55.54
    MPSNR/dB 29.19 26.40 31.04 27.92 35.20
    MSSIM 0.767 0.715 0.884 0.728 0.881
    SAM 0.090 0.127 0.074 0.093 0.036
    ERGAS 197.29 961.37 154.44 583.71 91.20
    下载: 导出CSV

    表  3   不同张量核范数方法的性能对比

    Table  3   Performance comparison of different tensor nuclear norm method

    data case method MPSNR/dB MSSIM SAM ERGAS
    Washington DC Mall SSDIP 31.69 0.931 0.087 86.25
    SSDSN 32.41 0.935 0.071 78.03
    SSDFS 31.90 0.926 0.096 89.81
    proposed 33.81 0.948 0.065 66.39
    SSDIP 26.53 0.849 0.124 157.01
    SSDSN 31.89 0.925 0.075 82.64
    SSDFS 27.21 0.850 0.124 155.91
    proposed 33.11 0.939 0.070 71.68
    Salinas SSDIP 33.87 0.863 0.054 537.44
    SSDSN 35.37 0.894 0.036 67.93
    SSDFS 34.38 0.869 0.048 196.51
    proposed 36.18 0.901 0.035 67.72
    SSDIP 28.78 0.777 0.083 499.51
    SSDSN 34.63 0.870 0.042 95.88
    SSDFS 29.22 0.787 0.082 511.01
    proposed 35.20 0.881 0.036 91.20
    下载: 导出CSV

    表  4   不同方法在4种方案下的平均运行时间对比/s

    Table  4   Comparison of average running times of different methods under four cases/s

    data LRMR LRTV HySuDeep DIP2D proposed
    Washington DC Mall 74.56 55.69 29.12 59.91 28.82
    Salinas 65.39 47.72 26.47 52.65 27.53
    下载: 导出CSV
  • [1] 赵佳乐, 王广龙, 周冰, 等. 基于边缘剔除的陆基高光谱图像噪声评估方法[J]. 激光技术, 2023, 47(1): 121-126. DOI: 10.7510/jgjs.issn.1001-3806.2023.01.019

    ZHAO J L, WANG G L, ZHOU B, et al. Noise evaluation method for land-based hyperspectral images based on edge elimination[J]. Laser Technology, 2023, 47(1): 121-126(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2023.01.019

    [2] 刘翠连, 陶于祥, 罗小波, 等. 混合卷积神经网络的高光谱图像分类方法[J]. 激光技术, 2022, 46(3): 355-361. DOI: 10.7510/jgjs.issn.1001-3806.2022.03.009

    LIU C L, TAO Y X, LUO X B, et al. Hyperspectral image classification based on hybrid convolutional neural network[J]. Laser Techno-logy, 2022, 46(3): 355-361(in Chinese). DOI: 10.7510/jgjs.issn.1001-3806.2022.03.009

    [3]

    WU Z B, SUN J, ZHANG Y, et al. Recent developments in parallel and distributed computing for remotely sensed big data processing[J]. Proceedings of the IEEE, 2021, 109(8): 1282-1305. DOI: 10.1109/JPROC.2021.3087029

    [4]

    XU X S, GAO R X, QING Y, et al. Hyperspectral image mixed noise removal via tensor robust principal component analysis with tensor-ring decomposition[J]. International Journal of Remote Sensing, 2023, 44(5): 1556-1578. DOI: 10.1080/01431161.2023.2187720

    [5]

    SUN H, ZHENG X T, LU X Q. A supervised segmentation network for hyperspectral image classification[J]. IEEE Transactions on Ⅰ-mage Processing, 2021, 30(4): 2810-2825.

    [6]

    RAO W Q, GAO L R, QU Y, et al. Siamese transformer network for hyperspectral image target detection[J]. IEEE Transactions on Geoscience And Remote Sensing 2022, 60(4): 1-19.

    [7]

    DAO P D, MANTRIPRAGADA K, HE Y, et al. Improving hyperspectral image segmentation by applying inverse noise weighting and outlier removal for optimal scale selection[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2021, 171(1): 348-366.

    [8]

    ZHANG P, NING J. Hyperspectral image denoising via group sparsity regularized hybrid spatio-spectral total variation[J]. Remote Sensing, 2022, 14(10): 2348-2372. DOI: 10.3390/rs14102348

    [9]

    XIONG F C, ZHOU J, ZHAO Q L, et al. MAC-Net: Model-aided nonlocal neural network for hyperspectral image denoising[J]. IEEE Transactions on Geoscience and Remote Sensing, 2021, 60(11): 1-14.

    [10]

    MA F, HUO S, YANG F X. Graph-based logarithmic low-rank tensor decomposition for the fusion of remotely sensed images[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2021, 14(4): 11271-11286.

    [11]

    ZENG Z Y, HUANG T Z, CHEN Y, et al. Nonlocal block-term decomposition for hyperspectral image mixed noise removal[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2021, 14(2): 5406-5420.

    [12]

    WANG Y, PENG J J, ZHAO Q, et al. Hyperspectral image restoration via total variation regularized low-rank tensor decomposition[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2017, 11(4): 1227-1243.

    [13]

    XU C G, WU Z M, LI F, et al. Spectral-spatial joint sparsity unmixing of hyperspectral images based on framelet transform[J]. Infrared Physics & Technology, 2021, 112(5): 103564.

    [14]

    KONG W F, SONG Y Y, LIU J. Hyperspectral image denoising via framelet transformation based three-modal tensor nuclear norm[J]. Remote Sensing, 2021, 13(19): 3829-3852.

    [15]

    XU X, LI J, LI S T, et al. Curvelet transform domain-based sparse nonnegative matrix factorization for hyperspectral unmixing[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2020, 13(3): 4908-4924. http://doc.paperpass.com/foreign/rgArti2020228667341.html

    [16]

    PAN E, MA Y, MEI X G, et al. D2net: deep denoising network in frequency domain for hyperspectral image[J]. IEEE/CAA Journal of Automatica Sinica, 2022, 10(3): 813-815.

    [17]

    ZHANG H Y, HE W, ZHANG L P, et al. Hyperspectral image restoration using low-rank matrix recovery[J]. IEEE Transactions on Geoscience and Remote Sensing, 2013, 52(8): 4729-4743.

    [18]

    HE W, ZHANG H Y, ZHANG L P, et al. Total-variation-regula-rized low-rank matrix factorization for hyperspectral image restoration[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 54(1): 178-188.

    [19]

    ZHUANG L N, NG M K, FU X Y. Hyperspectral image mixed noise removal using subspace representation and deep cnn image prior[J]. Remote Sensing, 2021, 13(20): 2071-2092.

    [20]

    SIDOROV O, HARDEBERG J. Deep hyperspectral prior: single-image denoising, inpainting, super-resolution[C]//Proceedings of the IEEE/CVF International Conference on Computer Vision Workshops. New York, USA: IEEE, 2019: 3844-3851.

图(9)  /  表(4)
计量
  • 文章访问数:  10
  • HTML全文浏览量:  0
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-06-25
  • 修回日期:  2023-07-26
  • 发布日期:  2024-05-24

目录

/

返回文章
返回