#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
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
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
Bonjour,
il y a combien de valeur nodata ?
Y.
Hors ligne
#6 Fri 15 April 2011 20:08
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