/*************************************************************************************************************************************/
/* */
/* E R O S I O N O U D I L A T A T I O N D ' U N E I M A G E : */
/* */
/* */
/* Definition : */
/* */
/* Cette commande genere une image */
/* dont le nom est le premier argument */
/* d'appel ; elle est le resultat de */
/* la dilation ou de l'erosion de l'image */
/* argument. */
/* */
/* */
/* Author of '$xci/dilate.01$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 1994??????????). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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
/* Introduit le 20101019135619... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* P A R A M E T R E S : */
/* */
/*************************************************************************************************************************************/
#define DILATER_L_IMAGE \
VRAI \
/* Precise s'il faut dilater l'image ('VRAI') ou bien l'eroder ('FAUX')... */
#define ATTENUER_LA_DILATATION \
FAUX
#define FACTEUR_D_ATTENUATION_DE_LA_DILATATION \
FU
/* Lorsque 'IL_FAUT(dilater_l_image)', precise s'il faut attnuer la dilatation ('VRAI') ou */
/* bien la conserver tel quel ('FAUX'). Ceci fut introduit le 20211024101937... */
#define NOMBRE_D_ITERATIONS \
UN \
/* Nombre d'iterations de l'erosion ou de la dilatation. Ceci fut introduit le */ \
/* 20190715151321 apres que j'ai decouvert que ce processus avit un defaut lorsqu'il */ \
/* y avait en quelque sorte des "canaux" etroits entre deux zones a traiter (cela s'est */ \
/* vu avec 'v $xiirv/IMPO.51.1')... */
#define TRAITER_TOUS_LES_NIVEAUX \
VRAI \
/* Indique s'il faut traiter tous les niveaux ('VRAI') ou uniquement le 'NOIR' ('FAUX'). */
#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 COMPATIBILITE_20230117 \
FAUX \
/* Permet de proceder a des calculs compatibles a ceux qui furent effectues anterieurement */ \
/* au 20230117184209. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* M A C R O S U T I L E S : */
/* */
/*************************************************************************************************************************************/
#include xci/convol.01.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* E R O S I O N O U D I L A T A T I O N D ' U N E I M 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(dilater_l_image,DILATER_L_IMAGE));
/* Precise s'il faut dilater l'image ('VRAI') ou bien l'eroder ('FAUX')... */
DEFV(Logical,INIT(attenuer_la_dilatation,ATTENUER_LA_DILATATION));
DEFV(Float,INIT(facteur_d_attenuation_de_la_dilatation,FACTEUR_D_ATTENUATION_DE_LA_DILATATION));
/* Lorsque 'IL_FAUT(dilater_l_image)', precise s'il faut attnuer la dilatation ('VRAI') ou */
/* bien la conserver tel quel ('FAUX'). Ceci fut introduit le 20211024101937... */
DEFV(Int,INIT(nombre_d_iterations,NOMBRE_D_ITERATIONS));
/* Nombre d'iterations de l'erosion ou de la dilatation. Ceci fut introduit le */
/* 20190715151321 apres que j'ai decouvert que ce processus avit un defaut lorsqu'il */
/* y avait en quelque sorte des "canaux" etroits entre deux zones a traiter (cela s'est */
/* vu avec 'v $xiirv/IMPO.51.1')... */
DEFV(Logical,INIT(traiter_tous_les_niveaux,TRAITER_TOUS_LES_NIVEAUX));
/* Indique s'il faut traiter tous les niveaux ('VRAI') ou uniquement le 'NOIR' ('FAUX'). */
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(Logical,DTb1(inhibition_du_noyau,TAILLE_MAXIMALE_D_UN_NOYAU_DE_CONVOLUTION));
/* Et sa liste d'activite. */
DEFV(Logical,INIT(compatibilite_20230117,COMPATIBILITE_20230117));
/* Permet de proceder a des calculs compatibles a ceux qui furent effectues anterieurement */
/* au 20230117184209. */
/*..............................................................................................................................*/
GET_ARGUMENTSi(nombre_d_arguments
,BLOC(GET_ARGUMENT_L("compatibilite_20230117=",compatibilite_20230117);
/* Options introduites le 20230117184209... */
GET_ARGUMENT_L("Pconvolution_____compatibilite_2014082110=""compatibilite_2014082110="
,Pconvolution_____compatibilite_2014082110
);
/* Options introduites le 20140821170005... */
GET_ARGUMENT_L("Pconvolution_____compatibilite_20141124=""compatibilite_20141124="
,Pconvolution_____compatibilite_20141124
);
/* Options introduites le 20141124185450... */
GET_ARGUMENT_L
("Fparcours_circulaire_d_une_spirale_carree_____compatibilite_2014082117=""compatibilite_2014082117="
,Fparcours_circulaire_d_une_spirale_carree_____compatibilite_2014082117
);
/* Options introduites le 20140821175732... */
GET_ARGUMENTS2_L(chain_Aconcaten3("Ierosion_ou_dilatation_____compatibilite_20140901="
,"IFerosion_ou_dilatation_____compatibilite_20140901="
,"compatibilite_20140901="
)
,Ierosion_ou_dilatation_____compatibilite_20140901
,IFerosion_ou_dilatation_____compatibilite_20140901
);
/* Options introduites le 20140901084800... */
GET_ARGUMENT_C("imageA=""A=",nom_imageA);
GET_ARGUMENT_C("imageR=""R=",nom_imageR);
GET_ARGUMENT_L("standard=",les_images_sont_standards);
/* Option introduite le 20101019135619... */
GET_ARGUMENT_L("dilater=""dilatation=",dilater_l_image);
/* Parametre introduit le 20051007151412... */
GET_ARGUMENT_N("eroder=""erosion=",dilater_l_image);
/* Parametre introduit le 20050921144609 et complete le 20051007151412... */
GET_ARGUMENT_L("attenuer_dilatation=""attenuer=""ad=",attenuer_la_dilatation);
GET_ARGUMENT_F("facteur_attenuation_dilatation=""fad=",facteur_d_attenuation_de_la_dilatation);
/* Parametres introduits le 20211024101937... */
GET_ARGUMENT_I("iterations=""iter=",nombre_d_iterations);
/* Option introduite le 20190715151321... */
GET_ARGUMENT_L("tous_niveaux=""tous=",traiter_tous_les_niveaux);
GET_ARGUMENT_N("NOIR_uniquement=""NOIR=",traiter_tous_les_niveaux);
/* Le 20090130183703, 'GET_ARGUMENT_N(...,traiter_tous_les_niveaux)' fut introduit... */
/* */
/* Le 20171024125656 je note que 'traiter_tous_les_niveaux=FAUX' n'est pas tres */
/* interessant comme cela s'est vu lors de la mise au point de 'v $xiirv/HYCO.11$R04'. */
/* En effet, soit une image toute NOIR avec quelques points BLANC que l'on cherche a */
/* agrandir. On pourrait esperer que : */
/* */
/* NOIR=VRAI */
/* */
/* optimiserait le processus. Il n'en est rien car, en effet, ce sont bien les points qui */
/* sont NOIR et dans le voisinnage de points BLANC qui deviennent BLANC afin de faire */
/* l'agrandissment attendu. Ainsi, "NOIR=VRAI" ne signifie pas ne traiter que les points */
/* qui ne sont pas NOIR, mais exactement le CONTRAIRE ! */
GET_ARGUMENT_L("noyau_carre=""ncarre=",le_noyau_de_convolution_est_carre);
/* Le 20060118092621, "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 20071003113832... */
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);
/* Option introduite le 20040703110138... */
GET_ARGUMENT_L("circulaire_optimale="
,Pconvolution_____adapter_le_nombre_de_points_lors_du_parcours_circulaire_des_spirales_carrees
);
/* Option introduite le 20061224122235... */
GET_ARGUMENT_I("points=""p=",nombre_de_points);
/* On rappelle le 20161031192719 que faire : */
/* */
/* points=1 */
/* */
/* rend ce programme neutre qu'il soit en mode dilatation ou bien en mode erosion... */
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 20191205114121... */
GET_ARGUMENT_F("ponderation=",Ierosion_ou_dilatation_____ponderation_de_dilatation);
GET_ARGUMENTS2_L("forcer_moyenne=""moyenne=""fm="
,Ierosion_ou_dilatation_____forcer_le_marquage_de_la_moyenne_dans_tous_les_cas
,IFerosion_ou_dilatation_____forcer_le_marquage_de_la_moyenne_dans_tous_les_cas
)
/* Options introduites le 20140831112842. Elles permettent de visualiser la moyenne telle */
/* qu'elle sera utilisee lors d'une erosion et donc de faciliter la determination du seuil */
/* '...seuil_de_la_moyenne_sur_la_spirale_lors_d_une_erosion' eventuel... */
GET_ARGUMENTS2_F("seuil_moyenne=""seuil=""sm="
,Ierosion_ou_dilatation_____seuil_de_la_moyenne_sur_la_spirale_lors_d_une_erosion
,IFerosion_ou_dilatation_____seuil_de_la_moyenne_sur_la_spirale_lors_d_une_erosion
)
/* Options introduites le 20140830180329... */
GET_ARGUMENTS2_N("minimum="
,Ierosion_ou_dilatation_____conserver_le_maximum_local_lors_d_une_erosion
,IFerosion_ou_dilatation_____conserver_le_maximum_local_lors_d_une_erosion
)
GET_ARGUMENTS2_L("maximum="
,Ierosion_ou_dilatation_____conserver_le_maximum_local_lors_d_une_erosion
,IFerosion_ou_dilatation_____conserver_le_maximum_local_lors_d_une_erosion
)
/* Options introduites le 20051207101624 et etendues le 20101019140153... */
)
);
Test(IL_NE_FAUT_PAS(compatibilite_20230117))
Bblock
Test(IFEQ(nombre_de_points,UN))
Bblock
EGAL(Pconvolution_____parcourir_circulairement_les_spirales_carrees,FAUX);
/* Introduit le 20230117184209... */
Eblock
ATes
Bblock
Eblock
ETes
Eblock
ATes
Bblock
Eblock
ETes
Test(IL_FAUT(attenuer_la_dilatation))
Bblock
Test(IFET(EST_VRAI(les_images_sont_standards),IL_FAUT(dilater_l_image)))
Bblock
Eblock
ATes
Bblock
PRINT_ATTENTION("l'attenuation ne peut s'appliquer qu'aux images standards subissant une dilatation");
Eblock
ETes
Eblock
ATes
Bblock
Eblock
ETes
CALi(gInettoyage(les_images_sont_standards,ImageR,IFmageR));
/* Initialisation de l'image Resultat. */
BoIn(niveau,NOIR,BLANC,PAS_COULEURS)
Bblock
Test(IFOU(IL_FAUT(traiter_tous_les_niveaux)
,IFET(IL_NE_FAUT_PAS(traiter_tous_les_niveaux)
,IFEQ(niveau,NOIR)
)
)
)
Bblock
EGAL(ITb1(niveaux_a_traiter,INDX(niveau,NOIR)),VRAI);
Eblock
ATes
Bblock
EGAL(ITb1(niveaux_a_traiter,INDX(niveau,NOIR)),FAUX);
Eblock
ETes
EGAL(ITb1(niveaux_cumulables,INDX(niveau,NOIR)),VRAI);
/* Initialisation telle que tous les niveaux soient a la fois "traitables" et "cumulables" */
/* sauf si l'on doit ne traiter que le 'NOIR'... */
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
DEFV(Int,INIT(index,UNDEF));
/* Index d'initialisation du noyau. */
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. */
/* */
/* Je note le 20191205155932 que c'est a cause de la mise a 'ACTIF' de tous les points */
/* du noyau qu'il est impossible de faire de la dilatation ou de l'erosion anisotrope... */
/* */
/* De plus, le 20191205102500 j'ai supprime le forcage des neuf premiers points du */
/* noyau (comme cela est fait dans 'v $xci/convol.01$K NOYAU.ADD2.PREMIER_POINT.' par */
/* exemple) : cela etait justifie par le fait que les valeurs du noyau ne jouent aucun */
/* role dans le processus de dilatation ou d'erosion... */
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 20111031140910 car, en effet, manquait cruellement... */
Eblock
ETes
Test(PAS_D_ERREUR(CODE_ERROR(gIload_image(les_images_sont_standards,ImageR,IFmageR,nom_imageA))))
/* Chargement de l'image a transformer. */
Bblock
Repe(nombre_d_iterations)
/* Possibilite introduite le 20190715151321... */
Bblock
CALS(gImove(les_images_sont_standards,ImageA,IFmageA,ImageR,IFmageR));
CALS(gIerosion_ou_dilatation(les_images_sont_standards
,ImageR,IFmageR
,ImageA,IFmageA
,niveaux_a_traiter,niveaux_cumulables
,nombre_de_points,noyau,inhibition_du_noyau
,dilater_l_image
)
);
/* Et dilatation (ou erosion) de l'image Argument. */
Test(EST_VRAI(les_images_sont_standards))
Bblock
Test(IL_FAUT(dilater_l_image))
Bblock
Test(IL_FAUT(attenuer_la_dilatation))
Bblock
BDEFV(image,distance_ImageA_ImageR);
BDEFV(image,distance_ImageA_ImageR_attenuee);
BDEFV(image,ImageR_attenuee);
CALS(Idistance_NOIR_des_points(distance_ImageA_ImageR,ImageA,ImageR));
CALS(Iscale(distance_ImageA_ImageR_attenuee
,facteur_d_attenuation_de_la_dilatation
,distance_ImageA_ImageR
,FZERO
)
);
CALS(Iproduit_de_masquage(ImageR_attenuee,distance_ImageA_ImageR_attenuee,ImageR));
CALS(Imove(ImageR,ImageR_attenuee));
/* Attenuation de la dilatation (introduite le 20211024101937)... */
/* */
/* On notera le 20211024112603 que : */
/* */
/* dilatation=VRAI attenuer_dilatation=VRAI facteur_attenuation_dilatation=1 */
/* */
/* n'est pas equivalent a : */
/* */
/* dilatation=VRAI attenuer_dilatation=FAUX */
/* */
/* et ce a cause apparemment du calcul de la distance... */
EDEFV(image,ImageR_attenuee);
EDEFV(image,distance_ImageA_ImageR_attenuee);
EDEFV(image,distance_ImageA_ImageR);
Eblock
ATes
Bblock
Eblock
ETes
Eblock
ATes
Bblock
Eblock
ETes
Eblock
ATes
Bblock
Eblock
ETes
Eblock
ERep
CALi(gIupdate_image(les_images_sont_standards,nom_imageR,ImageR,IFmageR));
Eblock
ATes
Bblock
Test__CODE_ERREUR__ERREUR07;
Eblock
ETes
RETU_Commande;
Eblock
ECommande