Discussion :

[le fab]

Hello

je suis sur une optimisation (fmincon) dont le résultat est une surface
j'ai modifier la PlotFcns optimplotx pour avoir un bel affichage pendant l'optimisation

jusque là tout va bien
après, comme mon optimisation dure longtemps, j'ai eu envie de l'enregistrer au format video pour la voir à posteriori ... et la ca coince
mon code (avec en rouge mon enregistrement video) :

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
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
function stop = optimplotfmc(x,optimValues,state,varargin)

persistent writerObj

stop = false;
switch state
    case 'init'

        set(gcf,'doublebuffer','on')
        
        writerObj = VideoWriter('optimfmc.avi');
        writerObj.FrameRate = 1;
        open(writerObj);
        
    case 'done'
        
        close(writerObj);
        
        
    case 'iter'
%         % Reshape if x is a matrix
%         x = x(:);
%         xLength = length(x);
%         xlabelText = sprintf('Number of variables: %g',xLength);

        vx = evalin('base','vEcartHbassin');
        vy = evalin('base','vAmplCycle');
        
        x = reshape(x,numel(vy),numel(vx));
        
        [vx,vy] = meshgrid(vx,vy);
        
        

        % Display up to the first 100 values
%         if length(x) > 100
%             x = x(1:100);
%             xlabelText = {xlabelText,sprintf('Showing only the first 100 variables')};
%         end
        
        if optimValues.iteration == 0
            % The 'iter' case is  called during the zeroth iteration,
            % but it now has values that were empty during the 'init' case
            %plotx = plot3(vx(:),vy(:),x(:),'bo');
            plotx = surf(vx,vy,x,'edgecolor','r','FaceAlpha',0.66);
            colormap winter

            view(3)
            grid on
            title(sprintf('Energie = %.2f GWh',- optimValues.fval/1000),'interp','none');
            zlabel(sprintf('Current point'),'interp','none');
            ylabel('amplitude marée','interp','none');
            xlabel('ecart bassin','interp','none');
%             set(plotx,'edgecolor','none')
%             set(gca,'xlim',[0,1 + xLength])
            set(plotx,'Tag','optimplotfmc');
        else
            plotx = findobj(get(gca,'Children'),'Tag','optimplotfmc');
            set(plotx,'Zdata',x);
            title(sprintf('Energie = %.2f GWh',- optimValues.fval/1000),'interp','none');
        end
        
        drawnow
        frame = getframe(gcf);
        writeVideo(writerObj,frame);
        
end

le problème c'est que ca me récupère des frames presque vide
presque car j'ai un système d'axe quand même, mais 2D et pas du tout à l'échelle de ce que j'affiche (en fait c'est le système d'axe à la création de la figure plotFcns)

en pieces jointes la figure que j'essaye de capturer et la video resultat

une idée ??

Merci
Fab

[membreComplexe12]

salut le fab,

je ne suis pas bien habitué à faire ce genre de truc mais la seule fois où j'ai fais une vidéo j'ai procédé comme ceci :

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
aviobj1 = avifile('FOne.avi','fps',6,'quality',100)
 numframes1=length(F1(:,1));
%A=moviein(numframes1); % create the movie matrix
 set(gca,'NextPlot','replacechildren')
 %axis equal
 i=0;
 for i=1:numframes1
	plot_mesh(node,IEN,F1(i,:),1,MaxColorFv,MinColorFv);
    A(:,i)=getframe(gcf);
	F=A(:,i);
	aviobj1 = addframe(aviobj1,F);
 end
 %movie(A,10,3) % Play the MATLAB movie
%save movie.mat A 
aviobj1 = close(aviobj1);

j'espère que ça pourra t'aider...

[le fab]

salut

merci de ta réponse
ta méthode est assez similaire à la mienne quant à la capture (getframe) sauf que tu utilise avifile au lieu de VideoWriter

en fait mon problème est que le getframe ne prend pas ma fenetre en "photo" mais autre chose ... (je le vois bien en mode débug en regardant ce que retourne getframe : j'ai bien un soucis avec cette fonction et non avec VideoWriter)
j'ai déjà fait des videos de plot animés sans soucis, (avec la même méthode que dans mon code ci dessus), mais là, j'ai l'impression d'avoir un problème lié au fait qu'il ne s'agisse pas d'une figure "standard" mais d'une PlotFcns

Fabien

[membreComplexe12]

salut,

ah OK... je ne peux donc pas t'aider dans ce cas je ne sais même pas ce que c'est "PlotFcns " mais d'après ce que tu dis il semblerai que se soit lié à cette fonction...

bon courage

[le fab]

en fait les PlotFcns sont des fonctions d'affichage graphique "temps réel" d'une optimisation (fmincon, ...)

[membreComplexe12]

super intéressant ça !
faudrait que je jete un coup d'oeil un de ces quatres

[FLB]

Salut Fabien,
Je ne connaissais pas les PlotFcns, mais peut être que ceci t'aidera :
tu as essayé de remplacer gcf par le handle que te renvoie le h=figure associée aux PlotFcns ?

[le fab]

salut

en fait il n'y a pas de h=figure explicite dans la fonction
j'ai essayé gca aussi, mais ca ne marche pas non plus
et je ne pense pas que le problème vienne de là, car un print(gcf,'toto.jpg') fonctionne

bon là en même temps je n'ai plus trop de consacré à ce problème, je verrai plus tard

[le fab]

je m'auto répond en citant une réponse de Jerome dans un autre post
pas le temps de tester là maintenant mais c'est probablement ca, alors j'en garde ici une trace pour plus tard

Bonjour,

Je pense que le problème vient de la fonction getframe qui retourne toujours la première image. Ce problème est assez récurrent sous Windows.
Voir ce thread.