Génération d'isosurface avec Gnuplot lente

guillaumech

Expert
Bonsoir,



Bonsoir à tous !

Je vous propose une petite question Gnuplot pour changer :) Voici mon soucis, j'ai le script ci-dessous qui me génère des isosurface. Tout fonctionne comme je l'aimerai, cependant le traitement est très lent. Mon fichier de données est assez conséquent, du coup je m'en remets aux experts si des fois il y avait des erreurs qui ralentiraient le traitement.

Code:
set term png transparent size 7200,3610 font "DejaVuSerif-Bold" 8
set output './01.png'
set dgrid3d 361,720,1
set isosamples 361,720
unset surf
set table './processing/isoline/PRMSLmeansealevelTable'
splot './processing/isoline/PRMSLmeansealevel'
unset table
set contour base
set cntrparam order 8
set cntrparam bspline
set cntrparam levels incremental 95000, 100, 150000
unset surf
set table './processing/isoline/PRMSLmeansealevelTable2'
splot './processing/isoline/PRMSLmeansealevel'
unset table
reset
unset key
set xrange [0:360]
set yrange [-90:90]
set contour base
set cntrparam bspline
unset xtics
unset ytics
unset ztics
unset surf
l '<./processing/config/./cont.sh ./processing/isoline/PRMSLmeansealevelTable2 0 10 0
p './processing/isoline/PRMSLmeansealevelTable' w ima, '<./processing/config/./cont.sh ./processing/isoline/PRMSLmeansealevelTable2 2 0 0' w l lt 8 lw 1

Le fichier de données initial comporte 360 000 lignes. Je peux le fournir pour ceux qui le souhaitent. Le fichier cont.sh dans le script est utile pour afficher les valeurs de chaque courbes de niveau trouvées. Le traitement global est de plus de 10 min sur un bon matériel. Cela me parait énorme, mais peut-être que je me trompe.
Les commandes lentes sont les splot.

Merci à vous tous en tout cas.
Guillaume.

PS : je cherche à tracer des isolines de pressions atmosphériques, données que j'ai dans mon fichier initial.
 

guillaumech

Expert
Hey, un des gentils dinosaures du forum :)

Oui, toujours, tombé dedans quand j'étais petit :) Tu n'es pas à fond Gnuplot mon cher Zeb ? Remarques, j'élargie la question, je comptes convertir une grille de données en image d'isoline, comprendre par là, chaque même valeur de la grille se voit couvert par une ligne (souvent avec une interpolation spline). J'ai regardé instinctivement du côté de Gnuplot, mais il existe peut-être d'autres outils, ou même un moyen de le coder rapidement (afin d'aller plus vite que Gnuplot).
 

zeb

Modérateur
Je ne suis pas spécialiste, mais si j'avais à le faire, je regarderais du côté de GNUplot :/

-----

/private
Tiens M. Météo, j'ai eu hier le plaisir d'observer hier en fin d'après-midi sur Paris un très joli phénomène : un parhélie !
 

guillaumech

Expert
Oui, moi aussi. J'ai même tenté du côté d'ImageMagick, avec un resize de la grille pour créer une interpolation correcte et ensuite faire de la détection de contour de zones, mais ce n'est pas très concluant ... Traiter 365 000 points en 15min ne te parait pas étrange ?
Je vais attendre Batchy qui a déjà utilisé cet outil, on verra.

Chanceux pour le parhélie, c'est peu fréquent !
 
Vous devez vous inscrire ou vous connecter pour répondre ici.
Derniers messages publiés
Statistiques globales
Discussions
730 125
Messages
6 717 769
Membres
1 586 361
Dernier membre
Florian3549
Partager cette page
Haut