Le portail francophone de la géomatique

MathieuR · Tue 16 February 2010 15:50

Bonjour,

Je tente la vectorisation d'une image cadastrale.
Je bloque à l'étape de vectorisation où il faut affiner le trait.

Voici mes différentes étapes:

Import de l'image dans GRASS

Code:

r.in.gdal in=C:/tests/GRASS/BDP_08_1234_6134_010_L93.tif out=raster

Datum <Reseau_Geodesique_Francais_1993> not recognised by GRASS and no parameters found
Projection of input dataset and current location appear to match
r.in.gdal complete. Raster map <raster2> created.

Une fois cela fait, je lance la commande r.thin:

Code:

r.thin in=raster out=rasterThin

Mais la fenêtre de sortie Log m'affiche le message suivant:

Code:

File raster2 -- 131957 rows X 108701 columns
ERROR: (null): Error writing temporary file
(Tue Feb 16 16:04:59 2010) Command finished (1057 sec)

Qu'est-ce qui peut bloquer?
un manque d'espace sur le disque dur pour l'écriture du fichier temporaire?
Faut-il alléger le fichier source?
Comment allouer plus d'espace à GRASS pour le stockage des fichiers temporaires?

La vectorisation de la couche "routes" du jeu de données spearfish s'est effectuée, elle, avec succès..

En vous remerciant par avance de votre aide

Dernière modification par baobazz (Tue 16 February 2010 16:13)

Sylvain_M · Tue 16 February 2010 16:35

essaies de relancer la commande r.thin en mode bavard (option --v) ... sinon tu peux aussi augmenter le niveau de verbosité de grass (DEBUG=1à5 dans le fichier .grassrc6)

Dernière modification par yarma (Tue 16 February 2010 16:36)

MathieuR · Tue 16 February 2010 17:26

ok, je te remercie pour ta réponse!

L'opération a pris son temps mais au final, j'ai pu parvenir à une vectorisation via la commande
r.to.vect

Dans le raster d'origine, les contours étaient saccadés. Y a-t-il moyen de lisser les lignes noires d'une image via GRASS?

MathieuR · Tue 16 February 2010 18:37

Avec la commande, j'ai toujours une erreur pour l'écriture du fichier temporaire:

Code:

r.thin in=raster out=rasterthin --v                                         
File farinole -- 131957 rows X 108701 columns
ERROR: (null): Error writing temporary file
(Tue Feb 16 18:33:03 2010) Command finished (0 sec)

Sinon, où se trouve le fichier .grassrc? Je ne l'ai pas trouvé..

Merci!

Dernière modification par baobazz (Tue 16 February 2010 18:39)

vblanque · Tue 16 February 2010 21:59

Vous pouvez peut être essayer avec
gdal_polygonize qui permet de
vectoriser en polygones des rasters.

MathieuR · Wed 17 February 2010 11:48

Oui, mais j'utilise FWTools - version 2.4.7 et l'outil gdal_polygonize ne fonctionne pas.
Il indique que la version de gdal est trop ancienne pour l'outil
le même souci qu'ici
Chez vous, ça marche sous Windows et si oui, l'utilisez-vous via FWtools ou via l'installeur oSGeo4W?

Sinon, où puis-je trouver le fichier .grassrc6 afin de la commande r.thin plus verbeuse?

Dernière modification par baobazz (Wed 17 February 2010 12:02)

Sylvain_M · Wed 17 February 2010 13:00

le fichier .grassrc6 doit se trouver dans le dossier "Mes Documents"

vu le genre d'erreur affichée, je pense qu'effectivement ca peut etre bien que tu passe grass en mode debug car ce genre d'erreur me semble plus lié à un souci de place ou de memoire qu'a la fonction r.thin en elle-meme ...

MathieuR · Wed 17 February 2010 13:02

après plusieurs heures, la commande r.thin m'affiche le résultat suivant:

Code:

r.thin in=farinole out=farinoleThin                                             
File farinole -- 140000 rows X 90000 columns
ERROR: (null): Unable to find bounding box for lines
(Wed Feb 17 12:53:51 2010) Command finished (9949 sec)

D'où cela peut-il venir?
l'image à l'origine n'était pas géoréférencée.
J'avais écrasé les définitions de projection par les définitions courantes de la région lors de l'import de l'image
r.in.gdal avec l'option -o
Peut-être cela vient-il de là?

arnaud17 · Mon 18 June 2012 15:53

Bonjour Baobazz,

