【图像加密】Arnold置乱图像加密解密【含GUI Matlab源码 1239期】

最新推荐文章于 2024-05-21 16:30:09 发布

Matlab武动乾坤

最新推荐文章于 2024-05-21 16:30:09 发布

阅读量536

点赞数 6

分类专栏： Matlab图像处理（进阶版）文章标签： matlab

本文链接：https://blog.csdn.net/KeepingMatlab/article/details/136048409

版权

Matlab图像处理（进阶版）专栏收录该内容

833 篇文章 15 订阅

订阅专栏

✅博主简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，Matlab项目合作可私信。
🍎个人主页：海神之光
🏆代码获取方式：
海神之光Matlab王者学习之路—代码获取方式
⛳️座右铭：行百里者，半于九十。

⛄一、Arnold置乱图像加密解密简介

0 前言
网络已经成为我们传递信息的主要平台, 为我们提供诸多便捷的同时, 也存在一些安全问题, 特别是一些重要信息的传递.如果在信息传递前先对其进行加密, 能够在一定程度上保护所传递的信息.数字图像作为重要的信息资源在人们的生活中发挥着越来越重要的作用[1], 因此, 数字图像的加密是一项值得研究的重要课题.本文介绍的就是一种基于Arnold变换的图像加解密算法。

1 Arnold变换
将N*N图像上的点 (x, y) 通过如下变换转成 (x′, y′) 如式 (1) , 该变换即称为Arnold变换.通过变换公式可发现, 其变换的本质是点的位置的变换, 并且这种变换保证变换前后的点保持一一对应的关系。
在这里插入图片描述
如果将一次变换的输出作为下一次变换的输入, 这就是迭代变换如式 (2) .当一次变换置乱效果不佳时, 往往需要迭代变换获得更好的置乱效果。

2 基于Arnold变换的图像加解密
图像加密也称图像置乱, 是对图像的像素进行混乱和扩散, 使加密后的图像在视觉上无法获得有效信息.空域加密是常用的方法, 分为空域置乱和序列加密.空域置乱是对像素坐标进行变换使其混乱, 解密时恢复原像素坐标.图像的加密既可以作为独立的信息隐藏方法, 也可以用来作为数字水印技术中图像水印的预处理。

加密图像可以是灰度图像也可以是彩色图像, 如果是灰度图像则需要将各点的灰度值带到新的坐标点, 如果是彩色图像则需要将各点的RGB值带到新的坐标点.本文所加密的图像为彩色图像, 因此首先保存原图点 (x, y) 处的RGB值, 然后使用Arnold变换改变点 (x, y) 位置到 (x′, y′) , 同时将原图点 (x, y) 处的RGB值带到 (x′, y′) 处.我们采用256*256的彩色图像作为原始图像如图1a所示, 在Visual Studio2008平台下运用式 (1) 对原图进行一次Arnold变换, 变换后的效果如图1b所示.

由图1可知, 对图像进行一次Arnold变换后仍会留下较多信息, 故在此基础上继续使用Arnold变换, 通过5次变换后 (见图2b) 其结果在视觉上就很难再获得什么有价值的信息, 在视觉上起到了图像加密的作用.同时变换的次数往往可以作为密钥, 从而在一定程度上提高安全性。
在这里插入图片描述
图1 原图和一次Arnold变换结果

图2 第2次和第5次Arnold变换结果
然而此时的加密并不完全可靠, 若已知采用Arnold变换作为加密方式, 则通过暴力求解法, 经过若干次的变换还是可以解密出原图.因此在一般情况下, 我们往往需要在Arnold变换中加入密钥 (ku, kv) , 以提高安全性, 其公式如下:
在这里插入图片描述
当需要对图像进行解密时, 则根据加密后的图像中的像素点与原图中的像素点是一一对应的关系进行解密, 因此图像解密的过程为:根据原图中像素点的坐标 (x, y) 按照Arnold变换算出新图中像素点的坐标 (x′, y′) , 获取新像素点 (x′, y′) 的RGB值, 将该值归还给 (x, y) .根据加密过程使用的变换次数, 解密时则要将以上过程重复相应的次数就可以恢复原图 (见图3) 。

3 基于三维Arnold变换的图像加密
前面提到的图像置乱是通过打乱像素点位置而达成目标, 图像的颜色分量值是保持不变的.因此为了提高加密效果, 我们可以在像素点位置已乱的前提下, 再通过改变RGB颜色分量的值对图像进行加密.采用如下三维Arnold变换 (如式4) 进行RGB分量值的改变。
在这里插入图片描述
我们将一次Arnold变换后的图像作为原始图像, 在此基础上对各像素点的RGB采用以上公式进行扰乱, 获得如图4所示结果.通过两种加密方式的综合, 则进一步提高了加密图像的安全性.式4中的值可以做相应的改变, 同时可以将变换的次数作为密钥提高安全性。
在这里插入图片描述
图3 第五次Arnold变换结果, 恢复后的加密图像

图4 第一次Arnold变换结果, RGB扰乱后的加密图像

⛄二、部分源代码

