function [Z,Imcluster] = isodata(image, k, L, I, ON, OC, OS, NO, min);
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  Parametros interos de la funcion  %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
data = double(imread(image));
% data=rgb2gray(data);
[a,b] = size(data);
ImVecteur = reshape(data,[a*b 1]);
r = size(ImVecteur);
cluster=zeros(r); % Espacio temporal interno de la funcion.


cluster=zeros(r); % Espacio temporal interno de la funcion.
iter=0;
final=0;
vuelve3=0;
A=1;     % Nous l'avons initialis