Le portail francophone de la géomatique

Bonjour,

> Je suis moi aussi a la recherche de solution pour exploiter les fichiers
srtm sour arcview 3, vous est il possible de donner un fichier header
exemple (hdr) a partir d'un srtm particulier (ex. N02E096.hgt)?

BYTEORDER M
LAYOUT BIL
NROWS 1201
NCOLS 1201
NBANDS 1
NBITS 16
BANDROWBYTES 2402
TOTALROWBYTES 2402
BANDGAPBYTES 0
NODATA -32768
ULXMAP 95.999583333333333
ULYMAP 3.0004166666666667
XDIM 0.000833333333333
YDIM 0.000833333333333

Sauf erreur de ma part.
Pour le lire, il vous faut changer les extensions (*.hgt en *.bil) et garder
le nom original des fichiers.

Ci-dessous, un bout de code pour automatiser la generation des headers sous
Arcview 3.x.

Cordialement,
Sylvain

--------

' -----
-------
' Nom .........: ecriture_header_srtm
' Version .....: 1
' Auteur ......: Sylvain Thery (thery@texte-a-enlever.ccr.jussieu.fr)
' Date .........: avril 2004

' -----
-------
' Self ........: nil
' Return ......: nil
' -----
-------
' Commentaires : genere les headers des fichiers SRTM (*.hgt renommes en
*.bil)

' -----
-------
' Completez la ligne ci-dessous avec vos fichiers....
' vous pouvez en mettre autant que vous le souhaitez.

myList = { d:N02E096.bil ,  d:N02E097.bil , ...  })

for each rec in myList
myHDR_f = rec.substitute( .bil ,  .hdr )
myHDR = myHDR_f.substitute( e ,  E )
myHDR = myHDR_f.substitute( n ,  N )
myHDR = myHDR_f.substitute( s ,  S )
myHDR = myHDR_f.substitute( w ,  W )
myHDR_f = myHDR_f.asfilename
myName = rec.asFilename.getbasename.substitute( .bil ,   )
myLat = myName.left(1)
myLat_v = (myName.left(3)).right(2).asNumber
myLong = (myName.right(4)).left(1)
myLong_v = myName.right(3).asNumber

' ****************
' l'instruction suivante est en une seule ligne !

myHDR_s =  BYTEORDER M +NL+ LAYOUT BIL +NL+ NROWS 1201 +NL+ NCOLS
1201 +NL+ NBANDS 1 +NL+ NBITS 16 +NL+ BANDROWBYTES 2402 +NL+ TOTALROWBYTES
2402 +NL+ BANDGAPBYTES 0 +NL+ NODATA -32768 +NL+ ULXMAP

' l'instruction precedente est en une seule ligne !
' ****************

if ((myLong =  e ) or (myLong =  E )) then
myVal = myLong_v-0.0004166666666667
myval.setformat( d.dddddddddddddddd )
else
myVal = 0-(myLong_v)+0.0004166666666667
myval.setformat( d.dddddddddddddddd )
end
myHDR_s = myHDR_s + myval.asString+nl+ ULYMAP

if ((myLat =  n ) or (myLat =  N )) then
myval = myLat_v + 1.0004166666666667
myval.setformat( d.dddddddddddddddd )
else
myval = 0-(myLat_v) - 1.0004166666666667
myval.setformat( d.dddddddddddddddd )
end
myHDR_s = myHDR_s + myval.asString +NL+ XDIM 0.000833333333333 +NL+ YDIM
0.000833333333333

lf = LineFile.Make(myHDR_f, #FILE_PERM_WRITE )
lf.WriteElt(myHDR_s)
lf.Close
end

' -----
-------
' EoS
' -----
-------

De: Raphael Kourdian
Bonjour,
Il n'est pas possible de lire ce format .hgt d'image brute (certainement Aster) a partir d'ArcView 3.1. Cependant des logiciels de teledetection permettent l'importation de ce fichier et son exportation en fichier image lisible (tif, jpg, etc) : TNT (Microimages) ou Erdas Imagine (Leica) permettent cette transformation.
En esperant avoir pu vous aider.

Cordialement,

Raphael Kourdian

De: Ihab Hawater
Bonjour,

En effet 'Grid Translator Pro for MapInfo' lit les donnes .hgt, mais si vous avez ArcGIS vous pouvez utiliser 'Grid Translator Pro for ArcGIS'.
http://www.geomaticssystems.com/Products/GTPArcGIS.htm

Une version gratuite temporaire de quelque jours peut etre telechargee a:
http://www.geomaticssystems.com/Product … STryIt.htm

Ihab