[GRASS-user] wxGUI Graphical Modeler - variables

2011-09-26 Thread Luís Ferreira
Hello.

I have some selected lines in a command line script:

#! /bin/sh -x

MYMAPSET=`g.gisenv -s get=MAPSET`
targetdir=/home/lf/TEMP/

LIST=`g.mlist type=rast separator=, pattern=clip_toar* mapset=
$MYMAPSET`
echo $LIST
i.group group=group input=$LIST

r.out.gdal -c -f input=group format=HFA type=Float32 nodata=0 output=
$targetdir/$MYMAPSET.clip_toar.img

g.remove group=group

#



In wxGUI Graphical Modeler, how can I  use the variables inside the
commands options? I'm trying using string data type variables 

declaring variables

MYMAPSET
default value `g.gisenv -s get=MAPSET`

targetdir 
default value /home/lf/TEMP/


then using in r.out.gdal options

r.out.gdal -c -f input=group format=HFA type=Float32 nodata=0 output=%
targetdir/%MYMAPSET.clip_toar.img
or 
$targetdir/$MYMAPSET.clip_toar.img


Without results.

Anyone have ideas how to solve this problem? Is it possible in Graphical
Modeller?

Thank you.

Luís Ferreira
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] 3D vectorial editing

2010-07-13 Thread Luís Ferreira
Hello,

At the moment, apparently, the unique implemented 3D vectorial GRASS
data type is the face (3D area).
We can transform vector geometry from raster using v.extrude or v.drape,
or transform from 2D to 3D using v.to.3d, v.to.db and v.in.db.
Is there any way to edit directly the z coordinate of vertices and nodes
in GRASS vectorial data?

Thank you,
Luís Ferreira



___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: [GRASS-dev] Fw: digitize toolbar error

2010-06-20 Thread Luís Ferreira
I'm not sure, but probably the problem is connected with the version of
swig you are using on compiling. I use swig version 1.3.36 and don't
have problems with wxpython gui vector digitizer.

Luís


On Dom, 2010-06-20 at 01:26 -0700, Mohammed Rashad wrote:
 
  
 --
 Rashad
 
 
 
 - Forwarded Message 
 From: Mohammed Rashad rasha...@yahoo.com
 To: grass-...@lists.osgeo.org
 Sent: Thu, June 17, 2010 8:53:22 PM
 Subject: digitize toolbar error
 
 I already send a mail regarding digitization tool error and didnt
 recieve any response
 so I am sending this mail again.
 
 When selecting digitize from mapdisplay window toolbars combobox
 and select a new vector map and given a name for newVectorMap
 I am getting this error
 
 Unable to initialize display driver of vector digitizer. See 'Command
 output' for details.
 
 Reason: message=_(Unable to initialize display driver of vector 
 
 Traceback (most recent call last):
   File
 /usr/local/grass-6.5.svn/etc/wxpython/gui_modules/toolbars.py, line
 1130, in StartEditing
 message=_(Unable to initialize display driver of vector 
 DigitError: unprintable DigitError object
 
 
 what i am missing.
  
 --
 Rashad
 
 
 
 
 
 ___
 grass-dev mailing list
 grass-...@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-dev


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] r.stream : vetor of horton order

2010-06-19 Thread Luís Ferreira
Take a look at http://grass.osgeo.org/wiki/GRASS_AddOns, on r.stream.*
collection.

Luís

On Sáb, 2010-06-19 at 15:33 -0300, Milton Cezar Ribeiro wrote:
 Dear all,
 
 I extracted the drainage and direction from a DEM using
 r.stream.extract.
 r.stream.extract --o elevation=aster_dem_wgs84_UTM threshold=50
stream_vect=aster_dem_wgs84_UTM_vect_th00050 
stream_rast=aster_dem_wgs84_UTM_rast_th00050 
direction=aster_dem_wgs84_UTM_rast_th00050_direction
 
 How can I obtain a vector of Horton order?!
 
 cheers
 
 milton
 
 
 ___
 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] i.rectify error

2010-06-08 Thread Luís Ferreira
On 

Ubuntu 10.04 x86_64
GRASS 6.4 SVN rev 42488

when trying to rectify an raster  group, the following error appears:

