请教matlab抵抗旋转攻击的水印提取代码
在pudn上下载了一个代码，发现提取水印的代码都是限定了图片的大小，如果图片的大小改变了代码就会报错。而进行旋转攻击后图片的大小是要改变的，请问这时候怎么提取水印啊？谁有能抵抗旋转攻击的水印植入和提取代码啊？讲一讲原理也好。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%function averagemark=putout(X,Y)
%函数作用:对含水印图象,用原始图象比较然后提取出水印
%输入变量:X,原始图象的像素值,Y,嵌入水印后的图象
%输出变量:averagemark,提取的水印图象求的平均
%实验不妨固定小波包变换结构为:
%1 2 5 6 13 14
%3 4 7 8 15 16
%9 10 17 18
%11 12 19 20
%21 22
%23 24
% 25 26
% 27 28
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function averagemark=putout(X,Y)
%对原始图象进行小波包变换（实验用固定结构的小波包变换，没有使用验证密码来随机改变这个小波包的%变换结构）
%把X进行三次小波变换，小波基为db1
%一级小波包变换
[ca1,ch1,cv1,cd1]=dwt2(X,'db1');
%二级小波包变换
[ca1_a,ca1_h,ca1_v,ca1_d]=dwt2(ca1,'db1');
[ch1_a,ch1_h,ch1_v,ch1_d]=dwt2(ch1,'db1');
[cv1_a,cv1_h,cv1_v,cv1_d]=dwt2(cv1,'db1');
[cd1_a,cd1_h,cd1_v,cd1_d]=dwt2(cd1,'db1');
%三级小波包变换
[ca1_aa,ca1_ah,ca1_av,ca1_ad]=dwt2(ca1_a,'db1');
[ca1_ha,ca1_hh, ca1_hv, ca1_hd]=dwt2(ca1_h,'db1');
[ca1_va,ca1_vh, ca1_vv, ca1_vd]=dwt2(ca1_v,'db1');
[ch1_ha,ch1_hh, ch1_hv, ch1_hd]=dwt2(ch1_h,'db1');
[ch1_va,ch1_vh, ch1_vv, ch1_vd]=dwt2(ch1_v,'db1');
[cv1_aa,cv1_ah, cv1_av, cv1_ad]=dwt2(cv1_a,'db1');
[cv1_da,cv1_dh, cv1_dv, cv1_dd]=dwt2(cv1_d,'db1');
%组织每块小波块，并标上序号
W1= ca1_aa;
W2=ca1_ah;
W3= ca1_av;
W4=ca1_ad;
W5=ca1_ha;
W6=ca1_hh;
W7=ca1_hv;
W8=ca1_hd;
W9=ca1_va;
W10=ca1_vh;
W11=ca1_vv;
W12=ca1_vd;
W13=ch1_ha;
W14=ch1_hh;
W15=ch1_hv;
W16=ch1_hd;
W17=ch1_va;
W18=ch1_vh;
W19=ch1_vv;
W20=ch1_vd;
W21=cv1_aa;
W22=cv1_ah;
W23=cv1_av;
W24=cv1_ad;
W25=cv1_da;
W26=cv1_dh;
W27=cv1_dv;
W28=cv1_dd;
W=[W1 W2 W3 W4 W5 W6 W7 W8 W9 W10 W11 W12 W13 W14 W15 W16 W17 W18 W19 W20 W21 W22 W23 W24 W25 W26 W27 W28];
%对小波块进行模糊分类
%首先对纹理强度模糊分类
%首先计算每个小波块中小波系数大于1的像素数目,记录在number中
for k=0:27
temp=W(1:32,32*k+1:32*k+32);
num=0;
for i=1:32
for j=1:32
if temp(i,j)>0
num=1+num;
end
end
end
%number(k+1)=num;
% end
%计算每块小波系数的隶属度
%计算对small的隶属度A_small
if num<=400
A_small=1;
end
if (num>400)&(num<=450)
A_small=(450-num)/(450-400);
end
if (num>450)
A_small=0;
end
%计算对middle的隶属度A_middle
if num<=400
A_middle=0;
end
if (num>400)&( num<=450)
A_middle =(num-400)/(450-400);
end
if (num>450)&( num<=500)
A_middle =1;
end
if (num>500)&( num<=550)
A_middle =(550-num)/(550-500);
end
if num>500
A_middle=0;
end
%计算对big的隶属度A_big
if num<=500
A_big=0;
end
if (num>500)&(num<=550)
A_big=(num-500)/(550-500);
end
if (num>550)
A_big=1;
end
%得到模糊分类的结果.返回的结果1代表是small，2代
%表是middle,3代表是big;
if A_small==max([A_small;A_middle;A_big])
fuzzy_result_Texture(k+1)=1;
end
if A_middle==max([A_small;A_middle;A_big])
fuzzy_result_Texture(k+1)=2;
end
if A_big==max([A_small;A_middle;A_big])
fuzzy_result_Texture(k+1)=3;
end
end
%对每个小波块进行能量大小的模糊分类。
%（1）计算每个小波块的能量，值存在e中
e=zeros(1,28);
for k=0:27
temp=W(1:32,32*k+1:32*k+32);
sum=0;
for i=1:32
for j=1:32
sum=temp(i,j)^2+sum;
end
end
e(1+k)=sum;
end
%(2)对每个小波块的能量按小到大排列得到排名号
%注：最大为28号，最小为1号。排名号存在变量paiming
%中
paiming=zeros(1,28);
for i=1:28
for j=i+1:28
if e(i)<=e(j)
s=0; t=1;
else
s=1; t=0;
end
paiming(i)=paiming(i)+s;
paiming(j)=paiming(j)+t;
end
end
paiming=paiming+1;
%计算每块小波块的隶属度
%计算对small的隶属度B_small
for i=1:28
if paiming(i)<=5
B_small(i)=1;
end
if (paiming(i)>5)&(paiming(i)<=10)
B_small(i)=(10-paiming(i))/(10-5);
end
if paiming(i)>10
B_small(i)=0;
end
end
%计算对middle的隶属度B_middle
for i=1:28
if paiming(i)<=5
B_middle(i)=0;
end
if (paiming(i)>5)&( paiming(i)<=10)
B_middle(i) =(paiming(i)-5)/(10-5);
end
if (paiming(i)>10)&( paiming(i)<=15)
B_middle(i) =1;
end
if (paiming(i)>15)&( paiming(i)<=20)
B_middle(i) =(20- paiming(i))/(20-15);
end
if paiming(i)>20
B_middle(i)=0;
end
end
%计算对big的隶属度B_big
for i=1:28
if paiming(i)<=15
B_big(i)=0;
end
if (paiming(i)>15)&( paiming(i)<=20)
B_big(i)=(paiming(i)-15)/(20-15);
end
if paiming(i)>20
B_big(i)=1;
end
end
%得到模糊分类的结果.返回的结果1代表是small，2代
%表是middle,3代表是big;
for i=1:28
if B_small(i)==max([B_small(i);B_middle(i);B_big(i)])
fuzzy_result_Energy(i)=1;
end
if B_middle(i)==max([B_small(i);B_middle(i);B_big(i)])
fuzzy_result_Energy(i) =2;
end
if B_big(i)==max([B_small(i);B_middle(i);B_big(i)])
fuzzy_result_Energy(i) =3;
end
end
%根据模糊规则得到最后的分类结果,存在变量B中,1代%表small,3代表big
for i=1:28
if (fuzzy_result_Texture(i)==3)&( fuzzy_result_Energy(i)==1)
B(i)=3;
end
if (fuzzy_result_Texture(i)==3)&( fuzzy_result_Energy(i)==2)
B(i)=2;
end
if (fuzzy_result_Texture(i)==3)&( fuzzy_result_Energy(i)==3)
B(i)=1;
end
if (fuzzy_result_Texture(i)==2)&( fuzzy_result_Energy(i)==3)
B(i)=1;
end
if (fuzzy_result_Texture(i)==2)&( fuzzy_result_Energy(i)==2)
B(i)=2;
end
if (fuzzy_result_Texture(i)==2)&( fuzzy_result_Energy(i)==1)
B(i)=2;
end
if (fuzzy_result_Texture(i)==1)&( fuzzy_result_Energy(i)==3)
B(i)=2;
end
if (fuzzy_result_Texture(i)==1)&( fuzzy_result_Energy(i)==2)
B(i)=1;
end
if (fuzzy_result_Texture(i)==1)&( fuzzy_result_Energy(i)==1)
B(i)=3;
end
end
%根据B中的值得到相应的水印嵌入强度,存放在 q中
for i=1:28
if B(i)==1
q(i)=0.06;
end
if B(i)==2
q(i)=0.11;
end
if B(i)==3
q(i)=0.15;
end
end
%对含水印图象也做三次小波变换,然后分出小波块
%把Y进行同样的三次小波变换，小波基为db1
%一级小波包变换
[hca1,hch1,hcv1,hcd1]=dwt2(Y,'db1');
%二级小波包变换
[hca1_a,hca1_h,hca1_v,hca1_d]=dwt2(hca1,'db1');
[hch1_a,hch1_h,hch1_v,hch1_d]=dwt2(hch1,'db1');
[hcv1_a,hcv1_h,hcv1_v,hcv1_d]=dwt2(hcv1,'db1');
%三级小波包变换
[hca1_aa,hca1_ah,hca1_av,hca1_ad]=dwt2(hca1_a,'db1');
[hca1_ha,hca1_hh,hca1_hv, hca1_hd]=dwt2(hca1_h,'db1');
[hca1_va,hca1_vh, hca1_vv, hca1_vd]=dwt2(hca1_v,'db1');
[hch1_ha,hch1_hh, hch1_hv, hch1_hd]=dwt2(hch1_h,'db1');
[hch1_va,hch1_vh, hch1_vv, hch1_vd]=dwt2(hch1_v,'db1');
[hcv1_aa,hcv1_ah, hcv1_av, hcv1_ad]=dwt2(hcv1_a,'db1');
[hcv1_da,hcv1_dh, hcv1_dv, hcv1_dd]=dwt2(hcv1_d,'db1');
%组织每块小波块，并标上序号
HW1= hca1_aa;
HW2=hca1_ah;
HW3= hca1_av;
HW4=hca1_ad;
HW5=hca1_ha;
HW6=hca1_hh;
HW7=hca1_hv;
HW8=hca1_hd;
HW9=hca1_va;
HW10=hca1_vh;
HW11=hca1_vv;
HW12=hca1_vd;
HW13=hch1_ha;
HW14=hch1_hh;
HW15=hch1_hv;
HW16=hch1_hd;
HW17=hch1_va;
HW18=hch1_vh;
HW19=hch1_vv;
HW20=hch1_vd;
HW21=hcv1_aa;
HW22=hcv1_ah;
HW23=hcv1_av;
HW24=hcv1_ad;
HW25=hcv1_da;
HW26=hcv1_dh;
HW27=hcv1_dv;
HW28=hcv1_dd;
HW=[HW1 HW2 HW3 HW4 HW5 HW6 HW7 HW8 HW9 HW10 HW11 HW12 HW13 HW14 HW15 HW16 HW17 HW18 HW19 HW20 HW21 HW22 HW23 HW24 HW25 HW26 HW27 HW28];
%比较每个小波块的差别,把每块的差别存在difference变量中
difference=HW-W;
%从difference中求出嵌入的每块水印的平均值
sum=zeros(32,32);
for i=0:27
temp=difference(1:32,32*i+1:32*i+32);
temp=temp/q(1,i+1);
sum=sum+temp;
end
averagemark=sum/28;
