/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        T R A N S F O R M E E   D E   F O U R I E R   C O M P L E X E   I N V E R S E  :                                           */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xci/fft.03$K' :                                                                                                */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 1989??????????).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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
#include  image_image_IMAGESJ_EXT
#include  maths_trans_FOURIER_EXT
#include  image_image_QUAD_IMAGE_EXT

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        P A R A M E T R E S  :                                                                                                     */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#define   DECENTRER_LA_TRANSFORMEE_DE_FOURIER                                                                                           \
                    FAUX                                                                                                                \
                                        /* Faut-il decentrer de {0.5,0.5} en {0,0} la transformee de Fourier ('VRAI') ou pas         */ \
                                        /* ('FAUX') ?                                                                                */
#define   FACTEUR_MULTIPLICATIF_DE_LA_TRANSFORMEE_DE_FOURIER                                                                            \
                    FU                                                                                                                  \
                                        /* Facteur multiplicatif de la transformee de Fourier Resultat introduit le 20030506141238.  */

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

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        T R A N S F O R M E E   D E   F O U R I E R   C O M P L E X E   I N V E R S E  :                                           */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
     Bblock
     DEFV(CHAR,INIC(POINTERc(nom_imageR),NOM_PIPE));
                                        /* Nom de l'image  Resultat complexe.                                                        */
     DEFV(CHAR,INIC(POINTERc(nom_transformee_de_Fourier),NOM_PIPE));
                                        /* Nom de la transformee de Fourier que l'on va inverser.                                    */
     DEFV(Logical,INIT(decentrer_la_transformee_de_Fourier,DECENTRER_LA_TRANSFORMEE_DE_FOURIER));
                                        /* Faut-il decentrer de {0.5,0.5} en {0,0} la transformee de Fourier ('VRAI') ou pas         */
                                        /* ('FAUX') ?                                                                                */
     DEFV(Float,INIT(facteur_multiplicatif_de_la_transformee_de_fourier,FACTEUR_MULTIPLICATIF_DE_LA_TRANSFORMEE_DE_FOURIER));
                                        /* Facteur multiplicatif de la transformee de Fourier Resultat introduit le 20030506141238.  */
     /*..............................................................................................................................*/
     GET_ARGUMENTSi(nombre_d_arguments
                   ,BLOC(GET_ARGUMENT_C("transformee_de_fourier=""transformee_de_Fourier=""tf=""ft=""TF=""FT=""A="
                                       ,nom_transformee_de_Fourier
                                        );
                         GET_ARGUMENT_L("decentrer=",decentrer_la_transformee_de_Fourier);
                         GET_ARGUMENT_C("imageR=""R=",nom_imageR);
                         GET_ARGUMENT_F("multiplicateur=""facteur=",facteur_multiplicatif_de_la_transformee_de_fourier);
                         )
                    );

     begin_nouveau_block
          Bblock
          BDEFV(imageJ,transformee_de_Fourier_a_inverser);
                                        /* Image complexe contenant la transformee de Fourier que l'on va inverser...                */
          Test(PAS_D_ERREUR(CODE_ERROR(IloadJ_image(transformee_de_Fourier_a_inverser,nom_transformee_de_Fourier))))
               Bblock
                                        /* Acces a l'image dont on veut la transformee de Fourier inverse.                           */

               Test(IL_FAUT(decentrer_la_transformee_de_Fourier))
                    Bblock
                    BDEFV(imageJ,transformee_de_Fourier_decentree);
                                        /* Image complexe contenant la transformee de Fourier decentree...                           */
                    DEFV(deltaF_2D,translation_de_decentrage);
                                        /* Translation de decentrage...                                                              */

                    INITIALISATION_ACCROISSEMENT_2D(translation_de_decentrage
                                                   ,NEGA(_____lNORMALISE_OX(MOIT(dimX)))
                                                   ,NEGA(_____lNORMALISE_OY(MOIT(dimY)))
                                                    );
                                        /* Definition de la translation de decentrage...                                             */

                    CALS(IJtranslation(transformee_de_Fourier_decentree
                                      ,transformee_de_Fourier_a_inverser
                                      ,ADRESSE(translation_de_decentrage)
                                      ,VRAI
                                      ,VRAI
                                       )
                         );
                    CALS(IJmove(transformee_de_Fourier_a_inverser
                               ,transformee_de_Fourier_decentree
                                )
                         );
                                        /* Et on decale la transformee de Fourier...                                                 */

                    EDEFV(imageJ,transformee_de_Fourier_decentree);
                                        /* Image complexe contenant la transformee de Fourier decentree...                           */
                    Eblock
               ATes
                    Bblock
                    Eblock
               ETes

               CALS(IJfourier_2D_inverse_complexe(IJmageR,transformee_de_Fourier_a_inverser));
                                        /* Puis la transformee inverse.                                                              */

               Test(IFNE(facteur_multiplicatif_de_la_transformee_de_fourier,FU))
                                        /* Ce test, a priori superflu, est destine a etre sur de la compatibilite par rapport aux    */
                                        /* versions anterieures au 20030506141238 (de plus, il optimise un peu...).                  */
                    Bblock
                    CALS(Icomplexe_reelle(IFmage_PR,IJmageR));
                    CALS(Icomplexe_imaginaire(IFmage_PI,IJmageR));
                                        /* Recuperation des parties reelles et imaginaires.                                          */
                    CALS(IFscale(IFmage_PR
                                ,facteur_multiplicatif_de_la_transformee_de_fourier
                                ,IFmage_PR
                                ,FZERO
                                 )
                         );
                    CALS(IFscale(IFmage_PI
                                ,facteur_multiplicatif_de_la_transformee_de_fourier
                                ,IFmage_PI
                                ,FZERO
                                 )
                         );
                                        /* Application du facteur multiplicatif...                                                   */
                    CALS(Ireelle_complexe(IJmageR,IFmage_PR));
                    CALS(Iimaginaire_complexe(IJmageR,IFmage_PI));
                                        /* Generation de l'image Resultat complexe avec (IFmage_PR,IFmage_PI).                       */
                    Eblock
               ATes
                    Bblock
                    Eblock
               ETes

               CALi(IupdateJ_image(nom_imageR,IJmageR));
               Eblock
          ATes
               Bblock
               Test__CODE_ERREUR__ERREUR07;
               Eblock
          ETes

          EDEFV(imageJ,transformee_de_Fourier_a_inverser);
                                        /* Image complexe contenant la transformee de Fourier que l'on va inverser...                */

          Eblock
     end_nouveau_block

     RETU_Commande;
     Eblock
ECommande



Copyright © Jean-François COLONNA, 2019-2024.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / École polytechnique, Institut Polytechnique de Paris, 2019-2024.