基于感兴趣区域的海岸带高光谱图像分布式有损编码和解码方法
【专利摘要】本发明公开一种基于感兴趣区域的海岸带高光谱图像分布式有损编码和解码方法，属于图像处理领域，该方法依据高光谱图像的谱间相关性将谱段帧分成若干谱段组，每个谱段组包含1个关键帧和若干个WZ帧；对每个谱段组的关键帧进行感兴趣区域提取后，对关键帧的感兴趣区域进行JPEG?LS无损编解码，对其背景区域进行基于SPIHT的编解码，而对每个谱段组的WZ帧则采用基于压缩感知的编解码。本发明有效地实现了编码复杂度低、较强的抗误码性以及码流具有质量可分级等的特性。
【专利说明】
基于感兴趣区域的海岸带高光谱图像分布式有损编码和解码 方法
技术领域
[0001] 本发明设及图像处理领域，尤其是一种编码效率高、抗误码性能好、质量可分级、 编码复杂度低的基于感兴趣区域的海岸带高光谱图像分布式有损编码和解码方法。
【背景技术】
[0002] 近年来，随着环勸海等多个沿海重要经济区发展规划上升为国家战略，基于高光 谱遥感的海岸带环境监测与管理的需求也在不断增加。然而，越来越高的光谱分辨率和空 间分辨率导致其信息量也飞速增加。同时，用于星上传输的信道的恶劣环境给高光谱遥感 图像的压缩和传输带来了巨大的困难。在运种情况下，如何有效地压缩和传输海岸带的高 光谱遥感图像信息就成为了当前高光谱海岸带图像应用的前沿问题之一。
[0003] 英国萨里空间中屯、和萨里卫星技术有限公司的研究人员对遥感图像编码算法进 行总结后指出，实际应用中目前最为流行的星载数据压缩编码算法主要包括肝EG算法(如 AL0S、北京一号小卫星、Cartosat-l/2、Pleiades-HR、Proba-2、Rapi 祀ye 和THE0S等）、 肝EG2000算法(如IMS-1、RASAT和X-SAT等)和CCSDS-LDC无损压缩算法(如PICA畑等）。然而， 上述方法都只是仅仅通过去除一个谱段的空间相关性实现压缩，其性能远不能满足高光谱 遥感技术的应用需求，至今为止尚未真正形成一种适应星载遥感图像处理的低复杂度要 求、强鲁棒性信息传输和高编码效率的高光谱图像压缩实用算法和标准。
[0004] -方面，2003年后出现的分布式信源编码架构因其编码端具有较低的编码复杂 度、高效的压缩性能和较高的抗误码性等优点成为针对星上图像编码的重要工具。但是，基 于分布式信源编码的高光谱遥感图像编码方案还不是很多，特别是面向传感器网络、遥感 系统、克服有限带宽和能量限制的分布式编码方法还处于不断发展阶段，尚未有被广泛认 可的、成熟稳定的算法出现。
[0005] 另一方面，大部分观察者在观看图像时，都只会对图像中感兴趣区域的信息进行 重点分析，而不是对全局信息进行分析。于是，研究人员逐渐达成了一个共识：由于图像的 图像感兴趣区域(Region of Interest, R0I)包含了整幅图像的主要信息和关键信息，R0I 可W近似代替原始图像。运样，在处理图像时选择重点对R0I进行处理和分析，可W大幅度 地降低计算量，从而能够更好地应对高光谱图像数据量大幅增长的问题。
【发明内容】
[0006] 本发明是为了解决现有技术所存在的上述技术问题，提供一种编码效率高、抗误 码性能好、质量可分级、编码复杂度低的基于感兴趣区域的海岸带高光谱图像分布式有损 编码和解码方法。
[0007] 本发明的技术解决方案是:一种基于感兴趣区域的海岸带高光谱图像分布式有损 编码方法，其特征在于按照如下步骤进行： 步骤1.输入待编码的海岸带高光谱图像I，将谱段帖的高度N和宽度Μ输出至压缩码流， 每个值固定占2个字节； 步骤2.按公式(1.)计算I的每一个谱段k与其相邻的两个较低谱段裘一I、套-逸的谱间 相关系数，若相关系数均小于0.9，则将谱段k设置为关键帖X，否则，将其设置为WZ帖Y;
纖 式（1)中，玄辨法游和?辩:补嵌轉分别为第k谱段和第教，補.谱段图像的i行、j列像素的 灰度值；餐和％^则分别为第k和采剌谱段所有像素点灰度值的平均值，运里r氏玄I；，此 夕h第1个谱段固定作为关键帖;将所选关键帖的数目G置于1个字节中，将关键帖的谱段号 置于其后面的G个字节中，然后将运:|資半||个字节输出至压缩码流； 步骤3.将每个关键帖作为一个谱段组，然后将每个WZ帖划分到与之最近的、且谱段号 小于该WZ帖谱段号的关键帖所在的谱段组，从而形成G个谱段组;令g二1; 步骤4.若客< G，则转入步骤5;否则，转入步骤10; 步骤5.对第g个谱段组的关键帖进行3级最大提升格形态小波变换，在最低分辨率下的高 频子带好抒中利用化nny算子进行边缘检测，并确定包含所有边缘点的最小包围盒，将与该 包围盒对应的空间域像素集合作为感兴趣区域*3汝，而其余区域作为背景^3，然后将包围 盒左上角和右下角的坐标写入压缩码流，每个坐标分量固定占用pQ霉 字节； 步骤6.对第g个谱段组关键帖的^3^进行肝EG-LS无损编码，对第g个谱段组关键帖的 -Yg采用可分级的编码方法SPIHT，二者的编码码流输出至压缩码流； 步骤7.对于第g个谱段组，将关键帖作为参考帖，为该组中的每一个WZ帖构造基于感兴 趣区域的1阶线性预测器，利用最小二乘法计算使预测均方差取得最小值的预测参数垃和 /3,并将其作为边信息经算术编码后输出到压缩码流，用W指导和矫正WZ帖的传输和解码： .?二輝'龄东激'乎洛)' 錢I 式(2)中，为与單?€>/相对应的WZ帖感兴趣区域，所述^为预测误差； 步骤8.对预测误差ζ进行典型的压缩感知采样，观测矩阵采用高斯随机矩阵，并将 采样值经算术编码后输出至压缩码流； 步骤9.令各二置手1，转入步骤4; 步骤10.将压缩码流通过下行信道传输至解码端，算法结束。
[0008] -种与上述基于感兴趣区域的海岸带高光谱图像分布式有损编码方法对应的解 码方法，其特征在于按照如下步骤进行： 步骤1.读入压缩码流的前4个字节，解码出谱段帖的高度N和宽度Μ; 步骤2.读入压缩码流的第5个字节，解码出关键帖的数目G; 步骤3.从压缩码流中继续读入G个字节，解码出关键帖的谱段号码； 步骤4.令器二1; 步骤5.若复< G，则转入步骤6;否则，算法结束； 步骤6.从压缩码流中继续读入
t字节，解码出第Ji:个谱 段组的感兴趣区域坐标； 步骤7.采用JPEG-LS无损解码重构出第g个谱段组的关键帖的感兴趣区域施/，采用 SPIHT解码重构出该关键帖的非感兴趣区域S ; 步骤8.用算术解码从压缩码流中提取出该谱段组的每个WZ帖的预测参数泣和舞，并 利用关键帖的-瑪9,讯1阶线性预测去胃夺获得第g个谱段组中每个WZ帖的预测； 步骤9.用算术解码从压缩码流中提取出每个WZ帖的预测误差的压缩采样值，然后利用 协同稀疏的压缩感知重构算法计算得到解码后的预测误差7。，并通过公式巧)得到每个WZ 帖的感兴趣区域的解码结果
鑛 步骤10.将该谱段组的关键帖非感兴趣区域作为每个WZ帖的非感兴趣区域，并将 每个WZ帖的与相加，获得每个WZ帖的解码结果； 步骤1].令复二复十至，转入步骤5。
[0009] 本发明与现有技术相比，具有W下优点：第一，采用分布式编码作为总体的编码框 架，并结合压缩感知技术，其编码复杂度低，对噪声的鲁棒性高，适合星载高光谱图像的编 码;第二，对感兴趣区域采用肝EG-LS无损编码方法可确保感兴趣区域的保真度W及解码图 像的人眼主观质量，对背景区域采用SPIHT编码方法使所获得的码流具有质量可分级特性； 第立，利用高光谱图像的谱间相关性对WZ帖进行线性预测，由于预测误差的稀疏性较高，采 用压缩感知编码在降低复杂度的同时，有利于提高图像感兴趣区域的重构质量。
【附图说明】
[0010] 图1是本发明实施例的流程图。
[0011] 图2是化/No=l.5地时，不同编码方法在不同码率下的解码图像对比图。
[0012] 图3是化/No=l. 5地、码率为0.4时，解码图像的放大结果对比图。 具体实施例
[0013]如图1所示:本发明实施例所提供的基于感兴趣区域的海岸带高光谱图像分布式 有损编码方法包括如下步骤： 步骤1.输入待编码的海岸带高光谱图像I，将谱段帖的高度N和宽度Μ输出至压缩码流， 每个值固定占2个字节； 步骤2.按公式α咐算I的每一个谱段k与其相邻的两个较低谱段冕一:|、煮交的谱间相 关系数，若相关系数均小于0.9，则将谱段k设置为关键帖X，否则，将其设置为WZ帖X;
a 式（1)中，巧策书凝和玄終+糾,巧分别为第k谱段和第患幸?.谱段图像的i行、j列像素的 灰度值;馬巧％:.则分别为第k和装给谱段所有像素点灰度值的平均值，运里# 2:1，此 夕h第1个谱段固定作为关键帖;将所选关键帖的数目G置于1个字节中，将关键帖的谱段号 置于其后面的G个字节中，然后将运I窝辛||个字节输出至压缩码流； 步骤3.将每个关键帖作为一个谱段组，然后将每个WZ帖划分到与之最近的、且谱段号 小于该WZ帖谱段号的关键帖所在的谱段组，从而形成G个谱段组;令咨二-i; 步骤4.若复-< G，则转入步骤5;否则，转入步骤10; 步骤5.对第g个谱段组的关键帖进行3级最大提升格形态小波变换，在最低分辨率下的高 频子带过中利用Canny算子进行边缘检测，并确定包含所有边缘点的最小包围盒，将与该包 围盒对应的空间域像素集合作为感兴趣区域^6^，而其余区域作为背景^，然后将包围 盒左上角和右下角的坐标写入压缩码流，每个坐标分量固定占用po豁??触文I黏i部! 字节； 步骤6.对第g个谱段组关键帖的?馈进行肝EG-LS无损编码，对第g个谱段组关键帖的 廉S采用可分级的编码方法SPIHT，二者的编码码流输出至压缩码流； 步骤7.对于第g个谱段组，将关键帖作为参考帖，为该组中的每一个WZ帖构造基于感兴 趣区域的1阶线性预测器，利用最小二乘法计算使预测均方差取得最小值的预测参数坊和 #，并将其作为边信息经算术编码后输出到压缩码流，用W指导和矫正WZ帖的传输和解码： .?.二玉:男級一.择迹纖.中媒 II) 其中，所述为与馬^^相对应的WZ帖感兴趣区域，所述为预测误差； 步骤8.对预测误差K进行典型的压缩感知采样，观测矩阵采用高斯随机矩阵，并将 采样值经算术编码后输出至压缩码流； 步骤9.令榮二祭去1，转入步骤4; 步骤10.将压缩码流通过下行信道传输至解码端，算法结束。
[0014] 如图1所示:本发明实施例所提供的基于感兴趣区域的海岸带高光谱图像分布式 有损解码方法包括如下步骤： 步骤1.读入压缩码流的前4个字节，解码出谱段帖的高度N和宽度M; 步骤2.读入压缩码流的第5个字节，解码出关键帖的数目G; 步骤3.从压缩码流中继续读入G个字节，解码出关键帖的谱段号码； 步骤4.令容二1; 步骤5.若含-< G，则转入步骤6;否则，算法结束； 步骤6.从压缩码流中继续读>
字节，解码出第g个谱 段组的感兴趣区域-坐标； 步骤7.采用JPEG-LS无损解码重构出第g个谱段组的关键帖的感兴趣区域-fsK，采用 SPIHT解码重构出该关键帖的非感兴趣区域jf's ; 步骤8.用算术解码从压缩码流中提取出该谱段组的每个WZ帖的预测参数货和y，并 利用关键帖的防和1阶线性预测（货JTg自J手获得第g个谱段组中每个WZ帖的预测； 步骤9.用算术解码从压缩码流中提取出每个WZ帖的预测误差的压缩采样值，然后利用 协同稀疏的压缩感知重构（Compressive Sensing Recovery via Collaborative Sparsity)算法计算得到解码后的预测误差1?，并通过公式得到每个WZ帖的感兴趣区 域的解码结果§^03:
鑛 步骤10.将该谱段组的关键帖非感兴趣区域?5作为每个WZ帖的非感兴趣区域，并将 每个WZ帖的f 厢加，获得每个WZ帖的解码结果； 步骤11.令置二咨^^，转入步骤5。
[0015] Eb/No=l. 5地时，不同编码方法在不同码率下的解码图像对比如图2所示。
[0016] Eb/No=l.5地、码率为0.4时，解码图像的放大结果对比如图3所示。
[0017] 恥/No=l. 5地时不同编码方法在不同码率下的解码图像的PSNR比较如表1。
[001引 表1
不同编码方法在不同码率下的运行时间比较如表2。 [0019]表2
【主权项】
1. 一种基于感兴趣区域的海岸带高光谱图像分布式有损编码方法，其特征在于按照如 下步骤进行： 步骤1.输入待编码的海岸带高光谱图像I，将谱段帧的高度N和宽度Μ输出至压缩码流， 每个值固定占2个字节； 步骤2.按公式_计算I的每一个谱段k与其相邻的两个较低谱段1一|、是一交的谱间相 关系数，若相关系数均小于〇. 9，则将谱段k设置为关键帧X，否则，将其设置为WZ帧Y;式⑴中，i構-祿和别为第k谱段和第谱段图像的?行、j列像素的 灰度值；尤和％4:则分别为第k和.的谱段所有像素点灰度值的平均值，这里|:名|:1;，此 外，第1个谱段固定作为关键帧;将所选关键帧的数目G置于1个字节中，将关键帧的谱段号 置于其后面的G个字节中，然后将这(G +1)个字节输出至压缩码流； 步骤3.将每个关键帧作为一个谱段组，然后将每个WZ帧划分到与之最近的、且谱段号 小于该WZ帧谱段号的关键帧所在的谱段组，从而形成G个谱段组;令g = 1; 步骤4.若g S G，则转入步骤5;否则，转入步骤10; 步骤5.对第g个谱段组的关键帧进行3级最大提升格形态小波变换，在最低分辨率下的高 频子带中利用Canny算子进行边缘检测，并确定包含所有边缘点的最小包围盒，将与该 包围盒对应的空间域像素集合作为感兴趣区域-％_，而其余区域作为背景@，然后将包围 盒左上角和右下角的坐标写入压缩码流，每个坐标分量固定占用 字~P ; 步骤6.对第g个谱段组关键帧的进行JPEG-LS无损编码，对第g个谱段组关键帧的 采用可分级的编码方法SPIHT，二者的编码码流输出至压缩码流； 步骤7.对于第g个谱段组，将关键帧作为参考帧，为该组中的每一个WZ帧构造基于感兴 趣区域的1阶线性预测器，利用最小二乘法计算使预测均方差取得最小值的预测参数心和 ?，并将其作为边信息经算术编码后输出到压缩码流，用以指导和矫正WZ帧的传输和解码：式(2)中，^沉为与相对应的WZ帧感兴趣区域，所述:?为预测误差； 步骤8.对预测误差_进行典型的压缩感知采样，观测矩阵采用高斯随机矩阵，并将采 样值经算术编码后输出至压缩码流； 步骤9.令:玄.=書+1.，转入步骤4; 步骤10.将压缩码流通过下行信道传输至解码端，算法结束。2. -种与权利要求1所述基于感兴趣区域的海岸带高光谱图像分布式有损编码方法对 应的解码方法，其特征在于按照如下步骤进行： 步骤1.读入压缩码流的前4个字节，解码出谱段帧的高度N和宽度M; 步骤2.读入压缩码流的第5个字节，解码出关键帧的数目G; 步骤3.从压缩码流中继续读入G个字节，解码出关键帧的谱段号码； 步骤4.令茗=1; 步骤5.若g S G，则转入步骤6;否则，算法结束； 步骤6.从压缩码流中继续读入|个字节，解码出第g个谱段 组的感兴趣区域坐标； 步骤7.采用JPEG-LS无损解码重构出第g个谱段组的关键帧的感兴趣区域^，采用 SPIHT解码重构出该关键帧的非感兴趣区域^ ; 步骤8.用算术解码从压缩码流中提取出该谱段组的每个WZ帧的预测参数祕和多，并利 用关键帧的和1阶线性预测+获得第g个谱段组中每个WZ帧的预测； 步骤9.用算术解码从压缩码流中提取出每个WZ帧的预测误差的压缩采样值，然后利用 协同稀疏的压缩感知重构算法计算得到解码后的预测误差，并通过公式微得到每个WZ 帧的感兴趣区域的解码结果_^^:步骤10.将该谱段组的关键帧非感兴趣区域^作为每个WZ帧的非感兴趣区域，并将 每个WZ帧的与相加，获得每个WZ帧的解码结果； 步骤11.令:倉，转入步骤5。
【文档编号】G06T9/00GK105825530SQ201610159167
【公开日】2016年8月3日
【申请日】2016年3月21日
【发明人】宋传鸣, 王相海, 张智迪, 李畅
【申请人】辽宁师范大学