Discussion :

[Supernul]

Bonjour à tous,

je débute sur MATLAB et je cherche à optimiser une fonction que j'ai écrite :

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
 
function [] = createjson (inputfile, outputdir)
    tic
    imp=load(inputfile,'input');
    cd('mat2html')
    mkdir(outputdir)
    cd('lib_init')
    copyfile('page.html', strcat('..', filesep, outputdir))
    cd('..')
    cd(outputdir)
 
    filename='buildings.js';
    XY=imp.input.XY;
    data=imp.input.data;
 
    dlmwrite(filename, ['var buildings={"crs":{"type":"name","properties":{"name":"urn:ogc:def:crs:EPSG:2154"}},"type":"FeatureCollection","features":['], 'delimiter', '')
 
    for i=1:numel(XY);
        x=XY{i}(:,1);
        y=XY{i}(:,2);
 
        dlmwrite(filename, ['{"type":"Feature","geometry":{"type":"Polygon","coordinates":['], '-append', 'delimiter', '')
        s = '[';
        for j=1:numel(x);
            if (~isnan(x(j)) & ~isnan(y(j)));
        	   s =  strcat(s, strcat('[', num2str(round(x(j), 1)), ',', num2str(round(y(j), 1)), '],'));
            end
        end
        s = strcat(s(1:end-1), ']');
        dlmwrite(filename, [s,']},"properties":{"solar":', num2str(data(i)), '}},'], '-append', 'delimiter', '')
    end
    dlmwrite(filename, [']}'], '-append', 'delimiter', '')
    dlmwrite(filename, ['var min_val=', num2str(min(data))], '-append', 'delimiter', '')
    dlmwrite(filename, ['var max_val=', num2str(max(data))], '-append', 'delimiter', '')
 
    toc
    t = toc;
   cd('../..') 
end

Le principe est d'écrire un fichier json à partir une structure qui contient des coordonnées XY et des données data qui sera affichée dans une page html.

Travaillant avec R et python je me pose quelques questions :
- éviter les boucles (comme dans R) peut-il faire gagner du temps ?
- j'ai essayé des versions favorisant dlmwrite (écriture dans le fichier à chaque fois) ou strcat (mise en mémoire puis écriture en une seule fois). Dans les deux cas je double le temps de traitement. Il y-a-t-il des fonctions plus rapides ?

Merci de votre aide

[Jerome Briot]

Essaie 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
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
function createjson(inputfile, outputdir)
 
outputdir = ['./mat2html/', outputdir];
 
mkdir(outputdir)
 
copyfile('./lib_init/page.html', outputdir)
 
imp = load(inputfile,'input');
 
fid = fopen(fullfile(outputdir, 'buildings.js'), 'wt');
 
fprintf(fid, 'var buildings={"crs":{"type":"name","properties":{"name":"urn:ogc:def:crs:EPSG:2154"}},"type":"FeatureCollection","features":[');
 
for i = 1:numel(imp.input.XY);
 
    x = imp.input.XY{i}(:,1);
    y = imp.input.XY{i}(:,2);
 
    fprintf(fid, '{"type":"Feature","geometry":{"type":"Polygon","coordinates":[[');
 
    idx = ~isnan(x) & ~isnan(y);
 
    fprintf(fid, '[%f,%f],', [x(idx(1:end-1)) y(idx(1:end-1))].');
 
    fprintf(fid, '[%f,%f]]', x(idx(end)), y(idx(end)));
 
    fprintf(fid, ']},"properties":{"solar":%f}}', imp.input.data(i));
 
end
 
fprintf(fid, ']}\nvar min_val=%f\nvar max_val=%f', min(imp.input.data), max(imp.input.data));
 
fclose(fid);

C'est un début de solution sans doute avec quelques corrections à effectuer.
À toi de l'améliorer en fonction des données réelles que tu utilises.

Si tu n'y arrives pas, donne nous un fichier mat (mis dans une archive zip ou rar) avec la structure imp

[Supernul]

Bonjour,

merci beaucoup. Je n'en attendais pas tant. J'ai modifié dans petites choses dans le code proposé en raison d'une imprécision dans mon premier message sur la structuration de mes dossiers et le code fonctionne parfaitement.

je vais approfondir de mon côté la différence entre fprintf et dlmwrite pour comprendre la différence (de vitesse surtout).

Il y a une ligne que je ne suis pas certain de comprendre : idx = ~isnan(x) & ~isnan(y);
Si je comprends bien, on affecte à idx une liste binaire qui permet ensuite de retenir les x et y pour lesquels idx vaut 1.

L'écart sur ma donnée est énorme : je passe de 181 ms à 0,4 ms !!!

Il va falloir que je change un peu ma façon de penser pour une approche plus matricielle.

Encore merci.

[Jerome Briot]

Il y a une ligne que je ne suis pas certain de comprendre : idx = ~isnan(x) & ~isnan(y);
Si je comprends bien, on affecte à idx une liste binaire qui permet ensuite de retenir les x et y pour lesquels idx vaut 1.

Voir la FAQ : Qu'est-ce que l'indexage logique ?