CN111681156B - 一种应用于无线传感器网络的深度压缩感知图像重构方法 - Google Patents

一种应用于无线传感器网络的深度压缩感知图像重构方法 Download PDF

Info

Publication number
CN111681156B
CN111681156B CN202010533495.XA CN202010533495A CN111681156B CN 111681156 B CN111681156 B CN 111681156B CN 202010533495 A CN202010533495 A CN 202010533495A CN 111681156 B CN111681156 B CN 111681156B
Authority
CN
China
Prior art keywords
image
wireless sensor
sensor network
image reconstruction
crop
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
CN202010533495.XA
Other languages
English (en)
Other versions
CN111681156A (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.)
Nankai University
Original Assignee
Nankai 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 Nankai University filed Critical Nankai University
Priority to CN202010533495.XA priority Critical patent/CN111681156B/zh
Publication of CN111681156A publication Critical patent/CN111681156A/zh
Application granted granted Critical
Publication of CN111681156B publication Critical patent/CN111681156B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T1/00General purpose image data processing
    • G06T1/0007Image acquisition
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30181Earth observation
    • G06T2207/30188Vegetation; Agriculture

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Image Processing (AREA)

Abstract

一种应用于无线传感器网络的深度压缩感知图像重构方法。本发明属于信号处理技术领域,具体涉及一种应用于无线传感器网络的农作物图像重构方法。本发明深入发掘农业物联网中农作物图像相似性,提出了一种基于生成模型的深度压缩感知农作物图像重构方法。该方法同时利用了农作物图像的稀疏性与低秩性,以及实际应用中图像的单一性,结合深度卷积生成对抗模型,利用随机梯度下降法对无线传感器网络中图像的重构问题进行求解。深度压缩感知模型与传统方法相比,可以更有效地构建出图像的低维表示。实验结果表明,本发明提出的深度压缩感知图像重构方法与现有方法相比,可以有效提升无线传感器网络中农作物图像的重构效果。

Description