GRASS 6.4.0svn (grass_32629):~  i.rectify
*** glibc detected *** i.rectify: munmap_chunk(): invalid pointer:
0x01c08280 ***
=== Backtrace: =
/lib/libc.so.6(+0x775b6)[0x7fa63fe905b6]
/lib/libc.so.6(+0x9b039)[0x7fa63feb4039]
/lib/libc.so.6(+0x9b1b9)[0x7fa63feb41b9]
/usr/lib/grass-6.4.0svn/lib/libgrass_gis.so(G_date+0x20)[0x7fa640a6a954]
/usr/lib/grass-6.4.0svn/lib/libgrass_gis.so(G_short_history
+0x19)[0x7fa640a77c2e]
/usr/lib/grass-6.4.0svn/lib/libgrass_gis.so(+0x1a95a)[0x7fa640a6295a]
/usr/lib/grass-6.4.0svn/lib/libgrass_gis.so(G_close_cell
+0x7c)[0x7fa640a61f10]
i.rectify(write_map+0x1ce)[0x406c1e]
i.rectify(rectify+0x498)[0x4061c4]
i.rectify(exec_rectify+0x30d)[0x404519]
i.rectify(main+0x7d8)[0x405328]
/lib/libc.so.6(__libc_start_main+0xfd)[0x7fa63fe37c4d]
i.rectify[0x402b19]
=== Memory map: 
0040-00408000 r-xp  08:01
400720 /usr/lib/grass-6.4.0svn/bin/i.rectify
00607000-00608000 r--p 7000 08:01
400720 /usr/lib/grass-6.4.0svn/bin/i.rectify
00608000-00609000 rw-p 8000 08:01
400720 /usr/lib/grass-6.4.0svn/bin/i.rectify
00609000-00619000 rw-p  00:00 0 
01bec000-01c36000 rw-p  00:00 0
[heap]
7fa63db42000-7fa63db4e000 r-xp  08:01
5345   /lib/libnss_files-2.11.1.so
7fa63db4e000-7fa63dd4d000 ---p c000 08:01
5345   /lib/libnss_files-2.11.1.so
7fa63dd4d000-7fa63dd4e000 r--p b000 08:01
5345   /lib/libnss_files-2.11.1.so
7fa63dd4e000-7fa63dd4f000 rw-p c000 08:01
5345   /lib/libnss_files-2.11.1.so
7fa63dd4f000-7fa63dd59000 r-xp  08:01
5349   /lib/libnss_nis-2.11.1.so
7fa63dd59000-7fa63df58000 ---p a000 08:01
5349   /lib/libnss_nis-2.11.1.so
7fa63df58000-7fa63df59000 r--p 9000 08:01
5349   /lib/libnss_nis-2.11.1.so
7fa63df59000-7fa63df5a000 rw-p a000 08:01
5349   /lib/libnss_nis-2.11.1.so
7fa63df5a000-7fa63df71000 r-xp  08:01
5339   /lib/libnsl-2.11.1.so
7fa63df71000-7fa63e17 ---p 00017000 08:01
5339   /lib/libnsl-2.11.1.so
7fa63e17-7fa63e171000 r--p 00016000 08:01
5339   /lib/libnsl-2.11.1.so
7fa63e171000-7fa63e172000 rw-p 00017000 08:01
5339   /lib/libnsl-2.11.1.so
7fa63e172000-7fa63e174000 rw-p  00:00 0 
...
etc

Any ideas? Bug?

Luís Ferreira
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Creation and managing Location/Mapset

2010-05-17 Thread Luís Ferreira


 
 1- Should I create 1 Location with 4 different mapsets (for each of
 the considered scales)? Or should I use 1 location, one mapset and 4
 different regions?

one location if all the analysis will be done with one coordinate system

one mapset = one region, one resolution

so, perhaps, is better you to use four different mapsets (four levels) 

 2- How can I define a region that is not a perfect square' Because in
 Set Region we have North, South, East and West Boundaries. And my
 regions are not perfect squares.

A region don't need to be square. And you can use an analysis mask.

Raster  Mask, or r.mask.


Luís Ferreira


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Using v.surf.rst

2010-04-08 Thread Luís Ferreira
Hello Eric

 
 there are points outside specified 2D/3D region--ignored 6 points
 
 ERROR2: zero points in the given region!
 
 ERROR: Input failed


 Keeping in mind I’m completely new to GRASS (and don’t understand a
 lot of the jargon), can you explain what I’m doing wrong?

