<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Oct 10, 2014 at 9:00 AM, Paulo van Breugel <span dir="ltr"><<a href="mailto:p.vanbreugel@gmail.com" target="_blank">p.vanbreugel@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Thu, Oct 9, 2014 at 11:49 PM, 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:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">On Wed, Oct 8, 2014 at 9:44 AM, Paulo van Breugel<br>
<span><<a href="mailto:p.vanbreugel@gmail.com" target="_blank">p.vanbreugel@gmail.com</a>> wrote:<br>
> On Wed, Oct 8, 2014 at 9:21 AM, Markus Neteler <<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>> wrote:<br>
>><br>
>> On Tue, Oct 7, 2014 at 11:20 PM, Paulo van Breugel<br>
>> <<a href="mailto:p.vanbreugel@gmail.com" target="_blank">p.vanbreugel@gmail.com</a>> wrote:<br>
>> > I am trying to upload raster values of the raster 'grid' to the<br>
>> > attribute<br>
>> > table of the gbif vector layer, as below:<br>
>> ><br>
>> > g.region rast=grid<br>
>><br>
>> can you please post<br>
>><br>
>> g.region rast=grid -p<br>
><br>
><br>
> GRASS 7.1.svn (latlon):~ > g.region rast=grid -p<br>
> projection: 3 (Latitude-Longitude)<br>
> zone:       0<br>
> datum:      wgs84<br>
> ellipsoid:  a=6378137 es=0.006694379990141316<br>
> north:      37:22:30N<br>
> south:      39S<br>
> west:       20:12:30W<br>
> east:       57:40E<br>
> nsres:      0:02:30<br>
> ewres:      0:02:30<br>
> rows:       1833<br>
> cols:       1869<br>
> cells:      3425877<br>
<br>
</span>Can you please post with "-g" which is decimal?<br></blockquote></div></div><div><br>n=37.375<br>s=-34.875<br>w=-17.5416666666667<br>e=51.4166666666667<br>nsres=0.0416666666666667<br>ewres=0.0416666666666667<br>rows=1734<br>cols=1655<br>cells=2869770<br> </div><span class=""><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<span><br>
>><br>
>><br>
>> > v.what.rast map=gbif raster=grid column=cell2arcm<br>
>> > Reading features from vector map...<br>
>> >  100%<br>
>> > ERROR: Reading raster map <grid@> request for row 1830 is outside region<br>
>><br>
>> ... and also the boundary box of gbif?<br>
><br>
> Projection: Latitude-Longitude<br>
><br>
> N: 37:19:58.08N<br>
> S: 38:52:30S<br>
> E: 57:38:24E<br>
> W: 20:02:49.92W<br>
<br>
</span>Also here, can you please post with "-g" which is decimal?<br></blockquote></span><div><br>n=37.3328<br>s=-38.875<br>w=-20.0472<br>e=57.64<br>nsres=0.0416663750683434<br>ewres=0.0416776824034335<br>rows=1829<br>cols=1864<br>cells=3409256<br></div><span class=""><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
The error message is generated in<br>
lib/raster/get_row.c<br>
<br>
static int compute_window_row(int fd, int row, int *cellRow)<br>
{<br>
    struct fileinfo *fcb = &R__.fileinfo[fd];<br>
    double f;<br>
    int r;<br>
<br>
    /* check for row in window */<br>
    if (row < 0 || row >= R__.rd_window.rows) {<br>
        G_fatal_error(_("Reading raster map <%s@%s> request for row %d<br>
is outside region"),<br>
                      fcb->name, fcb->mapset, row);<br>
    }<br>
<br>
Perhaps a vector point just hits the extents of the computational region or map?<br></blockquote><div><br><br></div></span><div>You are right that the extends do not match completely. But I was under the impression that points outside the region simply got skipped? For example, if I reduce the region to a small area, and run v.what.rast, I get the warning message "WARNING: 2287133 points outside current region were skipped", after which the it happily continues to run updating the points that are within the region bounds, as per below.<br><br>WARNING: 2287133 points outside current region were skipped<br>Update vector attributes...<br> 100%<br>v.what.rast complete. 57621 records updated.<br><br></div><div></div></div></div></div></blockquote><div><br></div><div>Sorry, I wasn't completely awake yet. It is about the extend of the raster layer to be sampled being smaller than the region... and yes, that was the case. Reducing the region to fit within the bounds of the raster layer solved the problem. Seems obvious now I know the error I made, sorry. Thanks a lot for your help, much appreciated!<br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<span><font color="#888888"><br>
Markus<br>
</font></span></blockquote></div><br></div></div>
</blockquote></div><br></div></div>