[GRASS-user] v.rast.stats error: "Unable to seek"
Micha Silver
tsvibar at gmail.com
Mon Feb 1 12:32:00 PST 2021
Hi Luis:
On 2/1/21 5:30 PM, Luí s Moreira de Sousa wrote:
> 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 |
> | |
> +----------------------------------------------------------------------------+
>
So the vector contains about 20 million polygons, and the raster about 8
billion cells. (~=32GB ?)
Does you computer have enough muscle for this?
>
> --
> 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 <mailto: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 <http://v.info> PSU
>> r.info <http://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 <mailto:grass-user at lists.osgeo.org>
>> https://lists.osgeo.org/mailman/listinfo/grass-user
>> <https://lists.osgeo.org/mailman/listinfo/grass-user>
>>
>>
>>
>> --
>> Micha Silver
>> Ben Gurion Univ
>> Sde-Boker Remote Sensing Lab
>> cell: +972 (52) 3665918
>
--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918
https://orcid.org/0000-0002-1128-1325
More information about the grass-user
mailing list