Probably, you aren't defining the working region. See command
g.region.

Example: g.region -a vect=second...@permanent res=10 -p

This command will define the working region by the spatial extent of
secondTry, defining and aligning a 10 meters resolution with the
region extent and printing the result. 

 Also, can you recommend a tutorial that covers this type of conversion
 oriented toward beginners?

I strongly suggest the book  OPEN SOURCE GIS - A GRASS GIS Approach, 3rd
edition.

Luís Ferreira


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Re: v.breach addon

2009-12-15 Thread Luís Ferreira
Yes, it worked with sqlite db connection.

Luís


On Mon, 2009-12-14 at 15:52 -0500, M S wrote:
 Luis:
 
 I noticed you mentioned you had run it with sqlite and dbf database
 drivers.  Mine is mapped to sqlite currently.  I see some v.to.db
 commands in the script, but do I understand correctly that you didnt
 have issues using sqlite as backend?


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: v.breach addon

2009-12-14 Thread Luís Ferreira

I used v.breach with both SQLite and dbf, with ETRS89 metric negative
coordinates, from 20m to 2m resolution and worked fine.

System: Ubuntu 9.10 Karmic Koala, 64 bits, GRASS 6.4SVN.

But I remember to see some warnings like yours. Are these warnings
critical? Visually, the interpolated DEMs with v.surf.rst are coherent.


Luís


On Mon, 2009-12-14 at 04:57 +0100, Maciej Sieczka wrote:

 @Luís: I applied most your of patch. Thanks! Are you sure that's all
 errors? I keep on getting warnings like:
 
 WARNING: Unable to get point on line: cat = 1 offset = 2.831711 (line
   length = 0)
 
 at v.segment calls in the script, and a crash at v.in.ascii.



___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: v.breach addon

2009-12-06 Thread Luís Ferreira
Aftermath...

I had to change some strings to make v.breach work on GRASS 6.4svn
revision 39873, Ubuntu/Debian 9.10.

Luís Ferreira


-
line 310
v.parallel input=${OUTP}_grid_addcat_bnd_patch_break_ln_2l_tmp_${UNQ}
output=${OUTP}_segm_2l_r_tmp_${UNQ} distance=0.1 --o  /dev/null
and
line 312
v.parallel input=${OUTP}_grid_addcat_bnd_patch_break_ln_2l_tmp_${UNQ}
output=${OUTP}_segm_2l_l_tmp_${UNQ} distance=-0.1 --o  /dev/null

change to
v.parallel

input=${OUTP}_grid_addcat_bnd_patch_break_ln_2l_tmp_${UNQ}
output=${OUTP}_segm_2l_r_tmp_${UNQ} side=right distance=0.1 --o
 /dev/null
v.parallel

input=${OUTP}_grid_addcat_bnd_patch_break_ln_2l_tmp_${UNQ}
output=${OUTP}_segm_2l_l_tmp_${UNQ} side=left distance=0.1 --o
 /dev/null
--
line 388
v.db.addcol map=${OUTP} layer=2 'columns=x double, y double, z double,
z_breach double'  /dev/null

change to
v.db.addcol map=${OUTP} layer=2
'columns=x double precision, y double
precision, z double precision, z_breach
double precision'  /dev/null
---

line 474
db.execute input=$TMP.${PROG}.sql database=$dbs driver=$drv  /dev/null
change to
db.execute input=$TMP.${PROG}.sql database=${database}
driver=
${driver}  /dev/null

---

change $tbl to $table

--
v.clean tool=prune 

defaulting to treshold=0.0

--
line 464
change to
table=`v.db.connect -g ${OUTP} | awk '{print $2}'`
database=`v.db.connect -g ${OUTP} | awk '{print $4}'`
driver=`v.db.connect -g ${OUTP} | awk '{print $5}'`



make explicit field separator for awk
awk -F |



