/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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  :                                                           */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        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.                                                                             */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xci/convol.01$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
#include  image_image_IMAGESF_EXT

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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 20060304184438 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 20060302163815...                                           */

#define   LE_NOYAU_DE_CONVOLUTION_EST_CARRE                                                                                             \
                    VRAI                                                                                                                \
                                        /* Indique si le noyau de convolution est carre ('VRAI') ou circulaire ('FAUX').             */
#define   NOMBRE_DE_POINTS_DU_NOYAU                                                                                                     \
                    EXP2(DOUP(UN))                                                                                                      \
                                        /* Nombre de points du noyau.                                                                */
#define   VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION                                                                                       \
                    FU                                                                                                                  \
                                        /* Facteur implicite des elements du noyau de convolution.                                   */
#define   FACTEUR_MULTIPLICATIF_DE_CONVOLUTION                                                                                          \
                    FU                                                                                                                  \
                                        /* Facteur multiplicatif du produit de convolution en chaque point {X,Y}.                    */

#define   IL_FAUT_RENORMALISER                                                                                                          \
                    FAUX                                                                                                                \
                                        /* Faut-il renormaliser l'image ?                                                            */

#define   CONFONDRE_LES_IMAGES_ARGUMENT_ET_RESULTAT                                                                                     \
                    FAUX                                                                                                                \
                                        /* Faut-il confondre les images "non standard" Argument et Resultat ('VRAI') ou les          */ \
                                        /* differencier ('FAUX'). Cette possibilite a ete introduite le 20000728165555.              */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        M A C R O S   U T I L E S  :                                                                                               */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
#include  xci/convol.01.I"

#define   EFFECTIVE_IFmageA                                                                                                             \
                    COND(IL_NE_FAUT_PAS(confondre_les_images_Argument_et_Resultat),IFmageA,IFmageR)                                     \
                                        /* Afin de permettre de confondre les images Argument et Resultat non standard...            */

