Vignes arrachées, Entre-deux-Mers

01 · Le problème

Un angle mort : les vignes qui sortent des déclarations

La crise du vignoble bordelais a conduit à des arrachages massifs, aidés ou non. Les aides imposent ensuite de renaturer ou de diversifier les parcelles, et collectivités, conseillers et opérateurs fonciers cherchent à savoir où sont ces terres et ce qu'elles deviennent. La source ouverte la plus complète, le Registre parcellaire graphique (RPG) des déclarations PAC, n'y répond qu'en partie.

≈ 21 000 ha

de vignes arrachés en Gironde d'août 2023 à mai 2026 : 12 284 ha avec une aide publique, environ 8 700 ha sans aide.

Agreste Nouvelle-Aquitaine, Essentiel n°80, juin 2026.

34,9 %

des parcelles déclarées en vigne en 2022 ne sont plus déclarées du tout à la PAC en 2025. 53,5 % sont toujours en vigne, 11,6 % déclarées autre chose.

Agreste Nouvelle-Aquitaine, Essentiel n°80, juin 2026 (données PAC 2025).

≈ 4 400 ha

proposés à la Foncière d'avenir en Gironde (430 dossiers, plus de 6 500 parcelles), qui cherche à identifier le foncier libéré par l'arrachage. Les friches viticoles font partie des premiers critères retenus.

Préfecture de Nouvelle-Aquitaine, 19 juin 2026 ; R. Lecocq, Pleinchamp, 15 juillet 2026.

Une parcelle sortie de la vigne dans la PAC peut cacher trois réalités

Cas 1

Reconvertie et déclarée

Jachère, prairie, culture : la déclaration PAC dit ce qu'elle est devenue.

Devenir connu

Cas 2

Toujours en vigne

Plus de demande d'aide, changement d'exploitant : la vigne est toujours là. La compter comme arrachée enverrait sur le terrain pour rien.

Fausse piste à écarter

Cas 3

Arrachée, plus rien de déclaré

Entretenue, en friche, en train de se boiser ? La PAC ne dit plus rien. C'est l'angle mort que l'outil éclaire.

Devenir à découvrir

Trois questions pour décider

La zone d'étude

Un carré de 40 km de côté (1 600 km²), centré sur l'Entre-deux-Mers (0,20° ouest, 44,75° nord). Il déborde sur le Libournais au nord de la Dordogne et sur les Graves à l'ouest de la Garonne. Le calcul tourne sur 81 dalles de 5 km, traitées une par une.

La méthode s'applique à une autre zone en changeant son centre et sa taille : rien n'est propre à ce carré.

Carte d'occupation du sol à 10 m de la zone d'étude — La zone vue par le modèle Sentinel-2 : forêt (vert foncé), vignes et vergers (violet), prairies (vert clair), cultures (jaune), eau (bleu), bâti (rouge). Carte produite par le pipeline.

02 · Les données

Uniquement des données ouvertes, croisées parcelle par parcelle

Le script télécharge tout seul chaque source (services de l'IGN, catalogue d'images Sentinel-2, API naturalistes, Géorisques). Les contrôles sont faits par photo-interprétation, pas sur le terrain.

