<br><br><div class="gmail_quote">2010/7/1 Daniel Victoria <span dir="ltr"><<a href="mailto:daniel.victoria@gmail.com">daniel.victoria@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Or how about converting the polygons to raster and then run v.what.rast to add a collumn with the polycon categorie value? Or even better, v.what.vect?<br>
<br>
Or did I misunderstand the question?<br></blockquote><div><br>Hi Daniel, hi Markus,<br><br>v.what.vector was the solution. Thank you.<br><br>I had been looking in completely the wrong place. <br><br>Previously I had successfully overlayed two vector maps and got all the intersecting polygons with all the attributes of both maps. Also somewhere I read that points are no longer treated as sites but also as vectors.<br>
<br>I concluded that points (imported with v.in.ascii) with attributes could also be overlayed with a vector map, resulting in a new map with the points having the attributes of the original points and all the attributes of the polygon in which they fall. That does not work (why not?). v.overlay simply returns the original polygons with existing but empty columns for the attributes of the points. It does not report any kind of error in the process.<br>
<br>So I tried to find a way to get v.overlay to work with points and never considered other solutions.<br><br>Thanks again,<br>stn<br> </div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<font color="#888888"><br>
Daniel<br>
</font><div><div></div><div class="h5"><br>
On Wed, Jun 30, 2010 at 6:31 PM, Markus Neteler <<a href="mailto:neteler@osgeo.org">neteler@osgeo.org</a>> wrote:<br>
> On Tue, Jun 29, 2010 at 9:25 AM, stn <<a href="mailto:stneumann@web.de">stneumann@web.de</a>> wrote:<br>
>> Hi,<br>
>><br>
>> how can I identify which point lies in which polygon?<br>
>><br>
>> The points are in an ascii-table, each with latitude, longitude and some<br>
>> other attributes like name, altitude and id-number.<br>
>> The vector-map is a map with political boundaries found in the vmap0-dataset<br>
>> at <a href="http://www.gdf-hannover.de/gis.php?id=5" target="_blank">http://www.gdf-hannover.de/gis.php?id=5</a><br>
><br>
> I don't know an elegant solution to this but you may<br>
> create the list of polygon categories with v.categories and print method,<br>
> then loop over the polygons to run v.extract + v.select on each<br>
> extracted polygon.<br>
><br>
> Likely there are better solutions!<br>
><br>
> Markus<br>
</div></div><div><div></div><div class="h5">> _______________________________________________<br>
> grass-user mailing list<br>
> <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
> <a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
><br>
</div></div></blockquote></div><br>