shiyan.m

来自「一种利用ARNOLD变换和DCT变换的图象处理方法,附带图象可以马上实验!」· M 代码 · 共 65 行

I = imread('lena.bmp');
subplot(2,2,1)
        imshow(I)
        fun = @dct2;
        J = blkproc(I,[8 8],fun);
        
        subplot(2,2,2)
        image(J)
            %对水印进行置乱
         Q=double(imread('200502.bmp'));
M=Q;
Size_Q=size(Q);
n=0;
Frequency=1;
K=Size_Q;
M1_t=Q;
 M2_t=Q;
for s=1:Frequency
   
    n=n+1;
    if mod(n,2)==0
        for i=1:K(1)
            for j=1:K(2)
                c=M2_t(i,j);
                M1_t(mod(i+j-2,K(1))+1,mod(i+2*j-3,K(2))+1)=c;
            end
        end
    else
         for i=1:K(1)
            for j=1:K(2)
                c=M1_t(i,j);
                M2_t(mod(i+j-2,K(1))+1,mod(i+2*j-3,K(2))+1)=c;
            end
         end
    end
end
if mod(s,2)==0
    %mod(Frequency,2)==0
    M=M1_t;
else
     M=M2_t;
end

        %提取中频系数
     
        for j=0:31;
            for i=0:31;
                m=2+8*i;
                n=3+8*j;
               B(i+1,j+1)=J(m,n);
           %N(i+1,j+1)=B(i+1,j+1)+M(i+1,j+1) 
            end
        end
 N=B+M;
 
 
 JJ=J;
 for i=0:31
     for j=0:31
         JJ(2+8*i,3+8*j)=JJ(2+8*i,3+8*j)+N(i+1,j+1);
     end
 end
 
%         imshow(JJ)