Re: [GRASS-user] empty map window/map rendering doesn't work (wxgui)
Hi, 2009/2/8 kries kr...@gmx.at: I have a very annoying problem. Maps are not rendered in wxgui anymore after recompiling grass7 (svn). Compiled cleanly (distclean/clean) ./configure --with-cairo --with-mysql --with-mysql-includes=/usr/include/mysql --with-freetype--with-freetype-includes=/usr/include/freetype2/ --with-fretype-lib=/usr/lib/ --with-python --with-sqlite --with-odcb --with-wxwidgets=/usr/bin/wx-config --with-readline --with-opengl --with-cxx --with-gdal --with-opengl --with-wxwidgets-include=/usr/include/wx-2.8/ gentoo kde4.2 python 2.5 wxwidgets 2.6 / 2.8 gdal-1.5.3 gcc-4.3.3 / 4.2.4 any ideas? No problem here. No exact error message? 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] empty map window/map rendering doesn't work (wxgui)
command output: Traceback (most recent call last): File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 443, in OnPaint self.pdc.DrawToDCClipped(dc, rgn) File /usr/local/grass-7.0.svn/etc/wxpython/vdigit/grass7_ wxvdigit.py, line 386, in DrawToDCClipped def DrawToDCClipped(self, *args): return _grass7_wxvdigit.PseudoDC_DrawToDCClipped(self, *args) TypeError : in method 'PseudoDC_DrawToDCClipped', argument 2 of type 'void *' Traceback (most recent call last): File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 535, in OnIdle self.UpdateMap(render=True) File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 660, in UpdateMap self.Draw(self.pdc, self.img, drawid=id) File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 263, in Draw pdc.SetBackground(bg) File /usr/local/grass-7.0.svn/etc/wxpython/vdigit/grass7_ wxvdigit.py, line 378, in SetBackground def SetBackground(self, *args): return _grass7_wxvdigit.PseudoDC_SetBackground(self, *args) TypeError : in method 'PseudoDC_SetBackground', argument 2 of type 'wxBrush const ' Traceback (most recent call last): File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp.py, line 560, in OnRender self.MapWindow.UpdateMap(render=True) File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 660, in UpdateMap self.Draw(self.pdc, self.img, drawid=id) File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 263, in Draw pdc.SetBackground(bg) File /usr/local/grass-7.0.svn/etc/wxpython/vdigit/grass7_ wxvdigit.py, line 378, in SetBackground def SetBackground(self, *args): return _grass7_wxvdigit.PseudoDC_SetBackground(self, *args) TypeError : in method 'PseudoDC_SetBackground', argument 2 of type 'wxBrush const ' No exact error message? 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 -- View this message in context: http://n2.nabble.com/empty-map-window-map-rendering-doesn%27t-work-%28wxgui%29-tp2294077p2296504.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] empty map window/map rendering doesn't work (wxgui)
Maris Nartiss wrote: Yes. Don't use GRASS 7 unless You are doing some GRASS 7 development. Probably we need to add HUGE warning for 7 - Not for general consumption! That's an exaggeration. It's less stable than 6.x, but most things should work most of the time. Anything that doesn't work will only get fixed as people try it and report bugs. Given the relatively small number of developers, GRASS' continued existence is heavily dependent upon users' willingness to test development versions. -- 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] v.colors rules
On Mon, Feb 9, 2009 at 12:12 PM, Paolo Craveri pcrav...@gmail.com wrote: Hi to all How v.colors with option rules works ? Paolo, in line 160 of the script, can you change if [ -n $OPTNAME ] ; then to if [ -z $OPTNAME ] ; then and report back? For me it works like this (maybe it breaks something else?!). Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.colors rules
On Mon, 2009-02-09 at 12:12 +0100, Paolo Craveri wrote: Hi to all How v.colors with option rules works ? I tried: cat color.txt 1 red 2 green 3 blue 4 orange 5 purple v.db.select mypoints cat 5 1 4 3 2 v.colors map=mypoints column=cat rules=color.txt ERROR: Pick one of color, rules, or raster options v.colors map=mypoints color=rules column=cat rules=color.txt ERROR: value rules out of range for parameter color v.colors map=mypoints color=rules.txt ERROR: value rules.txt out of range for parameter color Thanks -- Ciao Paolo. It doesn't work for me either! Nikos = # query samples_unburned6 for unique values in the 2nd column v.db.select samples_unburned6 | cut -d| -f2 | sort -u 1 3 332 5 value # v.colors with random (or other) colors works just fine v.colors samples_unburned6 column=value color=random # v.colors rules= or raster fails # first the *rules* cat samples_color.rules -333 255:0:0 -331 0:255:0 1 255:0:0 2 255:255:0 3 0:255:0 5 0:0:255 332 200:200:200 334 0:0:0 # try to use rules= v.colors samples_unburned6 column=value rules=samples_color.rules ERROR: Pick one of color, rules, or raster options # list the colr file for some raster map cat training_6_prefire % -333 334 -333:255:0:0 -331:0:255:0 -331:0:255:0 1:255:0:0 1:255:0:0 2:255:255:0 2:255:255:0 3:0:255:0 3:0:255:0 5:0:0:255 5:0:0:255 332:200 332:200 334:0 # try to use raster= v.colors samples_unburned6 column=value raster=training_6_prefire ERROR: Pick one of color, rules, or raster options == ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] empty map window/map rendering doesn't work (wxgui)
Hi 2009/2/9 kries kr...@gmx.at: [...] pdc.SetBackground(bg) File /usr/local/grass-7.0.svn/etc/wxpython/vdigit/grass7_ wxvdigit.py, line 378, in SetBackground def SetBackground(self, *args): return _grass7_wxvdigit.PseudoDC_SetBackground(self, *args) TypeError : in method 'PseudoDC_SetBackground', argument 2 of type 'wxBrush const ' r35819 (trunk) and r35820 (devbr6) - hopefully helps(?) 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] v.colors rules
On Mon, Feb 9, 2009 at 1:44 PM, Markus Neteler nete...@osgeo.org wrote: On Mon, Feb 9, 2009 at 12:12 PM, Paolo Craveri pcrav...@gmail.com wrote: Hi to all How v.colors with option rules works ? Paolo, in line 160 of the script, can you change forget that. The problem is that the color bcyr is predefined. Try to deleted line 66: #% answer: bcyr and report back. Then it seems to work. I tested color= and rules=. Also added a new example to the manual page. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.colors rules
2009/2/9 Markus Neteler nete...@osgeo.org: Try to deleted line 66: #% answer: bcyr ok, done, but: v.colors mypoints rules=color.txt column=cat ERROR: Pick one of color, rules, or raster options I commented lines from 163 to 166 (test over options) and works: v.colors mypoints rules=color.txt column=cat v.db.select mypoints cat|GRASSRGB 5|127:0:255 1|255:0:0 4|255:127:0 3|0:0:255 2|0:255:0 ciao -- -- Paolo C. Lat. 44° 39' 11.08'' N Long. 7° 23' 25.26'' E ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] empty map window/map rendering doesn't work (wxgui)
Hi, thank you, unfortunately it still doesn't work. kries Traceback (most recent call last): File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 443, in OnPaint self.pdc.DrawToDCClipped(dc, rgn) File /usr/local/grass-7.0.svn/etc/wxpython/vdigit/grass7_ wxvdigit.py, line 386, in DrawToDCClipped def DrawToDCClipped(self, *args): return _grass7_wxvdigit.PseudoDC_DrawToDCClipped(self, *args) NotImplementedError : Wrong number of arguments for overloaded function 'PseudoDC_DrawToDCClipped'. Possible C/C++ prototypes are: DrawToDCClipped(gwxPseudoDC *,wxBufferedPaintDC *,wxRect const ) DrawToDCClipped(gwxPseudoDC *,wxGCDC *,wxRect const ) Traceback (most recent call last): File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 535, in OnIdle self.UpdateMap(render=True) File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 660, in UpdateMap self.Draw(self.pdc, self.img, drawid=id) File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 263, in Draw pdc.SetBackground(bg) File /usr/local/grass-7.0.svn/etc/wxpython/vdigit/grass7_ wxvdigit.py, line 378, in SetBackground def SetBackground(self, *args): return _grass7_wxvdigit.PseudoDC_SetBackground(self, *args) TypeError : in method 'PseudoDC_SetBackground', argument 2 of type 'wxBrush const ' Traceback (most recent call last): File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 443, in OnPaint self.pdc.DrawToDCClipped(dc, rgn) File /usr/local/grass-7.0.svn/etc/wxpython/vdigit/grass7_ wxvdigit.py, line 386, in DrawToDCClipped def DrawToDCClipped(self, *args): return _grass7_wxvdigit.PseudoDC_DrawToDCClipped(self, *args) NotImplementedError : Wrong number of arguments for overloaded function 'PseudoDC_DrawToDCClipped'. Possible C/C++ prototypes are: DrawToDCClipped(gwxPseudoDC *,wxBufferedPaintDC *,wxRect const ) DrawToDCClipped(gwxPseudoDC *,wxGCDC *,wxRect const ) Traceback (most recent call last): File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 443, in OnPaint self.pdc.DrawToDCClipped(dc, rgn) File /usr/local/grass-7.0.svn/etc/wxpython/vdigit/grass7_ wxvdigit.py, line 386, in DrawToDCClipped def DrawToDCClipped(self, *args): return _grass7_wxvdigit.PseudoDC_DrawToDCClipped(self, *args) NotImplementedError : Wrong number of arguments for overloaded function 'PseudoDC_DrawToDCClipped'. Possible C/C++ prototypes are: DrawToDCClipped(gwxPseudoDC *,wxBufferedPaintDC *,wxRect const ) DrawToDCClipped(gwxPseudoDC *,wxGCDC *,wxRect const ) Traceback (most recent call last): File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 443, in OnPaint self.pdc.DrawToDCClipped(dc, rgn) File /usr/local/grass-7.0.svn/etc/wxpython/vdigit/grass7_ wxvdigit.py, line 386, in DrawToDCClipped def DrawToDCClipped(self, *args): return _grass7_wxvdigit.PseudoDC_DrawToDCClipped(self, *args) NotImplementedError : Wrong number of arguments for overloaded function 'PseudoDC_DrawToDCClipped'. Possible C/C++ prototypes are: DrawToDCClipped(gwxPseudoDC *,wxBufferedPaintDC *,wxRect const ) DrawToDCClipped(gwxPseudoDC *,wxGCDC *,wxRect const ) Traceback (most recent call last): File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp.py, line 560, in OnRender self.MapWindow.UpdateMap(render=True) File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 660, in UpdateMap self.Draw(self.pdc, self.img, drawid=id) File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 263, in Draw pdc.SetBackground(bg) File /usr/local/grass-7.0.svn/etc/wxpython/vdigit/grass7_ wxvdigit.py, line 378, in SetBackground def SetBackground(self, *args): return _grass7_wxvdigit.PseudoDC_SetBackground(self, *args) TypeError : in method 'PseudoDC_SetBackground', argument 2 of type 'wxBrush const ' Traceback (most recent call last): File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 443, in OnPaint self.pdc.DrawToDCClipped(dc, rgn) File /usr/local/grass-7.0.svn/etc/wxpython/vdigit/grass7_ wxvdigit.py, line 386, in DrawToDCClipped def DrawToDCClipped(self, *args): return _grass7_wxvdigit.PseudoDC_DrawToDCClipped(self, *args) NotImplementedError : Wrong number of arguments for overloaded function 'PseudoDC_DrawToDCClipped'. Possible C/C++ prototypes are: DrawToDCClipped(gwxPseudoDC *,wxBufferedPaintDC *,wxRect const ) DrawToDCClipped(gwxPseudoDC *,wxGCDC *,wxRect const ) Traceback (most recent call last): File /usr/local/grass-7.0.svn/etc/wxpython/gui_modules/ma pdisp_window.py, line 443, in OnPaint self.pdc.DrawToDCClipped(dc, rgn) File /usr/local/grass-7.0.svn/etc/wxpython/vdigit/grass7_ wxvdigit.py, line 386, in
Re: [GRASS-user] empty map window/map rendering doesn't work (wxgui)
oh, ok :) yes, /usr/bin/wx-config --version-full 2.8.9.1 kries Martin Landa wrote: Hi, 2009/2/8 kries kr...@gmx.at: [...] --with-wxwidgets=/usr/bin/wx-config --with-readline --with-opengl --with-cxx --with-gdal --with-opengl --with-wxwidgets-include=/usr/include/wx-2.8/ BTW, `--with-wxwidgets-include` doesn't exists - it's just ignored during configuration. Does /usr/bin/wx-config point to 2.8? 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 -- View this message in context: http://n2.nabble.com/empty-map-window-map-rendering-doesn%27t-work-%28wxgui%29-tp2294077p2297531.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] Creating turbidity raster file
Hi all, I am new in grass. Maybe that can be easy for you.. 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. What I have done so far, I created MIFs in mapInfo then imported grass using v.in.ogr. Then I used v.surf.rst to make them raster. But I think it's failed. I got that http://n2.nabble.com/file/n2297732/MapDisplay.png MapDisplay.png How can I do this right? Thanks for your help -- View this message in context: http://n2.nabble.com/Creating-turbidity-raster-file-tp2297732p2297732.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
On Mon, Feb 9, 2009 at 4:27 PM, buraq zak...@gmail.com wrote: Hi all, I am new in grass. Maybe that can be easy for you.. 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. What I have done so far, I created MIFs in mapInfo then imported grass using v.in.ogr. Then I used v.surf.rst to make them raster. But I think it's failed. I got that http://n2.nabble.com/file/n2297732/MapDisplay.png MapDisplay.png How can I do this right? You have to use the zcolumn parameter to indicate which column to use. Hope this helps, Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Calculating standard error of many maps?
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 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? Thanks, Rainer -- 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 -- 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] Re: Calculating standard error of many maps?
Hi Martin On Mon, Feb 9, 2009 at 6:12 PM, Martin Wegmann martin.wegm...@uni-wuerzburg.de wrote: hello Rainer, just look at r.series. That is exactly what I was looking for. Thanks a lot Rainer If the naming of the files is similar you can try something like: r.series input=`g.mlist pattern='insitu_data.*' sep=,` output=insitu_data.stddev method=stddev regards, Martin On Montag, 9. Februar 2009 16:59:33 Rainer M Krug 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 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? Thanks, Rainer -- 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 -- ** University of Wuerzburg Institute of Geography Department of Remote Sensing Remote Sensing and Biodiversity Unit Am Hubland 97074 Wuerzburg, Germany @ German Aerospace Center (DLR) German Remote Sensing Data Center (DFD) Phone: +49-(0)931-888-4797 Fax:+49-(0)931-888-4961 Email: martin.wegm...@uni-wuerzburg.de url: http://www.remote-sensing.uni-wuerzburg.de ** -- 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] empty map window/map rendering doesn't work (wxgui)
hi, swig-1.3.38 kries Martin Landa wrote: Hi, 2009/2/9 kries kr...@gmx.at: yes, /usr/bin/wx-config --version-full 2.8.9.1 which version of swig do you use? M. -- 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 -- View this message in context: http://n2.nabble.com/empty-map-window-map-rendering-doesn%27t-work-%28wxgui%29-tp2294077p2298010.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] Re: Calculating standard error of many maps?
On 09/02/09 16:59, Rainer M Krug 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 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... Moritz ___ 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 Feb 9, 2009, at 8:15 AM, Moritz Lennert wrote: On 09/02/09 16:59, Rainer M Krug 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 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 John___ 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?
hello Rainer, just look at r.series. If the naming of the files is similar you can try something like: r.series input=`g.mlist pattern='insitu_data.*' sep=,` output=insitu_data.stddev method=stddev regards, Martin On Montag, 9. Februar 2009 16:59:33 Rainer M Krug 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 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? Thanks, Rainer -- 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 -- ** University of Wuerzburg Institute of Geography Department of Remote Sensing Remote Sensing and Biodiversity Unit Am Hubland 97074 Wuerzburg, Germany @ German Aerospace Center (DLR) German Remote Sensing Data Center (DFD) Phone: +49-(0)931-888-4797 Fax:+49-(0)931-888-4961 Email: martin.wegm...@uni-wuerzburg.de url: http://www.remote-sensing.uni-wuerzburg.de ** ___ 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 09/02/09 17:19, John C. Tull wrote: On Feb 9, 2009, at 8:15 AM, Moritz Lennert wrote: On 09/02/09 16:59, Rainer M Krug 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. But please correct me if I'm wrong. Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.colors rules
2009/2/9 Markus Neteler nete...@osgeo.org: Are you sure that you ran make and, if needed, make install afterwards? Markus Ok now it works. Many thanks -- -- Paolo C. Lat. 44° 39' 11.08'' N Long. 7° 23' 25.26'' E ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Cut/Fill Analysis
Grass Community, Anybody in the community have experience using grass for creating cut-fill maps a la Arc/Info (http://webhelp.esri.com/arcgisdesktop/9.2/index.cfm?TopicName=Cut/Fill) I have xyz coordinates for my current site. I'd like to input a desired slope, and have the GIS output the areas where earth needs to be moved and areas where earth needs to be added, with perhaps volume estimates as well. Thanks in advance, -- Frank Aragona ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Cut/Fill Analysis
I've used GRASS to run preliminary quantities for a landfill design - It works quite well. Create a DEM that represents your proposed grades and then using mapcalc you can figure the cuts and fills relative to the existing topography. Then you can easily compute the volumes. If you want to balance the cut/fill, remember that there is a bulking factor on excavated material (say, 1.3X) and compaction of the fill (if any) needs to be considered. Rasterizing the proposed cut and fill surfaces can be a tedious process if there's a lot of detail. I would say that if you have access to a CAD package that does earthwork for a living it'll be a lot easier - and quicker. But at a pinch, this is a very useful application of GRASS's many talents. Stu On Feb 9, 2009, at 3:06 PM, Frank Aragona wrote: Grass Community, Anybody in the community have experience using grass for creating cut-fill maps a la Arc/Info (http://webhelp.esri.com/arcgisdesktop/9.2/index.cfm?TopicName=Cut/Fill ) I have xyz coordinates for my current site. I'd like to input a desired slope, and have the GIS output the areas where earth needs to be moved and areas where earth needs to be added, with perhaps volume estimates as well. Thanks in advance, -- Frank Aragona ___ 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] Importing ArcASCII raster and associated dbf table?
thanks Hamish, The data is an Australian soils atlas including hydraulic properties in the dbf. From the accompanying documentation Spatial data are available as ArcASCII format raster files. The associated database is a binary file in dBase or dbf format. The projection information is in ArcGIS format. This data set consists of nine sets of data files each containing the following: aas_*.asc An ArcASCII file containing a raster of the soil landscape index. aas_*.dbf A dBase file of soil properties linked to the soil landscape index. aas_*.prj A file containing projection information in the ArcGIS format. I don't know what arcviews capabilities are, so I'm not sure how it would handle the data. However it seems like I should be able to link the database in grass to the vector map I created. The .dbf has the following columns the first column is (int) second (varchar(12) and the rest are all integers or floating pt : VALUE , MAP_UNIT , DATASUPT, SOLDEPTH, SOLPAWHC, A_THICK, A_KSAT, A_KS_ERR, A_SAT, A_FCP, A_WP, A_PAWHC, A_WR_UNR, B_PERCNT, B_THICK, B_KSAT, B_KS_ERR, B_SAT, B_FCP, B_WP, B_PAWHC, B_WR_UNR () The dbf column Value corresponds to the CAT number of the vector map I created. The DBF contains additional value codes since not all the soils are present in the area I'm looking at. So cat 3 matches to value 3 etc. V.db.connect may be the right tool however I don't know how to drive it. All the tutorials I found were for old grass versions and seemed to assume I knew how to handle databases. I'm now using the native windows grass 6.3 Kindest regards Jarl -Original Message- From: Hamish [mailto:hamis...@yahoo.com] Sent: Friday, 6 February 2009 9:37 AM To: grass-user@lists.osgeo.org; Devereux, Jarl M Subject: Re: [GRASS-user] Importing ArcASCII raster and associated dbf table? Jarl Devereux wrote: I'm new to grass and I'm trying to import an ArcASCII raster file with an associated dbf table into a cygwin grass 6.2. I have been using r.in.gdal to import the raster without issue (r.in.arc might work too) however handling the attribute data is still a mystery to me. Can I connect a raster to a dbf table or add the attributes as layers? what is the form of the fields in the database? If simply category:label you can use the rules= option of the r.category module in GRASS 6.3+ to attach them. For GRASS 6.2 you would have to create the cats/ file by hand. (it's just a text file) see http://grass.osgeo.org/grass64/manuals/html64_user/r.category.html There is the possibility that you could pipe the dbf into r.category with the dbview utility. The db.* commands seemed to be only vector based right, so I converted the raster to an area vector using the raster values as the cat. I am not sure of the steps to copy or attach the arc dbf table. v.db.connect? my cat values are not consecutive that's ok and the .dbf table that accompanied the data doesn't have a cat column it has a consecutive value column instead. so is the Arc map categorical or floating point real numbers? what does it represent? is the DBF file like 0.0-500m: Lowlands 500.0-1000.0m Midlands 1000.0-3000.0m Highlands or..? How do I go about adding this table to the vector so it correctly matches up value and cat? need more info. Hamish DISCLAIMER The information contained in the above e-mail message or messages (which includes any attachments) is confidential and may be legally privileged. It is intended only for the use of the person or entity to which it is addressed. If you are not the addressee any form of disclosure, copying, modification, distribution or any action taken or omitted in reliance on the information is unauthorised. Opinions contained in the message(s) do not necessarily reflect the opinions of the Queensland Government and its authorities. If you received this communication in error, please notify the sender immediately and delete it from your computer system network. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Creating turbidity raster file
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 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? Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
RE: [GRASS-user] Importing ArcASCII raster and associated dbf table?
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 if not, is there a link to download the dataset? This data set consists of nine sets of data files each containing the following: aas_*.asc An ArcASCII file containing a raster of the soil landscape index. aas_*.dbf A dBase file of soil properties linked to the soil landscape index. aas_*.prj A file containing projection information in the ArcGIS format. It is an interesting problem to me. Are you interested in one of the DBF fields, or several/all of them? As a quick guess, I think the answer will be to create 1 map per DB column by some simple script. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Creating turbidity raster file
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 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? 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) with z and z' given in meters. This could be easily applied using r.mapcalc Jose Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- José A. Ruiz-Arias Solar Radiation and Atmosphere Modelling Group http://www.ujaen.es/investiga/tep220 Physics Department, University of Jaén Campus Lagunillas, Building A3 066 23071 Jaén Spain Tlf: +34 953 212 474 Fax: +34 953 212 838 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Creating turbidity raster file
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
RE: [GRASS-user] Importing ArcASCII raster and associated dbf table?
-Original Message- From: Hamish [mailto:hamis...@yahoo.com] Sent: Tuesday, 10 February 2009 9:34 AM To: grass-user@lists.osgeo.org; Devereux, Jarl M Subject: RE: [GRASS-user] Importing ArcASCII raster and associated dbf table? 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 if not, is there a link to download the dataset? It seems to be a different data set here is a link to the one I'm using you have to register to download it. http://www.toolkit.net.au/Tools/SHPA This data set consists of nine sets of data files each containing the following: aas_*.asc An ArcASCII file containing a raster of the soil landscape index. aas_*.dbf A dBase file of soil properties linked to the soil landscape index. aas_*.prj A file containing projection information in the ArcGIS format. It is an interesting problem to me. Are you interested in one of the DBF fields, or several/all of them? As a quick guess, I think the answer will be to create 1 map per DB column by some simple script. Hamish The way I was going to use the atlas was to query it with either coordinates to see what soil properties are at that point or overlay other areas of interest to determine how the soils change across the area. I can do this externally to grass by referencing the CAT number outside grass but it isn't very convenient. On a similar situation I was doing this attribute management manually by querying the vector and manually adding a new column and defining the attribute. This makes me think I should be able to have all the table information in the one vector. I could also add the columns individually and then populate them external to grass through excel or something? Regards Jarl DISCLAIMER The information contained in the above e-mail message or messages (which includes any attachments) is confidential and may be legally privileged. It is intended only for the use of the person or entity to which it is addressed. If you are not the addressee any form of disclosure, copying, modification, distribution or any action taken or omitted in reliance on the information is unauthorised. Opinions contained in the message(s) do not necessarily reflect the opinions of the Queensland Government and its authorities. If you received this communication in error, please notify the sender immediately and delete it from your computer system network. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
RE: [GRASS-user] Importing ArcASCII raster and associated dbf table?
Jarl wrote: The data is an Australian soils atlas including hydraulic properties in the dbf. here is a link to the one I'm using you have to register to download it. http://www.toolkit.net.au/Tools/SHPA ok, not too hard. The trick is to convert the raster coverage to vector areas to take advantage of the rich vector attribute/database tie-ins: # import ArcASCII grid to UTM55S location r.in.gdal input=aas_z55.asc output=aas_z55 r.colors aas_z55.gdal color=random r.info aas_z55 # convert to vector areas, using raster value for category g.region rast=aas_z55 r.to.vect -v in=aas_z55.gdal out=aas_z55 feature=area # get rid of empty DBF table created by r.to.vect v.db.droptable aas_z55 # copy real DBF file into $MAPSET/dbf/ directory eval `g.gisenv` cp /tmp/au/aas_z55.DBF $GISDBASE/$LOCATION_NAME/$MAPSET/dbf/aas_z55.dbf # connect full DBF database to vector map, use VALUE as the key column v.db.connect map=aas_z55 driver=dbf table=aas_z55 key=value done! The way I was going to use the atlas was to query it with either coordinates to see what soil properties are at that point # manually query coords v.what -a map=aas_z55 east_north=595633,6568838 # if coords are already in a vector points map v.what.vect or overlay other areas of interest to determine how the soils change across the area. v.overlay then v.univar or one of the several report modules (or custom SQL query reports). I can do this externally to grass by referencing the CAT number outside grass but it isn't very convenient. automation is the magic key On a similar situation I was doing this attribute management manually by querying the vector and manually adding a new column and defining the attribute. This makes me think I should be able to have all the table information in the one vector. Yup, you're on the right track Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
