clear;
%装入原图像1
load woman;
I=X;
%小波函数
type = 'db1';
% 2维离散Daubechies小波变换
[CA1, CH1, CV1, CD1] = dwt2(I,type);
C1 = [CH1 CV1 CD1];
%系数矩阵大小
[length1, width1] = size(CA1);
[M1, N1] = size(C1);
% 定义阈值T1
T1 =50;
alpha = 0.2;
%在图像中加入水印
for counter2 = 1: 1: N1
for counter1 = 1: 1: M1
if( C1(counter1, counter2) > T1 )
marked1(counter1,counter2) = randn(1,1);
NEWC1(counter1, counter2) = double( C1(counter1, counter2) ) +...
alpha * abs( double( C1(counter1, counter2) ) ) * marked1(counter1,counter2) ;
else
marked1(counter1, counter2) = 0;
NEWC1(counter1, counter2) = double( C1(counter1, counter2) );
end;
end;
end;
%重构图像
NEWCH1 = NEWC1(1:length1, 1:width1);
NEWCV1 = NEWC1(1:length1, width1+1:2*width1);
NEWCD1 = NEWC1(1:length1, 2*width1+1:3*width1);
R1 = double( idwt2(CA1, NEWCH1, NEWCV1, NEWCD1, type) );
%分离水印
watermark1 = double(R1) - double(I);
figure(1);
subplot(1,2,1);
image(I);
axis('square');
title('原始图像');
subplot(1,2,2);
imshow(R1/250);
axis('square');
title('Daubechies小波变换后图像');
figure(2);
imshow(watermark1*10^16);
axis('square');
title('水印图像');
% 水印检测
newmarked1 = reshape(marked1, M1*N1, 1);
% 检测阈值
T2 = 60;
for counter2 = 1: 1: N1
for counter1 = 1: 1: M1
if( NEWC1(counter1, counter2) >T2 )
NEWC1X(counter1, counter2) = NEWC1(counter1, counter2);
else
NEWC1X(counter1, counter2) = 0;
end;
end;
end;
NEWC1X = reshape(NEWC1X, M1*N1, 1);
correlation1 = zeros(1000,1);
for corrcounter = 1: 1: 1000
if( corrcounter == 500)
correlation1(corrcounter,1) = NEWC1X'*newmarked1 / (M1*N1);
else
rnmark = randn(M1*N1,1);
correlation1(corrcounter,1) = NEWC1X'*rnmark / (M1*N1);
end;
end;
% 计算阈值
originalthreshold = 0;
for counter2 = 1: 1: N1
for counter1 = 1: 1: M1
if( NEWC1(counter1, counter2) > T2 )
originalthreshold = originalthreshold + abs( NEWC1(counter1, counter2) );
end;
end;
end;
originalthreshold = originalthreshold * alpha / (2*M1*N1);
corrcounter = 1000;
originalthresholdvector = ones(corrcounter,1) * originalthreshold;
figure(3);
plot(correlation1, '-');
hold on;
plot(originalthresholdvector, '--');
title('原始的加水印图像');
xlabel('水印');
ylabel('检测响应');
没有合适的资源?快使用搜索试试~ 我知道了~
资源推荐
资源详情
资源评论
收起资源包目录
小波分析理论与Matlab实现 (153个子文件)
program_19_01.m 2KB
program_11_14.m 2KB
program_04_01.m 2KB
program_11_10.m 2KB
program_16_04.m 2KB
program_16_03.m 2KB
program_16_01.m 2KB
program_02_01.m 1KB
program_20_1.m 1KB
program_16_07.m 1KB
program_16_02.m 1KB
program_13_14.m 1KB
program_10_16.m 1KB
program_16_06.m 1KB
program_16_05.m 1KB
program_09_10.m 1KB
program_16_14.m 1KB
program_16_08.m 1KB
program_13_17.m 1KB
program_11_11.m 1KB
program_11_04.m 1KB
program_13_01.m 1KB
program_11_05.m 1KB
program_17_05_SPECIAL.m 1KB
program_10_13.m 1KB
program_16_09.m 1022B
program_09_09.m 991B
program_11_06.m 986B
program_17_04_SPECIAL.m 980B
program_18_02_SPECIAL.m 972B
program_18_01_SPECIAL.m 967B
program_10_02.m 953B
program_09_08.m 928B
program_16_15.m 911B
program_16_10.m 876B
program_18_07_SPECIAL.m 862B
program_18_03_SPECIAL.m 855B
program_13_16.m 833B
program_18_06_SPECIAL.m 813B
program_16_11.m 794B
program_10_09.m 778B
program_16_13.m 770B
program_09_05.m 767B
program_16_17.m 709B
program_10_03.m 696B
program_11_15.m 675B
program_13_20.m 672B
program_14_09.m 666B
program_09_06.m 656B
program_17_03_SPECIAL.m 650B
program_18_05_SPECIAL.m 641B
program_18_04_SPECIAL.m 636B
program_10_06.m 633B
program_13_10.m 627B
program_13_09.m 627B
program_11_07.m 617B
program_13_08.m 616B
program_13_03.m 610B
program_10_14.m 604B
program_17_02_SPECIAL.m 602B
program_11_13.m 590B
program_16_23.m 588B
program_13_22.m 581B
program_01_03.m 576B
program_14_11.m 572B
program_13_19.m 571B
program_13_15.m 567B
program_17_01_SPECIAL.m 559B
program_13_18.m 547B
program_13_24.m 546B
program_13_07.m 527B
program_10_12.m 522B
program_10_15.m 520B
program_14_10.m 518B
program_13_12.m 498B
program_13_02.m 496B
program_15_03.m 490B
program_15_02.m 490B
program_14_08.m 489B
program_14_04.m 440B
program_10_10.m 436B
program_13_25.m 433B
program_13_06.m 407B
program_14_07.m 404B
program_11_02.m 403B
program_16_16.m 403B
program_01_01.m 398B
program_11_03.m 396B
program_09_01.m 393B
program_16_18.m 383B
program_09_14.m 382B
program_09_20.m 376B
program_10_05.m 376B
program_14_06.m 372B
program_13_04.m 370B
program_13_21.m 370B
program_09_26.m 365B
program_13_13.m 365B
program_12_10.m 353B
program_09_13.m 353B
共 153 条
- 1
- 2
资源评论
frogwangwang
- 粉丝: 1
- 资源: 2
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功