Hi Luis, From the error message, it seems v.to.rast is failing, which may be due to a lack of resources - as Micha suggests.
Can you try a) v.to.rast input=PSU output=test use=cat if that does not work, we at least know it is not v.rast.stats and we might get a more specific error message... Not sure if the memory option in v.to.rast would have an effect here... b) try to subdivide your data into chunks / tiles (polygons within tiles)? Cheers Stefan -----Original Message----- From: grass-user <[email protected]> On Behalf Of Micha Silver Sent: mandag 1. februar 2021 21:32 To: Luí s Moreira de Sousa <[email protected]>; GRASS user list <[email protected]> Subject: Re: [GRASS-user] v.rast.stats error: "Unable to seek" 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://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fprotonmail.com%2F&data=04%7C01%7C%7C3a73488d58cd411cbc0f08d8c6f07342%7C6cef373021314901831055b3abf02c73%7C0%7C0%7C637478083673619175%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=4t7uF0b7BR9kcAqAReOPlaM6z8L7fwAJ94zMC66lzTc%3D&reserved=0> > Secure Email. > > ‐‐‐‐‐‐‐ Original Message ‐‐‐‐‐‐‐ > On Friday, January 29, 2021 3:31 PM, Micha Silver <[email protected]> > wrote: > >> >> >> On Fri, Jan 29, 2021 at 3:53 PM Luís Moreira de Sousa via grass-user >> <[email protected] <mailto:[email protected]>> 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 >> <https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fv.i >> nfo%2F&data=04%7C01%7C%7C3a73488d58cd411cbc0f08d8c6f07342%7C6cef3 >> 73021314901831055b3abf02c73%7C0%7C0%7C637478083673619175%7CUnknown%7C >> TWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXV >> CI6Mn0%3D%7C3000&sdata=Lc8A9SfmNwQKAMt3TTXYh%2FYBNurG6kUYAubGcxKh >> 3hw%3D&reserved=0> PSU r.info >> <https://eur01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fr.i >> nfo%2F&data=04%7C01%7C%7C3a73488d58cd411cbc0f08d8c6f07342%7C6cef3 >> 73021314901831055b3abf02c73%7C0%7C0%7C637478083673619175%7CUnknown%7C >> TWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXV >> CI6Mn0%3D%7C3000&sdata=GyWpg6cTSQfOc0LOFHKvEzLa4lxJXbinsTJlPm0lbt >> Y%3D&reserved=0> MAL_Mode_5x5 >> >> >> >> s4a=# select count(*) from mal.psu; >> count >> ---------- >> 19468516 >> (1 row) >> >> >> >> >> -- >> Luís >> >> _______________________________________________ >> grass-user mailing list >> [email protected] <mailto:[email protected]> >> >> https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.osgeo.org%2Fmailman%2Flistinfo%2Fgrass-user&data=04%7C01%7C%7C3a73488d58cd411cbc0f08d8c6f07342%7C6cef373021314901831055b3abf02c73%7C0%7C0%7C637478083673629172%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=2P01wK1Ak3X0QzZsMNLQa%2FLrQoVWrd8C7FgNezSRYRM%3D&reserved=0 >> >> <https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fli >> sts.osgeo.org%2Fmailman%2Flistinfo%2Fgrass-user&data=04%7C01%7C%7 >> C3a73488d58cd411cbc0f08d8c6f07342%7C6cef373021314901831055b3abf02c73% >> 7C0%7C0%7C637478083673629172%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAw >> MDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=2 >> P01wK1Ak3X0QzZsMNLQa%2FLrQoVWrd8C7FgNezSRYRM%3D&reserved=0> >> >> >> >> -- >> 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://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Forcid.org%2F0000-0002-1128-1325&data=04%7C01%7C%7C3a73488d58cd411cbc0f08d8c6f07342%7C6cef373021314901831055b3abf02c73%7C0%7C0%7C637478083673629172%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=PuJmTChNnqdqy%2F3CrSb3FbfoIMprgxLPwf%2B0j86VS8s%3D&reserved=0 _______________________________________________ grass-user mailing list [email protected] https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.osgeo.org%2Fmailman%2Flistinfo%2Fgrass-user&data=04%7C01%7C%7C3a73488d58cd411cbc0f08d8c6f07342%7C6cef373021314901831055b3abf02c73%7C0%7C0%7C637478083673629172%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000&sdata=2P01wK1Ak3X0QzZsMNLQa%2FLrQoVWrd8C7FgNezSRYRM%3D&reserved=0 _______________________________________________ grass-user mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/grass-user
