, Bien qu'encourageants ces résultats sontà considérer avec prudence, car ces travaux sont récents et en cours de validation

, Nous avons eu l'occasion d'appliquer la méthode de détection de ZCAsà des données dans un contexte d'épidémiologie, où il ne s'agissait pas de détecter des changements abrupts entre différentes zones, mais de détecter des clusters dûsà l'émergence d'une maladie. La méthode s'est montrée efficace sur ce type de données car elle est parfaitement adaptéeà la détection de valeurs aberrantes

, La méthode de détection de ZCAs devraitêtre appliquée (et adaptée)à d'autres types de données

R. Adams, Contour mapping and differential systematics of geographic variation, Systematic Zoology, vol.19, pp.385-390, 1970.

R. Adler, On excursion sets, tube formulas and maxima of random fields, The Annals of Applied Probability, vol.10, issue.1, pp.1-74, 2000.

R. Adler, Geometry of Random Fields, 1981.
URL : https://hal.archives-ouvertes.fr/hal-01975050

D. Aldous, Probability approximations via the poisson clumping heuristic, 1989.

B. Allemand, Comparaison de deux modélisations probabilistes pour la classification de données géostatistiques : construction et analyse des algorithmes d'estimation, 2004.

M. Aronowich and R. Adler, Extrema and level crossings of ? 2 processes, Advances in Applied Probability, vol.18, pp.901-920, 1986.

M. Aronowich and R. Adler, Sample path behaviour of ? 2 surfaces at extrema, Advances in Applied Probability, vol.20, pp.719-738, 1988.

J. Azaïs and C. Delmas, Asymptotic expansions for the distribution of the maximum of gaussian random fields, Extremes, vol.5, issue.2, pp.183-214, 2002.

G. Barbujani and N. Oden, Detecting regions of abrupt change in maps of biological variables, Systematic Zoology, vol.28, pp.376-389, 1989.

Y. Benjamini and Y. Hochberg, Controlling the false discovery rate : a practical and powerful approach to multiple testing, Journal of the Royal Statistical Society B, vol.57, issue.1, pp.289-300, 1995.

J. Besag, On the statistical analysis of dirty pictures (with discussion), Journal of the Royal Statistical Society, serie B, vol.48, pp.259-302, 1986.

J. Besag, P. Green, D. Higdon, and M. K. , Bayesian computation and stochastic systems (with discussion), Statistical Science, vol.10, pp.3-66, 1995.

J. P. Bocquet-appel and J. Bacro, Generalized wombling, Systematic Biology, vol.43, issue.3, pp.442-448, 1994.

J. Cao, Excursion sets of random fields with application to human brain mapping, 1997.

J. Cao, The size of the connected components of excursion sets of ? 2 , t and F fields, Advances in Applied Probability, vol.31, pp.579-595, 1999.

A. Castrignanò, M. Maiorana, M. Pisante, A. V. Vonella, and G. Fecondo, A methodology to define soil units in support to precision agriculture, 2003.

P. Chaudhuri and J. Marron, SiZer for exploration of structures in curves, Journal of the American Statistical Association, 1999.

P. Chaudhuri and J. Marron, Scale space view of curve estimation, The Annals of Statistics, vol.28, issue.2, pp.408-428, 2000.

J. Chilès and P. Delfiner, Geostatistics : modeling spatial uncertainty, 1999.

, Precision agriculture : spatial and temporal variability of environmental quality, Ciba Foundation Symposium, vol.210, 1997.

N. Cressie, Statistics for spatial data, 1993.

C. Delmas, Distribution du maximum d'un champ aléatoire et applications statistiques, 2001.

K. Friston, C. Frith, P. Liddle, and R. Frackowiak, Comparing functional (PET) images : the assessment of significance change, Journal of Cerebral Blood Flow and Metabolism, vol.11, pp.690-699, 1991.

K. Friston, K. Worsley, R. Frackowiak, J. Mazziotta, and A. Evans, Assessing the significance of focal activations using their spatial extent, Human Brain Mapping, vol.1, pp.214-220, 1994.

S. Geman and D. Geman, Stochastic relaxation, Gibbs distribution and the bayesian restoration of images, IEEE Trans. Pattern Analysis end Machine Intelligence, issue.6, pp.721-741, 1984.

