[Grass-de] r.in.xyz; Anzahl von Vektorpunkten pro Rasterzelle

[Markus Neteler]

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