/*************************************************************************************************************************************/
/* */
/* C O N V O L U T I O N D E L ' I M A G E A R G U M E N T P O U R B O U C H E R L E S T R O U S : */
/* */
/* */
/* Definition : */
/* */
/* Cette commande genere une image */
/* dont le nom est le premier argument */
/* d'appel ; elle est le resultat de */
/* la convolution de la deuxieme image. */
/* Elle permet, une image ayant ete */
/* calculee avec (pasX,pasY) # (PasX,PasY) */
/* de "boucher" les trous ainsi crees par */
/* une interpolation ; il suffit de prendre */
/* (pasX,pasY) egale a (pasX,pasY) qui a */
/* servi a la generer... */
/* */
/* */
/* Author of '$xci/convol.04$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 1987??????????). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* P A R A M E T R E S : */
/* */
/*************************************************************************************************************************************/
#define SEUIL_DE_FAUX_VRAI_POUR_LES_NIVEAUX_A_TRAITER \
GRIS_4 \
/* Seuil de discrimination {FAUX,VRAI} introduit le 20060304185627 lors de l'entree de la */ \
/* liste des niveaux a traiter via une palette de couleur. */
#define EXCLURE_UN_NIVEAU_DE_LA_CONVOLUTION \
FAUX
#define NIVEAU_A_EXCLURE_DE_LA_CONVOLUTION \
NOIR
/* Faut-il ('VRAI') ou pas ('FAUX') exclure un niveau de la convolution (par defaut le */
/* NOIR) ? Cela fut introduit le 20060303101605... */
#define NOMBRE_DE_POINTS_DU_NOYAU \
EXP2(DOUP(UN)) \
/* Nombre de points du noyau. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* M A C R O S U T I L E S : */
/* */
/*************************************************************************************************************************************/
#include xci/convol.01.I"
/* Introduits le 20060304190400 pour 'SELECTION_DES_NIVEAUX(...)'. */
#define VALEUR_INITIALE_DU_NOYAU \
FRA10(FU) \
/* Valeur initiale de tous les elements du noyau qui pourra etre ensuite retouchee... */
#redefine NOYAU(numero,valeur,incrementer_le_numero) \
Bblock \
EGAL(ITb1(inhibition_du_noyau,INDX(numero,PREMIER_POINT)),ACTIF); \
/* Activation systematique du noyau... */ \
Test(IFLE(numero \
,LSTX(PREMIER_POINT \
,NOMBRE_DE_POINTS_EFFECTIF_D_UN_NOYAU_DE_CONVOLUTION_PARCOURU_CIRCULAIREMENT(nombre_de_points) \
) \
) \
) \
Bblock \
INCR(facteur_multiplicatif,valeur); \
/* Le facteur multiplicatif est egal a la somme des ponderations... */ \
Eblock \
ATes \
Bblock \
Eblock \
ETes \
Test(IL_FAUT(incrementer_le_numero)) \
Bblock \
EGAL(ITb1(noyau,INDX(numero,PREMIER_POINT)),valeur); \
Test(IFLE(numero \
,LSTX(PREMIER_POINT \
,NOMBRE_DE_POINTS_EFFECTIF_D_UN_NOYAU_DE_CONVOLUTION_PARCOURU_CIRCULAIREMENT(nombre_de_points) \
) \
) \
) \
Bblock \
DECR(facteur_multiplicatif,VALEUR_INITIALE_DU_NOYAU); \
/* Correction du facteur multiplicatif due a la valeur qui a ete forcee precedemment... */ \
Eblock \
ATes \
Bblock \
Eblock \
ETes \
INCR(numero,I); \
Eblock \
ATes \
Bblock \
EGAL(ITb1(noyau,INDX(numero,PREMIER_POINT)),VALEUR_INITIALE_DU_NOYAU); \
Eblock \
ETes \
Eblock \
/* Initialisation element par element du noyau. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* C O N V O L U T I O N D E L ' I M A G E A R G U M E N T P O U R B O U C H E R L E S T R O U S : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
DEFV(CHAR,INIC(POINTERc(nom_imageR),NOM_PIPE));
DEFV(CHAR,INIC(POINTERc(nom_imageA),NOM_PIPE));
DEFV(genere_p,INIT(seuil_de_FAUX_VRAI_pour_les_niveaux_a_traiter,SEUIL_DE_FAUX_VRAI_POUR_LES_NIVEAUX_A_TRAITER));
DEFV(CHAR,INIC(POINTERc(nom_paletteA),NOM_PIPE));
/* Nom de la palette a utiliser eventuellement et seuil de discrimination {FAUX,VRAI} */
/* introduits le 20060304185627... */
DEFV(Logical,INIT(exclure_un_niveau_de_la_convolution,EXCLURE_UN_NIVEAU_DE_LA_CONVOLUTION));
DEFV(genere_p,INIT(niveau_a_exclure_de_la_convolution,NIVEAU_A_EXCLURE_DE_LA_CONVOLUTION));
/* Faut-il ('VRAI') ou pas ('FAUX') exclure un niveau de la convolution (par defaut le */
/* NOIR) ? Cela fut introduit le 20060303101605... */
DEFV(Logical,DTb1(niveaux_a_traiter,COULEURS));
/* Definit les niveaux sur lesquels on doit faire la convolution par 'Pconvolution()'. */
DEFV(Logical,DTb1(niveaux_cumulables,COULEURS));
/* Definit les niveaux cumulables lors du calcul de 'Pconvolution_____cumul_courant'. */
DEFV(Float,INIT(facteur_multiplicatif,FZERO));
/* Facteur multiplicatif du produit de convolution en chaque point {X,Y}. */
DEFV(Int,INIT(nombre_de_points,NOMBRE_DE_POINTS_DU_NOYAU));
/* Nombre de points du noyau. */
DEFV(Float,DTb1(noyau,TAILLE_MAXIMALE_D_UN_NOYAU_DE_CONVOLUTION));
/* Noyau de la convolution, */
DEFV(Logical,DTb1(inhibition_du_noyau,TAILLE_MAXIMALE_D_UN_NOYAU_DE_CONVOLUTION));
/* Et sa liste d'activite. */
DEFV(Int,INIT(index,UNDEF));
/* Index d'initialisation du noyau. */
/*..............................................................................................................................*/
GET_ARGUMENTSi(nombre_d_arguments
,BLOC(GET_ARGUMENT_C("imageA=""A=",nom_imageA);
GET_ARGUMENT_C("imageR=""R=",nom_imageR);
GET_ARGUMENT_C("paletteA=""palette=""pA=",nom_paletteA);
/* ATTENTION : le parametre "p=" ne peut etre defini ici a cause de 'nombre_de_points'... */
GET_ARGUMENT_P("sp=""seuil_palette=",seuil_de_FAUX_VRAI_pour_les_niveaux_a_traiter);
GET_ARGUMENT_L("exclure=",exclure_un_niveau_de_la_convolution);
GET_ARGUMENT_P("niveau_a_exclure=""ne=",niveau_a_exclure_de_la_convolution);
GET_ARGUMENT_L("circulaire=",Pconvolution_____parcourir_circulairement_les_spirales_carrees);
/* Option introduite le 20040703110028... */
GET_ARGUMENT_L("circulaire_optimale="
,Pconvolution_____adapter_le_nombre_de_points_lors_du_parcours_circulaire_des_spirales_carrees
);
/* Option introduite le 20061224122147... */
GET_ARGUMENT_L("tore=",Pconvolution_____placer_l_image_sur_un_tore);
GET_ARGUMENT_L("ponderations_utiles=",Pconvolution_____normaliser_uniquement_avec_les_ponderations_utilisees);
GET_ARGUMENT_I("pmX=",Pconvolution_____pasX_de_la_mosaique_de_convolution);
GET_ARGUMENT_I("pmY=",Pconvolution_____pasY_de_la_mosaique_de_convolution);
GET_ARGUMENT_I("tmX=",Pconvolution_____translationX_de_la_mosaique_de_convolution);
GET_ARGUMENT_I("tmY=",Pconvolution_____translationY_de_la_mosaique_de_convolution);
GET_ARGUMENT_I("saut=",Pconvolution_____nombre_de_points_sautes);
GET_ARGUMENT_F("transfert=",Pconvolution_____exposant_de_la_fonction_de_transfert);
GET_ARGUMENT_I("spirale_horizontal=",SPIRALE_DEFINITION_____facteur_delta_horizontal);
GET_ARGUMENT_I("spirale_vertical=",SPIRALE_DEFINITION_____facteur_delta_vertical);
GET_ARGUMENT_I("points=""p=",nombre_de_points);
)
);
CALi(Inoir(ImageA));
/* Initialisation de l'image Resultat avant le changement eventuel de (pasX,pasY), mais */
/* ATTENTION, cette initialisation a lieu dans 'ImageA'... */
SELECTION_DES_NIVEAUX(nom_paletteA,niveaux_a_traiter,niveaux_cumulables);
/* Introduit sous cette forme le 20060304185627... */
Test(IFLE(NOMBRE_DE_POINTS_EFFECTIF_D_UN_NOYAU_DE_CONVOLUTION_PARCOURU_CIRCULAIREMENT(nombre_de_points)
,TAILLE_MAXIMALE_D_UN_NOYAU_DE_CONVOLUTION
)
)
Bblock
DoIn(index
,PREMIER_POINT
,LSTX(PREMIER_POINT,NOMBRE_DE_POINTS_EFFECTIF_D_UN_NOYAU_DE_CONVOLUTION_PARCOURU_CIRCULAIREMENT(nombre_de_points))
,I
)
Bblock
NOYAU(index,VALEUR_INITIALE_DU_NOYAU,FAUX);
/* Initialisation du noyau de convolution ; tous les elements a une meme valeur... */
Eblock
EDoI
Eblock
ATes
Bblock
PRINT_ERREUR("la taille demandee pour le noyau de convolution est incompatible avec les definitions");
EGAL(nombre_de_points,NOMBRE_DE_POINTS_DU_NOYAU);
CAL1(Prer1("(la valeur %d par defaut est forcee)\n",nombre_de_points));
/* Introduit le 20111031140822 car, en effet, manquait cruellement... */
Eblock
ETes
EGAL(index,PREMIER_POINT);
/* Depart sur le premier element afin de corriger ponctuellement l'initialisation ci-dessus. */
NOYAU(index,FU,VRAI);
NOYAU(index,FDU,VRAI);
NOYAU(index,FQU,VRAI);
NOYAU(index,FDU,VRAI);
NOYAU(index,FQU,VRAI);
NOYAU(index,FDU,VRAI);
NOYAU(index,FQU,VRAI);
NOYAU(index,FDU,VRAI);
NOYAU(index,FQU,VRAI);
/* Initialisation du noyau de convolution : */
/* */
/* 0.25 0.50 0.25 */
/* */
/* 0.50 1.00 0.50 */
/* */
/* 0.25 0.50 0.25 */
/* */
Test(PAS_D_ERREUR(CODE_ERROR(Iload_image(ImageR,nom_imageA))))
/* Chargement de l'image a transformer, mais ATTENTION, le chargement a lieu dans 'ImageR'. */
Bblock
Tant(IFET(IFGT(pasX,PasX),IFGT(pasY,PasY)))
Bblock
Test(IFOU(IFNE(pasX,DOUB(MOIT(pasX))),IFNE(pasY,DOUB(MOIT(pasY)))))
Bblock
PRINT_ATTENTION("les pas en 'X' et 'Y' doivent etre des puissances de 2");
Eblock
ATes
Bblock
Eblock
ETes
SET_ECHANTILLONNAGE(MOIT(pasX),MOIT(pasY));
/* Reduction dichotomique du maillage... */
CALS(Imove(ImageA,ImageR));
/* Deplacement pour permettre l'iteration (par 'Tant'). */
CALS(Iconvolution(ImageR
,facteur_multiplicatif
,ImageA
,niveaux_a_traiter,niveaux_cumulables
,nombre_de_points,noyau,inhibition_du_noyau
)
);
/* Et convolution de l'image Argument. */
Eblock
ETan
CALi(Iupdate_image(nom_imageR,ImageR));
Eblock
ATes
Bblock
Test__CODE_ERREUR__ERREUR07;
Eblock
ETes
RETU_Commande;
Eblock
ECommande