Re: [GRASS-user] R, Moran index, getis e ord.....
Maybe http://www.asdar-book.org/ can help you. Have a look in Additional materials. Robbie 2010/9/30 Gabriele N. gis...@libero.it Hello everyone. I have a raster which is the result of calculating an index of remote sensing (with decimals). I would like to calculate the Moran's index and the index of Getis and Ord To do this I installed R (following http://grass.fbk.eu/statsgrass/grass6_r_install.html). Then I installed the packages install.packages(spdep,dependencies=TRUE) install.packages(maptools,dependencies=TRUE) Then from the shell: R Then load these libraries library(spgrass6) library(spdep) library(maptools) library(boot) library(Matrix) library(lattice) I read the raster index_water_2009-readRAST6(index_water_2009) At this point I stopped. I saw on the mailing list and some tutorials. I can not understand the syntax to use. It only works on vector? Someone can give me some help? Thanks -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/R-Moran-index-getis-e-ord-tp5588509p5588509.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
Re: [GRASS-user] v.net.path extension
/wxpython/wx.path/wx.path.py On Fri, Sep 17, 2010 at 7:01 PM, Robbie Heremans robbie.herem...@gmail.com wrote: Starting from a file with start and end positions (x_start, y_start, x_stop, y_stop) and a vector map of roads, I have to find places on the roads where most people pass, assuming they will use the shortest path or path with the least cost. Is there a command in Grass that can do this (ic v.net.path for each record and summary map of these shortest paths)? Any suggestions/examples? Robbie ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Rashad ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Displaying layers - GRASS 6.4.0RC6 with MSYS
Thanks Micha for putting me on the right track. GRASS_PNGFILE=myimage.png GRASS_PNG_READ=TRUE export GRASS_PNGFILE GRASS_PNG_READ before the d.vect commands or before the for-done-loop does the job. Robbie 2010/9/16 Micha Silver mi...@arava.co.il Robbie Heremans wrote: I use the Windows-version *GRASS* 6.4.0RC6 with *MSYS* and I want to display 2 *layers* using following commands for the Spearfish60 sample location. d.vect streams col=blue d.vect roads col=red The output file map.png only displays the last layer. I think what you need is to set the environment variable GRASS_PNG_READ to TRUE. Then each run of d.vect will start with the existing png file and add to it. How can I display multiple layers with MSYS? Using the wxGUI interface is not an option, I think, since I want to remake the example on http://casoilresource.lawr.ucdavis.edu/drupal/node/340 which uses a loop (map.png only gives the points of the 5th cluster). # there are 5 clusters: show them all, and compute convex hulls for x in `seq 1 5` do v.extract --o in=bclust where=cluster=$x out=bclust_$x v.hull --o in=bclust_$x out=bclust_hull_$x d.vect bclust_hull_$x type=boundary fcol=none width=2 col=white d.vect bclust icon=basic/box fcol=black col=black size=6 done Any solutions? Robbie This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.net.path extension
Starting from a file with start and end positions (x_start, y_start, x_stop, y_stop) and a vector map of roads, I have to find places on the roads where most people pass, assuming they will use the shortest path or path with the least cost. Is there a command in Grass that can do this (ic v.net.path for each record and summary map of these shortest paths)? Any suggestions/examples? Robbie ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Displaying layers - GRASS 6.4.0RC6 with MSYS
I use the Windows-version *GRASS* 6.4.0RC6 with *MSYS* and I want to display 2 *layers* using following commands for the Spearfish60 sample location. d.vect streams col=blue d.vect roads col=red The output file map.png only displays the last layer. How can I display multiple layers with MSYS? Using the wxGUI interface is not an option, I think, since I want to remake the example on http://casoilresource.lawr.ucdavis.edu/drupal/node/340 which uses a loop (map.png only gives the points of the 5th cluster). # there are 5 clusters: show them all, and compute convex hulls for x in `seq 1 5` do v.extract --o in=bclust where=cluster=$x out=bclust_$x v.hull --o in=bclust_$x out=bclust_hull_$x d.vect bclust_hull_$x type=boundary fcol=none width=2 col=white d.vect bclust icon=basic/box fcol=black col=black size=6 done Any solutions? Robbie ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] specify tunnel/bridge in roadmap
Dear, In order to perform (shortest path) network analysis, how can I specify that the crossing of 2 lines in a roadmap is not an intersection but a tunnel/bridge (impossible to turn left or right)? Is v.digit the right way to do this (deleting en changing lines/vertices)? Robbie ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] count vector points (quadrat count)
Dear, How can I count vector points by region resolution (quadrat count) and write the results in a raster map (in order to perform r.mapcalc on these observed frequencies). Is v.to.rast the command I should use (v.to.rast input=myvec output=myras use=??? value=???) Robbie ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] count vector points (quadrat count)
Thanks John. v.out.ascii input=myvec format=point | r.in.xyz input=- output=myras method=n type=FCELL x=1 y=2 z=3 does the job. Robbie 2010/4/8 John Tate john.t...@ntlworld.com Hi, You could try r.in.xyz using 'method=n' which will count points per cell (resolution chosen in g.region) and write a raster map allocating the number count to the cell. John On Thursday 08 April 2010 07:39:12 Robbie Heremans wrote: Dear, How can I count vector points by region resolution (quadrat count) and write the results in a raster map (in order to perform r.mapcalc on these observed frequencies). Is v.to.rast the command I should use (v.to.rast input=myvec output=myras use=??? value=???) Robbie ___ 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