/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        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  :                                                           */
/*                                                                                                                                   */
/*************************************************************************************************************************************/
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 20060304184438...                                                           */
     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 20060302163815...                                           */

     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(Logical,INIT(le_noyau_de_convolution_est_carre,LE_NOYAU_DE_CONVOLUTION_EST_CARRE));
                                        /* Indique si le noyau de convolution est carre ('VRAI') ou circulaire ('FAUX').             */
     DEFV(Float,INIT(facteur_multiplicatif,FACTEUR_MULTIPLICATIF_DE_CONVOLUTION));
                                        /* 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(Float,INIT(noyau_01,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_02,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_03,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_04,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_05,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_06,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_07,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_08,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_09,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_10,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_11,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_12,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_13,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_14,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_15,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_16,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_17,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_18,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_19,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_20,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_21,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_22,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_23,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_24,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
     DEFV(Float,INIT(noyau_25,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION));
                                        /* Definition des 9 premiers elements du noyau (meme si le noyau fait moins de 9 points),    */
                                        /* et qui sont redefinissables. A titre d'exemple :                                          */
                                        /*                                                                                           */
                                        /*   p=9 n1=0 n2=+1 n3=0 n4=0 n5=0 n6=-1 n7=0 n8=0 n9=0 normaliser=FAUX renormaliser=VRAI    */
                                        /*                                                                                           */
                                        /* donne la derivee de 'ImageR' par rapport a 'OX', alors que :                              */
                                        /*                                                                                           */
                                        /*   p=9 n1=0 n2=0 n3=0 n4=+1 n5=0 n6=0 n7=0 n8=-1 n9=0 normaliser=FAUX renormaliser=VRAI    */
                                        /*                                                                                           */
                                        /* donne la derivee de 'ImageR' par rapport a 'OY'.                                          */
                                        /*                                                                                           */
                                        /* Le passage de 9 a 25 a ete introduit le 20210307120745...                                 */
     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.                                                          */

     DEFV(Logical,INIT(confondre_les_images_Argument_et_Resultat,CONFONDRE_LES_IMAGES_ARGUMENT_ET_RESULTAT));
                                        /* Faut-il confondre les images "non standard" Argument et Resultat ('VRAI') ou les          */
                                        /* differencier ('FAUX'). Cette possibilite a ete introduite le 20000728165555.              */
     /*..............................................................................................................................*/
     EGAL(les_images_standards_sont_a_renormaliser,IL_FAUT_RENORMALISER);
                                        /* Faut-il renormaliser l'image ?                                                            */

     GET_ARGUMENTSi(nombre_d_arguments
                   ,BLOC(GET_ARGUMENT_L("Pconvolution_____compatibilite_2014082110=""compatibilite_2014082110="
                                       ,Pconvolution_____compatibilite_2014082110
                                        );
                                        /* Options introduites le 20140821165931...                                                  */
                         GET_ARGUMENT_L("Pconvolution_____compatibilite_20141124=""compatibilite_20141124="
                                       ,Pconvolution_____compatibilite_20141124
                                        );
                                        /* Options introduites le 20141124185423...                                                  */

                         GET_ARGUMENT_L
                             ("Fparcours_circulaire_d_une_spirale_carree_____compatibilite_2014082117=""compatibilite_2014082117="
                             ,Fparcours_circulaire_d_une_spirale_carree_____compatibilite_2014082117
                              );
                                        /* Options introduites le 20140821175813...                                                  */

                         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 20040703105913...                                                    */
                         GET_ARGUMENT_L("circulaire_optimale="
                                       ,Pconvolution_____adapter_le_nombre_de_points_lors_du_parcours_circulaire_des_spirales_carrees
                                        );
                                        /* Option introduite le 20061224114709...                                                    */
                                        /*                                                                                           */
                                        /* Je note le 20190208144516 qu'utiliser :                                                   */
                                        /*                                                                                           */
                                        /*                  circulaire=VRAI circulaire_optimale=FAUX                                 */
                                        /*                                                                                           */
                                        /* provoque le message 'v $xiipf/fonction.3$FON point.courant.est.hors.limite..2.' pour      */
                                        /* tous les nombre de "points=", sauf pour :                                                 */
                                        /*                                                                                           */
                                        /*                  {2,3,4,5}                                                                */
                                        /*                  {10,11,12,13,26,27,28,29}                                                */
                                        /*                                                                                           */
                                        /* et je n'ai pas trouve d'autres valeurs inferieures a 100.                                 */
                                        /*                                                                                           */
                                        /* C'est cela qui a provoque la modification 'v $xiii/di_image$FON 20190208144946'...        */

                         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_L("noyau_carre=""ncarre=",le_noyau_de_convolution_est_carre);
                                        /* Le 20060118092454, "carre=" a ete remplace par "noyau_carre=" et "ncarre=" (risque de     */
                                        /* double definition...).                                                                    */
                         GET_ARGUMENT_N("noyau_circulaire=""ncirculaire=",le_noyau_de_convolution_est_carre);
                                        /* Ces deux parametres furent introduits le 20071003113747...                                */
                         GET_ARGUMENT_F("facteur=""f=",facteur_multiplicatif);
                         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_L("normaliser=",Pconvolution_____normaliser_le_cumul_pondere);
                         GET_ARGUMENT_F("transfert=",Pconvolution_____exposant_de_la_fonction_de_transfert);
                         GET_ARGUMENT_I("points=""p=",nombre_de_points);
                         GET_ARGUMENT_F("n1=",noyau_01);
                         GET_ARGUMENT_F("n2=",noyau_02);
                         GET_ARGUMENT_F("n3=",noyau_03);
                         GET_ARGUMENT_F("n4=",noyau_04);
                         GET_ARGUMENT_F("n5=",noyau_05);
                         GET_ARGUMENT_F("n6=",noyau_06);
                         GET_ARGUMENT_F("n7=",noyau_07);
                         GET_ARGUMENT_F("n8=",noyau_08);
                         GET_ARGUMENT_F("n9=",noyau_09);
                         GET_ARGUMENT_F("n10=",noyau_10);
                         GET_ARGUMENT_F("n11=",noyau_11);
                         GET_ARGUMENT_F("n12=",noyau_12);
                         GET_ARGUMENT_F("n13=",noyau_13);
                         GET_ARGUMENT_F("n14=",noyau_14);
                         GET_ARGUMENT_F("n15=",noyau_15);
                         GET_ARGUMENT_F("n16=",noyau_16);
                         GET_ARGUMENT_F("n17=",noyau_17);
                         GET_ARGUMENT_F("n18=",noyau_18);
                         GET_ARGUMENT_F("n19=",noyau_19);
                         GET_ARGUMENT_F("n20=",noyau_20);
                         GET_ARGUMENT_F("n21=",noyau_21);
                         GET_ARGUMENT_F("n22=",noyau_22);
                         GET_ARGUMENT_F("n23=",noyau_23);
                         GET_ARGUMENT_F("n24=",noyau_24);
                         GET_ARGUMENT_F("n25=",noyau_25);
                                        /* Le passage de 9 a 25 a ete introduit le 20210307120745...                                 */
                                        /*                                                                                           */
                                        /* Disposition des 25 premiers points du noyau sur la spirale :                              */
                                        /*                                                                                           */
                                        /*                                                                                           */
                                        /*                  17---16---15---14---13                                                   */
                                        /*                   |                   |                                                   */
                                        /*                   |                   |                                                   */
                                        /*                  18    5----4----3   12                                                   */
                                        /*                   |    |         |    |                                                   */
                                        /*                   |    |         |    |                                                   */
                                        /*                  19    6    1----2   11                                                   */
                                        /*                   |    |              |                                                   */
                                        /*                   |    |              |                                                   */
                                        /*                  20    7----8----9---10                                                   */
                                        /*                   |                                                                       */
                                        /*                   |                                                                       */
                                        /*                  21---22---23---24---25                                                   */
                                        /*                                                                                           */

                         GET_ARGUMENT_L("editer_noyau_convolution=""enc=",Pconvolution_____editer_le_noyau_de_convolution);
                         GET_ARGUMENT_I("Xedition_noyau_convolution=""Xenc="
                                       ,Pconvolution_____coordonnee_X_d_edition_du_noyau_de_convolution
                                        );
                         GET_ARGUMENT_I("Yedition_noyau_convolution=""Yenc="
                                       ,Pconvolution_____coordonnee_Y_d_edition_du_noyau_de_convolution
                                        );
                                        /* Possibilite introduite le 20210309110450...                                               */

                         GET_ARGUMENT_I("spirale_horizontal=",SPIRALE_DEFINITION_____facteur_delta_horizontal);
                         GET_ARGUMENT_I("spirale_vertical=",SPIRALE_DEFINITION_____facteur_delta_vertical);
                         GET_ARGUMENT_F("Pcumul=""Pcum=",Pconvolution_____ponderation_de_cumul_courant_normalise_apres_transfert);
                         GET_ARGUMENT_F("Pminimum=""Pmin=",Pconvolution_____ponderation_de_Pconvolution_____minimum_sur_la_spirale);
                         GET_ARGUMENT_F("Pmaximum=""Pmax=",Pconvolution_____ponderation_de_Pconvolution_____maximum_sur_la_spirale);
                         GET_ARGUMENT_F("Pmoyenne=""Pmoy=",Pconvolution_____ponderation_de_Pconvolution_____moyenne_sur_la_spirale);
                                        /* Options {"Pcumul=","Pminimum=","Pmaximum="} introduites le 20070206123053...              */
                                        /*                                                                                           */
                                        /* L'option "Pmoyenne=" a ete introduite le 20161117134006, meme si cela ne sert pas a grand */
                                        /* chose, la moyenne ressemblant beaucoup au cumul si ce n'est que ce dernier est pondere... */
                                        /*                                                                                           */
                                        /* L'option "Pcum=" a ete introduite le 20161117135256 par "symetrie" avec les suivantes...  */

                         GET_ARGUMENT_L("conserver_minimum=""cmin="
                                       ,Pconvolution_____conserver_le_niveau_du_centre_de_la_spirale_s_il_est_le_minimum
                                        );
                         GET_ARGUMENT_L("conserver_maximum=""cmax="
                                       ,Pconvolution_____conserver_le_niveau_du_centre_de_la_spirale_s_il_est_le_maximum
                                        );
                                        /* Options {"conserver_minimum=","conserver_maximum="} introduites le 20130630185044...      */

                         GET_ARGUMENT_L("disperser=",Pconvolution_____disperser_le_noyau_de_convolution);

                         GET_ARGUMENT_F("disperser_X_a00=""dX00=",Pconvolution_____disperser_le_noyau_de_convolution_X_a00);
                         GET_ARGUMENT_F("disperser_X_a01=""dX01=",Pconvolution_____disperser_le_noyau_de_convolution_X_a01);
                         GET_ARGUMENT_F("disperser_X_a02=""dX02=",Pconvolution_____disperser_le_noyau_de_convolution_X_a02);
                         GET_ARGUMENT_F("disperser_X_a10=""dX10=",Pconvolution_____disperser_le_noyau_de_convolution_X_a10);
                         GET_ARGUMENT_F("disperser_X_a11=""dX11=",Pconvolution_____disperser_le_noyau_de_convolution_X_a11);
                         GET_ARGUMENT_F("disperser_X_a12=""dX12=",Pconvolution_____disperser_le_noyau_de_convolution_X_a12);
                         GET_ARGUMENT_F("disperser_X_a20=""dX20=",Pconvolution_____disperser_le_noyau_de_convolution_X_a20);
                         GET_ARGUMENT_F("disperser_X_a21=""dX21=",Pconvolution_____disperser_le_noyau_de_convolution_X_a21);
                         GET_ARGUMENT_F("disperser_X_a22=""dX22=",Pconvolution_____disperser_le_noyau_de_convolution_X_a22);

                         GET_ARGUMENT_F("disperser_Y_a00=""dY00=",Pconvolution_____disperser_le_noyau_de_convolution_Y_a00);
                         GET_ARGUMENT_F("disperser_Y_a01=""dY01=",Pconvolution_____disperser_le_noyau_de_convolution_Y_a01);
                         GET_ARGUMENT_F("disperser_Y_a02=""dY02=",Pconvolution_____disperser_le_noyau_de_convolution_Y_a02);
                         GET_ARGUMENT_F("disperser_Y_a10=""dY10=",Pconvolution_____disperser_le_noyau_de_convolution_Y_a10);
                         GET_ARGUMENT_F("disperser_Y_a11=""dY11=",Pconvolution_____disperser_le_noyau_de_convolution_Y_a11);
                         GET_ARGUMENT_F("disperser_Y_a12=""dY12=",Pconvolution_____disperser_le_noyau_de_convolution_Y_a12);
                         GET_ARGUMENT_F("disperser_Y_a20=""dY20=",Pconvolution_____disperser_le_noyau_de_convolution_Y_a20);
                         GET_ARGUMENT_F("disperser_Y_a21=""dY21=",Pconvolution_____disperser_le_noyau_de_convolution_Y_a21);
                         GET_ARGUMENT_F("disperser_Y_a22=""dY22=",Pconvolution_____disperser_le_noyau_de_convolution_Y_a22);
                                        /* Les options de "dispersion" ont ete introduites le 20201002112707...                      */

                         GET_ARGUMENT_L("rn=""renormaliser=""r=",les_images_standards_sont_a_renormaliser);
                         GET_ARGUMENT_L("standard=",les_images_sont_standards);
                         GET_ARGUMENT_L("confondre=",confondre_les_images_Argument_et_Resultat);
                         )
                    );

     Test(IFET(EST_VRAI(les_images_sont_standards),IL_FAUT(confondre_les_images_Argument_et_Resultat)))
          Bblock
          PRINT_ATTENTION("confondre les images Argument et Resultat ne fonctionne qu'en mode non standard des images");
                                        /* En effet, 'v $xiii/di_image$DEF gIconvolution' fait appel a la fonction                   */
                                        /* 'v $xiii/di_image$FON Iconvolution_avec_renormalisation' qui passe par une 'imageF'       */
                                        /* intermediaire pour pouvoir appeler 'v $xiii/di_image$FON IFconvolution'. Alors la         */
                                        /* "confusion" de 'ImageA' et de 'ImageR' est "invisible"...                                 */
          Eblock
     ATes
          Bblock
          Eblock
     ETes

     CALi(gInettoyage(les_images_sont_standards,ImageR,IFmageR));
                                        /* Initialisation de l'image Resultat.                                                       */

     SELECTION_DES_NIVEAUX(nom_paletteA,niveaux_a_traiter,niveaux_cumulables);
                                        /* Introduit sous cette forme le 20060304184438...                                           */

     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_DE_CONVOLUTION);
                                        /* Initialisation du noyau de convolution.                                                   */
               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 20111031140801 car, en effet, manquait cruellement...                        */
          Eblock
     ETes

     NOYAU(ADD2(PREMIER_POINT,ZERO),noyau_01);
     NOYAU(ADD2(PREMIER_POINT,UN),noyau_02);
     NOYAU(ADD2(PREMIER_POINT,DEUX),noyau_03);
     NOYAU(ADD2(PREMIER_POINT,TROIS),noyau_04);
     NOYAU(ADD2(PREMIER_POINT,QUATRE),noyau_05);
     NOYAU(ADD2(PREMIER_POINT,CINQ),noyau_06);
     NOYAU(ADD2(PREMIER_POINT,SIX),noyau_07);
     NOYAU(ADD2(PREMIER_POINT,SEPT),noyau_08);
     NOYAU(ADD2(PREMIER_POINT,HUIT),noyau_09);
     NOYAU(ADD2(PREMIER_POINT,NEUF),noyau_10);
     NOYAU(ADD2(PREMIER_POINT,DIX),noyau_11);
     NOYAU(ADD2(PREMIER_POINT,ONZE),noyau_12);
     NOYAU(ADD2(PREMIER_POINT,DOUZE),noyau_13);
     NOYAU(ADD2(PREMIER_POINT,TREIZE),noyau_14);
     NOYAU(ADD2(PREMIER_POINT,QUATORZE),noyau_15);
     NOYAU(ADD2(PREMIER_POINT,QUINZE),noyau_16);
     NOYAU(ADD2(PREMIER_POINT,SEIZE),noyau_17);
     NOYAU(ADD2(PREMIER_POINT,DIX_SEPT),noyau_18);
     NOYAU(ADD2(PREMIER_POINT,DIX_HUIT),noyau_19);
     NOYAU(ADD2(PREMIER_POINT,DIX_NEUF),noyau_20);
     NOYAU(ADD2(PREMIER_POINT,VINGT),noyau_21);
     NOYAU(ADD2(PREMIER_POINT,VINGT_ET_UN),noyau_22);
     NOYAU(ADD2(PREMIER_POINT,VINGT_DEUX),noyau_23);
     NOYAU(ADD2(PREMIER_POINT,VINGT_TROIS),noyau_24);
     NOYAU(ADD2(PREMIER_POINT,VINGT_QUATRE),noyau_25);
                                        /* Initialisation des 9 premiers elements du noyau de convolution (meme si la taille         */
                                        /* demandee est inferieure a 9) avec des valeurs eventuellement arguments d'appel...         */
                                        /* On notera que l'on ecrase ainsi une partie de l'initialisation faite ci-dessus dans la    */
                                        /* boucle 'DoIn(...)' que l'on ne modifie pas pour simplifier...                             */
                                        /*                                                                                           */
                                        /* Le passage de 9 a 25 a ete introduit le 20210307120745...                                 */

     Test(PAS_D_ERREUR(CODE_ERROR(gIload_image(les_images_sont_standards,ImageA,EFFECTIVE_IFmageA,nom_imageA))))
                                        /* Chargement de l'image a transformer.                                                      */
          Bblock
          CALS(gIconvolution(les_images_sont_standards
                            ,ImageR,IFmageR
                            ,facteur_multiplicatif
                            ,ImageA,EFFECTIVE_IFmageA
                            ,niveaux_a_traiter,niveaux_cumulables
                            ,nombre_de_points,noyau,inhibition_du_noyau
                            ,les_images_standards_sont_a_renormaliser
                             )
               );
                                        /* Et convolution de l'image Argument.                                                       */

          CALi(gIupdate_image(les_images_sont_standards,nom_imageR,ImageR,IFmageR));
          Eblock
     ATes
          Bblock
          Test__CODE_ERREUR__ERREUR07;
          Eblock
     ETes

     RETU_Commande;
     Eblock
ECommande



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