Suite

Utilisation de Gdalwarp pour modifier la résolution du fichier en degrés en mètres


J'utilise la commande : gdalwarp -t_srs WGS84 -te 144 -36 150 -32 -crop_to_cutline AustTopo250_SI55.tif

Je dois changer la résolution des pixels à 50 m x 50 m mais l'option est donnée sous la forme -tr xres yres : définir la résolution du fichier de sortie (en unités géoréférencées cibles).

L'unité cible est le degré. La question est donc : comment générer le fichier de sortie en pixels de 50 m x 50 m ?

N'hésitez pas à expliquer la relation entre les mètres et les degrés.

Fichier cible:

E:>gdalinfo AustTopo250_SI55.tif Pilote : GTiff/GeoTIFF Fichiers : AustTopo250_SI55.tif La taille est 18750, le système de coordonnées 12500 est : GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563 , AUTORITE["EPSG","7030"]], AUTORITE["EPSG","6326"]], PRIMEM["Greenwich",0], UNITÉ["degré",0.0174532925199433], AUTORITE["EPSG"," 4326"]] Origine = (144.00000000000000000,-32.0000000000000000) Taille de pixel = (0.0003200000000000,-0.0003200000000000) Métadonnées : AREA_OR_POINT=Zone COLORSPACE=RGB COMPRESSION_RATE_TARGET=9 VERSION=2 Métadonnées de structure d'image : INTERLEAVE=PIXEL Coordonnées d'angle : en haut à gauche (144.0000000, -32.0000000) (144d 0' 0.00"E, 32d 0' 0.00"S) En bas à gauche ( 144.0000000, -36.0000000) (144d 0' 0.00"E, 36d 0' 0.00"S) En haut à droite ( 150.0000000, -32.0000000) ( 150d 0' 0.00"E, 32d 0' 0.00"S) En bas à droite ( 150.0000000, -36.0000000) (150d 0' 0.00"E, 36d 0' 0.00"S) Centre ( 147.0000000, -34.0000000) (147d 0' 0.00" E, 34d 0' 0.00"S) Bande 1 Bloc=18750x1 Type=Byte, ColorInterp=Rouge Description = Bloc Bande Rouge 2=18750x1 Type=Byte, ColorInterp=Vert Description = Bloc Vert Bande 3=18750x1 Type=Byte, ColorInterp=Bleu Description = Bleu

Fichier source:

E:>gdalinfo 76943_Aust_gda94.ecw Pilote : ECW/ERDAS Compressed Wavelets (SDK 3.x) Fichiers : 76943_Aust_gda94.ecw La taille est 128490, 109375 Le système de coordonnées est : GEOGCS["GEOCENTRIC DATUM of AUSTRALIA", DATUM["GDA94", SPHEROID["GRS80",6378137,298.257222101]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433]] Origine = (112.883360000000040,-8.9999200000000039) Taille de pixel = (0.0003200000000000,-0.0003200000000000) Métadonnées : COLORSPACE =RGB COMPRESSION_RATE_TARGET=9 VERSION=2 Coordonnées du coin : en haut à gauche ( 112.8833600, -8.9999200) (112d53' 0.10"E, 8d59'59.71"S) en bas à gauche ( 112.8833600, -43.9999200) (112d53' 0.10"E, 43d59'59.71 "S) En haut à droite ( 154.0001600, -8.9999200) (154d 0' 0.58"E, 8d59'59.71"S) En bas à droite ( 154.0001600, -43.9999200) (154d 0' 0.58"E, 43d59'59.71"S) Centre ( 133.4417600 , -26.4999200) (133d26'30.34"E, 26d29'59.71"S) Bande 1 Bloc=256x256 Type=Byte, ColorInterp=Rouge Description = Rouge Aperçus : 64245x54687, 32122x27343, 16061x13671, 8030x6835, 4015x3417, 2007x1708, 1003x854, 501x427, 250x213 Band 2 Block=256x256 Type=Byte, ColorInterp=Green Description = Green Aperçus : 64245x54687, 32122x27343, 16061x13671, 8030x6835, 4015x3417, 2007x1708, 1003x854, 501x427, 250x213 Band 3 Block=256x256 Color Type=Byblue Description = Aperçus bleus : 64245x54687, 32122x27343, 16061x13671, 8030x6835, 4015x3417, 2007x1708, 1003x854, 501x427, 250x213

Si vous avez besoin d'un raster avec une résolution fixe en mètres, vous devez reprojeter vos données dans un CRS projeté, comme UTM, Lambert conforme conique ou ce qui pourrait être utilisé en Australie.

Alors que les degrés de latitude ont plus ou moins la même taille en mètres, un degré de longitude diminue de l'équateur vers les pôles, où il est presque nul. Il est donc peu logique de substituer une unité de degré par une unité de mètre sans reprojeter correctement.


Voir la vidéo: Teams: Kansion luominen, tiedoston siirtäminen kansioon, tiedoston linkittäminen keskusteluun (Octobre 2021).