【正文】 dding in a Medical Image[J]. Image Processing,2004,4(1):2315~2318.[28] , lossless data embedding scheme for medical images in application of Annual the IEEE EMBS[A],2005,Page(s):852~855.[29] . The reconstruction of a multidimensional sequence from the phase or magnitude of the FFT[A].IEEE Acoustics,Speech and Signal Processing,1992,(4):140~154[30] , , M, recovery of an original image containing embedded data[A]. US Patent, Docket No:7102~7105,19990923.附 錄（1）水印嵌入程序：clear all。 clc。start_time=cputime。I=imread(39。39。)。I=rgb2gray(I)。I=double(I)/255。 I=round(I)。figure(1)。subplot(1,2,1)。imshow(I),title(39。水印圖像39。)dimI=size(I)。rm=dimI(1)。cm=dimI(2)。mark=I。alpha=10。k1=randn(1,8)。k2=randn(1,8)。a0=imread(39。39。)。psnr_cover=double(a0)。a0=rgb2gray(a0)。[r,c]=size(a0)。figure(2)。subplot(1,2,1)。imshow(39。39。)。title(39。原始圖像39。)。subplot(1,2,2)。imshow(a0,[]),title(39?；叶葓D像39。)。cda0=blkproc(a0,[8,8],39。dct239。)。 %分塊處理圖片cda1=cda0。 % cda1 = 256_256for i=1:rm % i=1:32 for j=1:cm % j=1:32 x=(i1)*8。y=(j1)*8。 if mark(i,j)==1 k=k1。 else k=k2。 end cda1(x+1,y+8)=cda0(x+1,y+8)+alpha*k(1)。 cda1(x+2,y+7)=cda0(x+2,y+7)+alpha*k(2)。 cda1(x+3,y+6)=cda0(x+3,y+6)+alpha*k(3)。 cda1(x+4,y+5)=cda0(x+4,y+5)+alpha*k(4)。 cda1(x+5,y+4)=cda0(x+5,y+4)+alpha*k(5)。 cda1(x+6,y+3)=cda0(x+6,y+3)+alpha*k(6)。 cda1(x+7,y+2)=cda0(x+7,y+2)+alpha*k(7)。 cda1(x+8,y+1)=cda0(x+8,y+1)+alpha*k(8)。 endenda1=blkproc(cda1,[8,8],39。idct239。)。 a_1=uint8(a1)。imwrite(a_1,39。39。,39。jpg39。)。figure(3)。subplot(1,2,1)。imshow(a1,[]),title(39。嵌入水印后的圖像39。)。（2）水印提取程序：M1=imread(39。39。)。 psnr_watermarked=M1。dca1=blkproc(M1,[8,8],39。dct239。)。 p=zeros(1,8)。for i=1:rm % i=1:32 for j=1:cm % j=1:32 x=(i1)*8。y=(j1)*8。 p(1)=dca1(x+1,y+8)。 p(2)=dca1(x+2,y+7)。 p(3)=dca1(x+3,y+6)。 p(4)=dca1(x+4,y+5)。 p(5)=dca1(x+5,y+4)。 p(6)=dca1(x+6,y+3)。 p(7)=dca1(x+7,y+2)。 p(8)=dca1(x+8,y+1)。 %sd1=sum(sum(p.*k1))/sqrt(sum(sum(p.^2)))。 %sd2=sum(sum(p.*k2))/sqrt(sum(sum(p.^2)))。 %if sd1sd2 if corr2(p,k1)corr2(p,k2),warning off MATLAB:divideByZero。 mark1(i,j)=1。 else mark1(i,j)=0。 end endendfigure(5)。 subplot(1,2,1)。imwrite(mark1,39。39。,39。jpg39。)。imshow(mark1),title(39。提取的水印圖像39。)。I=imread(39。39。)。 I=rgb2gray(I)。I=double(I)/255。 I=round(I)。subplot(1,2,2)。imshow(I),title(39。原始水印圖像39。)。（3）噪聲攻擊程序：a1=imread(39。39。)。 Wimage2=a1。noise0=2*randn(size(Wimage2))。Wimage2=double(Wimage2)+noise0。 figure(4)。subplot(1,2,1)。imshow(uint8(Wimage2))。title(39。加入白噪聲后圖像39。)。M1=Wimage2。M_1=uint8(M1)。imwrite(M_1,39。39。,39。jpg39。)。 （4）高斯低通濾波攻擊： a1=imread(39。39。)。 WImage3=a1。H=fspecial(39。gaussian39。,[10,10],2)。WImage3=imfilter(WImage3,H)。 figure(4)。subplot(1,2,1)。imshow(WImage3,[])。title(39。高斯低通濾波后圖像39。)。M1=WImage3。M_1=uint8(M1)。imwrite(M_1,39。39。,39。jpg39。)。 （5）剪切攻擊程序：a1=imread(39。39。)。 Wimage2=a1。Wimage2(1:64,1:128,)=256。figure(4)。subplot(1,2,1),imshow(Wimage2,[]),title(39。剪切后的圖象39。)。M1=Wimage2。M_1=uint8(M1)。imwrite(M_1,39。39。,39。jpg39。)。 （6）旋轉(zhuǎn)攻擊程序：a1=imread(39。39。)。 Wimage2=a1。Wimage2=imrotate(Wimage2,90,39。bilinear39。,39。crop39。)。figure(4)。subplot(1,2,1)。imshow(uint8(Wimage2))。title(39。加入白噪聲后圖像39。)。M1=Wimage2。M_1=uint8(M1)。imwrite(M_1,39。39。,39。jpg39。)。 （7）求解峰值信噪比和歸一化系數(shù)程序：x1=a0。x2=imread(39。39。)。x1=double(x1)。x2=double(x2)。[m,n]=size(x1)。[m2,n2]=size(x2)。if m2~=m||n2~=n。error(39。圖像選擇錯(cuò)誤39。)。end。msevalue=0。for i=1:m。for j=1:n。msevalue=msevalue+(x1(i,j)x2(i,j))^2。end。end。msevalue=msevalue/(m*n)。disp([39。輸入數(shù)據(jù)的mse為:39。,num2str(msevalue)])。psnrvalue=255^2/msevalue。psnrvalue=10*log10(psnrvalue)。disp([39。輸入數(shù)據(jù)的psnr為:39。,num2str(psnrvalue)])。（8）計(jì)算NC（歸一化相關(guān)系數(shù)）%tmark_1=tmark_1(1:mo,1:no)。%%tmark_1g_mark=double(mark1)。o_mark=double(mark)。[p,q]=size(g_mark)。nc_0=0。nc_1=0。nc_2=0。for i=1:p for j=1:q nc_0=nc_0+g_mark(i,j)*o_mark(i,j)。 nc_1=nc_1+o_mark(i,j)*o_mark(i,j)。 nc_2=nc_2+g_mark(i,j)*g_mark(i,j)。 endendNC=nc_0/sqrt(nc_1*nc_2)。disp([39。輸入數(shù)據(jù)的NC為:39。,num2str(NC)])。畢業(yè)設(shè)計(jì)（論文）原創(chuàng)性聲明和使用授權(quán)說明原創(chuàng)性聲明本人鄭重承諾：所呈交的畢業(yè)設(shè)計(jì)（論文），是我個(gè)人在指導(dǎo)教師的指導(dǎo)下進(jìn)行的研究工作及取得的成果。盡我所知，除文中特別加以標(biāo)注和致謝的地方外，不包含其他人或組織已經(jīng)發(fā)表或公布過的研究成果，也不包含我為獲得 及其它教育機(jī)構(gòu)的學(xué)位或?qū)W歷而使用過的材料。對本研究提供過幫助和做出過貢獻(xiàn)的個(gè)人或集體，均已在文中作了明確的說明并表示了謝意。作 者 簽 名：　　 　 　　日　 期：　　 　 　　 指導(dǎo)教師簽名：　　 　 　　 日　　期：　　 　 　　 使用授權(quán)說明本人完全了解 大學(xué)關(guān)于收集、保存、使用畢業(yè)設(shè)計(jì)（論文）的規(guī)定，即：按照學(xué)校要求提交畢業(yè)設(shè)計(jì)（論文）的印刷本和電子版本；學(xué)校有權(quán)保存畢業(yè)設(shè)計(jì)（論文）的印刷本和電子版，并提供目錄檢索與閱覽服務(wù)；學(xué)?？梢圆捎糜坝?、縮印、數(shù)字化或其它復(fù)制手段保存論文；在不以贏利為目的前提下，學(xué)校可以公布論文的部分或全部內(nèi)容。作者簽名：　　 　 　　 日　 期：　　 　 　　 學(xué)位論文原創(chuàng)性聲明本人鄭重聲明：所呈交的論文是本人在導(dǎo)師的指導(dǎo)下獨(dú)立進(jìn)行研究所取得的研究成果。除了文中特別加以標(biāo)注引用的內(nèi)容外，本論文不包含任何其他個(gè)人或集體已經(jīng)發(fā)表或撰寫的成果作品。對本文的研究做出重要貢獻(xiàn)的個(gè)人和集體，均已在文中以明確方式標(biāo)明。本人完全意識到本聲明的法律后果由本人承擔(dān)。作者簽名： 日期： 年 月 日學(xué)位論文版權(quán)使用授權(quán)書本學(xué)位論文作者完全了解學(xué)校有關(guān)保留、使用學(xué)位論文的規(guī)定，同意學(xué)校保留并向國家有關(guān)部門或機(jī)構(gòu)送交論文的復(fù)印件和電子版，允許論文被查閱和借閱。本人授權(quán)　　 　 　大學(xué)可以將本學(xué)位論文的全部或部分內(nèi)容編入有關(guān)數(shù)據(jù)庫進(jìn)行檢索，可以采用影印、縮印或掃描等復(fù)制手段保存和匯編本學(xué)位論文。涉密論文按學(xué)校規(guī)定處理。作者簽名： 日期： 年 月 日導(dǎo)師簽名： 日期： 年 月 日注 意 事 項(xiàng)（論文）的內(nèi)容包括：1）封面（按教務(wù)處制定的標(biāo)準(zhǔn)封面格式制作）2）原創(chuàng)性聲明3）中文摘要（300字左右）、關(guān)鍵詞4）外文摘要、關(guān)鍵詞 5）目次頁（附件不統(tǒng)一編入）6）論文主體部分：引言（或緒論）、正文、結(jié)論7）參考文獻(xiàn)8）致謝9）附錄（對論文支持必要時(shí)）：理工類設(shè)計(jì)（論文）正文字?jǐn)?shù)不少于1萬字（不包括圖紙、程序清單等）。：任務(wù)書、開題報(bào)告、外文譯文、譯文原文（復(fù)印件）。、圖表要求：1）文字通順，語言流暢，書寫字跡工整，打印字體及大小符合要求，無錯(cuò)別字，不準(zhǔn)請他人代寫2）工程設(shè)計(jì)類題目的圖紙，要求部分用尺規(guī)繪制，部分用計(jì)算機(jī)繪制，所有圖紙應(yīng)符合國家技術(shù)標(biāo)準(zhǔn)規(guī)范。圖表整潔，布局合理，文字注釋必須使用工程字書寫，不準(zhǔn)用徒手畫3）畢業(yè)論文須用A4單面打印，論文50頁以上的雙面打印4）圖表應(yīng)繪制于無格子的頁面上5）軟件工程類課題應(yīng)有程序清單，并提供電子文檔1）設(shè)計(jì)（論文）2）附件：按照任務(wù)書、開題報(bào)告、外文譯文、譯文原文（復(fù)印件）次序裝訂3）其它 41