Re: [GRASS-user] Re: [Qgis-user] QGIS GRASS
Micha Silver wrote: So, I agree with Markus and Paulo that QGIS should stay, on the surface, an easy to use viewer, editor, and map designer for GIS data, with the heavy lifting done thru plugins to other software. Micha And in so doing would be much more in harmony with the age-old Unix philosophy of do one thing and do it well. -- Eric Patton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] QGIS GRASS for online use
jgarcia wrote: What is the way you would advise to operate online (for any user with Internet Access) against a common GIS server with the full (or at least a good deal of) capabilities of GRASS? perhaps PyWPS or OpenLayers could help? http://pywps.wald.intevation.org/ http://grass.osgeo.org/screenshots/web.php Is it necessary to use mapserver, or is it possible to use a GRASS a a geodata server, and directly access them through QGIS remote clients? If mapserver is still required, do you know if there are some future plans to directly access GRASS as a server from remote GRASS or QGIS clients? both GRASS and QGIS have WMS input ability, now GDAL does as well. additionally QGIS has an easy to use mapsever output plugin AFAIK. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Getting rows,cols of a region in a script
On Mon, 2009-05-04 at 16:23 +0200, Nikos Alexandris wrote: I need the number of rows and columns to use with v.mkgrid. Currently I use the following, rather long and complex command: g.region -g | grep o | cut -d= -f2 | tr \n , | cut -d, -f1,2 Em.. that's meant to be used in a loop within a script for several maps. Sorry for not making it clear. Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.whar.vect does not respond
Since last update of grass6_devel (a few days back) v.what.vect does not work anymore!? It starts counting some percentages and then... that was it. It just stays there... Anyone had/has issues with it before I spend time on trying to debug? Kind regards, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Getting rows,cols of a region in a script
Hi, 2009/5/4 Martin Landa landa.mar...@gmail.com: [...] ret = grass.read_command('g.region', flags = 'g') reg = grass.parse_key_val(ret) in trunk and devbr6 a bit easier to parse command output reg = grass.read_command('g.region', flags = 'g', parse = True) Martin -- 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] Error in sqlite3_step()
On Mon, May 4, 2009 at 3:20 PM, Nikos Alexandris nikos.alexand...@felis.uni-freiburg.de wrote: On Mon, 2009-05-04 at 15:07 +0200, Nikos Alexandris wrote: [...] Count another error here: Error in sqlite3_step(): database is locked All I am trying to do is: v.what.rast vector=points raster=FINAL column=test WoW :-O All of a sudden, I can't copy vector maps, drop or add columns anymore. I checked dbase connections and all looks fine. The problem started after trying to use v.what.rast and then sqlitebrowser to remove some column. Could it be a sqlitebrowser side-effect? Perhaps... I found http://www.sqlite.org/lockingv3.html With lsof you can find out who's using a file or a process. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.whar.vect does not respond
On Mon, May 4, 2009 at 5:19 PM, Nikos Alexandris nikos.alexand...@felis.uni-freiburg.de wrote: Since last update of grass6_devel (a few days back) v.what.vect does not work anymore!? It starts counting some percentages and then... that was it. It just stays there... Anyone had/has issues with it before I spend time on trying to debug? I tried in Spearfish: g.copy vect=archsites,myarchsites v.db.addcol myarchsites col=owner varchar(25) v.what.vect myarchsites qvect=fields column=owner qcolumn=label # verification: v.db.select myarchsites No problem... Add -x to the first line in the v.what.vect script to see the script execution stepwise. Then you'll see where is stops: http://grass.osgeo.org/wiki/GRASS_Debugging#Shell_script_debugging Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Creating LiDAR first pulse vector map
I am still fighting with the LiDAR tools, but I feel like I am close... I thought I had sorted the v.lidar.growing problem as it appears to run when using sqlite backend. But the resulting output was only a small portion of the output? For v.lidar.growing, is there any special way to create the first pulse vector map? I feel this may be where my problem comes from. (Currently I was just using the raw data from the lidar text file which includes X Y Z and intensity readings for all points.) Thanks -- Jack Lonsdale, ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] passing output from GRASS to R.
On Mon, May 4, 2009 at 8:57 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Dear all, I am running R inside GRASS 6.4 and I can see a list of all my raster files when I type (from the R session) system(g.list rast) But I would like to pass the list of rasters to a vector object in R. I tryed innocently this, but not worked: my.rast.list-system(g.list rast) Try this: rasterfiles - read.table(textConnection(system(g.mlist rast,intern = TRUE))) str(rasterfiles) 'data.frame': 42 obs. of 1 variable: $ V1: Factor w/ 42 levels aspect,basin_50K,..: 31 1 2 3 4 5 6 7 8 9 ... rasterfiles V1 1myfacility 2aspect 3 basin_50K 4 boundary_county_500m ... I also would like to define a pattern of raster names to the g.list command, like only return that raster names that have _output_ on the name. So: rasterfiles - read.table(textConnection(system(g.mlist rast pattern=*_output_*,intern = TRUE))) str(rasterfiles) By the way, I am running GRASS inside MSYS, under windows. Please report if above works (as it does on Linux). Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error in sqlite3_step()
On 04/05/09 15:20, Nikos Alexandris wrote: On Mon, 2009-05-04 at 15:07 +0200, Nikos Alexandris wrote: [...] Count another error here: Error in sqlite3_step(): database is locked All I am trying to do is: v.what.rast vector=points raster=FINAL column=test WoW :-O All of a sudden, I can't copy vector maps, drop or add columns anymore. I checked dbase connections and all looks fine. The problem started after trying to use v.what.rast and then sqlitebrowser to remove some column. Could it be a sqlitebrowser side-effect? IIRC, if you use sqlitebrowser, you need to save your changes, before being able to access the db from GRASS. BTW, has anyone compared sqlitebrowser to the sqlite manager add-on to firefox [1,2] ? After very rough first trials I have the feeling that the latter actually has a nicer feel to it, but need to do more testing... Moritz [1] https://addons.mozilla.org/fr/firefox/addon/5817 [2] http://code.google.com/p/sqlite-manager/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error in sqlite3_step()
[...] Moritz: IIRC, if you use sqlitebrowser, you need to save your changes, before being able to access the db from GRASS. The fact is that I couldn't do any changes to save them! BTW, has anyone compared sqlitebrowser to the sqlite manager add-on to firefox [1,2] ? After very rough first trials I have the feeling that the latter actually has a nicer feel to it, but need to do more testing... Not a real comparison but I tend to agree that it feels smoother. Nikos [1] https://addons.mozilla.org/fr/firefox/addon/5817 [2] http://code.google.com/p/sqlite-manager/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.whar.vect does not respond
On Mon, 2009-05-04 at 21:45 +0200, Markus Neteler wrote: On Mon, May 4, 2009 at 5:19 PM, Nikos Alexandris nikos.alexand...@felis.uni-freiburg.de wrote: Since last update of grass6_devel (a few days back) v.what.vect does not work anymore!? It starts counting some percentages and then... that was it. It just stays there... Anyone had/has issues with it before I spend time on trying to debug? I tried in Spearfish: g.copy vect=archsites,myarchsites v.db.addcol myarchsites col=owner varchar(25) v.what.vect myarchsites qvect=fields column=owner qcolumn=label # verification: v.db.select myarchsites No problem... Add -x to the first line in the v.what.vect script to see the script execution stepwise. Then you'll see where is stops: http://grass.osgeo.org/wiki/GRASS_Debugging#Shell_script_debugging Markus Just a quick feedback before doing real-testing: I have some point vectors (created last week) and v.what.vect indeed works fine on/with them. Now, if I take one of those point vectors (g.copy), drop a column just for fun (v.db.dropcol), add it back again (v.db.addcol column='test integer') v.what.vect does not work anymore on the _newly_ added column. Will test on spearfish data to see if I can replicate this. Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.whar.vect does not respond
Nikos: Since last update of grass6_devel (a few days back) v.what.vect does not work anymore!? It starts counting some percentages and then... that was it. It just stays there... Markus N: I tried in Spearfish: g.copy vect=archsites,myarchsites v.db.addcol myarchsites col=owner varchar(25) v.what.vect myarchsites qvect=fields column=owner qcolumn=label # verification: v.db.select myarchsites No problem... Yep, no problem at all with Spearfish. Add -x to the first line in the v.what.vect script to see the script execution stepwise. Then you'll see where is stops: http://grass.osgeo.org/wiki/GRASS_Debugging#Shell_script_debugging I added -x but doesn't look very helpful as it get's (if it really stucks) after counting the progress-percentages of whatever process they refer to. ( For some modules I think it would be nice, if it is possible, to have sort of a progress indicator... just to know it's running and have a rough estimation of how long it will take to complete. ) The problem appears in a (at least one) map which I recreated today since it was not working at all (due to database locked message). The steps were: # rasterise some map r.to.vect input=burned_area_landsat_sample_1.rgb output=burned_area_sample_1_points feature=point # some info v.info burned_area_sample_1_points ---%--- Number of points: 99449 ---%--- # drop/add some column(s) v.db.dropcol map=burned_area_sample_1_points column=value v.db.dropcol map=burned_area_sample_1_points column=label v.db.addcol burned_area_sample_1_points column=gridcell integer # carry over from an area-vector to the point-vector the _cat_ numbers v.what.vect burned_area_sample_1_points col=gridcell qvect=sample_1_grid qcol=cat + [ -z /usr/local/grass-6.5.svn ] + [ burned_area_sample_1_points != @ARGS_PARSED@ ] + exec g.parser /usr/local/grass-6.5.svn/scripts/v.what.vect burned_area_sample_1_points col=gridcell qvect=sample_1_grid qcol=cat + [ -z /usr/local/grass-6.5.svn ] + [ @ARGS_PARSED@ != @ARGS_PARSED@ ] + exec v.distance from=burned_area_sample_1_points to=sample_1_grid column=gridcell to_column=cat upload=to_attr dmax=0.0 from_layer=1 to_layer=1 100% 100% 100% ## Here it waits... waits... maybe just wait? I am _sure_ it did not take too long some last time when I tried this with the same map(s) ## Currently I am watching the same process ran via strace! It throws out numerous lseek/lread messages. I'll give it some time but something is wrong because it takes too much time. Was there change recently in v.db.addcol or something related to v.db.dropcol? Nikos Regards, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.whar.vect does not respond
On Tue, 2009-05-05 at 01:13 +0200, Nikos Alexandris wrote: Nikos: Since last update of grass6_devel (a few days back) v.what.vect does not work anymore!? It starts counting some percentages and then... that was it. It just stays there... [...] Currently I am watching the same process ran via strace! It throws out numerous lseek/lread messages. I'll give it some time but something is wrong because it takes too much time. I was wrong about v.what.vect not working. Via strace I was ensured that something is going-on and I have to wait. Finally the process completed successfully. But I wonder why it took so much more than last time. Anyhow, sorry for the noise and thanks for the feedback. I learned about strace this time :-) Cheers, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] passing output from GRASS to R.
Hi Markus, Under windows I needed only small changes rasterfiles - read.table(textConnection(system(g.mlist.exe rast pattern=*output*, intern = TRUE))) Thanks a lot, milton 2009/5/4 Markus Neteler nete...@osgeo.org On Mon, May 4, 2009 at 8:57 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Dear all, I am running R inside GRASS 6.4 and I can see a list of all my raster files when I type (from the R session) system(g.list rast) But I would like to pass the list of rasters to a vector object in R. I tryed innocently this, but not worked: my.rast.list-system(g.list rast) Try this: rasterfiles - read.table(textConnection(system(g.mlist rast,intern = TRUE))) str(rasterfiles) 'data.frame': 42 obs. of 1 variable: $ V1: Factor w/ 42 levels aspect,basin_50K,..: 31 1 2 3 4 5 6 7 8 9 ... rasterfiles V1 1myfacility 2aspect 3 basin_50K 4 boundary_county_500m ... I also would like to define a pattern of raster names to the g.list command, like only return that raster names that have _output_ on the name. So: rasterfiles - read.table(textConnection(system(g.mlist rast pattern=*_output_*,intern = TRUE))) str(rasterfiles) By the way, I am running GRASS inside MSYS, under windows. Please report if above works (as it does on Linux). Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: [GRASS-dev] Seeking old version grass_src421v21.tar.gz
Markus, I have grass42src.tar.Z; will that do? I have the pdf manuals from that time, too… Tom Markus Neteler wrote: Hi, for historical reasons (think copyright, nostalgia, ...) I am seeking the good old GRASS 4.2.1 source code which I unfortunately removed in 2002 due to disk space limitations. Anyone having a copy? Any 4.2.1 subversion is fine. thanks Markus ___ grass-dev mailing list grass-...@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-dev -- Thomas E Adams National Weather Service Ohio River Forecast Center 1901 South State Route 134 Wilmington, OH 45177 EMAIL: thomas.ad...@noaa.gov VOICE: 937-383-0528 FAX:937-383-0033 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.whar.vect does not respond
Nikos wrote: + exec v.distance from=burned_area_sample_1_points to=sample_1_grid column=gridcell to_column=cat upload=to_attr dmax=0.0 from_layer=1 to_layer=1 100% 100% 100% ## Here it waits... waits... maybe just wait? I am _sure_ it did not take too long some last time when I tried this with the same map(s) ## Currently I am watching the same process ran via strace! It throws out numerous lseek/lread messages. I'll give it some time but something is wrong because it takes too much time. ok so v.distance is taking a very long time. can you try with v.what.vect --vebose ? that should give some more messages from v.distance instead of just 100% 100% 100%. also, you can turn on the debug messages: g.gisenv set=DEBUG=1 turn it up to 2,3,.. to get more messages. set it back to 0 to turn them off. I see one in v.distance G_debug(3, SQL: %s, db_get_string(stmt)); maybe more DB troubles? Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user