FR 3069750 A1 20190208 FR 1757359 A 20170801 PROCEDE POUR AUGMENTER LA PRODUCTIVITE D'UNE PARCELLE AGRICOLE VIS-A-VIS D'UNE VARIETE VEGETALE D'UNE ESPECE DONNEE DISPOSITIFS POUR LA MISE EN OEUVRE DE CE PROCEDE Domaine de l'invention Le domaine de l'invention est celui de l'agriculture de précision. Plus précisément, l'invention concerne un procédé pour accroître la productivité d'une parcelle agricole (P) vis-à-vis d'une variété végétale d'une espèce donnée, en améliorant la gestion des intrants, tels que le semis par la modulation de sa densité. L'invention vise également: un procédé d'apport d'au moins un intrant sur la parcelle (P) destinée à produire une variété végétale d'une espèce donnée sur une parcelle (P), selon une carte de préconisation; un procédé de semis d'une variété végétale d'une espèce donnée sur une parcelle (P), selon une carte de préconisation de la densité de semis sur la parcelle (P); un procédé de réalisation d'une préconisation adaptée, de préférence d'une carte de préconisation, d'apport d'au moins un intrant sur la parcelle; un logiciel pour la mise en oeuvre du procédé selon l'invention ; un support de stockage d'une carte de préconisation d'apport d'un intrant, de préférence le semis par modulation de sa densité, sur (P); un support de stockage du logiciel, un dispositif d'apport de l'intrant, en particulier un semoir ; un module informatique comprenant le logiciel et destiné à équiper le dispositif d'apport de l'intrant, en particulier le semoir. Arrière-plan technologique - problème technique L'agriculture de précision consiste à mettre en oeuvre des moyens technologiques pour l'optimisation des rendements et de la gestion des intrants au sein d'une parcelle agricole, dans laquelle sont cultivées diverses variétés végétales appartenant à différentes espèces. Ces moyens technologiques ont pour objectif l'optimisation de la gestion des parcelles agricoles tant d'un point de vue agronomique pour rationaliser les apports d'intrants et contrôler ainsi les coûts de production, que d'un point de vue environnemental pour préserver les ressources naturelles. L'agriculture de précision peut permettre d'adapter les densités de semis sur une parcelle (P), d'une variété végétale dans une espèce donnée, selon les qualités agricoles du sol, la topographie, la climatologie,.., de la parcelle (P). 2 En effet, certaines zones de la parcelle (P) peuvent avoir un potentiel de production plus élevé et peuvent donc recevoir une densité de semis par unité de surface, plus importante de façon à maximiser le rendement des zones les plus favorables. En revanche, dans les zones de (P) moins favorables, la mise en oeuvre de densités de 5 semis faibles optimise un compromis entre apport d'intrants et rendements, en limitant la compétition vis-à-vis des ressources disponibles. En fait, l'agriculture de précision est une prise en compte de la variabilité intraparcellaire ainsi que du contexte agro-climatique, aux fins de gérer la conduite de la parcelle vis-à-vis de l'apport de la bonne dose d'intrants, au bon endroit, au bon 10 moment. Les paramètres qui déterminent la variabilité intra-parcellaire sont entre autres : la topographie, la pédologie, la climatologie, et l'historique agricole de la parcelle... Les indicateurs de biomasse obtenus par l'imagerie aérienne/satellite sont un moyen de mesurer la résultante de ces différents facteurs. 15 La demande de brevet US2003/0018431A1 décrit des méthodes et des systèmes de gestion de parcelles agricoles, dans lesquels sont réalisés des profils de différents paramètres tels que l'altitude, l'image satellite, la conductivité du sol, les rendements historiques de la parcelle. Ces profils sont combinés et analysés au moyen d'un logiciel 20 dénommé AGIS® de la société Delta Data System, pour générer un profil de gestion de la parcelle. Ainsi, il est par exemple créé 12 zones agronomiques pour la parcelle, chacune de ces zones agronomiques tenant compte de l'analyse des profils de paramètres subdivisés. La parcelle est découpée en portions, dont chacune est catégorisée dans l'une des 12 zones agronomiques. Par exemple, la zone agronomique 3 25 peut-être corrélée à un profil donné d'altitude (dépression) et à un profil donné d'images satellites (faible rendement), tandis que la zone agronomique 9 peut-être, quant à elle, corrélée à un autre profil d'images satellites (rendement moyen) et à un profil de conductivité du sol donné (faible conductivité). Selon cette demande de brevet U52003/0018431A1, le profil de gestion de la parcelle 30 peut être utilisé pour générer une carte de préconisation d'apport d'un intrant sur ladite parcelle. Au-delà de ces généralités, la demande de brevet U52003/0018431A1 ne donne aucune information précise permettant la réalisation du profil de gestion de parcelle et de la carte de préconisation d'apport d'un intrant sur ladite parcelle. En outre, elle ne propose pas la mise en oeuvre d'un modèle de culture décrivant la croissance et 35 le développement de cultures en interaction avec leurs conditions agroenvironnementales. 3 La demande de brevet US2003/0018431A1 ne divulgue donc pas une méthode opérationnelle complète dans le domaine de l'agriculture de précision. En effet, cette demande ne concerne que la caractérisation de la variabilité intra-parcellaire mais pas la création d'une préconisation, qui est pourtant au coeur de l'agriculture de précision. 5 La demande de brevet US2015/0206255A1 concerne une méthode de préconisation de densité de semis consistant à collecter des données constituées par des images satellites, à les convertir en les traitant pour donner un indice de biomasse (NDVI) d'une parcelle agricole pendant la période de végétation, à produire une carte de la parcelle 10 représentant le NDVI et à générer une carte de préconisation de densité de semis à partir de cette carte NDVI. Cette dernière étape est effectuée en obtenant la densité de semis maximale recommandée pour une espèce végétale donnée, en faisant correspondre une valeur de NDVI égale à un, à la densité de semis maximale recommandée, et une valeur de NDVI égale à zéro, à la densité de semis minimale recommandée. Les densités de 15 semis pour les valeurs intermédiaires sont interpolées en prenant pour base la mise à l'échelle "scaling" des données NDVI et la carte NDVI. Enfin, la carte de préconisation de densité de semis variable est générée sur la base de la densité de semis variable définie pour la parcelle agricole considérée. La densité de semis maximale recommandée pour une espèce végétale pour une variété végétale donnée, peut être 20 remplacée par la densité de semis maximale pour cette espèce végétale issue de l'expérience de l'agriculteur. Cette méthode d'amélioration de la productivité au moyen d'une carte de préconisation de densité de semis, repose sur un seul paramètre, le NDVI, propre aux plantes cultivées sur la parcelle, et, s'agissant de l'espèce végétale considérée, seulement sur les recommandations du semencier, lesquelles se fondent sur 25 l'expérience de l'agriculteur pour cette espèce végétale. Aucun modèle de culture n'est employé pour simuler la croissance des plantes et déterminer les densités optimales de semis de la variété considérée en fonction du contexte pédo-climatique local. Ces limitations quant aux paramètres pris en compte, affectent la fiabilité de la méthode proposée dans cette demande de brevet US2015/0206255A1. 30 En outre, du fait de la prise en compte d'un historique de culture restreint, cette méthode n'est ni suffisamment représentative, ni suffisamment précise. Force est donc de constater que l'existant ne propose pas de procédé et de moyens associés permettant d'augmenter la productivité d'une parcelle agricole, qui soit basée 35 sur des méthodes concrètes et efficaces (non virtuelles) multiparamétriques, simples, fiables, inter- opérables et compatibles avec les différents formats de fichiers utilisables par les engins agricoles délivrant des intrants sur les cultures. 4 Objectifs de l'invention La présente invention vise à satisfaire au moins l'un des objectifs énumérés ci-après. o fournir un procédé perfectionné pour accroître la productivité d'une parcelle agricole (P) vis-à-vis d'une variété végétale d'une espèce donnée, en améliorant 5 la gestion des intrants, tels que le semis par la modulation de sa densité; o fournir un procédé perfectionné pour accroître la productivité d'une parcelle agricole (P) vis-à-vis d'une variété végétale d'une espèce donnée, en améliorant la gestion des intrants, tels que le semis par la modulation de sa densité; ce procédé permettant en particulier une bonne modulation des intrants en 10 apportant la bonne densité de semis d'intrants, au bon moment et au bon endroit ; o fournir un procédé perfectionné pour accroître la productivité d'une parcelle agricole (P) vis-à-vis d'une variété végétale d'une espèce donnée, en améliorant la gestion des intrants, tels que le semis par la modulation de sa densité; ce procédé étant d'une grande spécificité vis-à-vis de la variété végétale considérée, 15 vis-à-vis des qualités intrinsèques du sol de la parcelle agricole concernée et vis- à-vis de l'environnement climatologique de ladite parcelle ; o fournir un procédé perfectionné, simple et fiable, pour accroître la productivité d'une parcelle agricole (P) vis-à-vis d'une variété végétale d'une espèce donnée; o fournir un procédé perfectionné, simple et fiable d'apport d'au moins un intrant 20 sur la parcelle (P) destinée à produire une variété végétale d'une espèce donnée sur une parcelle (P), selon une carte de préconisation; o fournir un procédé perfectionné, simple et fiable de semis d'une variété végétale d'une espèce donnée sur une parcelle (P), selon une carte de préconisation de la densité de semis sur la parcelle (P); 25 o fournir un procédé perfectionné, simple et fiable de réalisation d'une préconisation adaptée, de préférence d'une carte de préconisation, d'apport d'au moins un intrant sur la parcelle; o fournir un logiciel pour la mise en oeuvre du procédé selon l'invention ; o fournir un support de stockage d'une carte de préconisation d'apport d'un intrant, 30 de préférence le semis par modulation de sa densité, sur (P); un support de stockage du logiciel, o fournir un dispositif d'apport de l'intrant, en particulier un semoir ; o fournir un module informatique comprenant le logiciel pour la mise en oeuvre du procédé selon l'invention et destiné à équiper le dispositif d'apport de l'intrant, en 35 particulier le semoir. 5 Brève description de l'invention Ces objectifs, parmi d'autres, sont atteints par la présente invention qui concerne tout d'abord un procédé pour augmenter la productivité d'une parcelle agricole (P) vis-à-vis d'une variété végétale d'une espèce donnée, caractérisé en ce qu'il consiste 5 essentiellement : a. à collecter des données géoréférencées sur P qui rendent compte de la productivité de différentes zones de P; ces données appartenant à au moins un type Tn (n=entier positif), ces types Tn étant de préférence choisis dans le groupe comprenant: 10 Ti. des données géophysiques de P, de préférence conductivité/ré- sistivité, T2. des données topographiques de P; T3. des cartes de rendements de P; T4. des indicateurs de la biomasse végétale produite sur P (indices de 15 végétation) obtenus par imagerie aérienne, de préférence le Normalized Difference Vegetation Index (NID VI); b. à traiter ces données Tn pour les regrouper sur P en utilisant des méthodes géostatistiques appropriées; c. à décomposer P, pour tout ou partie des regroupements -de préférence 20 pour chaque regroupement- obtenus à l'étape (b), en différentes zones de gestion (Zg) et classifier ces différentes zones (Zg) en différentes classes Ck (k=entier positif), par classification automatique statistique; d. de préférence à regrouper les zones de gestion (Zg) de P, de manière à supprimer les zones de gestion (Zg) dont la surface ou la plus grande 25 dimension, est inférieure à une limite (L) fixée par les contraintes techni- ques liées au matériel agricole et aux pratiques agricoles mis en oeuvre sur la parcelle P; e. à déterminer un Indice de Productivité Adimensionnel (IPA) pour chacune des classes Ck affectées aux zones de gestion (Zg) de P, en 30 attribuant par pondération à chaque classe Ck, une valeur représentative de sa contribution à la productivité de la parcelle; f à fixer un Potentiel de Rendement sur P (PRP) de l'espèce considérée à partir d'une moyenne des rendements historiques sur P pour l'espèce considérée; 35 g. à calculer, à partir de l'IPA obtenu à l'étape (e) et du PRP obtenu à l'étape (f) et à l'aide d'une fonction mathématique donnée, le Potentiel de 6 Rendement (PR) de l'espèce considérée, pour chacune des zones de gestion (Zg); h. à déterminer une Densité de Semis (DS), pour la variété de l'espèce végétale considérée et pour chacune des classes Ck, à partir d'une 5 relation mathématique donnant DS en fonction de PR et dont au moins une partie des constantes est obtenue sur la base de simulations produites à partir d'un modèle écophysiologique: - qui est calibré à partir de données expérimentales, - qui reproduit le comportement de la variété de l'espèce végétale 10 considérée en fonction de conditions environnementales, de préférence sol et météorologie, - et qui permet, d'une part, de simuler les rendements de cette variété de l'espèce végétale en fonction de différentes conditions environnementales et, d'autre part, de déterminer les densités de 15 plantes optimales pour la variété considérée; i. à produire au moins une carte de P sur laquelle sont reportées les DS obtenues à l'étape (h) pour toutes les zones de gestion (Zg), cette carte étant une carte de préconisation de DS; j. à semer la variété végétale considérée sur la parcelle P en se conformant 20 à la carte de préconisation de DS de l'étape (i). Le procédé selon l'invention est avantageux en ce qu'il intègre: un ou plusieurs types différents Tn de données propres à la parcelle P et reflétant son potentiel de production par rapport à l'espèce végétale donnée, sur plusieurs 25 zones de gestion (Zg) hiérarchisées sur une échelle de productivité donnée par l'indice IPA; un potentiel de rendement (PR) pour chaque zone (Zg), issu de la combinaison de PIPA et d'un potentiel de rendement sur P (PRP) historique pour l'espèce considérée; 30 et des données de simulation de croissance de la variété végétale donnée, produites par un modèle éco-physiologique calibré à partir de données expérimentales. Cette approche nouvelle combinant des informations diverses et précises sur la variété 35 végétale et la parcelle, est particulièrement spécifique de la variété végétale donnée et de la parcelle. Elle donne accès à une préconisation intra-parcellaire localisée des plus 7 pertinentes, pour l'apport d'un intrant à la parcelle, en l'occurrence le semis par modulation de sa densité. Cela permet à l'agriculteur de procéder à une bonne et une fine modulation de la dose de semence, à délivrer parcelle par parcelle, pour optimiser le bilan économique, 5 agronomique et environnemental. Cette approche nouvelle permet de préconiser une densité adaptée à l'interaction historique sol-climat (sur plusieurs années), ce qui rend cette méthode précise et représentative. 10 Sur le même modèle de mise en oeuvre, il est envisageable de remplacer l'intrant : semis avec modulation de la densité de semis, par un autre intrant, tel que la dose de fertilisant(s) par unité de surface parcellaire. Selon un autre de ses aspects, l'invention concerne un support de stockage non 15 transitoire lisible par ordinateur, stockant au moins une version numérique d'au moins une carte de préconisation de DS pour une parcelle P, ladite carte étant telle que produite à l'étape (i) du procédé selon l'invention. Selon un autre de ses aspects, l'invention concerne un programme d'ordinateur 20 comprenant une série d'instructions qui, lorsqu'elles sont exécutées par un processeur, mettent en oeuvre un procédé selon l'invention. Selon un autre de ses aspects, l'invention concerne un support de stockage non transitoire lisible par ordinateur, stockant un programme d'ordinateur comprenant une 25 série d'instructions qui, lorsqu'elles sont exécutées par un processeur, mettent en oeuvre un procédé selon l'invention. Selon un autre de ses aspects, l'invention concerne un dispositif-semoir automatisé géolocalisable et équipé pour lire une carte numérique de préconisation de DS produite 30 à l'étape (i) du procédé selon l'invention. Selon un autre de ses aspects, l'invention concerne, un module informatique inclus ou susceptible d'être inclus dans le dispositif-semoir selon l'invention et permettant la géolocalisation dudit dispositif-semoir et le semis d'une variété végétale de l'espèce 35 donnée sur une parcelle P, en se conformant à une carte numérique de préconisation de DS produite à l'étape (i) du procédé selon l'invention. 8 Description détaillée de l'invention Définitions La terminologie employée dans le présent exposé est explicitée au travers des 5 définitions qui suivent et qui sont toutes des exemples non limitatifs: Un singulier désigne également un pluriel et réciproquement. " variété végétale" désigne une population artificielle obtenue en vue de son usage en agriculture, correspondant à la subdivision la plus basse du règne végétal. Elle doit être reproductible, homogène et stable dans ses caractéristiques génétiques et distincte des 10 autres variétés. "espèce végétale" désigne le groupe d'individus ayant des caractères morphologiques, physiologiques et chromosomiques semblables et qui peuvent se croiser entre eux "Modèle écophysiologique" : il s'agit d'algorithmes basés sur la connaissance écophysiologique des cultures, qui permettent de calculer, entres autres, à partir notamment 15 d'une description de la culture, du sol de la parcelle et de variables climatiques: -l'accumulation de biomasse par photosynthèse et donc le rendement de la culture, -les principaux besoins en éléments minéraux afférents à cette photosynthèse, dont ceux en eau et en éléments fertilisants (et donc de raisonner l'irrigation et la fertilisation). 20 -I- Procédé pour augmenter la productivité de la parcelle P Etape (a): Collecte des données de type(s) Tn Il s'agit de constituer une base de connaissances de la parcelle P en recueillant des 25 données de type(s) Tn, de natures diverses. Ces données de type(s) Tn sont géoréférencées, c'est-à-dire qu'à chaque donnée correspond avantageusement des coordonnées (x,y) sur P. Les données géophysiques de type Ti telles que la résistivité ou la conductivité du sol 30 de différentes unités de surface de la parcelle P, peuvent être mesurées à différentes profondeurs, par exemple entre 0 et 50 cm, entre 0 et 100 cm et entre 0 et 170 cm. Les données topographiques de type T2 peuvent être, par exemple, l'altitude centrée, la pente... de différentes unités de surface de la parcelle P. De telles données peuvent 35 notamment être acquises grâce à un système de positionnement par satellite de type GPS ou Galileo combiné avec la technique de positionnement dénommée "Cinématique temps réel (Real Time Kinematic, en anglais ou RTK)" ou grâce à un système de 9 télédétection par laser ou LIDAR, acronyme de l'expression en langue anglaise "light detection and ranging" ou "laser detection and ranging", est une technique de mesure à distance fondée sur l'analyse des propriétés d'un faisceau de lumière renvoyé vers son émetteur. 5 Les données de type T3 sont des cartes de rendements intraparcellaires de culture sur P de l'espèce végétale considérée ou d'autres espèces végétales, de préférence sur au moins 3 années de culture. 10 Sur une caractéristique remarquable de l'invention, les données T4, à savoir les indices de végétation, de préférence le NDVI, sont obtenus à partir d'au moins 2 - idéalement au moins 4 - images aériennes numérisées - de préférence satellites - de P, collectées sur au moins 2 - idéalement au moins 3 - années. Ces indices de végétation sont obtenus grâce à des capteurs de mesure embarqués dans 15 des satellites d'observation de la terre. Ces mesures sont par exemple celle du rayonnement solaire réfléchi dans le rouge et le proche infrarouge pour suivre l'activité photosynthétique et la surface foliaire ou la biomasse verte. Les plantes chlorophylliennes absorbent les longueurs d'onde du rouge pour alimenter la photosynthèse tandis que les structures cellulaires des parois réfléchissent fortement le 20 proche infrarouge. La combinaison de longueur d'onde permet de produire différents indices de végétation tels que le "Normalized Difference Vegetation Index" NDVI. Cet indice, parmi d'autres (tels que "Chlorophyll Vegetation Index" CVI, "Enhanced Vegetation Index" EVI, "Soil Adjusted Vegetation Index" SAVI, " Green Normalized Difference Vegetation Index " GNDVI, ....), peut être un témoin de la croissance de la 25 biomasse au sein de différentes parcelles P. Les données de type T4 étant, selon un mode préféré de mise en oeuvre de l'invention, des images, il est prévu, conformément à une caractéristique remarquable de l'invention, de retirer une zone tampon (buffer) des images concernées afin de limiter les effets de bordure, lors de l'interprétation desdites images. 30 La largeur de la zone tampon retirée est avantageusement déterminée en fonction de la résolution des images utilisées. Comme autres exemples de données de type(s) Tn sur la parcelle, on peut citer les données historiques sur le climat de la parcelle, les données historiques sur les 35 différentes cultures effectuées sur cette parcelle, les découpages et les remembrements historiques qui ont conduit à la délimitation de la parcelle, les observations visuelles permettant de détecter la nature géologique de la parcelle, les mauvaises herbes se 10 développant sur la parcelle, les réflectances mesurées à différentes longueurs d'ondes sur la parcelle à l'aide de capteurs placés à proximité du sol (N-sensore), la composition physico-chimique du sol de la parcelle, la structure, la texture, la granulométrie, l'humidité, la capacité d'échange cationique les taux d'éléments fertilisants, les éléments 5 traces et la teneur en matière organique ... Etape (b): Regroupement des différentes données de 0?pe(s) Tn sur P Les données de type(s) Tn sont sous forme de fichiers informatiques de différents 10 formats (e.g. rasters, vecteurs...) qui peuvent être stockés sur des supports différents selon leur nature, leur résolution ou leurs systèmes de projection géographiques. Ces données sont regroupées sur un repère géospatial unique en utilisant des méthodes géostatistiques appropriées, par exemple le krigeage ordinaire. Si ces données sont déjà représentées sur un repère géospatial unique, cette étape n'est plus nécessaire. 15 Ainsi, la position sur P d'une donnée de type Tn est avantageusement exprimée avec le même système de localisation que celui de toutes les autres données Tn exploitées. Le regroupement s'opère de préférence sur une grille Cette grille, qui se superpose au contour de la parcelle P présente un maillage régulier. Cette grille est constituée de points équidistants géoréférencés dans un système de projection en coordonnées 20 cartésiennes approprié à la localisation géographique de P. La distance entre les points est déterminée en fonction de la surface parcellaire ou de la résolution des données Tn. Idéalement cette distance est de 5 mètres, voire 10 mètres ou plus. Toutes les données Tn disponibles sont reprojetées dans le même système en coordonnées cartésiennes. Toutes les données Tn disponibles reprojetées sont ensuite interpolées sur la grille en 25 utilisant des méthodes géostatistiques appropriées, par exemple le krigeage ordinaire. La grille regroupant les données caractéristiques de P qui sont implémentées dans les étapes suivantes du procédé, peuvent être virtuelles et/ou peuvent faire l'objet d'un affichage sous forme d'images sur un support tel qu'un écran ou sur un élément matériel susceptible de recevoir une impression (feuille en papier). 30 Etape (c): Décomposition de P dans la grille obtenue à l'étape (b) en différentes zones de gestion Cette étape est l'une des étapes essentielles du procédé selon l'invention. Elle consiste à 35 générer une partition optimale de la parcelle P, pour tout ou partie des regroupements - de préférence pour chaque regroupement (e.g. sur la grille) de l'étape (b), au travers d'un 11 ou plusieurs types de données Tn, de façon à mettre en évidence l'hétérogénéité intraparcellaire. Pour ce faire, on procède à une décomposition de la parcelle en un certain nombre de zones de gestion (Zg) homogènes vis-à-vis d'un ou plusieurs types de données Tn. 5 Cette décomposition est réalisée par classification automatique statistique de P en (x) différentes zones de gestion (Zg) homogènes dans différentes classes Ck (k=entier positif). Chaque classe Ck se définit au travers d'une valeur ou d'un état sur lequel sont centrées l'ensemble des données Tn qui constitue une zone de gestion (Zg) homogène ou au travers d'un intervalle de valeurs ou d'états comprenant l'ensemble des données Tn qui 10 constitue une zone de gestion (Zg) homogène. Par exemple, si l'on considère des données Ti de résistivité Ro du sol de la parcelle P, comprises entre 0 et 500 am , on peut définir les 4 classes C1 C2 C3 C4 suivantes : Ci: 0 REVENDICATIONS1. Procédé pour augmenter la productivité d'une parcelle agricole (P) vis-à-vis d'une variété végétale d'une espèce donnée, caractérisé en ce qu'il consiste essentiellement : a. à collecter des données géoréférencées sur P et qui rendent compte de la productivité de différentes zones de P; ces données appartenant à au moins un type Tn (n=entier positif), ces types Tn étant de préférence choisis dans le groupe comprenant: T1. des données géophysiques de P, de préférence conductivité/ré- sistivité, T2. des données topographiques de P; T3. des cartes de rendements de P; T4. des indicateurs de la biomasse végétale produite sur P (indices de végétation) obtenus par imagerie aérienne, de préférence le Normalized Difference Vegetation Index (NDVI); b. à traiter ces données Tn pour les regrouper sur P en utilisant des méthodes géostatistiques appropriées; c. à décomposer P, pour tout ou partie des regroupements -de préférence pour chaque regroupement- obtenus à l'étape (b), en différentes zones de gestion (Zg) et classifier ces différentes zones (Zg) en différentes classes Ck Oc-entier positif), par classification automatique statistique; d. de préférence à regrouper les zones de gestion (Zg) de P, de manière à supprimer les zones de gestion (Zg) dont la surface ou la plus grande dimension, est inférieure à une limite (L) fixée par les contraintes techniques liées au matériel agricole et aux pratiques agricoles mis en oeuvre sur la parcelle P; e. à déterminer un Indice de Productivité Adimensionnel (IPA) pour chacune des classes Ck affectées aux zones de gestion (Zg) de P, en attribuant par pondération à chaque classe Ck, une valeur représentative de sa contribution à la productivité de la parcelle; f. à fixer un Potentiel de Rendement sur P (PRP) de l'espèce considérée à partir d'une moyenne des rendements historiques sur P pour l'espèce considérée; g. à calculer, à partir de PIPA obtenu à l'étape (e) et du PRP obtenu à l'étape (f) et à l'aide d'une fonction mathématique donnée, le Potentiel 30 de Rendement (PR) de l'espèce considérée, pour chacune des zones de gestion (Zg); h. à déterminer une Densité de Semis (DS), pour la variété de l'espèce végétale considérée et pour chacune des classes Ck, à partir d'une 5 relation mathématique donnant DS en fonction de PR et dont au moins une partie des constantes est obtenue sur la base de simulations produites à partir d'un modèle écophysiologique: - qui est calibré à partir de données expérimentales, - qui reproduit le comportement de la variété ou de l'espèce 10 végétale considérée en fonction de conditions environnementales, de préférence sou météorologie, - et qui permet, d'une part, de simuler les rendements de cette variété ou de cette espèce, en fonction de différentes conditions environnementales et, d'autre part, de déterminer les densités de 15 plantes optimales pour la variété considérée; j. à produire au moins une carte de P sur laquelle sont reportées les DS obtenues à l'étape (h) pour toutes les zones de gestion (Zg), cette carte étant une carte de préconisation de DS. J. 20 2. Procédé selon la revendication 1 caractérisé en ce qu'il comprend une étape (j) consistant à semer la variété végétale considérée sur la parcelle P en se conformant à la carte de préconisation de DS de l'étape (i). 25 I 3. Procédé selon la revendication 1 caractérisé en ce que les données collectées à l'étape a) sont de type T3 et sont des rendements historiques de culture sur P de l'espèce végétale considérée, de préférence sur au moins 3 années de culture. 4. Procédé selon l'une au moins des revendications précédentes caractérisé en ce 30 que les données collectées lors de l'étape a) sont de type T4, à savoir des indicateurs de la biomasse végétale produite sur P, de préférence le NDVI, sont obtenus à partir d'au moins 2 - idéalement au moins 4 - images aériennes numérisées - de préférence satellites - de P, collectées sur au moins 2 - idéalement au moins 3- années. 35 5. Procédé selon l'une au moins des revendications précédentes caractérisé en ce que l'étape (c) de décomposition de P en différentes zones de gestion (Zg), est 31 effectuée à partir d'au moins deux types Tn,Tn' de données, de préférence en combinant sur une même grille définissant un maillage régulier dans la parcelle P, les données d'au moins deux types Tn,Tn' de données, par interpolation géostatistique. 5 6. Procédé selon l'une au moins des revendications précédentes caractérisé en ce que la pondération de l'étape (e) qui permet l'attribution, à chaque classe Ck, d'un IPA, consiste essentiellement à soumettre tout ou partie des regroupements [de préférence chaque regroupement], comprenant les zones de 10 gestion (Zg) obtenues à partir d'un type ou de différents types de données Tn, à une transformation centrée réduite de la (des) donnée(s) Tn composant le regroupement, puis à calculer, pour chacune des classes Ck , la moyenne arithmétique de la (des) donnée(s) Tn centrée réduite, ladite moyenne constituant l'IPA. 15 7. Procédé selon l'une au moins des revendications précédentes caractérisé en ce que dans l'étape (g), le potentiel de rendement (PR) est déterminé pour chacune des classes par une fonction mathématique qui intègre PIPA et le PRP. 20 8. Procédé selon l'une au moins des revendications précédentes caractérisé en ce que, la relation mathématique mise en oeuvre dans l'étape (h) est la suivante: DS b,'- [(MaxPot) -PR) * a]; dans laquelle: a est la pente d'une droite de régression de DS de la variété considérée 25 en fonction du PR; MaxPot est le PR maximum donné par cette droite en conditions peu limitantes; bn, est l'ordonnée de la droite à MaxPot; a, b,', et MaxPot étant de préférence obtenu par simulation à partir d'un 30 modèle écophysiologique ou à partir de données expérimentales. 9. Procédé selon l'une au moins des revendications précédentes caractérisé en ce qu'une version numérique de la carte de préconisation de DS produite à l'étape (i) est transmise à un dispositif-semoir automatisé géolocalisable et équipé de 35 moyens informatiques pour lire cette carte numérique et pour effectuer le semis de la variété végétale considérée sur la parcelle P en se conformant à la carte de préconisation de DS. 32 10. Programme d'ordinateur comprenant une série d'instructions qui, lorsqu'elles 5 sont exécutées par un processeur, mettent en oeuvre un procédé selon l'une au moins des revendications 2 à 9. 11. Support de stockage non transitoire lisible par ordinateur, stockant un programme d'ordinateur selon la revendication 10. 10 12. Dispositif-semoir automatisé géolocalisable et équipé pour lire une carte numérique de préconisation de DS produite à l'étape (i) du procédé selon l'une au moins des revendications 1 à 9. 15 13. Module informatique inclus ou susceptible d'être inclus dans le dispositif- semoir selon la revendication 12 et permettant la géolocalisation dudit dispositif-semoir et le semis d'une variété végétale de l'espèce donnée sur une parcelle P en se conformant à une carte numérique de préconisation de DS produite à l'étape (i) du procédé selon l'une au moins des revendications 1 à 9. 20