« Clustering par K-means, segmentation d'image » : différence entre les versions

De Wiki du LAMA (UMR 5127)
Aller à la navigation Aller à la recherche
Ligne 40 : Ligne 40 :
import numpy as np
import numpy as np
</pre>
</pre>


== Transformation des données en tableaux ==

<pre>
def imageTab(fichier):
"""fonction qui récupère toutes les infos du fichier image, pour les convertir
en tableaux numpy."""
l = fichier.width
h = fichier.height

res = []

for x in range(0,l):
for y in range(0,h):
r,g,b = fichier.getpixel((x,y))

res += [ x, y, r, g, b]

tab = np.array(res, dtype=float)

nTab = tab.reshape(l*h,5)

return nTab


def coordonnees_alea(k, l, h, image):
"""donne k coordonnees et leurs couleurs"""
res = []
for i in range(0,k):

x = randint(0,l-1)

y = randint(0,h-1)

r,g,b = image.getpixel((x,y))
res = res + [[x, y, r, g, b]]

tab = np.array(res, dtype=float)

return tab
</pre>

== Calcul des distances ==

<pre>
def coefficiente_valeur(tab, l, h):
"""Entree : tableau numpy
Sortie : tableau numpy avec les coordonnées et les couleurs entre 0 et 1"""

tab2 = np.zeros((len(tab),5), dtype=float)
tab2[:,0] = tab[:,0] / l
tab2[:,1] = tab[:,1] / h
tab2[:,2] = tab[:,2] / 255
tab2[:,3] = tab[:,3] / 255
tab2[:,4] = tab[:,4] / 255
return tab2


def distancePoints(tab1, tab2, k, coef):
"""donne la distance en couleur entre 2 points
Entrée : deux tableaux numpy
Sortie : un tableau numpy avec les distances entre les points de l'image et les k"""

taille = len(tab1)

nTab = np.zeros((taille,k),dtype=float)

for i in range(0,taille):
for j in range(0,k):
nTab[i][j] = coef*( tab1[i][0] - tab2[j][0] )**2 + coef*( tab1[i][1] - tab2[j][1] )**2 + ( tab1[i][2] - tab2[j][2] )**2 + ( tab1[i][3] - tab2[j][3] )**2 + ( tab1[i][4] - tab2[j][4] )**2
return nTab


def distance_plus_courte(tab):
"""Entrée : tableau numpy avec les k distances pour chaque points
Sortie : tableau numpy des indices où la distance est la plus petite"""

indice = np.argmin(tab,axis=1)
return indice
</pre>

== Attribution aux clusters ==

<pre>
def kTab(k):
"""Entrée : un entier k
Sortie: un tableau de k tableaux vide"""
tab = []
for i in range(0,k):
tab = tab + [[]]
return tab


def attribution_aux_clusters(tabImage, tabIndice, k):
"""Entrée : 2 tableaux Numpy, un avec les info de l'image, et les indices des distances les plus courtes
Sortie : tableau de tableaux numpy des clusters"""

cluster = kTab(k)
for i in range(0,len(tabImage)):

indice = tabIndice[i]

cluster[indice] += [tabImage[i].tolist()]

return cluster


def clustersEnNumpy(tab,k):
"""Entrée : un tableau de tableaux de tableaux de points, un nb de cluster
Sortie : un tableau de tableaux Numpy"""

Ntab = kTab(k)

for i in range(0,k):
Ntab[i] = np.array(tab[i],dtype=float)

return Ntab
</pre>

== Nouveaux Centroïdes ==

<pre>
def NouveauxCentres(tab,k):
"""Entrée : tableau de tableaux numpy avec les coordonnées des points appartenant à un cluster
Sortie : un tableau Numpy avec les du nouveaux centres"""

Nclusters = np.zeros((k,5),dtype=float)

for i in range(0,k):

moy = ( np.sum(tab[i],axis=0) / len(tab[i]) )
moy = moy.astype(int)

Nclusters[i] = moy

return Nclusters


def calculCentres(tabImage, tabImageCoeff, tabCentre, k , l, h, coeff):
"""fonction qui reprend toutes celles précédentes, pour calculer un nouveau centre"""

coeffCentre = coefficiente_valeur(tabCentre, l , h)
distance = distancePoints(tabImageCoeff, coeffCentre,k,coeff)
indice = distance_plus_courte(distance)
clusters = attribution_aux_clusters(tabImage,indice,k)
NumpyCluster = clustersEnNumpy(clusters,k)

