<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">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 class=""><<a href="mailto:p.vanbreugel@gmail.com">p.vanbreugel@gmail.com</a>> wrote:<br>
> On Wed, Oct 8, 2014 at 9:21 AM, Markus Neteler <<a href="mailto:neteler@osgeo.org">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">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><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><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<span class=""><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><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><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><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><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<span class=""><font color="#888888"><br>
Markus<br>
</font></span></blockquote></div><br></div></div>