/*************************************************************************************************************************************/
/* */
/* R E C H E R C H E E T M A R Q U A G E D E V O I S I N A G E : */
/* */
/* */
/* Author of '$xci/MarqueVoi.01$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20070209101317). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 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_PREMIERS_VOISINS \
FU
#define VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION_SECONDS_VOISINS \
FZERO
/* Facteur implicite des elements du noyau de convolution. */
#define NIVEAU_INTERIEUR \
GRIS_8
#define NIVEAU_EXTERIEUR \
GRIS_0
#define NIVEAU_DE_MARQUAGE_DU_VOISINAGE \
GRIS_4
/* Si le point courant {X,Y} possede le niveau 'niveau_interieur' et si dans son voisinage */
/* il y a au moins un point de niveau 'niveau_exterieur', alors ce point courant {X,Y} est */
/* marque avec 'niveau_de_marquage_du_voisinage'... */
#define CONFONDRE_LES_IMAGES_ARGUMENT_ET_RESULTAT \
FAUX \
/* Faut-il confondre les images "non standard" Argument et Resultat ('VRAI') ou les */ \
/* differencier ('FAUX'). */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* M A C R O S U T I L E S : */
/* */
/*************************************************************************************************************************************/
#include xci/convol.01.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* R E C H E R C H E E T M A R Q U A G E D E V O I S I N A G E : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
DEFV(CHAR,INIC(POINTERc(nom_imageR),NOM_PIPE));
DEFV(CHAR,INIC(POINTERc(nom_imageA),NOM_PIPE));
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(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(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_SECONDS_VOISINS));
DEFV(Float,INIT(noyau_02,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION_PREMIERS_VOISINS));
DEFV(Float,INIT(noyau_03,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION_SECONDS_VOISINS));
DEFV(Float,INIT(noyau_04,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION_PREMIERS_VOISINS));
DEFV(Float,INIT(noyau_05,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION_SECONDS_VOISINS));
DEFV(Float,INIT(noyau_06,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION_PREMIERS_VOISINS));
DEFV(Float,INIT(noyau_07,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION_SECONDS_VOISINS));
DEFV(Float,INIT(noyau_08,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION_PREMIERS_VOISINS));
DEFV(Float,INIT(noyau_09,VALEUR_INITIALE_DU_NOYAU_DE_CONVOLUTION_SECONDS_VOISINS));
/* Definition des 9 premiers elements du noyau (meme si le noyau fait moins de 9 points), */
/* et qui sont redefinissables et initialise sur la 4-connexite... */
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(genere_p,INIT(niveau_interieur,NIVEAU_INTERIEUR));
DEFV(genere_p,INIT(niveau_exterieur,NIVEAU_EXTERIEUR));
DEFV(genere_p,INIT(niveau_de_marquage_du_voisinage,NIVEAU_DE_MARQUAGE_DU_VOISINAGE));
/* Si le point courant {X,Y} possede le niveau 'niveau_interieur' et si dans son voisinage */
/* il y a au moins un point de niveau 'niveau_exterieur', alors ce point courant {X,Y} est */
/* marque avec 'niveau_de_marquage_du_voisinage'... */
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'). */
/*..............................................................................................................................*/
GET_ARGUMENTSi(nombre_d_arguments
,BLOC(GET_ARGUMENT_C("imageA=""A=",nom_imageA);
GET_ARGUMENT_C("imageR=""R=",nom_imageR);
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 20060118093420, "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 20071003113619... */
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("circulaire=",Pconvolution_____parcourir_circulairement_les_spirales_carrees);
GET_ARGUMENT_L("circulaire_optimale="
,Pconvolution_____adapter_le_nombre_de_points_lors_du_parcours_circulaire_des_spirales_carrees
);
/* Option introduite le 20061224114709... */
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_I("spirale_horizontal=",SPIRALE_DEFINITION_____facteur_delta_horizontal);
GET_ARGUMENT_I("spirale_vertical=",SPIRALE_DEFINITION_____facteur_delta_vertical);
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 20210309110643... */
GET_ARGUMENT_L("standard=",les_images_sont_standards);
GET_ARGUMENT_L("confondre=",confondre_les_images_Argument_et_Resultat);
GET_ARGUMENT_P("niveau_interieur=""ni=",niveau_interieur);
GET_ARGUMENT_P("niveau_exterieur=""ne=",niveau_exterieur);
GET_ARGUMENT_P("nm=""niveau_de_marquage=",niveau_de_marquage_du_voisinage);
)
);
CALi(Inoir(ImageR));
/* Initialisation de l'image Resultat. */
BoIn(niveau,NOIR,BLANC,PAS_COULEURS)
Bblock
EGAL(ITb1(niveaux_a_traiter,INDX(niveau,NOIR)),VRAI);
EGAL(ITb1(niveaux_cumulables,INDX(niveau,NOIR)),VRAI);
/* Initialisation telle que tous les niveaux soient a la fois "traitables" et "cumulables". */
Eblock
EBoI
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,FZERO);
/* Initialisation du noyau de convolution de facon a evaluer a priori l'histogramme dans */
/* 'Irecherche_et_marquage_de_voisinage(...) nulle part. */
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 20111031140638 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);
/* 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... */
/* */
/* L'initialisation du noyau de convolution est faite de facon a n'evaluer l'histogramme */
/* dans 'Irecherche_et_marquage_de_voisinage(...) que pour les quatre premier voisins... */
Test(PAS_D_ERREUR(CODE_ERROR(Iload_image(ImageA,nom_imageA))))
/* Chargement de l'image a transformer. */
Bblock
CALS(Irecherche_et_marquage_de_voisinage(ImageR
,ImageA
,niveaux_a_traiter,niveaux_cumulables
,nombre_de_points,noyau,inhibition_du_noyau
,niveau_interieur
,niveau_exterieur
,niveau_de_marquage_du_voisinage
)
);
/* Et marquage des voisinages... */
/* */
/* Si le point courant {X,Y} possede le niveau 'niveau_interieur' et si dans son voisinage */
/* il y a au moins un point de niveau 'niveau_exterieur', alors ce point courant {X,Y} est */
/* marque avec 'niveau_de_marquage_du_voisinage'... */
/* */
/* On notera le 20070211101705 que cela ne donne pas (ainsi que je l'esperai naivement...) */
/* un unique ensemble de points 4-connexes, mais en fait evidemment un ensemble d'ensembles */
/* de points 4-connexes. Ce programme ne peut donc pas etre utilise pour extraire, par */
/* exemple, un front de diffusion (4-connexe)... */
CALi(Iupdate_image(nom_imageR,ImageR));
Eblock
ATes
Bblock
Test__CODE_ERREUR__ERREUR07;
Eblock
ETes
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.