res = NouveauxCentres(NumpyCluster,k)

return res


def meilleurs_centres(tabImage, tabImageCoeff,k,l,h,coeff,image):
"""Entrée : 2 tableaux Numpy avec les infos de l'image (1 coefficienté et l'autre non), 3 entiers, 1 flottant et une procédure d'ouverture d'image)
Sortie : un tableau Numpy, le centre optimal du cluster"""

centre = coordonnees_alea(k, l, h, image)
nouveauCentre = calculCentres(tabImage, tabImageCoeff, centre,k,l,h,coeff)
i = 0
print("Execution en cours : 0 %")
while i < 10:
centre = nouveauCentre
nouveauCentre = calculCentres(tabImage, tabImageCoeff, centre,k,l,h,coeff)
print("Execution en cours : "+str(i+1)+"0 %")
i = i +1

return nouveauCentre
</pre>

== Modification des couleurs de l'image ==

<pre>
def change_couleur_cluster(tabCluster, tabCentre, image):
"""Entrée : tableau numpy de points d'un cluster
Sortie : zone de couleur changée"""

numpyCluster = tabCluster.astype(int)
NtabCluster = numpyCluster.tolist()

numpyCentres = tabCentre.astype(int)
NtabCentres = numpyCentres.tolist()
for i in range(0,len(NtabCluster)):

image.putpixel( (NtabCluster[i][0] , NtabCluster[i][1] ), ( NtabCentres[2], NtabCentres[3], NtabCentres[4] ) )
</pre>


== Fonction finale, le clustering ==


<pre>
<pre>

Version du 11 mai 2021 à 19:27

Etudiant : Paul AUBRY

Tuteur : Jacques-Olivier LACHAUD


Introduction : Clustering par k-means

Le clustering ou algorithme des k moyennes a pour but de regrouper des populations en communautés disposant de critères communs proches, jusqu'à avoir des communautés homogènes qu'on appellera cluster et qui ont pour représentant un centroïde.
Pour déterminer ces clusters, on regroupera les différents éléments en fonction d'une distance. Cette notion de distance est différente en fonction des domaines d'applications.
Le clustering par k-means peut être utilisé pour faire de la segmentation de clientèle, du clustering en Data Mining ou encore sur des images.

Algorithme pour l'image

Nous allons voir ici, comment procéder de manière théorique, pour réaliser un clustering. Tout d'abord, nous devons choisir k points aléatoirement, qui seront les centroïdes, et un coefficient λ qui nous permettra de calculer les distances.
Ensuite, nous allons affecter chaque point de l'image à un cluster. Pour cela, il faut calculer la distance entre le point, et chaque k. La distance la plus courte nous permettra de déterminer quel cluster choisir.
Un point est représenté de la manière suivante : [ x , y , r , g , b ] .
Avec x et y les coordonnées et r, g, b les composantes de couleurs.

Pour calculer la distance on utilise la formule suivante :


Lorsque les clusters sont définis, on calcul la moyenne de chaque cluster, afin de récupérer de nouveaux centres.
On réitère les actions vus précédemment mais cette fois-ci avec les nouveaux centres.
Dès que l'on obtient des centres "stable", on peut modifier l'image.

Réalisation grâce à Python

Pour ce faire, il faut installer plusieurs bibliothèques.

  • "numpy", qui va nous servir à effectuer les calculs de manières bien plus rapide.
  • "PIL" pour le traitement des images.
from random import *
from PIL import Image
import numpy as np
def clustering(k,coef):
    """Entrée :1 entier un flottant
    Sortie : image modifiée"""
    
    image0 = Image.open("Kowloon-small-329x216.png")
    image = image0

    l = image.width
    h = image.height

    tabImage = imageTab(image)
    tabImageCoeff = coefficiente_valeur(tabImage, l, h)

    centroides = meilleurs_centres(tabImage, tabImageCoeff, k, l ,h, coef, image)
    coeffCentroides = coefficiente_valeur(centroides, l, h)

    distance = distancePoints(tabImageCoeff, coeffCentroides,k,coef)
    indice = distance_plus_courte(distance)

    clusters = attribution_aux_clusters(tabImage,indice,k)
    NumpyClusters = clustersEnNumpy(clusters,k)

    print("En cours de finalisation...")

    for i in range(0,k):
        change_couleur_cluster( NumpyClusters[i] , centroides[i] , image )
        

    image.save(str(k)+"-"+str(coef)+"Kowloon-small-329x216.png")    
    image.show()
    image0.close()
    image.close()

Résultats