Ecualización de una imagen Algoritmo de Mario D G

clear all
clc
A=imread('c:\moon.jpg');
%A=rgb2gray(A); %Si se necesita Ecualizar la imagen
B=double(A);
B=histeq(A);
minimo=min(min(B));
maximo=max(max(B));
dinamica=maximo-minimo;
cambio=255/dinamica;
posicion=minimo-1;
[fila,columna]=size(B);
Z=zeros(fila,columna);
for y=1:fila
    for x=1:columna
          Z(y,x)=((B(y,x)-posicion)*cambio);
    end
end
Z=uint8(Z);
imshow(Z)
figure
imhist(Z)
figure
imshow(A)
figure
imhist(A)
minimoZ=min(min(Z));
maximoZ=max(max(Z));
dinamicaZ=maximoZ-minimoZ;