[GRASS-de] Verschneidung von Rasterdaten

[Otto Dassau]

On Wed, 3 Aug 2005 09:25:01 +0200
Maik Trömel <maik.troemel at maitro.net> wrote:

> Hallo,
>  
> ich habe folgendes Problem.
> Ich habe Postleitzahlregionen als Vektordaten und Radarbilder mit
> Niederschlagsinformationen im Rasterformat vorliegen. Nun möchte ich
> anhand der Radarbilder herausfinden ob es in einer bestimmten PLZ-Region
> gehagelt hat. Ich habe nun ersteinmal die Vektordaten in Rasterdaten
> konvertiert und als Raster-Values die PLZ mitgegeben.
>  
> Mit dem Befehl r.mapcalc müsste man doch nun irgendwie herausbekommen
> welche Werte das Radarbild in den einzelnen PLZ-Regionen annimmt. Wie
> kann man so etwas realisieren?
> Oder wäre es sinnvoller mit einer Maske zu arbeiten?

Bestimmte Werte (Hagel) innerhalb eines bestimmten Gebietes könntes t du so herausfischen: 

r.mapcalc hagel30177="if(PLZkarte==30177 && \ regenkarte==hagel,regenkarte,null())"

Vielleicht hilft das?

  Gruesse
   Otto