Le portail francophone de la géomatique

C'est déjà fait pour Saga Gis : impossible d'importer des XYZ, à moins qu'il yait une astuce ? j'ai essayé avec "import tables" ...

Bonjour,

il est possible de faire ça avec Grass:
- import des fichiers avec v.in.ascii, ne pas oublier l'option -z "créer un vecteur 3D" ainsi que -b "ne pas construire la topologie" (l'import sera bien plus rapide)
- assembler les vecteurs avec v.patch
- filtrer tout ce qui dépasse du sol (bâtiments, arbres, ...) avec dans l'ordre v.lidar.edgedetection, puis v.lidar.growing, puis v.lidar.correction
- interpoler le raster mnt avec v.surf.bspline

les détails pour chaque commande sont expliqués dans le manuel (onglet "manual" dans la boite de dialogue de la commande).

Grass est capable de traiter de très gros jeux de données, il faut juste faire attention à ne pas être trop optimiste sur la quantité de mémoire diponible, et ne pas hésiter à effectuer les traitements par morceaux et faire l'assemblage des dalles mnt tout à la fin.

Bonjour,
J'ai fait l'acquisition d'un MNT au format XYZ à partir du CRIGE PACA. Le problème est que je n'arrive pas à l'ouvrir. En effet j'ai essayé diverses procédures :

Import sous QGIS (via Couche - Importer une couche raster). Cela donne le message d'erreur : 

Emplacement\fichier n'est pas une source de données raster gérées
Ungridded dataset: At line 2, X spacing was 0.000000. Expected >0 value

Import sous GRASS (interface QGIS) :

r.in.xyz "input=Emplacement\fichier" method=n type=FCELL fs=space x=1 y=2 z=3 output=couche

Le fichier raster semble bien importé mais on a un gros rectangle plein de valeurs 0 de la taille du secteur (France) sauf une quinzaine de pixels à l'emplacement normal de la couche avec des valeurs z autour de 50000.

Import sous IGN Map (Fichier : Import de MNT) : rien ne s'affiche.

Import sous IGN Map (Fichier : Import de CSV/TXT) : de nombreux points s'affichent, on peut interroger les points, on trouve l'altitude, mais, pour avoir les couches en TIF le travail est bien trop long (257 fichiers)


Pour info les données se présentent ainsi.
Pour les MNT.
La structure du MNT (XYZ) est la suivante (4 premières et 4 dernières lignes du fichier) :

1004003.010 6284997.500 331.330
1004003.010 6284992.500 332.320
1004003.010 6284987.500 333.300
1004003.010 6284982.500 334.290

1008998.010 6280017.500 159.060
1008998.010 6280012.500 160.290
1008998.010 6280007.500 161.580
1008998.010 6280002.500 162.920

Il y a 257 fichiers en tout, ayant la même structure (longitude latitude Lambert 93 et altitude).

Pour le shapefile.
Un shapefile est aussi livré, c'est une grille d'assemblage des rasters, avec les champs XMIN, XMAX, YMIN, YMAX (coordonnées Lambert 93), TILE_NAME et OLD_NAME (nom du fichier XYZ actuel). C'est peut-être la clé du problème.


Y a-t-il un moyen d'importer ce raster pour pouvoir le convertir en tif?
Je vous remercie d'avance.
Cordialement,
BAUDON Étienne

Bonjour,

As-tu essayé avec Gdal ?

gdal_translate -of GTiff ton_fichier.xyz ton_fichier.tif

Si tu trouves le fichier en sortie trop lourd, tu peux même en profiter pour le réechantillonner avec l'option -outsize

La doc est dispo ici

bonjour,
je dois dire que je reste étonné par la façon dont vous inventez des problèmes.
Je reste persuadé qu'avec Savgis on résoud le problème facilement.
Si l'on veut bien m'envoyer une copie des 300 fichiers, j'en ferais volontiers la démonstration.
Bon, je suppose que cela doit faire quelques gigas, c'est pour cela que c'est intéressant.
Cordialement.

Bonjour,

Personnellement j'ai eu le même souci avec des données XYZ, je vous écrit une méthode qui a bien marché pour moi, seul problème, ça s'adapte à une utilisation d'Arcgis avec toutes ses fonctionnalités, donc je ne réponds pas ici pour les logiciels libres mais pour les utilisateurs d'Arcgis ça peut vous faire économiser de bonnes heures.
Donc,
Préparation de la donnée
Etape 1:Sélection des dalles (par exemple Marseille, ou MPM etc.)
Etape 2:Ouverture des différents fichiers xyz avec le bloc note, conversion des points par des virgules.
Etape 3:Importation des fichiers xyz dans Excel, enregistrement sous format excel 2007, décaler les colonnes vers le bas, ajouter des entêtes F1, F2 et F3
Opérations dans Arcgis
Etape 4:Dans Arcgis, aller dans add data > add XY value, désigner F1 et F2 comme les valeurs X et Y et F3 comme Z. Réitérer l'opération sur l'ensemble des fichiers, normalement on doit avoir des séries de points
Etape 5:Reprojection des points sous format Lambert 93 (optionnel si vous travaillez en CC44)
Etape 6: interpolation des dalles de points avec la méthode de votre choix (personnellement, voisin naturel puisque les données sont très régulières et également réparties), j'ai utilisé model builder pour aller plus vite
Etape 7, normalement vous avez une série de dalles sous format GRID arcgis, vous n'avez plus qu'à créer un raster dataset (à ne pas confondre avec la mosaïque, qui peut aussi être une option selon ce que vous recherchez), rentrer les paramètres de sortie (par exemple 16 bit unsigned si vous voulez éviter les valeurs sous 0, définir -9999 comme le NoData etc.)
Etape 8: Load rasters  avec le clic droit sur votre icône raster dataset, chargez vos dalles raster, normalement vous avez un assemblage noir car les statistiques n'ont pas été calculées avec des valeurs hautes automatiques qui faussent l'échantillonnage de couleur, pour y remédier (ça m'a pris quelques temps avant de comprendre), simplement clic droit sur votre raster dataset > calculate statistics et votre raster est reclassé et visuellement utilisable.
Etape 9: dernier problème, normalement l'interpolation dalle par dalle donne des limites avec du nodata entre une dalle A et une dalle B(normal puisque le calcul d'interpolation n'a pas porté sur les points hors de la dalle), pour y remédier toujours sous Arcgis, j'ai rentré cette formule dans la calculatrice raster
Con(IsNull("myraster"), FocalStatistics("myraster", NbrRectangle(2,2, "CELL"), "MEAN"), "myraster")
2,2 signifie que la recherché autour de la cellule Nodata se fait dans un rayon de 2x2 carrés (cell), à vous de la faire varier selon vos données
Normalement vous avez un MNT tout beau tout propre et utilisable pour les tâches courantes, la méthode pourrait être encore perfectionnée mais ça a bien marché pour moi.

En espérant que ça vous aide