本发明涉及计算机技术与无线通信技术领域，尤其涉及一种基于压缩感知的无线图像传输方法。
背景技术：
在目前的无线环境下主要是以传统的数字传输方式进行图像传输。然而在实际的无线环境下，信道状态是是时变的，且对于广/多播的情况下，用户具有不同的信道状态。
传统的数字传输方式面临信道的时变或多样性会出现悬崖效应，影响用户的质量体验。近年来出现的去除了数字传输中的量化、熵编码和信道编码的传输方式，其通过去相关变换和功率失真优化可以实现用户重建质量随信道状态平滑衰落，避免了悬崖效应。这样一种传输方式，被称之为伪模拟传输方式，其更适合图像、视频无线传输。
然而，无线环境下经常面临带宽不足的问题，传统的数字方式可以通过改变量化参数来实现最优的带宽匹配。对于伪模拟传输，由于去除了量化和熵编码，如何实现最优的带宽匹配是一个开放的问题。目前的一种常用的简单策略是直接丢弃能量值低(重要性低)的系数块来实现带宽匹配，但是这并不是最优的带宽匹配策略。
技术实现要素：
本发明的目的是提供一种基于压缩感知的无线图像传输方法，可以提高带宽不足的无线环境下图像传输的性能。
本发明的目的是通过以下技术方案实现的：
一种基于压缩感知的无线图像传输方法，包括：
将图像进行小波变换，并将得到的系数均匀分块，获得一系列的系数块；
对于其中的低频系数块进行幅值偏移操作；
对于其中的高频系数块，结合高频系数块能量值，利用压缩感知来进行高频系数块的带宽压缩并实现带宽匹配；
将进行幅值偏移操作后的低频系数块与经过压缩感知处理后的高频系数块进行功率分配后向外传输。
由上述本发明提供的技术方案可以看出，根据系数块在变换域系数稀疏的特点，压缩感知的引入可以成功的对系数块进行压缩，从而可以传输更多的系数块。另外通过推导的针对压缩感知系数块的功率分配因子实现了最优的功率失真优化，相比于传统的丢弃策略，此方法可以实现更好的传输性能。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他附图。
图1为本发明实施例提供的一种基于压缩感知的无线图像传输方法的流程图。
具体实施方式
下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明的保护范围。
本发明实施例提供一种基于压缩感知的无线图像传输方法，在无线环境带宽不足的情况下，尽可能多的传输系数块，从而传输更多的信息量。本发明实施例中，图像传输是以系数块为单位，不同系数块具有不同的系数分布，即系数块能量值不同，稀疏度不同。为了实现在带宽匹配下传输更多的信息量，根据系数块能量值并利用压缩感知来进行系数块带宽压缩并实现带宽匹配。在进行带宽匹配后，我们推导出最优的基于压缩感知的伪模拟传输功率失真优化模型，确定系数块的功率分配因子。
如图1所示，本发明实施例提供的基于压缩感知的无线图像传输方法主要包括如下步骤：
1、将图像进行小波变换，并将得到的系数均匀分块，获得一系列的系数块。
本发明实施例中，为了减小计算复杂度，所有操作都是以系数块为单位进行。
本发明实施例中，系数块分为低频系数块与高频系数块，在对一张图片进行去相关变换之后，左上角的系数块属于低频系数块，包含有大部分的能量；右下角的系数属于高频部分，能量较少，而且是稀疏的，便于使用压缩感知。
2、对于其中的低频系数块进行幅值偏移操作。
本发明实施例中，对于小波变换后的低频系数，根据其能量值很高且相邻系数值相近的特性，可以对低频系数块进行幅值偏移操作并直接进行伪模拟传输(不采用压缩感知操作)，这样可以提高功率效率。
幅值偏移操作具体过程如下：将低频系数块划分为更小的系数块并计算划分后每一系数块的均值，再将划分后每一系数块内所有系数减去所对应的均值，得到幅值偏移后的系数块。
3、对于其中的高频系数块，结合高频系数块能量值，利用压缩感知来进行高频系数块的带宽压缩并实现带宽匹配。
本发明实施例中，对于高频系数块，根据高频系数块能量值权重进行带宽初分配，计算高频系数块的权重因子wi：
其中，λi表示高频系数块i的能量值；∑λj表示对所有高频系数块能量求和；
结合权重因子wi计算高频系数块i的初分配的带宽(采样率)bi：
bi＝b·wi；
其中，b为低频系数块占用部分带宽后剩余带宽。
本发明实施例中，假设每个高频系数块(所有待进行压缩感知的高频系数块，对应的它们的压缩比)所分配带宽为特定数值，特定数值可以设为1。因此每个高频系数块所分配的带宽等于其所对应的采样率。当bi大于1时，将其设置为1，然后将带宽小于1的高频系数块重新按照上述公式进行带宽计算。重复这个步骤直到所有高频系数块的带宽均不大于1。此处，对超过1的高频系数块置1，然后将剩余的带宽分配给小于1的高频系数块；这样做的目的是，一方面，保证对于所有的bi，有0<bi<＝1；另一方面，因为给定的带宽是固定的，将带宽大于1的系数块设定为1之后就会有部分带宽多余出来，从而将这部分多余出来的带宽再分配给带宽小于1的高频系数块。
筛选出带宽初分配结果小于特定数值的高频系数块，并通过下述方式进行带宽重分配：计算筛选出的每一高频系数块的稀疏度η，当η小于第一阈值时，即丢弃相应高频系数块；当η小于第二阈值时，相应高频系数块的带宽分配为a1；当η小第三阈值时，相应高频系数块的带宽分配为a2；当η小于第四阈值时，相应高频系数块的带宽分配为a3；当η大于第四阈值时，相应高频系数块的带宽分配为a4；其中，a4>a3>a2>a1；将筛选出的高频系数块按照能量值大小降序排列，并基于上述稀疏度与带宽的关系依次进行带宽分配，直至带宽用尽，并丢弃未分配到带宽的高频系数块。
本发明实施例中，稀疏度是指系数幅值大于预设值的个数占总个数的比值。
示例性的，预设值可以设为5；第一、第二、第三、第四阈值可对应的设为0.05、0.1、0.2、0.4。a1、a2、a3、a4可对应的设为0.45、0.6、0.75、0.9。需要说明的是，这些参数的具体数值仅为举例并非构成限制，在实际应用中，用户可以根据实际情况或者经验来设定各个参数的具体数值。
4、将进行幅值偏移操作后的低频系数块与经过压缩感知处理后的高频系数块进行功率分配后向外传输。
本发明实施例中，功率分配时，通过下式计算需要传输的各系数块(即进行幅值偏移操作后的低频系数块，以及压缩感知处理后的高频系数块)的功率分配因子：
其中，gk为当前所计算的系数块k的功率分配因子，系数块k可以是幅值偏移操作后的低频系数块，也可以是压缩感知处理后的高频系数块；p是总功率；分子、分母部分的ck、cu是指系数块k、系数块u的功率分配参数，系数块k与系数块u实际含义类似，只是在公式中系数块k对应为当前所计算的系数块，即可以理解为一个确定值，而u则对应于所有系数块中任一个，结合相关求和符号即表示针对所有系数块相关参数求和。
此外，假设每一系数块的大小为n×n，当系数块k或者系数块u，是低频系数块或者带宽为特定数值(如前文所述，特定数值可以设为1)的高频系数块时，功率分配参数为n2；当系数块k或者系数块u，是带宽小于特定数值的高频系数块时，结合其进行压缩感知处理后所分配的带宽大小m(即之前a1、a2、a3或a4)来计算功率分配参数ck与cu的计算公式表示为：ck与
类似的，λk'、λu'对应为系数块k、系数块u的能量值；当系数块k或者系数块u，是低频系数块或者带宽为特定数值的高频系数块时，能量值未发生变化，依然为原能量值；当系数块k或者系数块u，是带宽小于特定数值的高频系数块时，代表经过压缩感知处理之后的能量值。(由于经过压缩感知处理，此部分系数块系数放生变化，其能量值也随之变化)。
此外，向外传输时，功率分配后的各系数块通过qam调制后直接在无线信道上传输，所涉及的各系数块的能量值以及低频系数块进行幅值偏移操作所计算到的所有均值均作为元数据通过数字传输方式传输。
本发明实施例上述方案主要具有如下优点：根据系数块在变换域系数稀疏的特点，压缩感知的引入可以成功的对系数块进行压缩，从而可以传输更多的系数块。另外通过推导的针对压缩感知系数块的功率分配因子实现了最优的功率失真优化，相比于传统的丢弃策略，此方法可以实现更好的传输性能。
为了便于理解，下面结合示例做进行介绍。
本示例中，基于压缩感知的无线图像传输方法主要通过低频部分幅值偏移，高频部分带宽分配并压缩感知以及整体的功率分配组成。主要分为以下几个步骤：
1)图像处理。
将图像进行9/7小波变换处理(三级小波分解)，得到的系数被均匀划分成n×n的系数块(假设n＝32)。在本示例中并不限定图像的分辨率。
2)幅值偏移操作。
对划分后的低频系数块进行幅值偏移操作，我们将系数块进一步划分为更小的m×m系数块(假设m＝8)并计算每个系数块的均值。系数块中每个系数减去均值得到系数值更小的幅值偏移后的系数。本示例中划分后系数块的均值将被作为元数据。
3)压缩感知处理。
对于低频系数块，其带宽直接设置为1，即不进行压缩感知处理。
高频系数块根据前述操作进行带宽分配处理，得到每个系数块所对应的带宽值。这里计算稀疏度时所涉及的预设值设置为5，然后根据压缩感知操作进行线性映射。本示例中，选取高斯随机矩阵作为压缩感知的观测矩阵，矩阵大小为1024×1024。对于带宽为的bi系数块，将其排列成列并与观测矩阵的前l(l＝bi·n2)行进行矩阵乘积操作得到需要传输的测量值系数(即进行压缩感知处理之后得到的测量值系数)。
4)功率分配
计算需要传输的系数块的能量值(压缩感知系数块则计算线性映射后的系数块能量值)并根据前述功率因子公式计算每个系数块的功率分配因子。本示例中的系数块能量值作为元数据。
本领域技术人员可以理解，压缩感知之后能量值发生了变化，则重新计算高频系数块的能量值，对系数块中的系数值平方进行求和。
5)传输与重建
经过上述处理的系数块将通过qam调制直接在无线信道上传输，元数据因为其需要准确的传输将通过可靠的数字传输方式低速率传输。
在接收端接收到的系数块将通过最小均方误差重建(mmse)得到接收系数块，然后对经过压缩感知处理的系数块进行重建(nesta重建算法)。最后将得到的高频系数和低频系数进行合成，形成完整的小波系数。通过逆小波变换重建出接收的图像。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例可以通过软件实现，也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解，上述实施例的技术方案可以以软件产品的形式体现出来，该软件产品可以存储在一个非易失性存储介质(可以是cd-rom，u盘，移动硬盘等)中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。