J. Gilliot, A. Jullien, J. Michelin, and P. Huet, Essai de zonage d'un champà partir de cartes de rendement, 2004.

J. F. Gleyze, J. N. Bacro, and D. Allard, Detecting regions of abrupt change : wombling procedure and statistical significance, geoENV III : Geostatistics for Environmental Applications, 2001.
URL : https://hal.archives-ouvertes.fr/hal-02761023

F. Godtliebsen, J. Marron, and C. P. , Significance in scale space, 2001.

F. Godtliebsen, J. Marron, and S. Pizer, Spatial Clustering Modelling, chapter Significance in scale-space for clustering, 2002.

M. Guérif, N. Beaudoin, C. Durr, V. Houlès, J. Machet et al., Designing a field experiment for assessing soil and crop spatial variability and defining site-specific management strategies, Proceedings of the Third European Conference on Precision Agriculture, pp.677-682, 2001.

P. Hall and C. Rau, Local likelihood tracking of fault lines and boundaries, Journal of the Royal Statistical Society B, vol.63, issue.3, pp.569-582, 2001.

J. Hart, Central tendency in areal distributions, Economic Geographyy, issue.30, pp.48-59, 1954.

V. Houlès, Mise au point d'un outil de modulation intra-parcellaire de la fertilisation azotée du blé d'hiver basé sur la télédétection et un modèle de culture, 2004.

C. Lantuéjoul, Ergodicity and integral range, Journal of Microscopy, vol.161, pp.387-403, 1991.

A. Lawson and D. Denison, Spatial clustering modelling, 2002.

B. Mary, N. Beaudoin, J. M. Machet, and C. Bruchou, Characterization and analysis of soil variability within two agricultural fields : the case of water and mineral N profiles, Proceedings of the Third European Conference on Precision Agriculture, pp.431-436, 2001.

G. Matheron, Traité de géostatistique appliquée, Tome I. Number 14, 1962.

G. Matheron, Principles of geostatistics, Economic Geology, issue.58, pp.1246-1266, 1963.

K. Miller, Multidimensional Gaussian distribution, SIAM series in Applied Mathematics, 1964.

B. Nicoullaud, N. Beaudoin, J. Roque, A. Couturier, J. Maucorps et al., Relevés pédologiques et agriculture de précision, Journée de restitution de l'AIP, 2003.

N. Peyrard, Approximations de type champ moyen de modèles de champ de Markov pour la segmentation de données spatiales, 2001.

P. Qiu, Discontinuous regression surfaces fitting, The Annals of Statistics, vol.26, pp.2218-2245, 1998.

T. Shatar and A. Mcbratney, Subdividing a field into contiguous management zones using k-zones algorithm, Third European Conference on Precision Agriculture, vol.1, pp.115-120, 2001.

X. Shen, H. Huang, and N. Cressie, Nonparametric hypothesis testing for a spatial signal, Journal of the American Statistical Association, Theory and Method, vol.97, issue.460, pp.1122-1140, 2002.

M. Stein, Interpolation for spatial data : some theory for kriging, 1999.

Y. Wang, Change curve estimation via wavelets, Journal of the American Statistical Association, vol.93, pp.163-172, 1998.

B. Whelan and A. Mcbratney, Definition and interpretation of potential management zones in australia, Proceedings of the 11th Australian Agronomy Conference, 2003.

W. Womble, Differential systematics, Science, vol.114, pp.315-322, 1951.

K. Worsley, A. Evans, S. Marrett, and P. Neelin, A three dimensional statistical analysis for CBF activation studies in human brain, Journal of Cerebral Blood Flow and Metabolism, vol.12, pp.900-918, 1992.

K. Worsley, Local maxima and the expected Euler characteristic of excursion sets of ? 2 , F and t fields, Advances in Applied Probability, vol.26, pp.13-42, 1994.

K. Worsley, Estimating the number of peaks in a random field using the Hadwiger characteristic of excursion sets, with applications to medical images, The Annals of Statistics, vol.23, issue.2, pp.640-669, 1995.

K. Worsley, Testing for signals with unknown location and scale in a ? 2 random field, with an application to fMRI, Advances in Applied Probability, vol.33, issue.4, pp.773-793, 2001.

