[Grass-de] r.in.xyz; Anzahl von Vektorpunkten pro Rasterzelle
Christoph Heibl
christoph.heibl at gmx.net
Mo Aug 17 13:19:35 CEST 2009
Hallo Markus,
Mein Fehler: Ich hatte ' anstatt ` gesetzt. Meine shell verwendet
auch bash und mit der richtigen Syntax funktioniert es jetzt problemlos.
Vielen Dank für die schnelle und kompetente Hilfe!
Grüße,
Christoph
On Aug 17, 2009, at 1:00 PM, Markus Neteler wrote:
> Hallo Christoph,
>
> 2009/8/17 Christoph Heibl <christoph.heibl at gmx.net>:
>> Hallo Markus,
>>
>> Vielen Dank! Das hat wunderbar geklappt.
>>
>> Nur noch eine kleine Frage hätte ich noch dazu:
>>
>> eval `g.region -g` hat bei mir nicht funktioniert,
>
> Was ist denn passiert? Der Befehlt arbeitet "schweigend" da
> er ja Variablen setzt.
>
>> d.h. ich habe die
>> Spaltennummer im r.mapcalc-Befehl von Hand eingeben müssen. Woran
>> könnte das
>> liegen? (Ich verwende GRASS 6.4 in der 'Terminal' shell unter MAC OS
>> 10.4.11)
>
> Was sagt
> echo $SHELL
> ?
>
> Ich benutze /bin/bash, falls Du eine andere Shell hast, wird es eval
> vielleicht nicht geben (oder unter anderem Namen).
>
>
>>
>> Der Vollständigkeit halber hier noch kurz wie ich nach aus den
>> aktualisierten Vektordaten die Rasterkarte erzeugt habe:
>>
>> --------
>> echo "select x,y,count FROM random20 NATURAL INNER JOIN (SELECT
>> nummer,count(nummer) FROM random20 GROUP BY nummer ORDER BY
>> nummer) AS
>> random20_tmp" | db.select | v.in.ascii out=random20_s skip=1 x=1 y=2
>> columns='x double precision, y double precision, spec varchar(30),
>> sect
>> varchar(30), count int' --overwrite
>>
>> v.to.rast random20_s out=sampling use=attr column=count --o
>> ---------
>>
>> Eventuell geht das auch leichter ...
>
> Sieht gut aus - nur 2 Zeilen :)
>
> Viele Grüße,
> Markus