第八章圖象壓縮

第八章圖象壓縮第一頁，共三十四頁，2022年，8月28日

2．圖象一大特點是數(shù)據(jù)量大，給存貯、傳輸帶來困難，需壓縮

例：電話線傳輸速率一般為56Kbits/s（波特率）一幅彩色圖像512×512×24bit=6Mbits大小。傳一幅圖象需2分鐘左右。實時傳送更難：512×512×24bits×25幀/秒=150Mbits/S

如壓縮20倍，傳一幅圖6秒左右，可以接受，實用。實時，要專用信道（衛(wèi)星、微波網(wǎng)、專線網(wǎng)等技術）；另外，大量資料需存貯遙感、故宮、醫(yī)學CT、MR，需要壓縮。第二頁，共三十四頁，2022年，8月28日圖像數(shù)據(jù)可壓縮性分析：從統(tǒng)計觀點出發(fā)，簡單直觀地講，就是某一象素的灰度值，總是和周圍其他象素灰度值有某種關系，應用某種編碼方法提取并減少這些相關特性，便可實現(xiàn)圖象信息的數(shù)據(jù)壓縮。圖象數(shù)據(jù)壓縮的目的可以是節(jié)省圖象存貯器容量，也可以是減少傳輸信道容量，還可以是為了縮短圖象加工處理的時間。第三頁，共三十四頁，2022年，8月28日編碼例子：8個灰級，總計64*64=4096像素的圖像DAnkPr(rk)編碼1長度編碼2長度07900.1900031121/710230.2500130122/78500.2101031023/76560.16011300134/73290.081003000145/72450.0610130000156/71220.03110300000167/7810.0211130000006編碼1平均碼長：3bits;編碼2平均碼長：2.7bits第四頁，共三十四頁，2022年，8月28日（1）壓縮方法：按時間分：靜圖：靜止圖象（要求質量高）動圖：活動的序列圖象（相對質量要求低，壓縮倍數(shù)要高）

壓縮比＝未壓縮的圖象的存貯字節(jié)數(shù)／壓縮后圖象存貯字節(jié)數(shù)按失真與否分：無失真壓縮：經(jīng)壓縮后再恢復圖象與原圖象無任何區(qū)別，一般壓縮倍數(shù)<2

有限失真壓縮：單幀（靜）4～20倍。圖象序列（x、y、t）50～200倍

第五頁，共三十四頁，2022年，8月28日（2）壓縮失真度衡量標準：a．客觀：輸出與輸入之差：e(x,y)=g(x,y)－f(x,y)

整圖差：均方差：

量綱：灰度級。第六頁，共三十四頁，2022年，8月28日也可將輸出與輸入之差看作噪聲，此時：輸出圖象均方信噪比定義為：g2(x,y)/e2(x,y)第七頁，共三十四頁，2022年，8月28日常用：峰值信噪比

一般σ<2，或PSNR≥40dB

人眼看不出來

<30dB的圖象不能用～35dB可接受

[到目前為上，很難找出一個通用的評價圖象壓縮的客觀標準，仍然是很好的研究課題]第八頁，共三十四頁，2022年，8月28日b．主觀（人判別）－－專家投票的方法，實用。

人的視覺的主觀亮度是光強的對數(shù)函數(shù)。人眼對黑暗區(qū)誤差比明亮區(qū)更敏感。人眼對灰度突變邊緣比較敏感。第九頁，共三十四頁，2022年，8月28日

（3）圖象壓縮的技術指標：a．保真度－－與用途有關，例如偵察與體育比賽關心內容不同b．壓縮比：原圖象數(shù)據(jù)量／壓縮后數(shù)據(jù)量或0.3bit/pixelc．誤碼擴散程度d．實時性－－與壓縮算法、系統(tǒng)速度有關e．保密性－－傳輸中防止被盜，SPOT是數(shù)據(jù)需解碼（高明壓縮方法，壓縮后仍是一幅圖）數(shù)字水印技術（信息安全）－－watermarking第十頁，共三十四頁，2022年，8月28日8.2四叉樹（金字塔pyramid或Quadtree）適用于N=2n圖象四叉樹：

0層1層2層……n層象素數(shù)：

