Discussion :

[pootchi]

Bonjour,

Je suis sur du traitement d'images, plus précisément du face morphing.
Mon programme est un peu long à s'exécuter à mon goût. J'ai déjà
fait un PROFILE pour voir ce qui prenait le plus de temps, donc voilà
ce qui me pose problème :

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
 
coeff =    [1 oeil_left(1)+diff_pt1(1)  oeil_left(1)+diff_pt1(1);...
            1 1                         oeil_left(3)+diff_pt1(2);]
          /...
           [1 oeil_left(1)              oeil_left(1);...
            1 1                         oeil_left(3);...
            1 1                         1];
 
tform = maketform('affine', [coeff(1) coeff(3) 0; coeff(2) coeff(4) 0; 0 0 1]);
tri = src(1:oeil_left(3), 1:oeil_left(1), :);
tmp = imtransform(tri,tform);
tmp_size = size(tmp);
 
tmp_src(1:tmp_size(1) - 1, 1:tmp_size(2) - 1,:)=...
                    tmp(1:tmp_size(1) - 1, 1:tmp_size(2) - 1,:);

Ce qui prend le plus de temps est l'appel à la fonction imtransform,
comme c'est une fonction matlab, je doute que je puisse y faire grand
chose

Cependant, je dois répéter ce petit bout de code beaucoup de fois, et
les matrices que j'utilise pour calculer coeff sont toutes différentes.

Donc, au final, je me retrouve avec un fichier immense, avec ce petit bout
de code, quasi recopier, une trentaine de fois.

Est-ce qu'il y aurait moyen de réduire tout ça ?

* Est-ce que ça vaut le coup que je passe tout le fichier en C ?
* Je pourrais faire une grosse matrice regroupant tous mes "coeff"
et seulement ensuite faire un seul appel à imtransform, du genre :

Code :

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

1
2
3
4
tform = maketform('affine', liste_coeff(:));
tri = src(1:oeil_left(3), 1:oeil_left(1), :);
tmp = imtransform(tri,tform(:));

ça serait possible ? ou alors ça servirait à rien ?

Merci,

--
Aurélie Barbier

[Jerome Briot]

Mon programme est un peu long à s'exécuter à mon goût.

Combien de temps ?

[pootchi]

8.875sec en exécutant le bout de code 120 fois
c'est imtransform qui prend 85% du temps

[Jerome Briot]

Tu pourrais nous montrer le premier tableau récapitulatif de PROFILE ?

[pootchi]

Voilà (je l'ai fait tourné un peu plus longtemps...)



--
Aurélie Barbier

[Jerome Briot]

Je vois que REPMAT est appelé 2000 fois

Peux-tu montrer la (ou les) ligne(s) de code où cette fonction apparait.

Sinon, la fonction la plus couteuse (barre bleue foncée) est un fichier MEX (environ 7s). Ton code est donc quasi impossible à optimiser (sous MATLAB), à mon avis.

[pootchi]

Je vois que REPMAT est appelé 2000 fois

Peux-tu montrer la (ou les) ligne(s) de code où cette fonction apparait.

C'est imtransform qui l'appelle à chaque fois

Sinon, la fonction la plus couteuse (barre bleue foncée) est un fichier MEX (environ 7s). Ton code est donc quasi impossible à optimiser (sous MATLAB), à mon avis.

Ouaip, c'est ce que je me disais...Je me demandais si en C ça serait mieux...?

--
Aurélie Barbier

[Jerome Briot]

Ouaip, c'est ce que je me disais...Je me demandais si en C ça serait mieux...?

C'est déjà du C (a priori)

=> : Qu'est ce qu'un fichier MEX ?

[pootchi]

C'est déjà du C (a priori)

Yep, la fatigue, tout ça

Bien, merci, je vais laisser mon code comme ça alors...

--
Aurélie Barbier

[paradize3]

Hello,


tu pourrais appeler imtransform avec une 'Nearest-neighbor interpolation' a la place de bilineaire (defaut), si ca te convient.

[paradize3]

attention: les conseils suivants ne doivent etre applique que si a. les performances sont vraiment critiques b. tu es sur que ton algo est bien pense a un niveau "plus haut". c. tu as pas mal de temps a disposition!

Par rapport a l'implementation en C, il est tres probable que celle proposee ne soit pas multi-thredees... Donc il est peut-etre possible de faire une meilleure implementation sur proc multi-core. (mais a verifier avant de se lancer!).

Tu peux aussi aller voir les sources de resampsep (elles sont fournies avec la toolbox je crois). Il y a eventuellement des choses a optimiser dedans.... Ou a adapter a ton cas particulier.

Voila,

Bonne chance,

Gregoire

[pootchi]

attention: les conseils suivants ne doivent etre applique que si a. les performances sont vraiment critiques b. tu es sur que ton algo est bien pense a un niveau "plus haut". c. tu as pas mal de temps a disposition!

Par rapport a l'implementation en C, il est tres probable que celle proposee ne soit pas multi-thredees... Donc il est peut-etre possible de faire une meilleure implementation sur proc multi-core. (mais a verifier avant de se lancer!).

Tu peux aussi aller voir les sources de resampsep (elles sont fournies avec la toolbox je crois). Il y a eventuellement des choses a optimiser dedans.... Ou a adapter a ton cas particulier.

Ok, merci beaucoup pour ces indications. Comme c'est un projet de stage,
ce n'est pas grave si je n'optimise pas, mais quand même

Je vais donc réfléchir avant de me lancer dans cette optim

Merci en tout cas !

--
Aurélie Barbier