歡迎來(lái)到海神之光博客之家?博主簡(jiǎn)介熱愛(ài)科研的Matlab仿真開(kāi)發(fā)者修心和技術(shù)同步精進(jìn)個(gè)人主頁(yè)海神之光代碼獲取方式海神之光Matlab王者學(xué)習(xí)之路—代碼獲取方式??座右銘行百里者半于九十。更多Matlab圖像處理仿真內(nèi)容點(diǎn)擊①M(fèi)atlab圖像處理進(jìn)階版②付費(fèi)專(zhuān)欄Matlab圖像處理初級(jí)版??關(guān)注CSDN海神之光更多資源等你來(lái)?一、DCT數(shù)字水印嵌入與提取簡(jiǎn)介1 基本DCT變換目前基于DCT域的水印方法已經(jīng)成為數(shù)字水印算法研究的熱點(diǎn)它的核心思想就是通過(guò)離散傅立葉變換對(duì)圖像塊進(jìn)行處理后再選擇變換域中的一些系數(shù)值依據(jù)一定規(guī)則來(lái)嵌入水印。由于圖像塊中DCT系數(shù)頻帶分布由左上角的直流分量DC往下對(duì)應(yīng)的系數(shù)頻率由低頻升至高頻因此在不影響原圖質(zhì)量的前提下可將水印信息根據(jù)能量大小嵌入相應(yīng)系數(shù)頻帶中。通過(guò)圖像塊量化與水印嵌入結(jié)合的處理方法將水印信息均勻分布在圖像的整個(gè)空間域在圖像裁剪和濾波方面變換域的水印比在空間域的更能表現(xiàn)出一定的魯棒性。2 水印算法描述2.1 水印嵌入算法該算法采用加性嵌入的方式在經(jīng)過(guò)DCT變換后的子圖像塊的中頻域中選取隱秘位置嵌入水印信息具體的嵌入流程如下圖1所示圖1 分塊水印嵌入流程(1分塊處理設(shè)宿主圖像為P將其分塊處理為8*8的K個(gè)子塊。(2水印預(yù)處理設(shè)水印圖像為W對(duì)其進(jìn)行互補(bǔ)變換變換后的水印圖像和變換前的水印圖像相互補(bǔ)。(3對(duì)水印圖像進(jìn)行Arnold置亂變換并依據(jù)混沌映射規(guī)則選取密鑰混沌序列并與水印序列異或運(yùn)算將置換次數(shù)和異或運(yùn)算處理后的結(jié)果分別作為水印嵌入算法的密鑰1和密鑰2。(4)DCT變換對(duì)各子塊內(nèi)做DCT變換利用zig-zag對(duì)DCT系數(shù)進(jìn)行掃描得到第k塊子圖像塊的序列為Zk(i),i0,1,2…63.(5水印嵌入算法依據(jù)zig-zag排序在各子塊的中、低頻段選取特定系數(shù)x(m和x(n在系數(shù)坐標(biāo)a,b和c,d處嵌入水印信息圖像W并將其作為密鑰3。同理嵌入互補(bǔ)水印圖片W’并將嵌入的位置作為密鑰4。水印嵌入的方法如下(6)IDCT變換將每一個(gè)子圖像塊作二維DCT逆變換。(7子塊合并將每一個(gè)子塊合并成嵌入水印的圖像P’。2.2 水印提取算法將嵌入水印的圖像P’分塊處理并對(duì)各子塊進(jìn)行二維DCT變換由密鑰3和4推斷所選擇的水印系數(shù)若x(m≤x(n則水印信息為0若x(m)x(n則水印信息為1再利用密鑰1和2將初步水印的信息解密再進(jìn)行Arnold逆變換最終提取出水印信息。2.3 水印檢測(cè)算法本文通過(guò)計(jì)算峰值信噪比PSNR的值評(píng)價(jià)嵌入水印的宿主圖像的質(zhì)量一幅m和n的圖像PSNR度量標(biāo)準(zhǔn)定義為歸一化相關(guān)系數(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(‘平移后的頻譜’);其運(yùn)行結(jié)果如下圖  圖1 圖像傅里葉變換例 2.2 離散余弦變換DCT 離散余弦變換(DCT)是一組不同頻率和幅值的余弦函數(shù)和來(lái)近似一副圖像實(shí)際上是傅里葉變換的實(shí)數(shù)部分。由于離散余弦變量對(duì)于一副圖像其大部分可視化信息都集中在少數(shù)的變換系數(shù)上。因此離散余弦變量是數(shù)據(jù)壓縮常用的一個(gè)變換編碼方法它能將高相關(guān)數(shù)據(jù)能量集中使得它非常適用于圖像壓縮例如國(guó)際壓縮標(biāo)準(zhǔn)的JPEG格式中就采用了離散余弦變換。 在傅立葉變換過(guò)程中如果被展開(kāi)的函數(shù)是實(shí)偶函數(shù)那么其傅立葉變換中只包含余弦項(xiàng)基于傅立葉變換的這一特點(diǎn)人們提出了離散余弦變換。DCT變換先將圖像函數(shù)變換成偶函數(shù)形式再對(duì)其進(jìn)行二維離散傅立葉變換因此DCT變換可以看成是一種簡(jiǎn)化的傅立葉變換。離散余弦變換葉分為一維離散余弦變換和二維離散余弦變換在圖像處理中用到二維變化所以介紹下其二維定義其他定義可查閱《數(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ù)圖像);其運(yùn)行結(jié)果如下圖2圖2 離散余弦變換例2.3 小波變換DWT 小波變化是對(duì)傅里葉變換和短時(shí)傅里葉變換的一個(gè)突破其改變就在于將無(wú)限長(zhǎng)的三角函數(shù)基換成了有限長(zhǎng)的會(huì)衰減的小波。小波變化的理論較多更多內(nèi)容可查閱《數(shù)字圖像處理》——岡薩雷斯相關(guān)章節(jié)內(nèi)容以下從一個(gè)叫簡(jiǎn)單的角度來(lá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)?三、運(yùn)行結(jié)果?四、matlab版本及參考文獻(xiàn)1 matlab版本2014a2 參考文獻(xiàn)[1]萬(wàn)誼丹.基于Arnold和DCT的抗剪切攻擊圖像水印研究[J].網(wǎng)絡(luò)安全技術(shù)與應(yīng)用. 2021,(08)3 備注簡(jiǎn)介此部分摘自互聯(lián)網(wǎng)僅供參考若侵權(quán)聯(lián)系刪除 仿真咨詢(xún)1 各類(lèi)智能優(yōu)化算法改進(jìn)及應(yīng)用生產(chǎn)調(diào)度、經(jīng)濟(jì)調(diào)度、裝配線調(diào)度、充電優(yōu)化、車(chē)間調(diào)度、發(fā)車(chē)優(yōu)化、水庫(kù)調(diào)度、三維裝箱、物流選址、貨位優(yōu)化、公交排班優(yōu)化、充電樁布局優(yōu)化、車(chē)間布局優(yōu)化、集裝箱船配載優(yōu)化、水泵組合優(yōu)化、解醫(yī)療資源分配優(yōu)化、設(shè)施布局優(yōu)化、可視域基站和無(wú)人機(jī)選址優(yōu)化2 機(jī)器學(xué)習(xí)和深度學(xué)習(xí)方面卷積神經(jīng)網(wǎng)絡(luò)CNN、LSTM、支持向量機(jī)SVM、最小二乘支持向量機(jī)LSSVM、極限學(xué)習(xí)機(jī)ELM、核極限學(xué)習(xí)機(jī)KELM、BP、RBF、寬度學(xué)習(xí)、DBN、RF、RBF、DELM、XGBOOST、TCN實(shí)現(xiàn)風(fēng)電預(yù)測(cè)、光伏預(yù)測(cè)、電池壽命預(yù)測(cè)、輻射源識(shí)別、交通流預(yù)測(cè)、負(fù)荷預(yù)測(cè)、股價(jià)預(yù)測(cè)、PM2.5濃度預(yù)測(cè)、電池健康狀態(tài)預(yù)測(cè)、水體光學(xué)參數(shù)反演、NLOS信號(hào)識(shí)別、地鐵停車(chē)精準(zhǔn)預(yù)測(cè)、變壓器故障診斷3 圖像處理方面圖像識(shí)別、圖像分割、圖像檢測(cè)、圖像隱藏、圖像配準(zhǔn)、圖像拼接、圖像融合、圖像增強(qiáng)、圖像壓縮感知4 路徑規(guī)劃方面旅行商問(wèn)題TSP、車(chē)輛路徑問(wèn)題VRP、MVRP、CVRP、VRPTW等、無(wú)人機(jī)三維路徑規(guī)劃、無(wú)人機(jī)協(xié)同、無(wú)人機(jī)編隊(duì)、機(jī)器人路徑規(guī)劃、柵格地圖路徑規(guī)劃、多式聯(lián)運(yùn)運(yùn)輸問(wèn)題、車(chē)輛協(xié)同無(wú)人機(jī)路徑規(guī)劃、天線線性陣列分布優(yōu)化、車(chē)間布局優(yōu)化5 無(wú)人機(jī)應(yīng)用方面無(wú)人機(jī)路徑規(guī)劃、無(wú)人機(jī)控制、無(wú)人機(jī)編隊(duì)、無(wú)人機(jī)協(xié)同、無(wú)人機(jī)任務(wù)分配6 無(wú)線傳感器定位及布局方面?zhèn)鞲衅鞑渴饍?yōu)化、通信協(xié)議優(yōu)化、路由優(yōu)化、目標(biāo)定位優(yōu)化、Dv-Hop定位優(yōu)化、Leach協(xié)議優(yōu)化、WSN覆蓋優(yōu)化、組播優(yōu)化、RSSI定位優(yōu)化7 信號(hào)處理方面信號(hào)識(shí)別、信號(hào)加密、信號(hào)去噪、信號(hào)增強(qiáng)、雷達(dá)信號(hào)處理、信號(hào)水印嵌入提取、肌電信號(hào)、腦電信號(hào)、信號(hào)配時(shí)優(yōu)化8 電力系統(tǒng)方面微電網(wǎng)優(yōu)化、無(wú)功優(yōu)化、配電網(wǎng)重構(gòu)、儲(chǔ)能配置9 元胞自動(dòng)機(jī)方面交通流 人群疏散 病毒擴(kuò)散 晶體生長(zhǎng)10 雷達(dá)方面卡爾曼濾波跟蹤、航跡關(guān)聯(lián)、航跡融合