l3_4.m

来自「DigitalImageProcessing_base_on_Matlab 基」· M 代码 · 共 33 行

%程序L3_4.m:利用LGB训练三个不同大小与维度的码，本并分别进行向量量化%
load lena
th=0.05;
%训练大小为64、维度为4的码本cb_64x4%
nc=64;nd=4;
[cb_64x4]=LBG(X,nc,nd,th);
%训练大小为64、维度为16的码本cb_64x16%
nc=64;nd=16;
[cb_64x16]=LBG(X,nc,nd,th);
%训练大小为128、维度为4的码本cb_128x4%
nc=128;nd=4;
[cb_128x4]=LBG(X,nc,nd,th);

%使用cb_64x4的码本进行向量量化%
[M,N]=size(X);
[Y_64x4,bit_64x4]=VQ(X,cb_64x4);      %VQ函数在本程序范例中
%计算PSNR和CR%
MSE=(sum(sum((X-Y_64x4).^2)))/(M*N);
PSNR_64x4=20*log10(255/sqrt(MSE))
CR_64x4=M*N*8/bit_64x4

%使用cb_64x16的码本进行向量量化%
[Y_64x16,bit_64x16]=VQ(X,cb_64x16);
MSE=(sum(sum((X-Y_64x16).^2)))/(M*N);
PSNR_64x16=20*log10(255/sqrt(MSE))
CR_64x16=M*N*8/bit_64x16

%使用cb_128x4的码本进行向量量化%
[Y_128x16,bit_128x4]=VQ(X,cb_128x4;
MSE=(sum(sum((X-Y_128x4).^2)))/(M*N);
PSNR_128x4=20*log10(255/sqrt(MSE))
CR_128x4=M*N*8/bit_128x4