banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

Annonce

Rencontres QGIS 2025

L'appel à participation est ouvert jusqu'au 19 janvier 2025!

#1 Wed 13 April 2011 09:41

Sylvain PIERRE
Participant assidu
Lieu: Strasbourg
Date d'inscription: 6 Sep 2005
Messages: 170

Grass- perte de données à l'import (r.in.xyz)

Bonjour,

Voici mon problème.
J'importe une série de fichier xyz dans grass (dalles de 1000*1000, projection Lambert 1 nord) avec le script suivant :

Code:

for i in `ls -1 *.xyz`
do
                lay=${i/.xyz/}
                eval `cat $i | r.in.xyz -s -g input=- output=test fs=,`
                g.region n=$n s=$s w=$w e=$e res=1 -p
                g.region n=n+0.5 s=s-0.5 w=w-0.5 e=e+0.5 -p
                cat "$i" | r.in.xyz input=- output="$lay" fs=,
done

Tout marche bien sauf que j'ai perdu systématiquement les 128 premiers points de chaque fichier originel. Cela n'est visible que lorsqu'on zoom fortement sur la zone concernée, les 128 premiers pixels de la première ligne du raster au nord-ouest sont vides.
Aucun message d'erreur ...

Bizarre !?

Cordialement

Hors ligne

 

#2 Wed 13 April 2011 10:00

Yves
Membre du bureau
Lieu: Aix-les-Bains
Date d'inscription: 22 Mar 2006
Messages: 9869
Site web

Re: Grass- perte de données à l'import (r.in.xyz)

Sylvain,

Qu'est tu censé faire dans cette ligne ?

Code:

g.region n=n+0.5 s=s-0.5 w=w-0.5 e=e+0.5 -p

ne serait pas plutôt n=$n+0.5 ?

Y.


Yves Jacolin, bénévole de l'association GeoRezo.net, agit au nom et pour le compte de l'association - Partageons ce qui nous départage !!  - GeoRezo vous aide ? Aidez GeoRezo !

Hors ligne

 

#3 Thu 14 April 2011 09:26

Sylvain PIERRE
Participant assidu
Lieu: Strasbourg
Date d'inscription: 6 Sep 2005
Messages: 170

Re: Grass- perte de données à l'import (r.in.xyz)

D'abord merci de t'intéresser au sujet ;-)

Concernant ta remarque: non la syntaxe est bonne. Ca me redimensionne bien ma region en 1000 X 1000 (1000000 cells) alors qu'après le premier passage de la commande g.region la résolution est 999 X 999
Par contre c'est à l'import que les 128 premiers points disparaissent:

999872 points found in region


je penche pour un bug : 128 c'est un peu un chiffre "magique" en informatique...

Hors ligne

 

#4 Thu 14 April 2011 10:46

Yves
Membre du bureau
Lieu: Aix-les-Bains
Date d'inscription: 22 Mar 2006
Messages: 9869
Site web

Re: Grass- perte de données à l'import (r.in.xyz)

Bonjour,

il y a combien de valeur nodata ?

Y.


Yves Jacolin, bénévole de l'association GeoRezo.net, agit au nom et pour le compte de l'association - Partageons ce qui nous départage !!  - GeoRezo vous aide ? Aidez GeoRezo !

Hors ligne

 

#5 Fri 15 April 2011 20:07

Sylvain PIERRE
Participant assidu
Lieu: Strasbourg
Date d'inscription: 6 Sep 2005
Messages: 170

Re: Grass- perte de données à l'import (r.in.xyz)

128

Yves a écrit:

Bonjour,

il y a combien de valeur nodata ?

Y.

Hors ligne

 

#6 Fri 15 April 2011 20:08

Yves
Membre du bureau
Lieu: Aix-les-Bains
Date d'inscription: 22 Mar 2006
Messages: 9869
Site web

Re: Grass- perte de données à l'import (r.in.xyz)

Bonjour,

Ben voilà, tes valeurs nodata n'ont pas été importées, et je pense que c'est logique !

Y.


Yves Jacolin, bénévole de l'association GeoRezo.net, agit au nom et pour le compte de l'association - Partageons ce qui nous départage !!  - GeoRezo vous aide ? Aidez GeoRezo !

Hors ligne

 

#7 Tue 19 April 2011 09:23

Sylvain PIERRE
Participant assidu
Lieu: Strasbourg
Date d'inscription: 6 Sep 2005
Messages: 170

Re: Grass- perte de données à l'import (r.in.xyz)

Bonjour,

Je me suis trompé:
J'ai 0 valeurs nodata. Je perds bien les 128 premières lignes de mes fichiers...

Hors ligne

 

#8 Tue 19 April 2011 09:24

Sylvain PIERRE
Participant assidu
Lieu: Strasbourg
Date d'inscription: 6 Sep 2005
Messages: 170

Re: Grass- perte de données à l'import (r.in.xyz)

Trouvé!
    cat "$i" | r.in.xyz input=- output="$lay" fs=,  "oubli les 128 premieres lignes du fichier d'entrée
    r.in.xyz input="$i" output="$lay" fs=, marche correctement
Problème de compatibilité de commande unix sous windows?

Hors ligne

 

Pied de page des forums

Powered by FluxBB