/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        L O G A R I T H M E   D E C I M A L   ( O U   N E P E R I E N )   D ' U N E   I M A G E   F L O T T A N T E  :             */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Nota :                                                                                                                     */
/*                                                                                                                                   */
/*                    Ce programme s'appelle "decibel", mais                                                                         */
/*                  il s'agit evidemment d'un abus de langage                                                                        */
/*                  puisque le facteur '10' associe a la definition                                                                  */
/*                  du decibel est ici absent. Rappelons que le                                                                      */
/*                  nombre de decibels de deux puissances P1 et                                                                      */
/*                  P2 est egal a 10.log(P1/P2), le logarithme                                                                       */
/*                  etant decimal.                                                                                                   */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xci/decibel.01$K' :                                                                                            */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 1997??????????).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        I N T E R F A C E   ' listG '  :                                                                                           */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        :Debut_listG:                                                                                                              */
/*        :Fin_listG:                                                                                                                */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        D I R E C T I V E S   S P E C I F I Q U E S   D E   C O M P I L A T I O N  :                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        F I C H I E R S   D ' I N C L U D E S  :                                                                                   */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  INCLUDES_BASE
#include  image_image_IMAGESF_EXT

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        P A R A M E T R E S  :                                                                                                     */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#define   FORCER_LE_SEUIL_DES_VALEURS                                                                                                   \
                    VRAI                                                                                                                \
                                        /* Le seuil 'seuil_des_valeurs' doit-il etre force ('VRAI') ou etre le minimum des           */ \
                                        /* valeurs ('FAUX') ? Ceci fut introduit le 20070621140309...                                */

#define   SEUIL_DES_VALEURS                                                                                                             \
                    GRAND_EPSILON                                                                                                       \
                                        /* Seuil des valeurs en-deca duquel on ramene a ce meme seuil...                             */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        M A C R O S   U T I L E S  :                                                                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        L O G A R I T H M E   D E C I M A L   ( O U   N E P E R I E N )   D ' U N E   I M A G E   F L O T T A N T E  :             */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
     Bblock
     DEFV(CHAR,INIC(POINTERc(nom_imageR),NOM_PIPE));
     DEFV(CHAR,INIC(POINTERc(nom_imageA),NOM_PIPE));

     DEFV(Logical,INIT(forcer_le_seuil_des_valeurs,FORCER_LE_SEUIL_DES_VALEURS));
                                        /* Le seuil 'seuil_des_valeurs' doit-il etre force ('VRAI') ou etre le minimum des           */
                                        /* valeurs ('FAUX') ? Ceci fut introduit le 20070621140309...                                */
     DEFV(Float,INIT(seuil_des_valeurs,SEUIL_DES_VALEURS));
                                        /* Seuil des valeurs en-deca duquel on ramene a ce meme seuil...                             */
     /*..............................................................................................................................*/
     GET_ARGUMENTSi(nombre_d_arguments
                   ,BLOC(GET_ARGUMENT_C("imageA=""A=",nom_imageA);
                         GET_ARGUMENT_C("imageR=""R=",nom_imageR);

                         GET_ARGUMENT_L("forcer=",forcer_le_seuil_des_valeurs);
                         GET_ARGUMENT_N("minimum=",forcer_le_seuil_des_valeurs);
                         GET_ARGUMENT_F("seuil=",seuil_des_valeurs);

                         GET_ARGUMENT_L("log10=""LO1X=",IFdynamique_logarithmique_decimal_avec_seuillage_et_translation_____LO1X);
                         GET_ARGUMENT_N("log=""LO1G=",IFdynamique_logarithmique_decimal_avec_seuillage_et_translation_____LO1X);
                         )
                    );

     CALi(IFinitialisation(IFmageR,FZERO));
                                        /* Initialisation de l'image Resultat.                                                       */

     Test(PAS_D_ERREUR(CODE_ERROR(IloadF_image(IFmageA,nom_imageA))))
          Bblock
          Test(IL_NE_FAUT_PAS(forcer_le_seuil_des_valeurs))
                                        /* Test introduit le 20070621140309...                                                       */
               Bblock
               DEFV(genere_Float,INIT(niveau_minimum,FLOT__NIVEAU_UNDEF));
               DEFV(genere_Float,INIT(niveau_maximum,FLOT__NIVEAU_UNDEF));

               CALS(IFnivo_extrema(IFmageA
                                  ,ADRESSE(niveau_minimum)
                                  ,ADRESSE(niveau_maximum)
                                   )
                    );
                                        /* Recherche des extrema...                                                                  */
               EGAL(seuil_des_valeurs,niveau_minimum);
                                        /* Et forcage du seuil des valeurs...                                                        */
               Eblock
          ATes
               Bblock
               Eblock
          ETes

          Test(IZGT(seuil_des_valeurs))
               Bblock
               CALS(IFdynamique_logarithmique_decimal_avec_seuillage_et_translation
                        (IFmageR
                        ,IFmageA
                        ,seuil_des_valeurs
                        ,COND(EST_VRAI(IFdynamique_logarithmique_decimal_avec_seuillage_et_translation_____LO1X)
                             ,NEGA(LO1X(seuil_des_valeurs))
                             ,NEGA(LOGX(seuil_des_valeurs))
                              )
                         )
                    );
                                        /* Prise du logarithme decimal en ramenant tout a 0...                                       */
                                        /*                                                                                           */
                                        /* Le 20070621134728, 'LOGX(...)' fut ajoutee a 'LO1X(...)'...                               */
               Eblock
          ATes
               Bblock
               PRINT_ERREUR("le seuil des valeurs doit etre strictement positif");
               Eblock
          ETes

          CALi(IupdateF_image(nom_imageR,IFmageR));
          Eblock
     ATes
          Bblock
          Test__CODE_ERREUR__ERREUR07;
          Eblock
     ETes

     RETU_Commande;
     Eblock
ECommande