一种应用于无线传感器网络的深度压缩感知图像重构方法
【技术领域】本发明公开了一种应用于无线传感器网络的深度压缩感知图像重构方法,属于信号处理技术领域。
【背景技术】无线传感器网络(Wireless Sensor Network,WSN)以其组网自由、功耗低、鲁棒性强的优势,被广泛应用于辅助农业生产、地质灾害检测等领域。但是由于功耗的限制,无线传感器网络的图像采集和传输一直是一个难题。压缩感知(CompressiveSensing,CS)技术将稀疏信号的采样与压缩相结合,能够极大降低测量系统的采样率,突破奈奎斯特采样定律的限制。这使得压缩感知在无线传感器网络、核磁共振成像等领域有广泛应用。
压缩感知解决的问题是从m个观测中恢复未知
Figure BSA0000211084160000011
y=Ax+η
其中
Figure BSA0000211084160000012
是测量矩阵(m<<n),η是噪声。这是一个欠定方程组,除非x具有某种特性,否则无法得到唯一解。
幸运的是,现实中的信号都是稀疏的,或通过某种变换可以变成稀疏的。农业物联网中的图像本身不具有稀疏性,但是通过DCT或小波变换,能够成为稀疏信号。利用信号x的稀疏性,压缩感知问题的求解成为可能。
传统的压缩感知重构算法主要有贪婪算法、凸松弛算法和贝叶斯框架三类,但是这些算法在无线传感器网络图像压缩感知方面的表现不能令人满意。而卷积神经网络(Convolutional Neural Networks,CNN)在图像处理方面有非常优秀的表现。深度卷积生成网络(Deep Convolutional Generative Adversarial Networks,DCGAN)更是为图像的低维表示提供了新的思路。
对抗生成网络(Generative Adversarial Networks,GAN)是一种无监督的学习方法,它巧妙的利用了“对抗”的概念来完成生成器的学习,训练完成后可以生成与训练集相似的全新数据样本。DCGAN将GAN的概念引入到CNN中,能够生成高质量的图片样本。DCGAN的原理非常简单,它包括两个网络:生成网络G(Generator)和判别网络D(Discriminator)。G负责从随机的噪声中产生图片,D负责判断一张图片是不是真实的。在训练过程中,G的目的是产生尽可能骗过D的图片,D的目的是尽可能将生成图片与真实图片区分开。通过不断的学习,二者达到纳什均衡。G能够产生足够真实的图片。
引入深度生成网络的压缩感知模型被称为深度压缩感知(Deep CompressiveSensing,DCS)。经过充分训练的DCGAN网络的生成器能够以低维随机信号作为输入,生成高维图像,这与压缩感知的思想不谋而合。相较于传统的压缩感知图像重构方法,DCS模型能够在较低采样率的情况下较高质量的恢复植物幼苗图像,大大降低系统的采样率。这使得无线传感器网络传输图像信息的功耗大大降低。
【发明内容】本发明的目的在于将深度压缩感知模型应用于无线传感器网络中的农作物图像重建过程,提出一种应用于无线传感器网络的深度压缩感知图像重构方法。
本发明涉及方法的目标是通过观测值y=Ax+η和测量矩阵
Figure BSA0000211084160000021
寻找
Figure BSA00002110841600000214
接近x。
假设所有可能的农作物图像构成空间
Figure BSA0000211084160000022
由于X是农业物联网中特定种类的广义稀疏图像构成的集合,因此对于特定的x∈X,存在低维表示
Figure BSA0000211084160000023
所有可能的低维表示构成低维空间
Figure BSA0000211084160000024
由DCGAN网络训练得到的从低维空间Z到高维空间X的一个映射,即:
Figure BSA0000211084160000025
因此,求解问题变为寻找
Figure BSA00002110841600000215
使得
Figure BSA0000211084160000026
而z的维度通常是低于y的,这使得问题的求解变得简单。求解问题可以转化为一个非线性最小二乘求解问题:
Figure BSA0000211084160000027
这是一个凸优化问题,可以利用梯度下降法进行求解。实际应用时,添加正则项L(z)可以达到更好的恢复效果。即:
Figure BSA0000211084160000028
使用ADAM优化器进行问题上述问题进行求解。损失函数为:
Figure BSA0000211084160000029
梯度下降本质上是一种贪心算法,容易陷入局部最优解。因此引入智能搜索策略寻找全局最优解。
本发明的目的是通过下述技术方案来实现无线传感器网络中图像深度压缩感知问题的求解:
(1)定义:图像尺寸为(S,S,3),
Figure BSA00002110841600000210
为原图,其中n=S×S×3,测量矩阵为
Figure BSA00002110841600000211
Figure BSA00002110841600000212
观测值y=Ax+η,其中η为噪声,图像的低维表示为
Figure BSA00002110841600000213
(2)深度卷积生成网络模型构建:
生成器:输入为
Figure BSA0000211084160000031
经过全连接层投影并重塑为形状为(S/16,S/16,512)的张量,依次经过四个卷积核尺寸为5×5,步长为2的反卷积层,各层卷积核数量分别为256、128、64和3。除输出层外所有层使用批量归一化和Relu激活函数。输出层使用tanh激活函数;
判别器:输入为训练集图像或生成器产生的图像,依次经过三个卷积层和一个全连接层后输出一个判别值。卷积层的卷积核为5×5,步长为2,卷积核数量分别为64、128和256。除输出层外使用LeakyRelu激活函数、批量归一化并随机丢弃30%的神经元。
(3)使用交叉熵作为损失函数、Adam作为优化器训练深度卷积生成网络,最终的生成器为
Figure BSA0000211084160000032
(4)深度压缩感知重建图像:
①输入:观测值
Figure BSA0000211084160000033
测量矩阵
Figure BSA0000211084160000034
②初始化:最大搜索次数R,每次搜索最大迭代次数epoches,提前终止迭代阈值T1,重构成功判断阈值T2,学习率lr,正则项系数α、λ,早期判断阈值Ely;
③初始化z为服从高斯分布的随机向量,迭代次数置0;
④计算损失函数
Figure BSA0000211084160000035
⑤更新
Figure BSA0000211084160000036
⑥计算欧几里得距离
Figure BSA0000211084160000037
⑦如果迭代次数大于Ely且D>T1,则返回③;
⑧如果迭代次数小于epoches,则返回④;
⑨如果搜索次数小于R且D<T2,返回③;
⑩重建图像
Figure BSA0000211084160000038
【本发明的优点和积极效果】与现有技术相比,本发明具有如下优点和积极效果:
第一,本发明将深度压缩感知应用于无线传感器网络农作物图像的获取,利用农业物联网中目标农作物图像特征单一的特点,预先在植物图像数据集上学习该类图像的低维表示。将生成模型作为先验信息应用于农作物图像重构过程,显著降低了重构误差,大大减小了采样率。降低无线传感器网络图像传输的功耗,提高农作物图像的获取效率;
第二,针对随机梯度下降法容易陷入局部最优解的问题,引入智能搜索策略优化深度压缩感知重构模型,进一步提高重构成功率,降低重构误差,减少得到最优解消耗的时间。通过将提出方法与传统压缩感知重构方法进行比较,本发明效果大大优于传统方法;
第三,无线传感器网络实际应用于农业物联网时,采集到的农作物图像往往存在噪声。测试了本发明对于高斯噪声的抵抗能力后,实验结果表明,对于植物幼苗图像,本发明相较于传统方法,在采样率较低的情况下有着更强的抵抗噪声能力。
【附图说明】图1是本发明提出的一种应用于无线传感器网络的深度压缩感知图像重构方法流程图;
图2是本发明与Lasso算法在不同采样率下的重构误差比较;
图3是本发明和Lasso算法在采样率为0.1的情况下对农作物幼苗图像的重构样例;
图4是本发明和Lasso算法在采样率为0.1的情况下对存在高斯噪声农作物幼苗图像额重构样例。
【具体实施方式】为使本发明的实施方案与意义优势表述得更为清楚,下面结合后文附图及实施样例,对本发明进行更为详细的说明。
(1)定义:图像尺寸为(S,S,3),
Figure BSA0000211084160000041
为原图,其中n=S×S×3,测量矩阵为
Figure BSA0000211084160000042
Figure BSA0000211084160000043
观测值y=Ax+η,其中η为噪声,图像的低维表示为
Figure BSA0000211084160000044
(2)深度卷积生成网络模型构建:
生成器:输入为
Figure BSA0000211084160000045
经过全连接层投影并重塑为形状为(S/16,S/16,512)的张量,依次经过四个卷积核尺寸为5×5,步长为2的反卷积层,各层卷积核数量分别为256、128、64和3。除输出层外所有层使用批量归一化和Relu激活函数。输出层使用tanh激活函数;
判别器:输入为训练集图像或生成器产生的图像,依次经过三个卷积层和一个全连接层后输出一个判别值。卷积层的卷积核为5×5,步长为2,卷积核数量分别为64、128和256。除输出层外使用LeakyRelu激活函数、批量归一化并随机丢弃30%的神经元。
(3)使用交叉熵作为损失函数、Adam作为优化器训练深度卷积生成网络,最终的生成器为
Figure BSA0000211084160000051
(4)深度压缩感知重建图像:
①输入:观测值
Figure BSA0000211084160000052
测量矩阵
Figure BSA0000211084160000053
②初始化:最大搜索次数R,每次搜索最大迭代次数epoches,提前终止迭代阈值T1,重构成功判断阈值T2,学习率lr,正则项系数α、λ,早期判断阈值Ely;
③初始化z为服从高斯分布的随机向量,迭代次数置0;
④计算损失函数
Figure BSA0000211084160000054
⑤更新
Figure BSA0000211084160000055
⑥计算欧几里得距离
Figure BSA0000211084160000056
⑦如果迭代次数大于Ely且D>T1,则返回③;
⑧如果迭代次数小于epoches,则返回④;
⑨如果搜索次数小于R且D<T2,返回③;
⑩重建图像
Figure BSA0000211084160000057
本发明的仿真实验硬件配置为:Intel(R)Xeon(R)W-2145@3.70GHz CPU;64.0GBDDR4内存;NVIDIA Quadro RTX4000 GPU。
本发明的仿真实验软件配置为:Windows10 64位操作系统,仿真语言为Python,软件库为Tensorflow2.0。
在仿真实验中,使用的数据集为去除背景的Plant seedling dataset图片集。该图片集包括超过5000张植物幼苗图片,实验中,将100张图片作为测试集,其余图片用于DCGAN网络的训练。对比实验将本发明提出的一种应用于无线传感器网络的深度压缩感知图像重构方法与Lasso算法进行对比。采用每像素均方误差(Mean Squared Error,MSE)作为重构误差。噪声测试实验中噪声的添加方式为:
noisy(x)=x+σg
其中g的每个元素都服从标准正态分布。σ是噪声比例,即噪声与图片RGB数值的比值,实验中设置为0.2。
实验中将采样率分别设定为0.003、0.005、0.007、0.01、0.04、0.07、0.1、0.3、0.5,本发明的参数选择具有一定的普适性,本实验的所有图像均使用同一组参数。重构误差对比图中的竖线为95%置信区间。可以看出,本发明方法在采样率较低相比于传统算法能够达到更好的重构效果,并且对于图像噪声具有更强的地勘更努力。更具体地,在采样率为为0.01的情况下,本发明的重构误差仅为0.013,远优于Lasso算法的0.034。同样,注意到当采样率较高时,本发明重构误差高于Lasso算法,但是对于诸如无线传感器网络此类需要低采样率的应用场合,本发明所提方法更加适用。