On Wed, 2009-12-02 at 21:35 +0100, Maciej Sieczka wrote:
 Luís Ferreira pisze:
  I want to make a hydrological corrected DEM, using v.surf.rst for 
  interpolation/aproximation and v.breach to introduce aditional 
  elevation points at streams.
  
  My spatial reference system is EPSG:3763 ETRS89 / Portugal TM06, 
  (metrical).
  
  As required in v.breach document: - each input vect line have a 
  unique category in layer 1; - input lines not cross themselves; - 
  input lines not constitute loops; - the input vector lines not stand 
  out of the input raster DEM; - the input lines direction is from 
  upstream to downstream.
  
  After creating a DEM using v.surf.rst I want to correct the false 
  sinks in narrow valleys and sparse points areas. Using v.breach I 
  have as result two output vectors for streams lines and points but 
  without entities or attributes. I tried with both PostgreSQL and 
  Sqlite database connections.
  
  For PostgreSQL connection:
  
  DBMI-Postgres driver error: Cannot select: SELECT cat,lcat,z,z_breach
   FROM topografia.vbreachpoints_2 ERROR:  column z does not exist 
  LINE 1: SELECT cat,lcat,z,z_breach FROM topografia.vbreachpoints_2 ^
   Unable to open select cursor: 'SELECT cat,lcat,z,z_breach FROM 
  topografia.vbreachpoints_2' Copy table failed Table vbreach_2 does 
  not exist in database host=localhost,dbname=grass_pp_muda_pttm06 
  The table vbreach_2 is now part of vector map vbreach and may be
   deleted or overwritten by GRASS modules DBMI-Postgres driver error:
   Cannot create index: create unique index vbreach_2_cat on vbreach_2
  ( cat ) ERROR:  relation vbreach_2 does not exist Cannot create
  index DBMI-Postgres driver error: Cannot grant on table: grant select
  on vbreach_2 to public ERROR:  relation vbreach_2 does not exist 
  Cannot grant privileges on table vbreach_2 Select privileges were 
  granted on the table
  
  Any ideas? If wanted I'll send the files needed for replication.
 
 Luís,
 
 Glad to see someone trying to use my stuff :).
 
 What GRASS version do you use? What OS? What was your command line, exactly?
 
 Maybe the culprit is the DB backend. I have never used the script with
 anything else than DBF. Try to switch to DBF in GRASS (db.connect). Or a
 GRASS version too new. I think I used it only against 6.3. There
 *might* have been changes to underlying GRASS modules.
 
 I've written the script, and used it the last time, about 2 years ago.
 It did the job fine then. I don't have a personal interest in developing
 it at the moment. I could provide you a commercial support if you like.
 That would include testing the script against your data, in a GRASS
 setup similar to yours, and fixing the tool if necessary.
 
 Please contact me in priv if you are interested. I don't charge much :).
 Or, if there's a GRASS hacker willing to do it for free - please go
 ahead, anybody.
 
 Best,
 Maciek
 



v.breach_changed
Description: application/shellscript

[GRASS-user] v.breach addon

2009-12-02 Thread Luís Ferreira
I want to make a hydrological corrected DEM, using v.surf.rst for
interpolation/aproximation and v.breach to introduce aditional elevation
points at streams.

My spatial reference system is EPSG:3763 ETRS89 / Portugal TM06,
(metrical).

As required in v.breach document:
 - each input vect line have a unique category in layer 1;
 - input lines not cross themselves;
 - input lines not constitute loops;
 - the input vector lines not stand out of the input raster DEM;
 - the input lines direction is from upstream to downstream.

After creating a DEM using v.surf.rst I want to correct the false sinks
in narrow valleys and sparse points areas. Using v.breach I have as
result two output vectors for streams lines and points but without
entities or attributes. I tried with both PostgreSQL and Sqlite database
connections.

For PostgreSQL connection:

DBMI-Postgres driver error:
Cannot select:
SELECT cat,lcat,z,z_breach FROM topografia.vbreachpoints_2
ERROR:  column z does not exist
LINE 1: SELECT cat,lcat,z,z_breach FROM
topografia.vbreachpoints_2
^
Unable to open select cursor: 'SELECT cat,lcat,z,z_breach FROM
topografia.vbreachpoints_2'
Copy table failed
Table vbreach_2 does not exist in database
host=localhost,dbname=grass_pp_muda_pttm06
The table vbreach_2 is now part of vector map vbreach and may be
deleted or overwritten by GRASS modules
DBMI-Postgres driver error:
Cannot create index:
create unique index vbreach_2_cat on vbreach_2 ( cat )
ERROR:  relation vbreach_2 does not exist
Cannot create index
DBMI-Postgres driver error:
Cannot grant on table:
grant select on vbreach_2 to public
ERROR:  relation vbreach_2 does not exist
Cannot grant privileges on table vbreach_2
Select privileges were granted on the table

