Le portail francophone de la géomatique

Si personne ne sait pour les .DEM peut etre quelqu'un poura m'aider pour des MNT en .asc
J'ai un mnt en .asc que je voudrai ouvrir sur QGIS via le plugin grass afin de pouvoir réaliser des analyse dessus. Je n'arrive pas à l'importer, j'ai le même problème qu'avec le .dem . J'ai utilisé plusieur modul d'import dans les outils grass mais à chaque fois que j'essai d'ouvrir le raster créé, qgis se plante et se ferme sans message d'erreur

j'ai essayé de l'ouvrir directement dans grass, le fichier s'ouvre mais l'image ne s'affiche pas, enfin il y a une image blanche avec des valeurs de pixels un peu bizards.

est ce que quelqu'un voit le problème?

Merci d'avance

non rectification j'ai le même problème avec la couche raster "gtopo30" de la demo Qgis plante quand j'essai de l'ouvrir

j'ai un problème carement d'ouverture des couches raster grass ; de quoi ça peut venir?

Bonjour,

Merci beaucoup pour cette information !!

Y.

Bonjour,

Quels sont les messages d'erreurs éventuels ? La commande utilisée ?

Merci,

Y.

Bonjour,

La réponse à ta question est dans l'onglet "Manuel" du module GRASS.


DESCRIPTION r.in.ascii allows a user to create a (binary) GRASS raster map layer from an ASCII raster input file with (optional) TITLE.
The GRASS ASCII input file has a header section which describes the location and size of the data, followed by the data itself.
The header has 6 lines:
north:   xxxxxx.xx
south:   xxxxxx.xx
east:    xxxxxx.xx
west:    xxxxxx.xx
rows:    r
cols:    c
The north, south, east, and west field values entered are the coordinates of the edges of the geographic region. The rows and cols field values entered describe the dimensions of the matrix of data to follow. The data which follows is r rows of c integers.
Optionally the following parameters can be defined in the header section:
null: nn
type: float
multiplier: 2.
"null" defines a string or number to be converted to NULL value (no data).
"type" defines the data type (int, float double) and is not required.
"multiplier" is an optional parameter to multiply each cell value.


Si tu ouvres un fichier ASCII de la BD Alti tu verras que l'entête est écrit de la façon suivante :

ncols         5201
nrows         3201
xllcorner     339987.50
yllcorner     2419987.50
cellsize      25
NODATA_value  -9999


Par conséquent, tu dois réécrire l'entête avec les valeurs correspondantes en utilisant TextEdit ou l'équivalent sous Windows.

Exemple :

La BD Alti 250 m en L93 est écrite comme cela :

ncols         4641
nrows         4321
xllcorner     89875
yllcorner     6039875
cellsize      250
NODATA_value  -9999

Réécrit, cela donne :

north:        7120125
south:        6039875
east:        1250125
west:        89875
rows:        4321
cols:        4641

Et il ne te reste plus qu'à recommencer l'importation des données dans GRASS, cela fonctionnera sans aucun problème.

Rappel : dans le module GRASS de QGis, n'hésitez pas à consulter l'onglet "Manuel" si vous avez des problèmes avec GRASS. C'est peut-être parfois abscons, mais il y a toutes les réponses.

Bonne journée.

Dominique.

Dernière modification par Dominique Toquet (Wed 30 June 2010 08:45)

Réécrit seulement l'entête du fichier .asc, pas le fichier .txt joint avec.

Cela doit se présenter comme cela :

north:        7120125
south:        6039875
east:        1250125
west:        89875
rows:        4321
cols:        4641
0 0 0 0 0 0 0 0 0 0 0 0 0 0 et plein d'autres zéros.

Ensuite, dans le module r.in.asc tu précises bien à la ligne "String representing NULL value data cell" la valeur -9999.

Essaye avec le fichier modifié ci-joint (c'est le fichier MNT250_L93_FRANCE.ASC téléchargeable sur le site de l'IGN), le lien expirera dans 24 heures.

http://files.me.com/dominique.toquet/qlcjxb

