Suite

Ouvrez le DEM LiDAR du Royaume-Uni avec l'extension *.asc

Ouvrez le DEM LiDAR du Royaume-Uni avec l'extension *.asc


J'ai téléchargé une vignette d'un accès libreDEM LiDARde laRoyaume-Unià partir de http://environment.data.gov.uk/ds/survey#/download.

Mais je ne sais pas comment l'ouvrir. Les fichiers à l'intérieur ont le*.ascextension. Je peux voir que cette extension ne fait pas partie de celles prises en charge parGdal. Et, ni la plate-forme d'application SentinelSE CASSERniENVIréussi à l'ouvrir.

Quel logiciel est utilisé pour ouvrir cesDEMdes dossiers?


Il s'agit des fichiers de grille ESRI ASCII :

$ head -6 sd5755_DSM_2m.asc ncols 500 nrows 500 xllcorner 357000 yllcorner 455000 cellssize 2 NODATA_value -9999

avec les coordonnées dans les références de grille OSGB (code EPSG 27700). C'est à peu près la seule chose que GDAL ne sache pas à leur sujet :

$ gdalinfo sd5755_DSM_2m.asc Pilote : AAIGrid/Arc/Info Fichiers de grille ASCII : sd5755_DSM_2m.asc La taille est de 500, 500 Le système de coordonnées est "Origine = (357000.000000000000000,456000.0000000000000000000) Taille de pixel = (2.000000000000000,-2.000000000000000) Coordonnées du coin supérieur gauche 357000.000, 456000.000) En bas à gauche ( 357000.000, 455000.000) En haut à droite ( 358000.000, 456000.000) En bas à droite ( 358000.000, 455000.000) Centre ( 357500.000, 455500.000) Band 1 Block=500x1 Type=Float32, ColorInterp=Undefined NoData Value=-9999

vous devrez donc informer votre SIG du CRS si vous avez besoin de le projeter sur autre chose.


GDAL peut l'ouvrir pour vous, il doit s'agir d'un fichier ASCII. Vous pouvez facilement traduire en tif avec

gdal_translate source.asc dest.tif


J'opterais pour une approche en deux étapes pour créer un seul fichier à partir de tous les fichiers au format ASCII téléchargés à l'aide de deux utilitaires GDAL. Depuis le fichier décompressé :

  1. Créez un 'Virtual Raster' avec tous les fichiers :

    gdalbuildvrt merged_mosaic.vrt *asc

Je ne sais pas comment l'extension shell fonctionnera si vous utilisez Windows - vous pouvez également transmettre une liste de fichiers, voir http://www.gdal.org/gdalbuildvrt.html pour plus de détails.

  1. Convertissez le VRT en GeoTIFF (ou tout autre format pris en charge par GDAL) et attribuez la projection OSGB :

    gdal_translate -of GTiff -a_srs EPSG:27700 merged_mosaic.vrt merged_mosaic.tif

Si cela est utile, j'ai également écrit un script qui mosaïquera les fichiers à partir du fichier zip. Il utilise RSGISLib (http://rsgislib.org/) mais pourrait être adapté pour utiliser à la place gdalbuildvrt / gdal_translate. Voir : https://spectraldifferences.wordpress.com/2015/09/02/mosaic-environment-agency-dtmdsm-tiles-using-rsgislib/ pour plus de détails.

Mise à jour : j'ai également écrit un script qui fera des mosaïques en utilisant GRASS : https://gist.github.com/arsf/fcf51a0dd8b7dcde8243


Voir la vidéo: La frontière entre la France et le Royaume-Uni fermée pour 48 heures