Donnée	Producteur	Résolution, dates	Rôle dans l'outil
Registre parcellaire graphique (RPG)	IGN, ASP	Parcelles déclarées, 2018 à 2024	Repérer les parcelles déclarées en vigne puis plus en vigne ; lire ce qui est déclaré à la place.
Sentinel-2 L2A	Copernicus, via Microsoft Planetary Computer	10 m, 10 bandes, 2023 à 2025 ; scènes à moins de 70 % de nuages ; environ observations sans nuage par pixel	Carte d'occupation du sol, fauches, sol nu, évolution du NDVI.
BD ORTHO	IGN	20 cm (lue à 25 cm par pixel), prise de vue du 11 mai 2024	Voir les rangs de vigne.
BD TOPO, BD Forêt, OCS GE	IGN	Vecteur	Exemples d'apprentissage du modèle (forêt, eau, bâti, cultures), distance à la forêt.
Zonages et milieux	INPN (PatriNat), IGN	ZNIEFF 1 et 2, Natura 2000, arrêtés de biotope, réserves, zones humides et prairies sensibles BCAE, haies	Indice écologique et enjeu principal.
Espèces observées	GBIF, OpenObs (MNHN)	Mailles de 2 km, observations de 2000 à 2026 ; espèces quasi menacées à en danger critique (UICN)	Composante « espèces menacées » de l'indice.
Eau, sols, risques	BD TOPO, carte des sols INRAE, SoilGrids (ISRIC), Géorisques	Cours d'eau, hydromorphie, pH et carbone modélisés, zones inondables, argiles, cavités	Contraintes des projets, analyses de sol conseillées.
Urbanisme	Géoportail de l'urbanisme	Zonage des PLU, servitudes	Contexte réglementaire dans chaque fiche.
Contrôles à l'œil	Photo-interprétation	138 parcelles sur l'orthophoto de mai 2024, 131 sur des images de 2026	Caler le modèle, corriger l'estimation de surface.
Arrachage sanitaire, 1re vague	DDTM de la Gironde (data.gouv.fr)	Parcelles cadastrales, demandes de fin 2023, Licence Ouverte 2.0	Validation indépendante.

Six variables Sentinel-2 sur une dalle de 5 km — Une dalle de 5 km vue par Sentinel-2 : couleurs naturelles, NDVI d'été, part du temps en sol nu, fauches, tendance du NDVI, années depuis la dernière gestion. Six des 108 variables du modèle.

Extrait réel de l'orthophoto puis le même extrait en pixels de 10 m — À gauche, un extrait réel de l'orthophoto IGN de mai 2024 (340 m de côté, réduit) : une vigne striée de rangs en haut à droite, une parcelle arrachée devenue jachère dans le contour jaune. À droite, le même extrait moyenné en pixels de 10 m, la résolution de Sentinel-2 : les rangs disparaissent.

03 · La méthode

Cinq étapes, avec des règles explicites

Après le point de départ (repérer les parcelles sorties de la vigne), cinq étapes, numérotées comme dans la vidéo. Chacune écrit ses résultats sur le disque, se relance seule et se vérifie. Clique sur une étape pour le détail : ce qu'elle fait, avec quels paramètres, et pourquoi.

Étape 3 en pratique

Essayer le score de rangs

La démo refait le calcul du script en direct sur une image simulée de 64 m de côté. Des rangs réguliers donnent un point lumineux dans le spectre (à droite) : plus il est net par rapport au bruit voisin, plus le score monte.

Image simulée, 64 × 64 m

Spectre. Anneaux : périodes de 4 m et 1 m. Cercle vert : pic retenu.

Score de rangs

…

Écartement des rangs 2,0 m Orientation 30° Contraste du feuillage 1,0 Bruit, ombres 0,5

Étape 5 en détail

Les règles de qualification

Ce sont des règles expertes, écrites dans le script avec leurs seuils. Elles ne sont pas calibrées sur des données de terrain : elles servent à ouvrir la discussion, à valider avec des conseillers.

04 · Les résultats

Les réponses aux trois questions

Question 1

Lesquelles sont vraiment arrachées ?

L'entonnoir suit les étapes de la méthode. Chaque barre représente toute la surface sortie de la vigne PAC, découpée un peu plus finement à chaque étape.

Hectares arrondis à l'unité : une somme d'arrondis peut différer d'un hectare du total.

Ce que la carte compte, et ce qu'on estime

Pourquoi tant de vignes « disparues » sont encore là

Une parcelle peut sortir de la vigne dans la PAC sans être arrachée, par exemple si elle n'est simplement plus déclarée. L'outil ne dit pas pourquoi : il constate que, sur l'orthophoto de mai 2024, ces parcelles ont des rangs nets ou un signal de culture pérenne. Sans ce tri, plus de la moitié des parcelles d'une liste « vignes arrachées » tirée de la PAC seraient de fausses pistes ().

Question 2

Que sont-elles devenues ?

Avec le temps, les parcelles sortent de la PAC

Devenir de la surface arrachée selon le nombre d'années depuis la dernière vigne déclarée (état lu en 2024).

Question 3

Qu'en faire, et par où commencer ?

