Extraire données du résultat d'une commande

Version imprimable

Bonjour,
J'ai besoin d'extraire une donnée précise de ... 524 288 fichiers 8O

J'utilise la commande suivante :

Code:

1
2
 
gdalinfo MyFile.tiff

Ce qui me donne

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
 
Driver: GTiff/GeoTIFF
Files: map_11_0_0.tiff
Size is 256, 256
Coordinate System is:
GEOGCRS["WGS 84",
    DATUM["World Geodetic System 1984",
        ELLIPSOID["WGS 84",6378137,298.257223563,
            LENGTHUNIT["metre",1]]],
    PRIMEM["Greenwich",0,
        ANGLEUNIT["degree",0.0174532925199433]],
    CS[ellipsoidal,2],
        AXIS["geodetic latitude (Lat)",north,
            ORDER[1],
            ANGLEUNIT["degree",0.0174532925199433]],
        AXIS["geodetic longitude (Lon)",east,
            ORDER[2],
            ANGLEUNIT["degree",0.0174532925199433]],
    ID["EPSG",4326]]
Data axis to CRS axis mapping: 2,1
Origin = (-180.000000000000000,90.000000000000000)
Pixel Size = (0.001373291015625,-0.001373291015625)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (-180.0000000,  90.0000000) (180d 0' 0.00"W, 90d 0' 0.00"N)
Lower Left  (-180.0000000,  89.6484375) (180d 0' 0.00"W, 89d38'54.38"N)
Upper Right (-179.6484375,  90.0000000) (179d38'54.38"W, 90d 0' 0.00"N)
Lower Right (-179.6484375,  89.6484375) (179d38'54.38"W, 89d38'54.38"N)
Center      (-179.8242188,  89.8242188) (179d49'27.19"W, 89d49'27.19"N)
Band 1 Block=256x256 Type=Byte, ColorInterp=Palette
  Image Structure Metadata:
    NBITS=1
  Color Table (RGB with 2 entries)
    0: 0,0,0,255
    1: 255,255,255,255

j'ai besoin d'extraire Upper Left plus précisément -180.0000000, 90.0000000 et de mettre ça dans un fichier global (pour la totalité des 524 288 fichiers) ou une ligne serait du type
MyFile => -180.0000000, 90.0000000Comment je pourrais faire ça ?

Merci de votre aide ;)

Bonjour :coucou:

Code:

1
2
$ gdalinfo MyFile.tiff | sed '/^Upper Left/!d;s/[^(]*(\([^)]*\)).*/\1/' 
-180.0000000,  90.0000000

18/11/2021, 08h12
rlelamer

Merci beaucoup ça fonctionne parfaitement ;)

Je vais tenter de modifier la sortie avec sed pour obtenir une ligne comme ceci MyFile.tiff => -180.0000000, 90.0000000 avec retour à la ligne
Pas évident pour un non initié cette syntaxe mais man sed devrait être mon ami (ainsi que les exemples / explications sur le net)
18/11/2021, 16h52
jack-ft
Citation:
Envoyé par rlelamer

Bonjour,
J'ai besoin d'extraire une donnée précise de ... 524 288 fichiers 8O

J'utilise la commande suivante :

Code:

1 2 gdalinfo MyFile.tiff

Ce qui me donne

Code:

1 2 3 4 Corner Coordinates: Upper Left (-180.0000000, 90.0000000) (180d 0' 0.00"W, 90d 0' 0.00"N) Lower Left (-180.0000000, 89.6484375) (180d 0' 0.00"W, 89d38'54.38"N)

j'ai besoin d'extraire Upper Left plus précisément -180.0000000, 90.0000000 et de mettre ça dans un fichier global (pour la totalité des 524 288 fichiers) ou une ligne serait du type
MyFile => -180.0000000, 90.0000000
Citation:
Envoyé par Flodelarab

Bonjour :coucou:

Code:

1 2 $ gdalinfo MyFile.tiff | sed '/^Upper Left/!d;s/[^(]*($[^)]*$).*/\1/' -180.0000000, 90.0000000
Du coup, une 'tite boucle:

Code:

for f in *.tiff; do echo "$f => $(gdalinfo "$f" | sed '/^Upper Left/!d;s/[^(]*($[^)]*$).*/\1/')"; done > le_fichier_global.txt

à supposer que *.tiff ramène la liste des fichiers qui sont tous dans le même dossier.
Sinon y a find...
18/11/2021, 17h53
rlelamer

Citation:

Envoyé par jack-ft

Du coup, une 'tite boucle:

Code:

for f in *.tiff; do echo "$f => $(gdalinfo "$f" | sed '/^Upper Left/!d;s/[^(]*($[^)]*$).*/\1/')"; done > le_fichier_global.txt

à supposer que *.tiff ramène la liste des fichiers qui sont tous dans le même dossier.
Sinon y a find...

Pourquoi chercher compliqué quand on peut faire beaucoup plus simple ...
Moi qui comptais m'exciter les neurones pour faire ça qu'avec sed ...

Merci beaucoup a vous 2, vous venez de me faire gagner un temps ... monstrueux :D ;)

EDIT:
J'ai une erreur que je comprends pas très bien d'ailleurs puisque la commande "seule" fonctionne correctement

Code:

-bash: !d: event not found
18/11/2021, 18h12
N_BaH

gdalinfo n'accepte qu'un seul fichier ? il ne peux pas prendre une liste ?

Code:

gdalinfo *.tiff

?

Code:

1
2
3
4
5
6
7
8
9
10
11
12
 
MacBook-Air-de-Romain: ~/VRTiles/tiff $ gdalinfo *.tiff
-bash: /opt/homebrew/bin/gdalinfo: Argument list too long
 
MacBook-Air-de-Romain: ~/VRTiles/tiff $ gdalinfo map_11_0_0.tiff map_11_0_1.tiff
ERROR 6: Too many command options 'map_11_0_1.tiff'
Usage: gdalinfo [--help-general] [-json] [-mm] [-stats | -approx_stats] [-hist] [-nogcp] [-nomd]
                [-norat] [-noct] [-nofl] [-checksum] [-proj4]
                [-listmdd] [-mdd domain|`all`] [-wkt_format WKT1|WKT2|...]*
                [-sd subdataset] [-oo NAME=VALUE]* [-if format]* datasetname
 
MacBook-Air-de-Romain: ~/VRTiles/tiff $

Bah flute alors ...

18/11/2021, 19h40
N_BaH
Citation:
Envoyé par rlelamer

Code:

1 2 MacBook-Air-de-Romain: ~/VRTiles/tiff $ gdalinfo *.tiff

Bah flute alors ...
j'avais zappé "524 288 fichiers" ! 8O :calim2:

Alors je ne sais pas pourquoi mais si je le fais comme ceci
Code:

1 2 3 4 5 6 7 #!/bin/bash for f in *.tiff; do echo "$f => $(gdalinfo "$f" | sed '/^Upper Left/!d;s/[^(]*($[^)]*$).*/\1/')"; done > PosGPS.txt
et que je lance ce script, bah ça fonctionne correctement 🤷

18/11/2021, 20h33
N_BaH

le ! dans le terminal appelle le Développement d'historique.
or, par défaut, il n'y a pas d'historique dans un script.