Question

MATLAB，求大虾，帮忙给我个图像dct变换和逆变换的算法呗。。具体的dct是怎么实现的？谢谢了。。。

Answer 1

不知道你的dct变换 是用在哪里 数字水印or 图像压缩？？

%程序员：李鹏
%编写时间：2003.11.25%修改时间：2004.3.8
%函数功能：本函数用于DCT域的信息隐藏
function [count,msg,result]=hidedctadv1(image,imagegoal,msg,key,alpha) %按位读取秘密信息
frr=fopen(msg,'r');
[msg,count]=fread(frr,'ubit1');
fclose(frr); data0=imread(image); %将图象矩阵转为double型
data0=double(data0)/255; %取图象的一层做隐藏
data=data0(:,:,1); %对图象分块
T=dctmtx(8); %对分块图象做DCT变换
DCTrgb=blkproc(data,[8 8],'P1*x*P2',T,T'); DCTrgb0=DCTrgb; %产生随机的块选择,确定图像块的首地址
[row,col]=size(DCTrgb);
row=floor(row/8);
col=floor(col/8);
a=zeros([row col]);
[k1,k2]=randinterval1(a,count,key);
for i=1:count
k1(1,i)=(k1(1,i)-1)*8+1;
k2(1,i)=(k2(1,i)-1)*8+1;
end %信息嵌入
temp=0;
for i=1:count
if msg(i,1)==0
if DCTrgb(k1(i)+4,k2(i)+1)>DCTrgb(k1(i)+3,k2(i)+2)
temp=DCTrgb(k1(i)+4,k2(i)+1);
DCTrgb(k1(i)+4,k2(i)+1)=DCTrgb(k1(i)+3,k2(i)+2);
DCTrgb(k1(i)+3,k2(i)+2)=temp;
end
else
if DCTrgb(k1(i)+4,k2(i)+1)<DCTrgb(k1(i)+3,k2(i)+2)
temp=DCTrgb(k1(i)+4,k2(i)+1);
DCTrgb(k1(i)+4,k2(i)+1)=DCTrgb(k1(i)+3,k2(i)+2);
DCTrgb(k1(i)+3,k2(i)+2)=temp;
end
end
if DCTrgb(k1(i)+4,k2(i)+1)>DCTrgb(k1(i)+3,k2(i)+2)
DCTrgb(k1(i)+3,k2(i)+2)=DCTrgb(k1(i)+3,k2(i)+2)-alpha;%将原本小的系数调整得更小
else
DCTrgb(k1(i)+4,k2(i)+1)=DCTrgb(k1(i)+4,k2(i)+1)-alpha;
end
end %信息写回保存
DCTrgb1=DCTrgb;
data=blkproc(DCTrgb,[8 8],'P1*x*P2',T',T); r
esult=data0; result(:,:,1)=data;
imwrite(result,imagegoal);

追问

图像压缩。

追答

下次再发你

追问

20843878@qq.com、 谢谢

追答

q错了吧 注意啊

追问

sorry 。308453878@qq.com

Answer 2

我有求dct系数的程序 。。前天刚写的。