Conversion de fichiers XYZ (avec xy régulièrement espacés) au format ESRI GRID?

16

Normalement, les consultants (enquête LiDAR) me fournissent des données selon mes spécifications; déjà au format ESRI GRID ou ASCII GRID. De cette façon, les données LiDAR peuvent être facilement utilisées avec ArcGIS. Cette fois, je n'ai que des fichiers xyz (c'est-à-dire avec des coordonnées xy espacées régulières).

Auparavant, j'utilisais un utilitaire de Min-Lang Huang nommé GridBatch qui convertissait XYZ en ASCII GRID puis en fichiers ESRI GRID. Malheureusement, l'utilitaire se bloque depuis que j'ai mis à niveau ma machine vers W7 64 bits et ArcGIS vers la version 10.

Existe-t-il une autre façon de convertir par lots des données XYZ de surface LiDAR en fichiers ESRI GRID?

(Convertir en points -> Interpoler est une option que je voudrais éviter car il y a plus de 300 tuiles à convertir).

Jakub Sisak GeoGraphics
la source
Voici l'ancien fil de discussion du forum ESRI qui a un lien vers l'outil: forums.esri.com/Thread.asp?c=93&f=986&t=268865
Jakub Sisak GeoGraphics
Merci pour tout le monde. L'utilitaire d'origine était fourni avec le code source. J'ai obtenu une autorisation du créateur pour le mettre à jour. J'ai réussi à éliminer les erreurs, à simplifier pour mes besoins et à le mettre à jour pour et pour la version 9.3 sur Vista. resources.arcgis.com/gallery/file/geoprocessing/…
Jakub Sisak GeoGraphics
1
Merci pour la réponse, ce serait vraiment utile. J'ai cependant réussi à contourner cela. Pour plus d'informations: utilisez ASCII 3D vers entité (Outils 3D Analyst> Conversion> À partir d'un fichier) en veillant à générer une sortie vers le point, puis ajoutez des données de hauteur à l'aide de l'outil Ajouter des informations Z (3D Analyst> Fonctions 3D). Ensuite, j'ai interpolé en utilisant Natural Neighbour (3D Analyst> Raster Interpolation).
veedub
1
@Jakub Parce que le lien ne fonctionne plus, ce message n'est plus une réponse valide. Je le convertis en commentaire afin que les informations ne soient pas entièrement perdues pour tout le monde (juste au cas où ESRI restaure la cible).
whuber
1
@AndreSilva - en supposant des intervalles régulièrement espacés. Il s'agit d'un poste de 5 ans. Les arpenteurs fournissent maintenant LiDAR dans une variété de formats, donc ce n'est plus un problème pour moi. Vous ne voyez pas non plus un certain nombre de messages qui ont été supprimés ici. J'ai demandé à l'auteur de l'outil le code source et, avec sa permission, j'ai créé un complément mis à jour, mais je n'en ai pas eu besoin ni l'ai mis à jour depuis la version 10 et je l'ai retiré. Par conséquent, la solution n'étant plus disponible, les messages faisant référence à cette solution ont également été supprimés.
Jakub Sisak GeoGraphics

Réponses:

8

On dirait que vous voulez faire cela en lot (ne vous en voulez pas)

Comme l'a dit STH , il semble que Global Mapper effectuera effectivement des conversions par lots . Bon prix aussi.

FME Desktop peut utiliser le transformateur RasterDEMGenerator et un peu de liaison avec le lecteur / écrivain, mais vous aurez besoin de la version Pro, non disponible dans l'extension ArcGIS Data Interop.

Si la gratuité est plus dans votre budget, gdal_grid devrait le faire. Vous devrez passer un peu de temps à déterminer l'étendue de tous vos fichiers xyz, mais ogrinfo fera rapidement cela, puis rédigera un fichier / script par lots pour l'exécuter sur chaque fichier.

Ou vous pouvez utiliser ArcGIS ASCII 3D to Feature Class (nécessite 3D Analyst) chaîné avec Point to Raster , comme décrit ici , puis batchifier le processus avec ArcPy ou ModelBuilder. Si vous n'avez pas 3D Analyst, vous pouvez utiliser n'importe quel nombre d'arcscripts ou LASTools txt2Las et las2shp pour l'obtenir dans un format avec lequel ArcGIS peut fonctionner. Je sais que vous ne vouliez pas convertir les fichiers txt en classes d'entités intermédiaires, mais cela ne demande pas trop de travail et vous avez déjà accès au logiciel.

wwnick
la source
Merci. J'ai essayé ASCII 3D vers FC puis je pointe vers raster, aussi Topo vers raster etc. Les résultats sont corrects mais le processus prend trop de temps. Les fichiers XYZ bruts originaux vers ASCII GRID donnent les meilleurs résultats sans interpolation / conversion supplémentaire. L'outil que j'ai mentionné est livré avec le code source VB6, j'ai donc pu le mettre à jour (procédure ArcCatalog VBA) et il semble très bien fonctionner. J'essaierai cependant Global Mapper. Comme ce qu'il peut faire d'après ce que j'ai lu à ce sujet.
Jakub Sisak GeoGraphics du
Si vous allez essayer FME, je vous suggère d'utiliser FME2011 (www.safe.com/beta). La prise en charge de Lidar / Point Cloud est nettement améliorée. Vous pourriez probablement lire à partir de Lidar XYZ et écrire dans un format de grille sans utiliser de transformateurs, et cela fera du batch.
Mark Ireland
7

L'outil lasgrid.exe ( README ) de LAStools effectue une conversion "directe" avec maillage (sans triangulation) très rapidement. L'outil las2dem.exe ( README ) de LAStools rasterise les points via un TIN temporaire.

Martin Isenburg
la source
Merci, Martin. Avez-vous vu la réponse (de 9/2010) publiée par scw ou la référence à LAStools par wwnick (également le 9/2010)? Ils semblent couvrir exactement le même terrain, alors je me demande quelle distinction ou amélioration vous essayez de faire ici.
whuber
2
Je pense que le commentaire sur l'outil peut fonctionner sans triangulation clarifie une partie de la confusion.
mike
6

Vous pouvez essayer d'utiliser le pilote XYZ GDAL pour lire le format raster et le convertir en n'importe quel nombre de formats à l' aide gdal_translate. Vous avez besoin de GDAL 1.8 ou plus récent (par exemple, à partir d' OSGeo4W ).

