TA的每日心情 | 衰 2019-11-19 15:32 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
5 _3 B& [/ ^0 e. S$ `6 K+ Q3 h一、简介+ V2 x3 I9 U1 w; w% h6 G+ a
数字图像通常需要很大的比特数,这给图像的传输和存储带来相当大的困难。例如用8bit存储一幅512512的灰度图像的比特数为256K。而一部60分钟的彩色电影,如果每秒放映24帧,数字化后每帧包含512512像素,每像素的R、G、B分量分别占8bit,则这样一部电影的总比特数为64800M,若用一张600M的CD存储则需要100多张CD光盘来存储。由此可见对图像数据进行压缩显得非常必要。减少存储空间、缩短传输时间这成为促进图像压缩编码技术发展的主导因素。图像压缩是通过编码来实现的,所以通常将压缩与编码统称为图像的压缩编码。图像压缩编码从本质上来说就是对要处理的图像数据按照一定的规则进行变换和组合,从而达到以尽可能少的数据来表示尽可能多的数据信息。
+ p4 r6 {8 O1 W; j- _5 `, D
( O6 j3 x4 r8 r5 t1.数据的冗余
! \6 _1 h8 K% |8 |, }
8 A! [. G1 r# N) K大多数信息的表达都存在着一定的冗余度(相关性),通过采用一定的模型和编码方法,可以降低这种冗余度。0 h1 w1 E/ H* @; X- N
' l8 x9 z& L$ P& T
设:为表达一个信息,方法A用了n1个数据,改进方法B用了n2个数据(压缩数据):
$ A$ U4 w7 V8 }7 T e7 k9 ?. c- l) C# W; F3 B5 }* z
压缩率(压缩比):C = n1 / n2
6 y. p( W! n3 C& W" h e
4 H1 J8 ?8 ^: x) F& g) |$ D相对数据冗余:R = 1 – 1/C) u! s# W0 b; R4 @
- V1 p. J0 F9 O: B& D1 c2.图像数据中的冗余0 F' P8 h' A( w6 P' z
! l, W1 P# H t, ^' K+ ]% p N
三种数据冗余:编码冗余、像素冗余、视觉心理冗余
9 ^; f& Z/ V9 k' T8 p- ?: \( @- G3 l1 I: ?2 C* v
(1)编码冗余:如果一个图像的灰度级编码,使用了多于实际需要的编码符号,就称该图像包含了编码冗余。
( m# X, I: V& G; I# _4 l
: u% N. U! V% R7 I' I: ~: s(2)像素冗余(像素间相关)$ h" P% {* `, ~7 P0 D. W$ a: B
" }6 ], ? }6 N1 l: E% D4 T V任何给定位置像素值,原理上均与相邻像素相关,都可以通过它的邻居预测到,单个像素携带的信息相对较少。. C% H) _6 D K- U' ]$ t% B
8 ~; v) K4 Y: U7 E( t& R/ }! y% w
对于一幅图像,大量单个像素对视觉的贡献是冗余的。这是建立在对邻居值预测的基础上。. X5 y7 c6 y2 ^( ]- Z" }
' i% X6 i7 M4 ]$ y+ ^例:原图像数据:234 223 231 238 235
( |+ t2 W' V. w& X! @- v7 l" D6 S0 ^3 K3 }5 L
压缩后数据: 234 -11 8 7 -35 }, W. j" V' o5 P, T. d
+ w$ [2 ]% t: w2 e7 m/ Z
% D0 u& ]& A4 T8 B7 b6 P9 J; J(3)视觉心理冗余
4 o- [& h9 C: A* Y) j0 D! {3 b1 ~+ X3 g/ s- f" a1 M% n& m2 B$ c
一些信息在一般视觉处理中比其它信息的相对重要程度要小,这种信息就被称为视觉心理冗余。
, P5 A; y$ n- O1 x
6 n* S7 W. T8 ?" D7 t% @人类的视觉系统由于受生理特性的限制,对于图像场的注意是非均匀的,人对细微的颜色差异感觉不明显。 例如,人类视觉的一般分辨能力为26灰度等级,而一般的图像的量化采用的是28灰度等级,即存在视觉冗余。 人类的听觉对某些信号反映不太敏感,使得压缩后再还原有允许范围的变化,人也感觉不出来。
/ h8 u8 I" p8 x& s# c- {% n N( y
3.利用图像数据的冗余开展压缩
l* T; n, a' ^( O4 R' q) t( A5 ^6 `9 d
(1)因为有编码冗余和像素冗余,当我们将图像信息的描述方式改变之后,可以压缩掉这些冗余,进行“无损”压缩。" u. {7 j' Z" N4 k8 h& a
& r7 T4 _7 e" Y% t/ r( g(2)因为有主观视觉冗余,当我们忽略一些视觉不太明显的微小差异,可以进行所谓的“有损”压缩。+ K( j9 ^, X2 s; Y2 \& @9 v
8 v l6 X4 O1 M/ X; M
$ n2 b# k+ u+ a5 q4 c9 X二、源代码% [. V7 C) t q0 E# O
- close all; clear all; clc;
- I=imread('lena.jpg'); %读入图像
- I1=rgb2gray(I);
- I2=I1(:); %将原始图像写成一维的数据并设为 I2
- I2length=length(I2); %计算 I2 的长度
- I3=im2bw(I1,0.5); %将原图转换为二值图像,阈值为 0.5
- %以下程序为对原图像进行行程编码,压缩
- X=I3(:); %令 X 为新建的二值图像的一维数据组
- L=length(X);
- j=1;
- I4(1)=1;
- for z=1:1:(length(X)-1) %行程编码程序段
- if X(z)==X(z+1)
-
- I4(j)=I4(j)+1;
- else
- data(j)=X(z); % data(j)代表相应的像素数据
- j=j+1;
- I4(j)=1;
- end
- end
- data(j)=X(length(X)); %最后一个像素数据赋给 data
- I4length=length(I4); %计算行程编码后的所占字节数,记为 I4length
- CR=I2length/I4length; %比较压缩前于压缩后的大小
- %下面程序是行程编码解压
- l=1;
- for m=1:I4length
- for n=1:1:I4(m);
- decode_image1(l)=data(m);
- l=l+1;
- end
- end
- decode_image=reshape(decode_image1,256,256); %重建二维图像数组
- h=figure,set(h,'color','white')
- x=1:1:length(X);
- subplot(131),plot(x,X(x)); %显示行程编码之前的图像数据
- y=1:1:I4length ;
- subplot(132),plot(y,I4(y)); %显示编码后数据信息
- u=1:1:length(decode_image1);
- subplot(133),plot(u,decode_image1(u)); %查看解压后的图像数据
- h=figure;
- subplot(131);imshow(I); %显示原图的二值图像
- subplot(132);imshow(I3); %显示原图的二值图像
- subplot(133),imshow(decode_image); %显示解压恢复后的图像
- set(h,'color','white')
- disp('压缩比: ')
- disp(CR);
- disp('原图像数据的长度:')
- disp(L);
- disp('压缩后图像数据的长度:')
- disp(I4length);
- disp('解压后图像数据的长度:')
- disp(length(decode_image1));. X" _* e" J- Q- X, |1 [3 p
* Y( k7 o; h/ s: l! M; A; w( V: a3 ~. r9 r+ P& L3 D+ T
三、运行结果- y- N( o$ F# U c
& k Q0 s+ X+ [" U# I+ d/ b2 s0 C
, ?1 Y- c3 X% _ N
+ |: [+ ~/ w1 m. J1 L) {
& l7 R4 b+ w: o8 U P: Q ]8 }
, t4 N. j; b. j. t k1 X8 x: I* e! P+ S' u3 b6 x. J& i
|6 |6 j: y9 N7 @1 D" c* p! m# a
% B4 I1 X, ^0 F( W! ?. r+ D
|
|