Trois lectures complémentaires : ce qu'il faut préserver, ce qui est réaliste sur le terrain, et où agir d'abord.

Clique sur une ligne pour voir la règle exacte.

Ce qui est réaliste, par famille de projets

Ce qu'il faut retenir

Avant tout projet sur une ancienne vigne, chaque fiche conseille des analyses de sol, dont le cuivre hérité des traitements, et un sondage à la tarière pour l'hydromorphie.

Ce qui bloque chaque famille de projets — Ce qui bloque, par famille de projets. En vert, les freins qu'on peut lever : une parcelle trop petite se regroupe avec ses voisines, un sol hydromorphe (gorgé d'eau une partie de l'année) se draine, un point d'eau trop loin se compense par l'irrigation. En gris, les freins difficiles à lever : pente forte, sol engorgé en permanence, zone inondable. Une parcelle peut cumuler plusieurs freins.

Les 20 parcelles à regarder en premier

Score sur 100 : grande surface 35 %, enjeu écologique 25 %, années sans gestion 25 %, regroupable avec ses voisines 15 %. Clique pour voir la parcelle sur la carte.

Concentration de la surface et blocs regroupables — À gauche : la moitié des parcelles concentre 80 % de la surface arrachée. À droite : les blocs de parcelles distantes de moins de 100 m, utiles pour monter un projet à l'échelle d'un îlot.

Arrachages par commune et par année (première année sans vigne déclarée), et devenir selon l'ancienneté.

05 · Fiabilité et limites

Peut-on s'y fier ? Deux contrôles, des limites mesurées

Une liste de parcelles n'est utile que si l'on sait à quel point s'y fier. Premier contrôle : comparer l'outil à mon œil, sur les mêmes images. Second contrôle : le confronter à des arrachages connus par ailleurs.

1. Contrôle à l'œil, à deux dates

Sur un échantillon de parcelles, je note ce que montre l'image (vigne ou arrachée) et je compare au verdict de l'outil.

Comment on estime la surface sans croire la carte sur parole

Méthode d'Olofsson et al. (2014), standard en télédétection : les parcelles candidates sont rangées en 4 cases (rangs détectés ou non, signal pérenne du satellite ou non). Dans chaque case, la part réellement arrachée mesurée sur les vérifications est appliquée à toute la surface de la case. Les erreurs de l'outil sont ainsi corrigées. La marge vient d'un bootstrap : le calcul est refait 2 000 fois en rééchantillonnant les parcelles vérifiées, et on garde les 95 % de résultats du milieu.

Matrices de confusion et estimation de surface — Figure produite par le pipeline.

Trois parcelles types vues à quatre dates — À quoi ressemble un contrôle à l'œil : trois parcelles types, dans l'ordre des dates (SPOT 2023, SPOT 2024, orthophoto de mai 2024, SPOT 2025). Ligne 1 : arrachée, devenue jachère. Ligne 2 : sortie du dernier RPG, rangs encore nets. Ligne 3 : plus déclarée mais toujours en vigne. Imagerie IGN.

2. Test indépendant : les arrachages sanitaires de la DDTM

La DDTM de la Gironde publie les parcelles cadastrales arrachées avec l'aide sanitaire (première vague, demandes du 20 novembre au 20 décembre 2023, pour lutter contre la flavescence dorée). Ce sont des arrachages certains, recensés sans l'outil. On regarde ce que l'outil en dit.

Répartition des surfaces DDTM selon le verdict de l'outil

Catégorie	ha	Lecture

Ce que l'outil fait bien

Ses limites, mesurées

06 · À quoi ça sert

Une aide à la décision, pas un verdict

L'outil ne décide pas à la place des gens. Il dit où regarder, ce qui est probable, ce qui bloque et ce qui reste à vérifier.

Conseillers agricoles

Chambre d'agriculture, bureaux d'études

Préparer un diagnostic de diversification avec la fiche de la parcelle : pente, sol, eau, freins.
Repérer les îlots regroupables pour un projet collectif (élevage, maraîchage).

Collectivités

Communes, intercommunalités

Suivre les parcelles sans entretien et les friches viticoles de leur territoire.
Nourrir un PLU ou une stratégie foncière avec une carte mise à jour chaque année.

