Re: [GRASS-user] grass startup cmd line
Thank you Hamish, As you suggest, I wrote a shell script containing the instruction g.gui workspace=... After declaring the GRASS_BATCH_JOB variable, and launching grass in text mode, the wxpython GUI starts and my_workspace loads correctly, but the problem is I can't get the control on the process ; I guess grass is waiting for instructions from /within/ the shell script ! Otherwise I tried to pipe the command like this : echo 'g.gui gui=wxpython workspace=my_workspace.gxw' | grass64 -text my_gisbase/my_location/my_mapset There again gui and workspace load properly but the grass prompt does not appear. Anyway, it's a minor problem... but it could be interesting to improve the capabilities of this concept of workspace by allowing to launch grass only by e.g. double-clicking on a .gxw icon. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] menu items
Koen Hufkens wrote: I'm cleaning up my grass scripts and making them user friendly using the g.parser module. However, I'm wondering if it is possible to add a menu item to the gis manager window in a similar fashion. I found d.menu but this doesn't do what I expected it to do. Any tips and hints on how to add on an extra menu item in the gis.m gis manager? You need to modify $GISBASE/etc/gm/gmmenu.tcl. Or alternatively, create $GISBASE/etc/xtnmenu.dat file, or dir/xtnmenu.dat for any dir in $GRASS_ADDON_ETC. However, I'm not entirely sure of the format; examine the code at the top of gmmenu.tcl for clues. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Calculating standard error of many maps?
Hamish wrote: I am trying it at the moment with 1002 maps and, apart from the fact that it takes some time (which I am sure it would in GRASS as well), it seems to be working. I will be looking at R.series again later, when I run the simulations again. But I assume that there will be serious problems with 25000 maps - I might have to do it with a sample of maps. Hi, for the standard stats available from r.univar n minimum maximum range mean mean of absolute values standard deviation variance variation coefficient sum (ie not median, quartiles) you can use r.series in a loop, eg 1000 maps at a time, and then sum up the resulting 25 summary maps before doing a little math to get the overall answers. (or lots of 100 maps x250, ...) I think you would need to add a new sum_of_squares method to r.series to calculate variance/stdev. see r.univar/stats.c print_stats() Or 25000 runs of r.mapcalc $map.sq = $map^2 I can add sum of squares/cubes/etc to r.series easily enough if there is a use for them. Regarding command-line limits: is it worth extending G_parser() to allow arguments to be read from a file? E.g. r.series @args.txt. It was quite common for this feature to be added to DOS ports of Unix programs, as DOS had a much lower limit than Unix on the maximum length of a command line. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Appropriate use of v.surf.rst (lidar data)
Wesley Roberts wrote: Dear Grass Users, I am having trouble with interpolating a canopy height model using v.surf.rst and some lidar data. The lidar point density ranges from 5 - 8 points per meter squared. I know I can import the data using r.in.xyz but I would prefer to use an interpolation technique as opposed to a spatial average. I chose r.surf.xyz based on the paper written by Helena Mitasova, Lubos Mitas, and Russell S. Harmon. While the applications differ the use of a large data set with many points means dealing with the same issues. I ran the algorithm using default values and the result was a smooth moon-like surface not really approximating a plantation forest canopy (see link 2) v.surf.rst input=areaone_4...@wesley layer=1 zcolumn=dbl_3 elev=Area1_4_def tension=40. segmax=40 npmin=300 dmin=0.05 dmax=0.25 zmult=1.0 Given that I am working with very dense data with a high number of points I then ran the algorithm with a smaller tension value (roughly 2.5) with an incrase in the maximum number of points per segment (decrese computation time) and I find that the interpolation is blocky in nature. (see link 3) v.surf.rst input=areaone_4...@wesley layer=1 zcolumn=dbl_3 elev=Area1_4_def6 tension=2.5 smooth=0.0001 segmax=100 npmin=300 dmin=0.05 dmax=0.25 zmult=1.0 Is the blockiness related to the increase in the maximum number of points in a segment i.e. segmax value? I have also been looking at using R to do the interpolation using gstat or automap. Does anyone on the list have any advice with regards to methods to use for the interpolation of lidar canopy returns, with a view to creating a 3-D model of the canopy. The website spatial-analyst (http://spatial-analyst.net) suggests the use of regression-kriging but I don't have any ancillary data besides height and intensity. I have uploaded pics of the study area and interpolation results to flickr (everyone should have access to these) Link 1 http://www.flickr.com/photos/35273...@n07/3271646498/ Link 2 http://www.flickr.com/photos/35273...@n07/3270854029/ Link 3 http://www.flickr.com/photos/35273...@n07/3270904761/ Many thanks for taking the time to read my post and for your valuable contributions Wesley (I am using Grass 6.3.0 on Ubuntu Hardy Heron) Wesley Roberts MSc. Researcher: Earth Observation (Ecosystems) Natural Resources and the Environment CSIR Tel: +27 (21) 888-2490 Fax: +27 (21) 888-2693 To know the road ahead, ask those coming back. - Chinese proverb Hi Wesley, I've tried various methods for processing LiDAR data. I have a few comments: 1) Play with the resolution of your map and try different settings. 2) r.in.xyz does give good results, so don't rule it out for first approximation, and to help decide the resolution that you need. 3) r.surf.nnbathy might do the trick, but I haven't tried it. 4) If I want to interpolate to a higher resolution that my data covers, I usually krige in gstat. I found that easier than learning how to use R on top of everything else. The examples on p348-353 of the GRASS Book are easily adapted to your own data. (http://www.springerlink.com/content/j21632/) Note that the default method calculates a global variogram, which can take a long time with LiDAR data. The examples on p53-54 of the gstat manual describe how you can limit the analysis to a local neighborhood or set a maximum number of points to include in each calculation and thus dramatically speed up the interpolation. (http://www.gstat.org/gstat.pdf) Later John -- Dr John Stevenson Postdoctoral Research Associate School of Earth, Atmospheric and Environmental Sciences Williamson Building (Room 2.42) University of Manchester Manchester M13 9PL, UK tel. +44(0)161 306 6585; fax. +44(0)161 306 9361; john.steven...@manchester.ac.uk ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] use commands window
On 10/02/09 19:00, aude Valade wrote: Hi, I am very new at GRASS (I am using version 6.3.0 on windows) and I cannot type anything in the command window. When I am in mark mode any try to type a letter produces an error sound. Assuming that you used the native Windows installer, the command window that opens when you launch GRASS is just an artefact of the way the grass session functions in Windows and cannot be used for command line entry. See here: http://grass.osgeo.org/grass63/binary/mswindows/native/#Launching%20GRASS for details about the different options for launcing grass. Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass startup cmd line
Hello Moritz, unfortunately, I tried it too, but I got bailed out, with several error statements that don't really spark off anything to me : Traceback (most recent call last): File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 70, in module import gui_modules.preferences as preferences File /usr/local/grass-6.4.0svn/etc/wxpython/gui_modules/preferences.py, line 799, in module globalSettings = Settings() File /usr/local/grass-6.4.0svn/etc/wxpython/gui_modules/preferences.py, line 488, in __init__ self.ReadSettingsFile() File /usr/local/grass-6.4.0svn/etc/wxpython/gui_modules/preferences.py, line 558, in ReadSettingsFile mapset_file = os.path.join(gisdbase, location_name, mapset_name, self.fileName) File /usr/lib/python2.5/posixpath.py, line 60, in join if b.startswith('/'): AttributeError: 'NoneType' object has no attribute 'startswith' I did not go further in this direction, am I wrong ? Thank you, Vincent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Problemas con winGRASS 6.4.0RC3 del instalador OSGeo4W
Hi, please consider this is international ML (for those who cannot speak Spanish). Martin 2009/2/11 Marilena Yeguez myeg...@hotmail.com: Saludos, me anime a descargar el instalador OSGeo4W que trae GRASS para Windows entre muchas otras cosas, tal como me recomendaron. Lo instalé y al abrir GRASS me dio un error al no estar predefinido el directorio de datos, estuve revisando a ver si encontraba una carpeta que el instalador creara para tal fin, pero no la halle, realmente la crea?...Decidí crear mi propio directorio, copie una locación de prueba, pero GRASS me arroja el siguiente error: Erro setting region (Problem with g.region?): child lilled:unknown signal Agradezco a quien pueda ayudarme... Marilena Get news, entertainment and everything you care about at Live.com. Check it out! Invite your mail contacts to join your friends list with Windows Live Spaces. It's easy! Try it! ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] menu items
On Wed, Feb 11, 2009 at 10:59 AM, Glynn Clements gl...@gclements.plus.com wrote: Koen Hufkens wrote: I'm cleaning up my grass scripts and making them user friendly using the g.parser module. However, I'm wondering if it is possible to add a menu item to the gis manager window in a similar fashion. I found d.menu but this doesn't do what I expected it to do. Any tips and hints on how to add on an extra menu item in the gis.m gis manager? You need to modify $GISBASE/etc/gm/gmmenu.tcl. Or alternatively, create $GISBASE/etc/xtnmenu.dat file, or dir/xtnmenu.dat for any dir in $GRASS_ADDON_ETC. However, I'm not entirely sure of the format; examine the code at the top of gmmenu.tcl for clues. ...added as new FAQ: http://grass.osgeo.org/wiki/GIS_manager Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] executing multiple commands at once
Hi- I am new to GRASS (using version 6.3.0), and am trying to figure out how to execute multiple commands at once in the command line. I am currently tring to import many LANDSAT image files as rasters, and it would make it a lot easier/faster if I could import more than one at a time. Most of the online examples I have found show that people are able to do this, but I am still not aware how to do so. Any tips would be much appreciated. Thanks, Renae Mackas -- View this message in context: http://n2.nabble.com/executing-multiple-commands-at-once-tp2311081p2311081.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] db.select and SELECT COUNT query
Hi, has anybody ever had this kind of message ? Working on grass6.4.0svn, considering a database connection through the pg driver ; in order to count records matching an attribute (attr) value I execute this command on a table (table1) linked to a vector map : echo select count(*) from table1 where attr=201 | db.select The result should be 3. And db.select returns : ATTENTION : column 'count' : type int8 (bigint) is stored as integer (4 bytes) some data may be damaged count 3 In the end, it returns the result, but what does this message mean ? (Of course the same instruction given at the psql prompt raises no special warning). Any idea ? Thank you in advance, Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] menu items
Updated the FAQ... Should work for most systems as described in the FAQ Koen Op 11-feb-09, om 18:03 heeft Markus Neteler het volgende geschreven: On Wed, Feb 11, 2009 at 10:59 AM, Glynn Clements gl...@gclements.plus.com wrote: Koen Hufkens wrote: I'm cleaning up my grass scripts and making them user friendly using the g.parser module. However, I'm wondering if it is possible to add a menu item to the gis manager window in a similar fashion. I found d.menu but this doesn't do what I expected it to do. Any tips and hints on how to add on an extra menu item in the gis.m gis manager? You need to modify $GISBASE/etc/gm/gmmenu.tcl. Or alternatively, create $GISBASE/etc/xtnmenu.dat file, or dir/xtnmenu.dat for any dir in $GRASS_ADDON_ETC. However, I'm not entirely sure of the format; examine the code at the top of gmmenu.tcl for clues. ...added as new FAQ: http://grass.osgeo.org/wiki/GIS_manager Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] db.select and SELECT COUNT query
Sorry for this /stupid/ question (might be time to go to sleep !) I realize it must be related to a truncature operated between pg and grass or something like it. As it was annoying in a recursive query script, I turned the command in this : echo select count(*) from table1 where attr=201\ | db.select 2/dev/null Good evening ;-) Vincent Le mercredi 11 février 2009 à 21:36 +0100, Vincent Bain a écrit : Hi, has anybody ever had this kind of message ? Working on grass6.4.0svn, considering a database connection through the pg driver ; in order to count records matching an attribute (attr) value I execute this command on a table (table1) linked to a vector map : echo select count(*) from table1 where attr=201 | db.select The result should be 3. And db.select returns : ATTENTION : column 'count' : type int8 (bigint) is stored as integer (4 bytes) some data may be damaged count 3 In the end, it returns the result, but what does this message mean ? (Of course the same instruction given at the psql prompt raises no special warning). Any idea ? Thank you in advance, Vincent. ___ 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] Open source GIS cataloging software?
Sorry for the cross-posting. I was wondering if anyone has a unix (preferably) or windows program that can spider a directory, recursively searching for raster and vector data, and create bounding box polygons for each raster/vector it finds, with attributes indicating the path-to-file. Thanks! --j -- Jonathan A. Greenberg, PhD Postdoctoral Scholar Center for Spatial Technologies and Remote Sensing (CSTARS) University of California, Davis One Shields Avenue The Barn, Room 250N Davis, CA 95616 Cell: 415-794-5043 AIM: jgrn307, MSN: jgrn...@hotmail.com, Gchat: jgrn307 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-dev] Re: [GRASS-user] Re: Calculating standard error of many maps?
Markus Neteler wrote: Regarding command-line limits: is it worth extending G_parser() to allow arguments to be read from a file? E.g. r.series @args.txt. It was quite common for this feature to be added to DOS ports of Unix programs, as DOS had a much lower limit than Unix on the maximum length of a command line. Just to understand - does this allow to exceed the command limits? Then it would be pretty helpful (as text files are easy to generate). The idea is to have G_parser() recognise arguments of the form @filename, and read the file, treating each line the same way that it would treat an argv[i]. IOW, the commands: r.series @args.txt ... and: r.series `cat args.txt` ... would be equivalent[1], except that the former wouldn't run the risk of exceeding any OS limit on the maximum length of a command line. [1] Actually, they would interpret whitespace differently. The logical implementation of @filename would treat each line as a single argument. This is probably the most useful behaviour if you're generating the arguments with commands, and eliminates the need to provide a quoting mechanism (so long as you don't need to allow newlines within an argument). -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Thiessen Polygons
Hamish a écrit : Kurt Springs wrote: Does anyone know how to do Thiessen Polygons in GRASS? I need to use them around various types of megalithic tomb sites that are in point vector files. v.voronoi there is some replacement test code to look at in grass-addons as well: http://trac.osgeo.org/grass/browser/grass-addons/vector/voronoi While we are on this topic, is there a way to get a weigthed voronoi diagram using grass ? The ability to rank a point to tune the area's influence would be great, for that purpose I've been using an arcgis'extension* but with grass it is not possible**. Is there a way to get a similar result ? *http://www.geog.unt.edu/~pdong/software.htm **http://osdir.com/ml/gis.grass.user/2004-04/msg00036.html Regards, MORREALE Jean Roc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
RE: [GRASS-user] Open source GIS cataloging software?
Recently, I came across this: http://spatialguru.com/ideas/data_cataloguing_background Bob Moskovitz Seismic Hazard Evaluation Project California Geological Survey http://gmw.consrv.ca.gov/shmp CONFIDENTIALITY NOTICE: This communication is intended only for the use of the individual or entity to which it is addressed. This message contains information from the State of California, California Geological Survey, which may be privileged, confidential and exempt from disclosure under applicable law, including the Electronic Communications Privacy Act. If the reader of this communication is not the intended recipient, you are hereby notified that any dissemination, distribution, or copying of this communication is strictly prohibited. -Original Message- From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-boun...@lists.osgeo.org]on Behalf Of Jonathan Greenberg Sent: Wednesday, February 11, 2009 1:16 PM To: grass-user@lists.osgeo.org; r-sig-...@stat.math.ethz.ch Subject: [GRASS-user] Open source GIS cataloging software? Sorry for the cross-posting. I was wondering if anyone has a unix (preferably) or windows program that can spider a directory, recursively searching for raster and vector data, and create bounding box polygons for each raster/vector it finds, with attributes indicating the path-to-file. Thanks! --j -- Jonathan A. Greenberg, PhD Postdoctoral Scholar Center for Spatial Technologies and Remote Sensing (CSTARS) University of California, Davis One Shields Avenue The Barn, Room 250N Davis, CA 95616 Cell: 415-794-5043 AIM: jgrn307, MSN: jgrn...@hotmail.com, Gchat: jgrn307 ___ 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] executing multiple commands at once
On Wed, Feb 11, 2009 at 9:17 PM, Renae Mackas rmac...@sfu.ca wrote: Hi- I am new to GRASS (using version 6.3.0), and am trying to figure out how to execute multiple commands at once in the command line. I am currently tring to import many LANDSAT image files as rasters, and it would make it a lot easier/faster if I could import more than one at a time. Most of the online examples I have found show that people are able to do this, but I am still not aware how to do so. Any tips would be much appreciated. I usually use a shell script for this, in a for loop I get in all files. Example: http://www.grassbook.org/examples_menu3rd.php - Scripts to bulk import LANDSAT-TM5/LANDSAT-TM7 scenes from GLCF Maryland into GRASS - download You can easily adapt those scripts to your needs. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Thiessen Polygons
On Wednesday 11 February 2009, MORREALE Jean Roc wrote: Hamish a écrit : Kurt Springs wrote: Does anyone know how to do Thiessen Polygons in GRASS? I need to use them around various types of megalithic tomb sites that are in point vector files. v.voronoi there is some replacement test code to look at in grass-addons as well: http://trac.osgeo.org/grass/browser/grass-addons/vector/voronoi While we are on this topic, is there a way to get a weigthed voronoi diagram using grass ? The ability to rank a point to tune the area's influence would be great, for that purpose I've been using an arcgis'extension* but with grass it is not possible**. Is there a way to get a similar result ? *http://www.geog.unt.edu/~pdong/software.htm **http://osdir.com/ml/gis.grass.user/2004-04/msg00036.html Regards, MORREALE Jean Roc What exactly is a 'weighted voronoi diagram' (the link does not work)? If you are after a weighted network connecting points, here is one approach: http://casoilresource.lawr.ucdavis.edu/drupal/node/698 Dylan -- Dylan Beaudette Soil Resource Laboratory http://casoilresource.lawr.ucdavis.edu/ University of California at Davis 530.754.7341 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Thiessen Polygons
On Wed, Feb 11, 2009 at 5:03 PM, Dylan Beaudette debeaude...@ucdavis.edu wrote: On Wednesday 11 February 2009, MORREALE Jean Roc wrote: Hamish a écrit : Kurt Springs wrote: Does anyone know how to do Thiessen Polygons in GRASS? I need to use them around various types of megalithic tomb sites that are in point vector files. v.voronoi there is some replacement test code to look at in grass-addons as well: http://trac.osgeo.org/grass/browser/grass-addons/vector/voronoi While we are on this topic, is there a way to get a weigthed voronoi diagram using grass ? The ability to rank a point to tune the area's influence would be great, for that purpose I've been using an arcgis'extension* but with grass it is not possible**. Is there a way to get a similar result ? *http://www.geog.unt.edu/~pdong/software.htm **http://osdir.com/ml/gis.grass.user/2004-04/msg00036.html Regards, MORREALE Jean Roc Strike that. Link appears to work and RTFM-ing the user guide... Dylan ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Thiessen Polygons
Dylan Beaudette wrote: Does anyone know how to do Thiessen Polygons in GRASS? I need to use them around various types of megalithic tomb sites that are in point vector files. v.voronoi there is some replacement test code to look at in grass-addons as well: http://trac.osgeo.org/grass/browser/grass-addons/vector/voronoi While we are on this topic, is there a way to get a weigthed voronoi diagram using grass ? The ability to rank a point to tune the area's influence would be great, for that purpose I've been using an arcgis'extension* but with grass it is not possible**. Is there a way to get a similar result ? *http://www.geog.unt.edu/~pdong/software.htm **http://osdir.com/ml/gis.grass.user/2004-04/msg00036.html Regards, MORREALE Jean Roc Now that I have read about 'weighted voronoi diagrams', I wonder if a combination of r.cost + r.mapcalc would solve this problem. Something along those lines is demonstrated here: http://casoilresource.lawr.ucdavis.edu/drupal/node/288 This example isn't quite what is requested, although using r.cost with start=point_i, and stop=neighbor_points (derived from v.delaunay / v.distance?) may work. It would then be a little more work to convert the weighted-distance rasters into polygons, and link back to the original attribute tables... but (hopefully) not outside the realm of possibility via a script. The problem with using r.cost is that you would need to know the cost for each cell before you have created the polygons. I think that the simplest accurate approach would be to modify r.grow.distance. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user