sss.rar_image sharing _sss_z

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %Detecting Duplicated Image Regions % % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function [Du_image]=DupImRe(image) im=image; %读入图像 Du_image={}; if( ndims(im) == 3 ) %判断图像的维数，确定图像为彩色图像还是灰度图像 TYPE = 'c'; % color image else TYPE = 'g'; % grayscale image end b0=8; %设定分块边长 b=b0^2; %分块像素数 %e=0.01; %提取精度 Q=256; %量化常数 Nn=100; %相似行阈值 Nd=16; %分块位移阈值 %tem_D=0; fre_offset=[]; if TYPE=='g' [width,height]=size(im); %读取图像大小 N=width*height; %图像总像素数 Nb0=sqrt(N)-b0+1; Nb=Nb0^2; %分块个数 E0=zeros(b); X=zeros(b,Nb); X0=zeros(b,Nb); A0=zeros(Nb,b); p=1; C=zeros(b); Co=zeros(2,Nb); for i=1:Nb0 for j=1:Nb0 X(1:b,p)=div_block(i,j,b0,1,im); %转化分块像素为向量形式 Co(1,p)=i; %存储分块坐标 Co(2,p)=j; p=p+1; end end U=sum(X')'/Nb; for i=1:Nb %将向量组转化为零均值 X0(:,i)=X(:,i)-U; end for i=1:Nb %计算协方差矩阵 C=C+X0(:,i)*X0(:,i)'; end [E,D]=eig(C); [D0,d_sort]=sort(diag(D),'descend'); %排序特征向量组成基础解系 %sum_D0=sum(D0); %for i=1:b % tem_D=tem_D+D0(i); % if tem_D>(1-e)*sum_D0 % Nt=i; % break % end %end for i=1:b E0(:,i)=E(:,d_sort(i)); end A=X0'*E0; %产生系数矩阵 %[g,h]=size(A) for i=1:b A0(:,i)=fix(( A(:,i)-ones(Nb,1)*min(A(:,i)))/(max(A(:,i))-min(A(:,i)))*(Q-1))+1; end %量化系数矩阵 A1=[A0 Co']; Result=sortrows(A1,1:b); for i=1:Nb %累积位移频次 for j=i:min((i+Nn),Nb) tempx=abs(Result(i,b+1)-Result(j,b+1)); if tempx==0 tempy=abs(Result(i,b+2)-Result(j,b+2)); else tempy=Result(i,b+2)-Result(j,b+2); end distance=sqrt(tempx^2+tempy^2); if distance>Nd fre_offset=dis_save(tempx,tempy,fre_offset); end end end else if TYPE=='c' [width,height]=size(im); %读取图像大小 b=b*3; N=width*height/ndims(im); %图像总像素数 Nb0=sqrt(N)-b0+1; Nb=Nb0^2; %分块个数 E0=zeros(b); X=zeros(b,Nb); X0=zeros(b,Nb); A0=zeros(Nb,b); p=1; C=zeros(b); Co=zeros(2,Nb); for i=1:Nb0 for j=1:Nb0 X(1:b,p)=div_block(i,j,b0,3,im); %转化分块像素为向量形式 Co(1,p)=i; %存储分块坐标 Co(2,p)=j; p=p+1; end end U=sum(X')'/Nb; for i=1:Nb %将向量组转化为零均值 X0(:,i)=X(:,i)-U; end for i=1:Nb %计算协方差矩阵 C=C+X0(:,i)*X0(:,i)'; end [E,D]=eig(C); [D0,d_sort]=sort(diag(D),'descend'); %排序特征向量组成基础解 %sum_D0=sum(D0); %for i=1:b % tem_D=tem_D+D0(i); % if tem_D>(1-e)*sum_D0 % Nt=i; % break % end %end for i=1:b E0(:,i)=E(:,d_sort(i)); end A=X0'*E0; %产生系数矩阵 for i=1:b %量化系数矩阵 A0(:,i)=fix(( A(:,i)-ones(Nb,1)*min(A(:,i)))/(max(A(:,i))-min(A(:,i)))*(Q-1))+1; end A1=[A0 Co']; Result=sortrows(A1,1:b); for i=1:Nb %累积位移频次 i for j=i:min((i+Nn),Nb) tempx=abs(Result(i,b+1)-Result(j,b+1)); if tempx==0 tempy=abs(Result(i,b+2)-Result(j,b+2)); else tempy=Result(i,b+2)-Result(j,b+2); end distance=sqrt(tempx^2+tempy^2); if distance>Nd fre_offset=dis_save(tempx,tempy,fre_offset); end end end end end cd c:\my_dir save dup Result fre_offset;