数字水印 | Arnold 变换的 Python 代码实现

2024-05-14 03:04

本文主要是介绍数字水印 | Arnold 变换的 Python 代码实现,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

🥭 参考博客: Arnold 阿诺德置乱(猫脸变换)图像盲水印注入预处理(Python)



1 回顾:Arnold 公式

A r n o l d \mathsf{Arnold} Arnold 变换公式如下:

[ x n + 1 y n + 1 ] = [ 1 b a a b + 1 ] [ x n y n ] m o d ( N ) \begin{bmatrix} x_{n+1} \\ y_{n+1} \end{bmatrix}= \begin{bmatrix} 1 & b\\ a & ab+1 \end{bmatrix} \begin{bmatrix} x_{n} \\ y_{n} \end{bmatrix} mod(N) [xn+1yn+1]=[1abab+1][xnyn]mod(N)

其中, x n x_n xn 代表置乱 n n n 次后像素的行坐标, y n y_n yn 代表置乱 n n n 次后像素的列坐标。相应地, x n + 1 x_{n+1} xn+1 y n + 1 y_{n+1} yn+1 是在 x n x_{n} xn y n y_{n} yn 的基础上再置乱一次的结果。

A r n o l d \mathsf{Arnold} Arnold 逆变换公式如下:

[ x n y n ] = [ a b + 1 − b − a 1 ] [ x n + 1 y n + 1 ] m o d ( N ) \begin{bmatrix} x_{n} \\ y_{n} \end{bmatrix}= \begin{bmatrix} ab+1 & -b\\ -a & 1 \end{bmatrix} \begin{bmatrix} x_{n+1} \\ y_{n+1} \end{bmatrix} mod(N) [xnyn]=[ab+1ab1][xn+1yn+1]mod(N)

本质上就是乘置乱矩阵的逆矩阵。

个人理解:所谓的置乱,就是根据公式改变原始图像中每个像素的位置。



2 核心代码实现

A r n o l d \mathsf{Arnold} Arnold 变换的代码实现:

def arnold(img, shuffle_times, a, b):r, c, d = img.shapep = np.zeros(img.shape, np.uint8)for s in range(shuffle_times):for i in range(r):for j in range(c):x = (i + b * j) % ry = (a * i + (a * b + 1) * j) % cp[x, y, :] = img[i, j, :]img = np.copy(p)return p

A r n o l d \mathsf{Arnold} Arnold 逆变换的代码实现:

def de_arnold(img, shuffle_times, a, b):r, c, d = img.shapep = np.zeros(img.shape, np.uint8)for s in range(shuffle_times):for i in range(r):for j in range(c):x = ((a * b + 1) * i - b * j) % ry = (- a * i + j) % cp[x, y, :] = img[i, j, :]img = np.copy(p)return p

两个函数只有最里层的 for 循环体不一样。

参数说明:

  • i m g \mathsf{img} img 是待处理的图像;
  • s h u f f l e _ t i m e s \mathsf{shuffle\_times} shuffle_times 是变换的次数;
  • a , b \mathsf{a,b} a,b 是置乱矩阵的参数,可以自行指定;

代码说明:

x = (i + b * j) % r
y = (a * i + (a * b + 1) * j) % c

上述代码没有像公式中那样统一 m o d ( N ) mod(N) mod(N),而是根据行数 r \mathsf{r} r 和列数 c \mathsf{c} c 分别取余。

p[x, y, :] = img[i, j, :]

在变换得到的图像 P \mathsf{P} P ( x , y ) \mathsf{(x,y)} (x,y) 位置上,是原始图像 i m g \mathsf{img} img ( i , j ) \mathsf{(i,j)} (i,j) 位置上的像素。此外,: 代表 P \mathsf{P} P 复制了 i m g \mathsf{img} img 的所有通道,从而实现了彩色图像的置乱!



3 完整代码实现

import cv2
import numpy as np
from matplotlib import pyplot as pltdef arnold(img, shuffle_times, a, b):r, c, d = img.shapep = np.zeros(img.shape, np.uint8)for s in range(shuffle_times):for i in range(r):for j in range(c):x = (i + b * j) % ry = (a * i + (a * b + 1) * j) % cp[x, y, :] = img[i, j, :]img = np.copy(p)return pdef de_arnold(img, shuffle_times, a, b):r, c, d = img.shapep = np.zeros(img.shape, np.uint8)for s in range(shuffle_times):for i in range(r):for j in range(c):x = ((a * b + 1) * i - b * j) % ry = (- a * i + j) % cp[x, y, :] = img[i, j, :]img = np.copy(p)return pImg_path = 'white_bear.jpg'
Img = cv2.imread(Img_path)
Img = Img[:, :, [2, 1, 0]]Img_arnold = arnold(Img, 5, 2, 3)
Img_inverse_arnold = de_arnold(Img_arnold, 5, 2, 3)plt.subplot(1, 2, 1)
plt.title("arnold", fontsize=12, loc="center")
plt.axis('off')
plt.imshow(Img_arnold)plt.subplot(1, 2, 2)
plt.title("de_arnold", fontsize=12, loc="center")
plt.axis('off')
plt.imshow(Img_inverse_arnold)
plt.show()

