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