基于二维CS-SCHT变换和扩频方法的彩色图像水印嵌入和提取算法matlab仿真
1.算法运行效果图预览
2.算法运行软件版本
matlab2022a
3.算法理论概述
水印嵌入原理
水印提取原理:
将嵌入水印的图像再次进行二维CS-SCHT变换。
提取变换后的低频系数,并按照嵌入时的规则去除宿主图像内容的影响,恢复出水印信息 Wm′。
4.部分核心程序
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 | % figure; % subplot(121); % imshow(Irgb) % subplot(122); % imshow(I1) %RGB转换为RGB2YUV R = double (Irgb(:,:,1)); G = double (Irgb(:,:,2)); B = double (Irgb(:,:,3)); Y = 0.299*R + 0.587*G + 0.114*B; U = -0.147*R - 0.289*G + 0.436*B; V = 0.615*R - 0.515*G - 0.100*B; YUV = cat (3, Y, U, V); [I0,Marks1] = func_read_images(Y,I1); I0 = imresize(I0,[512,512]); %设置嵌入强度 Power = 5; %设置块的大小 Blksize = 8; RR = 60; CC = 60; %对水印进行置乱 Marks1s = Arnold(Marks1,1,0); % figure % subplot(131); % imshow(Irgb,[]); % title('原始图像'); % subplot(132); % imshow(Marks1,[]); % title('水印'); % subplot(133); % imshow(Marks1s,[]); % title('置乱后水印'); %对图像进行处理 [Mwk_1st] = func_wk_insert(I0,Marks1s,Power,Blksize,RR,CC); Y = imresize(Y,[512,512]); U = imresize(U,[512,512]); V = imresize(V,[512,512]); Mwk_1stYUV2 = cat (3,Mwk_1st,U,V); %YUV2RGB RGB1 = zeros ( size (Mwk_1stYUV2)); RGB1(:,:,1) = Mwk_1st + 1.14 * V; RGB1(:,:,2) = Mwk_1st - 0.39 * U - 0.58 * V; RGB1(:,:,3) = Mwk_1st + 2.03 * U; RGB1n(:,:,1)=awgn(RGB1(:,:,1),SNRS(ij), 'measured' ); RGB1n(:,:,2)=awgn(RGB1(:,:,2),SNRS(ij), 'measured' ); RGB1n(:,:,3)=awgn(RGB1(:,:,3),SNRS(ij), 'measured' ); %水印提取 Rs = double (RGB1n(:,:,1)); Gs = double (RGB1n(:,:,2)); Bs = double (RGB1n(:,:,3)); Ys = 0.299*Rs + 0.587*Gs + 0.114*Bs; Us = -0.147*Rs - 0.289*Gs + 0.436*Bs; Vs = 0.615*Rs - 0.515*Gs - 0.100*Bs; Mwk_1st2 = Ys(:,:,1); Msg1 = func_wk_desert( real (Mwk_1st2),Blksize,RR,CC,4); %对水印进行逆置乱 Msg1s = Arnold( uint8 (255*Msg1),1,1); % figure % subplot(231); % imshow(Irgb,[]); % title('原始图像'); % subplot(232); % imshow(uint8(YUV)); % title('RGB转为YUV图片'); % subplot(233); % imshow(Marks1); % title('水印'); % subplot(234); % imshow(uint8(RGB1)); % title('JPEG压缩攻击'); % subplot(236); % imshow(Msg1s,[]); % title('水印提取'); %% %水印提取,计算NC和PSNR Marks1 = imresize(Marks1,[RR,CC]); NC(ij,jk) = func_nc( uint8 (Marks1), uint8 (Msg1s)) ; end end figure ; plot (SNRS, mean (NC,2), 'b-o' ); xlabel ( 'SNR' ); ylabel ( '水印提取NC值' ); |
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 一个费力不讨好的项目,让我损失了近一半的绩效!
· 清华大学推出第四讲使用 DeepSeek + DeepResearch 让科研像聊天一样简单!
· 实操Deepseek接入个人知识库
· CSnakes vs Python.NET:高效嵌入与灵活互通的跨语言方案对比
· Plotly.NET 一个为 .NET 打造的强大开源交互式图表库