Istogramma equalizzazione per Python

Question

devo questa routine di fare equalizzazione di una foto:

def histeq(im,nbr_bins=256): 

    #get image histogram 
    imhist,bins = histogram(im.flatten(),nbr_bins,normed=True) 
    cdf = imhist.cumsum() #cumulative distribution function 
    cdf = 255 * cdf/cdf[-1] #normalize 

    #use linear interpolation of cdf to find new pixel values 
    im2 = interp(im.flatten(),bins[:-1],cdf) 

    return im2.reshape(im.shape), cdf 

#im = array(Image.open('AquaTermi_lowcontrast.jpg').convert('L')) 
im = array(Image.open('Unequalized.jpg').convert('L')) 
#Image.open('plant4.jpg').convert('L').save('inverted.jpg') 

im2,cdf = histeq(im) 

plt.imshow(im2) 
plt.savefig("outputhisto.jpg") 

Quando eseguo questo con il picture dalla pagina wiki di histogram equalization si traduce in questo:

Invece di regolare correttamente il contrasto dell'immagine su qualcosa lungo le linee di this. Che cosa sto facendo di sbagliato?

Answer 1

Sei sicuro di non eseguire il rendering con la mappa colori errata? Prova

plt.imshow(im2, cmap=plt.cm.gray) 

O

plt.imshow(im2, cmap=plt.get_cmap('gray'))