Opérateurs fonciers

SAFER, Foncière d'avenir en Gironde

Repérer le foncier libéré par l'arrachage, y compris celui qui ne figure plus à la PAC.
Écarter d'emblée les vignes encore en place.

Naturalistes

Conservatoires, associations, gestionnaires d'espaces

Cibler les parcelles à enjeu : zones humides, lisières forestières, végétation spontanée qui s'installe.
Prioriser les visites de terrain.

Comment s'en servir

Ouvrir la carte et chercher sa commune.
Lire la fiche de chaque parcelle : devenir, enjeu, projets possibles et leurs raisons, priorité.
Exporter : tableau des priorités en CSV, couche QGIS avec toutes les colonnes.
Vérifier sur le terrain et en parler avec les exploitants avant toute décision.
Mettre à jour chaque année : nouveau RPG, nouvelles images, puis relancer le notebook.

Pour aller plus loin

Voir les arrachages partiels : mesurer les rangs par fenêtres de 30 m à l'intérieur de chaque parcelle, plutôt que sur la parcelle entière.
Estimer par points : un échantillon de points plutôt que de parcelles, pour corriger aussi les arrachages partiels.
Suivre dans le temps : refaire la lecture des rangs sur chaque nouvelle orthophoto.
Valider les règles de projets et d'enjeux avec des conseillers et des naturalistes, puis les recalibrer.
Étendre à toute la Gironde : le script fonctionne aussi par département.

Annexe · Le code

Sous le capot

Un seul fichier Python, script_vignes_arrachees.py (7 795 lignes, 256 fonctions), et un notebook commenté qui l'appelle étape par étape. Conçu pour Onyxia (SSP Cloud) : il installe ses bibliothèques, reprend là où il s'est arrêté et s'applique à une autre zone en changeant deux paramètres.

Étape	Fonction	Ce qu'elle fait	Ce qu'elle écrit
1. Zone	etape_zone	Carré d'étude découpé en 81 dalles de 5 km, traitées une par une	0_zone/
2. Bases	etape_bd	Télécharge RPG 2018 à 2024, BD TOPO, BD Forêt, OCS GE, zonages, sols, risques, haies, PLU	1_references/
3. Sentinel-2	etape_s2	Cherche les images 2023 à 2025 (catalogue STAC) et calcule par pixel composites saisonniers, indices, séries, fauches	2_sentinel2/
4. Exemples	etape_etiquettes	Transforme les bases en exemples d'apprentissage (9 classes), en érodant les bords	3_apprentissage/
5. Modèle	etape_modele	Forêt aléatoire de 300 arbres, nettoyage des exemples douteux, validation croisée par blocs	modele/
6. Carte	etape_prediction	Prédit chaque dalle puis assemble la mosaïque à 10 m	4_cartes/
7. Parcelles	etape_parcelles	Repère les vignes sorties du RPG, découpe les îlots, calcule devenir, enjeux, projets, priorité	parcelles_arrachees.gpkg
7b. Correction	etape_correction	Mesure les rangs sur l'orthophoto, cale le modèle, corrige, estime la surface	chiffre_cle.json
8. Rapport	etape_rapport	Couches de vérification QGIS, carte web, figures, rapport HTML	5_resultats/
Validation	verite_ddtm	Confronte l'outil aux arrachages DDTM et corrige le rappel avec la vérification	validation/
Export	exporter_site	Rassemble les données de ce site	site_donnees.zip

Choix d'ingénierie

Reprise : avec REPRENDRE = True, toute étape déjà calculée est relue. Les scores de rangs sont mis en cache parcelle par parcelle : un « Run All » relit ce qui est fait au lieu de le recalculer.
Réseau fragile : chaque requête est réessayée jusqu'à 6 fois avec une attente qui double. Une requête WFS qui renvoie le maximum d'objets est redécoupée en quatre, récursivement, jusqu'à être complète.
Mémoire : la mosaïque est stockée en entiers 16 bits et lue bande par bande ; les dalles libèrent leurs images au fur et à mesure.
Vitesse : orthophotos téléchargées en parallèle (4 fils), calculs vectorisés avec NumPy et SciPy.
Honnêteté statistique : validations croisées par blocs spatiaux (2,5 km pour la carte, 2 km pour le modèle de vigne), pour ne jamais tester sur les voisines d'un exemple d'apprentissage.
Tests : autotest() fait tourner la correction, les indices, l'estimation et le rapport sur une petite zone fictive, sans réseau.
Traçabilité : chaque exécution ajoute une ligne au journal (seuil, méthode, chiffres clés).