A. Yaglom, Correlation theory of stationary and related random functions, 1986.

J. Ye, Measuring and correcting the effects of data mining and model selection, Journal of the American Statistical Association, vol.98, pp.120-131, 1998.

, Echantillon de 100 points répartis aléatoirement, issus d'un processus : a) stationnaire, b) avec une discontinuité le long des pointillés. Les ZCAs significatives sont colorées en noir et les non significatives en gris : c). La grille utilisée est de dimension 60 × 60

, ZCAs (significatives en noir et non significatives en gris) obtenues pour différentes valeurs de 1 ? ? : a) 1 ? ? = 0.995, b) 1 ? ? = 0, vol.9994

. Histogramme, . De-la-statistique-x-=-t-?-s-t?, and . Det, /2 /? comparéeà une exponentielle d'espérance 2, dans le cas d'une grille 30 × 30 (à gauche) et 60 × 60 (à droite), et pour 1 ? ? = 0.95 (première ligne) et 1 ? ? = 0.995 (deuxième ligne)

. Histogramme, . De-la-statistique-x-=-t-?-s-t?, and . Det, /2 /? comparéè a une exponentielle d'espérance 2, E(2), pour un ensemble de 100 simulations, dans le cas d'une grille 30 × 30 (à gauche) et 60 × 60 (à droite), et pour 1?? [i,j] (première ligne), 1?? moy (deuxième ligne), 1?? G (troisième ligne) et 1 ?? (dernière ligne)

, ZCAs potentielles obtenues pour la simulation utilisée pour illustrer la méthode, en fonction du niveau 1 ? ? : 1 ? ? [i,j] , 1 ? ? moy , 1 ? ? G et 1 ?? (de gaucheà droite), et de la discrétisation : 30 × 30 (première ligne) et 60 × 60 (deuxième ligne)

, Cconfigurations de points x i , x j pour lesquels les couples {Z(x i ), Z(x j )} ne sont pas pris en compte dans le calcul du variogramme expérimental, p.108

, Application de la procédure itérative sur l'exemple utilisé pour illustrer la méthode. a) Echantillon de points considéré avec une discontinuité a = 2.5 introduite en x 1 = 0.4. Variogrammes obtenusà la première itération, b), età la convergence, c)

, avec introduction d'une valeur abérante (cercle plein). b) ZCAs (significatives en noir) obtenuesà partir de cetéchantillon

, ZCAs (significatives en noir et non significatives en gris) au niveau 1 ??, obtenues en utilisant le krigeage simple (1 ?? = 0.9994), a), ou le krigeage ordinaire (1 ?? = 0.9993), b), comme interpolateur,à partir de l'exemple utilisé pour illustrer la méthode

, ZCAs significatives

, Variogrammes et ZCAsà la première itération : a) et f ),à la convergence pour 1 ?? : b) et g), pour 1 ? ? 0 = 0.99 : c) h)

, 90 : e) et j)

, 8}, pour le krigeage simple (deux premières lignes) avec a) (a 1 , a 2 ) = (1, ?1), b) (1.5, ?1.5) et c) (2, ?2), d) (a 1 , a 2 ) = (2, 0), e) (2, ?1) et f ) (3, ?1) et pour le krigrage ordinaire (dernière ligne) avec g) (a 1 , a 2 ) = (1, ?1), h) (1.5, ?1.5) et i) (2, ?2)

, 8}, pour le krigeage simple (deux premières lignes) avec a) (a 1 , a 2 ) = (1, ?1), b) (1.5, ?1.5) et c) (2, ?2), d) (a 1 , a 2 ) = (2, 0), e) (2, ?1) et f ) (3, ?1) et pour le krigrage ordinaire (dernière ligne) avec g) (a 1 , a 2 ) = (1, ?1), h) (1.5, ?1.5) et i) (2, ?2)

. Modèle and .. .. De-la-puissance-en-x-?-d-?-r-2, 2 a) Exemple d'échantillon de 100 points répartis aléatoirement dans un carré de côté un, issus d'un champ gaussien standard de covariance exponentielle de portée b = 0.1. Une discontinuité est introduite en y = 0.5. b) Puissance moyenne en un point sur 1000échantillons du type a), le long de la discontinuité

