[GRASS-de] vektormaske ?

Florian Jurgeit Florian.Jurgeit at uibk.ac.at
Mi Jun 16 10:23:29 CEST 2004


Jochen Schwarze wrote:

> Hallo!
> Mein Vorschlag:
> Die gewünschten Polygone anhand ihrer Katergorie selektieren und mit 
> v.to.rast in Rasterdaten konvertieren. Diese lassen sich ja dann mit 
> r.mask weiterverarbeiten. Das sollte im Batchbetrieb ganz gut 
> funktionieren.
>
> Gruss, Jochen
>
>  
>
Genau so hab ich das ganze dann erledigt - aber gleich mit r.in.shape 
... - dann die ID's der Flächen aus einer txt-Datei mit cat in einer 
for-Schleife geholt:
Hier der einfache Code:
# This Script writes all altitude-values of a defined region (e.g. Polys 
from a shape-file) in a text-file (a Perl-script makes a 
hypsomertic-curve ...)
# Florian Jurgeit, 2004
#************************************************************

# Don't forget to import the Shape-File (Poly-Features) to a GRASS-GRID 
with:
#   r.in.shape in=Export_gem_nt_st.shp out=gem_ntst_grid cat=ID 
lab=GEM_NAME
# It's important that something like an ID is defined as cell-values ...
#
# Don't forget to set the region ...

#Variables:
basedem=srtm;
featuregrid=gem_ntst_grid;
#id=170220;


for id in `cat gemids.txt`;
  do echo processing $id
  echo "Creating tmp-Grids ..."
  r.mapcalc "mask_tmp = if(($featuregrid != $id),null(),1)"

  r.mapcalc "dem_tmp = (mask_tmp*$basedem)"

  echo "Building stats..."
  r.stats -1 -n input=dem_tmp | sort -n -r > hypso/$id.txt

  echo "Removing tmp-grids..."
  g.remove rast=dem_tmp
  g.remove rast=mask_tmp
done


_______________________________________________
grass-de mailing list
grass-de at grass-verein.de
https://grass-verein.de/mailman/listinfo/grass-de
-------------- nächster Teil --------------
Ein Dateianhang mit Binärdaten wurde abgetrennt...
Dateiname   : Florian.Jurgeit.vcf
Dateityp    : text/x-vcard
Dateigröße  : 275 bytes
Beschreibung: nicht verfügbar
URL         : <https://lists.fossgis.de/pipermail/fossgis-talk-liste/attachments/20040616/2be8e8d5/Florian.Jurgeit.vcf>