4041424n四叉樹全部保留，總數(shù)據(jù)量為：全分辨率圖：

4n點數(shù)，數(shù)據(jù)量增加1/3，對搜索有好處，找到有關興趣區(qū)域后再提取感應趣的下一層圖象。第十一頁，共三十四頁，2022年，8月28日8.3信息有關的基本概念1．信息的度量：

7月份天氣報：“明天下雨”，信息量小

10月份天氣報：“明天下雪”－－小概率事件發(fā)生，信息量大即：概率低－－信息量大概率高－－信息量小怎么定量？戈壁灘一年下雨概率1/128，睛概率127/128。對天氣預報的信息量多少？引入-log2p代表信息量如：-log21/128=7

有雨

-log2127/128≈

0晴第十二頁，共三十四頁，2022年，8月28日信息具有可加性：一個9×9方陣，一個士兵在某一位置的概率1/81，信息=-log21/81=6.34

在特定行上的概率1/9：信息=-log21/9=3.17

在特定列上的概率1/9：信息=-log21/9=3.17

行、列信息量和：3.17+3.17=6.34引入一個信息度量I(P)

特點：0p<1信息量I(P)連續(xù)變化P=0，I(P)=P=1，I(P)=0P1>P2，I(P2)>I(P1)I(p?q)=I(p)+I(q)I(P)=-logbp—

一般b=2第十三頁，共三十四頁，2022年，8月28日2．信息熵（Entropy）

預期能得到的平均信息量，I(Pi)的數(shù)學期望值。對于M個灰階：Pii=1,2,….,M例1，當Pi=1/4,i=1,2,3,4H=2

當P1=1/2,P2=1/4,P3=P4=1/8H=1.75第十四頁，共三十四頁，2022年，8月28日例2，二值圖象，xi=0,10值概率：p1值概率：1-pH～p變化關系可以發(fā)現(xiàn)，當p=0.5時熵最大。

00.51pH1.0第十五頁，共三十四頁，2022年，8月28日一般講，對M個灰值：

求H的條件極值，即：

當

M=256時，Hmax=8注意熵H的特點：

H>0Hmax=log2MH<Hmax第十六頁，共三十四頁，2022年，8月28日8.4可逆保真編碼（無失真壓縮）熵：平均信息量一般圖象壓縮模式：代表：孤立地考察每個象素點出現(xiàn)概率。去相關（去冗余）量化無失真編碼S=8大S小0255灰度

百分比第十七頁，共三十四頁，2022年，8月28日注意：恰當做變換熵可以減小。例：（板書）熵變小，把信息轉移到變換模型中（高階熵概念）

（舉例：降低比特位問題）結論：

不作變換時無失真編碼壓縮極限是熵第十八頁，共三十四頁，2022年，8月28日一幅圖(8*8*3bits)4444444045555450456664504567645045666450455554504444444044444440灰度頻數(shù)概率熵080.1250.3751000200030004310.4840.5075160.250.50680.1250.375710.0160.096和641.01.853第十九頁，共三十四頁，2022年，8月28日1．Huffman碼－－緊湊編碼

n個由(m1…mn)信息，概率為P1…Pn，基本思想：對概率大的信息分配短碼，小的分配長碼。步驟：先作直方圖統(tǒng)計，按出現(xiàn)概率由大到小排隊。而后反向編碼：概率大的賦0，概率小的賦1。實現(xiàn)：可以按塊、按總體統(tǒng)計或自適應邊統(tǒng)計邊修正。特點：此碼變長，不易實現(xiàn)，當概率變，則編碼方式亦變化。

第二十頁，共三十四頁，2022年，8月28日2．B碼：每個碼分成兩部分，即延續(xù)比特和信息比特。分割碼字代表有效信息C×C×C×C×……(×：0或1表示)C相同的表示同一個碼字

C變號表示前一個碼字結束。

3．移位碼

2位碼可有00，01，10，而11作為移位碼，超過3的倍數(shù)用11的次數(shù)表之。

第二十一頁，共三十四頁，2022年，8月28日4．行程編碼(Run--LengthCoding)

由灰度長度對組成，即(灰度，長度)，根據(jù)前圖按垂直方向掃描,可編為：(4,9)(5,5)(4,3)(5,1)(6,3)(5,1)(4,3)(5,1)(6,1)(7,1)(6,1)(5,1)(4,3)(5,1)(6,3)(5,1)(4,3)(5,5)(4,10)(0,8)