, Echantillon de points issus d'un processus avec une discontinuité le long des pointillés ; b) : ZCAs (significatives en noir) détectées ; c) : les pointillés délimitent la zone permettant de définir les vraies ZCAs et la zone grise représente les bords que nouséliminons ; d) : ZCAs (significatives en noir) obtenues aprèsélimination des bords

, Histogrammes du cumul des coordonnées i des pixels [i, j] correspondants a des ZCAs significatives pour différents niveaux 1 ? ? (1 ?? en blanc et 0.995 en noir) pour 100 simulations. Le cumul des pixels correspondantà la discontinuité (a = 2.5) en i = 24 aété tracé de couleur différente. Les bords ontété conservés figure a) etéliminés figure b)

, Histogrammes des coordonnées i des pixels [i, j] constituant les ZCAs significatives au niveau 1 ?? pour 100 simulations, en fonction de la densité d'échantillonnage : a) n = 50, b) n = 100 et c) n = 200. Le cumul des pixels correspondantà la discontinuité (a = 2.5) en i = 24 aété tracé de couleur différente

, Effet de la discrétisation sur la détection des ZCAs au niveau 1 ?? (significatives en noir, non significatives en gris)

, Histogrammes des coordonnées i des pixels [i, j] constituant les ZCAs significatives au niveau 1?? pour 100 simulations, en fonction de la discrétisation : a) 15 × 15, b) 30 × 30 et c) 60 × 60. Le cumul des pixels correspondantà la discontinuité (a = 2.5) aété tracé de couleur différente, p.165

, Le cumul des pixels correspondantà la discontinuité (i = 24) aété tracé de couleur différente. Les bords ontété conservés figures a) et c) etéliminés figures b) et d)

. Carte-pédologique-;-nicoullaud, 176 13.2 a) Carte pédologique simplifiée ; b) Schéma d'échantillonnage, p.177, 2003.

, Points d'échantillonnage (proportionnelsà leur valeur) et interpolation de la teneur en eau par krigeage ordinaire (valeurs fortes en blanc et faibles en noir). La mêmeéchelle est utilisée pour illustrer les interpolations, p.178

, Points d'échantillonnage (proportionnelsà leur valeur) et interpolation de la teneur en azote minéral (valeurs fortes en blanc, valeurs faibles en noir)

, Deséchelles différentes sont utilisées pour illustrer les interpolations, p.180

, Interpolation des variables permanentes (valeursélevées en blanc et valeurs faibles en noir)

, Points d'échantillonnage (proportionnelsà leur valeur) de la teneur en eau en mars 2000. b) Interpolation sur une grille 62 × 98 par krigeage ordinaire. c) et d) (resp. e) et f ), g) et h), i) et j), k) et l)) Ajustement du variogramme et ZCAs détectéesà la première (resp. deuxième, troisième, quatrième, cinquième) itération. Les ZCAs significatives sont représentées en noir et les non significatives en gris

, ZCAs obtenuesà la convergence pour la teneur en eau en mars 2000, avec : a) 1 ? ? = 1 ? ? G , b) 1 ?? dans la procédure itérative, p.186

, ZCAs (significatives en noir, non significatives en gris) détectées pour la teneur en eauà la convergence

, ZCAs (significatives en noir, non significatives en gris) détectées pour les variables permanentesà la convergence

, ZCAs (significatives en noir, non significatives en gris) obtenues pour la teneur en eau en juillet 2002, avec un modèle de covariance a) exponentiel

, ZCAs (significatives en noir, non significatives en gris) obtenues pour la teneur en eau en juillet, avec différentes discrétisations, pp.62-98, 2002.

, et e)) et 155 × 245 (c) et f ))

, ZCAs (significatives en noir) obtenues pour la teneur en eau en juillet 2002 au niveau 1 ? ? G , en utilisant comme interpolateur a) le krigeage ordinaire, b) le krigeage simple

, Echantillon de points de la teneur en eau en juillet 2002 obtenus enéliminant m % des points de l'échantillon initial. a) m = 0, b) m = 10, c) m = 20, d) m = 30, e) m = 40, f ) m = 50, g) m = 60, h) m = 70. Les points gris correspondent aux pointséliminés et les points noirs aux points conservés, p.192