%函数功能—★图像加密与解密处理总界面★
%--------------------------------------------------------------------
clear all;
button_color=[0.529 0.808 0.922];
%创建一个窗口
h=figure(‘Color’,[0 0.4 0],‘MenuBar’,‘none’,‘name’,…
‘图像加密软件’,‘NumberTitle’,‘off’,…
‘Position’,[420 150 600 500]);
%'MenuBar’的值figure表示显示控制菜单none表示隐藏
%'NumberTitle’的值’off’标题栏中没有窗口序列号on表示有
%--------------------------------------------------------------------

%--------------------------------------------------------------------
%界面背景图案默认设置
h_b=axes(‘units’,‘normalized’,‘position’,[0 0 1 1]);
uistack(h_b,‘down’)
Q_bg=imread(‘bg.jpg’);
image(Q_bg);
colormap gray
set(h_b,‘handlevisibility’,‘off’,‘visible’,‘off’);
%--------------------------------------------------------------------

%--------------------------------------------------------------------
%菜单
%一级菜单
f1=uimenu(h,‘Label’,‘&File 文件’);
f2=uimenu(h,‘Label’,‘&Edit 编辑’);
f3=uimenu(h,‘Label’,‘图像处理’);
f4=uimenu(h,‘Label’,‘图像加密’);
f5=uimenu(h,‘Label’,‘&help 帮助’);

%二级菜单
m11=uimenu(f1,‘Label’,‘打开’,‘Accelerator’,‘O’,‘Callback’,‘picture_open’);
m12=uimenu(f1,‘Label’,‘保存’,‘Accelerator’,‘S’,‘Callback’,‘picture_save’);
m13=uimenu(f1,‘Label’,‘图像格式转换’);
m14=uimenu(f1,‘Label’,‘关闭所有窗口’,‘Callback’,‘close all;’);
m15=uimenu(f1,‘Label’,‘退出’,‘Accelerator’,‘Q’,‘Callback’,‘exit’,‘Separator’,‘On’);

% m21=uimenu(f2,‘Label’,‘改变按钮颜色’,‘Callback’,‘button_color=uisetcolor’);
m22=uimenu(f2,‘Label’,‘背景更换’,‘Callback’,‘bg_select;’);
m23=uimenu(f2,‘Label’,‘图像裁剪’,‘Callback’,‘imshow(Q_show);Q_imcrop=imcrop;Q_global=Q_imcrop;imshow(Q_imcrop);’);

m31=uimenu(f3,‘label’,‘图像旋转’,‘Callback’,‘picture_rotate;’);
m32=uimenu(f3,‘label’,‘图像灰化’,‘Callback’,‘picture_rgbtogray;’);
m33=uimenu(f3,‘label’,‘图像二值化’,‘Callback’,‘picture_graytotwo;’);
m34=uimenu(f3,‘label’,‘图像马赛克’,‘Callback’,‘picture_mosaic;’);

m41=uimenu(f4,‘label’,‘Arnold_jms’,‘Callback’,‘bu_trigger1;’);
m42=uimenu(f4,‘label’,‘Arnold_chaos’,‘Callback’,‘bu_trigger2;’);
%函数功能—★Arnold_jms加密/解密界面★
%--------------------------------------------------------------------
button_color=[0.8 0.68 0.604];
%创建一个窗口
h1=figure(‘Color’,button_color,‘MenuBar’,‘none’,‘name’,…
‘★Arnold_jms加密/解密★’,‘NumberTitle’,‘off’,…
‘Position’,[220 150 600 500]);
%'MenuBar’的值figure表示显示控制菜单none表示隐藏
%'NumberTitle’的值’off’标题栏中没有窗口序列号on表示有
axes(‘Position’,[.0 .3 .6 .6]);
Q_open=Q_global;
imshow(Q_open);title(‘待处理图像’); %在窗口中显示待处理的图像
%--------------------------------------------------------------------

%--------------------------------------------------------------------
%界面背景图案默认设置
h_b=axes(‘units’,‘normalized’,‘position’,[0 0 1 1]);
uistack(h_b,‘down’)
Q_bg=imread(‘bg1.jpg’);
image(Q_bg);
colormap gray
set(h_b,‘handlevisibility’,‘off’,‘visible’,‘off’);
%--------------------------------------------------------------------

%--------------------------------------------------------------------
%密钥输入窗口
EditBox=uicontrol(h1,‘Style’,‘edit’,‘Position’,[475 320 50 20],…
‘Callback’,…
[‘NumStr=get(EditBox,’‘String’‘);’…
‘Num=str2num(NumStr);’] );
Textpos=uicontrol(h1,‘Style’,‘text’,‘Position’,…
[400 320 60 20],‘FontSize’,12,‘String’,‘密钥 r’,…
‘BackgroundColor’,[1 0.894 0.710],‘FontSize’,13);
%按扭设置
button=uicontrol(h1,‘Style’,‘pushbutton’,‘String’,‘保存图像’,‘Position’,…
[425 150 100 40],‘BackgroundColor’,button_color,…
‘Callback’,‘picture_save;’);

button_bg=uicontrol(h1,‘Style’,‘pushbutton’,‘String’,‘背景更换’,‘Position’,…
[425 225 100 40],‘BackgroundColor’,button_color,…
‘Callback’,‘bg_select;’);

