CN103886542A - 量子Arnold图像置乱方法 - Google Patents

量子Arnold图像置乱方法 Download PDF

Info

Publication number
CN103886542A
CN103886542A CN201410109511.7A CN201410109511A CN103886542A CN 103886542 A CN103886542 A CN 103886542A CN 201410109511 A CN201410109511 A CN 201410109511A CN 103886542 A CN103886542 A CN 103886542A
Authority
CN
China
Prior art keywords
quantum
arnold
image
scramble
totalizer
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410109511.7A
Other languages
English (en)
Other versions
CN103886542B (zh
Inventor
姜楠
吴文亚
王珞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Kangao Medical Technology Co ltd
Original Assignee
Beijing University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201410109511.7A priority Critical patent/CN103886542B/zh
Publication of CN103886542A publication Critical patent/CN103886542A/zh
Application granted granted Critical
Publication of CN103886542B publication Critical patent/CN103886542B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Two-Way Televisions, Distribution Of Moving Picture Or The Like (AREA)

Abstract

本发明涉及量子Arnold图像置乱方法,其特征在于能够在量子计算机上实现图像的Arnold置乱。所述方法包括:给定原始量子图像,对原始图像的x轴坐标进行Arnold置乱变换,对原始图像的y轴坐标进行Arnold置乱变换。对x轴坐标的置乱变换还包括计算x+y和对2n取模;对y轴坐标的置乱变换还包括计算2y、计算x+2y和对2n取模。本发明算法以量子线路的形式给出,能够保证算法能在量子计算机上直接运行;与已有的量子Arnold置乱算法相比,本发明的取模操作通过直接忽略进位实现,大大降低了网络复杂度。

Description

量子Arnold图像置乱方法
技术领域
本发明涉及量子计算技术和图像置乱技术,涉及一种在量子计算机上实现Arnold图像置乱效果的方法。
背景技术
1.Arnold图像置乱
图像置乱是指通过移动像素的位置将图像变为不可读的形式,以达到保护图像信息的目的。为了保证置乱的置乱度,置乱按照某种规律进行。按照规律的不同,有多种置乱方法,Arnold置乱是其中之一。
在经典计算机上,数字图像可以看作是一个二维矩阵I(x,y)。假设图像尺寸为2n×2n,n为非负整数,置乱后图像为I′(x′,y′),则Arnold置乱可定义为:
x′=(x+y)mod2n         (1)
y′=(x+2y)mod2n
其中,mod为取模运算。
置乱后,像素从(x,y)位置移动到了(x′,y′)位置,即原来(x,y)位置存储的像素值被拷贝到(x′,y′)位置。
2.量子计算机和量子逻辑门
量子计算机利用量子力学的规律进行信息处理和逻辑操作,是不同于经典计算机的新型计算设备,有着经典计算机无可比拟的优势:具有与生俱来的并行计算能力,并能避免能量的额外损耗。
所有量子计算机上的信息处理和逻辑操作,即量子算法,均可以用量子线路的形式来表示。量子线路由量子逻辑门组成。量子逻辑门的输入和输出是量子比特,量子比特也是量子计算机上信息存储和处理的最小单位,用符号|·>表示。量子线路的本质是一个流程图,告诉量子计算机怎样经过一步一步的处理,将量子比特从输入态变为输出态,每一步处理对应一个量子逻辑门。量子线路又可以看成一种量子装置,可以将输入态变为输出态。
常见量子逻辑门及其符号表示和功能如表1所示。所有逻辑门默认左端是输入,右端是输出。
表1常见量子逻辑门
Figure BDA0000480791660000011
Figure BDA0000480791660000021
3.量子图像
量子计算机上存储和处理的图像是量子图像。量子图像的表示方法完全不同于经典计算机上的数字图像。一个2n×2n的图像用2n+1个量子比特表示,其中n个量子比特存储X方向的位置信息,另外n个量子比特存储Y方向的位置信息,1个量子比特存储灰度信息。即:
| I > = 1 2 n Σ i = 0 2 2 n - 1 | c i > ⊗ | i >
| c i > = cos θ i | 0 > + sin θ i | 0 > , θ i ∈ [ 0 , π 2 ] , i = 0,1 , . . . , 2 2 n - 1 - - - ( 2 )
|i>=|y>|x>=|yn-1yn-2…y0>|xn-1xn-2…x0>,|yj>|xj>∈{0,1},j=n-1,…,0
其中,|yn-1yn-2…y0>|xn-1xn-2…x0>是位置信息,|ci>是灰度信息,
Figure BDA0000480791660000024
是克罗内克积。
4.量子加法器
量子加法器是在量子计算机上实现两数相加的装置。假设|a>和|b>是两个量子数,均由n个量子比特构成,|a>=|an-1an-2…a0>,|b>=|bn-1bn-2…b0>,ai,bi∈{0,1},则量子加法器可以算出a和b的和,计算结果保留在|can-1an-2…a0>中。其中|c>,c∈{0,1},是初值为0的辅助位,用来存储进位。
5.量子Arnold图像置乱方法与经典Arnold图像置乱方法的不同
本发明要实现的量子Arnold图像置乱方法,在目的上与经典Arnold图像置乱方法是一致的,即要达到同样的置乱效果。但是两者又存在非常多的不同:
(1)适用设备不一样:经典方法在经典计算机上运行;量子方法在量子计算机上运行。相互之间无法通用。
(2)图像表示方法不一样:经典计算机上一个2n×2n图像,需要一个2n×2n数组来存储,如果数组中每个元素用8比特表示,即像素色深为8,则整个图像需要占用8×2n×2n比特,如果色深为24,RGB各占8比特,则整个图像需要占用24×2n×2n比特,空间复杂度为O(22n);量子计算机上,仅需2n+1个量子比特,空间复杂度为O(n)。
(3)像素移动原理不一样:经典计算机中进行像素移动时,真的需要将像素值从一个位置拷贝到另一个位置;量子计算机中是通过调整位置信息|yn-1yn-2…y0>和|xn-1xn-2…x0>中各个量子比特的值来实现,比如,如果量子线路能够把位置信息|11>变为|00>,那么原来|11>位置上的像素就被搬移到了|00>位置。
6.已有量子Arnold置乱方法
姜楠等发表在Quantum Information Processing上的论文“The quantum realization of Arnoldand Fibonacci image scrambling”,公开了一种量子Arnold置乱方法,该方法基于量子加法器实现,但对取模运算的处理过于复杂。
发明内容
本发明的目的是给出一种在量子计算机上实现Arnold图像置乱效果的方法。
为了实现上述目的,本发明采用以下技术方案。
图1所示为总体流程图。从公式(1)能够看出,x′和y′分别用两个公式计算,彼此之间不存在交叉。因此给出量子算法时,x′和y′可分别用两个量子线路实现。
图2所示为实现x′的流程图。根据公式(1),实现x′时,先将x和y相加,再实现取模运算。
图3所示为x′的量子实现方法。为实现x和y相加,采用量子加法器。将位置信息|yn-1yn-2…y0>和|xn-1xn-2…x0>作为加法器的输入,最下面的量子比特|0>用来存储进位。
为实现取模操作,取结果时,直接将进位忽略即可。这是因为两个n比特的二进制数相加,其和小于2n+1,所以直接舍弃进位位即是其和对2n取模的值。
线路图中的CARRY和SUM是两个量子计算模块,模块中的黑色竖条表明模块中量子逻辑门的摆放顺序。如果黑色竖条在右边,则模块的左边是输入,右边是输出;如果黑色竖条在左边,则模块的右边是输入,左边是输出,即将模块中的逻辑门顺序颠倒使用。
图4所示为上图中两个量子计算模块的分解。
图5所示为实现y′的流程图。根据公式(1),实现y′时,先将y乘以2,再将x和2y相加,最后实现取模运算。
图6所示为y′的量子实现方法。为实现y乘以2,将y中的所有量子比特左移1位,最右端补0,并忽略最高位,即将|yn-1yn-2…y0>变为|yn-2…y00>。左移1位并在最右端补0,是为了乘以2;忽略最高位是因为后面要进行取模运算,而(x+2y)mod2n≡(x+((2y)mod2n))mod2n,因此可以在此处先忽略掉一次最高位。在量子计算机中,将|yn-1yn-2…y0>变为|yn-2…y00>不需要用量子逻辑门处理,只需要将量子比特|yn-2…y00>输入到相应的位置即可。
为实现x和2y相加,采用量子加法器。将位置信息|yn-2…y00>和xn-1xn-2…x0>作为加法器的输入。
为实现取模操作,取结果时,直接将进位忽略即可。
与现有技术相比,本发明具有以下优点:
(1)本发明适用于量子计算机,能够在量子计算机上实现图像的Arnold置乱。
(2)本发明算法以量子线路的形式给出,保证算法能在量子计算机上直接运行。
(3)与已有的量子Arnold置乱算法相比,本发明的取模操作通过直接忽略进位实现,网络复杂度大大降低。
附图说明
图1为本发明所涉及方法的总体流程图;
图2为实现x′的流程图;
图3为x′的量子实现原理图;
图4为量子计算模块分解示意图;
图5为实现y′的流程图;
图6为y′的量子实现示意图;
图7为1个4×4图像的示意图;
图8为具体实施例中x′的量子实现原理图;
图9为具体实施例中y′的量子实现原理图。
具体实施方式
下面结合附图和实施例对本发明做进一步说明。
本发明所述的Arnold置乱方法的流程图如图1所示,包括以下步骤:
步骤1,给定原始量子图像|I>,其表达式为:
| I > = 1 2 n Σ i = 0 2 2 n - 1 | c i > ⊗ | i >
| c i > = cos θ i | 0 > + sin θ i | 0 > , θ i ∈ [ 0 , π 2 ] , i = 0,1 , . . . , 2 2 n - 1
|i>=|y>|x>=|yn-1yn-2…y0>|xn-1xn-2…x0>,|yj>|xj>∈{0,1},j=n-1,…,0
其中,2n×2n为经典图像尺寸,|yn-1yn-2…y0>|xn-1xn-2…x0>是位置信息,|ci>是灰度信息,
Figure BDA0000480791660000052
是克罗内克积。
步骤2,对原始图像的x轴坐标进行Arnold置乱变换。
(1)计算x+y。
采用量子加法器。将原始图像的位置信息|yn-1yn-2…y0>和|xn-1xn-2…x0>作为加法器的输入,加法器最下面的量子比特|0>用来存储进位。
(2)对2n取模。
因为两个n比特的二进制数相加,其和小于2n+1,所以直接舍弃进位位即是其和对2n取模的值。
Arnold置乱后的x轴坐标依然存储在量子比特|xn-1xn-2…x0>中,各个量子比特的值由xn-1xn-2…x0变为x′n-1x′n-2…x′0
步骤3,对原始图像的y轴坐标进行Arnold置乱变换。
(1)计算2y。
将y中的所有量子比特左移1位,最右端补0,并忽略最高位,即将|yn-1yn-2…y0>变为|yn-2…y00>。
(2)计算x+2y。
采用量子加法器。将左移1位后的y轴位置信息|yn-2…y00>和原始图像的x轴位置信息|xn-1xn-2…x0>作为加法器的输入,加法器最下面的量子比特0用来存储进位。
(3)对2n取模。
为实现取模操作,取结果时,直接将进位忽略即可。
Arnold置乱后的y轴坐标依然存储在量子比特|yn-1yn-2…y0>中,各个量子比特的值由yn-1yn-2…y0变为y′n-1y′n-2…y′0
其中步骤2和步骤3顺序可颠倒,也可并行执行。经过步骤2和步骤3给出的量子线路的处理,量子图像的位置信息由量子态|y>|x>,即|yx>,变为量子态|y′x′>,也就是将量子图像中处在|yx>位置的像素移动到|y′x′>位置,即完成了Arnold图像置乱。
下面给出一个对4×4像素大小的图片进行Arnold置乱的实例。
图7所示为1个4×4图像,这里n=2,即位置信息x和y各用2个量子比特表示,(a)是原图,(b)是Arnold置乱后的效果。
图8所示为具体实施例对应的实现x′的量子线路。由于n=2,其输入为位置信息|y1y0>、|x1x0>以及3个初值为0的辅助位。线路中每部分与图3的对应关系,在图中用虚线标出。
在该实施例中,注意到最下面的辅助位不做控制位,即该量子位对其他量子位没有影响,且它是最终被忽略的进位,因此该量子位对本实施例来讲是多余的,可从线路中删除,与之有关联从左往右数的第4个和第6个逻辑门随之也被删除。这样第5个、第7个、第8个逻辑门就成为有同样输入的并排在一起的控制非门1。3个控制非门1并列的作用相当于1个控制非门,可以将其中的2个删除。这样实现x′的量子线路被等价化简为等价符号后面的样子。
图9所示为具体实施例对应的实现y′的量子线路。该线路与实现x′的量子线路结构是完全相同的,差别仅在于输入不同。
根据图8和图9所示的量子线路,以及表1给出的各量子逻辑门的作用,可以得到Arnold置乱真值表,如表2所示。根据真值表,Arnold置乱后的图像为图7(b)。比如像素“E”原来所处的位置在第01行、00列,置乱后,该像素处在第10行、01列。
表2对4×4图片进行Arnold置乱处理的真值表
Figure BDA0000480791660000061

Claims (2)

1.量子Arnold图像置乱方法,其特征在于能够在量子计算机上实现图像的Arnold置乱;所述方法包括以下步骤:
步骤1,给定原始量子图像|I>,其表达式为:
| I > = 1 2 n Σ i = 0 2 2 n - 1 | c i > ⊗ | i >
| c i > = cos θ i | 0 > + sin θ i | 0 > , θ i ∈ [ 0 , π 2 ] , i = 0,1 , . . . , 2 2 n - 1
|i>=|y>|x>=|yn-1yn-2…y0>|xn-1xn-2…x0>,|yj>|xj>∈{0,1},j=n-1,…,0
其中,2n×2n为经典图像尺寸,|yn-1yn-2…y0>|xn-1xn-2…x0>是位置信息,|ci>是灰度信息,
Figure FDA0000480791650000013
是克罗内克积;
步骤2,对原始图像的x轴坐标进行Arnold置乱变换;
(1)计算x+y;
采用量子加法器;将原始图像的位置信息|yn-1yn-2…y0>和|xn-1xn-2…x0>作为加法器的输入,加法器最下面的量子比特|0>用来存储进位;
(2)对2n取模;
步骤(1)计算得到的x+y对2n取模;
取模后完成对原始图像x轴坐标的Arnold置乱,置乱后的x轴坐标依然存储在量子比特|xn-1xn-2…x0>中,各个量子比特的值由xn-1xn-2…x0变为x′n-1x′n-2…x′0
步骤3,对原始图像的y轴坐标进行Arnold置乱变换;
(1)计算2y;
将y中的所有量子比特左移1位,最右端补0,并忽略最高位,即将>yn-1yn-2…y0>变为|yn-2…y00>;
(2)计算x+2y;
采用量子加法器;将左移1位后的y轴位置信息|yn-2…y00>和原始图像的x轴位置信息|xn-1xn-2…x0>作为加法器的输入,加法器最下面的量子比特|0>用来存储进位;
(3)对2n取模;
步骤(2)计算得到的x+2y对2n取模;
取模后完成对原始图像y轴坐标的Arnold置乱,置乱后的y轴坐标依然存储在量子比特|yn-1yn-2…y0>中,各个量子比特的值由yn-1yn-2…y0变为y′n-1y′n-2…y′0
2.根据权利要求1所述的量子Arnold图像置乱方法,其特征在于,所述步骤2和步骤3所述的对2n取模运算通过在取结果时直接将进位忽略实现。
CN201410109511.7A 2014-03-24 2014-03-24 量子Arnold图像置乱方法 Active CN103886542B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410109511.7A CN103886542B (zh) 2014-03-24 2014-03-24 量子Arnold图像置乱方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410109511.7A CN103886542B (zh) 2014-03-24 2014-03-24 量子Arnold图像置乱方法

Publications (2)

Publication Number Publication Date
CN103886542A true CN103886542A (zh) 2014-06-25
CN103886542B CN103886542B (zh) 2017-01-25

Family

ID=50955419

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410109511.7A Active CN103886542B (zh) 2014-03-24 2014-03-24 量子Arnold图像置乱方法

Country Status (1)

Country Link
CN (1) CN103886542B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104217391A (zh) * 2014-08-26 2014-12-17 长春理工大学 一种多通道量子图像的加密方法
CN106454382A (zh) * 2016-09-20 2017-02-22 北京工业大学 一种量子图像制备方法
CN107832851A (zh) * 2017-10-19 2018-03-23 南京邮电大学 一种基于级联ghz态编码的逻辑量子比特量子门构造方法
CN109741236A (zh) * 2018-12-12 2019-05-10 四川大学 一种在ibm量子实验平台实现的量子图像阈值分割方法
CN110473134A (zh) * 2019-07-18 2019-11-19 广西师范大学 一种基于gneqr的量子图像置乱方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6190276A (ja) * 1984-10-09 1986-05-08 Yokogawa Hokushin Electric Corp 距離変換回路
CN101894359A (zh) * 2010-07-13 2010-11-24 北京信息科技大学 图像置乱加密方法
CN102833448A (zh) * 2012-07-28 2012-12-19 辽宁大学 一种基于线性链表的图像置乱方法
CN103617590A (zh) * 2013-11-28 2014-03-05 辽宁大学 一种基于魏尔斯特拉斯函数的图像置乱方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6190276A (ja) * 1984-10-09 1986-05-08 Yokogawa Hokushin Electric Corp 距離変換回路
CN101894359A (zh) * 2010-07-13 2010-11-24 北京信息科技大学 图像置乱加密方法
CN102833448A (zh) * 2012-07-28 2012-12-19 辽宁大学 一种基于线性链表的图像置乱方法
CN103617590A (zh) * 2013-11-28 2014-03-05 辽宁大学 一种基于魏尔斯特拉斯函数的图像置乱方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
宋国强: "数字图像加密技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
曹光辉等: "基于Logistic均匀分布图像置乱方法", 《物理学报》 *
郭文文: "基于经典变换域的数字图像水印算法的研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104217391A (zh) * 2014-08-26 2014-12-17 长春理工大学 一种多通道量子图像的加密方法
CN104217391B (zh) * 2014-08-26 2017-07-18 长春理工大学 一种多通道量子图像的加密方法
CN106454382A (zh) * 2016-09-20 2017-02-22 北京工业大学 一种量子图像制备方法
CN106454382B (zh) * 2016-09-20 2019-04-30 北京工业大学 一种量子图像制备方法
CN107832851A (zh) * 2017-10-19 2018-03-23 南京邮电大学 一种基于级联ghz态编码的逻辑量子比特量子门构造方法
CN107832851B (zh) * 2017-10-19 2020-02-14 南京邮电大学 一种基于级联ghz态编码的逻辑量子比特量子门构造方法
CN109741236A (zh) * 2018-12-12 2019-05-10 四川大学 一种在ibm量子实验平台实现的量子图像阈值分割方法
CN109741236B (zh) * 2018-12-12 2022-11-11 四川大学 一种在ibm量子实验平台实现的量子图像阈值分割方法
CN110473134A (zh) * 2019-07-18 2019-11-19 广西师范大学 一种基于gneqr的量子图像置乱方法
CN110473134B (zh) * 2019-07-18 2021-02-19 广西师范大学 一种基于gneqr的量子图像置乱方法

Also Published As

Publication number Publication date
CN103886542B (zh) 2017-01-25

Similar Documents

Publication Publication Date Title
CN103886542A (zh) 量子Arnold图像置乱方法
Guo et al. Threshold visual secret sharing by random grids with improved contrast
US9450938B1 (en) Information secure proactive multiparty computation (PMPC) protocol with linear bandwidth complexity
CN107579813A (zh) 信息加密、解密方法及装置
CN104737495A (zh) 校正量子密钥分配协议中的比特误差的排列方法
Liu et al. Image hiding scheme by use of rotating squared sub-image in the gyrator transform domains
Wang et al. A new chaotic image encryption algorithm based on L-shaped method of dynamic block
Song et al. Multi-image reorganization encryption based on SLF cascade chaos and bit scrambling
CN106127669A (zh) 基于保面积Baker映射的混沌图像加密方法
CN108416445B (zh) 一种量子实信号的存储与量子线路实现的设计方法
CN103368725B (zh) 一种g0类s盒构造方法及其电路
CN103780794A (zh) 一种基于混沌系统的图像加密改进方法
Nagpaul Topics in applied abstract algebra
CN103077498A (zh) 基于移位变换的图像加密算法
Sang et al. Quantum realization of Arnold scrambling for IFRQI
CN103236858A (zh) 基于循环左移的cmmb中准循环ldpc串行编码器
US20240187209A1 (en) Data encryption and decryption using screens and lfsr-generated logic blocks
Gates et al. N= 4 and N= 8 SUSY Quantum Mechanics and Klein's Vierergruppe
Huang et al. Design of Image Encryption Algorithm Based on Compound Two-dimensional Maps.
CN109450618A (zh) 一种基于md5的加密方法及系统
RU2485592C1 (ru) Способ формирования целочисленных неортогональных декоррелирующих матриц заданных размеров и устройство для его осуществления
CN103929207A (zh) Cmmb中部分并行输入的右移累加qc-ldpc编码器
Fowler et al. Z= p Z= p actions on Sn Sn
CN104363014A (zh) 一种基于四值逻辑电路的数据译码方法
CN103905055A (zh) 部分并行输入的右移累加qc-ldpc编码器

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20210127

Address after: Room 1424, Floor 4, Peony Pioneer Building, No. 2 Garden Road, Haidian District, Beijing, 100191

Patentee after: Beijing Zhonglian Technology Service Co.,Ltd.

Address before: 100124 No. 100 Chaoyang District Ping Tian Park, Beijing

Patentee before: Beijing University of Technology

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220302

Address after: 200000 room 5g11, No. 1, Lane 129, Datian Road, Jing'an District, Shanghai

Patentee after: Shanghai Huashi Investment Co.,Ltd.

Address before: Room 1424, Floor 4, Peony Pioneer Building, No. 2 Garden Road, Haidian District, Beijing, 100191

Patentee before: Beijing Zhonglian Technology Service Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240109

Address after: Room 1260, South Building of Film and Television Park, No. 4915 Beisong Road, Songjiang District, Shanghai, 200000

Patentee after: SHANGHAI KANGAO MEDICAL TECHNOLOGY Co.,Ltd.

Address before: 200000 room 5g11, No. 1, Lane 129, Datian Road, Jing'an District, Shanghai

Patentee before: Shanghai Huashi Investment Co.,Ltd.

TR01 Transfer of patent right