<html><head><style type="text/css">body {word-wrap: break-word; background-color:#ffffff;}</style></head><body><div style="font-family: sans-serif; font-size: 16px">Once you've added columns for the pixell values to the point vector, then you can, of course, display those attributes on your map.<br><br>--<br>Micha Silver<br><a href="tel:0523665918">052-366-5918</a></div><br><br>-----Original message-----<br><blockquote style="; border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><div style="font-family: sans-serif; font-size: 14px"><b>From: </b>Helena Herrera <helenaherrera1980@gmail.com><b><br>To: </b>Micha Silver <micha@arava.co.il><b><br>Cc: </b>grass-user@lists.osgeo.org<b><br>Sent: </b>Thu, 29 Mar 2012, 14:36:12 GMT+02:00<b><br>Subject: </b>Re: [GRASS-user] How to retrieve raster pixel value over a vector point layer<br><br></div><div>Hi<br><br>I was thinking in something much simpler like just printing Pixel values next to Point ID.<br><br>I will give a try on your suggestions.<br>Thanks<br><br><div class="gmail_quote">2012/3/28 Micha Silver <span dir="ltr"><<a href="mailto:micha@arava.co.il">micha@arava.co.il</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="direction:ltr" bgcolor="#FFFFFF" text="#000000"><div class="im">
On 03/28/2012 08:11 PM, Helena Herrera wrote:
<blockquote type="cite">
<div>Greetings</div>
<div> </div>
<div>I have a raster layer and a vector point layer and I need to
associate pixel values with points:</div>
</blockquote>
<br></div>
When you say "associate" If you mean to add the raster values to a
column in the point vector attribute table, then v.what.rast is what
you're looking for.<div class="im"><br>
<br>
<blockquote type="cite">
<div>1- pixel values that are over points;</div>
<div>2- 3x3 pixel neighbor pixel values over a point.</div>
</blockquote>
<br></div>
Create a new raster with r.neighbors then run v.what.rast again to
upload these 3x3 neighborhood values to another column in the point
vector. <br><div class="im">
<br>
<blockquote type="cite">
<div>Is this possible do to in a more automatic way than query
each by each point?</div>
</blockquote>
<br></div>
Here are the steps:<br>
v.db.addcol vector_map col="rast_val double, rast_val_3x3 double"<br>
v.what.rast vector_map rast=raster_map col=rast_val<br>
r.neighbors raster_map out=raster_map_3x3 # assumes averaging a 3x3
window <br>
v.what.rast vector_map rast=raster_map_3x3 col=rast_val_3x3<br>
<br>
<blockquote type="cite">
<div>Thanks</div>
<div>Regards</div>
<div>helena</div>
<br>
This mail was received via Mail-SeCure System.<br>
<br>
<fieldset></fieldset>
<br>
<pre>_______________________________________________
grass-user mailing list
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a>
This mail was received via Mail-SeCure System.
</pre><span class="HOEnZb"><font color="#888888">
</font></span></blockquote><span class="HOEnZb"><font color="#888888">
<p><br>
</p>
<br>
<pre cols="72">--
Micha Silver
GIS Consultant, Arava Development Co.
<a href="http://www.surfaces.co.il" target="_blank">http://www.surfaces.co.il</a></pre>
</font></span></div>
</blockquote></div><br>
<br>
This mail was received via Mail-SeCure System.<br>
</div></blockquote></body></html>