help to fix this

2 visualizaciones (últimos 30 días)
Hamrit
Hamrit el 10 de Jun. de 2012
clear all
clc
global Gdi
global c
O=imread('12.jpg');
D1= im2double(O);
for ng=1:3
XD=D1(:,:,ng);
[lg,cg]=size(XD);
lg1=lg/8;
lg2=fix(lg1);
g1=lg2*8;
cg1=cg/8;
cg2=fix(cg1);
g2=cg2*8;
ofs = ones(g1,g2)*128;
p=0;
t=1;
v=1;
while v<=(g1-7)
while t<=(g2-7)
pix=XD(v:(v+7),t:(t+7));
pig=pix*256;
pg=fix(pig);
p=p+1;
Gm(:,:,p)=pg(:,:,1);
%____________________________________________________________________
for i=1:1:1
for j=1:1:8
b(i,j)=1/sqrt(8);
end
end
for i=2:1:8
for j=1:1:8
b(i,j)=sqrt(.25)*cos((2*(j-1)+1)*(i-1)*pi/16);
end
end
%dtr=b';
c(:,:,p)=b*Gm(:,:,p);
d(:,:,p)=c(:,:,p)*b';
t=t+8;
end
t=1;
v=v+8;
end
for z=1:p
n=[3 5 7 9 11 13 15 17;
5 7 9 11 13 15 17 19;
7 9 11 13 15 17 19 21;
9 11 13 15 17 19 21 23;
11 13 15 17 19 21 23 25;
13 15 17 19 21 23 25 27;
15 17 19 21 23 25 27 29;
17 19 21 23 25 27 29 31];
QG(:,:,z)=d(:,:,z)./n;
DQG(:,:,z)=fix(QG(:,:,z));
%$$$$$$$$$ DECOMPRESSION $$$$$$$$$$$$$$$$$$$$$$$
GDCTi(:,:,z)=DQG(:,:,z).*n;
Go(:,:,z)= IDCT (GDCTi(:,:,z));
for i=1:1:1
for j=1:1:8
b(i,j)=1/sqrt(8);
end
end
for i=2:1:8
for j=1:1:8
b(i,j)=sqrt(.25)*cos((2*(j-1)+1)*(i-1)*pi/16);
end
end
%dtr=b';
c(:,:,z)=b*GDCTi(:,:,z);
Gdi(:,:,z)=c(:,:,z)*b;
end
%Gdi =uint8(Gdi(:,:,z));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%ùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùùù
p=1;
pas_sbloc=8;
for i=1:pas_sbloc:lg
for j=1:pas_sbloc:cg
C(i:pas_sbloc+i-1,j:pas_sbloc+j-1)=Gdi(:,:,p);
p=p+1;
end
end
if ng==1
Rh=C;
elseif ng==2
Gh=C;
else
Bh=C;
end
end
figure(2);
A=cat(3,Rh,Gh,Bh);
imshow(uint8(A));
  1 comentario
Walter Roberson
Walter Roberson el 10 de Jun. de 2012
In what manner is it broken? Do you get an error message? How does the output differ from your expectations?

Iniciar sesión para comentar.

Respuestas (0)

Etiquetas

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by