Claims (1)

1.一种应用于无线传感器网络的深度压缩感知图像重构方法,包括如下内容
(1)定义:图像尺寸为(S,S,3),
Figure FSA0000211084150000011
为原图,其中n=S×S×3,测量矩阵为
Figure FSA0000211084150000012
Figure FSA0000211084150000013
观测值y=Ax+η,其中η为噪声,图像的低维表示为
Figure DEST_PATH_IMAGE002
(2)深度卷积生成网络模型构建:
生成器:输入为
Figure FSA0000211084150000015
经过全连接层投影并重塑为形状为(S/16,S/16,512)的张量,依次经过四个卷积核尺寸为5×5,步长为2的反卷积层,各层卷积核数量分别为256、128、64和3, 除输出层外所有层使用批量归一化和Relu激活函数, 输出层使用tanh激活函数;
判别器:输入为训练集图像或生成器产生的图像,依次经过三个卷积层和一个全连接层后输出一个判别值,卷积层的卷积核为5×5,步长为2,卷积核数量分别为64、128和256,除输出层外使用LeakyRelu激活函数、批量归一化并随机丢弃30%的神经元;
(3)使用交叉熵作为损失函数、Adam作为优化器训练深度卷积生成网络,最终的生成器为
Figure DEST_PATH_IMAGE004
(4)深度压缩感知重建图像:
①输入:观测值
Figure FSA0000211084150000016
测量矩阵
Figure FSA0000211084150000017
②初始化:最大搜索次数R,每次搜索最大迭代次数epoches,提前终止迭代阈值T1,重构成功判断阈值T2,学习率lr,正则项系数α、λ,早期判断阈值Ely;
③初始化z为服从高斯分布的随机向量,迭代次数置0;
④计算损失函数
Figure FSA0000211084150000018
⑤更新
Figure FSA0000211084150000019
⑥计算欧几里得距离
Figure FSA00002110841500000110
⑦如果迭代次数大于Ely且D>T1,则返回③;
⑧如果迭代次数小于epoches,则返回④;
⑨如果搜索次数小于R且D<T2,返回③;
⑩重建图像
Figure FSA00002110841500000111
CN202010533495.XA 2020-06-16 2020-06-16 一种应用于无线传感器网络的深度压缩感知图像重构方法 Active CN111681156B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010533495.XA CN111681156B (zh) 2020-06-16 2020-06-16 一种应用于无线传感器网络的深度压缩感知图像重构方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010533495.XA CN111681156B (zh) 2020-06-16 2020-06-16 一种应用于无线传感器网络的深度压缩感知图像重构方法

