<div dir="ltr"><div><div>Thanks for the suggestion, but if I think that if I follow the procedure which you have indicated I will only be able to partially solve my problem. What about those parts of the polygons that are not boundaries? These might well overlap with certain cells of the raster grid...<br>
<br></div>thanks,<br><br></div>Paolo<br><div><div><br><br><div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/7/9 Moritz Lennert <span dir="ltr"><<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="HOEnZb"><div class="h5">On 09/07/13 13:21, Paolo Ruzza wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear GRASS users,<br>
<br>
I am confronted with a "non-standard" problem and I would appreciate any<br>
help that you can give me. I have several vector files (polygons) and an<br>
arbitrary raster grid, only with category numbers associated to it. I<br>
would like to overlay the vectors onto the raster grid in order to query<br>
it and extract the categories of the raster grid which overlap with the<br>
polygons of each vector file. I am not interested in getting any info<br>
from the vector polygons, I just want know where the vector and the<br>
raster overlap. I have already thought of a way of doing this-I guess it<br>
would naturally come to mind of every GRASS user. This would involve<br>
converting the vectors into a rasters and then query the raster grid and<br>
the rasterized vectors with r.stats. This would be very easy to do, but<br>
my vectors have very small islands, which can, of course, be converted<br>
into raster cells, if I set them my region to a very fine resolution<br>
for. However, if I went down this route, this would: i) massively slow<br>
down vector to raster conversion ii) massively slow down the r.stats<br>
querying process. Just to give you an idea of the magnitude of my<br>
problem I am using several vector layers (hundreds) at global level and<br>
a raster grid at a 2 arc-minutes resolution. I am kind of stuck and I am<br>
not sure if there is any computationally efficient strategy to do what I<br>
want. Any help would be greatly appreciated!<br>
</blockquote>
<br></div></div>
One possible option I see allowing you to keep your original resolution is to use v.to.points (type=boundary and dmax sufficiently small) to transform your polygon boundaries to points, then v.to.rast to transform that to raster.<span class="HOEnZb"><font color="#888888"><br>

<br>
Moritz<br>
</font></span></blockquote></div><br></div></div></div></div></div>