% Q：加密图像
%alignment：加密序列
%m：输入密钥
%a,b,c：图像三维
%加密k=0;解密k=1;
%u取3.5699456到4之间的数
%lq取100以上的数

function Q_global=chaos(Image,k,u,lq,mm)
if u>4 || u<3.5699456
msgbox(‘u 必须为3.5699456到4之间的数！’,‘错误警告’,‘warn’);
return;
end
if lq<100
msgbox(‘lq必须大于100！’,‘错误警告’,‘warn’);
return;
end
if mm>1 || mm<0
msgbox(‘mm必须大于0且小于1！’,‘错误警告’,‘warn’);
return;
end

if length(size(Image))2%灰度图像
if k0
%加密
ssq=lq;%取ssq位后的混沌序列
xs=u;
Q=Image;
[a,b]=size(Q);
N=a*b;
alignment=zeros(1,N);
alignment(1)=mm;

for yu=1:ssq
alignment(1)=xsalignment(1)-xsalignment(1)^2;
end
for k=1:N-1
alignment(k+1)=xsalignment(k)-xsalignment(k)^2;
end
alignment=mod(10^13*alignment,256);
alignment=uint8(alignment);
n=1;
for ii=1:a
for jj=1:b
Q(ii,jj)=bitxor(alignment(n),Q(ii,jj));
n=n+1;
end
end
Qtwochaos_jm=Q;
Q_global=Qtwochaos_jm;
imshow(Qtwochaos_jm);
title(‘已加密图像’);
msgbox(‘加密成功！’,‘提示’,‘warn’);

%解密
else
ssq=lq;%取ssq位后的混沌序列
xs=u;
Qcrypt=Image;
[aa,bb]=size(Qcrypt);
NN=aabb;
alignment=zeros(1,NN);
alignment(1)=mm;
for yu=1:ssq
alignment(1)=xsalignment(1)-xsalignment(1)^2;
end
for k=1:NN-1
alignment(k+1)=xsalignment(k)-xsalignment(k)^2;
end
alignment=mod(10^13alignment,256);
alignment=uint8(alignment);
nn=1;
for ii=1:aa
for jj=1:bb
Qcrypt(ii,jj)=bitxor(alignment(nn),Qcrypt(ii,jj));
nn=nn+1;
end
end
Qtwochaos_de=Qcrypt;
Q_global=Qtwochaos_de;
imshow(Qtwochaos_de);
title(‘已解密图像’);
msgbox(‘解密成功！’,‘提示’,‘warn’);
end

else %rgb
if k==0
%加密
ssq=lq;%取ssq位后的混沌序列
xs=u;
Q=Image;
[a,b,~]=size(Q);
N=a*b;
alignment=zeros(1,N);
alignment(1)=mm;

for yu=1:ssq
alignment(1)=xsalignment(1)-xsalignment(1)^2;
end
for k=1:N-1
alignment(k+1)=xsalignment(k)-xsalignment(k)^2;
end
alignment=mod(10^13*alignment,256);
alignment=uint8(alignment);
n=1;
for ii=1:a
for jj=1:b
Q(ii,jj,1)=bitxor(alignment(n),Q(ii,jj,1));
Q(ii,jj,2)=bitxor(alignment(n),Q(ii,jj,2));
Q(ii,jj,3)=bitxor(alignment(n),Q(ii,jj,3));
n=n+1;
end
end
Qrgbchaos_jm=Q;
Q_global=Qrgbchaos_jm;
imshow(Qrgbchaos_jm);
title(‘已加密图像’);
msgbox(‘加密成功！’,‘提示’,‘warn’);

%解密
else
ssq=lq;%取ssq位后的混沌序列
xs=u;
Qcrypt=Image;
[aa,bb,~]=size(Qcrypt);
NN=aabb;
alignment=zeros(1,NN);
alignment(1)=mm;
for yu=1:ssq
alignment(1)=xsalignment(1)-xsalignment(1)^2;
end
for k=1:NN-1
alignment(k+1)=xsalignment(k)-xsalignment(k)^2;
end
alignment=mod(10^13alignment,256);
alignment=uint8(alignment);
nn=1;
for ii=1:aa
for jj=1:bb
Qcrypt(ii,jj,1)=bitxor(alignment(nn),Qcrypt(ii,jj,1));
Qcrypt(ii,jj,2)=bitxor(alignment(nn),Qcrypt(ii,jj,2));
Qcrypt(ii,jj,3)=bitxor(alignment(nn),Qcrypt(ii,jj,3));
nn=nn+1;
end
end
Qrgbchaos_de=Qcrypt;
Q_global=Qrgbchaos_de;
imshow(Qrgbchaos_de);
title(‘已解密图像’);
msgbox(‘解密成功！’,‘提示’,‘warn’);
end

end
%--------------------------------------------------------------------

⛄三、运行结果

在这里插入图片描述

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]张义,宛楠.一种基于Arnold变换的数字图像加解密算法[J].安徽工程大学学报. 2013,28(03)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流人群疏散病毒扩散晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合