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 40 41 42 43
| clear all;
clc;
Ima1 = imread('C:\image1.tif');
FenSize = [500 500];
ImaSize= size(Ima1);
Class_Vect = unique(Ima1); % toutes les classes de l'image
Ima2 = imread('C:\image2.tif');
for i = 1:FenSize(1):size(Ima1,1)
for j = 1:FenSize(2):size(Ima1,2)
Fen = Ima1(i:i+FenSize(1)-1 , j:j+FenSize(2)-1);% Prélèvement du bloc
Fen1 = Ima2 (i:i+FenSize(1)-1 , j:j+FenSize(2)-1);% Prélèvement du bloc
MoyIm = mean(Fen1(:));
Fen1 = double(Fen1(:));
VarIm = var(Fen1(:));
Nbr = histc(Fen(:),Class_Vect); % Classes du bloc Fen
Nbr(Nbr==0)=[];
result = [Class_Vect',Nbr', MoyIm, VarIm];
folder = 'c:\SARAH\'
newimagename = [folder 'image_small_' num2str(i) '.bmp'];
imwrite(Fen,newimagename)
end
if exist('D:\images\results.txt','file')~=7
fid=fopen('D:\images\results.txt','a+');
if fid ==-1;
errordlg(['File: can not be recreated']);
return;
else
fprintf(fid,'%f %f %f %f \t\r\n',result');
fclose(fid);
end
end
end |
Partager