# extrait réel : le score de rangs (simplifié)
def score_rangs(gris, masque, res=0.25, pmin=1.0, pmax=4.0):
    a = gris - ndi.gaussian_filter(gris, 6.0 / res)   # passe-haut
    a *= np.outer(np.hanning(h), np.hanning(w))     # fenêtre
    P = np.abs(np.fft.fftshift(np.fft.fft2(a))) ** 2  # spectre
    bande = (F >= 1 / pmax) & (F <= 1 / pmin)        # 1 à 4 m
    i = np.argmax(np.where(bande, P, -1))            # pic
    # même direction (6° près), fréquences voisines
    bas  = (dt <= 6) & (F >= 0.55 * f0) & (F <= 0.8 * f0)
    haut = (dt <= 6) & (F >= 1.25 * f0) & (F <= 1.6 * f0)
    fond = np.sqrt(np.median(P[bas]) * np.median(P[haut]))
    return P.ravel()[i] / fond                   # score

# cellule de paramètres du notebook
import script_vignes_arrachees as fp
fp.POINT_LONLAT = (-0.20, 44.75)   # centre de la zone
fp.TEST_KM = 40                    # carré de 40 km
fp.REPRENDRE = True                # Run All rapide
res = fp.etape_correction()        # rangs, correction, estimation
fp.verite_ddtm()                   # validation indépendante
fp.exporter_site()                 # données de ce site

Bibliothèques : GeoPandas, Shapely, Rasterio, odc-stac et pystac, scikit-learn, SciPy, NumPy, pandas, Matplotlib. Exécution sur Onyxia, stockage MinIO.

Compétences mobilisées

Télédétection : séries Sentinel-2, indices, composites Apprentissage automatique : forêt aléatoire, régression logistique, validation croisée spatiale Traitement du signal : transformée de Fourier 2D Statistique : estimation de surface, bootstrap SIG : GeoPandas, Rasterio, services WFS et WMS, QGIS Agronomie : sols, eau, contraintes de projets Aide à la décision : enjeux, freins, priorités

Sources

D'où viennent les chiffres

Agreste Nouvelle-Aquitaine, Essentiel n°80, « La filière viticole girondine : chiffres clés », juin 2026. PDF
Préfecture de Nouvelle-Aquitaine, « Ouverture du dispositif d'arrachage sanitaire en Gironde », novembre 2023. Lien
Préfecture de Nouvelle-Aquitaine, « Foncière d'avenir en Gironde : avec 430 dossiers déposés, l'appel à manifestation d'intérêt… », 19 juin 2026. Lien
R. Lecocq, « Foncière d'avenir en Gironde : 4400 ha potentiellement concernés », Pleinchamp, 15 juillet 2026. Lien
DDTM de la Gironde, « Dispositif d'arrachage sanitaire des vignes en Gironde en renaturation et diversification, première vague », data.gouv.fr, Licence Ouverte 2.0, mis à jour le 4 juin 2026. data.gouv.fr
Olofsson P. et al. (2014), Good practices for estimating area and assessing accuracy of land change, Remote Sensing of Environment 148, p. 42 à 57. doi
Données de l'outil : IGN Géoplateforme (RPG, BD TOPO, BD Forêt, OCS GE, BD ORTHO, haies, couches BCAE, cadastre), Copernicus Sentinel-2 via Microsoft Planetary Computer, INPN (PatriNat), GBIF, OpenObs, ISRIC SoilGrids, Géorisques, Géoportail de l'urbanisme.

Un angle mort : les vignes qui sortent des déclarations

Une parcelle sortie de la vigne dans la PAC peut cacher trois réalités

Reconvertie et déclarée

Toujours en vigne

Arrachée, plus rien de déclaré