【正文】 ? clear ? close all ? I = imread(39。39。)。 ? tform = maketform(39。affine39。,[1 0 0。 .5 1 0。 0 0 1])。 ? J = imtransform(I,tform)。 ? imshow(I), figure, imshow(J) 圖像的二維空間變換 ? 透視投影變換示例： ? clc ? clear ? close all ? I = imread(39。39。)。 ? udata = [0 1]。 vdata = [0 1]。 % input coordinate system ? tform = maketform(39。projective39。,[ 0 0。 1 0。 1 1。 0 1],... ? [4 2。 8 3。 3 5。 6 3])。 ? [B,xdata,ydata] = imtransform(I, tform, 39。bicubic39。, ... ? 39。udata39。, udata,... ? 39。vdata39。, vdata,... ? 39。size39。, size(I),... ? 39。fill39。, 128)。 ? subplot(1,2,1), imshow(udata,vdata,I), axis on ? subplot(1,2,2), imshow(xdata,ydata,B), axis on 圖像的鄰域和塊操作 ? 塊操作 ? 鄰域操作 塊操作 ? blkproc ( ) ? 實現(xiàn)對圖像的分塊處理 ? B = blkproc(A,[m n],fun) ? 將 A分成 m*n的互不重疊矩形塊，若 A不能被 m和 n整除，則補充 0 ? 然后對每個塊按函數(shù) fun進行處理 ? 示例： ? I = imread(39。39。)。 ? fun = @dct2。 ? J = blkproc(I,[8 8],fun)。 ? imshow(J) 鄰域操作 ? nlfilter ( ) ? 實現(xiàn)對圖像的滑塊鄰域操作 ? B = nlfilter(A, [m n], fun) ? 一次運算只處理圖像的一個像素 ? 輸出圖像的像素值由其對應輸入圖像像素點的鄰域內的像素采用有關運算來決定 ? 當操作從一個元素移動到另一個元素時，矩陣塊也以相同的方向移動 鄰域操作 ? clc ? clear ? close all ? A = imread(39。39。)。 ? fun = @median。 ? B = nlfilter(A,[3 3],fun)。 ? imshow(A), figure, imshow(B) 提高 matlab的運行效率 ? 數(shù)組預分配 ? 循環(huán)向量化 ? 簡單函數(shù)內斂話 數(shù)組預分配 ? x=2。 ? x(3,6)=6 ? x = ? 0 0 0 0 0 0 ? 0 0 0 0 0 0 ? 0 0 0 0 0 6 ? 自動調整矩陣的能力，會極大的增加系統(tǒng)開銷，影響程序運行效率 數(shù)組預分配 ? 例如： ? clc ? clear ? tic ? a(1)=1。 ? b(1)=0。 ? for k=2:8000 ? a(k)=*a(k1)*b(k1)。 ? b(k)=*a(k1)+*b(k1)。 ? end ? toc ? t=toc ? 輸出： ? t = ? 數(shù)組預分配 ? 例如： ? clc； clear ? tic ? a=zeros(1,80000)。 ? b=zeros(1,80000)。 ? a(1)=1。 b(1)=0。 ? for k=2:8000 ? a(k)=*a(k1)*b(k1)。 ? b(k)=*a(k1)+*b(k1)。 ? end ? toc； ? t=toc ? 輸出： ? t = ? 循環(huán)向量化 ? MATLAB是矩陣語言為矢量和矩陣操作設計的。要盡可能把要進行的計算編寫成用矢量和矩陣運算的形式 ? clc； clear ? tic 簡單函數(shù)內斂化 ? MATLAB在調用任何一個函數(shù)時，除了函數(shù)本身執(zhí)行的開銷外，還要為函數(shù)調用增加額外開銷 ? 如果一段循環(huán)執(zhí)行的代碼中頻繁的調用函數(shù)，這種函數(shù)調用對程序的運行效率影響非常明顯 ? 建議對功能 非常簡單的函數(shù) ，如果在循環(huán)執(zhí)行代碼中需要反復調用它們的話， 可以將這些簡單函數(shù)的實現(xiàn)直接設計到自己代碼中 ，這樣可以提高運行效率 簡單函數(shù)內斂化 ? 例如： ? function a=ff(n) ? if (n==1) ? a=1。 ? elseif (n==2) ? a=0。 ? elseif (n==3) ? a=1。 ? else ? a=ff(n1)2*ff(n2)+ff(n3)。 ? end ? 主函數(shù) main： ? clc ? Clear ? tic ? for i=1:20 ? fac=zeros(1,20) ? fac(i)=ff(i) ? end ? toc ? t=toc ? 輸出： ? t = ? 簡單函數(shù)內斂化 ? 例如： ? clc ? clear ? tic ? f(1)=1。 ? f(2)=0。 ? f(3)=1。 ? for n=4:20 ? f(n)=f(n1)2*f(n2)+f(n3)。 ? end ? toc ? t2=toc ? 輸出： ? t 2= ? 小結 ? 掌握 matlab圖像的數(shù)據(jù)類型及轉換函數(shù) ? 掌握圖像文件的讀取和顯示方法 ? 掌握圖像的運算 ? 熟悉提高圖像運行效率的方法