本發(fā)明涉及遙感圖像加密及去噪
技術(shù)領(lǐng)域：
，具體地說(shuō)，是一種針對(duì)多幀遙感圖像的安全外包融合去噪方法。
背景技術(shù)：
：遙感圖像在獲取和傳輸過(guò)程中，受各種噪聲影響，使圖像的邊緣紋理等細(xì)節(jié)丟失，線條迷糊，從而導(dǎo)致圖像的質(zhì)量降低，增大了遙感圖像的分析與處理難度。為獲取清晰的、高質(zhì)量的遙感圖像必須對(duì)其進(jìn)行去噪或降噪的預(yù)處理。隨著空天地立體化遙感技術(shù)的發(fā)展，遙感圖像數(shù)據(jù)呈爆炸式增長(zhǎng)，傳統(tǒng)的基于單機(jī)模式的遙感圖像處理暴露出許多問(wèn)題，如處理速度慢、并發(fā)性能差等。因而傳統(tǒng)遙感圖像處理模式已經(jīng)逐漸無(wú)法滿足用戶海量高精度遙感圖像處理的需求，所以迫切需要尋求一種新的高效的遙感圖像處理模式與處理方法。云計(jì)算就是一種新的計(jì)算模式，它是一個(gè)典型的分布式、并行計(jì)算模型，該模型能夠大大縮短計(jì)算任務(wù)的執(zhí)行時(shí)間。隨著遙感圖像處理速度要求不斷地提高，遙感圖像處理越來(lái)越多地適合在分布式并行計(jì)算的云計(jì)算平臺(tái)上進(jìn)行。由于云平臺(tái)本身的開(kāi)放性特點(diǎn)，為遙感圖像的安全帶來(lái)新威脅，無(wú)法保證海洋遙感影像的安全性。因此，為了保護(hù)數(shù)據(jù)的安全性，很有必要對(duì)于敏感數(shù)據(jù)進(jìn)行加密。中國(guó)發(fā)明專利cn201610159959.9，公開(kāi)日為2016.08.10，公開(kāi)了一種雙目視覺(jué)的圖像超分辨率融合去噪方法，其步驟為輸入待超分辨率融合的雙目圖像，全局位置配準(zhǔn)圖像，局部位置匹配圖像，更新雙目圖像，圖像超分辨率融合，圖像去噪，輸出最終高分辨率圖像。該方法在圖像配準(zhǔn)中加入局部位置配準(zhǔn)方法，并采用拉普拉斯算子應(yīng)用于圖像超分辨率融合，再對(duì)融合圖像使用非局部均值濾波方法去噪，克服了現(xiàn)有技術(shù)無(wú)法解決圖像局部位置不匹配，細(xì)節(jié)增強(qiáng)不足，以及不具備抑制噪聲的缺陷，得到圖像細(xì)節(jié)增加，噪聲減少的超分辨率融合圖像。但是這種方法無(wú)法保證在外包過(guò)程中信息安全性。因此，亟需一種能夠提高多幀遙感圖像去噪效率，避免外包過(guò)程中出現(xiàn)信息泄漏的安全外包融合去噪方法，而目前關(guān)于這種融合去噪方法還未見(jiàn)報(bào)道。技術(shù)實(shí)現(xiàn)要素：本發(fā)明的目的是針對(duì)現(xiàn)有技術(shù)中的不足，提供一種針對(duì)多幀遙感圖像的安全外包融合去噪方法。為實(shí)現(xiàn)上述目的，本發(fā)明采取的技術(shù)方案是：一種針對(duì)多幀遙感圖像的安全外包融合去噪方法，包括本地平臺(tái)和云平臺(tái)，所述的本地平臺(tái)包括遙感圖像拆分模塊、遙感圖像置換模塊、遙感圖像轉(zhuǎn)換模塊、遙感圖像加密模塊、遙感圖像解密模塊、遙感圖像恢復(fù)排列模塊、遙感圖像合并模塊，所述的云平臺(tái)包括遙感圖像融合去噪模塊，所述的遙感圖像拆分模塊為對(duì)遙感圖像按照?qǐng)D像波段數(shù)進(jìn)行拆分，所述的遙感圖像置換模塊為對(duì)遙感圖像進(jìn)行隨機(jī)置換，所述的遙感圖像轉(zhuǎn)換模塊為對(duì)遙感圖像進(jìn)行johnson-lindenstrauss轉(zhuǎn)換，所述的遙感圖像加密模塊為對(duì)遙感圖像進(jìn)行paillier加密，所述的遙感圖像解密模塊為對(duì)遙感圖像進(jìn)行paillier解密，其中，所述的針對(duì)多幀遙感圖像的安全外包融合去噪方法的工作流程為，s01：遙感圖像拆分；s02：遙感圖像置換；s03：遙感圖像轉(zhuǎn)換；s04：遙感圖像加密；s05：遙感圖像融合去噪；s06：遙感圖像解密；s07：遙感圖像恢復(fù)排列；s08：遙感圖像合并。作為一種優(yōu)選的技術(shù)方案，在所述的步驟s01中，在本地平臺(tái)，讀取針對(duì)同一場(chǎng)景的多幀遙感圖像，遙感圖像的個(gè)數(shù)為n，遙感圖像的長(zhǎng)為im_h，遙感圖像的寬為im_w，遙感圖像的波段數(shù)為im_band_number，按照波段數(shù)把每個(gè)遙感圖像拆分成im_band_number個(gè)圖像，分別讀取個(gè)波段圖像的16位深灰度值，并對(duì)每一個(gè)波段的圖像進(jìn)行步驟s02-s07。作為一種優(yōu)選的技術(shù)方案，在所述的步驟s02中，在本地平臺(tái)，為打亂圖像像素的基本位置，對(duì)n個(gè)圖像的該波段圖像進(jìn)行隨機(jī)置換，以第t個(gè)圖像為例，其工作流程如下：s021：生成一個(gè)長(zhǎng)度為im_h×im_w的隨機(jī)數(shù)序列，記錄其索引順序i1；s022：對(duì)隨機(jī)數(shù)序列按照從小到大進(jìn)行排列，生成一個(gè)新序列，記錄其索引順序i2；s023：根據(jù)i1→i2，得到圖像置換數(shù)列，根據(jù)i2→i1，得到圖像恢復(fù)數(shù)列；s024：對(duì)遙感圖像進(jìn)行置換，圖像的16位深灰度值數(shù)組為im_array_t，以數(shù)組中第a行第b列的點(diǎn)為例，對(duì)于數(shù)組中的每一點(diǎn)做如下操作：s0241：獲取以im_array_t[a,b]為中心的d×d大小的矩形塊，并將該矩形塊轉(zhuǎn)化成一行，所有的轉(zhuǎn)化后的行組成一個(gè)新數(shù)組；s0242：利用圖像置換序列對(duì)新數(shù)組進(jìn)行行置換；s0243：對(duì)完成行置換的數(shù)組進(jìn)行變維操作，得到一個(gè)三維數(shù)組im_block_permutation_array_t，該數(shù)組的長(zhǎng)為im_h，寬為im_w，高為d2；s025：利用圖像置換序列對(duì)遙感圖像進(jìn)行像素置換，得到一個(gè)二維數(shù)組im_pixel_permutation_array_t，該數(shù)組的長(zhǎng)為im_h，寬為im_w。作為一種優(yōu)選的技術(shù)方案，在所述的步驟s03中，在本地平臺(tái)，對(duì)n個(gè)圖像的該波段圖像進(jìn)行johnson-lindenstrauss轉(zhuǎn)換，以第t個(gè)圖像為例，即對(duì)im_block_permutation_array_t進(jìn)行johnson-lindenstrauss轉(zhuǎn)換，將每一個(gè)1×d2的向量壓縮呈一個(gè)1×k的向量，得到一個(gè)三維數(shù)組im_jl_array_t，該數(shù)組的長(zhǎng)為im_h，寬為im_w，高為k，其工作流程如下：s031：隨機(jī)生成一個(gè)服從平均值為0，方差為1/k的高斯分布的d2×k的矩陣p；以數(shù)組im_array_t中第a行第b列這個(gè)點(diǎn)為例，對(duì)于數(shù)組中的每一個(gè)點(diǎn)重復(fù)步驟s032-s036；s032：將數(shù)組im_block_permutation_array_t中的im_block_permutation_array_t[a,b,0]，im_block_permutation_array_t[a,b,1]，…，im_block_permutation_array_t[a,b,d2-2]，im_block_permutation_array_t[a,b,d2-1]，組成一個(gè)行矩陣α；s033：隨機(jī)生成一個(gè)服從平均值為0，標(biāo)準(zhǔn)差為ζ的高斯分布1×k的矩陣δ；s034：計(jì)算αp，得到一個(gè)新矩陣y；s035：計(jì)算y+δ，得到一個(gè)新矩陣α′；s036：將新矩陣α′進(jìn)行變維操作，即可得到一個(gè)三維數(shù)組im_jl_array_t，im_jl_array_t[a,b,0]＝α′[0,0]，im_jl_array_t[a,b,1]＝α′[0,1]，…，im_jl_array_t[a,b,k-2]＝α′[0,k-2]，im_jl_array_t[a,b,k-1]＝α′[0,k-1]。作為一種優(yōu)選的技術(shù)方案，在所述的步驟s04中，在本地平臺(tái)，對(duì)n個(gè)圖像的該波段圖像進(jìn)行paillier加密，以第t個(gè)圖像為例，使用已知的公鑰，對(duì)數(shù)組im_block_permutation_array_t進(jìn)行paillier加密，得到一個(gè)二維數(shù)組im_paillier_array_t，該數(shù)組的長(zhǎng)為im_h，寬為im_w，其工作流程如下：s041：選取兩個(gè)隨機(jī)素?cái)?shù)p和q，并根據(jù)n＝pq，計(jì)算n的值；s042：計(jì)算出paillier加密參數(shù)λ和g，得到paillier加密的公鑰為(n,g)，私鑰為λ；s043：計(jì)算gm的值，m∈[0,65535]，并把結(jié)果保存在第一個(gè)表(結(jié)構(gòu)如下)中；mgms044：隨機(jī)生成65536個(gè)隨機(jī)數(shù)r，計(jì)算(gn)r的值，并把結(jié)果保存在第二個(gè)表(結(jié)構(gòu)如下)中；r(gn)r以數(shù)組im_pixel_permutation_array_t中第a行第b列的點(diǎn)位列，對(duì)數(shù)組im_pixel_permutation_array_t中的每一個(gè)點(diǎn)重復(fù)步驟s045-s047；s045：從第一個(gè)表中查找m的值等于im_pixel_permutation_array_t[a,b]對(duì)應(yīng)的gm的值；s046：從第二個(gè)表中隨機(jī)抽取5個(gè)不同的(gn)r值，即(gn)r1，(gn)r2，(gn)r3，(gn)r4，(gn)r5；s047：im_paillier_array_t[a,b]＝gm(gn)r1(gn)r2(gn)r3(gn)r4(gn)r5modn2。作為一種優(yōu)選的技術(shù)方案，在所述的步驟s05中，在云平臺(tái)對(duì)遙感圖像進(jìn)行融合去噪，即對(duì)im_jl_array_t和im_paillier_array_t進(jìn)行非局部去噪，得到二維數(shù)組im_paillier_result_array，該數(shù)組的長(zhǎng)為im_h，寬為im_w，其工作流程如下：s051：計(jì)算im_jl_array_1，im_jl_array_2，…，im_jl_array_n-1，im_jl_array_n這n個(gè)數(shù)組的疊加平均，得到一個(gè)三維數(shù)組，im_jl_mean_array，該數(shù)組的長(zhǎng)為im_h，寬im_w，數(shù)組的高為k，以行a列b高c的這個(gè)點(diǎn)為例，s052：選取h1，h1為非局部去噪的濾波參數(shù)；s053：計(jì)算a的值，a為放大倍數(shù)，計(jì)算公式為a＝27×im_h×im_w；以數(shù)組im_paillier_result_array中第a行第b列的點(diǎn)為例，對(duì)于數(shù)組中的每一個(gè)點(diǎn)重復(fù)步驟s054-s0514；s054：用im_jl_mean_array[a,b,0]，im_jl_mean_array[a,b,1]，…，im_jl_mean_array[a,b,k-2]，im_jl_mean_array[a,b,k-1]組成一個(gè)行向量，記為v(a,b)；以第t個(gè)圖像為例，對(duì)每一個(gè)圖像重復(fù)步驟s055-s057；以數(shù)組im_paillier_array_t中第i行第j列的點(diǎn)為例，對(duì)于數(shù)組im_paillier_array_t中的每一個(gè)點(diǎn)重復(fù)步驟s055-s057；s055：用im_jl_array_t[i,j,0]，im_jl_array_t[i,j,1]，…，im_jl_array_t[i,j,k-2]，im_jl_array_t[i,j,k-1]組成一個(gè)行向量，記為v(i,j)；s056：計(jì)算os_dis_square_t[i,j]，計(jì)算公式為s057：計(jì)算dis_t[i,j]，計(jì)算公式為s058：確定w_number的值；s059：確定數(shù)組os_dis_square_1，os_dis_square_2，…，os_dis_square_n中第w_number個(gè)大的值為dis_square，并計(jì)算h2的值，計(jì)算公式為以第t個(gè)圖像為例，對(duì)每一個(gè)圖像重復(fù)步驟s0510；以以數(shù)組dis_t中第i行第j列的點(diǎn)為例，對(duì)于數(shù)組dis_t中的每一個(gè)點(diǎn)重復(fù)步驟s0510；s0510：若os_dis_square_t[i,j]＞dis_square，則dis_t[i,j]＝0；若os_dis_square_t[i,j]≤dis_square，則s0511：計(jì)算z的值，計(jì)算公式為其中ω為整個(gè)圖像；以第t個(gè)圖像為例，對(duì)每一個(gè)圖像重復(fù)步驟s0512到步驟s0513:以數(shù)組w_t中第i行第j列的點(diǎn)為例，對(duì)于數(shù)組w_t中的每一個(gè)點(diǎn)重復(fù)步驟s0512-s0513；s0512：計(jì)算q_t[i,j]的值，計(jì)算公式為s0513：計(jì)算w_t[i,j]的值，計(jì)算公式為w_t[i,j]＝aq_t[i,j]；s0514：計(jì)算im_paillier_result_array[a,b]的值，計(jì)算公式為其中ω為整個(gè)圖像。作為一種優(yōu)選的技術(shù)方案，在所述的步驟s06中，在本地平臺(tái)，使用私鑰λ對(duì)數(shù)組im_paillier_result_array進(jìn)行paillier解密，得到一個(gè)二維數(shù)組im_permutation_result_array，該數(shù)組的長(zhǎng)為im_h，寬為im_w，其工作流程如下：s061：計(jì)算u的值，計(jì)算公式為以數(shù)組im_paillier_result_array中第a行第b列的點(diǎn)為例，對(duì)數(shù)組im_paillier_result_array中的每一個(gè)點(diǎn)重復(fù)步驟s062-s063；s062：計(jì)算nl的值，計(jì)算公式為其中c等于im_paillier_result_array[a,b]；s063：計(jì)算im_permutation_result_array[a,b]的值，計(jì)算公式為im_permutation_result_array[a,b]＝nl/a。作為一種優(yōu)選的技術(shù)方案，在所述的步驟s07中，在本地平臺(tái)，利用圖像恢復(fù)數(shù)列，對(duì)im_permutation_result_array進(jìn)行遙感圖像的像素置換恢復(fù)，得到一個(gè)二維數(shù)組im_result_array，該數(shù)組的長(zhǎng)為im_h，寬為im_w。作為一種優(yōu)選的技術(shù)方案，在所述的步驟s08中，在本地平臺(tái)，對(duì)完成步驟s02-s07的所有波段得到的結(jié)果進(jìn)行合并處理，得到去噪后的完整遙感圖像。本發(fā)明優(yōu)點(diǎn)在于：1、，在本地平臺(tái)對(duì)遙感圖像進(jìn)行拆分、加密，密鑰保存在本地平臺(tái)，不存在信息泄漏問(wèn)題；2、將針對(duì)同一場(chǎng)景的多幀遙感圖像放在云平臺(tái)進(jìn)行去噪處理，減少本地平臺(tái)數(shù)據(jù)處理量，提高工作效率；3、云平臺(tái)將去噪完成的圖像返回到本地平臺(tái)，本地平臺(tái)進(jìn)行解密、合并，即可得到融合去噪的遙感圖像；4、既高效的完成遙感圖像去噪，又防止遙感圖像信息在外包過(guò)程中出現(xiàn)泄密，有效結(jié)合了云平臺(tái)計(jì)算效率高和本地平臺(tái)安全的優(yōu)點(diǎn)，同時(shí)避免了云平臺(tái)泄密和本地平臺(tái)計(jì)算效率低的缺點(diǎn)。附圖說(shuō)明附圖1是本發(fā)明的架構(gòu)圖。附圖2是本發(fā)明的總流程圖。附圖3是本發(fā)明的具體工作流程圖。附圖4是本發(fā)明的johnson-lindenstrauss轉(zhuǎn)換流程圖。附圖5是本發(fā)明的paillier加密、解密流程圖。具體實(shí)施方式下面結(jié)合附圖對(duì)本發(fā)明提供的具體實(shí)施方式作詳細(xì)說(shuō)明。實(shí)施例1參照?qǐng)D1、2，圖1是本發(fā)明的一種針對(duì)多幀遙感圖像的安全外包融合去噪方法的架構(gòu)圖，圖2是本發(fā)明的一種針對(duì)多幀遙感圖像的安全外包融合去噪方法的流程圖。一種針對(duì)多幀遙感圖像的安全外包融合去噪方法，包括本地平臺(tái)和云平臺(tái)，所述的本地平臺(tái)包括遙感圖像拆分模塊、遙感圖像置換模塊、遙感圖像轉(zhuǎn)換模塊、遙感圖像加密模塊、遙感圖像解密模塊、遙感圖像恢復(fù)排列模塊、遙感圖像合并模塊，所述的云平臺(tái)包括遙感圖像融合去噪模塊，所述的遙感圖像拆分模塊為對(duì)遙感圖像按照?qǐng)D像波段數(shù)進(jìn)行拆分，所述的遙感圖像置換模塊為對(duì)遙感圖像進(jìn)行隨機(jī)置換，所述的遙感圖像轉(zhuǎn)換模塊為對(duì)遙感圖像進(jìn)行johnson-lindenstrauss轉(zhuǎn)換，所述的遙感圖像加密模塊為對(duì)遙感圖像進(jìn)行paillier加密，所述的遙感圖像解密模塊為對(duì)遙感圖像進(jìn)行paillier解密，其中，所述的針對(duì)多幀遙感圖像的安全外包融合去噪方法的工作流程為，s01：遙感圖像拆分；s02：遙感圖像置換；s03：遙感圖像轉(zhuǎn)換；s04：遙感圖像加密；s05：遙感圖像融合去噪；s06：遙感圖像解密；s07：遙感圖像恢復(fù)排列；s08：遙感圖像合并。實(shí)施例2參照?qǐng)D3，圖3是本發(fā)明的一種針對(duì)多幀遙感圖像的安全外包融合去噪方法的具體工作流程圖。本發(fā)明的一種針對(duì)多幀遙感圖像的安全外包融合去噪方法的具體工作流程如下所示。s01：遙感圖像拆分在本地平臺(tái)，讀取針對(duì)同一場(chǎng)景的多幀遙感圖像，遙感圖像的個(gè)數(shù)為n，遙感圖像的長(zhǎng)為im_h，遙感圖像的寬為im_w，遙感圖像的波段數(shù)為im_band_number，按照波段數(shù)把每個(gè)遙感圖像拆分成im_band_number個(gè)圖像，分別讀取個(gè)波段圖像的16位深灰度值，并對(duì)每一個(gè)波段的圖像進(jìn)行步驟s02-s07。s02：遙感圖像置換在本地平臺(tái)，為打亂圖像像素的基本位置，對(duì)n個(gè)圖像的該波段圖像進(jìn)行隨機(jī)置換，以第t個(gè)圖像為例。s021：生成一個(gè)長(zhǎng)度為im_h×im_w的隨機(jī)數(shù)序列，記錄其索引順序i1；s022：對(duì)隨機(jī)數(shù)序列按照從小到大進(jìn)行排列，生成一個(gè)新序列，記錄其索引順序i2；s023：根據(jù)i1→i2，得到圖像置換數(shù)列，根據(jù)i2→i1，得到圖像恢復(fù)數(shù)列；s024：對(duì)遙感圖像進(jìn)行置換，圖像的16位深灰度值數(shù)組為im_array_t，以數(shù)組中第a行第b列的點(diǎn)為例，對(duì)于數(shù)組中的每一點(diǎn)做如下操作：s0241：獲取以im_array_t[a,b]為中心的d×d大小的矩形塊，并將該矩形塊轉(zhuǎn)化成一行，所有的轉(zhuǎn)化后的行組成一個(gè)新數(shù)組；s0242：利用圖像置換序列對(duì)新數(shù)組進(jìn)行行置換；s0243：對(duì)完成行置換的數(shù)組進(jìn)行變維操作，得到一個(gè)三維數(shù)組im_block_permutation_array_t，該數(shù)組的長(zhǎng)為im_h，寬為im_w，高為d2；s025：利用圖像置換序列對(duì)遙感圖像進(jìn)行像素置換，得到一個(gè)二維數(shù)組im_pixel_permutation_array_t，該數(shù)組的長(zhǎng)為im_h，寬為im_w。s03：遙感圖像轉(zhuǎn)換參照?qǐng)D4，圖4是本發(fā)明的一種針對(duì)多幀遙感圖像的安全外包融合去噪方法的johnson-lindenstrauss轉(zhuǎn)換流程圖。在本地平臺(tái)，對(duì)n個(gè)圖像的該波段圖像進(jìn)行johnson-lindenstrauss轉(zhuǎn)換，以第t個(gè)圖像為例，即對(duì)im_block_permutation_array_t進(jìn)行johnson-lindenstrauss轉(zhuǎn)換，將每一個(gè)1×d2的向量壓縮呈一個(gè)1×k的向量，得到一個(gè)三維數(shù)組im_jl_array_t，該數(shù)組的長(zhǎng)為im_h，寬為im_w，高為k。s031：隨機(jī)生成一個(gè)服從平均值為0，方差為1/k的高斯分布的d2×k的矩陣p；以數(shù)組im_array_t中第a行第b列這個(gè)點(diǎn)為例，對(duì)于數(shù)組中的每一個(gè)點(diǎn)重復(fù)步驟s032-s036；s032：將數(shù)組im_block_permutation_array_t中的im_block_permutation_array_t[a,b,0]，im_block_permutation_array_t[a,b,1]，…，im_block_permutation_array_t[a,b,d2-2]，im_block_permutation_array_t[a,b,d2-1]，組成一個(gè)行矩陣α；s033：隨機(jī)生成一個(gè)服從平均值為0，標(biāo)準(zhǔn)差為ζ的高斯分布1×k的矩陣δ；s034：計(jì)算αp，得到一個(gè)新矩陣y；s035：計(jì)算y+δ，得到一個(gè)新矩陣α′；s036：將新矩陣α′進(jìn)行變維操作，即可得到一個(gè)三維數(shù)組im_jl_array_t，im_jl_array_t[a,b,0]＝α′[0,0]，im_jl_array_t[a,b,1]＝α′[0,1]，…，im_jl_array_t[a,b,k-2]＝α′[0,k-2]，im_jl_array_t[a,b,k-1]＝α′[0,k-1]。s04：遙感圖像加密參照?qǐng)D5，圖5是本發(fā)明的paillier加密、解密流程圖。在本地平臺(tái)，對(duì)n個(gè)圖像的該波段圖像進(jìn)行paillier加密。以第t個(gè)圖像為例，使用已知的公鑰，對(duì)數(shù)組im_block_permutation_array_t進(jìn)行paillier加密，得到一個(gè)二維數(shù)組im_paillier_array_t，該數(shù)組的長(zhǎng)為im_h，寬為im_w。s041：選取兩個(gè)隨機(jī)素?cái)?shù)p和q，并根據(jù)n＝pq，計(jì)算n的值；s042：計(jì)算出paillier加密參數(shù)λ和g，得到paillier加密的公鑰為(n,g)，私鑰為λ；s043：計(jì)算gm的值，m∈[0,65535]，并把結(jié)果保存在第一個(gè)表(結(jié)構(gòu)如下)中；mgms044：隨機(jī)生成65536個(gè)隨機(jī)數(shù)r，計(jì)算(gn)r的值，并把結(jié)果保存在第二個(gè)表(結(jié)構(gòu)如下)中；r(gn)r以數(shù)組im_pixel_permutation_array_t中第a行第b列的點(diǎn)位列，對(duì)數(shù)組im_pixel_permutation_array_t中的每一個(gè)點(diǎn)重復(fù)步驟s045-s047；s045：從第一個(gè)表中查找m的值等于im_pixel_permutation_array_t[a,b]對(duì)應(yīng)的gm的值；s046：從第二個(gè)表中隨機(jī)抽取5個(gè)不同的(gn)r值，即(gn)r1，(gn)r2，(gn)r3，(gn)r4，(gn)r5；s047：im_paillier_array_t[a,b]＝gm(gn)r1(gn)r2(gn)r3(gn)r4(gn)r5modn2。s05：遙感圖像融合去噪在云平臺(tái)對(duì)遙感圖像進(jìn)行融合去噪，即對(duì)im_jl_array_t和im_paillier_array_t進(jìn)行非局部去噪，得到二維數(shù)組im_paillier_result_array，該數(shù)組的長(zhǎng)為im_h，寬為im_w。s051：計(jì)算im_jl_array_1，im_jl_array_2，…，im_jl_array_n-1，im_jl_array_n這n個(gè)數(shù)組的疊加平均，得到一個(gè)三維數(shù)組，im_jl_mean_array，該數(shù)組的長(zhǎng)為im_h，寬im_w，數(shù)組的高為k；以行a列b高c的這個(gè)點(diǎn)為例，s052：選取h1，h1為非局部去噪的濾波參數(shù)；s053：計(jì)算a的值，a為放大倍數(shù)，計(jì)算公式為a＝27×im_h×im_w；以數(shù)組im_paillier_result_array中第a行第b列的點(diǎn)為例，對(duì)于數(shù)組中的每一個(gè)點(diǎn)重復(fù)步驟s054-s0514；s054：用im_jl_mean_array[a,b,0]，im_jl_mean_array[a,b,1]，…，im_jl_mean_array[a,b,k-2]，im_jl_mean_array[a,b,k-1]組成一個(gè)行向量，記為v(a,b)；以第t個(gè)圖像為例，對(duì)每一個(gè)圖像重復(fù)步驟s055-s057；以數(shù)組im_paillier_array_t中第i行第j列的點(diǎn)為例，對(duì)于數(shù)組im_paillier_array_t中的每一個(gè)點(diǎn)重復(fù)步驟s055-s057；s055：用im_jl_array_t[i,j,0]，im_jl_array_t[i,j,1]，…，im_jl_array_t[i,j,k-2]，im_jl_array_t[i,j,k-1]組成一個(gè)行向量，記為v(i,j)；s056：計(jì)算os_dis_square_t[i,j]，計(jì)算公式為s057：計(jì)算dis_t[i,j]，計(jì)算公式為s058：確定w_number的值；s059：確定數(shù)組os_dis_square_1，os_dis_square_2，…，os_dis_square_n中第w_number個(gè)大的值為dis_square，并計(jì)算h2的值，計(jì)算公式為以第t個(gè)圖像為例，對(duì)每一個(gè)圖像重復(fù)步驟s0510；以以數(shù)組dis_t中第i行第j列的點(diǎn)為例，對(duì)于數(shù)組dis_t中的每一個(gè)點(diǎn)重復(fù)步驟s0510；s0510：若os_dis_square_t[i,j]＞dis_square，則dis_t[i,j]＝0；若os_dis_square_t[i,j]≤dis_square，則s0511：計(jì)算z的值，計(jì)算公式為其中ω為整個(gè)圖像；以第t個(gè)圖像為例，對(duì)每一個(gè)圖像重復(fù)步驟s0512到步驟s0513:以數(shù)組w_t中第i行第j列的點(diǎn)為例，對(duì)于數(shù)組w_t中的每一個(gè)點(diǎn)重復(fù)步驟s0512-s0513；s0512：計(jì)算q_t[i,j]的值，計(jì)算公式為s0513：計(jì)算w_t[i,j]的值，計(jì)算公式為w_t[i,j]＝aq_t[i,j]；s0514：計(jì)算im_paillier_result_array[a,b]的值，計(jì)算公式為其中ω為整個(gè)圖像。s06：遙感圖像解密參照?qǐng)D5，圖5是本發(fā)明的paillier加密、解密流程圖。在本地平臺(tái)，使用私鑰λ對(duì)數(shù)組im_paillier_result_array進(jìn)行paillier解密，得到一個(gè)二維數(shù)組im_permutation_result_array，該數(shù)組的長(zhǎng)為im_h，寬為im_w。s061：計(jì)算u的值，計(jì)算公式為以數(shù)組im_paillier_result_array中第a行第b列的點(diǎn)為例，對(duì)數(shù)組im_paillier_result_array中的每一個(gè)點(diǎn)重復(fù)步驟s062-s063；s062：計(jì)算nl的值，計(jì)算公式為其中c等于im_paillier_result_array[a,b]；s063：計(jì)算im_permutation_result_array[a,b]的值，計(jì)算公式為im_permutation_result_array[a,b]＝nl/a。s07：遙感圖像恢復(fù)排列在本地平臺(tái)，利用圖像恢復(fù)數(shù)列，對(duì)im_permutation_result_array進(jìn)行遙感圖像的像素置換恢復(fù)，得到一個(gè)二維數(shù)組im_result_array，該數(shù)組的長(zhǎng)為im_h，寬為im_w。s08：遙感圖像合并在本地平臺(tái)，對(duì)完成步驟s02-s07的所有波段得到的結(jié)果進(jìn)行合并處理，得到去噪后的完整遙感圖像。本發(fā)明的一種針對(duì)多幀遙感圖像的安全外包融合去噪方法的優(yōu)點(diǎn)在于，在本地平臺(tái)對(duì)遙感圖像進(jìn)行拆分、加密，密鑰保存在本地平臺(tái)，不存在信息泄漏問(wèn)題；將針對(duì)同一場(chǎng)景的多幀遙感圖像放在云平臺(tái)進(jìn)行去噪處理，減少本地平臺(tái)數(shù)據(jù)處理量，提高工作效率；云平臺(tái)將去噪完成的圖像返回到本地平臺(tái)，本地平臺(tái)進(jìn)行解密、合并，即可得到融合去噪的遙感圖像；既高效的完成遙感圖像去噪，又防止遙感圖像信息在外包過(guò)程中出現(xiàn)泄密，有效結(jié)合了云平臺(tái)計(jì)算效率高和本地平臺(tái)安全的優(yōu)點(diǎn)，同時(shí)避免了云平臺(tái)泄密和本地平臺(tái)計(jì)算效率低的缺點(diǎn)。實(shí)施例3本實(shí)施例為本發(fā)明的一種針對(duì)多幀遙感圖像的安全外包融合去噪方法的具體實(shí)施案例。s01：遙感圖像拆分在本地平臺(tái)，讀取針對(duì)同一場(chǎng)景的多幀遙感圖像，遙感圖像的個(gè)數(shù)為3，遙感圖像的長(zhǎng)為512，遙感圖像的寬為512，遙感圖像的波段數(shù)為11，按照波段數(shù)把每個(gè)遙感圖像拆分成11個(gè)圖像，分別讀取個(gè)波段圖像的16位深灰度值，并對(duì)每一個(gè)波段的圖像進(jìn)行步驟s02-s07。s02：遙感圖像置換在本地平臺(tái)，為打亂圖像像素的基本位置，對(duì)3個(gè)圖像的該波段圖像進(jìn)行隨機(jī)置換，以第t個(gè)圖像為例。s021：生成一個(gè)長(zhǎng)度為512×512的隨機(jī)數(shù)序列，記錄其索引順序i1；s022：對(duì)隨機(jī)數(shù)序列按照從小到大進(jìn)行排列，生成一個(gè)新序列，記錄其索引順序i2；s023：根據(jù)i1→i2，得到圖像置換數(shù)列，根據(jù)i2→i1，得到圖像恢復(fù)數(shù)列；s024：對(duì)遙感圖像進(jìn)行置換，圖像的16位深灰度值數(shù)組為im_array_t，以數(shù)組中第a行第b列的點(diǎn)為例，對(duì)于數(shù)組中的每一點(diǎn)做如下操作：s0241：獲取以im_array_t[a,b]為中心的5×5大小的矩形塊，并將該矩形塊轉(zhuǎn)化成一行，所有的轉(zhuǎn)化后的行組成一個(gè)新數(shù)組；s0242：利用圖像置換序列對(duì)新數(shù)組進(jìn)行行置換；s0243：對(duì)完成行置換的數(shù)組進(jìn)行變維操作，得到一個(gè)三維數(shù)組im_block_permutation_array_t，該數(shù)組的長(zhǎng)為512，寬為512，高為25；s025：利用圖像置換序列對(duì)遙感圖像進(jìn)行像素置換，得到一個(gè)二維數(shù)組im_pixel_permutation_array_t，該數(shù)組的長(zhǎng)為512，寬為512。s03：遙感圖像轉(zhuǎn)換在本地平臺(tái)，對(duì)3個(gè)圖像的該波段圖像進(jìn)行johnson-lindenstrauss轉(zhuǎn)換，以第t個(gè)圖像為例，即對(duì)im_block_permutation_array_t進(jìn)行johnson-lindenstrauss轉(zhuǎn)換，將每一個(gè)1×25的向量壓縮呈一個(gè)1×18的向量，得到一個(gè)三維數(shù)組im_jl_array_t，該數(shù)組的長(zhǎng)為512，寬為512，高為18。s031：隨機(jī)生成一個(gè)服從平均值為0，方差為1/18的高斯分布的25×18的矩陣p；以數(shù)組im_array_t中第a行第b列這個(gè)點(diǎn)為例，對(duì)于數(shù)組中的每一個(gè)點(diǎn)重復(fù)步驟s032-s036；s032：將數(shù)組im_block_permutation_array_t中的im_block_permutation_array_t[a,b,0]，im_block_permutation_array_t[a,b,1]，…，im_block_permutation_array_t[a,b,23]，im_block_permutation_array_t[a,b,24]，組成一個(gè)行矩陣α；s033：隨機(jī)生成一個(gè)服從平均值為0，標(biāo)準(zhǔn)差為0.5的高斯分布1×18的矩陣δ；s034：計(jì)算αp，得到一個(gè)新矩陣y；s035：計(jì)算y+δ，得到一個(gè)新矩陣α′；s036：將新矩陣α′進(jìn)行變維操作，即可得到一個(gè)三維數(shù)組im_jl_array_t，im_jl_array_t[a,b,0]＝α′[0,0]，im_jl_array_t[a,b,1]＝α′[0,1]，…，im_jl_array_t[a,b,16]＝α′[0,16]，im_jl_array_t[a,b,17]＝α′[0,17]。s04：遙感圖像加密在本地平臺(tái)，對(duì)3個(gè)圖像的該波段圖像進(jìn)行paillier加密。以第t個(gè)圖像為例，使用已知的公鑰，對(duì)數(shù)組im_block_permutation_array_t進(jìn)行paillier加密，得到一個(gè)二維數(shù)組im_paillier_array_t，該數(shù)組的長(zhǎng)為512，寬為512。s041：選取兩個(gè)隨機(jī)素?cái)?shù)p和q，并根據(jù)n＝pq，計(jì)算n的值；s042：計(jì)算出paillier加密參數(shù)λ和g，得到paillier加密的公鑰為(n,g)，私鑰為λ；s043：計(jì)算gm的值，m∈[0,65535]，并把結(jié)果保存在第一個(gè)表(結(jié)構(gòu)如下)中；mgms044：隨機(jī)生成65536個(gè)隨機(jī)數(shù)r，計(jì)算(gn)r的值，并把結(jié)果保存在第二個(gè)表(結(jié)構(gòu)如下)中；r(gn)r以數(shù)組im_pixel_permutation_array_t中第a行第b列的點(diǎn)位列，對(duì)數(shù)組im_pixel_permutation_array_t中的每一個(gè)點(diǎn)重復(fù)步驟s045-s047；s045：從第一個(gè)表中查找m的值等于im_pixel_permutation_array_t[a,b]對(duì)應(yīng)的gm的值；s046：從第二個(gè)表中隨機(jī)抽取5個(gè)不同的(gn)r值，即(gn)r1，(gn)r2，(gn)r3，(gn)r4，(gn)r5；s047：im_paillier_array_t[a,b]＝gm(gn)r1(gn)r2(gn)r3(gn)r4(gn)r5modn2。s05：遙感圖像融合去噪在云平臺(tái)對(duì)遙感圖像進(jìn)行融合去噪，即對(duì)im_jl_array_t和im_paillier_array_t進(jìn)行非局部去噪，得到二維數(shù)組im_paillier_result_array，該數(shù)組的長(zhǎng)為512，寬為512。s051：計(jì)算im_jl_array_1，im_jl_array_2，…，im_jl_array_n-1，im_jl_array_n這n個(gè)數(shù)組的疊加平均，得到一個(gè)三維數(shù)組，im_jl_mean_array，該數(shù)組的長(zhǎng)為512，寬512，數(shù)組的高為18；以行a列b高c的這個(gè)點(diǎn)為例，s052：選取h1，h1為非局部去噪的濾波參數(shù)；s053：計(jì)算a的值，a為放大倍數(shù)，計(jì)算公式為a＝27×512×512；以數(shù)組im_paillier_result_array中第a行第b列的點(diǎn)為例，對(duì)于數(shù)組中的每一個(gè)點(diǎn)重復(fù)步驟s054-s0514；s054：用im_jl_mean_array[a,b,0]，im_jl_mean_array[a,b,1]，…，im_jl_mean_array[a,b,16]，im_jl_mean_array[a,b,17]組成一個(gè)行向量，記為v(a,b)；以第t個(gè)圖像為例，對(duì)每一個(gè)圖像重復(fù)步驟s055-s057；以數(shù)組im_paillier_array_t中第i行第j列的點(diǎn)為例，對(duì)于數(shù)組im_paillier_array_t中的每一個(gè)點(diǎn)重復(fù)步驟s055-s057；s055：用im_jl_array_t[i,j,0]，im_jl_array_t[i,j,1]，…，im_jl_array_t[i,j,16]，im_jl_array_t[i,j,17]組成一個(gè)行向量，記為v(i,j)；s056：計(jì)算os_dis_square_t[i,j]，計(jì)算公式為s057：計(jì)算dis_t[i,j]，計(jì)算公式為s058：確定w_number的值為128；s059：確定數(shù)組os_dis_square_1，os_dis_square_2，…，os_dis_square_n中第w_number個(gè)大的值為dis_square，并計(jì)算h2的值，計(jì)算公式為以第t個(gè)圖像為例，對(duì)每一個(gè)圖像重復(fù)步驟s0510；以以數(shù)組dis_t中第i行第j列的點(diǎn)為例，對(duì)于數(shù)組dis_t中的每一個(gè)點(diǎn)重復(fù)步驟s0510；s0510：若os_dis_square_t[i,j]＞dis_square，則dis_t[i,j]＝0；若os_dis_square_t[i,j]≤dis_square，則s0511：計(jì)算z的值，計(jì)算公式為其中ω為整個(gè)圖像；以第t個(gè)圖像為例，對(duì)每一個(gè)圖像重復(fù)步驟s0512到步驟s0513:以數(shù)組w_t中第i行第j列的點(diǎn)為例，對(duì)于數(shù)組w_t中的每一個(gè)點(diǎn)重復(fù)步驟s0512-s0513；s0512：計(jì)算q_t[i,j]的值，計(jì)算公式為s0513：計(jì)算w_t[i,j]的值，計(jì)算公式為w_t[i,j]＝aq_t[i,j]；s0514：計(jì)算im_paillier_result_array[a,b]的值，計(jì)算公式為其中ω為整個(gè)圖像。s06：遙感圖像解密在本地平臺(tái)，使用私鑰λ對(duì)數(shù)組im_paillier_result_array進(jìn)行paillier解密，得到一個(gè)二維數(shù)組im_permutation_result_array，該數(shù)組的長(zhǎng)為512，寬為512。s061：計(jì)算u的值，計(jì)算公式為以數(shù)組im_paillier_result_array中第a行第b列的點(diǎn)為例，對(duì)數(shù)組im_paillier_result_array中的每一個(gè)點(diǎn)重復(fù)步驟s062-s063；s062：計(jì)算nl的值，計(jì)算公式為其中c等于im_paillier_result_array[a,b]；s063：計(jì)算im_permutation_result_array[a,b]的值，計(jì)算公式為im_permutation_result_array[a,b]＝nl/a。s07：遙感圖像恢復(fù)排列在本地平臺(tái)，利用圖像恢復(fù)數(shù)列，對(duì)im_permutation_result_array進(jìn)行遙感圖像的像素置換恢復(fù)，得到一個(gè)二維數(shù)組im_result_array，該數(shù)組的長(zhǎng)為512，寬為512。s08：遙感圖像合并在本地平臺(tái)，對(duì)完成步驟s02-s07的所有波段得到的結(jié)果進(jìn)行合并處理，得到去噪后的完整遙感圖像。實(shí)施例4本實(shí)施例為本發(fā)明的一種針對(duì)多幀遙感圖像的安全外包融合去噪方法的具體實(shí)施案例。s01：遙感圖像拆分在本地平臺(tái)，讀取針對(duì)同一場(chǎng)景的多幀遙感圖像，遙感圖像的個(gè)數(shù)為5，遙感圖像的長(zhǎng)為256，遙感圖像的寬為256，遙感圖像的波段數(shù)為8，按照波段數(shù)把每個(gè)遙感圖像拆分成8個(gè)圖像，分別讀取個(gè)波段圖像的16位深灰度值，并對(duì)每一個(gè)波段的圖像進(jìn)行步驟s02-s07。s02：遙感圖像置換在本地平臺(tái)，為打亂圖像像素的基本位置，對(duì)5個(gè)圖像的該波段圖像進(jìn)行隨機(jī)置換，以第t個(gè)圖像為例。s021：生成一個(gè)長(zhǎng)度為256×256的隨機(jī)數(shù)序列，記錄其索引順序i1；s022：對(duì)隨機(jī)數(shù)序列按照從小到大進(jìn)行排列，生成一個(gè)新序列，記錄其索引順序i2；s023：根據(jù)i1→i2，得到圖像置換數(shù)列，根據(jù)i2→i1，得到圖像恢復(fù)數(shù)列；s024：對(duì)遙感圖像進(jìn)行置換，圖像的16位深灰度值數(shù)組為im_array_t，以數(shù)組中第a行第b列的點(diǎn)為例，對(duì)于數(shù)組中的每一點(diǎn)做如下操作：s0241：獲取以im_array_t[a,b]為中心的3×3大小的矩形塊，并將該矩形塊轉(zhuǎn)化成一行，所有的轉(zhuǎn)化后的行組成一個(gè)新數(shù)組；s0242：利用圖像置換序列對(duì)新數(shù)組進(jìn)行行置換；s0243：對(duì)完成行置換的數(shù)組進(jìn)行變維操作，得到一個(gè)三維數(shù)組im_block_permutation_array_t，該數(shù)組的長(zhǎng)為256，寬為256，高為9；s025：利用圖像置換序列對(duì)遙感圖像進(jìn)行像素置換，得到一個(gè)二維數(shù)組im_pixel_permutation_array_t，該數(shù)組的長(zhǎng)為256，寬為256。s03：遙感圖像轉(zhuǎn)換在本地平臺(tái)，對(duì)5個(gè)圖像的該波段圖像進(jìn)行johnson-lindenstrauss轉(zhuǎn)換，以第t個(gè)圖像為例，即對(duì)im_block_permutation_array_t進(jìn)行johnson-lindenstrauss轉(zhuǎn)換，將每一個(gè)1×9的向量壓縮呈一個(gè)1×4的向量，得到一個(gè)三維數(shù)組im_jl_array_t，該數(shù)組的長(zhǎng)為256，寬為256，高為4。s031：隨機(jī)生成一個(gè)服從平均值為0，方差為1/4的高斯分布的9×4的矩陣p；以數(shù)組im_array_t中第a行第b列這個(gè)點(diǎn)為例，對(duì)于數(shù)組中的每一個(gè)點(diǎn)重復(fù)步驟s032-s036；s032：將數(shù)組im_block_permutation_array_t中的im_block_permutation_array_t[a,b,0]，im_block_permutation_array_t[a,b,1]，…，im_block_permutation_array_t[a,b,7]，im_block_permutation_array_t[a,b,8]，組成一個(gè)行矩陣α；s033：隨機(jī)生成一個(gè)服從平均值為0，標(biāo)準(zhǔn)差為0.3的高斯分布1×4的矩陣δ；s034：計(jì)算αp，得到一個(gè)新矩陣y；s035：計(jì)算y+δ，得到一個(gè)新矩陣α′；s036：將新矩陣α′進(jìn)行變維操作，即可得到一個(gè)三維數(shù)組im_jl_array_t，im_jl_array_t[a,b,0]＝α′[0,0]，im_jl_array_t[a,b,1]＝α′[0,1]，im_jl_array_t[a,b,2]＝α′[0,2]，im_jl_array_t[a,b,3]＝α′[0,3]。s04：遙感圖像加密在本地平臺(tái)，對(duì)5個(gè)圖像的該波段圖像進(jìn)行paillier加密。以第t個(gè)圖像為例，使用已知的公鑰，對(duì)數(shù)組im_block_permutation_array_t進(jìn)行paillier加密，得到一個(gè)二維數(shù)組im_paillier_array_t，該數(shù)組的長(zhǎng)為256，寬為256。s041：選取兩個(gè)隨機(jī)素?cái)?shù)p和q，并根據(jù)n＝pq，計(jì)算n的值；s042：計(jì)算出paillier加密參數(shù)λ和g，得到paillier加密的公鑰為(n,g)，私鑰為λ；s043：計(jì)算gm的值，m∈[0,65535]，并把結(jié)果保存在第一個(gè)表(結(jié)構(gòu)如下)中；mgms044：隨機(jī)生成65536個(gè)隨機(jī)數(shù)r，計(jì)算(gn)r的值，并把結(jié)果保存在第二個(gè)表(結(jié)構(gòu)如下)中；r(gn)r以數(shù)組im_pixel_permutation_array_t中第a行第b列的點(diǎn)位列，對(duì)數(shù)組im_pixel_permutation_array_t中的每一個(gè)點(diǎn)重復(fù)步驟s045-s047；s045：從第一個(gè)表中查找m的值等于im_pixel_permutation_array_t[a,b]對(duì)應(yīng)的gm的值；s046：從第二個(gè)表中隨機(jī)抽取5個(gè)不同的(gn)r值，即(gn)r1，(gn)r2，(gn)r3，(gn)r4，(gn)r5；s047：im_paillier_array_t[a,b]＝gm(gn)r1(gn)r2(gn)r3(gn)r4(gn)r5modn2。s05：遙感圖像融合去噪在云平臺(tái)對(duì)遙感圖像進(jìn)行融合去噪，即對(duì)im_jl_array_t和im_paillier_array_t進(jìn)行非局部去噪，得到二維數(shù)組im_paillier_result_array，該數(shù)組的長(zhǎng)為256，寬為256。s051：計(jì)算im_jl_array_1，im_jl_array_2，…，im_jl_array_n-1，im_jl_array_n這5個(gè)數(shù)組的疊加平均，得到一個(gè)三維數(shù)組，im_jl_mean_array，該數(shù)組的長(zhǎng)為256，寬256，數(shù)組的高為4；以行a列b高c的這個(gè)點(diǎn)為例，s052：選取h1，h1為非局部去噪的濾波參數(shù)；s053：計(jì)算a的值，a為放大倍數(shù)，計(jì)算公式為a＝27×256×256；以數(shù)組im_paillier_result_array中第a行第b列的點(diǎn)為例，對(duì)于數(shù)組中的每一個(gè)點(diǎn)重復(fù)步驟s054-s0514；s054：用im_jl_mean_array[a,b,0]，im_jl_mean_array[a,b,1]，im_jl_mean_array[a,b,5]，im_jl_mean_array[a,b,3]組成一個(gè)行向量，記為v(a,b)；以第t個(gè)圖像為例，對(duì)每一個(gè)圖像重復(fù)步驟s055-s057；以數(shù)組im_paillier_array_t中第i行第j列的點(diǎn)為例，對(duì)于數(shù)組im_paillier_array_t中的每一個(gè)點(diǎn)重復(fù)步驟s055-s057；s055：用im_jl_array_t[i,j,0]，im_jl_array_t[i,j,1]，im_jl_array_t[i,j,3]，im_jl_array_t[i,j,3]組成一個(gè)行向量，記為v(i,j)；s056：計(jì)算os_dis_square_t[i,j]，計(jì)算公式為s057：計(jì)算dis_t[i,j]，計(jì)算公式為s058：確定w_number的值為128；s059：確定數(shù)組os_dis_square_1，os_dis_square_2，…，os_dis_square_n中第w_number個(gè)大的值為dis_square，并計(jì)算h2的值，計(jì)算公式為以第t個(gè)圖像為例，對(duì)每一個(gè)圖像重復(fù)步驟s0510；以以數(shù)組dis_t中第i行第j列的點(diǎn)為例，對(duì)于數(shù)組dis_t中的每一個(gè)點(diǎn)重復(fù)步驟s0510；s0510：若os_dis_square_t[i,j]＞dis_square，則dis_t[i,j]＝0；若os_dis_square_t[i,j]≤dis_square，則s0511：計(jì)算z的值，計(jì)算公式為其中ω為整個(gè)圖像；以第t個(gè)圖像為例，對(duì)每一個(gè)圖像重復(fù)步驟s0512到步驟s0513:以數(shù)組w_t中第i行第j列的點(diǎn)為例，對(duì)于數(shù)組w_t中的每一個(gè)點(diǎn)重復(fù)步驟s0512-s0513；s0512：計(jì)算q_t[i,j]的值，計(jì)算公式為s0513：計(jì)算w_t[i,j]的值，計(jì)算公式為w_t[i,j]＝aq_t[i,j]；s0514：計(jì)算im_paillier_result_array[a,b]的值，計(jì)算公式為其中ω為整個(gè)圖像。s06：遙感圖像解密在本地平臺(tái)，使用私鑰λ對(duì)數(shù)組im_paillier_result_array進(jìn)行paillier解密，得到一個(gè)二維數(shù)組im_permutation_result_array，該數(shù)組的長(zhǎng)為256，寬為256。s061：計(jì)算u的值，計(jì)算公式為以數(shù)組im_paillier_result_array中第a行第b列的點(diǎn)為例，對(duì)數(shù)組im_paillier_result_array中的每一個(gè)點(diǎn)重復(fù)步驟s062-s063；s062：計(jì)算nl的值，計(jì)算公式為其中c等于im_paillier_result_array[a,b]；s063：計(jì)算im_permutation_result_array[a,b]的值，計(jì)算公式為im_permutation_result_array[a,b]＝nl/a。s07：遙感圖像恢復(fù)排列在本地平臺(tái)，利用圖像恢復(fù)數(shù)列，對(duì)im_permutation_result_array進(jìn)行遙感圖像的像素置換恢復(fù)，得到一個(gè)二維數(shù)組im_result_array，該數(shù)組的長(zhǎng)為256，寬為256。s08：遙感圖像合并在本地平臺(tái)，對(duì)完成步驟s02-s07的所有波段得到的結(jié)果進(jìn)行合并處理，得到去噪后的完整遙感圖像。以上所述僅是本發(fā)明的優(yōu)選實(shí)施方式，應(yīng)當(dāng)指出，對(duì)于本
技術(shù)領(lǐng)域：
的普通技術(shù)人員，在不脫離本發(fā)明方法的前提下，還可以做出若干改進(jìn)和補(bǔ)充，這些改進(jìn)和補(bǔ)充也應(yīng)視為本發(fā)明的保護(hù)范圍。當(dāng)前第1頁(yè)12