Code Python pour créer les Moyades en 3D
# -*- coding: utf-8 -*- """ Visualisation interactive du MNT autour de l’îlot des Moyades Auteur : fcart Date : juillet 2025 """ # === IMPORTS === import pyvista as pv import numpy as np import rasterio # === LECTURE DU FICHIER === path = r"C:\Users\bureau\geotif\0890_6235_Moyades\0890_6235\LITTO3D_FRA_0892_6234_20150529_LAMB93_RGF93_IGN69\MNT1m\LITTO3D_FRA_0892_6234_MNT_20150128_LAMB93_RGF93_IGN69.asc" with rasterio.open(path) as src: data = src.read(1) transform = src.transform nodata = src.nodata # === MASQUAGE DES DONNÉES INEXISTANTES === data = np.where(data == nodata, np.nan, data) # === CONSTRUCTION DES GRILLES X ET Y === nrows, ncols = data.shape x0 = transform.c y0 = transform.f dx = transform.a dy = transform.e # attention : négatif x = np.arange(ncols) * dx + x0 y = np.arange(nrows) * dy + y0 xx, yy = np.meshgrid(x, y) # === DÉFINITION MANUELLE DE LA ZONE D'INTÉRÊT (coordonnées Lambert 93) === # Ajustage des bornes s...