Pour obtenir les valeurs North, South, East and West de tes données, tu vas dans les métadonnées.

"The geographic coordinates north, south, east, and west describe the outer edges of the geographic region. They run along the edges of the cells at the edge of the geographic region and not through the center of the cells at the edges. The NW value occurs at the beginning of the first line of data, and the SW value occurs at the beginning of the last line of data. "

Si cela ne fonctionne pas, le plus simple, sauf si tu veux aller au bout du problème, c'est d'utiliser la commande r.in.arc avec le fichier d'origine.

Je viens de faire la manipulation trois fois de suite sans aucun problème.

Et si, en fin de compte, tu n'obtiens rien, vérifie ta région et ton mapset.

Dominique.

Re,

C'est cool. ça avance de mon côté.

J'ai réussi à charger votre fichier, puis à réécrire l'entête du fichier. Les couches se chargent sans aucun problème (un gros carré gris est dessiné).

Cependant, je n'arrive pas à transposer ce gris en une palette de couleur qui me différencierait les altitudes. Sous l'onglet convention des signes/symbologie, le bouton colorisé par 3 bandes n'est pas en surbrillance.

Concernant les valeurs North, South, East and West de mes données : je vais dans n'importe quelle couche vectorielle "fonctionnelle" (par exemple une couche sous-bassin versant), vais sous l'onglet métadonnées, et prend les Xmin, Ymin, Xmax, Ymax pour remplir le nord sud est et ouest?

Cyrille

Il m'affiche le message ci-dessous quand j'essaye de charger le fichier.

C:/Documents and Settings/user/Bureau/Dossier sommaire de candidature/SIG/BDAlti/cotedor/France_250_ASC_L93/MNT250_L93_FRANCE.ASC n'est pas une source de données raster gérées

Symbologie

Le MNT est un raster 1 bande.

C'est ok, cela fonctionne bien pour la visualisation.

-----------------

Transformation du raster en vecteur

Afin de pourvoir travailler la couche, je l'ai converti en vecteur (grâce à la boîte à outils GRASS). Sans m'en rendre compte, il me l'a convertie sur une zone un peu plus large que mon bassin versant. C'est du à la création de mon Mapset?

-----------------

Après avoir créer cette couche vecteur, je n'ai plus accès aux outils GRASS. A quoi cela est-il dû?

-----------------

J'aurais aimé, à partir de cette couche vecteur :
- dessiner des sous-bassin versant,
- dessiner une vue 3D de mon bassin versant.

Quelles sont les manipulations à faire?

-----------------

Quelle est l'extension la plus pertinente pour travailler des rasters? est ce que cela vaut la peine de se "prendre le chou" à comprendre GRASS?

-----------------

Au fait, un premier merci pour ce soutien technique sous GRASS et Qgis!

Cyrille

Bonjour, tu peux dessiner un vecteur polygonal correspondant à ta zone d'étude, tu la transforme en raster grâce au module v.to.rast puis créer un masque grâce à la fonction r.mask en entrée ta couche que tu viens de rasteriser, et normalement tu ne travailleras que sur cette zone.

Concernant la création de sous bassin versants, des vues 3D,....tout est dans r.watershed.

Maintenant yapluka!

Re,

Je continue ma découverte de Qgis.

Je suis déjà arrivé à créer mon MNT, et c'est déjà pas mal! je remercie d'ailleurs ceux qui m'ont accompagner!

Diophaminou, c'est bien le fichier g.lock qui bloquait. A quoi correspond il?

J'aurais maintenant voulu utiliser ma Bd Alti pour faire des modèles hydrologiques.

J'ai donc testé r.watershed avec ma bdalti, mais comme d'hab, j'ai quelques soucis. Le message d'erreur est en pièce jointe.

Que faut il respecter pour utiliser le module hydrologic modelling?

Il n'est pas possible d'utiliser des données vecteurs pour faire des analyses spatiales?

Merci. Cyrille.

Dernière modification par Cyrille BOUVIER (Wed 07 July 2010 11:31)