Logistic映射在图像加密中的应用

Logistic映射

Logistic映射是一种可产生的非线性系统,模型表示如下:
x n + 1 = μ x n ( 1 − x n ) x_{n+1}=\mu x_n(1-x_n) xn+1=μxn(1xn)
式中: 0 < μ ≤ 4 ; 0 < x < 1 ; n ∈ N ; x n ∈ [ 0 , 1 ] 0<\mu\leq4;0<x<1;n∈N;x_n∈[0,1] 0<μ4;0<x<1;nN;xn[0,1] μ \mu μ是分岔参数, 当 3.5699456... < μ ≤ 4 当3.5699456...<\mu\leq4 3.5699456...<μ4时,映射进入混沌(chaos)区域。

Logistic映射分岔图像如下图所示。

在这里插入图片描述

注:

Logistic模型分叉图生成过程及代码见文章尾部参考链接1

混沌方法的数字图像加密实现

摘要

目前混沌系统与加密技术相结合是现如今最热门的一个课题之一,虽然有大量的加密算法面世,但是这些加密算法并不成熟,仍然需要进一步的研究。而混沌系统具有许多重要的性质,如对初始条件和系统参数的敏感依赖性,使得基于混沌系统的图像加密得到了广泛的研究与应用。

虫口模型-Logistic混沌映射

传统的混沌映射加密可分为置换和扩散两个阶段,像素置乱变换与灰度值替换都由混沌系统产生的混沌序列所控制,扩散过程如下:

  • 步骤1:利用混沌序列获取密钥流。
  • 步骤2:加密图像矩阵的像素值。
  • 步骤3:重复步骤2,直到达到密文图像。

算法

设图像(i,j)处的灰度值为I(i,j),满足 1 ≤ i ≤ M 、 1 ≤ j ≤ N 1\leq i \leq M、1\leq j \leq N 1iM1jN, I ′ ( i , j ) I'(i,j) I(i,j)表示替换后 I ( i , j ) I(i,j) I(i,j)在(i,j)出的灰度值。

加密公式
I ′ ( i , j ) = ( ( ( r 1 ( i , j ) + r 2 ( i , j ) ) ⊕ r 3 ( i , j ) ) + I ( i , j ) ) m o d 256 I'(i,j)=(((r_1(i,j)+r_2(i,j))\oplus r_3(i,j))+I(i,j))mod256 I(i,j)=(((r1(i,j)+r2(i,j))r3(i,j))+I(i,j))mod256
其中:mod表示求模运算; ⊕ \oplus 表示按位异或运算。r1,r2,r3表示的是混沌序列值与255的乘积(采用多混沌序列产生的密钥空间大于单一的混沌序列所产生的密钥空间),替换变换的密钥由r1,r2,r3对应的混沌系统提供,可以进行多次变换寻求更好加密效果。

