Le portail francophone de la géomatique

tu as bien raison de poser des questions, c'est également comme ca qu'on apprend

je te propose donc de regarder le script que j'ai ecris pour importer automatiquement les dalles ASTER

for list in `find -type f  -iname "*_dem.tif" | sed 's/^\.\///'`
do
#folder=$(dirname $list);
file=$(basename $list);
name=$(basename $file _dem.tif);

echo `r.in.gdal -e --o input=${list} output=${name} title=${name}`;
echo `r.colors map=${name} color=srtm`;

done;

l'explication de texte:
- la boucle for liste de manière récursive le dossier courant à la rechercge de tous les fichiers contenant "_dem.tif"
- on recupere le nom du fichier (basename) duquel on supprimer la chaine _dem.tif pour obtenir le nom final
- on importe via r.in.gdal
- on colorise selon la palette srtm qui va est tres bien pour le relief

Bonjour,

Je suppose que le "echo" est là pour débuguer

Y.

Merci pour les explications supplémentaires, j'ai testé et ça marche. C'est d'ailleurs assez impressionnant. Je n'ai pas encore compris à quoi servent les symboles { }.

En tout cas ça m'a fait gagné du temps ; mais surtout j'ai pu mettre un premier pied dans la programmation en shell et ça c'est génial.
J'ai tenté de poursuivre en affichant tous les rasters importés avec la commande

for list in `g.list type=rast mapset=ASTER`
do
echo $list
echo `d.rast -o map=${list}`
done;

Mais l'affichage du raster c'est-à-dire la commande d.rast ne fonctionne pas dans le terminal on dirait. Après je dois peut-être créer un autre sujet si je veux continuer cette discussion.

En tout cas merci beaucoup pour le coup de main

Dernière modification par mygeomatic (Tue 01 December 2009 22:22)

Bonjour,

Bien merci pour l'explication.

Es-ce que je dois lancer un nouveau sujet pour la commande d.rast dans le terminale ou je peux continuer d'en parler ici ?

Voilà, quand j'exécute mon script pour afficher tous les rasters du mapset ASTER le terminal bloque à la fonction d.rast et me demande choisir un "monitor".
Seulement je ne sais s'il en existe un qui corresponde au Display x de ma fenêtre Layer Manager ou de ma fenêtre GRASS GIS Map Display x.  J'ai bien essayé avec un d.mon start=x1 , d.mon select=x1 ça lance bien une nouvelle fenêtre pour afficher les rasters mais je voudrais les voir dans le Layer Manager et le Map Display.

Voilà si quelqun connait la manip à faire.

Sinon je me permet encore une question : je fais ces manips chez moi sous ubuntu mais je me demande quel est le moyen de les faire sous windows (pour le travail). J'ai bien essayé d'utiliser le OSGeo Shell mais sans succé.

Bonne soirée

@mygeomatic :  Tu peux aussi émuler un système unix avec cygwin http://www.cygwin.com/ et travailler indifféremment en bash sous linux et windows

Bon je reviens sur la discussion pour dire que je me suis lancé dans l'apprentissage du Python.
En parallèle je test cygwin mais je n'arrive pas à lancer Grass depuis le Shell. Pour info j'ai installé GRASS avec OSGeo4W donc tous les dossier /bin /etc... sont dans C:\OSGeo4W.
J'ai regardé dans les fichiers grass64.bat et grass.bat pour essayer de comprendre comment le lancer mais je n'y arrive toujours pas.
Voilà je m'en remet encore à vos conseils et à votre gentillesse

Ok merci.

Je crois que je vais discrètement installer une machine virtuelle et ubuntu sur mon pc du boulot pour pouvoir travailler sous linux.

Bon week end