[GRASS-user] v.rast.stats error: "Unable to seek"
Veronica Andreo
veroandreo at gmail.com
Mon Feb 1 07:52:23 PST 2021
Does the example in v.rast.stats manual [1] work for you, Luis? Is your
polygon vector map "healthy"? The error you get seems to be while
converting the vector to raster, hence I wonder if boundaries of all areas
in your vector are right.
my 0.02 cents
Vero
[1] https://grass.osgeo.org/grass78/manuals/v.rast.stats.html
El lun, 1 feb 2021 a las 16:31, Luís Moreira de Sousa via grass-user (<
grass-user at lists.osgeo.org>) escribió:
> Hi Micha, thank you for replying.
>
> I recreated the PSU layer in another mapset with the default SQLite
> back-end, v.rast.stats fails with the exact same error. So this is not
> related to the back-end.
>
> Below are the outputs you requested. Thank you.
>
> > v.info PSU
>
> +----------------------------------------------------------------------------+
> | Name:
> PSU |
> | Mapset:
> MAL |
> | Location:
> S4AHomolosine |
> | Database:
> /home/duque004/Work/GRASSDATA |
> |
> Title: |
> | Map scale:
> 1:1 |
> | Name of creator:
> duque004 |
> |
> Organization: |
> | Source date: Fri Jan 29 08:13:08
> 2021 |
> | Timestamp (first layer):
> none |
>
> |----------------------------------------------------------------------------|
> | Map format:
> native |
>
> |----------------------------------------------------------------------------|
> | Type of map: vector (level:
> 2) |
> |
> |
> | Number of points: 0 Number of centroids:
> 19468516 |
> | Number of lines: 0 Number of boundaries:
> 38945853 |
> | Number of areas: 19468516 Number of islands:
> 1 |
> |
> |
> | Map is 3D:
> No |
> | Number of dblinks:
> 1 |
> |
> |
> | Projection:
> unknown |
> |
> |
> | N: 4289569.27224353 S:
> -4452930.72775647 |
> | E: 6679312.25515029 W:
> -2226387.74484971 |
> |
> |
> | Digitization threshold:
> 0 |
> |
> Comment: |
> |
> |
>
> +----------------------------------------------------------------------------+
>
> > r.info MAL_Mode_5x5
>
> +----------------------------------------------------------------------------+
> | Map: MAL_Mode_5x5 Date: Thu Jan 28 09:08:14
> 2021 |
> | Mapset: MAL Login of Creator:
> duque004 |
> | Location:
> S4AHomolosine |
> | DataBase:
> /home/duque004/Work/GRASSDATA |
> | Title: 5x5 neighborhood: mode of
> MAL_AFRICA |
> | Timestamp:
> none |
>
> |----------------------------------------------------------------------------|
> |
> |
> | Type of Map: raster Number of Categories:
> 0 |
> | Data Type:
> CELL |
> | Rows:
> 87425 |
> | Columns:
> 89057 |
> | Total Cells:
> 7785808225 |
> | Projection:
> unknown |
> | N: 4289569.27224353 S: -4452930.72775647 Res:
> 100 |
> | E: 6679312.25515029 W: -2226387.74484971 Res:
> 100 |
> | Range of data: min = 0 max =
> 1 |
> |
> |
> | Data
> Description: |
> | generated by
> r.neighbors |
> |
> |
> |
> Comments: |
> | r.neighbors input="MAL_AFRICA" output="MAL_Mode_5x5" method="mode"
> s\ |
> |
> ize=5 |
> |
> |
>
> +----------------------------------------------------------------------------+
>
>
> --
> Luís
>
> Sent with ProtonMail <https://protonmail.com> Secure Email.
>
> ‐‐‐‐‐‐‐ Original Message ‐‐‐‐‐‐‐
> On Friday, January 29, 2021 3:31 PM, Micha Silver <tsvibar at gmail.com>
> wrote:
>
>
>
> On Fri, Jan 29, 2021 at 3:53 PM Luís Moreira de Sousa via grass-user <
> grass-user at lists.osgeo.org> wrote:
>
>> Dear all,
>>
>> I am getting the error "Unable to seek" with v.stats.error. There is not
>> much information that could point the cause, just a warning saying that
>> some data base files are not found. I checked the database connection and
>> everything looks in order (see below). Any hints on what may be causing
>> this error?
>>
>> Thank you.
>>
>> > v.rast.stats map=PSU raster=MAL_Mode_5x5 method=number
>> column_prefix=mal
>> ERROR: Unable to seek: Invalid argument
>> ERROR: An error occurred while converting vector to raster
>> WARNING: No data base element files found
>>
>> > db.connect -p
>> driver: pg
>> database: s4a
>> schema: mal
>> group:
>>
>> > psql -d s4a -p 5434
>> psql (12.5 (Ubuntu 12.5-0ubuntu0.20.04.1))
>> Type "help" for help.
>>
>> s4a=# \dt mal.*
>> List of relations
>> Schema | Name | Type | Owner
>> --------+--------------------------+-------+----------
>> mal | psu | table | duque004
>> (1 rows)
>>
>>
> Could it be simply CAPS in the vector name vs small letters in the
> Postgres table?
> Can you post the outputs from:
> v.info PSU
> r.info MAL_Mode_5x5
>
>
>
>> s4a=# select count(*) from mal.psu;
>> count
>> ----------
>> 19468516
>> (1 row)
>>
>>
>>
>>
>> --
>> Luís
>>
>> _______________________________________________
>> grass-user mailing list
>> grass-user at lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/grass-user
>>
>
>
> --
> Micha Silver
> Ben Gurion Univ
> Sde-Boker Remote Sensing Lab
> cell: +972 (52) 3665918
>
>
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20210201/9988a09b/attachment.html>
More information about the grass-user
mailing list