Chapitre IV : Etude de bruit de fond
37
Figure 12: Emplacement des enregistrements H/V sur fond d’une image satellitaire
IV.5. Traitement des données
Nous avons récupéré les données des enregistrements, via la carte flash, sur PC avec le
logiciel ReadCity.
Les courbes H/V
obtenus ont été calculées à l’aide du logiciel libre Geopsy
[www.geopsy.org]. Ce logiciel (figure 13) a été développé dans le cadre du projet européen
SESAME [2004]. L’algorithme de calcul du rapport spectral H/V est le suivant :
1) sélection des fenêtres de temps stationnaires, de durée au moins égale à 10 fois la période
minimale que l'on souhaite analyser.
2) calcul sur chaque fenêtre de temps n, les spectres de Fourier des signaux enregistrés sur les
trois composantes (deux horizontales H1 et H2, et une verticale V).
3) lissage des spectres avec la procédure de fenêtre glissante à largeur adaptée à la fréquence
d'analyse.
4) calcul de la moyenne quadratique des spectres horizontaux lissés SH(Ni) (composantes
Nord-Sud et Est-Ouest) pour chaque fenêtre :
( ) √
( )
( )
(2)
5) calcul des rapports spectraux H/V(Ni) pour chaque fenêtre :
( ) (
( )
( )
)
(3)
37
Figure 12: Emplacement des enregistrements H/V sur fond d’une image satellitaire
IV.5. Traitement des données
Nous avons récupéré les données des enregistrements, via la carte flash, sur PC avec le
logiciel ReadCity.
Les courbes H/V
obtenus ont été calculées à l’aide du logiciel libre Geopsy
[www.geopsy.org]. Ce logiciel (figure 13) a été développé dans le cadre du projet européen
SESAME [2004]. L’algorithme de calcul du rapport spectral H/V est le suivant :
1) sélection des fenêtres de temps stationnaires, de durée au moins égale à 10 fois la période
minimale que l'on souhaite analyser.
2) calcul sur chaque fenêtre de temps n, les spectres de Fourier des signaux enregistrés sur les
trois composantes (deux horizontales H1 et H2, et une verticale V).
3) lissage des spectres avec la procédure de fenêtre glissante à largeur adaptée à la fréquence
d'analyse.
4) calcul de la moyenne quadratique des spectres horizontaux lissés SH(Ni) (composantes
Nord-Sud et Est-Ouest) pour chaque fenêtre :
( ) √
( )
( )
(2)
5) calcul des rapports spectraux H/V(Ni) pour chaque fenêtre :
( ) (
( )
( )
)
(3)
