Suite

Fichier asc de rééchantillonnage Gdal/QGIS

Fichier asc de rééchantillonnage Gdal/QGIS


J'ai un fichier DEM (.asc) avec une résolution de 1000 mètres par pixel. Par conséquent, le fichier a beaucoup de colonnes et de lignes. Je voudrais hasarder la résolution par pixel de 1000 à 3000 mètres par pixel. J'aimerais donc le rééchantillonner.

quelle commande gdal dois-je utiliser pour cela? je sais que c'est gdalwarp mais comment est la ligne de commande exacte avec les options supplémentaires ?

si vous savez comment le faire dans QGIS, j'en suis également satisfait.


La tâche peut sembler triviale en lisant la documentation de gdalwarp http://www.gdal.org/gdalwarp.html et GDAL AAIGrid -- Arc/Info ASCII Grid driver documentation http://www.gdal.org/frmt_various.html. La taille du pixel cible est trois fois plus grande que la résolution native 0,008333333333 degrés/pixel (pas 1000 m/pixel, voir les commentaires).

gdalwarp -of AAIGrid -tr 0,024999 0,024999 input.asc output.asc

Cependant, c'est un peu plus difficile que cela.

  • gdalwarp ne prend pas en charge le format AAIGrid comme format de sortie directe
  • La méthode de rééchantillonnage par défaut de gdalwarp est le voisin le plus proche, ce qui ne convient pas bien au DEM

Par conséquent, la conversion doit être effectuée en deux étapes et avec une meilleure méthode de rééchantillonnage.

La première étape consiste à créer une sortie intermédiaire en tant que raster virtuel GDAL (.VRT) avec un rééchantillonnage "moyen"

gdalwarp -of VRT -r moyenne -tr 0,024999 0,024999 input.asc interim.vrt

La deuxième étape consiste à convertir le DEM provisoire en un nouveau fichier de grille ASCII avec gdal_translate

gdal_translate -of AAIGrid interim.vrt output.asc

Noter:

Les formats pris en charge pour la sortie peuvent être répertoriés avec la commande gdalwarp --formats. Les formats qui prennent en charge la sortie directe sont marqués du caractère "+". Cependant, tous ne fonctionneront pas vraiment (fichier NTv2 Datum Grid Shift par exemple).

Sous Windows

gdalwarp --formats|find "+" FITS -raster- (rw+): Système de transport d'images flexible HDF4Image -raster- (rw+): Jeu de données HDF4 netCDF -raster- (rw+s): Network Common Data Format VRT -raster- ( rw+v) : Virtual Raster GTiff -raster- (rw+vs) : GeoTIFF NITF -raster- (rw+vs) : National Imagery Transmission Format HFA -raster- (rw+v) : Erdas Imagine Images (.img) ELAS -raster- (rw+v) : ELAS MEM -raster- (rw+) : en mémoire Raster BMP -raster- (rw+v) : MS Windows Device Independent Bitmap PCIDSK -raster,vector- (rw+v) : PCIDSK Database Fichier ILWIS -raster- (rw+v): ILWIS Raster Map SGI -raster- (rw+): SGI Image File Format 1.0 Leveler -raster- (rw+): Leveler heightfield Terragen -raster- (rw+): Terragen heightfield ISIS2 -raster - (rw+v) : USGS Astrogeology ISIS cube (Version 2) ERS -raster- (rw+v) : ERMapper .ers Étiqueté RMF -raster- (rw+v) : Raster Matrix Format RST -raster- (rw+v ): Idrisi Raster A.1 INGR -raster- (rw+v): Intergraph Raster GSBG -raster- (rw+v): Golden Software Binary Grid (.grd) GS7BG -raster- (rw+v) : Golden Software 7 Binary Grid (.grd) PNM -raster- (rw+v) : Portable Pixmap Format (netpbm) ENVI -raster- (rw+v) : ENVI .hdr Étiqueté EHdr -raster- ( rw+v) : ESRI .hdr Labeled Paux -raster- (rw+) : PCI .aux Labeled MFF -raster- (rw+) : Vexcel MFF Raster MFF2 -raster- (rw+) : Vexcel MFF2 (HKV) Raster BT -raster- (rw+v) : VTP .bt (Terrain binaire) 1.3 Format LAN -raster- (rw+v) : Erdas .LAN/.GIS IDA -raster- (rw+v) : Image Data and Analysis GTX -raster- ( rw+v : NOAA Vertical Datum .GTX NTv2 -raster- (rw+vs): NTv2 Datum Grid Shift CTable2 -raster- (rw+v): CTable2 Datum Grid Shift KRO -raster- (rw+v): KOLOR Raw ADRG -raster- (rw+vs) : ARC Digitized Raster Graphics SAGA -raster- (rw+v) : SAGA GIS Binary Grid (.sdat) PDF -raster,vector- (rw+vs) : PDF géospatial

Voir la vidéo: QGIS Python PyQGIS - Get raster band data with GDAL