[GRASS-user] Which colour palette is currently applied to a raster
Hi Users, How can I tell which colour palette set is currently applied to a raster in Grass 7.8.3? nothing is stated in metadata and checked in 'PERMANENT/colr' for colour palette 'name' but nothing. Many thanks, *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Lost attribute connection - grass 7.6.1
Hi Users, I often create vector data using 'r.to.vect' which vectorises with no issues but If I then export it to PostGIS using 'v.out.postgis' the first layer will export fully but consecutive layers fail with the following error:- GRASS_INFO_WARNING(4924,1): Table not found in database using driver GRASS_INFO_END(4924,1) After which the table does not open again in grass. I've checked to ensure both 'sqlite' and 'pg' connections are present. When I get the above error I am connected to 'pg' db.tables -p (shows me i'm connecetd to 'pg') db.connect -p or -g driver: pg database: ## schema: grass_gis group: postgis_reader (Tue Aug 06 15:56:34 db.connect -g This is probably something obvious i'm missing!! Thanks, Paul Shapley --- Traceback (most recent call last): File "C:\OSGeo4W\apps\grass\grass76\gui\wxpython\lmgr\frame.py", line 2018, in OnShowAttributeTable selection=selection) File "C:\OSGeo4W\apps\grass\grass76\gui\wxpython\dbmgr\man ager.py", line 74, in __init__ **kwargs) File "C:\OSGeo4W\apps\grass\grass76\gui\wxpython\dbmgr\base.py", line 749, in __init__ self.dbMgrData['mapDBInfo'] = VectorDBInfo(self.dbMgrData['vectName']) File "C:\OSGeo4W\apps\grass\grass76\gui\wxpython\dbmgr\vinfo.py", line 87, in __init__ VectorDBInfoBase.__init__(self, map) File "C:\OSGeo4W\apps\grass\grass76\gui\wxpython\gui_core\ gselect.py", line 803, in __init__ self._DescribeTables() # -> self.tables File "C:\OSGeo4W\apps\grass\grass76\gui\wxpython\gui_core\ gselect.py", line 835, in _DescribeTables database=self.layers[layer]["database"])['cols']: File "C:\OSGeo4W\apps\grass\grass76\etc\python\grass\scrip t\db.py", line 43, in db_describe s = read_command('db.describe', flags='c', table=table, **args) File "C:\OSGeo4W\apps\grass\grass76\etc\python\grass\scrip t\core.py", line 478, in read_command return handle_errors(returncode, stdout, args, kwargs) File "C:\OSGeo4W\apps\grass\grass76\etc\python\grass\scrip t\core.py", line 334, in handle_errors returncode=returncode) grass.exceptions . CalledModuleError : Module run None ['db.describe', '-c', 'table=grass_gis.nameofmylayer', 'driver=pg', 'database=webgis'] ended with error Process ended with non-zero return code 1. See errors in the (error) output. ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.external
Hi Users, I've upgraded to Grass 7.6.0 (r 73958) and i've tried using 'v.external' to link to our Postgresql/PostGIS tables in the usual way but when I try to link no tables appear in the list. I've selected 'Database' 'ODBC' then entered PG:dbname=MyDatabase then I get the error:- *No Layers Selected. Operation Cancelled* If i enter 'db.tables -p' i can see all the tables and schemas listed and the 'VAR' file shows the connections so i think it is connecting OK. Thanks, *Paul J. Shapley* ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Image signature preview
Hi Users, I want to isolate a single signature created from an area polygon with i.class/g.gui.iclass. I can do this and preview the result i need but what I would like to do is create a full region map from the 'SIG' in order to display just the one category type of land (isolated in the training area and applied to the remaining map). How could I achieve this in Grass 7.4.1 Thanks, -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r,in.gdal error when importing a jpeg2000 image
Hi Nikos, Would i be able to use the 'i.segment' module on subsets of the image if i use 'r.external' instead of importing? Can you point me to where the differences are between importing and linking? Thank You, Paul Shapley On 18 July 2018 at 13:48, Nikos Alexandris wrote: > * Paul Shapley [2018-07-18 09:32:56 +0100]: > > Hi Users, >> >> Trying to import an 8.7gb (jpeg2000) image into Grass 7.4.1 64 bit. I get >> the error below. Tested the image in another OS GIS package and it opens >> without issue. Not sure why Grass wont import it. >> >> Thanks everyone! >> >> Paul Shapley >> >> > If it's a driver issue (like Daniel mentions in his post) and GDAL can > read it, do you really need to import it? You know that you can link to > it via r.external. > > Nikos > > [..] > -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r,in.gdal error when importing a jpeg2000 image
Hi Users, Trying to import an 8.7gb (jpeg2000) image into Grass 7.4.1 64 bit. I get the error below. Tested the image in another OS GIS package and it opens without issue. Not sure why Grass wont import it. Thanks everyone! Paul Shapley (Wed Jul 18 09:19:16 2018) r.in.gdal -o input=J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2 output=aerial_25cm_master memory=3000 Over-riding projection check Importing 3 raster bands... Importing raster map ... ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0 ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0 ERROR 1: Marker is not compliant with its position ERROR 1: opj_decode() failed ERROR 1: J:\Spatial Data\Aerial\APGB_received_20180416\aerial_25cm.jp2, band 1: IReadBlock failed at X offset 0, Y offset 0
Re: [GRASS-user] Switching databases for holding attributes in Grass 7.4.1
Hi Marcus, In the end i just removed the details in the 'VAR' file and saved it and it now defaults back to sqlite. Many thanks everyone for your help as always. On 3 July 2018 at 07:48, Markus Metz wrote: > > > On Tue, Jul 3, 2018 at 8:09 AM, Moritz Lennert < > mlenn...@club.worldonline.be> wrote: > > > > > > > > Am 2. Juli 2018 21:06:15 MESZ schrieb Markus Metz < > markus.metz.gisw...@gmail.com>: > > >On Mon, Jul 2, 2018 at 4:54 PM, Moritz Lennert > > > > > >wrote: > > >> > > >> Le Mon, 2 Jul 2018 16:36:30 +0200, > > >> Markus Metz a écrit : > > >> > > >> > On Mon, Jul 2, 2018 at 4:04 PM, Paul Shapley > > >> > wrote: > > >> > > > > >> > > Hi Markus, > > >> > > > > >> > > Used 'db.connect' driver=sqlite > > >> > database=$GISDBASE/$LOCATION_NAME/$MAPSET/sqlite/sqlite.db > > >> > > > > >> > > I had to delete 'database schema=grass_gis' > > >> > > > > >> > > I had to delete 'default group=postgis_reader' > > >> > > > >> > That seems to be a bug in db.connect, apparently it is not possible > > >to > > >> > reset schema and group. Can you please open a ticket? > > >> > > >> AFAIR, schema= and group= (i.e. without any value) should reset these > > >> parameters. > > > > > >That means no answer to schema and group would be a reasonable default > > >answer in order to easily change the current connection parameters? > > > > No, should was the wrong word: does would have been better. AFAIR, > explicitly setting them to no answer already did reset them. > > True. However, the manual says > > *schema*=*name* Database schema Do not use this option if schemas are not > supported by driver/database server > Markus M > > > > > Moritz > > > > > -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Switching databases for holding attributes in Grass 7.4.1
Hi Markus, Used 'db.connect' driver=sqlite database=$GISDBASE/$LOCATION_NAME/$MAPSET/sqlite/sqlite.db *I had to delete* 'database schema=grass_gis' *I had to delete* 'default group=postgis_reader' Then 'RUN' output:- (Mon Jul 02 14:59:21 2018) db.connect (Mon Jul 02 14:59:22 2018) Command finished (1 sec) now attempted 'r.to.vect' r.to.vect input=APGB_aerial_2_i_segment_drained_peat_with_exposed_haggs@PERMANENT output=test1 type=area DBMI-SQLite driver error: Error in sqlite3_prepare(): unknown database grass_gis DBMI-SQLite driver error: Error in sqlite3_prepare(): unknown database grass_gis ERROR: Unable to create table: create table grass_gis.test1 ( cat integer, value integer, label varchar(10)) WARNING: Table linked to vector map does not exist Repeat as above!! Used 'db.connect' driver=sqlite database=$GISDBASE/$LOCATION_NAME/$MAPSET/sqlite/sqlite.db *I had to delete* 'database schema=grass_gis' *I had to delete* 'default group=postgis_reader' Then 'RUN' On 2 July 2018 at 14:41, Markus Metz wrote: > > > On Mon, Jul 2, 2018 at 3:25 PM, Paul Shapley wrote: > > > > Hi Markus, > > > > I used 'db.connect' > > which flags/options did you use? > > > followed by 'db.login' > > you don't need db.login for sqlite > > > > Please see below:- > > > > db.connect -p > > > driver: sqlite > > database: C:\\TempLocation\PERMANENT\sqlite\sqlite.db > > schema: grass_gis > > group: postgis_reader > > (Mon Jul 02 14:17:04 2018) Command finished (0 sec) > > > > db.connect -d > > > Default driver / database set to: > > driver: sqlite > > database: C:\\TempLocation\PERMANENT\sqlite\sqlite.db > > (Mon Jul 02 14:19:40 2018) Command finished (0 sec) > > > > Still holding 'schema' and 'group' details from PostGIS. Shouldn't these > have been dropped when changing to sqlite which is perhaps why the > connection is failing? > > yes. > > It seems that db.connect -g is better than db.connect -p to verify > connection settings, e.g. I get in the North Carolina sample dataset > > > db.connect -g > driver=sqlite > database=$GISDBASE/$LOCATION_NAME/$MAPSET/sqlite/sqlite.db > schema= > group= > > The database is not a real path but contains GRASS variables that are > evaluated on the fly. This is the default connection. > > Markus M > > > > > On 2 July 2018 at 14:10, Markus Metz > wrote: > >> > >> > >> > >> On Mon, Jul 2, 2018 at 2:00 PM, Paul Shapley > wrote: > >> > > >> > Hi Users, > >> > > >> > I want to use 'r.to.vect' but it seems to default to a 'Postgis' > table that has no geometry column just a table with a 'cat' id. I would > like to export to 'sqlite'. I have both postgres an sqlite login details > stored in 'db.login'. How do i switch between the two or at least make > sqlite the default? If I change settings with 'db.connect' to 'sqlite' > >> > >> what is the exact command you used to set the default connection to > sqlite? > >> > >> You can also try db.connect -d, otherwise you might need to specify the > default database to be used. > >> > >> In any case, please verify first with db.connect -p to see the > connection details used when creating a new vector. > >> > >> Markus M > >> > >> > it fails (error below) but still keeps holding onto the 'grass_gis' > schema form the Postgis connection! > >> > > >> > r.to.vect input=APGB_aerial_2_i_segment_drained_peat_with_exposed_ > haggs@PERMANENT output=test1 type=area > >> > DBMI-SQLite driver error: > >> > Error in sqlite3_prepare(): > >> > unknown database grass_gis > >> > DBMI-SQLite driver error: > >> > Error in sqlite3_prepare(): > >> > unknown database grass_gis > >> > ERROR: Unable to create table: create table grass_gis.test1 ( cat > integer, value integer, label varchar(10)) > >> > WARNING: Table linked to vector map does > not exist > >> > > >> > Is there a preferred option to change databases more cleanly. > >> > > >> > Many thanks, > >> > > >> > > >> > -- > >> > Paul J. Shapley MSc CGeog (GIS) FRGS > >> > > >> > > >> > ___ > >> > grass-user mailing list > >> > grass-user@lists.osgeo.org > >> > https://lists.osgeo.org/mailman/listinfo/grass-user > >> > > > > > > > > -- > > Paul J. Shapley MSc CGeog (GIS) FRGS > > > -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Switching databases for holding attributes in Grass 7.4.1
Hi Markus, I used 'db.connect' followed by 'db.login' Please see below:- db.connect -p driver: sqlite database: C:\\TempLocation\PERMANENT\sqlite\sqlite.db schema: grass_gis group: postgis_reader (Mon Jul 02 14:17:04 2018) Command finished (0 sec) db.connect -d Default driver / database set to: driver: sqlite database: C:\\TempLocation\PERMANENT\sqlite\sqlite.db (Mon Jul 02 14:19:40 2018) Command finished (0 sec) Still holding 'schema' and 'group' details from PostGIS. Shouldn't these have been dropped when changing to sqlite which is perhaps why the connection is failing? On 2 July 2018 at 14:10, Markus Metz wrote: > > > On Mon, Jul 2, 2018 at 2:00 PM, Paul Shapley wrote: > > > > Hi Users, > > > > I want to use 'r.to.vect' but it seems to default to a 'Postgis' table > that has no geometry column just a table with a 'cat' id. I would like to > export to 'sqlite'. I have both postgres an sqlite login details stored in > 'db.login'. How do i switch between the two or at least make sqlite the > default? If I change settings with 'db.connect' to 'sqlite' > > what is the exact command you used to set the default connection to sqlite? > > You can also try db.connect -d, otherwise you might need to specify the > default database to be used. > > In any case, please verify first with db.connect -p to see the connection > details used when creating a new vector. > > Markus M > > > it fails (error below) but still keeps holding onto the 'grass_gis' > schema form the Postgis connection! > > > > r.to.vect input=APGB_aerial_2_i_segment_drained_peat_with_exposed_ > haggs@PERMANENT output=test1 type=area > > DBMI-SQLite driver error: > > Error in sqlite3_prepare(): > > unknown database grass_gis > > DBMI-SQLite driver error: > > Error in sqlite3_prepare(): > > unknown database grass_gis > > ERROR: Unable to create table: create table grass_gis.test1 ( cat > integer, value integer, label varchar(10)) > > WARNING: Table linked to vector map does not > exist > > > > Is there a preferred option to change databases more cleanly. > > > > Many thanks, > > > > > > -- > > Paul J. Shapley MSc CGeog (GIS) FRGS > > > > > > ___ > > grass-user mailing list > > grass-user@lists.osgeo.org > > https://lists.osgeo.org/mailman/listinfo/grass-user > > -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Switching databases for holding attributes in Grass 7.4.1
Hi Users, I want to use 'r.to.vect' but it seems to default to a 'Postgis' table that has no geometry column just a table with a 'cat' id. I would like to export to 'sqlite'. I have both postgres an sqlite login details stored in 'db.login'. How do i switch between the two or at least make sqlite the default? If I change settings with 'db.connect' to 'sqlite' it fails (error below) but still keeps holding onto the 'grass_gis' schema form the Postgis connection! r.to.vect input=APGB_aerial_2_i_segment_drained_peat_with_exposed_haggs@PERMANENT output=test1 type=area DBMI-SQLite driver error: Error in sqlite3_prepare(): unknown database grass_gis DBMI-SQLite driver error: Error in sqlite3_prepare(): unknown database grass_gis ERROR: Unable to create table: create table grass_gis.test1 ( cat integer, value integer, label varchar(10)) WARNING: Table linked to vector map does not exist Is there a preferred option to change databases more cleanly. Many thanks, -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.out.postgis
Hi Moritz, My apologies... it appears that the result of 'r.to.vect was indeed sending the attributes to postgres but i was expecting them to appear in sqlite had i 'refreshed' the posgres schema i would have seen them. Best Wishes and thank you, Paul Shapley On 29 June 2018 at 14:50, Moritz Lennert wrote: > On 29/06/18 12:38, Paul Shapley wrote: > >> Hi Moritz, >> >> Here is the result:- >> >> v.db.select APGB_aerial_2_i_segment_drained_peat_with_exposed_haggs_ >> vector >> DBMI-PostgreSQL driver error: >> Unable to select: >> SELECT * FROM grass_gis.APGB_aerial_2_i_segment_drained_peat >> _with_exposed_haggs_vector >> ERROR: relation "grass_gis.apgb_aerial_2_i_segment_drained_ >> peat_with_exposed_haggs_vector" does not exist >> LINE 1: SELECT * FROM >> grass_gis.APGB_aerial_2_i_segment_drained_peat... >>^ >> DBMI-PostgreSQL driver error: >> Unable to select: >> SELECT * FROM grass_gis.APGB_aerial_2_i_segment_drained_peat >> _with_exposed_haggs_vector >> ERROR: relation "grass_gis.apgb_aerial_2_i_segment_drained_ >> peat_with_exposed_haggs_vector" does not exist >> LINE 1: SELECT * FROM >> grass_gis.APGB_aerial_2_i_segment_drained_peat... >>^ >> ERROR: Unable to open select cursor >> (Fri Jun 29 11:24:18 2018) Command finished (0 sec) >> > > So, apparently your vector map is linked to a table that does not exist. > Or it exists but not in the schema grass_gis. You will have to correct the > link between the map and its attribute table before you can export it. > > Moritz > > -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.out.postgis
Hi Moritz, Here is the result:- v.db.select APGB_aerial_2_i_segment_drained_peat_with_exposed_haggs_vector DBMI-PostgreSQL driver error: Unable to select: SELECT * FROM grass_gis.APGB_aerial_2_i_segment_drained_peat _with_exposed_haggs_vector ERROR: relation "grass_gis.apgb_aerial_2_i_segment_drained_ peat_with_exposed_haggs_vector" does not exist LINE 1: SELECT * FROM grass_gis.APGB_aerial_2_i_segment_drained_peat... ^ DBMI-PostgreSQL driver error: Unable to select: SELECT * FROM grass_gis.APGB_aerial_2_i_segment_drained_peat _with_exposed_haggs_vector ERROR: relation "grass_gis.apgb_aerial_2_i_segment_drained_ peat_with_exposed_haggs_vector" does not exist LINE 1: SELECT * FROM grass_gis.APGB_aerial_2_i_segment_drained_peat... ^ ERROR: Unable to open select cursor (Fri Jun 29 11:24:18 2018) Command finished (0 sec) - I have two databases that i use in grass, the (sqlite.db) for vector data (e.g using r.to.vect) then i usually export the results to postgresql/postgis (webgis). The 'db.login' file holds the PostGIS details above the sqlite details. If i try to open the attributes (of sqlite) from within Grass it is trying to find them in postgres/postgis (webgis). see cmdline error:- GRASS_INFO_WARNING<1408,1>:Table not found in database using driver GRASS_INFO_END<1408,1> Is the order of the connection details in 'pg.login' important? Paul Shapley On 29 June 2018 at 11:17, Moritz Lennert wrote: > On 29/06/18 10:39, Paul Shapley wrote: > >> Hi Markus, >> >> OK here is the result of v.db.connect -p >> >> v.db.connect -p APGB_aerial_2_i_segment_drained_peat_with_exposed_haggs_ >> vector >> Vector map >> is connected by: >> layer >> <1/grass_gis.APGB_aerial_2_i_segment_drained_peat_with_exposed_haggs_vector> >> table >> >> in database through driver with key >> > > So it's PostgreSQL as attribute database backend. > > What is the output of > > v.db.select APGB_aerial_2_i_segment_drained_peat_with_exposed_haggs_vector > > ? > > Moritz > > >> I recall some years ago you had to delete the 'db.login' file and >> recreate it but this would lose the connection to the vector layer >> attribute table (sqlite.db). >> >> Thanks, >> >> Paul Shapley >> >> On 28 June 2018 at 13:22, Markus Metz > <mailto:markus.metz.gisw...@gmail.com>> wrote: >> >> >> >> On Thu, Jun 28, 2018 at 1:40 PM, Moritz Lennert >> mailto:mlenn...@club.worldonline.be>> >> wrote: >> > >> > On 28/06/18 09:59, Paul Shapley wrote: >> >> >> >> Hi, >> >> >> >> I am trying to export a vector map to PostGIS. Has anyone come >> across the 'Unable to describe table' error? I suspect my config >> could be wrong. >> >> >> >> v.out.postgis --overwrite >> input=APGB_aerial_2_i_segment_drained_peat_with_exposed_hagg >> s_vct@PERMANENT >> output=PG:dbname=webgis >> output_layer=APGB_aerial_2_i_segment_drained_peat_with_expos >> ed_haggs_vct >> options=FID=cat, GEOMETRY_NAME=wkb_geometry, SPATIAL_INDEX=YES, >> PRIMARY_KEY=YES >> >> DBMI-PostgreSQL driver error: >> >> select * from grass_gis.APGB_aerial_2_i_segment_drained_peat >> >> _with_exposed_haggs_vct where 1 = 0 >> >> ERROR: relation "grass_gis.apgb_aerial_2_i_segment_drained_ >> >> peat_with_exposed_haggs_vct" does not exist >> >> LINE 1: select * from >> >> grass_gis.APGB_aerial_2_i_segment_drained_peat... >> >>^ >> >> DBMI-PostgreSQL driver error: >> >> select * from grass_gis.APGB_aerial_2_i_segment_drained_peat >> >> _with_exposed_haggs_vct where 1 = 0 >> >> ERROR: relation "grass_gis.apgb_aerial_2_i_segment_drained_ >> >> peat_with_exposed_haggs_vct" does not exist >> >> LINE 1: select * from >> >> grass_gis.APGB_aerial_2_i_segment_drained_peat... >> >>^ >> >> *ERROR: Unable to describe table* >> > _haggs_vct> >> > >> > >> > It looks like it is trying to read the input attrubute table from >> postgresql (schema grass_gis). Do you use PostgreSQL as attribute >> table backend ? Could you provide the output of db.connect -p i
Re: [GRASS-user] v.out.postgis
Hi Markus, OK here is the result of v.db.connect -p v.db.connect -p APGB_aerial_2_i_segment_drained_peat_with_exposed_haggs_vector Vector map is connected by: layer <1/grass_gis.APGB_aerial_2_i_segment_drained_peat_with_exposed_haggs_vector> table in database through driver with key I recall some years ago you had to delete the 'db.login' file and recreate it but this would lose the connection to the vector layer attribute table (sqlite.db). Thanks, Paul Shapley On 28 June 2018 at 13:22, Markus Metz wrote: > > > On Thu, Jun 28, 2018 at 1:40 PM, Moritz Lennert < > mlenn...@club.worldonline.be> wrote: > > > > On 28/06/18 09:59, Paul Shapley wrote: > >> > >> Hi, > >> > >> I am trying to export a vector map to PostGIS. Has anyone come across > the 'Unable to describe table' error? I suspect my config could be wrong. > >> > >> v.out.postgis --overwrite input=APGB_aerial_2_i_segment_ > drained_peat_with_exposed_haggs_vct@PERMANENT output=PG:dbname=webgis > output_layer=APGB_aerial_2_i_segment_drained_peat_with_exposed_haggs_vct > options=FID=cat, GEOMETRY_NAME=wkb_geometry, SPATIAL_INDEX=YES, > PRIMARY_KEY=YES > >> DBMI-PostgreSQL driver error: > >> select * from grass_gis.APGB_aerial_2_i_segment_drained_peat > >> _with_exposed_haggs_vct where 1 = 0 > >> ERROR: relation "grass_gis.apgb_aerial_2_i_segment_drained_ > >> peat_with_exposed_haggs_vct" does not exist > >> LINE 1: select * from > >> grass_gis.APGB_aerial_2_i_segment_drained_peat... > >>^ > >> DBMI-PostgreSQL driver error: > >> select * from grass_gis.APGB_aerial_2_i_segment_drained_peat > >> _with_exposed_haggs_vct where 1 = 0 > >> ERROR: relation "grass_gis.apgb_aerial_2_i_segment_drained_ > >> peat_with_exposed_haggs_vct" does not exist > >> LINE 1: select * from > >> grass_gis.APGB_aerial_2_i_segment_drained_peat... > >>^ > >> *ERROR: Unable to describe table* segment_drained_peat_with_exposed_haggs_vct> > > > > > > It looks like it is trying to read the input attrubute table from > postgresql (schema grass_gis). Do you use PostgreSQL as attribute table > backend ? Could you provide the output of db.connect -p in that mapset ? > > ore precise would be > v.db.connect -p APGB_aerial_2_i_segment_drained_peat_with_exposed_ > haggs_vct@PERMANENT > > because there may be a bug in v.out.postgis, not properly reading the > database connection of the input vector. > > Markus Metz > > > > > Moritz > > > > ___ > > grass-user mailing list > > grass-user@lists.osgeo.org > > https://lists.osgeo.org/mailman/listinfo/grass-user > > -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.out.postgis
Hi, I am trying to export a vector map to PostGIS. Has anyone come across the 'Unable to describe table' error? I suspect my config could be wrong. v.out.postgis --overwrite input=APGB_aerial_2_i_segment_drained_peat_with_exposed_haggs_vct@PERMANENT output=PG:dbname=webgis output_layer=APGB_aerial_2_i_segment_drained_peat_with_exposed_haggs_vct options=FID=cat, GEOMETRY_NAME=wkb_geometry, SPATIAL_INDEX=YES, PRIMARY_KEY=YES DBMI-PostgreSQL driver error: select * from grass_gis.APGB_aerial_2_i_segment_drained_peat _with_exposed_haggs_vct where 1 = 0 ERROR: relation "grass_gis.apgb_aerial_2_i_segment_drained_ peat_with_exposed_haggs_vct" does not exist LINE 1: select * from grass_gis.APGB_aerial_2_i_segment_drained_peat... ^ DBMI-PostgreSQL driver error: select * from grass_gis.APGB_aerial_2_i_segment_drained_peat _with_exposed_haggs_vct where 1 = 0 ERROR: relation "grass_gis.apgb_aerial_2_i_segment_drained_ peat_with_exposed_haggs_vct" does not exist LINE 1: select * from grass_gis.APGB_aerial_2_i_segment_drained_peat... ^ *ERROR: Unable to describe table* Thanks, *Paul J. Shapley* ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] G_malloc errors in i.segment
Hi Marcus, Thanks again for your patience with this issue. It seems that whatever i set in the memory of 'i.segment' to it throws up the same 'G_malloc' error with the same 3 gb shortfall. I had another idea which is to make a smaller subset image (but if i do this) will the new subset of 'id' segments be applicable if applied to the entire image. For example if a value of '5134' represents 'peat' could i use 'r.mapcalc' to query the entire image for those values of '5134'. Supervised and Unsupervised Classification would not yield the differences in the landscape i'm looking for which 'i.segment' succeeds. Thanks, Paul Shapley On 21 June 2018 at 08:52, Markus Metz wrote: > [please keep it in the list] > > On Thu, Jun 21, 2018 at 9:12 AM, Paul Shapley wrote: > > > > Thank you Markus...here are the results of g.region -p and r.info > > > > g.region -p > > > projection: 99 (OSGB 1936 / British National Grid) > > zone: 0 > > datum: osgb36 > > ellipsoid: airy > > north: 242815 > > south: 200862.25 > > west: 263134.25 > > east: 298833.75 > > nsres: 0.25 > > ewres: 0.25 > > rows: 167811 > > cols: 142798 > > cells: 23963075178 > > (Thu Jun 21 08:03:51 2018) Command finished (0 sec) > > > > (Thu Jun 21 08:08:02 2018) > > > r.info map=APGB_aerial_1.1@PERMANENT > > > +-- > --+ > > | Map: APGB_aerial_1.1@PERMANENT Date: Tue May 29 09:11:31 > 2018| > > | Mapset: PERMANENT Login of Creator: > administrator | > > | Location: TempLocation > | > > | DataBase: C:\ > | > > | Title:APGB_aerial_1.1 > | > > | Timestamp: none > | > > |-- > --| > > | > | > > | Type of Map: raster Number of Categories: 0 > | > > | Data Type:CELL > | > > | Rows: 167811 > | > > | Columns: 142798 > | > > | Total Cells: 23963075178 > | > > |Projection: OSGB 1936 / British National Grid > | > > |N: 242815S: 200862.25 Res: 0.25 > | > > |E: 298833.75W: 263134.25 Res: 0.25 > | > > | Range of data:min = 0 max = 255 > | > > | > | > > | Data Description: > | > > |generated by r.in.gdal > | > > | > | > > | Comments: > | > > |r.in.gdal -o -k input="J:\Spatial Data\Aerial\APGB_received_20180416\ > | > > |\source\RGB_25cm\output\test\JPEG.png" output="APGB_aerial_1" > memory\ | > > |=300 offset=0 num_digits=0 > | > > | > | > > +-- > --+ > > (Thu Jun 21 08:08:02 2018) Command finished (0 sec) > > OK, so you really have a group of rasters with 24 billion cells. Still, do > you really need to run i.segment on the full extents with such a high > resolution? If yes, try i.segment with memory=25000. Even if the system has > more RAM, the operating system and any other programs also need some RAM, > therefore it is safer to not use all free memory. Be aware that i.segment > will take a long time (several days) to process such a large region. > > > > > Should the output image from 'r.in.gdal' have had more memory allocated > than the 300mb? > > The memory option of r.in.gdal applies to the input image, not the output > image. Sometimes more memory can help, sometimes it does not have any > effect. However, it does not affect subsequent processing steps. > > Markus M > > > > > Thanks again Markus...I appreciate your help. > > > > Paul Shapley > > > > > > > > > > On 20 June 2018 at 16:15, Markus Metz > wrote: > >> > >> > >> > >> On Wed, Jun 20, 2018 at 1:30 PM, Paul Shapley > wrote: > >> > > >> > Hi Users, > >> > > >> > Apologies if i'm posting this error again but i need to get pass the > memory allocation issue of trying to allocate 3gb memory. The Grass version > is (using file r.contour). Processor has 32 gb ram and i've allocated 3 > mb under memory option of 'i.segment'. > >> > > >> > System Info > > >> > GRASS version: 7.4.1 > > >> > GRASS SVN revision: r72807 > > >> > Build date: 2
[GRASS-user] G_malloc errors in i.segment
Hi Users, Apologies if i'm posting this error again but i need to get pass the memory allocation issue of trying to allocate 3gb memory. The Grass version is (using file r.contour). Processor has 32 gb ram and i've allocated 3 mb under memory option of 'i.segment'. System Info GRASS version: 7.4.1 GRASS SVN revision: r72807 Build date: 2018-06-13 Build platform: i386-w64-mingw32 GDAL: 2.2.4 PROJ.4: 4.9.3 GEOS: 3.5.0 SQLite: 3.17.0 Python: 2.7.4 wxPython: 2.8.12.1 Platform: Windows-2008ServerR2-6.1.7601-SP1 (OSGeo4W) ERROR: G_malloc: unable to allocate 2995426350 bytes of memory at imagery/i.segment/flag.c:25 -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] i.segment
Hi Users, I am getting an error when i run i.segment on a small subset of an image:- 'ERROR: Insufficient number of non-NULL cells in current region' Googled the error and found this (if it helps) https://gis.stackexchange.com/questions/256431/error-in-grass-using-i-segment-uspo-function Thanks, -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] i.gensig - multiple signatures into a single file for i.maxlik
Hi Users, What is the correct format for placing multiple signatures into a single ascii file. Is a '#' required between each signature along with a blank line? or can i start the new signature with line 2 e.g. #grass and is there a '#' after the last line to end the file. 1 # 2 #water 3 4186 4 67.9508 48.7346 37.8915 15.3129 13.8473 12.0855 5 1.74334 6 0.439504 2.07267 7 0.662523 1.63501 4.21189 8 0.530339 2.40757 5.52857 22.433 9 0.561184 2.30762 5.18846 20.5364 20.4926 10 0.393218 1.2184 2.63628 9.61528 9.36025 5.8531411 11 #grass 12 2132 13 etc *OR* 1 # 2 #water 3 4186 4 67.9508 48.7346 37.8915 15.3129 13.8473 12.0855 5 1.74334 6 0.439504 2.07267 7 0.662523 1.63501 4.21189 8 0.530339 2.40757 5.52857 22.433 9 0.561184 2.30762 5.18846 20.5364 20.4926 10 0.393218 1.2184 2.63628 9.61528 9.36025 5.8531411 11 # 12 #grass 13 etc Many thanks, -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] i.maxlik yellow map
Hi Users, I am re-posting this to see if anyone else may have some suggestions. My thanks to Veronica who suggested checking 'g.region' but this appears to be fine. I've also tried the same on Grass 7.2.2/7.4 and nightly but all produce the same result. signature file content:- # #Bare_peat 2667621 60.4424 75.4845 87.2587 716.802 831.829 1086.7 869.829 1105.97 1179.92 0 0 0 -- Original message:- I have created a signature file in the g.gui.iclass (supervised classification module) and it looks correct in 'preview' and when i inspected the 'sig' file. When the signature is applied using 'i.maxlik' the result is a solid yellow colour as a layer. What could have gone wrong? -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Fwd: Supervised Classification Tool - g.gui.iclass
Hi Users, I want to know if there is any good sources of information on the above wxlClass. I want to know roughly how large the training areas should be to prevent the module from taking several hours to process a signature. Also in the Training Area Display must i use a composite RGB image to digitise a training area from or will a single channel be okay? Many thanks, -- *Paul J. Shapley *MSc CGeog (GIS) FRGS -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] i.maxlik yellow map
Hi Users, I have created a signature file in the g.gui.iclass (supervised classification module) and it looks correct in 'preview' and when i inspected the 'sig' file. When the signature is applied using 'i.maxlik' the result is a solid yellow colour as a layer. What could have gone wrong? Thanks, -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Supervised Classification Tool - g.gui.iclass
Hi Users, I want to know if there is any good sources of information on the above wxlClass. I want to know roughly how large the training areas should be to prevent the module from taking several hours to process a signature. Also in the Training Area Display must i use a composite RGB image to digitise a training area from or will a single channel be okay? Many thanks, -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.import/r.in.gdal jpeg2000 image
Hi Users, Trying to import a jpeg2000 image (size 8.78 gb) with both r.import and r.in.gdal but get the error below for each module. I am using the 64 bit version of Grass 7.4.0 on Windows Server 2008 Current region rows: 111307, cols: 111307 "ERROR: G_malloc: unable to allocate 1142385 bytes of memory at lib/raster/put_row.c:343" Any ideas how i can import it? -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] 3D Grid
Hi, I want to permanently display the grey x/y/z (3D grid) in a 3D view (it actually appears briefly when you apply any movement tool) in the 3D view. I can use a constant surface and a transparency value but it's not as useful as a grid you can adjust the resolution of. Is there a way to do it in either 7.3svn or 7.0.5? Thanks, -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.in.ogr - Grass 7.0.4RC1
Hi Users, Anyone else having issues getting 'v.in.ogr' to bring in Postgres/PostGIS layers. This happens when you select 'Database' then 'Format: PostgreSQL' then 'Name:your_database_name' *no* list of vector layers appear. Have used the new 'db.login' settings for 'Hostname' and 'Port' so at least i'm not getting the old error about incorrect password. -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Grass7 GUI does not start
Hi Laurent, Have you tried: g.gui wxpython Paul Shapley Sent from my iPad > On 25 Mar 2016, at 22:28, Laurent C. <lrn...@gmail.com> wrote: > > Dimitris, > > How did you install GRASS ? From the Suse repositories? > What happen if you try to launch the GUI manually with g.gui (possibly > with --v flag)? > > Regards, > Laurent > > 2016-03-18 7:42 GMT-06:00 ddarlac <dimn...@rutgers.edu>: >> I am a Qgis user, and wanted to do some 3D videos - I read Nviz might be the >> way; I wanted to work on it through Grass (not through the Qgis plugin), and >> grass fails to start. >> the start up window starts fine, I created the location and the mapset - >> when I hit >> 'Start GRASS session', I am getting: >> /Launching GUI in the background, please wait... >> GRASS 7.0.3 (Manhattan):~ > / >> and nothing happens after that. >> >> My system is opensuse leap, and the grass installation is 7.03; >> wxpython is installed >> >> any help is appreciated >> Dimitris >> >> >> >> -- >> View this message in context: >> http://osgeo-org.1560.x6.nabble.com/Grass7-GUI-does-not-start-tp5257181.html >> Sent from the Grass - Users mailing list archive at Nabble.com. >> ___ >> grass-user mailing list >> grass-user@lists.osgeo.org >> http://lists.osgeo.org/mailman/listinfo/grass-user > ___ > grass-user mailing list > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] GRASS Revision 64733
Hi users, I'm trying to get a working QGIS/GRASS plugin on Ubuntu 15.10 but with little luck. I seem to need the r64733 in order to be compatible with QGIS 2.14. My question is what is the easiest way of installing this revision of GRASS. Thanks, -- *Paul J. Shapley* ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] file formats from Ordnance Survey Survey Maps
I'm attempting to get some maps and use them in GRASS, which gives me a choice of formats to download from. Can anyone tell me if any one of these formats is better and will allow direct import into GRASS 7.0.2. without transformation. For one type of topography map I can get either GML2, DWG or File Geodatabase. For another type of Terrain type map I can get either Shape, GML3 or DWG. I have googled and looked at the choices provided by Hi Jim, We use either shapefile (not ideal for using annotation) or import directly into a Postgres/PostGIS database first. -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] band combinations using GRASS
Hi Josep, The images probably need a 'Histogram Stretch' using 'r.colors' use the 'grey' (color) and tick Histogram Stretch' but make sure you set 'g.region' to the full image first otherwise it wont work. Paul Shapley On 10 December 2015 at 11:03, Josep Sitjar <jo...@sigte.udg.edu> wrote: > Dear all, > I'm working with r.composite to create a natural color band combination > using Landsat 8 image. The result is a very dark image, so I would like to > know if there's any way to improve it. I'm using GRASS 6.X and also GRASS > 7. > > Thanks in advance, > > -- > *Josep Sitjar Suñer* > Tècnic SIG / Servei de Sistemes d'Informació Geogràfica i Teledetecció > - > Universitat de Girona > *SIGTE* > - > Pl. Ferrater Mora 1 > 17071 Girona > Tel +34 972 418 039 (7026 intern) > jo...@sigte.udg.edu > > http://www.sigte.udg.edu > Twitter http://twitter.com/SIGTE_UDG > > > ___ > grass-user mailing list > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user > -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] band combinations using GRASS
Hi Josep, My apologies but i should have read your email thoroughly. Please check the link below hope this helps you. http://www.paulshapley.com/2014/04/processing-landsat-8-bands-in-grass-gis.html Regards, Paul Shapley On 10 December 2015 at 11:03, Josep Sitjar <jo...@sigte.udg.edu> wrote: > Dear all, > I'm working with r.composite to create a natural color band combination > using Landsat 8 image. The result is a very dark image, so I would like to > know if there's any way to improve it. I'm using GRASS 6.X and also GRASS > 7. > > Thanks in advance, > > -- > *Josep Sitjar Suñer* > Tècnic SIG / Servei de Sistemes d'Informació Geogràfica i Teledetecció > - > Universitat de Girona > *SIGTE* > - > Pl. Ferrater Mora 1 > 17071 Girona > Tel +34 972 418 039 (7026 intern) > jo...@sigte.udg.edu > > http://www.sigte.udg.edu > Twitter http://twitter.com/SIGTE_UDG > > > ___ > grass-user mailing list > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user > -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Link external vector data - grass 7.0.2RC2
Hi Users, Getting issues trying to connect to a postgres database on Win7. -- "Error in db.databases unable to list databases. Try to define correct connection settings by db.login" -- The connection settings are correct and work with other software except Grass. If i use 'db.tables -p' the tables are all visible. Are the db.databases and db.login modules still got issues?. Any suggestions welcome. Thanks -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] db.test
Hi Grass Users, Trying to connect to an existing postgres/postgis database in 7.0.2RC1 but getting older connection details reported which may be preventing a newer connection. Here is the output from 'db.test' below. Note the 'VAR' file is showing the newest details (not revealed in full for the list). References in the 'search_path' to 'public' and 'schema sde' details are from connections going back some three years and creating a new mapset does not change this:- - WARNING: invalid value for parameter "search_path": ""$user", public, sde" DETAIL: schema "sde" does not exist - (Mon Oct 26 12:39:56 2015) db.test test=test1 Using DB driver: pg WARNING: invalid value for parameter "search_path": ""$user", public, sde" DETAIL: schema "sde" does not exist EXECUTE: OK WARNING: invalid value for parameter "search_path": ""$user", public, sde" DETAIL: schema "sde" does not exist EXECUTE: OK WARNING: invalid value for parameter "search_path": ""$user", public, sde" DETAIL: schema "sde" does not exist EXECUTE: OK DBMI-PostgreSQL driver error: Unable to select: SELECT * FROM grass_test1 ERROR: relation "grass_test1" does not exist LINE 1: SELECT * FROM grass_test1 ^ WARNING: invalid value for parameter "search_path": ""$user", public, sde" DETAIL: schema "sde" does not exist DBMI-PostgreSQL driver error: Unable to select: SELECT * FROM grass_test1 ERROR: relation "grass_test1" does not exist LINE 1: SELECT * FROM grass_test1 ^ ERROR: EXECUTE: ERROR 0a1,2 > 1|123.456|abcd > 2||xxx ERROR: RESULT: ERROR DBMI-PostgreSQL driver error: Unable to select: SELECT c1 FROM grass_test1 WHERE d1 < 500 / 2 AND i1 <> 2 AND c1 LIKE '%bc%' ERROR: relation "grass_test1" does not exist LINE 1: SELECT c1 FROM grass_test1 WHERE d1 < 500 / 2 AND i1 <> 2 A... ^ WARNING: invalid value for parameter "search_path": ""$user", public, sde" DETAIL: schema "sde" does not exist DBMI-PostgreSQL driver error: Unable to select: SELECT c1 FROM grass_test1 WHERE d1 < 500 / 2 AND i1 <> 2 AND c1 LIKE '%bc%' ERROR: relation "grass_test1" does not exist LINE 1: SELECT c1 FROM grass_test1 WHERE d1 < 500 / 2 AND i1 <> 2 A... ^ ERROR: EXECUTE: ERROR 0a1 > abcd ERROR: RESULT: ERROR WARNING: invalid value for parameter "search_path": ""$user", public, sde" DETAIL: schema "sde" does not exist EXECUTE: OK WARNING: invalid value for parameter "search_path": ""$user", public, sde" DETAIL: schema "sde" does not exist EXECUTE: OK WARNING: invalid value for parameter "search_path": ""$user", public, sde" DETAIL: schema "sde" does not exist EXECUTE: OK DBMI-PostgreSQL driver error: Unable to select: SELECT * FROM grass_test1 ERROR: relation "grass_test1" does not exist LINE 1: SELECT * FROM grass_test1 ^ WARNING: invalid value for parameter "search_path": ""$user", public, sde" DETAIL: schema "sde" does not exist DBMI-PostgreSQL driver error: Unable to select: SELECT * FROM grass_test1 ERROR: relation "grass_test1" does not exist LINE 1: SELECT * FROM grass_test1 ^ 0a1,3 > 1|123.456|abcd| > 2|18.6|xxx|987 > 3|0|_\'_| ERROR: EXECUTE: ERROR ERROR: RESULT: ERROR CREATE TABLE grass_test1 (i1 INTEGER, d1 DOUBLE PRECISION, c1 VARCHAR(20)) INSERT INTO grass_test1 VALUES ( 1, 123.456, 'abcd' ) INSERT INTO grass_test1 VALUES ( 2, null, 'xxx' ) SELECT * FROM grass_test1 SELECT c1 FROM grass_test1 WHERE d1 < 500 / 2 AND i1 <> 2 AND c1 LIKE '%bc%' INSERT INTO grass_test1 VALUES ( 3, 0.0, '_\''_' ) ALTER TABLE grass_test1 ADD COLUMN i2 INTEGER UPDATE grass_test1 SET d1 = 18.6, i2 = 987 WHERE i1 = 2 SELECT * FROM grass_test1 DROP TABLE grass_test1 WARNING: invalid value for parameter "search_path": ""$user", public, sde" DETAIL: schema "sde" does not exist EXECUTE: OK (Mon Oct 26 12:39:58 2015) Command finished (2 sec) -- 'VAR' File contents below DB_DRIVER: pg DB_DATABASE: host=?.?-??.???.??,dbname=??? DB_SCHEMA: grass_gis DB_GROUP: postgres - -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.colors
Hi All, Just testing 7.0.2RC1 with some landsat 8 images. I've imported and changed the DN values to TOAR values but when i try to use a 'Histogram Stretch' using the 'grey' table the maps display no color (white display) whereas in version 6 it worked. Also i need to increase the 'contrast' when displaying the individual bands but dont know if there is a way to do this. I have tried opening the 'i.fusion.hpf' script but it seems to hang. Many Thanks, -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Migrating to a different database engines in Grass 7.0.1RC2
Hello Grass Users, I would like to know why and where grass stores old database connection parameters set with 'db.connect' and 'db.login'. this is not in the 'VAR' file in the PERMANENT directory, I often switch between PostGIS and SQLite data but even if i add a new connection it seems to retrieve a three year old connection parameter and i dont want to keep having to create a new 'mapset' every time i change. Should i go back to 6.4 stable or is there another solution? Many thanks, -- *Paul J. Shapley *MSc CGeog (GIS) FRGS ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Projection errors
Hi Users, I'm having the problem of projecting a source map (landsat 8 image) into a target location (Ordnance Survey UK) in Grass 7. I have tried the r.pack/unpack method but i get the following error:- 'WARNING: Projection information does not match. Proceeding...' The resulting view shows both images but the landsat image is many kilometres to the north of where it should overlap. Have also made 'PROJ2'(see below) with the same n,s,e,w as 'PROJ1'(see below) but still the same result as above. As i understand it the underlying Landsat images are in the WGS84/lon-lat format. I'm wondering should i use m.proj to first convert the co-ordinates to OrdSurvey? I seem to recall just using 'r.proj' in the past with no problems at all but now i get error:- 'WARNING: pj_transform() failed: latitude or longitude exceeded limits' I'm a little confused as to which module to use to get around this problem. g.region -b north latitude: 89:59:47.524206N south latitude: 49:45:56.915978N west longitude: 18:03:37.69897W east longitude: 0:02:23.60354W center longitude: 9:03:00.651255W center latitude: 69:52:52.220092N 'PROJ1' - DEFAULT_WIND proj: 99 zone: 0 north: 5848507.5 south: 5605492.5 east: 540907.5 west: 300892.5 cols: 540915 rows: 5848815 e-w resol: 1 n-s resol: 1 top:1.000 bottom: 0.000 cols3: 1 rows3: 1 depths: 1 e-w resol3: 1 n-s resol3: 1 t-b resol: 1 'PROJ2' - DEFAULT_WIND proj: 99 zone: 0 north: 123 south: -123 east: 66 west: -66 cols: 108220 rows: 1169781 e-w resol: 5 n-s resol: 5 top:1.000 bottom: 0.000 cols3: 1 rows3: 1 depths: 1 e-w resol3: 5 n-s resol3: 5 t-b resol: 1 Thanks, -- *Paul J. Shapley* ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.proj (grass 7)
Hi Users, Having some problems moving some landsat 8 composite images into a different location/mapset with the same projection. Have tried r.proj (this produces the 'ERROR: Input raster map is outside current region') and v.in.region method (this produces a 'blank' map). Not sure if the resolutions need changing with g.region? Anyone else had this? Source Location g.region -p projection: 99 (Transverse Mercator) zone: 0 datum: osgb36 ellipsoid: airy north: 5848507.5 south: 5605492.5 west: 300892.5 east: 540907.5 nsres: 15 ewres: 15 rows: 16201 cols: 16001 cells: 259232201 Target Location g.region -p projection: 99 (Transverse Mercator) zone: 0 datum: osgb36 ellipsoid: airy north: 1466741.93548387 south: -315435.48387097 west: -309963.70967742 east: 1037084.67741935 nsres: 14.7828 ewres: 15.3772 rows: 118812 cols: 89803 cells: 10669674036 *Paul J. Shapley * ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Landsat 8 Workflow for Grass 7.0.0RC2
Hi Grass Users, I want to use the following Landsat 8 workflow to obtain some pansharpened images but have a problem when i convert the DN values to 8bit for pansharpening. Here is my preferred workflow in Grass 7.0.0RC2. · r.in.gdal (success) · i.landsat.toar (success) · r.rescale bands 432 8 (convert reflectance 'DN' values to 0-255, 8bit) · i.pansharpen bands 432 with 8 · r.rescale (convert 0-255 values back to reflectance 'DN) · i.colors.enhance (pansharp bands) · d.rgb My question...is this a valid workflow? because i'm getting the same two errors below and the images when displayed look solid grey when viewed with d.rgb:- *displaying the pansharpened images error log:-* (Tue Feb 17 14:57:26 2015) g.list rast LC82040242014206LGN00_B1 LC82040242014206LGN00_B10 LC82040242014206LGN00_B11 LC82040242014206LGN00_B2 LC82040242014206LGN00_B3 LC82040242014206LGN00_B4 LC82040242014206LGN00_B5 LC82040242014206LGN00_B6 LC82040242014206LGN00_B7 LC82040242014206LGN00_B8 LC82040242014206LGN00_B9 LC82040242014206LGN00_BQA LC82040242014206LGN00_refl1 LC82040242014206LGN00_refl10 LC82040242014206LGN00_refl11 LC82040242014206LGN00_refl2 LC82040242014206LGN00_refl3 LC82040242014206LGN00_refl4 LC82040242014206LGN00_refl5 LC82040242014206LGN00_refl6 LC82040242014206LGN00_refl7 LC82040242014206LGN00_refl8 LC82040242014206LGN00_refl9 LC82040242014318LGN00_B1 LC82040242014318LGN00_B10 LC82040242014318LGN00_B11 LC82040242014318LGN00_B2 LC82040242014318LGN00_B3 LC82040242014318LGN00_B4 LC82040242014318LGN00_B5 LC82040242014318LGN00_B6 LC82040242014318LGN00_B7 LC82040242014318LGN00_B8 LC82040242014318LGN00_B9 LC82040242014318LGN00_BQA LC82040242014318LGN00_refl1 LC82040242014318LGN00_refl10 LC82040242014318LGN00_refl11 LC82040242014318LGN00_refl2 LC82040242014318LGN00_refl2_rescale LC82040242014318LGN00_refl3 LC82040242014318LGN00_refl3_rescale LC82040242014318LGN00_refl4 LC82040242014318LGN00_refl4_rescale LC82040242014318LGN00_refl5 LC82040242014318LGN00_refl6 LC82040242014318LGN00_refl7 LC82040242014318LGN00_refl8 LC82040242014318LGN00_refl8_rescale LC82040242014318LGN00_refl9 (Tue Feb 17 14:57:27 2015) Command finished (0 sec) Command 'd.rast map=LC82040242014318LGN00_pansharp@PERMANENT' failed Details: Raster map LC82040242014318LGN00_pansharp@PERMANENT not found *i.colors.enhance error log:-* (Tue Feb 17 13:06:40 2015) i.colors.enhance red=LC82040242014318LGN00_B4@PERMANENT green=LC82040242014318LGN00_B3@PERMANENT blue=LC82040242014318LGN00_refl2@PERMANENT strength=96 Processing... Process Process-3: Traceback (most recent call last): File C:\OSGeo4W\apps\Python27\lib\multiprocessing\process.py, line 258, in _bootstrap self.run() File C:\OSGeo4W\apps\Python27\lib\multiprocessing\process.py, line 114, in run self._target(*self._args, **self._kwargs) File C:\OSGeo4W\apps\grass\grass-7.0.0RC2\scripts\i.color s.enhance.py, line 103, in get_percentile_mp result = get_percentile(map, percentiles) File C:\OSGeo4W\apps\grass\grass-7.0.0RC2\scripts\i.color s.enhance.py, line 89, in get_percentile percentiles = values, quiet = True) File C:\OSGeo4W\apps\grass\grass-7.0.0RC2\etc\python\gras s\script\core.py, line 427, in read_command return handle_errors(returncode, stdout, args, kwargs) File C:\OSGeo4W\apps\grass\grass-7.0.0RC2\etc\python\gras s\script\core.py, line 310, in handle_errors returncode=returncode) CalledModuleError: Module run None ['r.quantile', '--q', 'input=LC82040242014318LGN00_refl2@PERMANENT', 'percentiles=2,96'] ended with error Process ended with non-zero return code -1073741571. See errors in the (error) output. Process Process-2: Traceback (most recent call last): File C:\OSGeo4W\apps\Python27\lib\multiprocessing\process.py, line 258, in _bootstrap self.run() File C:\OSGeo4W\apps\Python27\lib\multiprocessing\process.py, line 114, in run self._target(*self._args, **self._kwargs) File C:\OSGeo4W\apps\grass\grass-7.0.0RC2\scripts\i.color s.enhance.py, line 103, in get_percentile_mp result = get_percentile(map, percentiles) File C:\OSGeo4W\apps\grass\grass-7.0.0RC2\scripts\i.color s.enhance.py, line 89, in get_percentile percentiles = values, quiet = True) File C:\OSGeo4W\apps\grass\grass-7.0.0RC2\etc\python\gras s\script\core.py, line 427, in read_command return handle_errors(returncode, stdout, args, kwargs) File C:\OSGeo4W\apps\grass\grass-7.0.0RC2\etc\python\gras s\script\core.py, line 310, in handle_errors returncode=returncode) CalledModuleError: Module run None ['r.quantile', '--q', 'input=LC82040242014318LGN00_B3@PERMANENT', 'percentiles=2,96'] ended with error Process ended with non-zero return code -1073741571. See errors in the (error) output. Process
[GRASS-user] Reconnecting vector attribution with Postgres/PostGIS
Hi Grass 7 users. I am trying to connect some vector layers with thier respective attribute tables in postgres/postgis. The tables all connected under version 6.4.4 but not since upgrading to version 7RC1 (on winxp in this instance) they do not. It looks like Grass cannot find the database password but it is present in 'db.login'. Has anyone else had this happen? Here are the outputs:- *Grass Command Prompt* Cleaning up temporary files... Starting GRASS GIS... ERROR: Location demolocation not found WARNING: Concurrent mapset locking is not supported on Windows __ ___ _____ / / __ \/ | / ___/ ___/ / / _/ ___/ / / __/ /_/ / /| | \__ \\_ \ / / __ / / \__ \ / /_/ / _, _/ ___ |___/ /__/ / / /_/ // / ___/ / \/_/ |_/_/ |_/// \/___/// Welcome to GRASS 7.0.0RC1 GRASS homepage: http://grass.osgeo.org This version running through:Command Shell (C:\WINDOWS\system32\cmd. exe) Help is available with the command: g.manual -i See the licence terms with: g.version -c If required, restart the GUI with: g.gui wxpython When ready to quit enter:exit Launching wxpython GUI in the background, please wait... Microsoft Windows XP [Version 5.1.2600] (C) Copyright 1985-2001 Microsoft Corp. H:\DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis GRASS_INFO_WARNING(2500,1): Unable to open database host=gisdb.neath-porttalbot .gov.uk, dbname=postgis by driver pg GRASS_INFO_END(2500,1) GRASS_INFO_WARNING(2500,2): Unable open database host=gisdb.neath-porttalbot.go v.uk, dbname=postgis by driver pg GRASS_INFO_END(2500,2) DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis GRASS_INFO_ERROR(2500,3): Unable to open database host=gisdb.neath-porttalbot.g ov.uk, dbname=postgis GRASS_INFO_END(2500,3) DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis GRASS_INFO_WARNING(4928,1): Unable to open database host=gisdb.neath-porttalbot .gov.uk, dbname=postgis by driver pg GRASS_INFO_END(4928,1) GRASS_INFO_WARNING(4928,2): Unable open database host=gisdb.neath-porttalbot.go v.uk, dbname=postgis by driver pg GRASS_INFO_END(4928,2) DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis GRASS_INFO_ERROR(4928,3): Unable to open database host=gisdb.neath-porttalbot.g ov.uk, dbname=postgis GRASS_INFO_END(4928,3) DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis GRASS_INFO_WARNING(3988,1): Unable to open database host=gisdb.neath-porttalbot .gov.uk, dbname=postgis by driver pg GRASS_INFO_END(3988,1) GRASS_INFO_WARNING(3988,2): Unable open database host=gisdb.neath-porttalbot.go v.uk, dbname=postgis by driver pg GRASS_INFO_END(3988,2) DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis GRASS_INFO_ERROR(3988,3): Unable to open database host=gisdb.neath-porttalbot.g ov.uk, dbname=postgis GRASS_INFO_END(3988,3) DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis GRASS_INFO_WARNING(4156,1): Unable to open database host=gisdb.neath-porttalbot .gov.uk, dbname=postgis by driver pg GRASS_INFO_END(4156,1) GRASS_INFO_WARNING(4156,2): Unable open database host=gisdb.neath-porttalbot.go v.uk, dbname=postgis by driver pg GRASS_INFO_END(4156,2) DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis GRASS_INFO_ERROR(4156,3): Unable to open database host=gisdb.neath-porttalbot.g ov.uk, dbname=postgis GRASS_INFO_END(4156,3) DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis DBMI-PostgreSQL driver error: Unknown option in database definition for PostgreSQL: dbname=postgis
[GRASS-user] r.proj
I am trying to re-project some landsat 8 bands in from one location into another location with the same Ordnance Survey Projection (epsg:27700) using r.proj. They should overlap but currently do not. Am i using the correct module? The error message is: -- (Tue Apr 29 11:30:48 2014) r.proj input=LC82040242013363LGN00_refl1 location=Landsat8 mapset=Landsat dbase=\\msastore\decsgroup\GIS Admin support share output=LC82040242013363LGN00_refl1 Input Projection Parameters: +proj=tmerc +lat_0=23 +lon_0=-96 +k=1 +x_0=0 +y_0=0 +no_defs +a=6377563.396 +rf=299.3249646 +towgs84=368.000,-120.000,425.000 Input Unit Factor: 1 Output Projection Parameters: +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=40 +y_0=-10 +no_defs +a=6377563.396 +rf=299.3249646 +towgs84=368.000,-120.000,425.000 Output Unit Factor: 1 WARNING: pj_transform() failed: latitude or longitude exceeded limits WARNING: pj_transform() failed: latitude or longitude exceeded limits WARNING: pj_transform() failed: latitude or longitude exceeded limits -- r.proj parameters From Projection projection: 99 (Transverse Mercator) zone: 0 datum: osgb36 ellipsoid: airy north: 5848830 south: 5605470 west: 296070 east: 536430 nsres: 30 ewres: 30 rows: 8112 cols: 8012 cells: 64993344 to Projection projection: 99 (Transverse Mercator) zone: 0 datum: osgb36 ellipsoid: airy north: 123 south: 0 west: 0 east: 66 nsres: 50 ewres: 50 rows: 24600 cols: 13200 cells: 32472 -- *Paul J. Shapley* ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Grass 7 ubuntu 13.10
Hi, I am trying to start Grass 7 svn on Ubuntu 13.10. I used the following command in terminal: sudo apt-get update; sudo apt-get install grass70 grass70-gui When i start Grass with 'grass70' in terminal the Grass 'splash' screen opens and then promptly crashes. Here is the error: username@ServerName localhost:~$ grass70 WELCOME TO GRASS 7.0.svn 1) Have at your side all available GRASS tutorials 2) When working on your location, the following materials are extremely useful: - A topo map of your area - Current catalog of available computer maps 3) Check the GRASS webpages for feedback mailinglists and more: http://grass.osgeo.org http://www.grass-gis.org Hit RETURN to continue Starting GRASS GIS... Fontconfig warning: /etc/fonts/conf.d/50-user.conf, line 14: reading configurations from ~/.fonts.conf is deprecated. __ ___ _____ / / __ \/ | / ___/ ___/ / / _/ ___/ / / __/ /_/ / /| | \__ \\_ \ / / __ / / \__ \ / /_/ / _, _/ ___ |___/ /__/ / / /_/ // / ___/ / \/_/ |_/_/ |_/// \/___/// Welcome to GRASS 7.0.svn GRASS homepage: http://grass.osgeo.org This version running through:Bash Shell (/bin/bash) Help is available with the command: g.manual -i See the licence terms with: g.version -c If required, restart the GUI with: g.gui wxpython When ready to quit enter:exit Launching wxpython GUI in the background, please wait... GRASS 7.0.svn (dtm_location):~ Fontconfig warning: /etc/fonts/conf.d/50-user.conf, line 14: reading configurations from ~/.fonts.conf is deprecated. Traceback (most recent call last): File /usr/lib/grass70/etc/gui/wxpython/wxgui.py, line 143, in module sys.exit(main()) File /usr/lib/grass70/etc/gui/wxpython/wxgui.py, line 136, in main app = GMApp(workspaceFile) File /usr/lib/grass70/etc/gui/wxpython/wxgui.py, line 50, in __init__ wx.App.__init__(self, False) File /usr/lib/python2.7/dist-packages/wx-2.8-gtk2-unicode/wx/_core.py, line 7981, in __init__ self._BootstrapApp() File /usr/lib/python2.7/dist-packages/wx-2.8-gtk2-unicode/wx/_core.py, line 7555, in _BootstrapApp return _core_.PyApp__BootstrapApp(*args, **kwargs) File /usr/lib/grass70/etc/gui/wxpython/wxgui.py, line 84, in OnInit workspace = self.workspaceFile) File /usr/lib/grass70/etc/gui/wxpython/lmgr/frame.py, line 110, in __init__ self._menuTreeBuilder = LayerManagerMenuData() File /usr/lib/grass70/etc/gui/wxpython/lmgr/menudata.py, line 40, in __init__ MenuTreeModelBuilder.__init__(self, filename, expandAddons=expandAddons) File /usr/lib/grass70/etc/gui/wxpython/core/menutree.py, line 69, in __init__ xmlTree = etree.parse(filename) File /usr/lib/python2.7/xml/etree/ElementTree.py, line 1182, in parse tree.parse(source, parser) File /usr/lib/python2.7/xml/etree/ElementTree.py, line 657, in parse self._root = parser.close() File /usr/lib/python2.7/xml/etree/ElementTree.py, line 1654, in close self._raiseerror(v) File /usr/lib/python2.7/xml/etree/ElementTree.py, line 1506, in _raiseerror raise err xml.etree.ElementTree.ParseError: no element found: line 1, column 0 -- *Paul Shapley* ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Fwd: FW: connecting to postgresql and postgis backend database - Grass 6.4.4svn
Hi, I'm having some problems with password authentication when trying to connect to postgresql. Here is the generated error: (Fri Mar 21 10:25:16 2014) g.remove vect=building@PERMANENT Removing vector building@PERMANENT DBMI-Postgres driver error: Cannot connect to Postgres: FATAL: password authentication failed for user pshapley WARNING: Unable to open database postgis by driver pg WARNING: Unable open database postgis by driver pg WARNING: Unable to find table npt_planning.building linked to vector map building WARNING: couldn't be removed WARNING: building nothing removed (Fri Mar 21 10:25:17 2014) Command finished (0 sec) Is this problem linked with Ticket #1274 because if i use 'v.in.ogr' (with the -t) option it imports all the geometries without the attributes (not much use)?. When the import is finished I am asked to 'try to rebuild topologies with v.build' but this does not work. Any ideas would be great. Thanks, *Paul Shapley* ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] grass 6.4.3 and db.connect
Hi, I am using 'db.connect' to get access to a sqlite database. I have renamed the database extension from '.sqlite' to '.db' and placed it in the 'PERMANENT' directory but not entirely sure if i've correctly set this sqlite access up correctly. When i use 'db.test' i get the following error. Unable to open database: unable to open database file ERROR: Unable to open database '/home/username/dtm_location/PERMANENT/gisdata.db' ERROR: EXECUTE: ERROR I'm sure i'm probably making some simple errors in the configuration somewhere. Can anyone help with setup? Many thanks, -- *Paul J. Shapley * ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.hazard.flood.py - GRASS 6.4.3RC2 WinXP
Hi, This flood module is working in the latest stable release, however there is an issue in that it does not produce the flood 'test1' map only 'test2' (the MTI map). Running the module seems to generate three maps 'r.clump', 'r.flood' and 'test2'. Just need confirmation that this is how it is meant to be and is the error below normal or does the module require fine tuning. Also is there any detailed information on this module? Where can i find it?. r.hazard.flood.py map=mosaic@PERMANENT flood=test1 mti=test2 Cellsize : 5.0 SECTION 1a (of 4): Initiating Memory. SECTION 1b (of 4): Determining Offmap Flow. SECTION 2: A * Search. SECTION 3: Accumulating Surface Flow with MFD. SECTION 4: Closing Maps. WARNING: Writing out only positive flow accumulation values. WARNING: Cells with a likely underestimate for flow accumulation can no longer be identified. Flow accumulation done. Slope raster map r_slope complete Slope map done. Exponent : 0.0335464203571 MTI threshold : 2.64732051769 Calculating mti raster map.. Cleaning up.. Removing raster r_accumulation Removing raster r_slope Calculating flood raster map.. Running r.clump.. Pass 1... Pass 2... r.clump complete. 133745 clumps. Deleting areas of less than 4 cells.. 'r.area' is not recognized as an internal or external command, operable program or batch file. Invalid map r_flood_th Invalid map r_flood_th Parse error ERROR: An error occurred while running r.mapcalc (Thu Jan 17 11:54:55 2013) Command finished (7 min 8 sec) -- Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.hazard.flood.py Grass 6.4.3RC1
Hi All, I'm getting the following error trying to execute the r.hazard.flood python script. The module will load and accept the parameters but does not seem to execute. Any advice appreciated. Thanks, Paul Shapley - **kwargs) wx._core . PyAssertionError : C++ assertion m_pendingEvents failed at ..\..\src\common\event.cpp(1172) in wxEvtHandler::ProcessPendingEvents(): Please call wxApp::ProcessPendingEvents() instead Traceback (most recent call last): File C:\Program Files\GRASS GIS 6.4.3RC1\etc\wxpython\gui_core\goutput.py, line 980, in OnProcessPendingOutputWindowEvents self.ProcessPendingEvents() File C:\Program Files\GRASS GIS 6.4.3RC1\Python27\lib \site-packages\wx-2.8-msw-unicode\wx\_core.py, line 3871, in ProcessPendingEvents return _core_.EvtHandler_ProcessPendingEvents(*args, **kwargs) wx._core . PyAssertionError : C++ assertion m_pendingEvents failed at ..\..\src\common\event.cpp(1172) in wxEvtHandler::ProcessPendingEvents(): Please call wxApp::ProcessPendingEvents() instead -- Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Python Scripts
Hi All, I'm trying to run a python script 'r.hazard.flood.py' in Grass 6.4.3RC1 (WinXP). I'm not sure where the script should be placed within the grass installation so i can use the GUI 'File Launch Script'. Should it be placed here: C:\Program Files\GRASS GIS 6.4.3RC1\Python27\Scripts or C:\Program Files\GRASS GIS 6.4.3RC1\etc\gui\scripts or somewhere else. If it is run from the 'gui\scripts' location it just crashes. Regards, -- Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Graphical Modeler Output
Hi Users, I'm experimenting with the Graphical Modeler in Grass 6.4.3 but seem to be having some trouble getting the output to select a monitor and display the data. Here is the output from a raster to vector conversion. r.to.vect -s --overwrite input=grass_15degrees_or_below@PERMANENToutput=grass_15degrees_or_below_modeltest@PERMANENTfeature=area WARNING: Vector map grass_15degrees_or_below_modeltest already exists and will be overwritten Extracting areas... Building topology for vector map grass_15degrees_or_below_modeltest... Registering primitives... 31632 primitives registered 435724 vertices registered Building areas... 18382 areas built 9878 isles built Attaching islands... Attaching centroids... Number of nodes: 23128 Number of primitives: 31632 Number of points: 0 Number of lines: 0 Number of boundaries: 24422 Number of centroids: 7210 Number of areas: 18382 Number of isles: 9878 Number of areas without centroid: 11172 r.to.vect complete. (Tue Nov 06 14:01:27 2012) Command finished (16 sec) (Tue Nov 06 14:01:27 2012) d.monsize setmonitor=x2 setwidth=601 setheight=391 Starting monitor x2 at 601X391 ... ERROR: no such monitor 'x2' ERROR: No such monitor as x2 ERROR: No such monitor as x2 Problem selecting x2. Will try once more Many thanks, Paul Shapley -- Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] G_malloc error
Hi Users, I'm getting 'G_malloc' errors when using large images in Grass 6.4.3rc1. I need to use much larger images than the one causing the error. Is there a way around this problem? i have 3gb of RAM on this PC (standard MS Windows installer not osgeo). Is this enough RAM or should i invest in much more... or is there another cause of this?. Error message: Current region rows: 19998, cols: 15999 ERROR: G_malloc: unable to allocate 1279792008 bytes of memory at gsds.c:575 -- Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Memory allocation in Grass 6.4.2 Windows
HI I'm getting the following error when using the r.to.vect module: ERROR: G_realloc: unable to allocate 2544 bytes at areas.c:678 I understand that this is due to the standard memory installation of 64mb of RAM allocation. How can i overcome this limitation.I also need to solve this on my Ubuntu 12.04 install -- Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Grass 6.4.2 on Ubuntu 11.10
Hi, I'm having some trouble installing the latest stable Grass 6.4.2 on Ubuntu. I've un-installed 6.4.1 but if i use the 'apt-get install Grass' command it re-installs 6.4.1. I think i may be missing something. -- Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Handling 700 Dtm files in Grass
Hi, I need to merge 700 arcinfo ascii files so they can be imported into Grass then saved as a single geotif, How can this best be done?. It is possible to import them one at a time but this is not practical. -- Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] NVIZ Wire Mesh
Hello, I want to have the NVIZ wire mesh opaque but don't know how best to do this. I don't want to be transparent. Does anyone know how this can be done? Thanks, -- Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.in.wms failure
Hi, I'm getting the following error telling me gdalwarp has failed. I'm using the windows version of Grass 6.4.2RC1. Any help appreciated as usual. Paul Shapley (Thu Nov 10 11:49:41 2011) r.in.wms output=grass_wms_request mapserver=mapserver= http://120.100.200.96:8080/geowebcache/service/wms?tiled=truelayers=os_imagery_2010 format=jpeg maxcols=1024 maxrows=1024 folder=folder=\\innogistic\mapping\test method=nearest Calculating tiles Requesting 1 tiles. mkdir: cannot create directory `folder=\\innogistic\\mapping\\test': No such file or directory C:/Program Files/GRASS 6.4RC1/etc/r.in.wms/wms.request: line 221: folder=\innogistic\mapping\test/grass_wms_request_.wget: No such file or directory C:/Program Files/GRASS 6.4RC1/etc/r.in.wms/wms.request: line 223: folder=\innogistic\mapping\test/grass_wms_request_.proj4: No such file or directory C:/Program Files/GRASS 6.4RC1/etc/r.in.wms/wms.request: line 240: folder=\innogistic\mapping\test/grass_wms_request__0.jgw: No such file or directory C:/Program Files/GRASS 6.4RC1/etc/r.in.wms/wms.request: line 241: folder=\innogistic\mapping\test/grass_wms_request__0.jgw: No such file or directory C:/Program Files/GRASS 6.4RC1/etc/r.in.wms/wms.request: line 242: folder=\innogistic\mapping\test/grass_wms_request__0.jgw: No such file or directory C:/Program Files/GRASS 6.4RC1/etc/r.in.wms/wms.request: line 243: folder=\innogistic\mapping\test/grass_wms_request__0.jgw: No such file or directory C:/Program Files/GRASS 6.4RC1/etc/r.in.wms/wms.request: line 244: folder=\innogistic\mapping\test/grass_wms_request__0.jgw: No such file or directory C:/Program Files/GRASS 6.4RC1/etc/r.in.wms/wms.request: line 245: folder=\innogistic\mapping\test/grass_wms_request__0.jgw: No such file or directory C:/Program Files/GRASS 6.4RC1/etc/r.in.wms/wms.request: line 248: folder=\innogistic\mapping\test/grass_wms_request_.wget: No such file or directory Downloading tiles cat: folder=\innogistic\mapping\test/grass_wms_request_.wget: No such file or directory All tiles downloaded successfully cat: folder=\\innogistic\mapping\test/grass_wms_request_.wget: No such file or directory file: could not find any magic files! ERROR: Required parameter input not set: (Raster file or files to be imported. If multiple files are specified they will be patched together.) Description: Warps and imports GDAL supported raster file complete with correct NULL values Keywords: raster, rotate, reproject Usage: r.in.gdalwarp [-peck] input=string[,string,...] output=string s_srs=string method=string [warpoptions=string] [--verbose] [--quiet] Flags: -p Don't reproject the data, just patch it. -e Extend location extents based on new dataset -c Make color composite image if possible -k Keep band numbers instead of using band color names --v Verbose module output --q Quiet module output Parameters: input Raster file or files to be imported. If multiple files are specified they will be patched together. output Name for resultant raster map. Each band will be name output.bandname s_srs Source projection in gdalwarp format method Reprojection method to use options: nearest,bilinear,cubic,cubicspline default: nearest warpoptions Additional options for gdalwarp ERROR: r.in.gdalwarp failed (Thu Nov 10 11:50:10 2011) Command finished (29 sec) -- Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] NVIZ Raster resolution
I'm using GRASS 6.4.2RC1-1 but cannot get a high resolution from an input dem/raster surface that i've been using with previous versions of Grass. NVIZ will display the raster okay but will not allow a high res render. I understand NVIZ is in transition to python but which Grass version (windows installer) does NVIZ work fully. Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Grass 6.4.2 RC1-1 Startup error on Windows XP
I keep getting the following error when attempting to start grass on Windows XP:- GRASS 6.4 Traceback (most recent call last): File c:/Program Files/GRASS 6.4RC1/etc/wxpython/wxgui.py, line 65, in module from gui_modules import layertree File C:\Program Files\GRASS 6.4RC1\etc\wxpython\gui_modules\layertree.py, line 44, in module import toolbars File C:\Program Files\GRASS 6.4RC1\etc\wxpython\gui_modules\toolbars.py, line 51, in module from icon import Icons File C:/Program Files/GRASS 6.4RC1\etc\wxpython\icons\icon.py, line 433, in module 'scriptSave' : MetaIcon(img = iconSet['script-save'], KeyError: 'script-save' GRASS 6.4 g.gui wxpython Traceback (most recent call last): File c:/Program Files/GRASS 6.4RC1/etc/wxpython/wxgui.py, line 65, in module from gui_modules import layertree File C:\Program Files\GRASS 6.4RC1\etc\wxpython\gui_modules\layertree.py, line 44, in module import toolbars File C:\Program Files\GRASS 6.4RC1\etc\wxpython\gui_modules\toolbars.py, line 51, in module from icon import Icons File C:/Program Files/GRASS 6.4RC1\etc\wxpython\icons\icon.py, line 433, in module 'scriptSave' : MetaIcon(img = iconSet['script-save'], KeyError: 'script-save' Launching 'wxpython' GUI in the background, please wait ... GRASS 6.4 I downloaded the following binaries from;- http://grass.fbk.eu/grass64/binary/mswindows/native/ Any ideas would be grateful -- Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Using vector data in Grass 6.4.0 on WinXP Pro
Hi, I have a vector layer of 10 polys and have noticed some odd behaviour. There is a long delay of a couple of minutes on opening the attribute table (i'm accessing (Postgresql/postgis locally at present). Also when the layer is open applying a simple query in the attribute table manager fails to get any results it just appears to freeze (not responding). It is does not zoom to attribute like other layers but zooms out to map extent instead. I have checked topology with v.build and it seems all okay. Does anyone else get this with larger vector layers? -- Paul J. Shapley ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user