1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39
| %%%%%%%%%%%%%%%%%
clc
clear all
[filename, pathname] = uigetfile('*.jpg;*.tif;*.gif;*.png;*.bmp','MultiSelect', 'On');
jj=length(filename);
a=0;
for k=1:jj
[c,map]=imread(filename{k});
[m,n]=size(c);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
ws = 9; % windws size
[nrow,ncol] = size(c);
s = (ws+1)/2;
imag = padarray(c, [s s], 'symmetric', 'both'); % Pad image array
imag(s+nrow,:) = []; imag(s,:) = [];
imag(:,s+ncol) = []; imag(:,s) = [];
ss = (ws-1)/2;
for i = ss+1:nrow+ss
for j = ss+1:ncol+ss
win = imag(i-ss:i+ss,j-ss:j+ss);
meanWin(i,j) = mean2(win);
stdWin(i,j) = std2(win);
end
end
meanWin = meanWin(ss+1:nrow+ss,ss+1:ncol+ss);
stdWin = stdWin(ss+1:nrow+ss,ss+1:ncol+ss);
imwrite(meanWin(k).CData,sprintf('im%d.gif',k),'gif')
end |
Partager