Discussion :

[Nana11]

Bonjour, j'ai un projet à réaliser et je dois créer un programme pour traiter des images pour relever les valeurs d'intensité pour chaque pixels.

Mais j'obtiens des résultats pas très pratiques à manier j'aimerais que les intensités soient affichées en colonnes et à côté le numéro d'image pour ensuite réaliser des graphiques.

Quelqu'un pourrait-il m'aider ?

Voici mon programme :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
import Image
im = Image.open("nov0311a1_8_8bits.tif")
F0= 158
for i in range(0,60):
    im.seek(i)
    region = im.crop((52,20,62,28))
    cc = region.getcolors()
    print (cc,i)

[josmiley]

ho ! mon code

comme ça ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
import Image
im = Image.open("nov0311b1_10_8bits.tif")
F0= 158
for i in range(0,60):
    im.seek(i)
    region = im.crop((52,20,62,28))
    for n,cc in region.getcolors():
        print "image num :",i,"nombre de pixel d'intensite",cc,":",n

[Nana11]

Ah super merci !
Par contre quand je veux enregistrer les résultats je fais :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
fout=open( "test.csv", "w")
for i in range (0,60):
    output = str(cc)+"\t"+str(i)+"\n"
    fout.write (output)
fout.close()

ça m'affiche : " fout=open( "test.csv", "w")
IOError: [Errno 13] Permission denied: 'test.csv' "

[josmiley]

tu n'as pas le droit d'écrire dans le répertoire courant ...

[Nana11]

Ah mince !
Pourtant avant sur le code de départ j'arrivais à enregistrer les résultats avec ce code là. Comment puis-je faire ?

[utopman]

Salut,
moi je te propose ça pour tes images

il faut Pil et matplotlib

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
 
import Image
import matplotlib.pyplot as plt
import sys
 
filename = sys.argv[1]
 
im = Image.open(filename)
data = list(im.getdata())
r,g,b = im.split()
 
 
red 	= r.getdata() 
green 	= g.getdata() 
blue 	= b.getdata() 
 
 
if '-s' in sys.argv:
	size = range(256)
	plt.title('intensity')
	plt.plot(size,r.histogram(),'r')
	plt.plot(size,g.histogram(),'g')
	plt.plot(size,b.histogram(),'b')
	plt.show()
 
if '-o' in sys.argv:
	o = open(filename + '.csv','w')
	o.write('r\tg\tb\n')
	hr = r.histogram()
	hg = g.histogram()
	hb = b.histogram()
 
	for d in xrange(255):
		o.write('%s\t%s\t%s\n' % (hr[d],hg[d],hb[d]))
	o.close()

tu peux afficher les resultats avec l'option -s et les sortir en csv avec -s
exemple : python monimage.jpg -o -s
affiche et sort un fichier csv avec les données de l'histogramme

je n'ai testé qu'avec des jpg...
cordialement

[Nana11]

On ne doit utiliser que la bibliothèque PIL
Merci quand même pour ta réponse =)

[PauseKawa]

tu n'as pas le droit d'écrire dans le répertoire courant ...

Ou le fichier existe en lecture seule pour vous.
Un test simple pour cela:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

fout=open( "test1.csv", "w")

Ou avec le module os:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
if os.access("test.csv", os.W_OK):
    fout=open( "test.csv", "w")
    for i in range (0,60):
        output = str(cc)+"\t"+str(i)+"\n"
        fout.write (output)
    fout.close()

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
if os.access("test.csv", os.W_OK):
    with open( "test.csv", "w") as fout:
        for i in range (0, 60):
            fout.write(str(cc)+"\t"+str(i)+"\n")

Il est toujours possible de passer par try plutôt qu'os:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
try:
    fout = open("test.csv", "w")
except IOError:
    print('Interdit')
else:
    with fout:
        for i in range (0, 180):
            fout.write(str(cc)+"\t"+str(i)+"\n")

@+

[utopman]

si tu avais pris la peine de lire un peu mieux ma réponse et donc mon message, tu aurais compris qu'il y a deux possibles. tu peux enlever l'utilisation de matplotlib et simplement sortir les valeurs en csv de tes intensités en omettant completement l'import de matplotlib et de la condition sur l'argument -s.

[Nana11]

Ah désolé utopman je pensais qu'il fallait utiliser les deux ( je ne suis tellement pas douée en Python que j'ai même pas su le reconnaître qu'on pouvait faire qu'avec un ^^)
Mais mon image est en noir et blanc du coup je ne peux pas la découper par bandes r,v,b si ?


PauseKawa : quand j'utilise ce que vous me proposez je n'ai plus de message d'erreur mais j'obtiens un fichier avec deux colonnes une qui varie entre 1 et 60 et l'autre qui reste bloquée à 158 alors que mes résultats sont complètement différents, comment ça se fait ?

[PauseKawa]

Bonjour,

L'erreur ne viens pas de mon code (qui ne concerne que l'écriture du fichier) mais du votre:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
for i in range (0,60):
    output = str(cc)+"\t"+str(i)+"\n"
    fout.write (output)

Ici seul i change, cc reste à la dernière valeur que vous lui avez donner (158).

@+

[Nana11]

Ah oui c'est bon ça marche merci pour votre aide =)