【图像加密】RSA图像加密解密【含Matlab源码 1442期】

Matlab武动乾坤

于 2024-01-24 12:38:56 发布

阅读量1k

点赞数 23

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

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

版权

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

797 篇文章 5 订阅

订阅专栏

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

⛄一、RSA图像加密简介

图像数据的获取、传输和处理已经遍及了数字时代的各个角落,图像的安全问题日益严重,尤其是在军事、商业和医疗等特殊领域。RSA是一种目前被公开应用的加密算法,它采用非对称的加密体制。下面,笔者在研究图像加密算法的基础上,将RSA算法应用到图像加密技术中,重点研究了算法优化、密钥生成、密钥分发和文件传输等内容,提出了一套完整的数字图像加密和解密解决方案。

1 RSA基本原理
RSA算法采用一种非对称密码加密体制,在整个加密过程中密钥和算法独立分开,使得密钥更能有效的等到分配。其采用大素数因子分解的难度来保障安全性,目前,大素数分解问题仍然没有很好的解决方法。
RSA算法加密的基本步骤如下:
Step1选取不同的2个大素数a,b,然后计算2个数乘积m=a×b。
Step2选取大整数加密密钥d,要求d与(a-1)×(b-1)互质。
Step3求解密密钥e,e×d=1mod(a-1)×(b-1)。
Step4将明文P加密为密文C(C=Pdmod m),密文C解密为明文P(P=Cemod m)。
根据以上步骤可以看出,只知道m和d不能计算出解密密钥e,可见,任何人都可以对明文P进行加密,但只有获得授权许可后才能对密文C进行解密。

2 RSA图像加密流程设计
RSA图像加密流程设计是先对图像进行读取并转化成十六进制数据流;然后生成RSA算法所需密钥,将密钥与图像进行幂乘及取模运算,生成十进制数据;最后将数据转换为字符串数据流进行保存。其中,最重要的是密钥的生成[ ,其决定最后图像加密效果。密钥的生成过程包括自动生成大素数并存储、对大素数进行算术运算、大数幂模与乘模运算和素数的自动生成4部分,具体加密流程如图1所示。

3 RSA图像加密功能设计
根据RSA图像加密流程,设计图像加密功能,主要包括RSA图像加密、解密和密钥生成等功能。系统功能模块如图2所示。
1)数字图像加密。对数字图像进行字节流的读取并转换为十六进制流,应用RSA算法对十六进制流进行加密,将加密后的数据转化为文本输出。
2)数字图像解密。加载加密后图像文件,利用密钥对其进行解密,对加密的图像进行还原。
在这里插入图片描述
图1 RSA图像加密流程图

图2 系统功能模块图
3)加密解密预览。在数字图像加密解密过程中,确保待加密及解密后的数字图像可视化。
4)密钥自动生成。RSA算法需要用户自定义输入2个大素数,主要实现自动生成大素数及密钥以减少用户操作,并能确保密钥使用的素数足够大。
5)密钥长度设定。用户可根据不同情况设定密钥的长度,能够灵活控制加密和解密速度。
6)密钥文件导出。以文件形式导出密钥文件,确保传输或存储过程的安全性。
7)密钥文件导入。对导出的密钥文件进行读入,具有识别功能。可以识别本软件导出的密钥文件,并对内容检查,确保导入文件安全。
8)密钥输入。解密时支持手动输入加密处理的密钥串,确保解密操作安全。
9)文件打印。实现对密钥文件的打印,可以更好的保存私钥。
10)文件传输。实现加密后的文件及图像直接传输,方便快捷且能更好的确保安全。

4 RSA图像加密关键技术实现
4.1 RSA密钥的大数存储
RSA算法生成密钥需要2个大素数为21024或更大比较安全[7] 。但编程语言中unsigned int类型最多存储2个字节,远远小于RSA安全密钥长度。笔者设计单元线性数组实现大素数存储,解决编程中大素数的存储问题。首先设置一个以unsigned为单元的线性数组用来存储大素数,定义2个无符号整数z和n来控制存储单元数,z是分配空间的单元数,如果大素数的长度超过unsigned数组的预定义数组长度,z会随着数字变大不断增大;n表示当前存储大素数已占用单元数。每个大素数最大可以达到232*占用单元数满足RSA的各种运算。

4.2 RSA密钥的大数运算
由于生成的大素数超过21024,原有数据运算方式不再适用。在大数存储基础上采用类间派生与关联方式实现大数运算。定义flex_unit派生得到vlong_value类实现新的运算函数,将原vlong_value类关联到新类vlong中,在新vlong类实现运算符重载。本类运算是按一定数制对数字的计算,乘除和取余也都按照竖式运算的原理实现。加运算的核心代码如下:
在这里插入图片描述
4.3 大素数幂模与乘模运算优化
幂模运算是RSA算法中比重最大的计算,其主要决定生成最后的公钥和私钥,直接地决定了RSA算法的性能。依据乘模的性质,把幂模运算转换成乘模运算,实现思想是指数不断的对分,具体流程如图3所示。

乘模运算能够提高运算速度,通常对于千位以上的二进制整数n,利用普通除法求模运算速度很慢[9] 。笔者采用Montgomery算法实现幂乘运算的优化来提高加密速度,改进思想是先选取与模数m互素的基数Y(Y=2k),m为奇数满足2k-1≤m<2k,再选取Y-1(0<Y-1<m)和m′(0<m′<m),使得YY-1-mm′=1。然后选取任意整数n(0≤n<Ym),最后实现求模乘法nY-1mod m的快速算法,来提高加密运算速度。

