多光谱遥感影像星上压缩影响评价.pdf
《多光谱遥感影像星上压缩影响评价.pdf》由会员分享,可在线阅读,更多相关《多光谱遥感影像星上压缩影响评价.pdf(7页珍藏版)》请在文库网上搜索。
1、2097-3012(2023)03-0328-07 Journal of Spatio-temporal Information 时空信息学报 收稿日期:2022-08-30;修订日期:2023-09-03 基金项目:战略性国际科技创新合作重点专项项目(2016YFE0205300)作者简介:林瑶瑶,研究方向为遥感图像处理与分析。E-mail: 通信作者:田永明,研究方向为遥感图像处理与分析。E-mail: 多光谱遥感影像星上压缩影响评价 林瑶瑶1,唐新明2,薛玉彩2,田永明1,齐麟3 1.济南市勘察测绘研究院,济南 250013;2.自然资源部国土卫星遥感应用中心,北京 100048;3.山
2、东高速工程检测有限公司,济南 250013 摘 要:开展卫星多光谱遥感影像星上压缩影响评价,是验证卫星影像是否满足应用指标的重要内容。高分七号卫星突破了我国 11 万立体测绘技术,可满足高精度住建应用、统计应用等。当前,仍缺少高分七号星上压缩指标论证的相关研究,因此,基于高分七号模拟多光谱影像,本文从影响影像灰度、纹理、相关性等构像质量特征和影像应用能力,包括影像分割一致性及影像分类精度两方面进行影像压缩评价,并以临朐县为例进行了验证。结果表明,对于星上 JPEG-LS 压缩算法提供的两种压缩比:压缩比为 31 时,影像灰度、纹理信息保持好;多光谱遥感影像分割后图斑总数变化率和平均图斑面积变化
3、率绝对值不超过 1%,影像分割一致性较高;影像分类后 Kappa 系数达 0.91。压缩比为 61 时,单波段 61 压缩后影像相关系数最小降低到 0.891,重建影像峰值信噪比小于 50,影像分类后 Kappa 系数为 0.87。多光谱影像宜采用 31 压缩比。本研究可为高分七号星上压缩指标论证提供参考。关键词:高分七号模拟影像;JPEG-LS 算法;影像特征;分割精度;分类精度 引用格式:林瑶瑶,唐新明,薛玉彩,田永明,齐麟.2023.多光谱遥感影像星上压缩影响评价.时空信息学报,30(3):328-334 Lin Y Y,Tang X M,Xue Y C,Tian Y M,Qi L.20
4、23.Evaluation of on-board compression effects on multispectral remote sensing images.Journal of Spatio-temporal Information,30(3):328-334,doi:10.20117/j.jsti.202303002 1 引 言 2012 年 1 月 9 日,资源三号卫星成功发射,它是我国首颗民用高分辨率光学立体测图卫星,主要应用于我国 15 万测绘产品制作,12.5 万地形图的修补与更替(王鸿燕等,2020;王煌文等,2022;刘翔等,2023)。高分七号测绘卫星的全色分辨率
5、达 0.7 m,多光谱影像分辨率也达到 2.8 m,是实现我国 11 万卫星测绘能力的关键卫星工程(赵文普等,2022;张太煜和冯献霆,2023)。随着图像分辨率的提高,卫星星上数据量亦是成倍增长,高分七号原始数据量已是资源三号的 5 倍。星上海量的遥感数据需要经过压缩后才能及时下传到地面站。在遥感领域,采取的压缩算法有基于小波变换的 JPEG2000 算法(连佳佳等,2019)、多级 树集合分裂排序(set partitioning in hierarchical trees,SPIHT)算法(叶辰等,2022)、基于离散余弦变换的 JPEG 压缩算法(刘婷和梁文龙,2022)、JPEG-L
6、S 压缩算法(尤传亮,2019)等。在星上进行影像压缩的实时性要求较高,同时星上缓存巨大数据量的能力也有限。为能够及时获得较高压缩倍率的数据,满足使用精度和要求,须采用硬件压缩方式来实现星载高分辨率遥感影像的缓存和下传。其中,JPEG-LS 压缩算法复杂度低,压缩性能良好,更适用于硬件系统(陈聪,2018;张人愉,2018)。遥感影像分类是多光谱影像应用的重要内容之一。经过压缩后的影像质量下降,部分影像信息丢失,影响到遥感影像在土地和地质应用、变化检测、专题信息提取与分析等方面的应用能力和应用范畴。Xue 等(2014)利用 JPEG-LS 压缩算法对资 林瑶瑶 等:多光谱遥感影像星上压缩影响
7、评价 329 源三号多光谱影像进行不同倍率压缩,分析影像融合效果,结果表明 3 倍压缩后的遥感影像可以满足测绘应用需求。Danoedoro(2013)比较分析了基于像素和基于对象的分类方法对 ALOS-AVNIR2 多光谱遥感影像分类精度的影响,结果表明基于对象的分类方法结果更稳定、分类精度更高。谭可清(2018)研究了 JPEG2000 压缩对 Landsat 8 与资源三号监督分类结果的影响。Elkholy 等(2019)对 Geo-eye 卫星遥感影像研究了离散小波变换(discrete wavelet transformation,DWT)压缩对 Geo-eye 卫星遥感影像的神经网络
8、及最大似然分类结果的影响。胡晔(2020)选择 JPEG2000 压缩算法对高分二号、Landsat 8 及资源三号卫星影像压缩,利用最大似然和支持向量机分类算法对影像进行分类,并计算 Kappa 系数、总体分类精度等精度指标。综上所述,目前对多光谱遥感影像压缩影响分析多集中在不同卫星影像受压缩算法影响程度,以及不同分类方法对多光谱影像应用能力影响等方面的研究,缺少对高分辨率遥感影像星上压缩质量影响的相关研究。因此,本文对高分七号模拟多光谱影像进行了系统的星上压缩影响评价实验。利用JPEG-LS算法对高分七号模拟多光谱影像进行不同倍率压缩,比较压缩后影像与未压缩影像灰度特征、纹理特征的差异;对
9、遥感影像进行多尺度分割,选取分割后的图斑单元作为样本,继而采用面向对象的分类方法对遥感影像进行分类。同时,以临朐县为例,将未压缩影像的分割和分类结果作参照,对压缩后影像进行分割一致性和分类精度的验证。2 多光谱影像压缩质量评价方法 2.1 构像质量评价 根据遥感影像构像质量分析内容、方向的不同,从灰度特征、纹理特征分析和影像对照差异分析两个方面进行构像质量分析评价(林瑶瑶,2019)。前者通过分别分析未压缩影像和重建影像的特征,研究不同压缩倍率后影像特征的变化;后者采取未压缩影像和重建影像直接进行关联的方法,研究压缩导致的影像质量的差别。影像特征分析从影像光谱特征和纹理特征两方面进行。其中,光
10、谱特征包括影像灰度平均值和灰度标准方差;纹理特征包括基于灰度共生矩阵提取的对比度和信息熵。影像光谱特征可以反映影像亮度值变化,影像纹理特征反映了纹理沟纹、影像信息量变化等。影像对比特征分析通过比较原始未压缩影像和重建影像差异来进行。相关系数与峰值信噪比是常用的两种评价因子。其中:11001111220000()()CC()()MNgfjiMNMNfgjijif i,jg i,jf i,jg i,j (1)式中,CC 为影像相关系数;(,)f i j、f分别为未压缩影像在第 i 行、j 列的灰度值及整幅影像的灰度平均值;(,)i jg、g分别为重建影像在第 i 行、j列的灰度值及整幅影像的灰度平
11、均值;M、N 分别为影像总行数、总列数。影像相关系数用于表示两幅影像灰度相似程度。另:102MSEPSNR10log(21)B (2)112001MSE(,)(,)MNjif i jg i jMN (3)式中,B为影像位数;PSNR为峰值信噪比,与均方误差(mean square error,MSE)相关,将压缩前后影像相同行、列号处的像元灰度差当作噪声,将未压缩影像当作信号,计算其信噪比。峰值信噪比大,影像受压缩影响小,影像质量高。2.2 分割一致性评价 本文利用基于区域生长的分割技术对多光谱影像数据进行分割。eCognition的多尺度分割方法将影像的空间及光谱异质性综合考虑,有效防止了只
12、顾及光谱异质性最小造成的分割后图斑形状较为零碎现象,在面向对象分类中广为应用(陈蕊和张继超,2020;林荣清等,2020)。本文采用控制变量的方法,通过设置分割尺度参数、形状因子及紧致度因子分析实验影像不同地物的分割效果,继而选取合适的分割参数对多光谱影像进行分割,实验中,各个波段权重设为1。结合遥感影像视觉分割效果和最佳分割尺度图,确定本文采用的分割尺度为28、形状因子为0.2、紧致度因子为0.5。采用确定的分割参数对原始未压缩影像和压缩后重建影像进行多尺度分割,以未压缩影像分割330 Journal of Spatio-temporal Information 时空信息学报 2023,30
13、(3)图斑数据为参考,分析不同压缩倍率对重建影像分割结果的影响,如图1所示。图 1 影像分割影响评价 Fig.1 Evaluation of image segmentation accuracy 2.3 分类精度评价 不同地物目标拥有不同的对象特征,关系到影像分析与信息提取;其中,光谱信息是最基础的特征,随着高分辨率影像的涌现,对象的形状和纹理信息得到充分表达。因此,本文运用遥感影像的光谱信息、形状信息和纹理信息进行多光谱影像的分类(潘俊虹和梁明,2022;王书珩,2022;武萌华,2023)。首先,利用监督分类方法对原始未压缩影像及压缩后重建影像分别进行分类,得到各自的分类结果;其次,将重
14、建影像分类结果与原始影像分类结果进行对比分析,研究压缩算法对多光谱影像专题信息分类精度及整体分类精度影响,评价路线如图2所示。图 2 影像分类精度评价 Fig.2 Evaluation of image classification accuracy 实验中,将地物类型分为建设用地、水体、林地、耕地、阴影、其他土地共六类;采用像元平均亮度值、每个波段平均值、相异性、最大差异、形状指数、长宽比及归一化植被指数(normalized difference vegetation index,NDVI)作为提取样本的分类指标。利用混淆矩阵可以计算出评价分类效果的制图精度、用户精度、总体分类精度(邵亚琴
15、等,2022)和Kappa系数(徐辉等,2019)。3 评价实验与结果分析 3.1 实验数据 遥感影像压缩评价实验与测试影像有着直接关系,综合考虑压缩算法与卫星研制工程,开展高分七号卫星影像模拟。进行多光谱遥感影像压缩质量评价,需要综合考虑不同地物地貌的遥感影像受压缩算法的影响,实验选取影像大小为1700像素1700像素、分辨率为2.8 m的高分七号模拟多光谱影像作为测评影像,如图3所示。多光谱实验影像区域选取山东省潍坊市临朐县,影像包括大面积的裸土、山林,以及成片的河流和城市建筑,地物类型丰富。图 3 多光谱影像 Fig.3 Multispectral images 3.2 压缩算法 JPE
16、G-LS的核心是LOCO-I(low complexity lossless compression for images)算法,是一种基于自回归上下文来预测残差的编码方法,其编码方式由被编码像素点的失真控制参数和局部梯度的关 林瑶瑶 等:多光谱遥感影像星上压缩影响评价 331 系决定:如果被编码像素点的局部梯度在失真控制参数范围外,编码进入常规模式;否则,进入游程模式。在JPEG-LS游程编码模式中,编码器以x为起点计算游程长度,遇到不同像素值(对于有损压缩则为其值超出误差限)的图斑或当前行的行尾时,游程终止。JPEG-LS常规编码方式包含重建值预测、预测误差计算和对预测误差进行Golomb
17、编码三个部分(刘向增等,2021)。其原理如 min()max(,)max()min()abcabxabcababcx,xxxxPx,xxx,xxxx,其他 (4)式中,xP为预测值;ax、bx、cx分别为位于a、b、c处的像素值。当x上方存在水平边缘时,选取像素a作为x的预测值;若x左侧存在垂直边缘,选用像素b作为x的预测值;如果算法没有检测到边 缘,则用abc来作为x处的预测值,如图4所示。预测误差为x处像素预测值与实际值的差值,利用Golomb对预测误差进行编码。图 4 预测值选取示意图 Fig.4 Diagram of prediction value selection 3.3 影像
18、特征评价结果 计算未压缩影像与压缩后影像各个波段及多波段灰度平均值和标准方差,如表1所示。31压缩后单个波段和多波段差值统计不大,多波段统计均值差值为0.939,方差差值为0.018;61压缩后多波段均值差值为1.694,方差差值为0.725,各个波段和多波段差值统计较大。表 1 影像灰度平均值、标准方差统计 Tab.1 Statistics of image gray average and standard deviation 未压缩 31 压缩 61 压缩 波段 灰度平均值 灰度标准方差 灰度平均值 灰度标准方差 灰度平均值 灰度标准方差 红 183.454 69.636 184.388
19、 69.780 184.984 70.648 绿 342.804 63.316 343.738 63.514 344.186 64.535 蓝 302.997 31.614 303.965 31.892 305.820 32.290 近红外 331.780 118.549 332.701 118.663 332.821 119.752 多波段 290.259 63.861 291.198 63.843 291.953 63.136 影像压缩前后纹理特征,如表2所示。31压缩后影像纹理差异不明显,多波段对比度差值和信息熵差值为0.042和0.001;61压缩后多波段对比度差值和信息熵差值为0.3
20、27和0.007,各个波段和多波段差值较明显。表 2 影像对比度、信息熵统计 Tab.2 Statistics of image contrast and information entropy 未压缩 31 压缩 61 压缩 波段 对比度 信息熵 对比度 信息熵 对比度 信息熵红 12.374 1.865 12.817 1.888 13.9031.743绿 12.979 1.964 13.004 1.974 15.2101.812蓝 3.588 1.595 3.817 1.606 4.9451.343近红外 47.371 2.121 47.693 2.124 48.4292.106多波段 1
21、4.520 2.011 14.478 2.012 14.8472.004 影像压缩后与原始未压缩影像的差异特征,如表3所示。31压缩后影像相关系数高,峰值信噪比大于50,压缩算法对重建影像影响较小,影像质量良好;单波段61压缩后影像相关系数最小降低到0.891,峰值信噪比小于50,重建影像质量较差。表 3 影像对比特征统计 Tab.3 Statistics of image contrast features 31 压缩 61 压缩 波段 相关性 峰值信噪比 相关性 峰值信噪比红 0.997 51.349 0.971 41.636 绿 0.996 51.318 0.965 41.562 蓝 0
22、.986 51.628 0.891 42.596 近红外0.999 51.22 0.988 41.023 多波段0.999 56.984 0.990 46.990 332 Journal of Spatio-temporal Information 时空信息学报 2023,30(3)3.4 影像分割一致性评价结果 设置分割尺度为28,形状因子为0.2、紧致度因子为0.5,未压缩影像及31、61压缩比情况下的分割变化情况,如图5、图6所示。在房屋、耕地此类地物边界明显的地区,影像执行分割后图斑形状变化较小;在林地、水体此类地物边界模糊的地区,影像执行分割后图斑形状变化相对明显;整体上,61重建影
23、像分割图斑形状较31重建影像分割图斑形状变化明显。图 5 原始影像和重建影像的分割结果示例一 Fig.5 Example of segmentation results for the original image and the reconstructed imageExample 1 图 6 原始影像和重建影像的分割结果示例二 Fig.6 Example of segmentation results for the original and the reconstructed imageExample 2 影像进行不同倍率压缩后,影像光谱特征及纹理特征会发生不同程度的变化,影像分割后图斑
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 光谱 遥感 影像 压缩 影响 评价