, ZCAs (significatives en noir, non significatives en gris) obtenuesà la convergence du cas I pour la teneur en eau en juillet 2002 enéliminant m % des points de l'échantillon initial. a) m = 0, b) m = 10, c) m = 20, d) m = 30, e) m = 40, f ) m = 50, g) m = 60, h) m = 70

, Liste des tableaux

. .. Al.-;-), 26 1.2 Configurations des courbures en fonction des valeurs propres, Configurations de test possibles proposées par Gleyze et, p.31, 2001.

. .. , Surface S et valeur critique p des ZCAs détectées (figure 6.2c), p.93

. .. Nombre-d'erreurs-commises-lors-du-test-de-n-hypothèses, 98 7.2 Valeurs de 1 ? ? [i,j] obtenues en chaque pixel de la grille d'interpolation pour un ensemble de 100 simulations

, Valeurs de 1?? pour une discrétisation 30×30 et 60×60 pour un ensemble de 100 simulations utilisant un modèle de covariance exponentiel de portée b = 0.1

, Portée b et variance ? 2 estiméesà la première itération età la convergence pour différentes valeurs de 1 ? ? 0

, x), gradient de Z * (x) : W (x) et variance de W (x) : ?(x), en fonction de l'interpolateur de Z(x

, Dimension de la matrice de covariance ? V ? en fonction de la taille de la fenêtre

, Nombre de simulations (sur 100) ayant des vraies ZCAs et/ou des autres ZCAs pour différentes valeurs de 1 ? ? (grille 60 × 60, n = 100, covariance exponentielle de portée b = 0.1)

, Nombre de simulations, N , sur 100 présentant au moins une vraie ZCA et taille moyenne,s, des ZCAs détectées, en fonction de la densité de points n et de la valeur de la discontinuité a

N. Nombre and D. Simulations, sur 100) ayant des ZCAs significatives et taille moyenne,s, des ZCAs détectées, en considérant comme modèle vrai un modèle exponentiel de portée b = 0.1, pour une grille 60×60 et n = 100. Les ZCAs ontété calculées pour cinq autres modèles, p.162

N. Nombre and D. Simulations, sur 100) ayant des ZCAs significatives au niveau 1 ?? et taille moyenne,s, des ZCAs détectées, pour différentes discrétisation. La méthode utilise les paramètres d'une fonction de covariance exponentielle de portée b = 0.1 et deséchantillons de n = 100 points, p.164

N. Nombre and D. Simulations, sur 100) ayant des ZCAs et taille moyenne, s, des ZCAs détectées pour deux interpolateurs du champ : krigeage simple et krigeage ordinaire (grille 60 × 60, n = 100, covariance exponentielle de portée b = 0.1)

, Evolution culturale et nombre de points d'échantillonnage des variables dynamiques en fonction de la date de prélèvement

. Résumé-statistique-de-la-teneur-en and . Eau,

. .. , Résumé statistique de la teneur en azote minéral (en kg/ha), p.181

. .. , Résumé statistique des variables permanentes, p.181

, variance (en mm 2 ) et niveau 1 ? ? Gé valués lors de la procédure itérative pour la teneur en eau en mars, Portée b (en m), p.185, 2000.

, Niveaux de confiance associés au test local de détection obtenus aux différentes itérations pour la teneur en eau en mars, 2000.

, Portée b (en m), variance ? 2 (en mm 2 ) et niveau local de confiance 1 ? ? G obtenusà la première itération età la convergence pour la teneur en eau, p.186

S. Surfaces, en m 2 ) et valeurs critiques, p, des ZCAs obtenues lors de la convergence pour la teneur

, Portée b (en m) et variance ? 2 (en mm 2 ) des fonctions de covariance estimées aux différentes itérations, en fonction du modèle de covariance initial

B. Portée, (en m), variance, ? 2 (en mm 2 ) et niveau 1 ? ? G obtenus lors de la procédure itérative appliquéeà la teneur en eau en juillet 2002, pour de? echantillons obtenus enéliminant m % des points de l'échantillon initial, p.193