Je suis face au même problème que toi, avec ce message d'erreur "ERROR: (null) : Unable to find bounding box for lines"
Le plus surprenant dans mon cas, c'est que je lance le module sur plusieurs raster à la fois et je n'ai ce problème que sur un raster. J'ai essayé le mode bavard comme indiqué dans la discussion, sans résultat également.
Comment as tu fait pour régler ce problème ?

Merci d'avance.

MathieuR · Mon 18 June 2012 16:39

en fait, en réglant correctement la région: résolution + étendue.

gene · Mon 18 June 2012 18:33

Comme toujours, dans GRASS, la notion de région est fondamentale. On a beau le répéter, l'écrire, le graver sur la roche, ça ne rentre pas...

MathieuR · Tue 19 June 2012 08:30

en fait, lors de l'opération de thinning, que j'avais tentée de réaliser sur la bdparcellaire image soit à une résolution de 0.5 m, la mémoire demandée était top élevée. Dans ce cas, soit on diminue la résolution: 1m, etc... ou bien on ajuste l'étendue. Or, il est préférable d'ajuster la région dans le cas d'une tentative de vectorisation de la bdparcellaire. Il faut donc sélectionner une étendue plus petite. En gros, cela implique de diviser l'étendue globale du secteur à vectoriser en entités/dalles plus petites. Assez fastidieux quand le nombre de dalles est élevé. D'où nécessité de scripter tout cela. Et oui, la région est fondamentale.

Sylvain_M · Tue 19 June 2012 11:02

Bonjour,

en dehors de la problématique de l'utilisation de GRASS (que j'adore), est-ce qu'il ne serait pas plus rapide pour vous d'utiliser directement le Registre Parcellaire Graphique en format vectoriel ?

MathieuR · Tue 19 June 2012 13:13

Le RPG ici dont vous parlez est agricole.

En ce qui concerne la bdparcellaire vectorisée, elle n'est pas pas dispo sur toutes les communes

Dernière modification par baobazz (Tue 19 June 2012 13:13)

gene · Tue 19 June 2012 22:14

Si vous utilisez la version 6.4.2 tout est scriptable en Python, ce qui rend les choses plus faciles pour les utilisateurs de Windows. Pour un exemple de boucles, regardez le script que j'ai fourni sur les listes de GRASS

http://osgeo-org.1560.n6.nabble.com/Aut … 78801.html

MathieuR · Wed 20 June 2012 13:19

j'avais scripté aussi pour la vectorisation du cadastre mais bon,...pas mal de retouches à effectuer au final. je n'ai pas trouvé de méthode optimale.

Sylvain_M · Thu 21 June 2012 09:52

ah oui, autant pour moi, je n'ai pas fait attention à la nuance ...

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Tue 16 February 2010 15:50

[GRASS][vectorisation] r.thin bloque

Code:

Code:

Code:

#2 Tue 16 February 2010 16:35

Re: [GRASS][vectorisation] r.thin bloque

#3 Tue 16 February 2010 17:26

Re: [GRASS][vectorisation] r.thin bloque

#4 Tue 16 February 2010 18:37

Re: [GRASS][vectorisation] r.thin bloque

Code:

#5 Tue 16 February 2010 21:59

Re: [GRASS][vectorisation] r.thin bloque

#6 Wed 17 February 2010 11:48

Re: [GRASS][vectorisation] r.thin bloque

#7 Wed 17 February 2010 13:00

Re: [GRASS][vectorisation] r.thin bloque

#8 Wed 17 February 2010 13:02

Re: [GRASS][vectorisation] r.thin bloque

Code:

#9 Mon 18 June 2012 15:53

Re: [GRASS][vectorisation] r.thin bloque

#10 Mon 18 June 2012 16:39

Re: [GRASS][vectorisation] r.thin bloque

#11 Mon 18 June 2012 18:33

Re: [GRASS][vectorisation] r.thin bloque

#12 Tue 19 June 2012 08:30

Re: [GRASS][vectorisation] r.thin bloque

#13 Tue 19 June 2012 11:02

Re: [GRASS][vectorisation] r.thin bloque

#14 Tue 19 June 2012 13:13

Re: [GRASS][vectorisation] r.thin bloque

#15 Tue 19 June 2012 22:14

Re: [GRASS][vectorisation] r.thin bloque

#16 Wed 20 June 2012 13:19

Re: [GRASS][vectorisation] r.thin bloque

#17 Thu 21 June 2012 09:52

Re: [GRASS][vectorisation] r.thin bloque

Pied de page des forums