Par exemple, pour traduire un fichier XYZ my_raster.xyzau format de grille Esri ASCII, utilisez la commande shell suivante (par exemple, à partir d'un shell OSGeo4W):

gdal_translate -of AAIGrid my_raster.xyz my_raster.asc
Mike T
la source
4

Je viens d'avoir ce même problème - le plus simple pour moi est de concaténer tous les fichiers en dehors du SIG.
Comme j'utilise Windows, j'ai utilisé la commande DOS

pour% f dans (* .xyz) tapez "% f" >> agrégat.xyz

puis je peux exécuter mes tâches de traitement normales (convertir de l'ASCII en grille en DEM ou autre) un seul pas 155 fois. Bien sûr, mon fichier texte résultant était volumineux (> 1 Go) mais ESRI pouvait le gérer.

dérose
la source
3

Jetez un œil à l' XYZ2DTMoutil de Fusion . Selon son manuel:

XYZ2DTM convertit les modèles de surface stockés sous forme de fichiers de points ACSII XYZ au format PLANS DTM. Les fichiers de points d'entrée incluent un enregistrement pour chaque point de grille avec les valeurs X, Y et d'élévation séparées par des virgules, des espaces ou des tabulations. En général, cet utilitaire n'est utilisé que lorsque les modèles de surface sont livrés dans ce format.

toto
la source
0

wwnick a mentionné LASTools , qui contient également une routine de conversion de nuage de points en DEM très rapide avec LAS2DEM . Si les performances sont importantes, je ferais un essai car je suppose que ce sera beaucoup plus rapide que certains des autres packages.

scw
la source
1
Merci. las2dem.exe triangule d'abord les points en TIN puis produit DEM. Je recherche une conversion directe.
Jakub Sisak GeoGraphics
-2

Je suggère d'essayer GlobalMapper pour cela (www.globalmapper.com). La version de démonstration peut être téléchargée à partir de leur page Web. Vous pouvez obtenir une licence de démonstration afin que vous puissiez essayer la fonctionnalité complète pendant une durée limitée si vous envoyez un e-mail au support.

Si vous voulez un outil gratuit pour manipuler les fichiers LIDAR, vous pouvez jeter un œil à l'outil de commande (LasTOOLS) http://www.cs.unc.edu/~isenburg/lastools/

STH
la source
Merci. Je vais essayer Global Mapper. En ce qui concerne les LasTools, je l'ai regardé il y a quelque temps et le seul que j'aurais pu utiliser était le las2dem.exe Cependant, cet outil triangule d'abord puis convertit en dem: "Cet outil lit les points LIDAR du format LAS, les triangule temporairement dans un TIN, puis raster le TIN sur un DEM. "
Jakub Sisak GeoGraphics
1
Je pense que votre réponse pourrait être améliorée en la modifiant pour inclure un lien ou une autre référence qui indique que GlobalMapper ne peut pas simplement utiliser des données LiDAR mais peut réellement exporter au format Esri ASCII Raster. Sinon, il semble que ce ne soit que spéculation.
PolyGeo
-2

SAGA GIS est gratuit et portable et peut gérer les données LIDAR. Je ne travaille pas avec ce type de données, mais SAGA prend en charge plusieurs formats à exporter (y compris le format de grille Esri ASCII), vous devriez donc essayer. Et il a une interface graphique relativement conviviale.

Juhele
la source
4
Je pense que votre réponse pourrait être améliorée en la modifiant pour inclure un lien ou une autre référence qui indique que SAGA ne peut pas simplement utiliser des données LiDAR mais peut réellement exporter au format Esri ASCII Raster. Sinon, il semble que ce ne soit que spéculation.
PolyGeo