Any ideas?
If wanted I'll send the files needed for replication.

Thank you

Luís Ferreira


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] v.db.join - created fields with only one character long

2009-11-22 Thread Luís Ferreira
When trying to join a vector table with another table, the created
fields are only type character(1). The original fields to import vary
from 1 to 100 characters.

Is this a bug or I'm not doing something previously necessary?

GRASS 6.4 rev 39766 from SVN
postgresql 8.4.1 from repository
Ubuntu 9.10 

Luís Ferreira




v.db.join map=cos90_cl...@cos90 column=leg otable=grass.cos90rcl
ocolumn=codigo

DBMI-Postgres driver error:
Cannot execute:
UPDATE grass.cos90_clean SET codigo=(SELECT codigo
FROM grass.cos90rcl WHERE
grass.cos90rcl.codigo=grass.cos90_clean.leg)
ERROR:  value too long for type character(1)
ERROR: Error while executing: 'UPDATE grass.cos90_clean SET
codigo=(SELECT codigo
FROM grass.cos90rcl WHERE grass.cos90rcl.codigo=grass.cos90_clean.leg)'

DBMI-Postgres driver error:
Cannot execute:
UPDATE grass.cos90_clean SET leg2=(SELECT leg2
FROM grass.cos90rcl WHERE
grass.cos90rcl.codigo=grass.cos90_clean.leg)
ERROR:  value too long for type character(1)
ERROR: Error while executing: 'UPDATE grass.cos90_clean SET leg2=(SELECT
leg2
FROM grass.cos90rcl WHERE grass.cos90rcl.codigo=grass.cos90_clean.leg)'

DBMI-Postgres driver error:
Cannot execute:
UPDATE grass.cos90_clean SET nivel1=(SELECT nivel1
FROM grass.cos90rcl WHERE
grass.cos90rcl.codigo=grass.cos90_clean.leg)
ERROR:  value too long for type character(1)
ERROR: Error while executing: 'UPDATE grass.cos90_clean SET
nivel1=(SELECT nivel1
FROM grass.cos90rcl WHERE grass.cos90rcl.codigo=grass.cos90_clean.leg)'

etc.

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] v.dissolve and SQLite

2009-11-21 Thread Luís Ferreira
After:

# v.patch -e input=x...@cos90,x...@cos90,x...@cos90,x...@cos90,x...@cos90
output=cos90_patch

# v.clean input=cos90_pa...@cos90 output=cos90_clean
type=point,line,boundary,centroid,area tool=break,snap,rmdupl,bpol
thresh=0.01

# v.dissolve input=cos90_cl...@cos90 output=cos90_diss column=LEG


...gives me this error:

DBMI-SQLite driver error:
Cannot step:
SQL logic error or missing database
ERROR: Unable to fetch data from table cos90_clean
ERROR: Unable to open vector map cos90_diss_11...@cos90 on level 2.
Try to rebuild vector topology by v.build.
ERROR: Vector map cos90_diss not found in the current mapset


An similar error is documented. 
http://trac.osgeo.org/grass/ticket/548

Any idea?
I'm going to try the same operations with a postgres connection.

Luís Ferreira

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] v.dissolve and SQLite

2009-11-21 Thread Luís Ferreira
Yes, connecting to postgresql solved the problem.

Luís

On Sat, 2009-11-21 at 20:48 +0100, Martin Landa wrote:
 Hi,
 
 2009/11/21 Luís Ferreira lferreira7...@gmail.com:
  An similar error is documented.
  http://trac.osgeo.org/grass/ticket/548
 
  Any idea?
 
 right, the bug still waits to be be fixed. PostgreSQL backend should work.
 
 Martin
 


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] mosaic raster files

2009-07-02 Thread Luís Ferreira
30 minutes for 35 tiles.

I tried to use some commercial software: 2 hours until 10 files with
subsequent error :D. 

g.region rast=`g.mlist pat=AST*dem sep=,`

r.series input=`g.mlist pat=AST*dem sep=,` output=aster_dem
method=average

