3. traitement d'image

3.1 lecture des fichiers

La première étape par laquelle nous avons dû passer est la lecture des fichiers contenant les données et leur conversion en images exploitables par le logiciel que nous utilisions, Aphelion.

3.1.1 description théorique

Le fichier fournit est de type .BDIR. Il consiste en 

    - un en-tête de deux fois quatre octets ( deux données de type unsigned long ) :

            - le premier contient le nombre de données présentes dans le fichier

             - le second donne le nombre de colonnes de l'image à obtenir.

    - données image de type float ( quatre octets ) que l'on doit associer par paire:

             - un float correspondant à la partie réelle

             - l'autre à la partie imaginaire d'un même nombre complexe.

Au lieu de fournir en sortie de l'opérateur de lecture une image complexe ou deux images, réelle et imaginaire, nous avons préféré calculer l'amplitude et la phase, c'est à dire le module et l'argument du nombre complexe. En effet, en imagerie radar, ces informations sont beaucoup plus intéressantes.

Amplitude

Phase

Cependant l'inverse de la  tangente d'un angle n'est compris qu'entre -p/2 et +p/2, modulo 2p. On lève cette indétermination en examinant le signe de la partie réelle.

3.1.2  description algorithmique

La création des fichiers .BDIR s'est faite sur des stations UNIX qui range les octets dans un sens alors que nous travaillions sur des PC qui utilisent la convention inverse. Nous avons donc créé une fonction, nommée "convert", qui inverse le sens des octets.

Nous avons choisi de présenter la phase sur un intervalle égal à [ 0 ; 2p] ainsi nous sommes certains qu'aucune image de départ n'a de niveaux de gris inférieurs à zéro. 

3.1.3  présentation des résultats obtenus

Pour créer chaque image radar, nous lisons deux fichiers: un d'amplitude et un d'interférogramme. On obtient donc quatre images. 

. Ci-dessous se trouve l'image synthétisée d'un seul et même bâtiment. Sur l'image d'amplitude et l'interférogramme on peut discerner trois parties:

        - le sol où l'amplitude est homogène et les bandes de phase se discernent bien

        - la zone de recouvrement où rien n'est structuré

        - le toit plat du bâtiment qui a les mêmes caractéristiques que le sol.

On va chercher la région correspondante au recouvrement. 

figure 1

. Ci-dessous, se trouve synthèse d'une zone urbaine. Elle peut être décomposer en plusieurs régions:

        - le sol où les caractéristiques sont les mêmes que précédemment

        - les zones de recouvrement correspondant aux pixels négatifs

        - quatre bâtiments mais dont le toit n'est pas forcément horizontal.

Ici le but sera de repérer les lignes de discontinuité.