#!/bin/csh

#######################################################################################################################################
#                                                                                                                                     #
#         T R A N S F O R M E E   E N   O N D E L E T T E   D E   M O R L E T                                                         #
#         D ' U N E   I M A G E   F L O T T A N T E  :                                                                                #
#                                                                                                                                     #
#                                                                                                                                     #
#         Utilisation :                                                                                                               #
#                                                                                                                                     #
#                   $xci/morlet.02$Z    <nomA> <nomR> [<Dilatation> [<Rotation>]]                                                     #
#                                                                                                                                     #
#                                                                                                                                     #
#         Author of '$xci/morlet.02$Z' :                                                                                              #
#                                                                                                                                     #
#                     Jean-Francois Colonna (LACTAMME, 20030415161242).                                                               #
#                                                                                                                                     #
#######################################################################################################################################

set       nomA=$1
set       nomR=$2
set       Dilatation=$3
set       Rotation=$4

setArgum  Dilatation                    0.1

setArgum  Rotation                      $pis2

$xci/morlet.02$X                                                                                                                        \
                    dilatation=$Dilatation                                                                                              \
                    rotation=$Rotation                                                                                                  \
                    R=$xT/ONDELETTE$J                                                                                                   \
                                                                                                                        $formatI
                                        # Generation de l'ondelette de Morlet demandee.                                               #
$xci/fftd$Z                                                                                                                             \
                    $xT/ONDELETTE                                                                                                       \
                    FAUX                                                                                                                \
                    $NEXIST
                                        # Calcul de la transformee de Fourier directe de l'ondelette de Morlet.                       #
$xci/init$X                                                                                                                             \
                    standard=FAUX                                                                                                       \
                    v=0                                                                                                                 \
                    R=$xT/ZERO                                                                                                          \
                                                                                                                        $formatI
$xci/fftdf$Z                                                                                                                            \
                    $nomA                                                                                                               \
                    $xT/ZERO                                                                                                            \
                    $xT/ARGUMENT
                                        # Calcul de la transformee de Fourier directe de l'image complexe {ImageArgument,Zero}.       #
$xci/fftp$Z                                                                                                                             \
                    $xT/ONDELETTE                                                                                                       \
                    $xT/ARGUMENT                                                                                                        \
                    $xT/PRODUIT
                                        # Produit de convolution des deux transformees de Fourier directes precedentes.               #
$xci/ffti$Z                                                                                                                             \
                    $xT/PRODUIT                                                                                                         \
                    FAUX                                                                                                                \
                    $NEXIST
                                        # Calcul de la transformee de Fourier inverse du produit precedent.                           #
$xci/J_RIMP$X                                                                                                                           \
                    A=$xT/PRODUIT$J                                                                                                     \
                    standard=FAUX                                                                                                       \
                    RR=$nomR$PR                                                                                                         \
                    RI=$nomR$PI                                                                                                         \
                    module_phase=VRAI                                                                                                   \
                    RM=$nomR$Mo                                                                                                         \
                    RP=$nomR$Ph                                                                                                         \
                                                                                                                        $formatI
                                        # Extraction des parties Reelle et Imaginaire de la transformee en ondelettes. Les images     #
                                        # {$Mo,$Ph} ont ete introduites le 20030415164841...                                          #

$DELETE   $xT/ONDELETTE*
$DELETE   $xT/ZERO
$DELETE   $xT/ARGUMENT*
$DELETE   $xT/PRODUIT*
                                        # Nettoyage...                                                                                #