Qui, 2009-07-02 às 08:28 -0700, Jamie Adams escreveu:
 If you have GDAL 1.6 installed, you can use gdalbuildvrt to create a
 VRT of the files, and then link to this using r.external.  This way
 you have a mosaic without duplicating the data.  
 
 On Thu, Jul 2, 2009 at 7:53 AM, Milton Cezar Ribeiro
 miltinho.astrona...@gmail.com wrote:
 Hi Guys,
  
 I will try first download the latest SVN. In fact I am running
 a compiled version of 6.4R4, and I remember, on that time, few
 programs (like r.patch if I am not wrong) not compiled
 properly. Case it not work, I try r.series.
  
 Thanks a lot,
  
 milton
 
 
 2009/7/2 Hamish hamis...@yahoo.com
 
 
 Milton Cezar Ribeiro wrote:
  I just downloaded some tiles from GDEM/NASA and
 imported on GRASS.
  Now I would like to generate a mosaic with the
 tiles,
  but unfortunatelly r.patch
 ( http://grass.itc.it/gdp/html_grass63/r.patch.html )
  is not compiled for GRASS 6.4R4 (windows / vista).
 
 
 Which installer?
 
 really there is nothing special about r.patch so it is
 a bad sign that
 it would be missing. what if you install RC5 or the
 SVN snapshot binary
 from June 2?
 
 
  So is there a way of I join a list of rasters on a
 unique mosaic file?
 
 
 r.series and r.mapcalc should both be able to do the
 job, but not as
 simply.
 
 use 'g.region rast=tile1,tile2,tile3,tile5' first to
 set the computational
 region to the outside bounds of all the tiles listed.
 
 
 Hamish
 
 
 
 
 
 
 
 ___
 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


Re: [GRASS-user] Call script automatically when launching Grass-GIS]

2008-09-21 Thread Luís Ferreira

 Dom, 2008-09-21 às 16:53 +1000, Richard Chirgwin escreveu:
  Hi all,
  
  I have a few user machines with lots of maps to deal with, and I want to
  centralize the map catalogues. So what I would like to do is:
  
  - start Grass-GIS in text mode in a location (/home/user/grassdata/location)
  - automatically launch the script (which is a simple bash loop on
  g.mlist, calling on v.info and r.info to pull the metadata into a text file)
  - exit at the end.
  
  The script works fine if I execute it in Grass-GIS but I would like to
  have it run without intervention. That way, I can schedule the catalogue
  to run overnight on all machines. So is there a way to pass the script
  to Grass-GIS to run 'hands-off'?

If you are using GNU/Linux see cron, gnome-schedule, kcron and bcron for 
scheduled tasks.

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] r.fillnulls (GRASS BOOK 3rd edition)

2008-09-18 Thread Luís Ferreira
Hi

On page 123 of the book Open Source Gis - A GRASS GIS Aproach, when
performing

GRASS 6.4.svn (nc_spm_08):~  r.fillnulls elev_srtm_30m_null
out=elev_srtm_30mfilled

it gives me the errors

/usr/lib/grass-6.4.svn/scripts/r.fillnulls: eval: line 82: unexpected
EOF while looking for matching `''
/usr/lib/grass-6.4.svn/scripts/r.fillnulls: eval: line 83: syntax error:
unexpected end of file
/usr/lib/grass-6.4.svn/scripts/r.fillnulls: line 83: GISDBASE: parameter
null or not set

I'm starting using GRASS. Any suggestion to solve this?

Thank you




___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] r.fillnulls (GRASS BOOK 3rd edition)

2008-09-18 Thread Luís Ferreira
Qui, 2008-09-18 às 15:40 -0700, Hamish escreveu:
  GRASS 6.4.svn (nc_spm_08):~  r.fillnulls \
elev_srtm_30m_null out=elev_srtm_30mfilled
  
  it gives me the errors
  
  /usr/lib/grass-6.4.svn/scripts/r.fillnulls: eval: line 82:
  unexpected EOF while looking for matching `''
 
 
 goes g.gisenv --help or just g.gisenv work from the command prompt?
 
 it sounds like GRASS is not fully installed correctly.

 Hamish

I had compiled GRASS 6.4 from 
https://svn.osgeo.org/grass/grass/branches/develbranch_6 grass6_devel.

