Discussion :

[Bionics]

Bonjour,

Après plusieurs recherches je n'ai pas trouvé ce que je voulais, seulement des pistes mais je n'arrive pas à faire le code que je veux. Voilà le problème :

J'ai un dossier contenant 8100 images, je veux toutes les traiter mais par groupe de 180.
Ce que je pensais faire c'est :

Une boucle qui sélectionne les 180 premières images, les lise, fasse des calculs dans cette boucle (moyenne etc).
Jusque là je n'ai pas de problèmes mais là où je bute c'est pour passer aux 180 images suivantes..
Il faudrait une boucle englobant le tout, qui partage d'abord les 8100 images en groupes de 180 images? Et chaque groupe serait indexé? Et donc je traiterai groupe par groupe?

Je suis perdu mais j'espère que vous avez compris là où je veux en venir, merci d'avance.

[Jean Dumoncel]

Bonjour,

il faudrait nous montrer la structure de ton code pour que l'on puisse t'aider efficacement. Tu devrais pouvoir faire une boucle du genre :

Code :

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

1
2
3
for k=1:180:8100
… traitement par groupe de 180
end

[Bionics]

Je ne comprends pas ce que fait le 8100 ici? Je suis pas très doué en boucle, pourtant c'est la base

Code :

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

for i=1:180:8100;

Pour en venir à mon code, prenons seulement 180 images pour commencer.
A l'heure actuelle la boucle est sensée lire une image appelée nuit02a03-i.fits avec i variant de 1 à 180. Avant de recommencer un tour de boucle, elle stock l'image (matrice) indexée i dans une grande matrice appelée matricefinale. Et a chaque tour la nouvelle matrice img contenant l'image nuit02a03-i.fits est stockée dans la matricefinale à la suite de la précédente. Je ne sais pas si cela fonctionne car je ne peux pas l'essayer.. Voilà le code :

Code :

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

1
2
3
4
5
6
7
8
9
10
repertoire = 'C:\Users\blablabla..';
 
for i=1:180;
 
    filenam = sprintf('nuit02a03-%01d.fits',i);
    chemin = fullfile(repertoire, filenam);
 
    img = fitsread(chemin);
    matricefinale = [ matricefinale img ]
end;

Ainsi, j’obtiens une grande matrice (ou plutôt un tenseur de rang 3) contenant mes 180 images empilées suivant Z. Je peux donc faire la moyenne de toutes les valeurs de ce tenseur et enregistrer cette moyenne dans un vecteur. Mais ici c'est seulement pour 180 images.. Et j'ai vraiment du mal à comprendre comment réitérer cette opération mais en partant d'un i = 181;360 et ce jusqu'à 8100..

Merci d'avance

[Gooby]

Bonjour

Un truc dans ce genre:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
repertoire = 'C:\Users\blablabla..';
 
infos = dir(repertoire);
 
for q=1:180:length(infos.name)-2-180
     for i=q:q+180;
         chemin = fullfile(repertoire, infos.name(i+2);
          img = fitsread(chemin);
          matricefinale{q,1} = cat(3,matricefinale,img);
     end
end

[Bionics]

Merci, ton code semble être exactement ce qu'il me faut! Cependant je ne comprends le infos.name(i+2). (Pour le code suivant j'ai mis i au lieu de i+2 du coup, est-ce une erreur de ma part?)

Admettons que je n'ai plus besoin de ces sous-séquences, que je traite les 8100 images d'un coup.
Dans une boucle for je fais la lecture de l'image, ensuite je modifie l'image et je calcule sa moyenne (moyenne de la valeur des pixels). Je veux ensuite entrer cette valeur dans une matrice à la position 1x1, et au second tour de boucle mettre la seconde moyenne à la position 1x2 etc etc (toujours une ligne). Est-ce que la dernière ligne me permet de faire ça?

Code :

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

1
2
3
4
5
6
7
8
9
10
repertoire = 'C:\Users\blablabla..';
infos = dir(repertoire);
 
for i=1:num;
    chemin = fullfile(repertoire, infos.name(i));
    img = fitsread(chemin);
    pupill = fonction_obtenir_pupille(img,1500);
    moy = fonction_intensite_moyenne(pupill);
    allmoy = [ allmoy , moy ];
end

Merci pour votre aide en tout cas!

[Bionics]

J'ai finalement compris pourquoi le i+2 était important, car autrement avec i et i+1 on traite la première image 2 ou 3 fois. Par contre dans la syntaxe tu t'es trompé :

Code :

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

chemin = fullfile(repertoire, infos(i+2).name);

Mais sinon ça a l'air de vouloir fonctionner sauf pour la partie que j'ai modifié :

Code :

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

allmoy = [ allmoy , moy ];

Merci beaucoup

[Gooby]

Le i+2 c'est parce que la fonction dir te renvoie . et .. en deux premières sorties.

Est-ce que la dernière ligne me permet de faire ça?

Oui, mais n'oublie pas d'initialiser ta variable allmoy

[Bionics]

Comme ceci?

Code :

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

allmoy = [];

avec ce code en dehors de la boucle, tout fonctionne c'est super, thanks!