<div dir="ltr"><div>Hi Markus,<br><br></div>Good idea. I've downloaded the NC dataset and tried it using the following commands (resizing the region to save time):<br><br>g.region res=30<br>g.region n=228500 s=228000 e=630500 w=630000<br>
r.to.vect --overwrite input=elev_srtm_30m@PERMANENT output=elev_srtm_30m_vec feature=point<br>v.db.addcol map=elev_srtm_30m_vec@user1 columns=srtm double precision<br>v.what.rast vector=elev_srtm_30m_vec@user1 raster=elev_srtm_30m@PERMANENT column=srtm<br>
<div><div><br></div><div>Checking just the NW corner, I'm getting the right results, so the issue doesn't appear to be with v.what.rast. Are there any issues with importing ascii grids or shape files that could cause this?<br>
<br></div><div>My example dataset is available <a href="http://ubuntuone.com/24Gb4JLTlXO2JHPPTj5GCQ">here</a>, made up of a shp file of points and an ascii raster and also these imported into a grass location and mapset. The commands I used are as follows:<br>
<br>r.in.arc input=/home/filepath/MeanTemp_Apr_1961-1990_LTA.txt output=Temp_apr<br>v.in.ogr -o dsn=/home/filepath/MO_grid_pts_GB.shp output=MO_grid_pts_GB<br>v.db.addcol map=MO_grid_pts_GB@PE columns=apr double precision<br>
v.what.rast vector=MO_grid_pts_GB@PE raster=Temp_apr@PE column=apr<br><br></div><div>When I then use the query raster/vector maps tool in the map display at, for example, grid ref 142500,962500 (top of the Western Isles) and the two points to the east I find the vector data attribute is taken from the raster cell 5 km to the east. As a quick check, this is the same in Edinburgh too.<br>
<br>Help greatly appreciated!<br></div><div>Michael<br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 10 July 2013 00:56, Markus Neteler <span dir="ltr"><<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Michael,<br>
<div><div class="h5"><br>
On Mon, Jul 8, 2013 at 6:43 PM, Michael Spencer<br>
<<a href="mailto:spencer.mike.r@gmail.com">spencer.mike.r@gmail.com</a>> wrote:<br>
> Hi,<br>
><br>
> I've a problem with the v.what.rast tool. I imported a vector grid of points<br>
> with a spacing of 5 km (v.in.ogr) to GRASS, along with 12 rasters<br>
> (r.in.arc). There is a point in the centre of each raster cell, with some<br>
> extra points as a buffer around the edge (still on a 5 km grid).<br>
><br>
> Having added extra columns to the point vector file I ran v.what.rast to<br>
> extract the value of the 12 rasters, in turn, at each point location.<br>
><br>
> Checking the result in the GRASS viewer the tool has appended data from the<br>
> raster cell 5 km to the east at all locations.<br>
><br>
> Is this a bug? Has anyone else experienced it? And most importantly - does<br>
> anyone know how to fix it?!<br>
<br>
</div></div>We would need an example with e.g. the North Carolina sample dataset to<br>
easily reproduce (or not) the issue:<br>
<br>
<a href="http://grass.osgeo.org/download/sample-data/" target="_blank">http://grass.osgeo.org/download/sample-data/</a><br>
<span class="HOEnZb"><font color="#888888"><br>
Markus<br>
</font></span></blockquote></div><br></div>