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 <tsvi...@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 <mailto: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 <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@lists.osgeo.org <mailto:grass-user@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
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user