Annonce
L'appel à participation est ouvert jusqu'au 19 janvier 2025!
#1 Thu 06 March 2014 17:39
- BSarah
- Participant occasionnel
- Date d'inscription: 24 Feb 2014
- Messages: 18
GRASS r.sunmask
Bonjour à tous,
J'essaye de travailler avec le logiciel Grass pour modéliser des ombres portées à partir d'un MNS. J'utilise donc l'outil r.sunmask.
Or, après 2h de lancement de cet outil, la boite de dialogue affiche toujours: "calcul des ombres depuis le MNT". Il n'y a pas de barre d'état m'indiquant la progression du traitement. De plus, aucune carte n'a été créée...
Je souhaiterais donc savoir si il s'agit d'un bug ou si le calcul est long, tout simplement.
Je suis débutante sur ce logiciel. J'aimerai donc savoir si quelqu'un a déjà utilisé cet outil? Est-ce qu'il y a des pré-traitements nécessaire ou je peux l'utiliser directement à partir de mon MNS? Je n'ai pas trouver d'indications sur cela...
Merci d'avance,
Cordialement,
Sarah
Hors ligne
#2 Fri 07 March 2014 14:58
Re: GRASS r.sunmask
Bonjour r.sunmask est réputé pour être très lent.
r.sun avec l'option -s permet également de calculer les ombres portées. Il prend plus de paramètres en entrée mais est beaucoup plus performant.
Par ailleurs, il permet d'autres calculs - de radiation notamment. Je te laisse explorer la doc : http://grass.osgeo.org/grass64/manuals/r.sun.html
Ayant déjà utilisé à plusieurs reprises cet outil, je pourrai au besoin te donner un coup de main.
Bien à toi,
Simon
Hors ligne
#3 Mon 10 March 2014 10:11
- BSarah
- Participant occasionnel
- Date d'inscription: 24 Feb 2014
- Messages: 18
Re: GRASS r.sunmask
Bonjour,
Je viens d'essayer de tester r.sun, malheureusement deux erreurs me parviennent:
1/ G_malloc, problème de mémoire
2/ J'ai coché "économiser la mémoire" mais après survient d'horizon pré-calculés avant d'avoir de nouveau un problème de mémoire.
Sur la page d'information de Grass j'ai vu que quand le problème de mémoire se présente, il faut utiliser la numpartition. Je ne comprend pas à quoi cela correspond et comment cela fonctionne...
Aurais-tu une idée?
Géo amicalement,
Sarah
Hors ligne
#4 Mon 10 March 2014 10:44
Re: GRASS r.sunmask
Cela vient sans doute de la région sur laquelle vous travaillez :
- l'étendue de celle-ci correspond-elle bien à votre zone d'étude ou est-elle plus large ?
- la résolution adoptée n'est-elle pas trop fine ?
g.region est l'outil permettant d'avoir les informations sur la région GRASS.
geodata au cerema et petits billets en géomatique
Hors ligne
#5 Mon 10 March 2014 11:03
- BSarah
- Participant occasionnel
- Date d'inscription: 24 Feb 2014
- Messages: 18
Re: GRASS r.sunmask
Pour l'étendue de la région, j'ai vérifiée cela correspond bien à ma zone d'étude.
Par contre, la taille de mes cellule de mon MNT es de 0.25 x 0.25m. Pensez-vous que ma résolution est trop fine?
Pour créer mon MNT, j'ai des données Lidar à la base. Je souhaite donc avoir une résolution fine pour garder la précision de mes données d'origines.
Hors ligne
#6 Mon 10 March 2014 14:00
Re: GRASS r.sunmask
Pour l'étendue de la région, j'ai vérifiée cela correspond bien à ma zone d'étude.
Par contre, la taille de mes cellule de mon MNT es de 0.25 x 0.25m. Pensez-vous que ma résolution est trop fine?
c'est fin, mais peut-être nécessaire à tes besoins (toi seule, le sais). Quelle est l'emprise ?
L'option numpartitions est utile mais inopérante pour le calcul des ombres (option -s).
Quelle version de GRASS utilises-tu?
Un peu de lecture en attendant : http://grasswiki.osgeo.org/wiki/Large_r … processing
Hors ligne
#7 Sun 13 September 2015 20:25
- ploicq
- Participant occasionnel
- Lieu: Tours
- Date d'inscription: 30 Sep 2009
- Messages: 26
Re: GRASS r.sunmask
Bonjour,
Je me permets de déterrer ce topic car je suis sur la même problématique. J'aimerais optimiser le temps de calcul avec r.sun sur Grass 7 pour déterminer les zones ombragées. Actuellement je fais mes tests sur les données échantillons de Caroline du Nord (données "elevation"):
Si j'applique directement r.sun:
r.sun elevation=elevation day=88 time=7 incidout=incidoutD88T07 --overwrite
ça prend 28 secondes
Si je crée des horizons tous les 30° comme conseillé dans le manuel:
r.horizon elevation=elevation@PERMANENT step=30 maxdistance=5000 output=horizon (pour info cette commande met 17min sur un bi-coeur)
puis que je les mets en input:
r.sun elevation=elevation horizon_basename=horizon horizon_step=30 day=88 time=7 incidout=incidoutD88T07 --overwrite
ça prend 51 secondes
Or il est mis partout que le pré-calcul des horizons est censé accélérer le calcul.
En fait je comprend pas (c'est noté nulle part) comment il fait le join entre l'azimuth solaire du 88ème jour à 7h00 et celui que j'ai calculé avant ? Prend-il le plus proche ?
Merci pour votre aide
Pierre de Tours
Dernière modification par ploicq (Sun 13 September 2015 20:26)
Hors ligne
#8 Sat 19 September 2015 16:03
Re: GRASS r.sunmask
Bonjour Pierre,
Pour ce type de question, je t'invite à consulter directement la liste GRASS sur laquelle tu trouveras les contributeurs du module r.sun : http://osgeo-org.1560.x6.nabble.com/Gra … 84509.html
Tu peux aussi fouiller les archives : http://osgeo-org.1560.x6.nabble.com/tem … ;n=3681316
Merci de nous tenir informé
simo
Hors ligne
#9 Mon 21 September 2015 11:27
- ploicq
- Participant occasionnel
- Lieu: Tours
- Date d'inscription: 30 Sep 2009
- Messages: 26
Re: GRASS r.sunmask
Merci de ton conseil,
Je viens de poster, voici le lien: http://osgeo-org.1560.x6.nabble.com/Spe … 25165.html
J'ai aussi posté sur StackExchange mais pas de réponses pour le moment: http://gis.stackexchange.com/questions/ … -gis-r-sun
Hors ligne
#10 Mon 21 September 2015 21:08
Re: GRASS r.sunmask
Je te renvoie, en passant, vers un message que j'avais moi-même posté il y a quelques années : http://osgeo-org.1560.x6.nabble.com/r-s … 7682.html. Peut-être y trouveras -tu des infos utiles ... A l'époque, je travaillais avec la v6.4
Potentiellement, si j'ai un jeu de données à portée de main et un peu de temps, je ferai également un test sur les temps de traitement avec la v7.
Hors ligne
#11 Tue 22 September 2015 09:21
- ploicq
- Participant occasionnel
- Lieu: Tours
- Date d'inscription: 30 Sep 2009
- Messages: 26
Re: GRASS r.sunmask
Merci. Ton lien ne marche pas chez moi mais j'ai retrouvé le sujet ici: http://osdir.com/ml/grass-user-gis/2013 … 00313.html
Je devrais pouvoir faire tourner la GRASS 7 sur linux bientôt, sur 16 coeurs, avec l'extension r.sun.hourly (r.sun.daily existe aussi)
Hors ligne
