Annonce
Pages: 1
- Sujet précédent - Passage de données au format jpeg2000 vers le format DTED1 - Sujet suivant
#1 Tue 30 March 2010 18:34
- smillet
- Juste Inscrit !
- Date d'inscription: 30 Mar 2010
- Messages: 1
Passage de données au format jpeg2000 vers le format DTED1
Bonjour,
Etant novice dans tout ce qui touche à la manipulation de fichiers de données géographiques, j'aurais besoin d'aide pour convertir des données en JP2 vers du DTED1. J'ai trouvé gdal_translate et gdalwarp mais suis un peu perdu qd à leurs utilisations, le jpeg2000 en entrée a les caractéristiques suivantes (gdalinfo):
Driver: JP2ECW/ERMapper JPEG2000
Files: srtm.jp2
Size is 119999, 74401
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (-31.999583000000008,60.000416333276860)
Pixel Size = (0.000833333333333,-0.000833333333333)
Corner Coordinates:
Upper Left ( -31.9995830, 60.0004163) ( 31d59'58.50"W, 60d 0'1.50"N)
Lower Left ( -31.9995830, -2.0004170) ( 31d59'58.50"W, 2d 0'1.50"S)
Upper Right ( 67.9995837, 60.0004163) ( 67d59'58.50"E, 60d 0'1.50"N)
Lower Right ( 67.9995837, -2.0004170) ( 67d59'58.50"E, 2d 0'1.50"S)
Center ( 18.0000003, 28.9999997) ( 18d 0'0.00"E, 29d 0'0.00"N)
Band 1 Block=119999x1 Type=Int16, ColorInterp=Gray
Overviews: arbitrary
Merci d'avance...
Hors ligne
#2 Wed 31 March 2010 15:24
- mygeomatic
- Participant occasionnel
- Date d'inscription: 14 Oct 2009
- Messages: 24
- Site web
Re: Passage de données au format jpeg2000 vers le format DTED1
Bonjour,
Peux-tu me dire sous quel système tu travail s'il te plait ?
Il s'agit de changer le format des données donc normalement gdalwarp devrait faire l'affaire. Après il me semble que le format jpeg2000 n'est pas systématiquement implémenté (je n'ai jamais eu l'occasion de tester), je pense qu'avec FWTools ça ne devrait pas poser de problèmes. Il va donc falloir travailler en ligne de commande. La commande devrait ressembler à quelque chose comme ça
Code:
gdalwarp -of "DTED" "adresse complète du fichier jp2" "adresse complète du fichier en sortie"
Après tu devrais trouver quelque ressources sur le forum pour utiliser les utilitaires gdal : http://georezo.net/wiki/main:logiciels:gdal_ogr
Et je te conseil également la page du projet gdal bien sûr : http://www.gdal.org/index.html
Cordialement,
Charles
Hors ligne
#3 Wed 31 March 2010 22:52
- rouault
- Participant assidu
- Date d'inscription: 26 Apr 2009
- Messages: 166
Re: Passage de données au format jpeg2000 vers le format DTED1
Il faut en fait découper ton JP2 en autant de tuiles de 1 degré x 1 degré de taille 1201 x 1201 pixels (ou 601 x 1201 pixels pour lat >= 50).
Installe Python et essaye le script suivant qui met en oeuvre gdalwarp et gdal_translate :
Code:
import os
src_file = 'srtm.jp2'
for x in range(-32,68):
for y in range(-2,60):
if y >= 50:
xres = 2. / 1200.
else:
xres = 1. / 1200.
yres = 1. / 1200.
xmin = x - xres / 2
ymin = y - yres / 2
xmax = x + 1 + xres / 2
ymax = y + 1 + yres / 2
try:
os.remove('dst.vrt')
except:
pass
cmd = 'gdalwarp -of VRT %s dst.vrt -tr %.16f %.16f -te %.16f %.16f %.16f %.16f' % (src_file, xres, yres, xmin, ymin, xmax, ymax)
print cmd
os.popen(cmd).read()
if y < 0:
ns = 'S'
else:
ns = 'N'
if x < 0:
ew = 'W'
else:
ew = 'E'
cmd = 'gdal_translate -of DTED dst.vrt %c%02d%c%03d.dt1' % (ns, abs(y), ew, abs(x))
print cmd
os.popen(cmd).read()Pour un autre fichier d'entrée, il suffirait d'adapter les range() d'itération de x et y à l'étendue spatiale du fichier d'entrée.
Hors ligne
Pages: 1
- Sujet précédent - Passage de données au format jpeg2000 vers le format DTED1 - Sujet suivant