Publications (2)

Publication Number Publication Date
CN111681156A CN111681156A (zh) 2020-09-18
CN111681156B true CN111681156B (zh) 2022-09-23

Family

ID=72454739

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010533495.XA Active CN111681156B (zh) 2020-06-16 2020-06-16 一种应用于无线传感器网络的深度压缩感知图像重构方法

Country Status (1)

Country Link
CN (1) CN111681156B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113052925A (zh) * 2021-04-02 2021-06-29 广东工业大学 一种基于深度学习的压缩感知重构方法和系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521879A (zh) * 2012-01-06 2012-06-27 肖华 一种2d转3d的方法
CN108510464A (zh) * 2018-01-30 2018-09-07 西安电子科技大学 基于分块观测的压缩感知网络及全图重构方法
CN109102461A (zh) * 2018-06-15 2018-12-28 深圳大学 低采样分块压缩感知的图像重构方法、装置、设备及介质

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108171762B (zh) * 2017-12-27 2021-10-12 河海大学常州校区 一种深度学习的压缩感知同类图像快速重构系统与方法
CN110361778B (zh) * 2019-07-04 2020-10-13 华北电力大学 一种基于生成对抗网络的地震数据重建方法
CN110751049B (zh) * 2019-09-20 2022-04-29 浙江工业大学 一种面向信号采样梯度攻击的防御方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521879A (zh) * 2012-01-06 2012-06-27 肖华 一种2d转3d的方法
CN108510464A (zh) * 2018-01-30 2018-09-07 西安电子科技大学 基于分块观测的压缩感知网络及全图重构方法
CN109102461A (zh) * 2018-06-15 2018-12-28 深圳大学 低采样分块压缩感知的图像重构方法、装置、设备及介质

