CN106875362B - 一种遥感图像的安全外包去噪方法 - Google Patents
一种遥感图像的安全外包去噪方法 Download PDFInfo
- Publication number
- CN106875362B CN106875362B CN201710089139.1A CN201710089139A CN106875362B CN 106875362 B CN106875362 B CN 106875362B CN 201710089139 A CN201710089139 A CN 201710089139A CN 106875362 B CN106875362 B CN 106875362B
- Authority
- CN
- China
- Prior art keywords
- array
- remote sensing
- sensing image
- paillier
- image
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000012946 outsourcing Methods 0.000 title claims abstract description 22
- 238000012545 processing Methods 0.000 claims abstract description 21
- 230000002688 persistence Effects 0.000 claims description 34
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000011084 recovery Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 4
- 238000003491 array Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000007792 addition Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000002131 composite material Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种遥感图像的安全外包去噪方法，所述去噪方法包括以下流程：遥感图像拆分处理；遥感图像随机排列；遥感图像Johnson–Lindenstrauss转换；遥感图像Paillier加密；遥感图像非局部去噪；遥感图像Paillier解密；遥感图像恢复排列；遥感图像合并处理。其优点在于，将在本地拆分并加密后的遥感图像放在云端进行图像去噪处理；不出现信息泄漏问题；只需在本地进行解密和合并后就可以得到该遥感图像完成去噪后的图像；可以提高原图像的PSNR的值。
Description
技术领域
本发明涉及图像去噪技术领域，具体地说，是一种遥感图像的安全外包去噪方法。
背景技术
遥感图像在获取和传输过程中，受各种噪声影响，使图像的边缘纹理等细节丢失，线条迷糊，从而导致图像的质量降低，增大了遥感图像的分析与处理难度。为获取清晰的、高质量的遥感图像必须对其进行去噪或降噪的预处理。
随着空天地立体化遥感技术的发展，遥感图像数据呈爆炸式增长，传统的基于单机模式的遥感图像处理暴露出许多问题，如处理速度慢、并发性能差等。因而传统遥感图像处理模式已经逐渐无法满足用户海量高精度遥感图像处理的需求，所以迫切需要寻求一种新的高效的遥感图像处理模式与处理方法。
云计算就是一种新的计算模式，它是一个典型的分布式、并行计算模型，该模型能够大大缩短计算任务的执行时间。随着遥感图像处理速度要求不断地提高，遥感图像处理越来越多地适合在分布式并行计算的云计算平台上进行。
由于云平台本身的开放性特点，为遥感图像的安全带来新威胁，无法保证海洋遥感影像的安全性。因此，为了保护数据的安全性，很有必要对于敏感数据进行加密。
中国发明专利CN201510220372.X，公开日为2016.12.07，公开了一种结合非局部先验的复合正则化图像去噪方法。但是该方法无法对图像进行加密及去噪。
因此，亟需一种能够对遥感图像加密并去噪的方法，既完成去噪又保证遥感图像安全，而目前关于这种去噪方法还未见报道。
发明内容
本发明的目的是针对现有技术中的不足，提供一种。
为实现上述目的，本发明采取的技术方案是：
一种遥感图像的安全外包去噪方法，所述去噪方法包括以下流程：
S1：遥感图像拆分处理；
S2：遥感图像随机排列；
S3：遥感图像Johnson–Lindenstrauss转换；
S4：遥感图像Paillier加密；
S5：遥感图像去噪；
S6：遥感图像Paillier解密；
S7：遥感图像恢复排列；
S8：遥感图像合并处理。
其中，所述的步骤S5中的遥感图像去噪为遥感图像非局部去噪，所述的步骤S1-S4、S6-S8在本地客户端进行，所述的步骤S5在云端进行。
所述的步骤S1遥感图像拆分处理中，其步骤为：读取遥感图像，设遥感图像的长为im_h，宽为im_w，波段数为im_band_number，按照波段将遥感图像拆分为im_band_number个图像，分别读取各个波段的16位深的灰度值。
所述的步骤S2遥感图像随机排列中，选取一个波段的遥感图像，对其进行随机排列，其步骤如下：
S21：生成一个长度为im_h×im_w的随机数序列，并记录其索引顺序为索引顺序A1；
S22：将这个随机数序列按照从小到大的顺序进行排列，并记录其索引顺序为索引顺序A2；
S23：根据索引顺序A1和索引顺序A2，得到图像置换数列和图像恢复数列；
S24：对图像进行置换，其具体过程如下，
S241：利用S23的图像进行随机排列的数列对遥感图像进行块置换；
S242：图像的16位深灰度值数组为im_array，获取以im_array[m,n]为中心的d×d大小的矩阵块，并将矩阵块转化为一行；
S243：将数组im_array中的所有点重复步骤S242，最终获得一个新数组；
S244：利用图像置换数列对S243中的新数组进行行置换，再次得到新数组；
S245：对S244中的新数组，即完成行置换的数组，进行变维操作，得到一个三维数组im_block_permutation_array，该数组的长为im_h，宽为im_w，高为d2；
S25：利用图像置换数列对遥感图像进行像素置换，得到一个二维数组im_pixel_permutation_array，该数组的长为im_h，宽为im_w。
所述的步骤S3遥感图像Johnson–Lindenstrauss转换中，对完成块置换的结果进行Johnson–Lindenstrauss转换，即将每一个1×d2的向量压缩成一个1×k的向量，得到一个三维数组im_JL_array，该数组的长为im_h，宽为im_w，高为k，其步骤如下，
S31：随机生成一个服从平均值为0，方差为1/k的高斯分布的d2×k的矩阵P；
S32：使用三维数组im_block_permutation_array，利用数组中的[m,n,0]、[m,n,1]、……、[m,n,d2-2]、[m,n,d2-1]，组成一个行矩阵α；
S33：随机生成一个服从平均值为0，方差为ζ的高斯分布的1×k的矩阵Δ；
S34：计算αP得到矩阵Y；
S35：计算Y+Δ得到α′；
S36：得到三维数组im_JL_array，即im_JL_array[m,n,0]＝α′[0,0]，im_JL_array[m,n,1]＝α′[0,1]，……，im_JL_array[m,n,k-2]＝α′[0,k-2]，im_JL array[m,n,k-1]＝α′[0,k-1]。
所述的步骤S4遥感图像Paillier加密中，对二维数组im_pixel_permutation_array进行paillier加密，得到一个二维数组im_paillier_array，该数组的长为im_h，宽为im_w，其具体步骤如下，
S41：选取两个随机素数p、q，并根据计算公式N＝pq，得出N的值；
S42：计算出paillier加密参数λ和g，得到paillier加密的公钥是(N,g)，私钥是λ；
S43：计算gm的值，m∈[0,65535]，并把所有结果保存在第一个表；
S44：随机生成65536个随机数r，计算(gN)r的值，并把所有结果保存在第二个表；
S45：对于二维数组的点im_pixel_permutation_array[i,j]，从第二个表中查找m的值等于该点对应的gm的值；
S46：从第二个表中随机抽取5个不同的(gN)r的值，即(gN)r1，(gN)r2，(gN)r3，(gN)r4，(gN)r5；
S47：im_paillier_array[i,j]＝gm(gN)r1(gN)r2(gN)r3(gN)r4(gN)r5mod N2；
S48：将二维数组im_pixel_permutation_array所有的点重复步骤S45-S47，得到数组im_paillier_array。
所述的步骤S5遥感图像去噪中，对数组im_JL array和im_paillier_array进行非局部去噪，得到一个二维数组im_paillier_result_array，该数组的长为im_h，宽为im_w，其具体步骤如下：
S51：选取h1的值，h1为进行非局部去噪的滤波参数；
S52：计算A的值，A为放大倍数，计算公式为A＝27×im_h×im_w；
S53：用im_JL array[m,n,0]，im_JL_array[m,n,1]，……，im_JL_array[m,n,k-2]，im_JL array[m,n,k-1]组成一个行向量，记为V(m,n)；
S54：用im_JL_array[i,j,0]，im_JL_array[i,j,1]，……，im_JL_array[i,j,k-2]，im_JL_array[m,n,k-1]组成一个行向量，记为V(i,j)；
S55：计算os_dis_square[i,j]的值，os_di_square[i,j]为像素点i和像素点j之间的欧式距离的平方，计算公式为
S56：计算dis[i,j]的值，计算公式为
S57：计算Z的值，计算公式为
其中Ω为整个图像；
S58：计算Q[i,j]的值，计算公式为
S59：计算W[i,j]的值，计算公式为
S510：确定W_number的值；
S511：确定os_dis_square数组中第W_number的值为dis_square，并计算h2的值，h2为进行非局部去噪的滤波参数，计算公式为
S512：将os_dis_square[i,j]与square进行比较，若os_dis_square[i,j]>square，则W[i,j]＝0，若os_dis_square[i,j]<square，则
S513：计算im_paillier_result_array[i]的值，计算公式为
其中Ω为整个图像；
S514：将数组im_paillier_result_array中的所有点重复步骤S53-S513，将数组im_paillier_array中的所有点重复步骤S54-S59、S512；
S515：得到数组im_paillier_result_array。
所述的步骤S6遥感图像Paillier解密中，利用私钥对二维数组im_paillier_result_array进行Paillier解密，得到一个二维数组im_permutation_result_array，该数组的长为im_h，宽为im_w，其具体步骤如下：
S61：计算U的值，计算公式为
S62：对于二维数组中的点im_paillier_result_array[i,j]，计算其对应像素点c的NL的值，计算公式为
S63：
S64：将数组im_paillier_result_array的所有点重复S62-S63，得到数组im_permutation_result_array。
所述的步骤S7遥感图像恢复排列中，其步骤为：利用图像恢复数列，对数组im_permutation_result_array进行遥感图像的像素置换恢复，得到一个二维数组im_result_array，该数组的长为im_h，宽为im_w。
所述的步骤S8遥感图像合并处理中，其步骤为：将每一个波段的遥感图像重复S2-S7，进行合并处理，得到完成去噪后的完整遥感图像。
本发明优点在于：
1、将在本地拆分并加密后的遥感图像放在云端进行图像去噪处理；
2、不出现信息泄漏问题；
3、只需在本地进行解密和合并后就可以得到该遥感图像完成去噪后的图像；
4、可以提高原图像的PSNR的值。
附图说明
附图1是本发明的一种遥感图像的安全外包去噪方法的总工作流程图。
附图2是本发明的一种遥感图像的安全外包去噪方法的具体工作流程图。
附图3是本发明的一种遥感图像的安全外包去噪方法的单波段遥感图像去噪示意图。
附图4是本发明的一种遥感图像的安全外包去噪方法的遥感图像Johnson–Lindenstrauss转换流程图。
具体实施方式
下面结合附图对本发明提供的具体实施方式作详细说明。
实施例1
参照图1，本发明的一种遥感图像的安全外包去噪方法的总工作流程如下：
S1：遥感图像拆分处理；
S2：遥感图像随机排列；
S3：遥感图像Johnson–Lindenstrauss转换；
S4：遥感图像Paillier加密；
S5：遥感图像去噪；
S6：遥感图像Paillier解密；
S7：遥感图像恢复排列；
S8：遥感图像合并处理。
其中，所述的步骤S5中的遥感图像去噪为遥感图像非局部去噪，所述的步骤S1-S4、S6-S8在本地客户端进行，所述的步骤S5在云端进行。
实施例2
参照图2以及图3，本发明的一种遥感图像的安全外包去噪方法的具体工作流程如下：
S1：在本地客户端进行遥感图像拆分处理
遥感图像拆分处理的步骤为：读取遥感图像，设遥感图像的长为im_h，宽为im_w，波段数为im_band_number，按照波段将遥感图像拆分为im_band_number个图像，分别读取各个波段的16位深的灰度值。
S2：在本地客户端进行遥感图像随机排列
选取一个波段的遥感图像，对其进行随机排列，其步骤如下：
S21：生成一个长度为im_h×im_w的随机数序列，并记录其索引顺序为索引顺序A1；
S22：将这个随机数序列按照从小到大的顺序进行排列，并记录其索引顺序为索引顺序A2；
S23：根据索引顺序A1和索引顺序A2，得到图像置换数列和图像恢复数列；
S24：对图像进行置换，其具体过程如下，
S241：利用S23的图像进行随机排列的数列对遥感图像进行块置换；
S242：图像的16位深灰度值数组为im_array，获取以im_array[m,n]为中心的d×d大小的矩阵块，并将矩阵块转化为一行；
S243：将数组im_array中的所有点重复步骤S242，最终获得一个新数组；
S244：利用图像置换数列对S243中的新数组进行行置换，再次得到新数组；
S245：对S244中的新数组，即完成行置换的数组，进行变维操作，得到一个三维数组im_block_permutation_array，该数组的长为im_h，宽为im_w，高为d2；
S25：利用图像置换数列对遥感图像进行像素置换，得到一个二维数组im_pixel_permutation_array，该数组的长为im_h，宽为im_w。
S3：在本地客户端进行遥感图像Johnson–Lindenstrauss转换
参照图4，对完成块置换的结果进行Johnson–Lindenstrauss转换，即将每一个1×d2的向量压缩成一个1×k的向量，得到一个三维数组im_JL_array，该数组的长为im_h，宽为im_w，高为k，其步骤如下，
S31：随机生成一个服从平均值为0，方差为1/k的高斯分布的d2×k的矩阵P；
S32：使用三维数组im_block_permutation_array，利用数组中的[m,n,0]、[m,n,1]、……、[m,n,d2-2]、[m,n,d2-1]，组成一个行矩阵α；
S33：随机生成一个服从平均值为0，方差为ζ的高斯分布的1×k的矩阵Δ；
S34：计算αP得到矩阵Y；
S35：计算Y+Δ得到α′；
S36：得到三维数组im_JL_array，即im_JL array[m,n,0]＝α′[0,0]，im_JL array[m,n,1]＝α′[0,1]，……，im_JL array[m,n,k-2]＝α′[0,k-2]，im_JL array[m,n,k-1]＝α′[0,k-1]。
S4：在本地客户端进行遥感图像Paillier加密
对二维数组im_pixel_permutation_array进行paillier加密，得到一个二维数组im_paillier_array，该数组的长为im_h，宽为im_w，其具体步骤如下，
S41：选取两个随机素数p、q，并根据计算公式N＝pq，得出N的值；
S42：计算出paillier加密参数λ和g，得到paillier加密的公钥是(N,g)，私钥是λ；
S43：计算gm的值，m∈[0,65535]，并把所有结果保存在第一个表；
S44：随机生成65536个随机数r，计算(gN)r的值，并把所有结果保存在第二个表；
S45：对于二维数组的点im_pixel_permutation_array[i,j]，从第二个表中查找m的值等于该点对应的gm的值；
S46：从第二个表中随机抽取5个不同的(gN)r的值，即(gN)r1，(gN)r2，(gN)r3，(gN)r4，(gN)r5；
S47：im_paillier_array[i,j]＝gm(gN)r1(gN)r2(gN)r3(gN)r4(gN)r5mod N2；
S48：将二维数组im_pixel_permutation_array所有的点重复步骤S45-S47，得到数组im_paillier_array。
S5：在云端进行遥感图像非局部去噪
对数组im_JL array和im_paillier_array进行非局部去噪，得到一个二维数组im_paillier_result_array，该数组的长为im_h，宽为im_w，其具体步骤如下，
S51：选取h1的值，h1为进行非局部去噪的滤波参数；
S52：计算A的值，A为放大倍数，计算公式为A＝27×im_h×im_w；
S53：用im_JL array[m,n,0]，im_JL_array[m,n,1]，……，im_JL_array[m,n,k-2]，im_JL array[m,n,k-1]组成一个行向量，记为V(m,n)；
S54：用im_JL_array[i,j,0]，im_JL_array[i,j,1]，……，im_JL_array[i,j,k-2]，im_JL_array[m,n,k-1]组成一个行向量，记为V(i,j)；
S55：计算os_dis_square[i,j]的值，os_di_square[i,j]为像素点i和像素点j之间的欧式距离的平方，计算公式为
S56：计算dis[i,j]的值，计算公式为
S57：计算Z的值，计算公式为
其中Ω为整个图像；
S58：计算Q[i,j]的值，计算公式为
S59：计算W[i,j]的值，计算公式为
S510：确定W_number的值；
S511：确定os_dis_square数组中第W_number的值为dis_square，并计算h2的值，h2为进行非局部去噪的滤波参数，计算公式为
S512：将os_dis_square[i,j]与square进行比较，若os_dis_square[i,j]>square，则W[i,j]＝0，若os_dis_square[i,j]<square，则
S513：计算im_paillier_result_array[i]的值，计算公式为
其中Ω为整个图像；
S514：将数组im_paillier_result_array中的所有点重复步骤S53-S513，将数组im_paillier_array中的所有点重复步骤S54-S59、S512；
S515：得到数组im_paillier_result_array。
S6：在本地客户端进行遥感图像Paillier解密
利用私钥对二维数组im_paillier_result_array进行Paillier解密，得到一个二维数组im_permutation_result_array，该数组的长为im_h，宽为im_w，其具体步骤如下，
S61：计算U的值，计算公式为
S62：对于二维数组中的点im_paillier_result_array[i,j]，计算其对应为c的NL的值，计算公式为
S63：
S64：将数组im_paillier_result_array的所有点重复S62-S63，得到数组im_permutation_result_array。
S7：在本地客户端进行遥感图像恢复排列
遥感图像恢复排列的步骤为：利用图像恢复数列，对数组im_permutation_result_array进行遥感图像的像素置换恢复，得到一个二维数组im_result_array，该数组的长为im_h，宽为im_w。
S8：在本地客户端进行遥感图像合并处理
遥感图像合并处理的步骤为：将每一个波段的遥感图像重复S2-S7，进行合并处理，得到完成去噪后的完整遥感图像。
本发明的一种遥感图像的安全外包去噪方法的优点在于：将在本地拆分并加密后的遥感图像放在云端进行图像去噪处理；不出现信息泄漏问题；只需在本地进行解密和合并后就可以得到该遥感图像完成去噪后的图像；可以提高原图像的PSNR的值。
实施例3
本发明的一种遥感图像的安全外包去噪方法的具体实施过程如下：
S1：在本地客户端进行遥感图像拆分处理
读取Landsat8卫星的遥感图像，该遥感图像的长为512，宽为512，波段数为11，按照波段将遥感图像拆分为11个图像，并读取各个波段的16位深的灰度值。
S2：在本地客户端进行遥感图像随机排列
选取一个波段的遥感图像，对其进行随机排列，其步骤如下：
S21：生成一个长度为512×512的随机数序列，并记录其索引顺序为索引顺序A1；
S22：将这个随机数序列按照从小到大的顺序进行排列，并记录其索引顺序为索引顺序A2；
S23：根据索引顺序A1和索引顺序A2，得到图像置换数列和图像恢复数列；
S24：对图像进行置换，其具体过程如下，
S241：利用S23的图像进行随机排列的数列对遥感图像进行块置换；
S242：图像的16位深灰度值数组为im_array，获取以im_array[m,n]为中心的5×5大小的矩阵块，并将矩阵块转化为一行；
S243：将数组im_array中的所有点重复步骤S242，最终获得一个新数组；
S244：利用图像置换数列对S243中的新数组进行行置换，再次得到新数组；
S245：对S244中的新数组，即完成行置换的数组，进行变维操作，得到一个三维数组im_block_permutation_array，该数组的长为512，宽为512，高为25；
S25：利用图像置换数列对遥感图像进行像素置换，得到一个二维数组im_pixel_permutation_array，该数组的长为512，宽为512。
S3：在本地客户端进行遥感图像Johnson–Lindenstrauss转换
对完成块置换的结果进行Johnson–Lindenstrauss转换，即将每一个1×25的向量压缩成一个1×k的向量，得到一个三维数组im_JL_array，该数组的长为512，宽为512，高为k，其步骤如下，
S31：随机生成一个服从平均值为0，方差为1/k的高斯分布的d2×k的矩阵P；
S32：使用三维数组im_block_permutation_array，利用数组中的[m,n,0]、[m,n,1]、……、[m,n,23]、[m,n,24]，组成一个行矩阵α；
S33：随机生成一个服从平均值为0，方差为ζ的高斯分布的1×k的矩阵Δ；
S34：计算α×P得到矩阵Y；
S35：计算Y+Δ得到α′；
S36：得到三维数组im_JL_array，即im_JL_array[m,n,0]＝α′[0,0]，im_JL_array[m,n,1]＝α′[0,1]，……，im_JL_array[m,n,k-2]＝α′[0,k-2]，im_JL_array[m,n,k-1]＝α′[0,k-1]。
S4：在本地客户端进行遥感图像Paillier加密
对二维数组im_pixel_permutation_array进行paillier加密，得到一个二维数组im_paillier_array，该数组的长为512，宽为512，其具体步骤如下，
S41：选取两个随机素数p、q，并根据计算公式N＝pq，得出N的值；
S42：计算出paillier加密参数λ和g，得到paillier加密的公钥是(N,g)，私钥是λ；
S43：计算gm的值，m∈[0,65535]，并把所有结果保存在第一个表；
S44：随机生成65536个随机数r，计算(gN)r的值，并把所有结果保存在第二个表；
S45：对于二维数组的点im_pixel_permutation_array[i,j]，从第二个表中查找m的值等于该点对应的gm的值；
S46：从第二个表中随机抽取5个不同的(gN)r的值，即(gN)r1，(gN)r2，(gN)r3，(gN)r4，(gN)r5；
S47：im_paillier_array[i,j]＝gm(gN)r1(gN)r2(gN)r3(gN)r4(gN)r5mod N2；
S48：将二维数组im_pixel_permutation_array所有的点重复步骤S45-S47，得到数组im_paillier_array。
S5：在云端进行遥感图像非局部去噪
对数组im_JL array和im_paillier_array进行非局部去噪，得到一个二维数组im_paillier_result_array，该数组的长为512，宽为512，其具体步骤如下，
S51：选取h1的值为384，h1为进行非局部去噪的滤波参数；
S52：计算A的值，A为放大倍数，计算公式为A＝27×512×512；
S53：用im_JL array[m,n,0]，im_JL_array[m,n,1]，……，im_JL_array[m,n,k-2]，im_JL array[m,n,k-1]组成一个行向量，记为V(m,n)；
S54：用im_JL_array[i,j,0]，im_JL_array[i,j,1]，……，im_JL_array[i,j,k-2]，im_JL_array[m,n,k-1]组成一个行向量，记为V(i,j)；
S55：计算os_dis_square[i,j]的值，os_di_square[i,j]为像素点i和像素点j之间的欧式距离的平方，计算公式为
S56：计算dis[i,j]的值，计算公式为
S57：计算Z的值，计算公式为
其中Ω为整个图像；
S58：计算Q[i,j]的值，计算公式为
S59：计算W[i,j]的值，计算公式为
S510：确定W_number为128；
S511：确定os_dis_square数组中第W_number的值为dis_square，并计算h2的值，h2为进行非局部去噪的滤波参数，计算公式为
S512：将os_dis_square[i,j]与square进行比较，若os_dis_square[i,j]>square，则W[i,j]＝0，若os_dis_square[i,j]<square，则
S513：计算im_paillier_result_array[i]的值，计算公式为
其中Ω为整个图像；
S514：将数组im_paillier_result_array中的所有点重复步骤S53-S513，将数组im_paillier_array中的所有点重复步骤S54-S59、S512；
S515：得到数组im_paillier_result_array。
S6：在本地客户端进行遥感图像Paillier解密
利用私钥对二维数组im_paillier_result_array进行Paillier解密，得到一个二维数组im_permutation_result_array，该数组的长为512，宽为512，其具体步骤如下，
S61：计算U的值，计算公式为
S62：对于二维数组中的点im_paillier_result_array[i,j]，计算其对应像素点c的NL的值，计算公式为
S63：
S64：将数组im_paillier_result_array的所有点重复S62-S63，得到数组im_permutation_result_array。
S7：在本地客户端进行遥感图像恢复排列
遥感图像恢复排列的步骤为：利用图像恢复数列，对数组im_permutation_result_array进行遥感图像的像素置换恢复，得到一个二维数组im_result_array，该数组的长为512，宽为512。
S8：在本地客户端进行遥感图像合并处理
遥感图像合并处理的步骤为：将每一个波段的遥感图像重复S2-S7，进行合并处理，得到完成去噪后的完整遥感图像。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员，在不脱离本发明方法的前提下，还可以做出若干改进和补充，这些改进和补充也应视为本发明的保护范围。
Claims (9)
1.一种遥感图像的安全外包去噪方法，其特征在于，所述的安全外包去噪方法包括以下流程：
S1：遥感图像拆分处理；
S2：遥感图像随机排列；
S3：遥感图像Johnson–Lindenstrauss转换；
S4：遥感图像Paillier加密；
S5：遥感图像去噪；
S6：遥感图像Paillier解密；
S7：遥感图像恢复排列；
S8：遥感图像合并处理；
其中，所述的步骤S5中的遥感图像去噪为遥感图像非局部去噪，所述的步骤S1-S4、S6-S8在本地客户端进行，所述的步骤S5在云端进行。
2.根据权利要求1所述的安全外包去噪方法，其特征在于，所述的步骤S1遥感图像拆分处理中，其步骤为：读取遥感图像，设遥感图像的长为im_h，宽为im_w，波段数为im_band_number，按照波段将遥感图像拆分为im_band_number个图像，分别读取各个波段的16位深的灰度值。
3.根据权利要求1所述的安全外包去噪方法，其特征在于，所述的步骤S2遥感图像随机排列中，选取一个波段的遥感图像，对其进行随机排列，其步骤如下：
S21：生成一个长度为im_h×im_w的随机数序列，并记录其索引顺序为索引顺序A1；
S22：将这个随机数序列按照从小到大的顺序进行排列，并记录其索引顺序为索引顺序A2；
S23：根据索引顺序A1和索引顺序A2，得到图像置换数列和图像恢复数列；
S24：对图像进行置换，其具体过程如下，
S241：利用S23的图像进行随机排列的数列对遥感图像进行块置换；
S242：图像的16位深灰度值数组为im_array，获取以im_array[m,n]为中心的d×d大小的矩阵块，并将矩阵块转化为一行；
S243：将数组im_array中的所有点重复步骤S242，最终获得一个新数组；
S244：利用图像置换数列对S243中的新数组进行行置换，再次得到新数组；
S245：对S244中的新数组，即完成行置换的数组，进行变维操作，得到一个三维数组im_block_permutation_array，该数组的长为im_h，宽为im_w，高为d2；
S25：利用图像置换数列对遥感图像进行像素置换，得到一个二维数组im_pixel_permutation_array，该数组的长为im_h，宽为im_w。
4.根据权利要求1所述的安全外包去噪方法，其特征在于，所述的步骤S3遥感图像Johnson–Lindenstrauss转换中，对完成块置换的结果进行Johnson–Lindenstrauss转换，即将每一个1×d2的向量压缩成一个1×k的向量，得到一个三维数组im_JL_array，该数组的长为im_h，宽为im_w，高为k，其步骤如下：
S31：随机生成一个服从平均值为0，方差为1/k的高斯分布的d2×k的矩阵P；
S32：使用三维数组im_block_permutation_array，利用数组中的[m,n,0]、[m,n,1]、……、[m,n,d2-2]、[m,n,d2-1]，组成一个行矩阵α；
S33：随机生成一个服从平均值为0，方差为ζ的高斯分布的1×k的矩阵Δ；
S34：计算αP得到矩阵Y；
S35：计算Y+Δ得到α′；
S36：得到三维数组im_JL_array，即im_JL_array[m,n,0]＝α′[0,0]，
im_JL_array[m,n,1]＝α′[0,1]，……，im_JL_array[m,n,k-2]＝α′[0,k-2]，
im_JL_array[m,n,k-1]＝α′[0,k-1]。
5.根据权利要求1所述的安全外包去噪方法，其特征在于，所述的步骤S4遥感图像Paillier加密中，对二维数组im_pixel_permutation_array进行paillier加密，得到一个二维数组im_paillier_array，该数组的长为im_h，宽为im_w，其具体步骤如下：
S41：选取两个随机素数p、q，并根据计算公式N＝pq，得出N的值；
S42：计算出paillier加密参数λ和g，得到paillier加密的公钥是(N,g)，私钥是λ；
S43：计算gm的值，m∈[0,65535]，并把所有结果保存在表1；
S44：随机生成65536个随机数r，计算(gN)r的值，并把所有结果保存在表2；
S45：对于二维数组的点im_pixel_permutation_array[i,j]，从表2中查找m的值等于该点对应的gm的值；
S46：从表2中随机抽取5个不同的(gN)r的值，即(gN)r1，(gN)r2，(gN)r3，(gN)r4，(gN)r5；
S47：im_paillier_array[i,j]＝gm(gN)r1(gN)r2(gN)r3(gN)r4(gN)r5mod N2；
S48：将二维数组im_pixel_permutation_array所有的点重复步骤S45-S47，
得到数组im_paillier_array。
6.根据权利要求1所述的安全外包去噪方法，其特征在于，所述的步骤S5遥感图像去噪中，对数组im_JL array和im_paillier_array进行非局部去噪，得到一个二维数组im_paillier_result_array，该数组的长为im_h，宽为im_w，其具体步骤如下：
S51：选取h1的值，h1为进行非局部去噪的滤波参数；
S52：计算A的值，A为放大倍数，计算公式为A＝27×im_h×im_w；
S53：用im_JL array[m,n,0]，im_JL_array[m,n,1]，……，im_JL_array[m,n,k-2]，im_JL array[m,n,k-1]组成一个行向量，记为V(m,n)；
S54：用im_JL_array[i,j,0]，im_JL_array[i,j,1]，……，im_JL_array[i,j,k-2]，im_JL_array[m,n,k-1]组成一个行向量，记为V(i,j)；
S55：计算os_dis_square[i,j]的值，os_di_square[i,j]为像素点i和像素点j之间的欧式距离的平方，计算公式为
S56：计算dis[i,j]的值，计算公式为
S57：计算Z的值，计算公式为
其中Ω为整个图像；
S58：计算Q[i,j]的值，计算公式为
S59：计算W[i,j]的值，计算公式为
S510：确定W_number的值；
S511：确定os_dis_square数组中第W_number的值为dis_square，并计算h2的值，h2为进行非局部去噪的滤波参数，计算公式为
S512：将os_dis_square[i,j]与square进行比较，若os_dis_square[i,j]>square，则W[i,j]＝0，若os_dis_square[i,j]<square，则
S513：计算im_paillier_result_array[i]的值，计算公式为
其中Ω为整个图像；
S514：将数组im_paillier_result_array中的所有点重复步骤S53-S513，将数组im_paillier_array中的所有点重复步骤S54-S59、S512；
S515：得到数组im_paillier_result_array。
7.根据权利要求1所述的安全外包去噪方法，其特征在于，所述的步骤S6遥感图像Paillier解密中，利用私钥对二维数组im_paillier_result_array进行Paillier解密，得到一个二维数组im_permutation_result_array，该数组的长为im_h，宽为im_w，其具体步骤如下：
S61：计算U的值，计算公式为
S62：对于二维数组中的点im_paillier_result_array[m,n]，计算其对应像素点c的NL的值，计算公式为
S63：
S64：将数组im_paillier_result_array的所有点重复S62-S63，得到数组im_permutation_result_array。
8.根据权利要求1所述的安全外包去噪方法，其特征在于，所述的步骤S7遥感图像恢复排列中，其步骤为：利用图像恢复数列，对数组im_permutation_result_array进行遥感图像的像素置换恢复，得到一个二维数组im_result_array，该数组的长为im_h，宽为im_w。
9.根据权利要求1所述的安全外包去噪方法，其特征在于，所述的步骤S8遥感图像合并处理中，其步骤为：将每一个波段的遥感图像重复S2-S7，进行合并处理，得到完成去噪后的完整遥感图像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710089139.1A CN106875362B (zh) | 2017-02-20 | 2017-02-20 | 一种遥感图像的安全外包去噪方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710089139.1A CN106875362B (zh) | 2017-02-20 | 2017-02-20 | 一种遥感图像的安全外包去噪方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106875362A CN106875362A (zh) | 2017-06-20 |
| CN106875362B true CN106875362B (zh) | 2020-04-10 |
Family
ID=59167180
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710089139.1A Active CN106875362B (zh) | 2017-02-20 | 2017-02-20 | 一种遥感图像的安全外包去噪方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106875362B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107194895B (zh) * | 2017-05-27 | 2020-07-28 | 上海海洋大学 | 一种针对多帧遥感图像的安全外包融合去噪方法 |
| GB2567697B (en) * | 2017-10-23 | 2022-02-16 | Thales Holdings Uk Plc | Data encryption and processing |
| CN110032877A (zh) * | 2019-03-01 | 2019-07-19 | 阿里巴巴集团控股有限公司 | 图像存取方法及其系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102081513A (zh) * | 2011-01-24 | 2011-06-01 | 山东大学 | Aes加密算法中列混淆过程指令优化方法及其指令集处理器 |
| CN102903082A (zh) * | 2012-09-14 | 2013-01-30 | 重庆大学 | 一种基于开关中值滤波的加密图像去噪方法 |
| CN105323212A (zh) * | 2014-06-12 | 2016-02-10 | 江苏真云计算科技有限公司 | 一种安全外包的基于属性的加密技术 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8843762B2 (en) * | 2009-09-04 | 2014-09-23 | Gradiant, Centro Tecnolóxico de Telecomunicacións de Galicia | Cryptographic system for performing secure iterative computations and signal processing directly on encrypted data in untrusted environments |
-
2017
- 2017-02-20 CN CN201710089139.1A patent/CN106875362B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102081513A (zh) * | 2011-01-24 | 2011-06-01 | 山东大学 | Aes加密算法中列混淆过程指令优化方法及其指令集处理器 |
| CN102903082A (zh) * | 2012-09-14 | 2013-01-30 | 重庆大学 | 一种基于开关中值滤波的加密图像去噪方法 |
| CN105323212A (zh) * | 2014-06-12 | 2016-02-10 | 江苏真云计算科技有限公司 | 一种安全外包的基于属性的加密技术 |
Non-Patent Citations (1)
| Title |
|---|
| "PRIVACY PRESERVING APPROXIMATION OF L1 DISTANCE";Shantanu Rane.et al;《2010 IEEE International Conference on Multimedia and Expo》;20100723;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106875362A (zh) | 2017-06-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103778593B (zh) | 基于图像关联分解的量子图像加密方法 | |
| CN106875362B (zh) | 一种遥感图像的安全外包去噪方法 | |
| CN108648133B (zh) | 一种结合块旋转和马赛克的无嵌入伪装方法 | |
| CN101044761A (zh) | 安全处理输入图像的方法 | |
| CN107194895B (zh) | 一种针对多帧遥感图像的安全外包融合去噪方法 | |
| CN101044760A (zh) | 安全处理一系列输入图像的方法和系统 | |
| CN104851071A (zh) | 一种基于三维混沌系统的数字图像加密方法 | |
| Askari et al. | A novel visual secret sharing scheme without image size expansion | |
| CN111064558A (zh) | 一种基于云计算的同态加密矩阵连乘安全外包方法 | |
| Sun et al. | An image encryption algorithm utilizing Mandelbrot set | |
| Khizrai et al. | Image encryption using different techniques for high security transmission over a network | |
| CN112330551A (zh) | 一种基于秘密共享的遥感图像外包降噪方法 | |
| CN108717683B (zh) | 一种结合密钥和随机正交张量基的密图伪装恢复方法 | |
| Ergun | Privacy preserving face recognition in encrypted domain | |
| CN110569945A (zh) | 基于单强度记录相位恢复过程的qr码加密方法 | |
| CN116309164A (zh) | 一种图像处理方法及装置 | |
| Yang et al. | Secret image sharing without keeping permutation key | |
| CN106683030B (zh) | 基于量子多图像模型与三维变换的量子多图加密算法 | |
| Jiang et al. | Compressive-sensing-based double-image encryption algorithm combining double random phase encoding with Josephus traversing operation | |
| CN112184841A (zh) | 块替换生成式信息隐藏及恢复方法、设备及介质 | |
| Hosseini et al. | A parallel image encryption based on elementary cellular automata using two processors | |
| Dhane et al. | A new framework for secure biometric data transmission using block-wise reversible data hiding through encryption | |
| Ding et al. | Participants Increasing for Shamir’s Polynomial-Based Secret Image Sharing Scheme | |
| CN108154460B (zh) | 一种基于局部感知思想的混沌图像加密方法和装置 | |
| Hussien | Encryption of stereo images after estimated the motion using spatially dependent algorithms |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |