13.07.2015 Views

Rapport DEA Philippe Buhr - INSA de Lyon

Rapport DEA Philippe Buhr - INSA de Lyon

Rapport DEA Philippe Buhr - INSA de Lyon

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

<strong>DEA</strong> Images et SystèmesOptimisation <strong>de</strong> prédiction <strong>de</strong> couverture radiola source ailleurs au bas <strong>de</strong> la pyrami<strong>de</strong> puis la propager vers le haut pour enfin propagerl’ensemble <strong>de</strong>s flux vers le bas. Ceci est très avantageux car la majorité du temps <strong>de</strong> calcul estutilisé pour calculer les matrices <strong>de</strong> passage (nombreuses inversions matricielles sur <strong>de</strong>smatrices <strong>de</strong> très gran<strong>de</strong> taille) d’un niveau à un autre.Nous reviendrons plus loin sur le détail du fonctionnement du moteur <strong>de</strong> propagation.Ce moteur a été développé sous Matlab par J.-M.Gorce et a été porté sous environnementJava pas E. Jullo. La recherche <strong>de</strong> l’implantation <strong>de</strong>s algorithmes <strong>de</strong> création <strong>de</strong> la structurepyramidale, et <strong>de</strong> propagation <strong>de</strong>s sources ont été réalisées par E. Jullo. La structure <strong>de</strong>sdonnées a été choisie <strong>de</strong> manière à optimiser la place mémoire et le temps <strong>de</strong> calcul.E. Jullo montre qu’une telle approche permet la réduction du temps <strong>de</strong> calcul d’un facteurasymptotique significatif : pour un espace bidimensionnel <strong>de</strong> taille NxN, on a une complexité<strong>de</strong> calcul proportionnelle à N 2 .log 2 (N) au lieu <strong>de</strong> N 3 pour la métho<strong>de</strong> initiale dans le domainetemporel.Détail du déroulement d’un calcul <strong>de</strong> couverture :1. Pré-traitement : lecture dans un fichier <strong>de</strong>l’environnement modélisé, découpage <strong>de</strong>l’environnement <strong>de</strong> manière irrégulière selon laconstitution (matériau) <strong>de</strong> chaque pixel puiscréation <strong>de</strong> la pyrami<strong>de</strong> hiérarchique2. Calcul <strong>de</strong>s matrices <strong>de</strong> passage pour passer d’unniveau <strong>de</strong> la pyrami<strong>de</strong> au niveau directementsupérieur ou directement inférieur3. Lecture <strong>de</strong> la position <strong>de</strong> la source et remontée <strong>de</strong>la source jusqu’en haut <strong>de</strong> la pyrami<strong>de</strong> (calcul <strong>de</strong>la source équivalente au sommet)4. Re<strong>de</strong>scente le long <strong>de</strong> la pyrami<strong>de</strong> vers <strong>de</strong>s pixels(points <strong>de</strong> l’environnement) donnés ou versl’ensemble <strong>de</strong> l’environnement, et ce jusqu’à unerésolution donnée (résolution fixe ou blocshomogènes)Figure 3 : étapes 3 et 4 <strong>de</strong>propagation <strong>de</strong> la source vers lehaut puis re<strong>de</strong>scente vers la cible2.2.2 Limites <strong>de</strong> cette métho<strong>de</strong>Le moteur travaille à partir <strong>de</strong> matériaux (mur, air, vitre, etc.) dont il faut connaître lescaractéristiques, à savoir son coefficient d’atténuation et l’indice du milieu. Souvent cescaractéristiques ne peuvent pas être évaluées précisément <strong>de</strong> manière simple et risquent <strong>de</strong>changer fortement d’un bâtiment à un autre. Des facteurs assez variables influent cesparamètres : épaisseur <strong>de</strong>s murs, matériau du mur (plâtre, béton, brique), épaisseur <strong>de</strong>s vitreset type <strong>de</strong> verre,… On en vient donc tout naturellement à la conclusion que ces paramètresdoivent être i<strong>de</strong>ntifiés avant <strong>de</strong> calculer une carte <strong>de</strong> couverture dans un bâtiment. C’estexactement ce point qui fera l’objet principal <strong>de</strong> ce <strong>DEA</strong> : optimisation et i<strong>de</strong>ntification <strong>de</strong>sparamètres <strong>de</strong>s matériaux.La principale limitation <strong>de</strong> ce moteur <strong>de</strong> propagation est le temps <strong>de</strong> calcul ainsi que la placemémoire occupée par le programme. Les blocs et leurs matrices associées sont enregistrés surle disque dur <strong>de</strong> l’ordinateur. Ce sont principalement les matrices <strong>de</strong> propagation quiconsomment du temps <strong>de</strong> calcul. En effet, une inversion matricielle est réalisée.<strong>Philippe</strong> <strong>Buhr</strong> 12/42 2002/2003

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!