Kernel  : Linux 2.6.26-1-vserver-amd64 (x86_64)
Compiled: #1 SMP Thu Aug 28 13:09:10 UTC 2008
C Library   : GNU C Library version 2.7 (stable)
Distribution: Debian GNU/Linux lenny/sid

Code lines in /usr/lib/grass-6.4.svn/scripts/r.fillnulls
82 eval `g.gisenv`
83 : ${GISBASE?} ${GISDBASE?} ${LOCATION_NAME?} ${MAPSET?}
84 LOCATION=$GISDBASE/$LOCATION_NAME/$MAPSET

Using g.gisenv and g.gisenv --help:

GRASS 6.4.svn (nc_spm_08):~  g.gisenv
GISDBASE=/home/lf/GRASSDATA
LOCATION_NAME=nc_spm_08
MAPSET=user1
GRASS_HTML_BROWSER=/usr/lib/iceweasel/iceweasel
GRASS_DB_ENCODING=iso8859-1
GRASS_GUI=wxpython


GRASS 6.4.svn (nc_spm_08):~  g.gisenv --help

Descrição:
 Sai e modifica as atuais variáveis de usuário do GRASS.

Palavras chave:
 general

Uso:
g.gisenv [get=VARIABLE] [set=VARIABLE=value] [store=string]
   [--verbose] [--quiet]

Flags:
 --v   Saída do módulo verbosa
 --q   Saída do módulo silenciosa

Parâmetros:
get   Variável do GRASS para capturar
set   Variável do GRASS para ajustar
  store   Onde as variáveis do GRASS são guardadas
  opções: gisrc,mapset
  padrão: gisrc

 Using gis.m gives:


/usr/lib/grass-6.4.svn/scripts/r.fillnulls: eval: line
82: unexpected EOF while looking for matching `''
/usr/lib/grass-6.4.svn/scripts/r.fillnulls: eval: line
83: syntax error: unexpected end of file
Locating and isolating NULL areas...

Reading input raster map [EMAIL PROTECTED]...


Finding buffer zones...


Writing output raster map r_fillnulls_12974.buf...


Creating interpolation points...

Extracting points...

Building topology for vector map
vecttmp_fillnulls_12974...
egistering primitives: 
Building areas: 
Número de nós :   2187
Número de primitivas:   2187
Número de pontos:   2187
Número de linhas :   0
Número de fronteiras:   0
Número de centróides :   0
Número de áreas :   0
Número de ilhas :   0

Topology was built

0 areas built  
0 isles built
Anexando ilhas:
Anexando centróides:
r.to.vect completo. 

2187 primitives registered
2187 vertices registered
Interpolating 2187 points

Note: The following warnings may be ignored.

Removing raster MASK

Using segmentation for interpolation...

Porcentagem completada: 

Carregando dados da tabela de atributos ... 

Lendo linhas do mapa vetorial ... 

Lendo nós do mapa vetorial


Processing all selected output files
Porcentagem completada: 

will require 1818544 bytes of disk space for temp files

The number of points from vector map is 2187
The number of points outside of 2D/3D region 0
The number of points being used is 2187

há uma faixa com dados insuficientes

demorando demais para encontrar pontos para

Re: [GRASS-user] r.fillnulls (GRASS BOOK 3rd edition)

2008-09-18 Thread Luís Ferreira
Qui, 2008-09-18 às 18:33 -0700, Hamish escreveu:
 Luís Ferreira wrote:
  Using g.gisenv and g.gisenv --help:
  
 GRASS 6.4.svn (nc_spm_08):~  g.gisenv
 GISDBASE=/home/lf/GRASSDATA
 LOCATION_NAME=nc_spm_08
 MAPSET=user1
 GRASS_HTML_BROWSER=/usr/lib/iceweasel/iceweasel
 GRASS_DB_ENCODING=iso8859-1
 GRASS_GUI=wxpython
 
 
 edit ~/.grassrc6 and remove the GRASS_HTML_BROWSER line.
 
 Besides not needing the quotes it doesn't belong in that file as it is
 an environment variable not a GIS system variable.
 
 running this before starting GRASS should get that set correctly:
  export GRASS_HTML_BROWSER=iceweasel
 
 (assuming you use bash)
 
 Hamish

Editing ~/.grassrc6 solved all the r.fillnulls error warnings in bash,
gis.m and g.gui.

Thank you very much

Luís Ferreira

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user