Also Published As

Publication number Publication date
CN111681156A (zh) 2020-09-18

Similar Documents

Publication Publication Date Title
CN110135366B (zh) 基于多尺度生成对抗网络的遮挡行人重识别方法
CN110728192B (zh) 一种基于新型特征金字塔深度网络的高分遥感图像分类方法
CN110348399B (zh) 基于原型学习机制和多维残差网络的高光谱智能分类方法
CN109002848B (zh) 一种基于特征映射神经网络的弱小目标检测方法
CN109472837A (zh) 基于条件生成对抗网络的光电图像转换方法
Ye et al. CSformer: Bridging convolution and transformer for compressive sensing
CN113327231B (zh) 一种基于空谱联合的高光谱异常目标检测方法和系统
CN112614070B (zh) 一种基于DefogNet的单幅图像去雾方法
CN111932452B (zh) 基于可见光图像增强的红外图像卷积神经网络超分辨方法
CN114170088A (zh) 一种基于图结构数据的关系型强化学习系统及方法
CN113793261A (zh) 一种基于3d注意力机制全通道融合网络的光谱重建方法
Nakashima et al. Learning to drop points for lidar scan synthesis
CN111681156B (zh) 一种应用于无线传感器网络的深度压缩感知图像重构方法
Liu et al. Modal regression-based graph representation for noise robust face hallucination
CN111798531B (zh) 一种应用于植株监测的图像深度卷积压缩感知重构方法
CN112184552B (zh) 一种基于高频特征学习的子像素卷积图像超分辨方法
CN113705340A (zh) 一种基于雷达遥感数据的深度学习变化检测方法
CN117036901A (zh) 一种基于视觉自注意力模型的小样本微调方法
CN110728728A (zh) 一种基于非局部正则的压缩感知网络图像重建方法
CN107330912B (zh) 一种基于多特征融合的稀疏表示的目标追踪方法
CN116385281A (zh) 一种基于真实噪声模型与生成对抗网络的遥感图像去噪方法
CN115760589A (zh) 一种用于运动模糊图像的图像优化方法及装置
CN115439669A (zh) 基于深度学习的特征点检测网络及跨分辨率图像匹配方法
Ramasinghe et al. A learnable radial basis positional embedding for coordinate-mlps
CN110503061B (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