Annonce
#1 Fri 10 April 2015 15:50
QGIS 2.8 / GDAL2XYZ : probleme sur un GRD multi-bande
Bonjour,
Je soumets à vos lumière le fichier joint à ce message. A première vue, c'est un fichier GRD (Surfer Grid) basique qui contient des valeurs altimétriques de type MNT.
Le hic, c'est qu'au lieu d'être mono-bande comme on s'attendrait à ce qu'il le soit, il est... multi-bande (4 en tout).
Je voudrais le convertir en mono-bande et retrouver mes altimétries, mais je ne vois pas trop comment faire sur ce coup.
Pour info :
* Global Mapper le lit sans problème, l'interprète, et je ne sais comment me donne une altimétrie (entre 204 et 219).
* Les valeurs ne sont pas sur la bande 4
* RGB2PCT en convient pas > ça ne donne pas de Z...
Des idées ?
Pour complément
FILENAME=Aconvertir.grd
DESCRIPTION=Aconvertir.grd
UPPER LEFT X=747000.500
UPPER LEFT Y=6621999.500
LOWER RIGHT X=747999.500
LOWER RIGHT Y=6621000.500
WEST LONGITUDE=3° 36' 53.9489" E
NORTH LATITUDE=46° 41' 47.1928" N
EAST LONGITUDE=3° 37' 41.3789" E
SOUTH LATITUDE=46° 41' 14.5564" N
PROJ_DESC=Lambert-93 / GRS80 / meters
PROJ_DATUM=GRS80 ELLIPSOID
PROJ_UNITS=meters
COVERED AREA=0.998 sq km
NUM COLUMNS=1000
NUM ROWS=1000
NUM BANDS=1
PIXEL WIDTH=1 meters
PIXEL HEIGHT=1 meters
MIN ELEVATION=204.19 meters
MAX ELEVATION=219.67 meters
ELEVATION UNITS=meters
BIT DEPTH=16
SAMPLE TYPE=Unsigned 16-bit Integer
gdalinfo Aconvertir.grd
Driver: NWT_GRD/Northwood Numeric Grid Format .grd/.tab
Files: Aconvertir.grd
Size is 1000, 1000
Coordinate System is:
PROJCS["unnamed",
GEOGCS["unnamed",
DATUM["GRS_80",
SPHEROID["GRS 80",6378137,298.257222101],
TOWGS84[0,0,0,-0,-0,-0,0]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",44],
PARAMETER["standard_parallel_2",49],
PARAMETER["latitude_of_origin",46.5],
PARAMETER["central_meridian",3],
PARAMETER["false_easting",700000],
PARAMETER["false_northing",6600000],
UNIT["Meter",1]]
Origin = (747000.000000000000000,6622000.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Corner Coordinates:
Upper Left ( 747000.000, 6622000.000) ( 3d36'54.29"E, 46d41'47.21"N)
Lower Left ( 747000.000, 6621000.000) ( 3d36'53.93"E, 46d41'14.80"N)
Upper Right ( 748000.000, 6622000.000) ( 3d37'41.40"E, 46d41'46.95"N)
Lower Right ( 748000.000, 6621000.000) ( 3d37'41.03"E, 46d41'14.54"N)
Center ( 747500.000, 6621500.000) ( 3d37'17.66"E, 46d41'30.88"N)
Band 1 Block=1000x1 Type=Byte, ColorInterp=Red
Band 2 Block=1000x1 Type=Byte, ColorInterp=Green
Band 3 Block=1000x1 Type=Byte, ColorInterp=Blue
Band 4 Block=1000x1 Type=Float32, ColorInterp=Undefined
NoData Value=-9.9999999338158125e+036
Offset: 204.190002441406, Scale:0.000236213198149648
Dernière modification par Robin (Thu 16 April 2015 15:49)
Hors ligne
#2 Fri 10 April 2015 19:32
Re: QGIS 2.8 / GDAL2XYZ : probleme sur un GRD multi-bande
Salut Robin,
Avec gdal_translate et le l'option -b, peut-être pour extraire la bande souhaitée..
geodata au cerema et petits billets en géomatique
Hors ligne
#3 Wed 15 April 2015 09:38
Re: QGIS 2.8 / GDAL2XYZ : probleme sur un GRD multi-bande
Merci de ta réponse. Le problème c'est qu'-e l'extract de bande avec -b .... m'extrait une bande uniquement.
Je n'ai pas résolu le mystère mais Global Mapper m'a sauvé la mise. Suite à mon message et pour complément, on a découvert que Vertical Mapper aussi lisait bien l'altimétrie. Je suppose qu'il doit y avoir un entête de fichier que gdal ne comprend pas (?)
Dernière modification par Robin (Wed 15 April 2015 09:38)
Hors ligne
#4 Wed 15 April 2015 11:27
Re: QGIS 2.8 / GDAL2XYZ : probleme sur un GRD multi-bande
oups, désolé, je pensais que tu voulais extraire une des 4 bandes alros qu'en fait, tu souhaitais une sorte de combinaison des 4 en une ?
geodata au cerema et petits billets en géomatique
Hors ligne
#5 Thu 16 April 2015 15:47
Re: QGIS 2.8 / GDAL2XYZ : probleme sur un GRD multi-bande
Oui ou plutôt, je suppose que c'est ce qu'il faudrait que je fasse pour pouvoir extraire mes Z via gdal2xyz. En fait, au fond de cette affaire, je voudrais surtout comprendre pourquoi QGIS ne comprend pas mon GRD aussi bien que Vertical Mapper ou Global Mapper 
Dernière modification par Robin (Thu 16 April 2015 15:48)
Hors ligne
