[GRASS-user] Site in Polygon Query

[Martin Wegmann]

Hello Michael,

On Wednesday, 23. May 2007 12:29:41 Michael Malicky wrote:
> Dear All!
>
> I am using GRASS for some time now, mainly for plotting distribution maps
> of plants and animals. Now a user came up with a different request:
> We have a map consisting of polygons and another one with sites. We need a
> query to filter all polygons that include at least one site. Is there a
> special command for this or a script ? Any help is welcome.

just an idea:

v.to.rast of polygons and sites (use=val)

r.statistics method=sum 
and then 
r.mapcalc "output = if(map_sum > 1,1,null())"

regards, Martin



> sincerely yours
>
> Michael
>
> DI Michael Malicky
> EDV Leiter
> Oberösterreichische Landesmuseen/Biologiezentrum
> Johann-Wilhelm-Klein Str. 73
> A-4040 Linz
> Österreich
> Tel.: +43/732/759733-33
> FAX.: +43/732/759733-99
> Mobil: +43/664/60072/52333
> E-Mail: m.malicky at landesmuseum.at
>
>
> _______________________________________________
> grassuser mailing list
> grassuser at grass.itc.it
> http://grass.itc.it/mailman/listinfo/grassuser
-- 
Martin Wegmann

DLR - German Aerospace Center
German Remote Sensing Data Center
@
Dept.of Geography
Remote Sensing and Biodiversity Unit
University of Wuerzburg
Am Hubland
97074 Würzburg

phone: +49-(0)931 - 888 4797
mobile: +49-(0)175 2091725
fax:   +49-(0)931 - 888 4961