/*************************************************************************************************************************************/
/* */
/* D I S T A N C E E N T R E D E U X I M A G E S E V A L U E E P A R L A R E C H E R C H E */
/* D E N I V E A U X I D E N T I Q U E S S U R U N E S P I R A L E : */
/* */
/* */
/* Definition : */
/* */
/* Cette commande genere une image */
/* dont le nom est le premier argument */
/* d'appel ; elle est le resultat du */
/* calcul des "distances" des points de */
/* meme niveau des deux images Arguments */
/* sur une spirale. */
/* */
/* */
/* Author of '$xci/distance.02$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 1996??????????). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* 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 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 FACTEUR_DU_NUMERO \
FZERO
#define FACTEUR_DE_LA_COORDONNEE_X \
FZERO
#define FACTEUR_DE_LA_COORDONNEE_Y \
FZERO
#define FACTEUR_DE_LA_DISTANCE_EUCLIDIENNE \
FU
#define FACTEUR_DU_MAXIMUM_DE_X_ET_DE_Y \
FZERO
#define FACTEUR_DU_NIVEAU \
FZERO
/* Facteurs multiplicatifs permettant de selectioner ce qui est reellement calcule... */
/* */
/* Le 20060125095944 fut introduit 'facteur_du_niveau' destine a permettre la realisation */
/* de pavages colores... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* M A C R O S U T I L E S : */
/* */
/*************************************************************************************************************************************/
#include xci/convol.01.I"
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D I S T A N C E E N T R E D E U X I M A G E S E V A L U E E P A R L A R E C H E R C H E */
/* D E N I V E A U X I D E N T I Q U E S S U R U N E S P I R A L E : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
DEFV(CHAR,INIC(POINTERc(nom_imageR),NOM_PIPE));
DEFV(CHAR,INIC(POINTERc(nom_imageA1),NOM_PIPE));
DEFV(CHAR,INIC(POINTERc(nom_imageA2),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(Float,INIT(facteur_multiplicatif,FACTEUR_MULTIPLICATIF_DE_CONVOLUTION));
/* Facteur multiplicatif du produit de convolution en chaque point {X,Y}. */
DEFV(Float,INIT(facteur_du_numero,FACTEUR_DU_NUMERO));
DEFV(Float,INIT(facteur_de_la_coordonnee_X,FACTEUR_DE_LA_COORDONNEE_X));
DEFV(Float,INIT(facteur_de_la_coordonnee_Y,FACTEUR_DE_LA_COORDONNEE_Y));
DEFV(Float,INIT(facteur_de_la_distance_euclidienne,FACTEUR_DE_LA_DISTANCE_EUCLIDIENNE));
DEFV(Float,INIT(facteur_du_maximum_de_X_et_de_Y,FACTEUR_DU_MAXIMUM_DE_X_ET_DE_Y));
DEFV(Float,INIT(facteur_du_niveau,FACTEUR_DU_NIVEAU));
/* Facteurs multiplicatifs permettant de selectioner ce qui est reellement calcule... */
/* */
/* Le 20060125095944 fut introduit 'facteur_du_niveau' destine a permettre la realisation */
/* de pavages colores... */
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 20210308132435... */
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. */
/*..............................................................................................................................*/
EGAL(Pconvolution_____adapter_le_nombre_de_points_lors_du_parcours_circulaire_des_spirales_carrees,FAUX);
/* Introduit le 20230407135543 ('v $xiii/di_image$FON 20230407133928')... */
GET_ARGUMENTSi(nombre_d_arguments
,BLOC(GET_ARGUMENT_C("imageA1=""A1=",nom_imageA1);
GET_ARGUMENT_C("imageA2=""A2=",nom_imageA2);
GET_ARGUMENT_C("imageR=""R=",nom_imageR);
GET_ARGUMENT_L("standard=",les_images_sont_standards);
GET_ARGUMENT_L("circulaire=",Pconvolution_____parcourir_circulairement_les_spirales_carrees);
/* L'option : */
/* */
/* GET_ARGUMENT_L("circulaire_optimale=" */
/* ,Pconvolution_____adapter_le_nombre_de_points_lors_du_parcours_... */
/* ); */
/* */
/* introduite le 20061224122329, a ete supprimee le 20230407135543 et ce a cause de */
/* 'v $xiii/di_image$FON 20230407133928'... */
GET_ARGUMENT_L("centre=",Pconvolution_____tester_le_niveau_du_centre_de_la_spirale);
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 20060118092726, "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 20071003113913... */
GET_ARGUMENT_I("rang=",Pconvolution_____rang_du_n_ieme_point_de_meme_niveau);
GET_ARGUMENT_L("remplacer=""inverser=""different="
,Pconvolution_____remplacer__meme_niveau__par__niveau_different
);
GET_ARGUMENT_N("meme=",Pconvolution_____remplacer__meme_niveau__par__niveau_different);
GET_ARGUMENT_F("facteur=""f=",facteur_multiplicatif);
GET_ARGUMENT_F("facteurN=""fN=",facteur_du_numero);
GET_ARGUMENT_F("facteurY=""fY=",facteur_de_la_coordonnee_X);
GET_ARGUMENT_F("facteurX=""fX=",facteur_de_la_coordonnee_Y);
GET_ARGUMENT_F("facteurD=""fD=",facteur_de_la_distance_euclidienne);
GET_ARGUMENT_F("facteurM=""fM=",facteur_du_maximum_de_X_et_de_Y);
GET_ARGUMENT_F("fn=""facteurn=",facteur_du_niveau);
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 20210308132435... */
/* */
/* 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 20210309110337... */
GET_ARGUMENT_I("spirale_horizontal=",SPIRALE_DEFINITION_____facteur_delta_horizontal);
GET_ARGUMENT_I("spirale_vertical=",SPIRALE_DEFINITION_____facteur_delta_vertical);
GET_ARGUMENT_I("sinus_spirale=",SPIRALE_PARCOURS_____sinus_de_l_angle_de_la_rotation);
GET_ARGUMENT_P("minimum=""min=""NOIR=",Ifloat_std_____niveau_NOIR);
GET_ARGUMENT_P("maximum=""max=""BLANC=",Ifloat_std_____niveau_BLANC);
/* Parametres introduits le 20060125111320... */
)
);
CALi(IFinitialisation(IFmageR,FZERO));
/* Initialisation de l'image Resultat. */
Test(EST_VRAI(les_images_sont_standards))
Bblock
/* Cas d'une image a convertir en "standard" : */
CALi(Inoir(ImageR));
/* Initialisation de l'image Resultat. */
Eblock
ATes
Bblock
/* Autres cas : */
Eblock
ETes
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,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");
CAL1(Prer1("(%d points etaient demandes, ",nombre_de_points));
/* Introduit le 20171003131755... */
EGAL(nombre_de_points,NOMBRE_DE_POINTS_DU_NOYAU);
CAL1(Prer1("la valeur %d par defaut est forcee)\n",nombre_de_points));
/* Introduit le 20111031141115 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 20210308132435... */
Test(PAS_D_ERREUR(CODE_ERROR(Iload_image(ImageA1,nom_imageA1))))
Bblock
/* Chargement de la premiere image Argument, */
Test(PAS_D_ERREUR(CODE_ERROR(Iload_image(ImageA2,nom_imageA2))))
Bblock
/* Chargement de la deuxieme image Argument. */
CALS(IFdistance_aux_niveaux_identiques(IFmageR
,facteur_multiplicatif
,ImageA1,ImageA2
,niveaux_a_traiter,niveaux_cumulables
,nombre_de_points,noyau,inhibition_du_noyau
,facteur_du_numero
,facteur_de_la_coordonnee_X
,facteur_de_la_coordonnee_Y
,facteur_de_la_distance_euclidienne
,facteur_du_maximum_de_X_et_de_Y
,facteur_du_niveau
)
);
/* Calcul de la distance aux niveaux identiques. On notera que {A1,A2} ne commutent pas ; */
/* en general, on prendra : */
/* */
/* A1 = une image uniforme (par exemple 'BLANC'), */
/* A2 = une image "pointillee" (par exemple un certain nombre de points */
/* sur fond 'NOIR). */
/* */
/* */
/* si 'IL_NE_FAUT_PAS(Pconvolution_____remplacer__meme_niveau__par__niveau_different)' */
/* et : */
/* */
/* A1 = une image uniforme (par exemple 'NOIR'), */
/* A2 = une image "pointillee" (par exemple un certain nombre de points */
/* sur fond 'NOIR). */
/* */
/* si 'IL_FAUT(Pconvolution_____remplacer__meme_niveau__par__niveau_different)'. */
/* */
/* Le 20060125095944 fut introduit 'facteur_du_niveau' destine a permettre la realisation */
/* de pavages colores... */
Test(EST_VRAI(les_images_sont_standards))
Bblock
/* Cas ou 'nom_imageR' est une image standard : */
CALS(Ifloat_std_avec_renormalisation(ImageR,IFmageR));
CALi(Iupdate_image(nom_imageR,ImageR));
Eblock
ATes
Bblock
/* Cas ou 'nom_imageR' n'est pas une image standard : */
CALi(IupdateF_image(nom_imageR,IFmageR));
Eblock
ETes
Eblock
ATes
Bblock
Test__CODE_ERREUR__ERREUR07;
Eblock
ETes
Eblock
ATes
Bblock
Test__CODE_ERREUR__ERREUR07;
Eblock
ETes
RETU_Commande;
Eblock
ECommande