一种基于分段线性混沌映射的医学图像加密算法_秦秋霞.pdf

第 25 卷第 1 期2023 年 1 月大 连 民 族 大 学 学 报Journal of Dalian Minzu UniversityVol25，No1January 2023收稿日期:20221002;最后修回日期:20221122基金项目:国家自然科学基金资助项目(61672121)。作者简介:秦秋霞(1996)，女，山东济南人，大连民族大学计算机科学与工程学院硕士研究生，主要从事图像加密研究。文章编号:20961383(2023)01005707一种基于分段线性混沌映射的医学图像加密算法秦秋霞a，梁仲月a，徐毅b(大连民族大学 a计算机科学与工程学院;b理学院 预科教育学院，辽宁 大连 116650)摘要:为了加强对医学图像信息的保护，提出一种基于分段线性混沌映射的医学图像加密算法，采用多轮加密方式，提高加密算法的安全性。首先对原始图像进行随机数插入操作，使用分段线性混沌映射对随机数插入后的矩阵进行像素值替换，实现置乱，然后将替换后的矩阵与未经过替换前的矩阵进行按位异或操作，实现扩散;最后将得到的矩阵逆时针旋转 90 度，继续进行下一轮加密，经过四轮加密到最终的加密图像。经过实验仿真和安全性分析表明:本算法各项安全性能优越，可以有效抵抗噪声攻击等攻击方式，具有一定的鲁棒性，适用于医学图像加密。关键词:医学图像加密;分段线性混沌映射;按位异或中图分类号:TP312文献标志码:AA Medical Image Encryption Algorithm Based on Piecewise LinearChaotic MappingQIN Qiuxiaa，LIANG Zhongyuea，XU Yib(aSchool of Computer Science and Engineering;bSchool of Science School of Preuniversity，Dalian Minzu University，Dalian Liaoning 116650，China)Abstract:In order to strengthen the protection of medical image information，this paper proposesa medical image encryption algorithm based on piecewise linear chaotic mapping This article u-ses multiple rounds of encryption to improve the security of the encryption algorithm First，arandom number is inserted to the original image After that，piecewise linear chaotic mapping isused to replace the pixel values of the matrix after the random number has been inserted to a-chieve scrambling Then，a bitwise XO operation is performed on the matrix after the replace-ment and the matrix before the replacement to achieve diffusion Finally，the obtained matrix isrotated 90 degrees counter clockwise，and the next round of encryption is continued After fourrounds of encryption，the final encrypted image is obtained，and the encryption process is com-pleted Experimental simulation and security analysis show that this algorithm has superior secur-ity performance，can effectively resist noise attacks and other attacks，has a certain degree of ro-bustness，and is suitable for medical image encryptionKey words:medical image encryption;piecewise linear chaotic mapping;bitwise XO随着医疗水平高速发展，人口老龄化的急剧增长。医学图像信息得到大范围的使用，面对海量、冗余度高的医学图像信息，寻找更加安全、有效的医学图像信息加密算法是十分重要的13。DOI:10.13744/21-1431/g4.2023.01.010为了保证医学图像信息传输的安全性，防止出现个人、相关机构及国家医疗信息的泄露，本文提出了一种基于分段线性混沌映射的医学图像加密算法，通过使用多轮加密的方式，实现了医学图像的加密，提高了医学图像加密算法的安全性和鲁棒性。混沌系统一般指的是一种非常复杂、并且不能依据初始状态用来预测将来状态的具有很大随机性的系统45。同时，混沌系统敏感性的特征十分符合加密算法的要求，适用于扩散和置乱过程。纪等人6 提出了一种基于混沌的新型舰船安全图像加密算法，使用混沌置乱实现置乱操作，提高加密安全性;乔等人7 为了解决传统加密步骤繁琐的问题，提出了一种基于二维 Logistic 混沌系统的图像加密算法，使用混沌序列改变像素值的位置，该算法加密效果良好;张等人8 提出了结合S 盒与混沌映射的图像加密算法，使用二维 Logis-tic 映射和 Chen 混沌系统生成的混沌序列参与置乱、扩散过程，提高了加密效果;杜等人9 提出了基于混沌压缩感知和 DNA 编码的多图像加密算法，混沌系统使用 Lorenz 超混沌系统和 Chen 超混沌系统，产生随机性更强的混沌序列，实现了一次一密的加密方式。从仿真结果和安全性分析结果可以看出，将混沌系统引入图像加密算法可以使得加密效果更加优越，加密算法的随机性和健壮性更有保证。本文使用的混沌映射是 PWLCM 混沌映射，PWLCM 混沌映射是一种分段线性混沌映射，它具有良好的遍历性和很少的周期窗口等特点1012。混沌现象是发生在确定系统中的不确定行为。使用分段线性混沌映射可以实现医学图像中稳定状态和不稳定状态之间的转换。PWLCM 混沌映射的初值敏感性符合加密算法的扩散操作，而随机特性符合加密算法的置乱操作，因此分段线性混沌系统十分符合加密算法的要求。由于 PWLCM混沌映射在加密算法中具有较高的应用效果，因此本文将其引入到医学图像加密算法中。医学图像加密算法不仅可以对普通图像进行加密，还可以对医学图像进行加密。医学图像加密需要尽可能的保证加密过程的安全性，保护病患的隐私信息;另外，还需要保证解密后的得到的图像可以达到医生诊断的正常需求。依据上述问题以及分析，本文提出了一种基于分段线性混沌映射的医学图像加密算法。本算法采用四轮加密的方式，增强算法的安全性，使用分段线性混沌映射参与每一轮的置乱和扩散过程，提高算法的随机性。每一轮的置乱过程都采取随机数插入的方式进行，扩散过程采用按位异或的方式来实现，每一轮结束后进行矩阵旋转然后再继续进行下一轮加密。当四轮加密结束后，得到加密图像，完成整个加密过程。本文剩余章节的结构如下:第一部分主要介绍本算法使用的基本理论;第二部分介绍加、解密算法的流程和详细步骤;第三部分给出加密算法的仿真实验和安全性分析结果;第四部分给出结论。1基本原理本章主要介绍本文用到的分段线性混沌映射，医学图像相关原理。11分段线性混沌映射分段线性混沌映射(PWLCM)是一种常见的分段混沌映射。由于其具有周期窗口小、遍历性好、生成混沌序列分布均匀等特性，研究人员将其应用于图像加密领域1011，13。PWLCM 混沌映射定义如下:xn+1=F(xn，r)=xn/r，0 xn r(xn r)/(05 r)，r xn 05F(1 xn，r)，05xn 1。(1)式中:x 是系统参数;r 是控制参数，x (0，1)。当 r (0，05)时，系统处于混沌状态。PWLCM混沌映射具有一致的变量分布和较好的遍历性，可以提供满足加密算法随机性的随机序列，并且简单容易实现，适合用于图像加密领域。12医学图像医学图像是医护人员进行医疗诊断工作的主要参考数据之一，它是通过医疗器材获取的人体某一部位的图像信息1415。医学图像信息能够直观地反映人体某一部位是否存在病变，而且医学图像在不同的场景下应用时对图像分辨率有着不同的要求。当医生使用医学图像用来判断病人的身体状况时，此时要求医学图像必须具有高分辨率;当医学图像用于病历分析时，可以容忍医学图像存在一定程度的损失压缩，以达到节省存储空间的作用。但不论哪种情况，在使用医学图像时，必须保证医学图像具有较高的分辨率，并且在传输过程中要做到足够的安全性，以避免病人隐私的泄露。在临床应用中，医学图像主要有为 X射线图85大连民族大学学报第 25 卷像、CT、核磁共振成像等类型1617 如图 1。a)Brainb)Chestc)Tibial图 1不同类型的医学图像展示2加解密算法本文提出了一种基于分段线性混沌映射的医学图像加密算法。首先，对原始图像进行随机数的插入。其次，使用分段线性混沌映射替换图像矩阵的像素值，实现置乱操作。然后，将替换后的矩阵与未进行替换前的矩阵进行按位异或操作，实现像素值的改变。最后，将第一轮结束后的矩阵逆时针旋转 90 度，继续进行下一轮加密操作，经过四轮加密后生成加密图像，加密结束。多轮加密和分段线性映射的使用确保加密算法的安全性，可用于医学图像加密。具体的加密步骤如2123 节所示。加密流程图如图 2。图 2加密流程图21随机数插入在原始图像 P 中插入一列随机值，为后续替换和扩散操作做准备，具体如下:第一步:读取原始图像 P，原始图像的长是M，宽是 N。第二步:采用双层循环的方式生成随机数插入后的矩阵 ，的大小是 M (N+1)。其中，矩阵的第 1 列元素是通过随机数函数 randi()生成的，其他列的元素是原始图像的各个像素值。具体如下:(i，j)=randi(i)，j=1P(i，j 1)，j 1。i=1，2，M(2)其中，randi()函数用于生成均匀分布的伪随机整数。22第一轮加密本加密算法总共采用四轮加密，每一轮加密都是类似的，下面以第一轮加密为例进行详细的介绍。首先，对随机数插入后的矩阵 使用分段线性映射进行替换操作，然后将替换后的矩阵进行按位异或操作，最后将得到的矩阵进行矩阵的旋转，继续进行下一轮的加密，直到第四轮加密结束，完成加密。221替换操作整个替换过程分为两个部分，先对第 1 列的元素进行替换，再对其他元素进行替换，替换完成后，生成 x1矩阵。具体如下:第一步:对随机数插入后生成的 矩阵的第1 列进行替换操作。对于 矩阵的第 1 列，第1 行对应的元素，分段线性混沌映射使用已知的控制变量 r0，初值 x0生成的随机值将 矩阵的元素进行替换，作为 x1矩阵的第 1 列，第 1 行的元素。第二步:对于矩阵 的第 1 列的其他元素(i，j)，分段线性混沌映射将已知的 r0作为控制参数，生成的 x1(i 1，1)作为初值，生成随机值将 矩阵的元素进行替换。第三步:对于矩阵 的第 1 列之外的元素(i，j)，分段线性混沌映射将 r1作为控制参数，x1(i，j 1)作为初值生成随机值用来替代 矩阵中的元素。具体如下:x1(i，j)=pwlcm(r0，x0);i=1，j=1pwlcm(r0，x1(i 1，1);i 1，j=1pwlcm(r1，x1(i，j 1);i 1，j 1。(3)其中，pwlcm()函数是分段线性混沌映射函数。222扩散操作将替换操作后生成的 x1矩阵与未进行替换前生成的 矩阵进行按位异或操作，实现像素值的改变，具体如下:第一步:将 221 节生成的 x1矩阵进行处理，95第 1 期秦秋霞，等:一种基于分段线性混沌映射的医学图像加密算法让其值保持在 0，255 之间，以便进行按位异或操作，具体如下:x=floor(mod(x1(i，j)1015，256);i=1，2，M;j=1，2，N+1;(4)其中，floor()是向下取整函数，mod()是取余函数。第二步:将 x 矩阵、矩阵进行按位异或操作，生成大小为 M (N+1)的矩阵 B。具体如下:B(i，j)=B(i，j 1)(i，j)x;i=1，2，M;j=2，N+1;(5)其中，是按位异或操作。23加密算法具体的加密步骤如下:(1)原始图像 P 插入随机数生成矩阵 ，提高了置乱和扩散的效果。如 21 节所示。(2)分段线性混沌映射使用不同的控制参数和初值生成随机数，使用随机数将矩阵 对应的各个元素进行替换，生成矩阵 x1。如 221 节所示。(3)将 x1矩阵以及 矩阵进行按位异或操作，生成矩阵 B，实现扩散过程。如 222 节所示。(4)将矩阵 B 以逆时针的方式旋转 90 度，继续进行第二轮加密过程，每一轮与上述的(1)、(2)、(3)步类似。(5)当第四轮加密结束后，生成加密图像 C，完成加密过程。24解密算法解密算法是加密算法的逆过程，所以不再进行详细的介绍。在加密算法的每轮加密过程中，先进行随机数插入，其次进行混沌值替换，然后进行按位异或扩散，最后进行逆时针旋转 90 度操作。在解密算法中，要注意各个操作的顺序，以及各个操作对应的逆操作的准确性，加解密算法对于各种操作都特别敏感，存在一点错误都是无法正确解密的。3安全性分析使用个人电脑对三组不同形状、不同类型的的医学图像以及一组普通图像进行加密、解密仿真实验。具体环境:win10 系统，MATLAB 2020a软件。其中，本实验中使用的具体的加密密钥见表 1。675901 像素的 Brain 医学图像的原图、加密图以及解密图如图3。256256 像素的 Chest 医学图像的原图、加密图以及解密图如图 4。512512 像素的 Tibial 医学图像的原图、加密图以及解密图如图 5。512512 大小的 Lena 图像的原图、加密图以及解密图如图 6。表 1密钥表类型密钥PWLCM 的初值x0=047PWLCM 的控制参数r0=052r1=031r2=026r3=074r4=043a)Brain 原图b)Brain 密图c)Brain 解密图图 3Brain 医学图像的加密、解密效果a)Chest 原图b)Chest 密图c)Chest 解密图图 4Chest 医学图像的加密、解密效果a)Tibial 原图b)Tibial 密图c)Tibial 解密图图 5Tibial 医学图像的加密、解密效果b)Lena 原图b)Lena 密图c)Lena 解密图图 6Lena 图的加密、解密效果通过图 36 可以看出，本加密算法不仅能对医学图像进行加密、解密，而且也能够对普通图像进行加密、解密。从实验效果可以看出，使用本加06大连民族大学学报第 25 卷密算法不管是医学图像、还是普通图像得到的加密图是完全不同于原始图像的，丝毫看不出原始图像的任何痕迹，说明加密算法能够达到隐藏原始图像信息的目的。通过原始图像和解密图像的效果来看，从肉眼上看不出两者之间有任何区别，说明从直观上可以分析出本算法能够完全满足解密效果。31密钥空间通常，当加密算法的密钥长度大于 2100，就能够抵抗穷举攻击1819。本文提出的算法中，安全性是由四轮加密保障的。每一轮加密所使用的密钥都是 PWLCM 混沌映射1011 给定的的初值和控制参数。从密钥表中可以看出，一共有 6 个密钥值。每一个密钥值在计算机中都是 64 位，可以表示成具有 15 位有效数字的十进制数，所以本加密算法的密钥空间是(1015)6=1090，该密钥空间远远大于2100 20，所以说明本算法可以有效抵抗暴力攻击，保证算法的安全性。32信息熵信息熵是用来表征随机性的。信息熵可以作为衡量像素分布随机性的标准，信息熵的值越大，说明随机性越强2122。它的数学表达式定义如下:H(m)=Li=0P(mi)log2P(mi)。(6)式中:m 表示图像;H(m)表示该图像计算得出的信息熵;mi表示图像 m 的第 i 个灰度值;P(mi)表示 mi出现的概率;L 表示图像 m 的灰度级。当图像的灰度级为 8 时，图像的最大信息熵值为 8，也就是理想的随机图像的信息熵为 823。医学图像 Chest 的原图的信息熵和经过加密算法后得到密图的信息熵情况，另外还展示了本算法与几个最新算法的信息熵的对比情况，具体如见表 2。表 2原图、密图信息熵对比表信息熵对比原图信息熵密图信息熵本算法6494 87997 4文献20 7568 37996 8文献24 7486 57965 5文献25 7465 17997 3通过表 2 可以看出，原始图像的信息熵值与理想值 826 相差较大，说明随机性比较低。而经过本加密算法得到的信息熵值为 7997 4，与理想值十分接近，说明经过本算法得到的加密图像随机性很强，抗攻击的能力也强，说明本加密算法的性能比较好。同时为了更加客观的说明本算法的信息熵值，对比了现有文献的信息熵值，对比结果说明本算法的信息熵值优于现有算法的信息熵值，进一步说明本算法能够达到现有算法的安全级别。33相关性原始图像相邻像素间存在非常高的相关性，可以通过一个像素就能反映出相邻像素的信息27。所以对于图像加密来说最重要的是破坏相邻像素间这种强相关性，进而提高抵抗统计攻击能力。为了验证加密算法可以打破原始图像相邻像素之间的相关性，并且拥有抵抗统计攻击的能力，通过以下公式计算相关性:rxy=1NNi=1(xi E(x)(yi E(y)D(x)D(y);(7)E(x)=1NNi=1xi;(8)D(x)=1NNi=1(xi E(x)2)。(9)式中:x 和 y 都是图像的相邻像素值;rxy是计算得出的相邻像素点的相关性系数;E(x)是期望值;D(x)是方差。加密算法应该有效减低图像相关性，理想情况下，加密图的相关性为 028。本实验对 Chest 医学图像的原图、加密图、解密图分别随机选取 1000 对相邻像素进行相关系系数计算，另外与最新算法的相关性系数进行对比，具体情况见表 3。从表 3 可以看出，原图的三个方向的相关系系数都较大，说明相关性特性明显，被统计性攻击破解的风险极大。经过加密后得到的加密图的三个方向的相关系系数接近于 0，也就是说明加密图的像素值之间几乎不存在关系。本算法与最新算法的相关性比较，通过数值看出相关性基本上低于最新算法的相关性。因此，可以看出本加密算法降低图像相关性的能力较强。16第 1 期秦秋霞，等:一种基于分段线性混沌映射的医学图像加密算法表 3原图、密图相关性系数对比表相关性系数原图密图水平垂直对角线水平垂直对角线本算法0999 00999 30998 30006 00009 80005 3文献 200920 80955 90881 90005 70012 70007 2文献 240999 00999 30998 30004 20006 10016 7文献 290974 90954 60917 50011 20007 10010 2另外，本实验还通过观察图像的相关性散点图的方式展示相关性情况。其中，Chest 原图的水平、垂直、对角线方向的相关性情况如图 7。Chest加密图的水平、垂直、对角线方向的相关性情况如图 8。a)原图水平b)原图垂直c)原图对角线图 7原图相关性散点图a)密图水平b)密图垂直c)密图对角线图 8密图相关性散点图通过原图和密图的相关性散点图可以看出，原图的大多数相邻像素对集中分布在对角线或者对角线两侧的范围内，而经过本算法加密后的相邻像素对能够均匀、随机地分布在整个坐标系的范围内，说明本加密算法能够有效打破明文像素间的强相关性，并且能够有效的抵抗统计攻击。34直方图一般像素的直方图可以反映出像素值的分布情况9，30，攻击者通过对像素值进行分析即可破解图像。本小节得到的测试图像加、解密的灰度直方图如图 9。中直方图的横坐标用来表示其灰度值，纵坐标用来表示像素的个数。a)原图直方图b)密图直方图c)解密图直方图图 9原图、密图、解密图的直方图分布从图 9 可以看出，原图像素灰度值出现的概率非常不均匀，最多的出现次数高达 9 000 多次，而最少的几乎为 0，这就说明图像的安全性非常弱，是很容易被攻击者分析出图像特点的。而经过本算法加密后得到的加密图像直方图变得非常均匀，从而可以有效的防止出现原始图像信息的泄露，进而有效的抵抗统计攻击。35鲁棒性分析本算法对医学图像进行加密，除了需要保证加密效果良好，还需要有一定的鲁棒性。鲁棒性包括抵抗噪声攻击的能力，裁剪攻击的能力等25，31。本实验以先对加密图进行一定程度的椒盐噪声干扰，然后再进行解密的方式，分析本算法抵抗噪声攻击的能力。对 Chest 加密图像分别进行 0000 5、0005 以及 005 的椒盐噪声干扰，然后进行解密过程得到解密图像，具体如图 10。a)0000 5 噪声b)0005 噪声c)005 噪声图 10不同程度椒盐噪声后生成的解密图从图 10 可以看出，经过不同程度的椒盐噪声污染后生成的解密图虽然清晰度受到了一定程度的干扰，但是仍然看出图像的有效信息，说明本算法抵抗噪声攻击能力较强，具有一定的鲁棒性。4结语本文提出了一种基于分段线性混沌映射的医学图像加密算法。该算法采用四轮加密的方式，确保加密算法的有效性和安全性。每一轮加密中，首先原始图像进行随机数插入，其次使用分段线性混沌映射进行像素值的替换，然后对生成矩阵以及替换前的矩阵进行按位异或操作，改变像素值的大小，最后将矩阵逆时针旋转进行下一轮加密，持续进行四轮加密，得到加密图像。从仿真实验和安全性分析结果可以看出，该算法不仅可以对医学图像进行加密而且可以对普通图像进行26大连民族大学学报第 25 卷加密，与目前算法相比安全性能更高，并且具有一定的鲁棒性，适合用于医学图像加密。参考文献:1高国静，吕庆文基于复合混沌的鲁棒型医学图像加密算法 J 计算机与现代化，2021，2021(6):107112 2孟维鸣 基于超混沌的医学图像加密算法研究D 洛阳:河南科技大学，2020 3王成宇 基于新型压缩感知技术的医学图像安全传输 D 北京:北京邮电大学，2021 4陈楠楠 基于混沌映射与 S 盒的图像加密算法研究 D 哈尔滨:东北林业大学，2020 5符翔龙 基于 DNA 编码和压缩感知的混沌图像加密算法的设计与实现 D 郑州:河南大学，2018 6纪芳混沌理论的舰船图像安全传输加密算法 J 舰船科学技术，2021，43(18):187189 7乔建平基于二维 Logistic 混沌系统的图像加密算法分析 J 江苏科技信息，2021，38(28):3436 8张雷，陈川，谭淇匀，等结合 S 盒与混沌映射的图像加密算法 J 北京邮电大学学报，2021，44(6):4047 9杜鑫昌，高瑜翔，曹远杰，等基于混沌压缩感知和DNA 编码的多图像加密算法J 无线电工程，2022，52(3):476483 10彭静静 基于混沌系统的图像加密算法研究 D 郑州:河南大学，2020 11谭琳基于 DNA 序列和混沌的图像加密算法J 信息系统工程，2014，2014(11):9092 12田军锋，彭静静，左宪禹，等基于循环移位和多混沌映射的图像加密算法J 计算机科学，2020，47(10):327331 13王海珍，廉佐政，谷文成基于 LP 混沌映射和 AES的图像加密算法J 计算机仿真，2021，38(10):217221 14王春娥，崔岩，赵少卿，等一种新四维超混沌系统的分岔分析及应用 J 河南科技大学学报(自然科学版)，2022，43(1):9198 15王倩改进 Arnold 算法和超混沌系统的医学图像加密研究 J 计算机时代，2018，2018(2):4347 16王倩基于位分解和超混沌映射的医学图像加密研究 J 计算机仿真，2019，36(1):209212 17禹峰 基于混沌的医疗图像加密算法研究D 北京:北京邮电大学，2019 18单泽润，龚涛基于 DCT 变换的 JPEG 图像加密算法 J 现代计算机，2021，27(25):6973 19黎学军，苗苗，李庆宇，等一种基于神经网络的图像加密算法J 大连工业大学学报，2021，40(5):358366 20刘杨曦，张烨基于压缩感知和位平面多重置乱的图像加密算法J 南昌大学学报(工科版)，2021，43(3):291300 21吕翔，陈利强，陆岳峰，等采用小波变换的抗剪切图像加密算法 J 浙江师范大学学报(自然科学版)，2021，44(4):381388 22唐旭，舒立鹏，范继，等基于 ossler 超混沌系统的军事图像加密算法研究J 火炮发射与控制学报，2022，43(2):3238 23王晓辉基于矩阵变换的图像加密方法J 微型电脑应用，2021，37(9):13 24周红亮，刘洪娟结合 DNA 编码的快速混沌图像加密算法J 东北大学学报(自然科学版)，2021，42(10):13911399 25乔建平，关来德一种基于混沌系统的分段式图像加密算法 J 南方农机，2021，52(11):161163 26韦丞婧，李国东基于超混沌系统和 Logistic 映射的视频图像加密设计J 计算机工程，2022，48(5):263271 27许向亮，李国东，戴婉莹新混沌系统与神经网络结合的图像加密算法 J 西华大学学报(自然科学版)，2021，40(5):4252 28张艳鹏，侯冬梅，杨倩，等基于混沌同步技术的图像加密算法设计研究J 现代电子技术，2021，44(19):3942 29张小康，曾国辉基于压缩感知的图像混沌加密算法 J 内蒙古科技与经济，2021，2021(12):8790 30蒋东华，朱礼亚，沈子懿，等结合二维压缩感知和混沌映射的双图像视觉安全加密算法 J 西安交通大学学报，2022，56(2):139148 31王纪一种基于分块置乱和混合混沌的图像加密算法 J 长江信息通信，2021，34(10):6466(责任编辑王楠楠)36第 1 期秦秋霞，等:一种基于分段线性混沌映射的医学图像加密算法