解密公式:
I ( i , j ) = ( I ′ ( i , j ) − ( ( r 1 ( i , j ) + r 2 ( i , j ) ) ⊕ r 3 ( i , j ) ) ) m o d 256 I(i,j)=(I'(i,j)-((r_1(i,j)+r_2(i,j))\oplus r_3(i,j))) mod 256 I(i,j)=(I(i,j)((r1(i,j)+r2(i,j))r3(i,j)))mod256
解密是加密的逆,可以容易得到验证。

实验环境

  • 语言:python3.8
  • 编辑器:Sublime Text 3
  • python库:numpy、matplotlib、cv2
  • 实验图像下载

实验代码(python3.8)

import cv2
import numpy as np
import matplotlib.pyplot as plt

'''
加密函数
img_encrypt:原始图像(二维ndarray)
key:密钥列表,大小为7(1、2、3为混沌系统初始条件;4、5、6为分岔参数u,7为重复加密次数)
return:返回加密后的图像
'''
def encrypt(img,key):
  #图像的宽高
  [w,h]=img.shape
  #混沌系统初始条件
  x1=key[0]
  x2=key[1]
  x3=key[2]
  #分岔参数u
  u1=key[3]
  u2=key[4]
  u3=key[5]
  #加密次数
  n=key[6]
  #一个临时数组用于返回加密后的图像,可以不影响原始图像
  img_tmp=np.zeros((w,h))
  #对原始图像的每个像素都处理n次
  for k in range(n):
    for i in range(w):
      for j in range(h):
        #计算混沌序列值
        x1=u1*x1*(1-x1)
        x2=u2*x2*(1-x2)
        x3=u3*x3*(1-x3)
        ##混沌值位于[0,1]区间内,所以可以看做是一个系数,乘以最大灰度值并转成整数用于异或运算即可
        r1=int(x1*255)
        r2=int(x2*255)
        r3=int(x3*255)
        img_tmp[i][j]=(((r1+r2)^r3)+img[i][j])%256
    #下一轮加密重新初始化混沌系统
    x1=key[0]
    x2=key[1]
    x3=key[3]
  return img_tmp


'''
解密函数
img:加密后图像(二维ndarray)
key:密钥列表,大小为7(1、2、3为混沌系统初始条件;4、5、6为分岔参数u,7为重复加密次数)
return:返回解密后的图像
'''
def decrypt(img,key):
    #图像的宽高
  [w,h]=img.shape
  #混沌系统初始条件
  x1=key[0]
  x2=key[1]
  x3=key[2]
  #分岔参数u
  u1=key[3]
  u2=key[4]
  u3=key[5]
  #加密次数
  n=key[6]
  #一个临时数组用于返回加密后的图像,可以不影响传入的加密图像
  img_tmp=np.zeros((w,h))
  #对原始图像的每个像素都处理n次
  for k in range(n):
    for i in range(w):
      for j in range(h):
        x1=u1*x1*(1-x1)
        x2=u2*x2*(1-x2)
        x3=u3*x3*(1-x3)
        r1=int(x1*255)
        r2=int(x2*255)
        r3=int(x3*255)
        img_tmp[i][j]=(img[i][j]-((r1+r2)^r3))%256
        # img[i][j]=(img[i][j]-((r1[i*w+j]+r2[i*w+j])^r3[i*w+j]))%256
    #下一轮加密重新初始化混沌系统
    x1=key[0]
    x2=key[1]
    x3=key[3]
  return img_tmp


def main():
  #原始图像路径
  path='./lena.jpg'
  #加密密钥参数列表
  key=[0.343,0.432,0.63,3.769,3.82,3.85,1]
  #读取原始图像
  img=cv2.imread(path)
  #img.shape=>(512,512,3),宽、高、通道数量

  #图像灰度化
  img_gray = cv2.cvtColor(img,cv2.COLOR_RGB2GRAY)
  #混沌加密原始图像
  img_encrypt=encrypt(img_gray,key)
  #错误密钥解密图像
  wrong_key=[0.342,0.432,0.61,3.769,3.82,3.85,1]
  wrong_decrypt=decrypt(img_encrypt,wrong_key)
  #正确密钥解密图像
  img_decrypt=decrypt(img_encrypt,key)

  #结果展示
  plt.rcParams['font.sans-serif'] = ['SimHei']  # 中文乱码
  #子图1,原始图像
  plt.subplot(221)
  #plt默认使用三通道显示图像,所以需要制定cmap参数为gray
  #imshow()对图像进行处理,画出图像,show()进行图像显示
  plt.imshow(img_gray,cmap='gray')
  plt.title('原始图像(灰度化)')
  #不显示坐标轴
  plt.axis('off')

  #子图2,加密后图像
  plt.subplot(222)
  plt.imshow(img_encrypt,cmap='gray')
  plt.title('加密图像(密钥{}'.format(key))
  plt.axis('off')

  #子图3,错误密钥解密结果
  plt.subplot(223)
  plt.imshow(wrong_decrypt,cmap='gray')
  plt.title('解密图像(密钥{})'.format(wrong_key))
  plt.axis('off')

  #子图4,正确密钥解密结果
  plt.subplot(224)
  plt.imshow(img_decrypt,cmap='gray')
  plt.title('解密图像(密钥{})'.format(key))
  plt.axis('off')

  #设置子图默认的间距
  # plt.tight_layout()
  #保存图像
  # cv2.imwrite('./lean_original.jpg',img)
  # cv2.imwrite('./lean_gray.jpg',img_gray)
  # cv2.imwrite('./lean_encrypt.jpg',img_encrypt)
  # cv2.imwrite('./lean_gray.jpg',wrong_decrypt)
  # cv2.imwrite('./lean_decrypt.jpg',img_decrypt)

  #显示图像
  plt.show()


if __name__ == '__main__':
  main()

实验结果

在这里插入图片描述

从图中可以看出,哪怕错误密钥与正确密钥只是微小的差别,也难以从解密结果中看出一丝原始图像的信息。采用混沌系统来加密图像取得了较好的效果。

参考:

  • 9
    点赞
  • 63
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

z2bns

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值