实现效果:

在这里插入图片描述



这篇关于数字水印 | Arnold 变换的 Python 代码实现的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/987573

相关文章

java之Objects.nonNull用法代码解读

《java之Objects.nonNull用法代码解读》:本文主要介绍java之Objects.nonNull用法代码,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐... 目录Java之Objects.nonwww.chinasem.cnNull用法代码Objects.nonN

Python如何使用__slots__实现节省内存和性能优化

《Python如何使用__slots__实现节省内存和性能优化》你有想过,一个小小的__slots__能让你的Python类内存消耗直接减半吗,没错,今天咱们要聊的就是这个让人眼前一亮的技巧,感兴趣的... 目录背景:内存吃得满满的类__slots__:你的内存管理小助手举个大概的例子:看看效果如何?1.

Python+PyQt5实现多屏幕协同播放功能

《Python+PyQt5实现多屏幕协同播放功能》在现代会议展示、数字广告、展览展示等场景中,多屏幕协同播放已成为刚需,下面我们就来看看如何利用Python和PyQt5开发一套功能强大的跨屏播控系统吧... 目录一、项目概述:突破传统播放限制二、核心技术解析2.1 多屏管理机制2.2 播放引擎设计2.3 专

Python中随机休眠技术原理与应用详解

《Python中随机休眠技术原理与应用详解》在编程中,让程序暂停执行特定时间是常见需求,当需要引入不确定性时,随机休眠就成为关键技巧,下面我们就来看看Python中随机休眠技术的具体实现与应用吧... 目录引言一、实现原理与基础方法1.1 核心函数解析1.2 基础实现模板1.3 整数版实现二、典型应用场景2

Python实现无痛修改第三方库源码的方法详解

《Python实现无痛修改第三方库源码的方法详解》很多时候,我们下载的第三方库是不会有需求不满足的情况,但也有极少的情况,第三方库没有兼顾到需求,本文将介绍几个修改源码的操作,大家可以根据需求进行选择... 目录需求不符合模拟示例 1. 修改源文件2. 继承修改3. 猴子补丁4. 追踪局部变量需求不符合很

idea中创建新类时自动添加注释的实现

《idea中创建新类时自动添加注释的实现》在每次使用idea创建一个新类时,过了一段时间发现看不懂这个类是用来干嘛的,为了解决这个问题,我们可以设置在创建一个新类时自动添加注释,帮助我们理解这个类的用... 目录前言:详细操作:步骤一:点击上方的 文件(File),点击&nbmyHIgsp;设置(Setti

SpringBoot实现MD5加盐算法的示例代码

《SpringBoot实现MD5加盐算法的示例代码》加盐算法是一种用于增强密码安全性的技术,本文主要介绍了SpringBoot实现MD5加盐算法的示例代码,文中通过示例代码介绍的非常详细,对大家的学习... 目录一、什么是加盐算法二、如何实现加盐算法2.1 加盐算法代码实现2.2 注册页面中进行密码加盐2.

MySQL大表数据的分区与分库分表的实现

《MySQL大表数据的分区与分库分表的实现》数据库的分区和分库分表是两种常用的技术方案,本文主要介绍了MySQL大表数据的分区与分库分表的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有... 目录1. mysql大表数据的分区1.1 什么是分区?1.2 分区的类型1.3 分区的优点1.4 分

一文详解如何从零构建Spring Boot Starter并实现整合

《一文详解如何从零构建SpringBootStarter并实现整合》SpringBoot是一个开源的Java基础框架,用于创建独立、生产级的基于Spring框架的应用程序,:本文主要介绍如何从... 目录一、Spring Boot Starter的核心价值二、Starter项目创建全流程2.1 项目初始化(

Mysql删除几亿条数据表中的部分数据的方法实现

《Mysql删除几亿条数据表中的部分数据的方法实现》在MySQL中删除一个大表中的数据时,需要特别注意操作的性能和对系统的影响,本文主要介绍了Mysql删除几亿条数据表中的部分数据的方法实现,具有一定... 目录1、需求2、方案1. 使用 DELETE 语句分批删除2. 使用 INPLACE ALTER T