Re: [GRASS-user] v.colors rules
On Mon, Feb 9, 2009 at 8:14 PM, Paolo Craveri pcrav...@gmail.com wrote: 2009/2/9 Markus Neteler nete...@osgeo.org: Ok now it works. Many thanks Perfect. Fixed in all relevant GRASS versions (SVN). Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Creating turbidity raster file
Markus Neteler wrote: You have to use the zcolumn parameter to indicate which column to use. Hope this helps, Markus It worked!! thanks a lot -- View this message in context: http://n2.nabble.com/Creating-turbidity-raster-file-tp2297732p2301479.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
Re: [GRASS-user] Creating turbidity raster file
Dylan Beaudette wrote:
On Monday 09 February 2009, Hamish wrote:
buraq wrote:
I have formatted csv files for every month which is like
latitude;longitude;turbidityValue. So I want to create
raster files from these csvs for using in r.sun.
v.in.ascii - convert .csv to a vector points file
g.region # set up raster grid
v.surf.rast - interpolate vector points to raster surface
Minor correction:
v.surf.rast --- v.surf.rst
I have worried about r.sun's Linke Turbidity factor values in areas with
big changes in elevation. Is turbidity value heavily dependent on
altitude?
It should, as it is based on a measure of optical thickness (air mass):
Extracted from a document --
m = \frac{1}{sin(\alpha) + 0.15(\alpha + 3.885)^{-1.253}} e^{-0.0001184
\times
A}
\noindent
where $\alpha$ is the solar elevation angle at (12:00 pm local time) on
each
day of the year, and $A$ is the weather station elevation in meters above
mean sea level (MSL).
Extracted from a document --
This stuff is also on pages 40-41 of this document:
http://169.237.35.250/~dylan/d_e_beaudette-ms_thesis.pdf
Cheers,
Dylan
Hamish
--
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
I took the turbidity values from soda-is.com. I entered no altitude value
for the latitude and longitude.
--
View this message in context:
http://n2.nabble.com/Creating-turbidity-raster-file-tp2297732p2301484.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] grass startup cmd line
Hi list, as I'm quite a lazy boy, and am currently working on a single grass project, I would appreciate to launch grass from a terminal so that it starts with a given GUI, a given gisbase/location/mapset and a given workspace... This command grass64 -wxpython my_gisbase/my_location/my_mapset works, but is it possible to specify a default workspace to load i.e. grass64 -wxpython my_gisbase/my_location/my_mapset my_workspace.gxw Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass startup cmd line
Hi, 2009/2/10 Vincent Bain b...@toraval.fr: as I'm quite a lazy boy, and am currently working on a single grass project, I would appreciate to launch grass from a terminal so that it starts with a given GUI, a given gisbase/location/mapset and a given workspace... This command grass64 -wxpython my_gisbase/my_location/my_mapset works, but is it possible to specify a default workspace to load i.e. grass64 -wxpython my_gisbase/my_location/my_mapset my_workspace.gxw this is currently possible only from GRASS cmd, e.g. g.gui workspace=my_workspace.gxw (assuming that wxpython is default GUI) 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] grass startup cmd line
Hi, 2009/2/10 Martin Landa landa.mar...@gmail.com: this is currently possible only from GRASS cmd, e.g. g.gui workspace=my_workspace.gxw (assuming that wxpython is default GUI) Please update your contact: grassu...@grass.itc.it - grass-user@lists.osgeo.org Thanks, 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] Re: Calculating standard error of many maps?
Moritz Lennert wrote: Just one clarification: I would like to calculate these descriptive stats for each cell, to obtain variability maps. Rainer On Mon, Feb 9, 2009 at 5:57 PM, Rainer M Krug r.m.k...@gmail.com mailto:r.m.k...@gmail.com wrote: Hi I have 25000 maps, generated by simulation predictions, and I would like to calculate some descriptive stats, like mean, standard deviation, median, quartiles. Is there an easy way of doing this in GRASS? r.series ? You will probably have to use xargs to collate that many map names... Or you can use g.mlist if the names have a common prefix, suffix, etc. This assumes you want a single map that displays the result of all the input maps. E.g. r.series input=`g.mlist rast pattern='prefix_*' sep=,` output=prefix_stddev method=stddev Yes, but if I'm not mistaken, this won't work with 25000 maps as it would create a too long command line. 25000 maps will probably also exceed the resource limit for the number of open descriptors. It's common for user processes to be limited to 1024 descriptors; you need to modify /etc/security/limits.conf to increase the limit. -- 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] grass startup cmd line
Thank you, and I fear it is not possible to launch this command in the child process from the parent (e.g. grass64 -wxpython ... g.gui workspace=my_workspace.gxw) (OK, I purged the former grass-list address) Vincent Le mardi 10 février 2009 à 10:56 +0100, Martin Landa a écrit : Hi, 2009/2/10 Vincent Bain b...@toraval.fr: as I'm quite a lazy boy, and am currently working on a single grass project, I would appreciate to launch grass from a terminal so that it starts with a given GUI, a given gisbase/location/mapset and a given workspace... This command grass64 -wxpython my_gisbase/my_location/my_mapset works, but is it possible to specify a default workspace to load i.e. grass64 -wxpython my_gisbase/my_location/my_mapset my_workspace.gxw this is currently possible only from GRASS cmd, e.g. g.gui workspace=my_workspace.gxw (assuming that wxpython is default GUI) Martin ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] problems starting grass
*Hello, I just installed grass 6.2.3 and every time I start the graphical interface I get the following error message:* GRASS 6.2.3 (RM):~ *** glibc detected *** g.region: double free or corruption (fasttop): 0x08a48ee0 *** === Backtrace: = /lib/tls/i686/cmov/libc.so.6[0xb78e53f4] /lib/tls/i686/cmov/libc.so.6(cfree+0x96)[0xb78e7456] /usr/lib/libstdc++.so.6(_ZdlPv+0x21)[0xb69d00b1] /usr/lib/libstdc++.so.6(_ZNSs4_Rep10_M_destroyERKSaIcE+0x1d)[0xb69ac8bd] /usr/lib/libstdc++.so.6(_ZNSsD1Ev+0x51)[0xb69ae271] /lib/tls/i686/cmov/libc.so.6(__cxa_finalize+0xc1)[0xb78a5101] /usr/local/lib/libgeos-3.1.0.so[0xb6820e64] /usr/local/lib/libgeos-3.1.0.so[0xb68ca070] /lib/ld-linux.so.2[0xb80f0153] /lib/tls/i686/cmov/libc.so.6(exit+0x89)[0xb78a4d69] g.region(main+0x1893)[0x804b5f3] /lib/tls/i686/cmov/libc.so.6(__libc_start_main+0xe5)[0xb788c685] g.region[0x8049b21] === Memory map: 08048000-0804f000 r-xp 08:06 14665504 /usr/lib/grass/bin/g.region 0804f000-0805 r--p 6000 08:06 14665504 /usr/lib/grass/bin/g.region 0805-08051000 rw-p 7000 08:06 14665504 /usr/lib/grass/bin/g.region 08a48000-08a69000 rw-p 08a48000 00:00 0 [heap] b580-b5821000 rw-p b580 00:00 0 b5821000-b590 ---p b5821000 00:00 0 b599d000-b59dc000 r--p 08:06 14042007 /usr/lib/locale/es_ES.utf8/LC_CTYPE b59dc000-b59e rw-p b59dc000 00:00 0 b59e-b59e3000 r-xp 08:06 12542018 /lib/libgpg-error.so.0.3.0 b59e3000-b59e4000 rw-p 2000 08:06 12542018 /lib/libgpg-error.so.0.3.0 b59e4000-b59e6000 r-xp 08:06 12542023 /lib/libkeyutils-1.2.so b59e6000-b59e8000 rw-p 1000 08:06 12542023 /lib/libkeyutils-1.2.so b59e8000-b59e9000 rw-p b59e8000 00:00 0 b59e9000-b59f r-xp 08:06 14010499 /usr/lib/libkrb5support.so.0.1 b59f-b59f1000 r--p 6000 08:06 14010499 /usr/lib/libkrb5support.so.0.1 b59f1000-b59f2000 rw-p 7000 08:06 14010499 /usr/lib/libkrb5support.so.0.1 b59f2000-b5a02000 r-xp 08:06 12559183 /lib/tls/i686/cmov/ libresolv-2.8.90.so b5a02000-b5a03000 r--p f000 08:06 12559183 /lib/tls/i686/cmov/ libresolv-2.8.90.so b5a03000-b5a04000 rw-p 0001 08:06 12559183 /lib/tls/i686/cmov/ libresolv-2.8.90.so b5a04000-b5a06000 rw-p b5a04000 00:00 0 b5a06000-b5a1c000 r-xp 08:06 14010719 /usr/lib/libsasl2.so.2.0.22 b5a1c000-b5a1d000 r--p 00015000 08:06 14010719 /usr/lib/libsasl2.so.2.0.22 b5a1d000-b5a1e000 rw-p 00016000 08:06 14010719 /usr/lib/libsasl2.so.2.0.22 b5a1e000-b5a2a000 r-xp 08:06 14009570 /usr/lib/liblber-2.4.so.2.1.0 b5a2a000-b5a2b000 r--p b000 08:06 14009570 /usr/lib/liblber-2.4.so.2.1.0 b5a2b000-b5a2c000 rw-p c000 08:06 14009570 /usr/lib/liblber-2.4.so.2.1.0 b5a2c000-b5b5e000 r-xp 08:06 14041092 /usr/lib/i686/cmov/libcrypto.so.0.9.8 b5b5e000-b5b5f000 ---p 00132000 08:06 14041092 /usr/lib/i686/cmov/libcrypto.so.0.9.8 b5b5f000-b5b67000 r--p 00132000 08:06 14041092 /usr/lib/i686/cmov/libcrypto.so.0.9.8 b5b67000-b5b74000 rw-p 0013a000 08:06 14041092 /usr/lib/i686/cmov/libcrypto.so.0.9.8 b5b74000-b5b79000 rw-p b5b74000 00:00 0 b5b79000-b5bbb000 r-xp 08:06 14041601 /usr/lib/i686/cmov/libssl.so.0.9.8 b5bbb000-b5bbc000 r--p 00041000 08:06 14041601 /usr/lib/i686/cmov/libssl.so.0.9.8 b5bbc000-b5bbf000 rw-p 00042000 08:06 14041601 /usr/lib/i686/cmov/libssl.so.0.9.8 b5bbf000-b6695000 r--p 08:06 14010458 /usr/lib/libicudata.so.38.1 b6695000-b6696000 r--p 00ad5000 08:06 14010458 /usr/lib/libicudata.so.38.1 b6696000-b67b3000 r-xp 08:06 14010470 /usr/lib/libicuuc.so.38.1 b67b3000-b67bb000 r--p 0011d000 08:06 14010470 /usr/lib/libicuuc.so.38.1 b67bb000-b67bc000 rw-p 00125000 08:06 14010470 /usr/lib/libicuuc.so.38.1 b67bc000-b67be000 rw-p b67bc000 00:00 0 b67be000-b6901000 r-xp 08:06 14108383 /usr/local/lib/ libgeos-3.1.0.so b6901000-b6905000 r--p 00142000 08:06 14108383 /usr/local/lib/ libgeos-3.1.0.so b6905000-b6907000 rw-p 00146000 08:06 14108383 /usr/local/lib/ libgeos-3.1.0.so b6907000-b6908000 rw-p b6907000 00:00 0 b6908000-b6915000 r-xp 08:06 12541973 /lib/libgcc_s.so.1 b6915000-b6916000 r--p c000 08:06 12541973 /lib/libgcc_s.so.1 b6916000-b6917000 rw-p d000 08:06 12541973 /lib/libgcc_s.so.1 b6917000-b69fa000 r-xp 08:06 14008662 /usr/lib/libstdc++.so.6.0.10 b69fa000-b69fe000 r--p 000e3000 08:06 14008662 /usr/lib/libstdc++.so.6.0.10 b69fe000-b69ff000 rw-p 000e7000 08:06 14008662 /usr/lib/libstdc++.so.6.0.10 b69ff000-b6a05000 rw-p b69ff000 00:00 0 b6a05000-b6a1a000 r-xp 08:06 12559174 /lib/tls/i686/cmov/ libnsl-2.8.90.so b6a1a000-b6a1b000 r--p 00014000 08:06 12559174 /lib/tls/i686/cmov/ libnsl-2.8.90.so b6a1b000-b6a1c000 rw-p 00015000 08:06 12559174 /lib/tls/i686/cmov/ libnsl-2.8.90.so b6a1c000-b6a1e000 rw-p b6a1c000 00:00 0 b6a1e000-b6a27000 r-xp 08:06 12559170 /lib/tls/i686/cmov/ libcrypt-2.8.90.so b6a27000-b6a28000 r--p 8000 08:06 12559170 /lib/tls/i686/cmov/
Re: [GRASS-user] Re: Calculating standard error of many maps?
On Tue, Feb 10, 2009 at 12:06 PM, Glynn Clements
gl...@gclements.plus.com wrote:
Moritz Lennert wrote:
Just one clarification: I would like to calculate these descriptive
stats for each cell, to obtain variability maps.
Rainer
On Mon, Feb 9, 2009 at 5:57 PM, Rainer M Krug r.m.k...@gmail.com
mailto:r.m.k...@gmail.com wrote:
Hi
I have 25000 maps, generated by simulation predictions, and I would
like to calculate some descriptive stats, like mean, standard
deviation, median, quartiles.
Is there an easy way of doing this in GRASS?
r.series ?
You will probably have to use xargs to collate that many map names...
Or you can use g.mlist if the names have a common prefix, suffix, etc.
This assumes you want a single map that displays the result of all the
input maps.
E.g.
r.series input=`g.mlist rast pattern='prefix_*' sep=,`
output=prefix_stddev method=stddev
Yes, but if I'm not mistaken, this won't work with 25000 maps as it
would create a too long command line.
25000 maps will probably also exceed the resource limit for the number of
open descriptors. It's common for user processes to be limited to 1024
descriptors; you need to modify /etc/security/limits.conf to increase the
limit.
I asked the same question on the R mailing list, and got the following
promising response:
###
This is how can you can do it with the raster package
# install.packages(raster, repos=http://R-Forge.R-project.org;)
require(raster)
# Try it for a few files first..
n - 10
# create a list (or vector) of file names, e.g. :
fn - list()
for (i in 1:n) { fn[i] - paste('myfile', i, '.tif', sep='') }
# make a RasterStack
s - stack(fn)
r1 - mCalc(s, fun=mean)
r2 - mCalc(s, fun=sd)
#r can be plotted, coerced to sp objects, etc.
plot(r1)
# or saved to file
r1 - setFilename(r1, 'cellmeans.tif')
r1 - writeRaster(r1, format='GTiff')
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.
Rainer
--
Glynn Clements gl...@gclements.plus.com
--
Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation
Biology, UCT), Dipl. Phys. (Germany)
Centre of Excellence for Invasion Biology
Faculty of Science
Natural Sciences Building
Private Bag X1
University of Stellenbosch
Matieland 7602
South Africa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] trace channel from accummulation raster
Hi everybody, Is there a way to extract a channels from the accumulation raster, e.g. accumulation from r.watershed? Somehow r.drain delivers something which is not following the accumulation raster. In theory the accumulation along the drainage should always increase and never decrease. My results show a nice increasing trend, which is interrupted with very low values. Did anybody came across similar problems? thanks for your help! Christoff - Christoff Andermann Geosciences Rennes geomorphology and remote sensing www.rsg.tu-freiberg.de/twiki/bin/view/Main/ChristoffAndermann -- View this message in context: http://n2.nabble.com/trace-channel-from-accummulation-raster-tp2301893p2301893.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] trace channel from accummulation raster
Hi everybody, Is there a way to extract a channels from the accumulation raster, e.g. accumulation from r.watershed? Somehow r.drain delivers something which is not following the accumulation raster. In theory the accumulation along the drainage shoult always increase and never decrease. My results show a nice increasing trend, which is interrupted with very low values. Did anybody came across similar problems? thanks for your help! Christoff - Christoff Andermann Geosciences Rennes geomorphology and remote sensing www.rsg.tu-freiberg.de/twiki/bin/view/Main/ChristoffAndermann -- View this message in context: http://n2.nabble.com/trace-channel-from-accummulation-raster-tp2301896p2301896.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] menu items
Hi list, 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? Kind regards, Koen ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] g.findfile
Hi, Using grass6.4.0svn (revision 35838), the result of g.findfile appears not to run properly for raster query : create a raster map test, run g.findfile element=raster file=test. Looks like the file does not exist : name= mapset= fullname= file= (Tested successfully with a vector map) What am I doing wrong ? I remember the command worked in previous versions. Thanks, Vincent. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] g.findfile
Hi, 2009/2/10 Vincent Bain b...@toraval.fr: Using grass6.4.0svn (revision 35838), the result of g.findfile appears not to run properly for raster query : create a raster map test, run g.findfile element=raster file=test. Looks like the file does not exist : name= mapset= fullname= file= (Tested successfully with a vector map) What am I doing wrong ? I remember the command worked in previous versions. use element=cell to find raster map 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] g.findfile
Thank you Martin ! (for my info, is it a recent change or was the previous keyword really 'raster' ?) Vincent Le mardi 10 février 2009 à 14:39 +0100, Martin Landa a écrit : Hi, 2009/2/10 Vincent Bain b...@toraval.fr: Using grass6.4.0svn (revision 35838), the result of g.findfile appears not to run properly for raster query : create a raster map test, run g.findfile element=raster file=test. Looks like the file does not exist : name= mapset= fullname= file= (Tested successfully with a vector map) What am I doing wrong ? I remember the command worked in previous versions. use element=cell to find raster map Martin ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] trace channel from accummulation raster
Christoff wrote: Hi everybody, Is there a way to extract a channels from the accumulation raster, e.g. accumulation from r.watershed? You can either set the threshold option in r.watershed and get stream segments or use flow accumulation output and do something like r.mapcalc channels = if(flow_acc threshold, flow_acc, null()), threshold would here be the desired number of cells draining through a given cell when using r.watershed flow accumulation or size of the area draining through a given cell for r.terraflow flow accumulation. See also Further processing of output layers in the manual of r.watershed. Somehow r.drain delivers something which is not following the accumulation raster. In theory the accumulation along the drainage shoult always increase and never decrease. My results show a nice increasing trend, which is interrupted with very low values. As I understand, r.drain works on elevation and uses a different algorithm, so the results are likely different, but I don't know r.drain that well. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] g.findfile
Hi, 2009/2/10 Vincent Bain b...@toraval.fr: Thank you Martin ! (for my info, is it a recent change or was the previous keyword really 'raster' ?) no, 'cell' is historical artefact. One of topics for GRASS7 - clean up list of GIS elements used in GRASS 5/6 [1] Martin [1] http://trac.osgeo.org/grass/browser/grass/branches/develbranch_6/lib/gis/element_list -- 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] wxpython: vdigit crash
Hi, 2009/2/3 Craig Leat craig.l...@gmail.com: Hi Martin Martin Landa: right, here it works;-) Try to enable also vdigit debuging export GRASS_WX_DEBUG=3 g.gisenv set=DEBUG=3 last 50 lines of debug output should be enough... 50 lines of debug output attached. This snippet looks strange (the last two values are very LARGE): D3/3: Box(N,S,E,W,T,B): -2.891490e+01, -3.014472e+01, 3.091949e+01, 2.948216e+01, 1.797693e+308, -1.797693e+308 this is OK, T/B is vertical extent. From your debug, vector digitizer crashes after http://trac.osgeo.org/grass/browser/grass/branches/develbranch_6/gui/wxpython/vdigit/driver_draw.cpp#L52 Try to disable drawing areas and but some debug messages around http://trac.osgeo.org/grass/browser/grass/branches/develbranch_6/gui/wxpython/vdigit/driver_draw.cpp#L166 ? 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
[GRASS-user] updating vector map using another vector map as criteria
Is there a way I can update the attributes of one vector map, given another vector map(vm)? The selection vector map would be subset of original vm. As an example, using output from v.net.path to update a column in original vm, say just updating some arbitrary integer field to signify that there was a match between original and other vm. Thanks. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] use commands window
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. Am I missing something? Thank you very much for your help. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Tips on how to delete areas with no categories?
Hello Grass-Users, After using r.to.vect, I have areas with no categories. These areas appear to be holes. This is what the map look like with info from v.what: http://bobdebm.googlepages.com/example_area.png/example_area-full;init:.png How do I remove these areas? Thanks in advance. Bob Bob Moskovitz Research Analyst I Seismic Hazard Zonation Project California Geological Survey http://www.conservation.ca.gov/cgs/shzp?xml:namespace prefix = o ns = urn:schemas-microsoft-com:office:office / 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. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] DAFF/CSIRO soils map
Hamish, You mentioned this dataset in a previous discussion ... Jarl wrote: The data is an Australian soils atlas including hydraulic properties in the dbf. this? http://www.daff.gov.au/brs/data-tools/daas-download Being curious I took a look ... the data uncompresses as csiro2mgvfc.e00 but r.in.arc fails on this: Illegal line in header EXP 1 /REFERENCE/NATIONAL/ARCINFO/SOIL/COVERS/CSIRO2MG.E00 yllcorner field missing from header xllcorner field missing from header nrows field missing from header ncols field missing from header cellsize field missing from header nodata_value field missing from header Can't get cell header What would cause this? Richard C ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] DAFF/CSIRO soils map
Moskovitz, Bob wrote: Richard, I believe you need to use v.in.e00. -Bob /blush Sheepishly I thank you! RC -Original Message- From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-boun...@lists.osgeo.org]on Behalf Of Richard Chirgwin Sent: Tuesday, February 10, 2009 12:36 PM To: grass-user@lists.osgeo.org Subject: [GRASS-user] DAFF/CSIRO soils map Hamish, You mentioned this dataset in a previous discussion ... Jarl wrote: The data is an Australian soils atlas including hydraulic properties in the dbf. this? http://www.daff.gov.au/brs/data-tools/daas-download Being curious I took a look ... the data uncompresses as csiro2mgvfc.e00 but r.in.arc fails on this: Illegal line in header EXP 1 /REFERENCE/NATIONAL/ARCINFO/SOIL/COVERS/CSIRO2MG.E00 yllcorner field missing from header xllcorner field missing from header nrows field missing from header ncols field missing from header cellsize field missing from header nodata_value field missing from header Can't get cell header What would cause this? Richard C ___ 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] OGRS2009 Call For Research Papers, Showcases and Labs - Deadline extension
-- OGRS2009 Call For Research Papers, Showcases and Labs - Deadline extension First International Opensource Geospatial Research Symposium 2009 Dates place: July 8th to 10th, 2009, Ecole Centrale de Nantes (France) Website: http://www.ogrs2009.org -- Dear colleagues, (our apologies for cross postings) From July 8th to 10th, 2009, during the tenth edition of Libre Software Meeting, the Research Institute on Urban Sciences and Techniques (IRSTV) in partnership with the University Of Applied Sciences Western Switzerland, Prefecture and Council of Region Pays de La Loire will organize a Research and Innovation Symposium about free and open source geospatial methods and technologies. Main goals of this symposium are : * to build a panel of new scientific works built on open source models or using open source tools * to build and discuss a scientific framework about open source technology usage (benefits and limitations) * to provide a platform to network and develop ideas for future collaborative works between the academic research world and the everyday operational world (companies, local authorities ...) The symposium will integrate scientific conferences (research papers), companies and authority testimonies (showcases) and training sessions (labs). Thus, we invite to submit contributions which can be : * a full paper or extended abstract for the research session, * an abstract for the showcase session, * an abstract for the lab session. Contributions should address the following overall topics but are not limited to : * Cartography and advanced styling ; * Cultural heritage ; * Earth Observation ; * From Webmapping to WebGIS ; * Integrated risk and disaster management for society ; * GIS and SDI interoperability (data, knowledge, language, ...) ; * GIS and SDI for participatory learning and action ; * GIS concept and theory (semantics and ontology) ; * Landscape and ecology ; * Spatial Analysis and Integration ; * Spatial and Spatio-temporal data and analysis ; * Sustainable Development and Governance ; * Urban and Environmental Planning ; * Visual Languages and Querying ; * Web and desktop convergence in the geospatial field. Keynote speakers : * Câmara Gilberto, Brazil's National Institute for Space Research * Neteler Markus, Fondazione Mach - Centro di Ecologia Alpina, Italia * Steiniger Stefan, University of Calgary, Alberta, Canada See committee members on OGRS2009 website : http://committee.ogrs2009.org Important : submission deadlines has been extended due to several requests from authors Research paper submission (see http://research.ogrs2009.org) : - Full paper : * Extended submission deadline : April, 3th 2009 * Submit a 500-1000 words abstract and a 12 pages full paper * Notification to authors: May, 11th 2009 - Extended abstract : * Extended Submission deadline : April, 3th 2009 * Submit a 500-1000 words abstract * Notification to authors: April, 29th 2009 Lab abstract submission (see http://lab.ogrs2009.org) : * Extended submission deadline : April, 3th 2009 * Submit a 500-1000 words abstract * Notification to authors: April, 29th 2009 Showcase abstract submission (see http://showcase.ogrs2009.org) : * Extended submission deadline : April, 3 th 2009 * Submit a 500-1000 words abstract * Notification to authors: April, 29th 2009 Contacts : * for research submission : resea...@ogrs2009.org * for lab submission : l...@ogrs2009.org * for other questions : i...@ogrs2009.org We would appreciate if you could kindly distribute this call to other interested parties of your acquaintance. Best regards, OGRS2009 program committee ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Thiessen Polygons
Hi folks, 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. Thanks for any help you can give. Kurt ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Thiessen Polygons
Try v.surf.idw Cheers Harald Kurt Springs wrote: Hi folks, 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. Thanks for any help you can give. Kurt ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- --- Harald Hofmann PhD Candidate - Hydrogeology School of Geosciences Bldg 28 Monash University Clayton Campus Wellington Road VIC 3800 Australia harald.hofm...@sci.monash.edu.au Office: +61 3 9905 5786 Fax:+61 3 9905 4903 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Thiessen Polygons
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 Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Tips on how to delete areas with no categories?
Bob Moskovitz wrote: After using r.to.vect, I have areas with no categories. These areas appear to be holes. How do I remove these areas? Richard Chirgwin wrote: What if you use v.category to attach cats to the empty areas, and then use v.edit to delete the appropriate category range? (or v.dissolve instead of v.edit) another old trick: 'v.extract cat=1-9' will only extract actual data, ie will skip holes. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Thiessen Polygons
That did the trick. Thanks Hamish Kurt On Feb 10, 2009, at 6:47 PM, Hamish wrote: 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 Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass startup cmd line
Vincent Bain wrote: I would appreciate to launch grass from a terminal so that it starts with a given GUI, a given gisbase/location/mapset and a given workspace... Martin Landa: this is currently possible only from GRASS cmd, e.g. g.gui workspace=my_workspace.gxw (assuming that wxpython is default GUI) Vincent: I fear it is not possible to launch this command in the child process from the parent (e.g. grass64 -wxpython ... g.gui workspace=my_workspace.gxw) maybe GRASS_BATCH JOB could do it: http://grass.osgeo.org/wiki/GRASS_and_Shell#GRASS_Batch_jobs Also setting up an alias for the long command line in .alias is handy: alias spear='grass64 /home/user/grassdata/spearfish60/user1' or a little script in ~/bin/ to setup run the batch job. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Problemas con winGRASS 6.4.0RC3 del instalador OSGeo4W
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! http://spaces.live.com/spacesapi.aspx?wx_action=createwx_url=/friends.aspxmkt=en-us___ 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?
On Wed, Feb 11, 2009 at 4:06 AM, Hamish hamis...@yahoo.com wrote: Rainer Krug wrote: I have 25000 maps, generated by simulation predictions, and I would like to calculate some descriptive stats, like mean, standard deviation, median, quartiles. 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() Thanks a lot Hamish - I will look into that and see what can do. In the meantime, it works in R, as described above. Cheers Rainer Hamish -- Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation Biology, UCT), Dipl. Phys. (Germany) Centre of Excellence for Invasion Biology Faculty of Science Natural Sciences Building Private Bag X1 University of Stellenbosch Matieland 7602 South Africa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Creating turbidity raster file
Hamish:
I have worried about r.sun's Linke Turbidity factor values in areas with
big changes in elevation. Is turbidity value heavily dependent on
altitude?
Jose:
Linke turbidity is certainly height-dependent given the higher optical
path length at lower elevations. A simple approximated pressure correction
might be applied if the turbidity at a given altitude z' is known:
TL(z) = TL(z') exp( -(z-z') / 8435.2)
Dylan:
It should, as it is based on a measure of optical thickness (air mass):
m = \frac{1}{sin(\alpha) + 0.15(\alpha + 3.885)^{-1.253}} e^{-0.0001184
\timesA}
(yay LyX, ctrl-m, paste selection, view PostScript)
buraq wrote:
I took the turbidity values from soda-is.com. I entered no altitude
value for the latitude and longitude.
I did the same, I attempted to match their model grid (approx 1x1-deg IIRC)
for an array of lat/lon extractions over my study area, then v.surf.rst to
produce a Linke layer for r.sun. I can't recall off the top of my head if
we entered an altitude or not, I think not.
I work in fjords with 1000m vertical drops (including one of the top
10 tallest waterfalls in the world). With little idea of the elevation
values used for the SODA data, I just guess that it will be very rough
and subject to a somewhat arbitrary sampling error.
I assume I'll have to use one of the above formulas with r.mapcalc to
create the Linke layer from the SODA data a local DEM for delta-z. But
that just makes me wonder if turbidity really belongs in a 3D voxel grid,
not a 2D coverage map? i.e. I can derive a value for Linke at the
ground-surface for each raster cell in the DEM easily enough using the
above formulas, but isn't the important Linke value(s) what the beam
encounters on it's path through the atmosphere high above the ground,
not just the value at ground level at the beam's terminus?
If it helps, I do have quite a bit of PAR light meter data we collected
at sea level throughout the region over full years; and do know the air
is nearly as clear as it gets. What instrumentation could I add to our
met stations to get a better record? Deploying light loggers on mountain
peaks for the summer months may be an option (granite is beautiful stuff
to climb :).
hope that makes sense,
Hamish
ps- If anyone is interested, at one point I wrote a little C program
that fit a sine curve to the monthly data and gave you a per-day value
from per-month data, avoiding big jumps at the month changes. The same
could in theory be adapted to work with r.mapcalc to get daily maps from
the monthly raster layers, at considerable computational cost. (but you
just need to run it once(ie x365) in a dedicated mapset)
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
