/*************************************************************************************************************************************/ /* */ /* 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