結果：64個數(shù)據(jù)壓成到40個。通常掃描走向：第二十二頁，共三十四頁，2022年，8月28日另一種掃描方向線選取，Helbert曲線－－行程編碼。遞歸算法：

A:DAABB:CBBAC:BCCDD:ADDCDDDAC1級2級第二十三頁，共三十四頁，2022年，8月28日對前圖用Helbert曲線掃描結果：

4454445566765544455444444444556665554444440004400044566555444400得：(4,2)(5,1)(4,3)(5,2)(6,2)(7,1)(6,1)(5,2)(4,3)(5,2)(4,9)(5,2)(6,3)(5,3)(4,6)(0,3)(4,2)(0,3)(4,2)(5,1)(6,2)(5,3)(4,4)(0,2)24×2=48[目前算術編程碼效率最高,計算量大]第二十四頁，共三十四頁，2022年，8月28日8.5有限失真編碼1．插值編碼量化：灰度分辨率采樣：空間分辨率僅保留*部分灰值其它部分值求得：*o*o*o*o*o*o*o*o*o*o*o*o*o*o*s3*o*s1s0s2

*s4第二十五頁，共三十四頁，2022年，8月28日2．按塊取齊編碼將N×N圖像分成M×M大小之塊，M≤N然后分別對每個M×M塊內之數(shù)據(jù)壓成一位長，按力矩不變原則選閾值并二值化而得到，同時還需保留恢復力矩必須的信息。步驟：由中值分割成二值圖，由一、二次力矩計算出分別之A，B值，保留二值圖像及f,σ。

例：

原圖

114564737200247255160121694357251

恢復后圖204204171717204204204171717204171717204

二值圖

1100011100010001第二十六頁，共三十四頁，2022年，8月28日小范圍內恢復后的灰度值近似取為一樣，看不出差別，主要特征保留。

優(yōu)點：按塊等長編碼，防誤碼擴散。缺點：四行才能計算，需要行存貯。第二十七頁，共三十四頁，2022年，8月28日3．變換域編碼

概念：將空間域里的圖象，經(jīng)過某種變換（通常用正交變換，如付氏變換，沃爾什變換等），在某種變換域里描述，一般來說，在變換域里描述要比空間域里簡單，而且圖象相關性明顯下降，這樣變換本身就能實現(xiàn)編碼率的壓縮。若再對變換域圖象進行某種處理（如頻域中的二維濾波）以及熵編碼，則可進一步壓縮圖象的編碼比特率。對變換處理后的圖像信號施以上述的反變換，即可獲得空間域圖象。變換目的是去相關、相鄰象素相關性強。變換運算量與N2成正比。一般N取8×8塊或16×16塊（塊小壓縮倍數(shù)少，沒意義；塊大，相關性弱）。第二十八頁，共三十四頁，2022年，8月28日4．預測編碼：由前m點預測當前點，存預測值與實際值的差值△典型：相鄰象素灰值相減，存差值。

預測結果：多米諾骨牌，誤碼擴散嚴重。

分塊后邊緣損失。

第二十九頁，共三十四頁，2022年，8月28日5．矢量編碼簡化描述：全圖n個點，有A0A1A2…An，有m個組，近似歸納為有限個模式，如k組，全圖有n/m組，n/m》k,

存k組+n/m個組號。

第三十頁，共三十四頁，2022年，8月28日8.6圖象幀間編碼序列圖象（動圖）幀間相關性強1．靜止運動部分的關系：

Bell實驗室研究成果，人對靜止部分分辨率強，對運動部分分辨率弱－－要求不高。即：空間分辨率高，時間分辨率低一些。2．傳送幀間差－－對幀間差壓縮傳輸，兩幅之間對應象素的灰度差小于某一閾值，取為零。3．運動檢測：圖象分成一定大小的塊（M×N—16×16）子塊，在(M+2L,N+2L)范圍內到前一幀圖上搜索與某相關性最大的子塊。

相關函數(shù)：

第三十一頁，共三十四頁，2022年，8月28日8.7JPEG與MPEG（J