Discussion :

[geoffroy31]

Bonjour !
[desole pour l ecriture, je nai pas mon clavier francais]

Jai deja beaucoup cherche sur Internet et sur le Forum mais n est malheuresement pas trouve mon bonheur...

Je suis actuellement en stage et m occupe d un projet de post traitement d une Simulation numerique de fluide en 3D.

Les donnees d entree sont des Points de coordonnees X,Y,Z caracterises par des proprietes (Pression partielle, Pression totale...).
Au lieu d exploiter l Ensemble du domaine 3D, on ne s occupe que de tranches decoupees pour une valeur de X, Y ou Y donnee.

L objectif est de comparer les tranches de deux simulations differentes.
Les 2 tranches representent donc le meme plan, mais les points qui y sont contenus ont des coordonnees differentes (donc impossible de simplement comparer les points 2 a 2).

Pour parer cela, voila comment je procede:
1. ouverture des deux fichiers de points, lecture des fichiers
2. creation d un meshgrid (qui sera commun aux deux jeux de donnees)
3. interpollation sur le meshgrid du premier Ensemble de Point [et du second Ensemble] Avec la fonction GridData
[4. calcule la difference pour chaque Point]
5. on affiche la difference Avec contourf et on voit ou les deux ecoulement different

J essaie deja d appliquer le procede a la premiere Serie de donnees, donc les parties [...] ne sont pas implementees.

Voila le Resultat :


Le Probleme: dans les Ensembles de points, il y a des espaces vides, ou "creux", sans points. Ceux ci correspondent a la position de lobjet dans la Simulation numerique. Il ne doit normalement donc pas y avoir de valeurs interpolees a cette endroit, puisqu il n y a pas de fluide.
mais gridData interpole quand meme des valeurs dans ces creux (necessairement car pour matlab pas de difference entre l Interieur et l Exterieur de l object)

Ma solution pour detecter ces creux: utilisation de la fonction alphavol ( voir le sujet : http://www.developpez.net/forums/d13...on-circulaire/ )

La fonction Alpha volume me cree donc des polynones et extrude la Zone interne vide.
J utilise ensuite la fonction inpolygon pour detecter les points de gridData a linterieur et a l Exterieur de ce polygone.



Pour les points a lexterieur du polygone, je remplace leurs valeurs dans GridData par NaN pour avoir une case blanche.

[ le fichier ne se charge pas, mais cest simplement la premiere Image avec la tache blanche de la seconde Image au Milieu]

Les nouveaux problemes qui se posent :
- la fonction inpolygon est beaucoup trop longue car il y a plusieurs milliers de points et plusieurs milliers de polygones issus de AlphaVol. Serait il possible de fusionner tous ces polygones en un seul et unique ? jai cherche et je crois que non
- meme Avec cela, il est impossible de detecter les tres Petits creux (mais je ne vois pas ici de solution etant donne que cest creux sont parfois plus Petit que la distance entre deux points dans les zones ou le maillage de la Simulation CFD est peu raffine...)
- gridData interpole toujours "n importe comment" mais cest simplement caché

le premier objectif du code est de tourner tres vite (il y a beaucoup de plans a comparer) et je pense donc que ma Methode Avec alphavol et inpolygon nest pas adapte...

je me suis donc demande sil etait possible de limiter le Rayon d Action de linterpollation appliquee par GridData (je veux dire: si l objet fait 10 Pixels, limiter l interpollation a 3 Pixels pour qu elle ne puisse pas se faire entre le points du haut de l objet et le Point du bas) ?

si quelquun aurait d autres idees auxquelles je n aurais pas pensé, merci

merci beaucoup pour l Aide que vous pourriez apporter
Geoffroy

[le fab]

salut

tu peux essayer de mettre nan pour les valeurs absentes (avant griddata)

Fabien

[Jerome Briot]

Tu trouveras des alternatives optimisées à inpolygon sur le File Exchange

[geoffroy31]

@ le fab : jai deja essaye de remplacer mes valeurs par des NaN avant le griddata, mais GridData na pas l´air de le supporter car elle me renvoyait une erreur.

@ Dut : merci pour le lien ! je ne pensais pas que d autres fonctions inpolygon existait !

[geoffroy31]

EDIT :

Ne pouvant pas utiliser les fonctions proposees par Dut car je navais pas Internet sur l ordinateur Avec Matlab, jai continué et mon Probleme sest résolu tout seul :

A la dernière étape, j affiche les valeurs de grid data Avec la fonction "contourf".
En parametrant correctement, on peut supprimer certains les Levels de basses valeurs et ainsi effacer les zones ou gridData a interpolé pour rien (car leurs valeurs sont très Faible comme il y a peu de Points).
Il y a toujours une erreur au paroi, mais Avec un suffisamment Grand nombre de Points, elle est limitée.

Je vous remercie cependant pour l Aide et la rapidité Avec laquelle elle a été apportée !

Bon Weekend