歡迎來到海神之光博客之家?博主簡介熱愛科研的Matlab仿真開發(fā)者修心和技術同步精進個人主頁海神之光代碼獲取方式海神之光Matlab王者學習之路—代碼獲取方式??座右銘行百里者半于九十。更多Matlab圖像處理仿真內(nèi)容點擊①Matlab圖像處理進階版②付費專欄Matlab圖像處理初級版??關注CSDN海神之光更多資源等你來?一、DCT數(shù)字水印嵌入與提取簡介1 基本DCT變換目前基于DCT域的水印方法已經(jīng)成為數(shù)字水印算法研究的熱點它的核心思想就是通過離散傅立葉變換對圖像塊進行處理后再選擇變換域中的一些系數(shù)值依據(jù)一定規(guī)則來嵌入水印。由于圖像塊中DCT系數(shù)頻帶分布由左上角的直流分量DC往下對應的系數(shù)頻率由低頻升至高頻因此在不影響原圖質(zhì)量的前提下可將水印信息根據(jù)能量大小嵌入相應系數(shù)頻帶中。通過圖像塊量化與水印嵌入結(jié)合的處理方法將水印信息均勻分布在圖像的整個空間域在圖像裁剪和濾波方面變換域的水印比在空間域的更能表現(xiàn)出一定的魯棒性。2 水印算法描述2.1 水印嵌入算法該算法采用加性嵌入的方式在經(jīng)過DCT變換后的子圖像塊的中頻域中選取隱秘位置嵌入水印信息具體的嵌入流程如下圖1所示圖1 分塊水印嵌入流程(1分塊處理設宿主圖像為P將其分塊處理為8*8的K個子塊。(2水印預處理設水印圖像為W對其進行互補變換變換后的水印圖像和變換前的水印圖像相互補。(3對水印圖像進行Arnold置亂變換并依據(jù)混沌映射規(guī)則選取密鑰混沌序列并與水印序列異或運算將置換次數(shù)和異或運算處理后的結(jié)果分別作為水印嵌入算法的密鑰1和密鑰2。(4)DCT變換對各子塊內(nèi)做DCT變換利用zig-zag對DCT系數(shù)進行掃描得到第k塊子圖像塊的序列為Zk(i),i0,1,2…63.(5水印嵌入算法依據(jù)zig-zag排序在各子塊的中、低頻段選取特定系數(shù)x(m和x(n在系數(shù)坐標a,b和c,d處嵌入水印信息圖像W并將其作為密鑰3。同理嵌入互補水印圖片W’并將嵌入的位置作為密鑰4。水印嵌入的方法如下(6)IDCT變換將每一個子圖像塊作二維DCT逆變換。(7子塊合并將每一個子塊合并成嵌入水印的圖像P’。2.2 水印提取算法將嵌入水印的圖像P’分塊處理并對各子塊進行二維DCT變換由密鑰3和4推斷所選擇的水印系數(shù)若x(m≤x(n則水印信息為0若x(m)x(n則水印信息為1再利用密鑰1和2將初步水印的信息解密再進行Arnold逆變換最終提取出水印信息。2.3 水印檢測算法本文通過計算峰值信噪比PSNR的值評價嵌入水印的宿主圖像的質(zhì)量一幅m和n的圖像PSNR度量標準定義為歸一化相關系數(shù)NC的值判斷嵌入水印的圖像與宿主圖像的相似度其定義為?二、部分源代碼close all;clear all;clc;I imread(‘pepper.bmp’);J rgb2gray(I);K_1 fft2(J);L_1 abs(K_1/256);K_2 fftshift(K_1);L_2 abs(K_2/256);figure;subplot(2,2,1);imshow(I); title(‘原始圖像’);subplot(2,2,2);imshow(J); title(‘灰度圖像’);subplot(2,2,3);imshow(uint8(L_1)); title(‘灰度圖像傅里葉頻譜’);subplot(2,2,4);imshow(uint8(L_2)); title(‘平移后的頻譜’);其運行結(jié)果如下圖  圖1 圖像傅里葉變換例 2.2 離散余弦變換DCT 離散余弦變換(DCT)是一組不同頻率和幅值的余弦函數(shù)和來近似一副圖像實際上是傅里葉變換的實數(shù)部分。由于離散余弦變量對于一副圖像其大部分可視化信息都集中在少數(shù)的變換系數(shù)上。因此離散余弦變量是數(shù)據(jù)壓縮常用的一個變換編碼方法它能將高相關數(shù)據(jù)能量集中使得它非常適用于圖像壓縮例如國際壓縮標準的JPEG格式中就采用了離散余弦變換。 在傅立葉變換過程中如果被展開的函數(shù)是實偶函數(shù)那么其傅立葉變換中只包含余弦項基于傅立葉變換的這一特點人們提出了離散余弦變換。DCT變換先將圖像函數(shù)變換成偶函數(shù)形式再對其進行二維離散傅立葉變換因此DCT變換可以看成是一種簡化的傅立葉變換。離散余弦變換葉分為一維離散余弦變換和二維離散余弦變換在圖像處理中用到二維變化所以介紹下其二維定義其他定義可查閱《數(shù)字圖像處理》——岡薩雷斯。其二維離散余弦變換定義如下式  c close all; clear all; clc; I imread (boats.bmp); J dct2(I); figure(Name,離散余弦變換); subplot(1,2,1); imshow(I); title(原始圖像); subplot(1,2,2); imshow(log(abs(J))); title(離散余弦變換系數(shù)圖像);其運行結(jié)果如下圖2圖2 離散余弦變換例2.3 小波變換DWT 小波變化是對傅里葉變換和短時傅里葉變換的一個突破其改變就在于將無限長的三角函數(shù)基換成了有限長的會衰減的小波。小波變化的理論較多更多內(nèi)容可查閱《數(shù)字圖像處理》——岡薩雷斯相關章節(jié)內(nèi)容以下從一個叫簡單的角度來解釋小波變換。 表1 小波變換母小波function varargout manit1(varargin)% MANIT1 M-file for manit1.fig% MANIT1, by itself, creates a new MANIT1 or raises the existing% singleton*.%% H MANIT1 returns the handle to a new MANIT1 or the handle to% the existing singleton*.%% MANIT1(‘CALLBACK’,hObject,eventData,handles,…) calls the local% function named CALLBACK in MANIT1.M with the given input arguments.%% MANIT1(‘Property’,‘Value’,…) creates a new MANIT1 or raises the% existing singleton*. Starting from the left, property value pairs are% applied to the GUI before manit1_OpeningFcn gets called. An% unrecognized property name or invalid value makes property application% stop. All inputs are passed to manit1_OpeningFcn via varargin.%% *See GUI Options on GUIDE’s Tools menu. Choose “GUI allows only one% instance to run (singleton)”.%% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help manit1% Last Modified by GUIDE v2.5 14-Jun-2015 10:32:30% Begin initialization code - DO NOT EDITgui_Singleton 1;gui_State struct(‘gui_Name’, mfilename, …‘gui_Singleton’, gui_Singleton, …‘gui_OpeningFcn’, manit1_OpeningFcn, …‘gui_OutputFcn’, manit1_OutputFcn, …‘gui_LayoutFcn’, [] , …‘gui_Callback’, []);if nargin ischar(varargin{1})gui_State.gui_Callback str2func(varargin{1});endif nargout[varargout{1:nargout}] gui_mainfcn(gui_State, varargin{:});elsegui_mainfcn(gui_State, varargin{:});end% End initialization code - DO NOT EDIT% — Executes just before manit1 is made visible.function manit1_OpeningFcn(hObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% varargin command line arguments to manit1 (see VARARGIN)% Choose default command line output for manit1handles.output hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes manit1 wait for user response (see UIRESUME)% uiwait(handles.figure1);% — Outputs from this function are returned to the command line.function varargout manit1_OutputFcn(hObject, eventdata, handles)% varargout cell array for returning output args (see VARARGOUT);% hObject handle to figure% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Get default command line output from handles structurevarargout{1} handles.output;% — Executes on button press in browsecoverimg.function browsecoverimg_Callback(hObject, eventdata, handles)% hObject handle to browsecoverimg (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)[filename, pathname] uigetfile({‘.jpg’;.bmp’;‘.gif’;.*’}, ‘Pick an Image File’);S imread([pathname,filename]);Simresize(S,[512,512]);axes(handles.axes1)imshow(S)title(‘Cover Image’)set(handles.text3,‘string’,filename)handles.SS;guidata(hObject,handles)% — Executes on button press in browsemsg.function browsemsg_Callback(hObject, eventdata, handles)% hObject handle to browsemsg (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)[filename, pathname] uigetfile({‘.bmp’;.jpg’;‘.gif’;.*’}, ‘Pick an Image File’);msg imread([pathname,filename]);axes(handles.axes2)imshow(msg)title(‘Input Message’)set(handles.text4,‘string’,filename)handles.msgmsg;guidata(hObject,handles)% — Executes on button press in DWT.function DWT_Callback(hObject, eventdata, handles)% hObject handle to DWT (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)messagehandles.msg;cover_objecthandles.S;k10;[watermrkd_img,PSNR,IF,NCC,recmessage]dwt(cover_object,message,k);axes(handles.axes3)imshow(watermrkd_img)title(‘Watermarked Image’)axes(handles.axes4)imshow(recmessage)title(‘Recovered Message’)a[PSNR,NCC,IF]‘;thandles.uitable1;set(t,‘Data’,a)handles.aa;guidata(hObject,handles)% — Executes on button press in DWTDCT.function DWTDCT_Callback(hObject, eventdata, handles)% hObject handle to DWTDCT (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)ahandles.a;messagehandles.msg;cover_objecthandles.S;k10;[watermrkd_img,recmessage,PSNR,IF,NCC1] dwtdct(cover_object,message,k);axes(handles.axes3)imshow(watermrkd_img)title(‘DWTDCT Watermarked Image’)axes(handles.axes4)imshow(recmessage)title(‘DWTDCT Recovered Message’)b[PSNR,NCC1,IF]’;thandles.uitable1;set(t,‘Data’,[a b])handles.bb;guidata(hObject,handles)% — Executes on button press in DWTDCTBFO.function DWTDCTBFO_Callback(hObject, eventdata, handles)% hObject handle to DWTDCTBFO (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)ahandles.a;bhandles.b;messagehandles.msg;cover_objecthandles.S;[watermrkd_img,recmessage,PSNR,IF,NCC,pbest] BG(cover_object,message);axes(handles.axes3)imshow(watermrkd_img)title(‘DWTDCTBFO Watermarked Image’)axes(handles.axes4)imshow(recmessage)title(‘DWTDCTBFO Recovered Message’)c[PSNR,NCC,IF]‘;thandles.uitable1;set(t,‘Data’,[a b c])handles.cc;guidata(hObject,handles)% — Executes on button press in pushbutton6.function pushbutton6_Callback(hObject, eventdata, handles)% hObject handle to pushbutton6 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)ahandles.a;bhandles.b;chandles.c;dhandles.d;PSNR[a(1),b(1),c(1),d(1) ];NCC[a(2),b(2),c(2),d(2)];IF[a(3),b(3),c(3),d(3)];figurebar(PSNR)set(gca,‘XTickLabel’,’ ‘);ylabel(‘PSNR’)text(1,0,‘DWT’,‘Rotation’,260,‘Fontsize’,8);text(2,0,‘DWTDCT’,‘Rotation’,260,‘Fontsize’,8);text(3,0,‘DWTDCTBFO’,‘Rotation’,260,‘Fontsize’,8);text(4,0,‘DWTDCTPBFO’,‘Rotation’,260,‘Fontsize’,8);[t]get(gca,‘position’);set(gca,‘position’,[t(1) 0.31 t(3) 0.65])title(‘Bar Graph Comparison of PSNR’)%%%%figurebar(NCC)set(gca,‘XTickLabel’,’ ‘);ylabel(‘NCC’)text(1,0,‘DWT’,‘Rotation’,260,‘Fontsize’,8);text(2,0,‘DWTDCT’,‘Rotation’,260,‘Fontsize’,8);text(3,0,‘DWTDCTBFO’,‘Rotation’,260,‘Fontsize’,8);text(4,0,‘DWTDCTPBFO’,‘Rotation’,260,‘Fontsize’,8);[t]get(gca,‘position’);set(gca,‘position’,[t(1) 0.31 t(3) 0.65])title(‘Bar Graph Comparison of NCC’)% — Executes on button press in DWTDCTPBFO.function DWTDCTPBFO_Callback(hObject, eventdata, handles)% hObject handle to DWTDCTPBFO (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)ahandles.a;bhandles.b;chandles.c;messagehandles.msg;cover_objecthandles.S;[watermrkd_img,recmessage,PSNR,IF,NCC,pbest] BG_PSO(cover_object,message);axes(handles.axes3)imshow(watermrkd_img)title(‘DWTDCTPBFO Watermarked Image’)axes(handles.axes4)imshow(recmessage)title(‘DWTDCTPBFO Recovered Message’)d[PSNR,NCC,IF]’;thandles.uitable1;set(t,‘Data’,[a b c d])handles.dd;guidata(hObject,handles)?三、運行結(jié)果?四、matlab版本及參考文獻1 matlab版本2014a2 參考文獻[1]萬誼丹.基于Arnold和DCT的抗剪切攻擊圖像水印研究[J].網(wǎng)絡安全技術與應用. 2021,(08)3 備注簡介此部分摘自互聯(lián)網(wǎng)僅供參考若侵權(quán)聯(lián)系刪除 仿真咨詢1 各類智能優(yōu)化算法改進及應用生產(chǎn)調(diào)度、經(jīng)濟調(diào)度、裝配線調(diào)度、充電優(yōu)化、車間調(diào)度、發(fā)車優(yōu)化、水庫調(diào)度、三維裝箱、物流選址、貨位優(yōu)化、公交排班優(yōu)化、充電樁布局優(yōu)化、車間布局優(yōu)化、集裝箱船配載優(yōu)化、水泵組合優(yōu)化、解醫(yī)療資源分配優(yōu)化、設施布局優(yōu)化、可視域基站和無人機選址優(yōu)化2 機器學習和深度學習方面卷積神經(jīng)網(wǎng)絡CNN、LSTM、支持向量機SVM、最小二乘支持向量機LSSVM、極限學習機ELM、核極限學習機KELM、BP、RBF、寬度學習、DBN、RF、RBF、DELM、XGBOOST、TCN實現(xiàn)風電預測、光伏預測、電池壽命預測、輻射源識別、交通流預測、負荷預測、股價預測、PM2.5濃度預測、電池健康狀態(tài)預測、水體光學參數(shù)反演、NLOS信號識別、地鐵停車精準預測、變壓器故障診斷3 圖像處理方面圖像識別、圖像分割、圖像檢測、圖像隱藏、圖像配準、圖像拼接、圖像融合、圖像增強、圖像壓縮感知4 路徑規(guī)劃方面旅行商問題TSP、車輛路徑問題VRP、MVRP、CVRP、VRPTW等、無人機三維路徑規(guī)劃、無人機協(xié)同、無人機編隊、機器人路徑規(guī)劃、柵格地圖路徑規(guī)劃、多式聯(lián)運運輸問題、車輛協(xié)同無人機路徑規(guī)劃、天線線性陣列分布優(yōu)化、車間布局優(yōu)化5 無人機應用方面無人機路徑規(guī)劃、無人機控制、無人機編隊、無人機協(xié)同、無人機任務分配6 無線傳感器定位及布局方面?zhèn)鞲衅鞑渴饍?yōu)化、通信協(xié)議優(yōu)化、路由優(yōu)化、目標定位優(yōu)化、Dv-Hop定位優(yōu)化、Leach協(xié)議優(yōu)化、WSN覆蓋優(yōu)化、組播優(yōu)化、RSSI定位優(yōu)化7 信號處理方面信號識別、信號加密、信號去噪、信號增強、雷達信號處理、信號水印嵌入提取、肌電信號、腦電信號、信號配時優(yōu)化8 電力系統(tǒng)方面微電網(wǎng)優(yōu)化、無功優(yōu)化、配電網(wǎng)重構(gòu)、儲能配置9 元胞自動機方面交通流 人群疏散 病毒擴散 晶體生長10 雷達方面卡爾曼濾波跟蹤、航跡關聯(lián)、航跡融合