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.
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 Secure Email.
‐‐‐‐‐‐‐ Original Message ‐‐‐‐‐‐‐
On Friday, January 29, 2021 3:31 PM, Micha Silver <tsvibar@gmail.com> wrote:
On Fri, Jan 29, 2021 at 3:53 PM Luís Moreira de Sousa via grass-user <grass-user@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@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@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user