4.4 算法中素数自动生成及优化
笔者采用Eratosthenes筛选法对素数筛选进行优化,优化策略是对各小素数因子求模,得到当前a在素数搜索范围内的最小倍数在b[]中的对应位置,继续后移a个位置,直到将a在搜索范围内的所有倍数全部找到;在完成对所有小素数因子的类似操作后,其倍数在搜索范围内的位置标记b[r]被全部标记为0。其流程如图4所示。
在这里插入图片描述
图3 幂模转换乘模运算流程图

图4 素数搜索除去小素数因子倍数流程图
笔者采用费马小定理对素数测试进行优化,化策略是选取一个整数P,要求与a互素,并满足关系式为Pa-1mod a=1,输入大整数a满足关系式可能不是素数,需要改变P完成多次测试,测试后这个数很大概率为素数。

⛄二、部分源代码

clc
clear all
close all

%% Initialization
Test=0; % for right key encryption
% Test=1; % for key sensitivity test
addpath functions
Images_Path=‘Images\Orignal’;
fname={‘Lena’,‘baboon’,‘cameraman’,‘panda’,‘contact_lens’,‘checkerboard’,‘Black1’,‘Onion’,‘Football’}; % filename
ext=‘.jpg’;
fid=1; % file ID 1 for lena
IS =256; % Image size
Data=imread(strcat(Images_Path,fname{fid},ext));
if (size(Data,3)==3)
Data=rgb2gray(Data);
end
Data=imresize(Data,[IS IS]); % Image Size

[row,col]=size(Data);
[Data,padding]=Scalling(Data,8);
Data_binary=convert2bin(Data);

hex_key = ‘AAAAAAAAAAAAAAAA’;
[bin_key] = Hex2Bin( hex_key );
[K1,K2,K3,K4,K5]=SF_Key_Gen(bin_key);

orignal_msg=[];
encrypt_msg=[];
decrypt_msg=[];

%% Encryption Process
for kk=1:2
for i=1:size(Data_binary,1)
orignal=Data_binary(i,:);
tic
[cipher]=SF_Encrypt(orignal,K1,K2,K3,K4,K5);
encryption_time(i)=toc;
tK1=[K1(1:8),orignal(1:8)];tK2=orignal(9:24);tK3=orignal(25:40);tK4=orignal(41:56);tK5=[orignal(57:64),K5(9:16)];
K1=tK1;K2=tK2;K3=tK3;K4=tK4;K5=tK5;

encrypt_msg(:,i)=Binary2Dec(cipher);
cipher_data(i,:)=double(cipher);
if(kk<2)
Data_binary(i,:)=cipher_data(i,:);
end

end
function O = P(I)
if(Ilogical([0 0 0 0]))
O = logical([0 0 1 1]); end
if(Ilogical([0 0 0 1]))
O = logical([1 1 1 1]); end
if(Ilogical([0 0 1 0]))
O = logical([1 1 1 0]); end
if(Ilogical([0 0 1 1]))
O = logical([0 0 0 0]); end
if(Ilogical([0 1 0 0]))
O = logical([0 1 0 1]); end
if(Ilogical([0 1 0 1]))
O = logical([0 1 0 0]); end
if(Ilogical([0 1 1 0]))
O = logical([1 0 1 1]); end
if(Ilogical([0 1 1 1]))
O = logical([1 1 0 0]); end
if(Ilogical([1 0 0 0]))
O = logical([1 1 0 1]); end
if(Ilogical([1 0 0 1]))
O = logical([1 0 1 0]); end
if(Ilogical([1 0 1 0]))
O = logical([1 0 0 1]); end
if(Ilogical([1 0 1 1]))
O = logical([0 1 1 0]); end
if(Ilogical([1 1 0 0]))
O = logical([0 1 1 1]); end
if(Ilogical([1 1 0 1]))
O = logical([1 0 0 0]); end
if(Ilogical([1 1 1 0]))
O = logical([0 0 1 0]); end
if(Ilogical([1 1 1 1]))
O = logical([0 0 0 1]); end
end
function O = Q(I)
if(Ilogical([0 0 0 0]))
O = logical([1 0 0 1]); end
if(Ilogical([0 0 0 1]))
O = logical([1 1 1 0]); end
if(Ilogical([0 0 1 0]))
O = logical([0 1 0 1]); end
if(Ilogical([0 0 1 1]))
O = logical([0 1 1 0]); end
if(Ilogical([0 1 0 0]))
O = logical([1 0 1 0]); end
if(Ilogical([0 1 0 1]))
O = logical([0 0 1 0]); end
if(Ilogical([0 1 1 0]))
O = logical([0 0 1 1]); end
if(Ilogical([0 1 1 1]))
O = logical([1 1 0 0]); end
if(Ilogical([1 0 0 0]))
O = logical([1 1 1 1]); end
if(Ilogical([1 0 0 1]))
O = logical([0 0 0 0]); end
if(Ilogical([1 0 1 0]))
O = logical([0 1 0 0]); end
if(Ilogical([1 0 1 1]))
O = logical([1 1 0 1]); end
if(Ilogical([1 1 0 0]))
O = logical([0 1 1 1]); end
if(Ilogical([1 1 0 1]))
O = logical([1 0 1 1]); end
if(Ilogical([1 1 1 0]))
O = logical([0 0 0 1]); end
if(Ilogical([1 1 1 1]))
O = logical([1 0 0 0]); end
end

⛄三、运行结果

在这里插入图片描述

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]高新成,王莉利,李苏龙.基于RSA算法的图像加密系统设计与实现[J].长江大学学报(自科版). 2016,13(25)

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 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合