Re: [GRASS-user] Fonts missing
On Tue, May 27, 2008 at 6:44 PM, Paul Kelly [EMAIL PROTECTED] wrote: On Tue, 27 May 2008, Markus Neteler wrote: ... Despite the g.mkfontcap from Hamish, we need a fix in d.font to reject non-existing font names. Currently they are silently accepted (with zero effect of course). I just committed a fix to SVN to check the font names along the lines of what Glynn suggested, and also to check that the path existed and is a valid file in the case of the full path to a Freetype font. Please test - if there aren't any issues I think it would be a good idea to backport to 6.4. Great, thanks, Paul. I have polished the message and backported to 6.4. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Fonts missing
On Tue, May 27, 2008 at 11:13 PM, Richard Chirgwin [EMAIL PROTECTED] wrote: ... Fixed first time, thanks. sigh The number of things I still don't know about Grass! Richard, whenever possible, the software should suggest what to do in case it is clearly identifiable like here. The updated d.font now suggests to check for typo in the font name or to run g.mkfontcap :) Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] wxpython start fails
Hello, I start on a recent grass64svn (Ubuntu 8.04) with wxpython and get the followin message: [EMAIL PROTECTED]:~/src/grass6_devel$ grass64 -wxpython Cleaning up temporary files ... Starting GRASS ... Traceback (most recent call last): File /usr/local/grass-6.4.svn/etc/wxpython/gis_set.py, line 31, in module from gui_modules import globalvar File /usr/local/grass-6.4.svn/etc/wxpython/gui_modules/globalvar.py, line 48, in module CheckForWx() File /usr/local/grass-6.4.svn/etc/wxpython/gui_modules/globalvar.py, line 39, in CheckForWx except (ImportError, ValueError, wxversion.VersionError), e: UnboundLocalError: local variable 'wxversion' referenced before assignment Error in GUI startup. If necessary, please report this error to the GRASS developers. Switching to text mode now. Hit RETURN to continue... What is wrong? Thanks Manuel -- ___ Dr. Manuel Seeger Wiss. Assistent Scientific Assistant Physische GeographieDpt. of Physical Geography FB VI - Geographie/GeowissenschaftenGeography/Geosciences Universität Trier University of Trier D - 54286 Trier Tel.: +49-651-201 4557 Fax:+49-651-201 3976 Web:http://www-neu.uni-trier.de/index.php?id=9607 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] GUI equivalent of command
On Wed, May 28, 2008 at 6:23 AM, Nikos Alexandris [EMAIL PROTECTED] wrote: On Wed, 2008-05-28 at 08:58 +0530, H.S.Rai wrote: Is there any reference, where I can find GUI based equivalent of commands, as I want to introduce GRASS to beginners (avoiding command terminal to make them comfortable, to start with), like how to activate r.report with mouse. If there is such information, then I can use that, otherwise I will produce it and share with all. This dead link sounds very promising: http://wwwamb.bologna.enea.it/forgrass/documents/Grass-6-Visual-Tutorial.pdf Grab it here: http://web.archive.org/web/20060209031109/wwwamb.bologna.enea.it/forgrass/download.htm - VISUAL TUTORIAL FOR GRASS 6: 30 pages of tutorial (pdf) Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] d.thematic.area
On Tue, May 27, 2008 at 7:20 PM, kries [EMAIL PROTECTED] wrote: Thank you Markus. Sorry, I wasn't really very clear in my first posting. I have a freshly compiled grass, did distclean/clean beforehand. Two days ago I submitted a number of changes to run GRASS 7.svn (trunk) independently from GRASS 6. There were lots of GRASS 6 indications in trunk which I updated. Please be sure to have updated recently (yesterday or today). I do mean d.thematic.area: http://www.nabble.com/file/p17495633/dthematicarea01.jpg I see that Moritz didn't yet enable d.thematic.area for default compilation (I forgot that it exists...). I have now activated it in trunk for easier testing. Something got mixed up there, as you can see from my screenshots. Input and output must have got confused somehow. It's not just a GUI problem. It's the same with the command line. My old scripts don't work anymore. Sure - the GUI is autogenerated here. I have just compiled and it looks ok. Possible, that you run make distclean in grass64 instead of grass70? Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] RST - how does it chooses from near points?
I cant answer the question of how the v.surf.rst module choses points, but I can say that with the use of r.in.xyz you can control this quite well. For example, (and relating to the last question) you can use different methods to assign values to the cells from ASCII import. Meaning if you have multiple points in a cell, and you only want the values that are the minimum in the cell, use the 'method=min' option. Alternatively, you can use 'mean' and 'max' for other applications as well. then you would r.to.vect the data, as 3d points and then interpolate with v.surf.rst. If I recall right, by doing this r.in.xyz, it cut out about 35% of the processing time for v.surf.rst in my application. Especially for LiDAR data, I found the r.in.xyz module to be nothing short of spectacular. With the 'method=n' option, and an r.report raster unit=p, you can see the percent distribution of points per cell to help set cell size. These are some great tools! Mark On 5/27/08, Carlos Guâno Grohmann [EMAIL PROTECTED] wrote: It just came to me: if I have a vector map to be interpolated by RST, how does the algorithm chooses between points that are very close to each other? And if I have more than one point with the same coordinates, but different values? thks all Carlos -- +---+ Carlos Henrique Grohmann - Guano Geologist M.Sc - Doctorate Student at IGc-USP - Brazil Linux User #89721 - carlos dot grohmann at gmail dot com +---+ _ Good morning, doctors. I have taken the liberty of removing Windows 95 from my hard drive. --The winning entry in a What were HAL's first words contest judged by 2001: A SPACE ODYSSEY creator Arthur C. Clarke Can't stop the signal. ___ 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] wxpython start fails
Hi, which version of python and wxpython do you use? Martin 2008/5/28 Manuel Seeger [EMAIL PROTECTED]: Hello, I start on a recent grass64svn (Ubuntu 8.04) with wxpython and get the followin message: [EMAIL PROTECTED]:~/src/grass6_devel$ grass64 -wxpython Cleaning up temporary files ... Starting GRASS ... Traceback (most recent call last): File /usr/local/grass-6.4.svn/etc/wxpython/gis_set.py, line 31, in module from gui_modules import globalvar File /usr/local/grass-6.4.svn/etc/wxpython/gui_modules/globalvar.py, line 48, in module CheckForWx() File /usr/local/grass-6.4.svn/etc/wxpython/gui_modules/globalvar.py, line 39, in CheckForWx except (ImportError, ValueError, wxversion.VersionError), e: UnboundLocalError: local variable 'wxversion' referenced before assignment Error in GUI startup. If necessary, please report this error to the GRASS developers. Switching to text mode now. Hit RETURN to continue... What is wrong? Thanks Manuel -- ___ Dr. Manuel Seeger Wiss. Assistent Scientific Assistant Physische GeographieDpt. of Physical Geography FB VI - Geographie/GeowissenschaftenGeography/Geosciences Universität Trier University of Trier D - 54286 Trier Tel.: +49-651-201 4557 Fax:+49-651-201 3976 Web:http://www-neu.uni-trier.de/index.php?id=9607 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa * ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] wxpython start fails
Hi, well, when I uninstall wxPython I get a bit more informative error message Cleaning up temporary files ... Starting GRASS ... ERROR: Requested version of wxPython not found. wxPython = 2.8.1.1 is required. Detailed information in README file. Error in GUI startup. If necessary, please report this error to the GRASS developers. Switching to text mode now. Hit RETURN to continue... Martin 2008/5/28 Manuel Seeger [EMAIL PROTECTED]: Uups, I thought I hat installed wxPython... but it was not. Now it works Thanks a lot Martin! Manuel Martin Landa schrieb: Hi, which version of python and wxpython do you use? Martin 2008/5/28 Manuel Seeger [EMAIL PROTECTED]: Hello, I start on a recent grass64svn (Ubuntu 8.04) with wxpython and get the followin message: [EMAIL PROTECTED]:~/src/grass6_devel$ grass64 -wxpython Cleaning up temporary files ... Starting GRASS ... Traceback (most recent call last): File /usr/local/grass-6.4.svn/etc/wxpython/gis_set.py, line 31, in module from gui_modules import globalvar File /usr/local/grass-6.4.svn/etc/wxpython/gui_modules/globalvar.py, line 48, in module CheckForWx() File /usr/local/grass-6.4.svn/etc/wxpython/gui_modules/globalvar.py, line 39, in CheckForWx except (ImportError, ValueError, wxversion.VersionError), e: UnboundLocalError: local variable 'wxversion' referenced before assignment Error in GUI startup. If necessary, please report this error to the GRASS developers. Switching to text mode now. Hit RETURN to continue... What is wrong? Thanks Manuel -- ___ Dr. Manuel Seeger Wiss. Assistent Scientific Assistant Physische GeographieDpt. of Physical Geography FB VI - Geographie/GeowissenschaftenGeography/Geosciences Universität Trier University of Trier D - 54286 Trier Tel.: +49-651-201 4557 Fax:+49-651-201 3976 Web:http://www-neu.uni-trier.de/index.php?id=9607 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- ___ Dr. Manuel Seeger Wiss. Assistent Scientific Assistant Physische GeographieDpt. of Physical Geography FB VI - Geographie/GeowissenschaftenGeography/Geosciences Universität Trier University of Trier D - 54286 Trier Tel.: +49-651-201 4557 Fax:+49-651-201 3976 Web:http://www-neu.uni-trier.de/index.php?id=9607 -- 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 start fails
Uups, I thought I hat installed wxPython... but it was not. Now it works Thanks a lot Martin! Manuel Martin Landa schrieb: Hi, which version of python and wxpython do you use? Martin 2008/5/28 Manuel Seeger [EMAIL PROTECTED]: Hello, I start on a recent grass64svn (Ubuntu 8.04) with wxpython and get the followin message: [EMAIL PROTECTED]:~/src/grass6_devel$ grass64 -wxpython Cleaning up temporary files ... Starting GRASS ... Traceback (most recent call last): File /usr/local/grass-6.4.svn/etc/wxpython/gis_set.py, line 31, in module from gui_modules import globalvar File /usr/local/grass-6.4.svn/etc/wxpython/gui_modules/globalvar.py, line 48, in module CheckForWx() File /usr/local/grass-6.4.svn/etc/wxpython/gui_modules/globalvar.py, line 39, in CheckForWx except (ImportError, ValueError, wxversion.VersionError), e: UnboundLocalError: local variable 'wxversion' referenced before assignment Error in GUI startup. If necessary, please report this error to the GRASS developers. Switching to text mode now. Hit RETURN to continue... What is wrong? Thanks Manuel -- ___ Dr. Manuel Seeger Wiss. Assistent Scientific Assistant Physische GeographieDpt. of Physical Geography FB VI - Geographie/GeowissenschaftenGeography/Geosciences Universität Trier University of Trier D - 54286 Trier Tel.: +49-651-201 4557 Fax:+49-651-201 3976 Web:http://www-neu.uni-trier.de/index.php?id=9607 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- ___ Dr. Manuel Seeger Wiss. Assistent Scientific Assistant Physische GeographieDpt. of Physical Geography FB VI - Geographie/GeowissenschaftenGeography/Geosciences Universität Trier University of Trier D - 54286 Trier Tel.: +49-651-201 4557 Fax:+49-651-201 3976 Web:http://www-neu.uni-trier.de/index.php?id=9607 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] d.thematic.area
Hi, :clap: it works! Thanks for your help! kries -- View this message in context: http://www.nabble.com/d.thematic.area-tp17487551p17509812.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] syncing two locations
I completely agree to use your approach, with the shared location with multiple mapsets. Using something like Unison sync's the locations well, but has problems with single file reconciliation such as an sqlite.db file for all the vector attributes. Now I have some vector files that dont have attributes. :-( It seems to work fine for rasters, but cant say for absolute. I would stick with your approach, or use the linux clustering approach like you mentioned. thanks for the feedback, this will be a great help to distribute the processing load. Mark On 5/22/08, Markus Neteler [EMAIL PROTECTED] wrote: (I take liberty to cc the list again for discussion) On Thu, May 22, 2008 at 3:59 PM, M S [EMAIL PROTECTED] wrote: That is cool. I will definitely check out the wiki page. Perhaps I am not using locations or more specifically mapsets as efficiently as I could be. If I were to use a single location, but with different mapsets, would there be issues with changing the region boundary and/or the cell resolution? not at all, all mapsets are independent. For example, in my main area of interest (basin) there were areas I needed to do finer detailed analysis, going from 25 foot cells to 3 foot cells (r.in.xyz + r.report helped me find this optimal cell size, great stuff!), would this be an issue within the same location, but using a different mapset? I had thought that region settings applied to a location. no :) to a mapset! I have added http://grass.osgeo.org/wiki/Parallel_GRASS_jobs#Background I had (perhaps not the best choice, copied the location to another workstation) and changed the region resolution and boundary, and then ran more detailed analysis. If you do it serially, no problem. Otherwise it will conflict. Under this scenario, does it still make sense to use a single shared location? Most likely. One problem I foresee is with vectors using a single sqlite.db file, That's true. You cannot write simultaneously to the same sqlite.db file. http://en.wikipedia.org/wiki/SQLite A write access can only be satisfied if no other accesses are currently being serviced, otherwise the write access fails with an error code (or can automatically be retried until a configurable timeout expires). This concurrent access situation would change when dealing with temporary tables. But just use different names, so multiple sqlite.db files? in two different locations and then getting the vector database out of sync, and unable to rectify this with syncing the locations. It seems that reprojecting the data into the location is one way to deal with this single DB issue. I feel that it's making things worse. Most all of what I'm doing is time-consuming raster analysis though, and on the first sync, things seemed to work as expected. OK fine. Please let's collect ideas in the wiki how to deal with vector data and related DBs when processing in parallel. Markus Mark On 5/21/08, Markus Neteler [EMAIL PROTECTED] wrote: On Wed, May 21, 2008 at 5:58 PM, M S [EMAIL PROTECTED] wrote: I have two GRASS workstations performing various lidar and basin operations. My idea was to distribute the processing on various machines. As new data gets created in each location on each computer, when all processes are done, I would like to re-sync the locations to be the same. I did massive MODIS map processing on a cluster and used one location only in a shared network directory. Therein multiple mapsets with a final g.copy job to the mapset keeping the results. I have documented it here: http://grass.osgeo.org/wiki/Parallel_GRASS_jobs#PBS Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Open Source Geospatial Foundation http://www.osgeo.org/ http://www.grassbook.org/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] help
y need to make a poligon, but i make the boundary then the centroid and i dont know how to tell the program to make it a poligon.. can sombody help me? thanks a lot!! ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] syncing two locations || default DB - SQLite
One problem I foresee is with vectors using a single sqlite.db file, That's true. You cannot write simultaneously to the same sqlite.db file. ... But just use different names, so multiple sqlite.db files? IIRC you can make SQLite create per-vector map DBs by setting the [v.]db.connect database= string the appropriate way. I am not totally certain about that. If not it would be nice to have. By using a single file for the entire mapset's DB needs you much more quickly reach 2gb file sizes, make it less flexible for backups, have longer seek times, more chance for a single bit of data corruption to cascade etc. I don't know which way is better, for now I hope we can have both ways and race them against each other to see which one wins. I'll even press the issue. I just made SQLite the default DB in grass7/trunk. (using a single sqlite.db file per mapset) We'll see what happens.. Mark: I completely agree to use your approach, with the shared location with multiple mapsets. Using something like Unison sync's the locations well, but has problems with single file reconciliation such as an sqlite.db file for all the vector attributes. Now I have some vector files that dont have attributes. :-( It seems to work fine for rasters, but cant say for absolute. I would stick with your approach, or use the linux clustering approach like you mentioned. other ideas: NFS mount common dir; openMosix's distributed filesystem (seemingly built for this exact task); common Samba share; ... Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] help
On Thu, 2008-05-29 at 09:10 +, Isabel Alfaro Cardozo wrote: y need to make a poligon, but i make the boundary then the centroid and i dont know how to tell the program to make it a poligon.. can sombody help me? thanks a lot!! Hello! Some ideas: .Use v.info on your map and see how many areas are reported? Maybe you need to run v.clean to check for broken polygons (?). .Be sure that your boundary is closed. You will have to close it manually. .Add a centroid and Save Exit .In your GIS-Manager uncheck all except area to see whether you have area(s) or not Greetings, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
RE: [GRASS-user] help
y need to make a poligon, but i make the boundary then the centroid and i dont know how to tell the program to make it a poligon.. can sombody help me? thanks a lot!! You probably need to enable 'snapping' within the v.digit module to make the polygon vertices coincident. See the attached screenshot showing where to set this on the v.digit Settings. Once you have snapping enabled, you may still have to manually move vertices that were digitzed prior to enabling snapping. Use the 'Move point, line, boundary, or centroid' tool on the v.digit interface (see attached screenshot). Once you move a vertex on top of another with this tool, the two should snap an become coincident. I usually make the snapping parameter around a few tens of meters for projected (UTM-style) projections. ~ Eric. attachment: Screenshot-settings.pngattachment: Move_Boundary.png___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] GUI equivalent of command
Nikos: This dead link sounds very promising: http://wwwamb.bologna.enea.it/forgrass/documents/Grass-6-Visual-Tutorial.pdf Markus: Grab it here: http://web.archive.org/web/20060209031109/wwwamb.bologna.enea.it/forgrass/download.htm - VISUAL TUTORIAL FOR GRASS 6: 30 pages of tutorial (pdf) I really liked that tutorial and had hoped we could host it in the addons SVN tutorial/ dir*. Original Mac textedit source to the PDF is/was distributed with Lorenzo's GRASS 6.0 package so a good starting point is available. It would need to be updated for GRASS 6.3.0, maybe reimplement it using LaTeX or LyX for cross platform collaboration. I wait for Lorenzo's blessing before starting on that. [*] http://trac.osgeo.org/grass/browser/grass-addons/tutorials Lorenzo now links to this Italian 6.2 tutorial: http://www.geo.unipr.it/%7Eclerici/GRASS6.2.pdf (and also now has a 6.3.0 package ready) Also OpenOSX provides a quickstart: http://openosx.com/grass/OpenOSX-Grass-ReadMe-6.2.1.pdf Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] again help
i dont understand how to let the program know if im working with a point layer, a line layer or a poligon layer. the program allows me to have all kinds of elements in one file. how do i deal with this? thanks a lot for the help ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
RE: [GRASS-user] help
thank you a lot, im new in grass, im learning and i cant get the poligon area, i think the procedure is to use the boundary button, make it close then use the centroid option and then when i save the file and check the info it says i have a boundary and a centriod but no area, im i missing some steps? my snapping is the same as yours Please keep all responses on the Grass mailing list so that others may participate/respond. Thanks... Can you post a screenshot of your v.digit session? What happens when you manually move two vertices together, do they snap? Can you post the output from v.info? ~ Eric. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] again help
GRASS is not ArcGIS. So it doesn't have the same restrictions that ArcGIS does about points, lines, and polygons being in separate layers. Of course, YOU can make a vector file (i.e., layer) that is only points or lines or polygons, but you are not REQUIRED to do so. If you want to export vector objects from GRASS to shapefiles or something else, you can select only one feature type to export. GRASS also allows you to convert between vector types and converts from raster to vector. Why is this a problem? I always was annoyed about the requirement of ArcView and ArcGIS to have only one feature type per shapefile. Michael On 5/28/08 7:22 AM, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Date: Thu, 29 May 2008 09:15:05 +1930 From: Isabel Alfaro Cardozo [EMAIL PROTECTED] Subject: [GRASS-user] again help To: grass-user@lists.osgeo.org Message-ID: [EMAIL PROTECTED] Content-Type: text/plain; charset=iso-8859-1 i dont understand how to let the program know if im working with a point layer, a line layer or a poligon layer. the program allows me to have all kinds of elements in one file. how do i deal with this? thanks a lot for the help __ Michael Barton, Professor of Anthropology Director of Graduate Studies School of Human Evolution Social Change Center for Social Dynamics Complexity Arizona State University phone: 480-965-6213 fax: 480-965-7671 www: http://www.public.asu.edu/~cmbarton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] wxpython interface and r.lake
Hello everybody! Don't know if this is a bug, but the GUI for r.lake doesn't show any raster data from the drop-down menu. All is fine with tcltk. This is a screenshoot of the problem: http://img159.imageshack.us/my.php?image=rlakeinputob3.jpg Thank you very much Luca ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] GUI equivalent of command
On Wed, 2008-05-28 at 22:21 +0530, H.S.Rai wrote: On Wed, May 28, 2008 at 9:47 AM, Nikos Alexandris [EMAIL PROTECTED] wrote: [...] Just see what happened with me (before posting of my original email), when I clicked Raster Reports and statistics, i got pop up window for r.univar, making me wonder and confused. Now I found that it depends on position of GS-Manager on screen and position of mouse on Report and statistics. If I click near left edge of Report and statistics, then it works fine, otherwise it take click automatically on next window at the position of Univariate statistics. It may be sensitivity of my mouse or just a bad day for me :-( Hmmm...? It does not happen to me. Could it really rely on your mouse? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] segmentation fault r.li on 6.3.1svn (subject errata corrige)
For the list record: I could reproduce the problem with long file/path names on their machine. Buffer lengths are now corrected and the problem apparently fixed (6.3.svn, 6.4.svn and 7.trunk). Markus On Thu, May 22, 2008 at 12:31 AM, Markus Neteler [EMAIL PROTECTED] wrote: Hi Leo, I have tested with your data which you sent offlist to me in GRASS 6.4.svn: no problem - works. Markus On Wed, May 21, 2008 at 12:02 PM, Leonardo Lami [EMAIL PROTECTED] wrote: Hi all, I'm testing r.li.x on 6.3.1svn but I have some problems. I run r.li.setup setting a moving windows 60x60 cells in a region with a resolution 5x5 meters as sampling area. The region have this values: projection: 99 (Transverse Mercator) zone: 0 datum: rome40 ellipsoid: a=6378388 es=0.00672267002222 north: 4817400 south: 4812130 west: 1644645 east: 1651180 nsres: 5 ewres: 5 rows: 1054 cols: 1307 cells: 1377578 The sampling area file have this contenents: SAMPLINGFRAME 0|0|1|1 SAMPLEAREA -1|-1|0.0569259962049|0.0459066564652 MOVINGWINDOW When I run r.li.shannon I have this output: GRASS 6.3.1svn (Toscana):~ r.li.shannon [EMAIL PROTECTED] conf=vicarello output=uso_primavera_2005_vicarello_shannon [5]+ Doner.li.setup Segmentation fault WARNING: Can't open header file for [uso_primavera_2005 in (null)] CHILD[pid = 6124]: unable to open 2+ mask ... continue without!!! GRASS 6.3.1svn (Toscana):~ WARNING: Can't open header file for [uso_primavera_2005 in (null)] CHILD[pid = 6123]: unable to open 2+ mask ... continue without!!! WARNING: Can't open header file for [uso_primavera_2005 in (null)] CHILD[pid = 6122]: unable to open 2+ mask ... continue without!!! WARNING: Can't open header file for [uso_primavera_2005 in (null)] CHILD[pid = 6121]: unable to open 2+ mask ... continue without!!! WARNING: Can't open header file for [uso_primavera_2005 in (null)] CHILD[pid = 6120]: unable to open 2+ mask ... continue without!!! WARNING: Can't open header file for [uso_primavera_2005 in (null)] CHILD[pid = 6119]: unable to open 2+ mask ... continue without!!! WARNING: Can't open header file for [uso_primavera_2005 in (null)] CHILD[pid = 6118]: unable to open 2+ mask ... continue without!!! WARNING: Can't open header file for [uso_primavera_2005 in (null)] CHILD[pid = 6117]: unable to open 2+ mask ... continue without!!! WARNING: Can't open header file for [uso_primavera_2005 in (null)] CHILD[pid = 6116]: unable to open 2+ mask ... continue without!!! WARNING: Can't open header file for [uso_primavera_2005 in (null)] CHILD[pid = 6115]: unable to open 2+ mask ... continue without!!! The problem seems to be in the header raster file...but it doesn't seem particular strange: proj: 99 zone: 0 north: 4817400 south: 4812130 east: 1651180 west: 1644645 cols: 1307 rows: 1054 e-w resol: 5 n-s resol: 5 format: 2 compressed: 1 Have someone any ideas? Thank you very much for your help! Leonardo ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Open Source Geospatial Foundation http://www.osgeo.org/ http://www.grassbook.org/ -- Open Source Geospatial Foundation http://www.osgeo.org/ http://www.grassbook.org/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How many modules has GRASS?
On Wed, May 28, 2008 at 7:52 PM, Nikos Alexandris [EMAIL PROTECTED] wrote: Can they be counted automatically? Sure: GRASS 6.4.svn (spearfish60): ls -la $GISBASE/bin/*.* | wc -l 325 GRASS 6.4.svn (spearfish60): ls -la $GISBASE/scripts/*.* | wc -l 80 Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] d.thematic.area
On 28/05/08 10:15, Markus Neteler wrote: I see that Moritz didn't yet enable d.thematic.area for default compilation (I forgot that it exists...). I have now activated it in trunk for easier testing. Thanks! I had always meant to improve it a bit more and write the documentation (code and module) before making it compile by default, but I guess it's better that it gets more tester, and thus more potential coders, since I'm currently just too busy to do any work on that... :-( Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How many modules has GRASS?
On Wed, 2008-05-28 at 22:52 +0200, Nikos Alexandris wrote: ls -la $GISBASE/bin/*.* | wc -l Something like ls -la $GISBASE/bin/*.* | grep -F r. | wc -l #for raster ls -la $GISBASE/bin/*.* | grep -F v. | wc -l #for vector ls -la $GISBASE/bin/*.* | grep -F g. | wc -l #for generic etc. All in a script and a nice report. Or there other methods (better than grep -F) doing this? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] populate a table with values from another one
On 27/05/08 20:59, Glynn Clements wrote: Manuel Sangiao wrote: someone knows why this statement doesn't work? echo UPDATE grenzen gr SET gr.cleft=(SELECT adm.NAME FROM myadm0 amd WHERE gr.left= adm.cat) | db.execute Which database back-end are you using? The DBF driver only understands a fairly limited subset of SQL; in particular, it doesn't support subselects. And just as importantly, AFAIK, it does not support the use of more than one table in a statement, or ? Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How many modules has GRASS?
Hi, On 29.05.2008 00:04, Nikos Alexandris wrote: Or there other methods (better than grep -F) doing this? ls -1 $GISBASE/bin/r.* | wc -l #for raster etc. No need for grep or anything other then file names --Wolf -- :3 ) Wolf Bergenheim ( 8: ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How many modules has GRASS?
Nikos: Something like ls -la $GISBASE/bin/*.* | grep -F r. | wc -l #for raster ls -la $GISBASE/bin/*.* | grep -F v. | wc -l #for vector ls -la $GISBASE/bin/*.* | grep -F g. | wc -l #for generic etc. All in a script and a nice report. Or there other methods (better than grep -F) doing this? Wolf: ls -1 $GISBASE/bin/r.* | wc -l #for raster etc. No need for grep or anything other then file names and even the -1 is optional, it automatically switches to that if stdout is piped. also remember modules like xganim and nviz without a dot (if list method is not regex) and the few redundant wrapper scripts (r.cats, r.terraflow.short, ..) which are not seperate modules. I don't really see the worth in a counting script, but if you really want that the place to put it would be in the tools/ dir along with build_html_index.sh and module_synopsis.sh. another way, using tab completion: :) GRASS r.tabtab Display all 154 possibilities? (y or n) Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Extracting the external boundary of a vector file
Edward Tomlinson: I have a whole bunch of data from an RTK GPS. File exists as a MapInfo file. Imports into grass fine with area and centroids. My problem is however, i want to create another vector file that contains the mask or the external (an internal island) boundary. I have been able to create the external boundary fine through the following. v.category (delete categories) v.category (add category as 1 for all) v.dissolve (combine all boundaries with common category) This is great and works fine for the external boundary. However there is meant to be an island in there. Somewhere along the way (i think after the second call to v.category) it fills that island in. right. Has anyone got any hints or ideas? Essentially i am trying to create this vector to generate a mask for further interpolation of the point heights. I am assuming that i need this island removed (as there is no point data for it) otherwise it is going to effect the interpolation. Use v.overlay for the new and original map, operator=not, etc. or other combination of operator=and etc until you get what you want. another trick is v.extract list=1- to only extract features with a cat number. (then v.overlay) v.hull is another module to look at for this instead of v.dissolve. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How many modules has GRASS?
On Wed, 2008-05-28 at 20:09 -0700, Hamish wrote: Nikos: Something like ls -la $GISBASE/bin/*.* | grep -F r. | wc -l #for raster ls -la $GISBASE/bin/*.* | grep -F v. | wc -l #for vector ls -la $GISBASE/bin/*.* | grep -F g. | wc -l #for generic etc. All in a script and a nice report. Or there other methods (better than grep -F) doing this? Wolf: ls -1 $GISBASE/bin/r.* | wc -l #for raster etc. No need for grep or anything other then file names and even the -1 is optional, it automatically switches to that if stdout is piped. also remember modules like xganim and nviz without a dot (if list method is not regex) and the few redundant wrapper scripts (r.cats, r.terraflow.short, ..) which are not seperate modules. I don't really see the worth in a counting script, but if you really want that the place to put it would be in the tools/ dir along with build_html_index.sh and module_synopsis.sh. I was asked (some times) how many modules has GRASS and was curious to count. another way, using tab completion: :) GRASS r.tabtab Display all 154 possibilities? (y or n) Wowww! This is the fastest I suppose ;-) Cool! Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Triangles polygons with Arc*
Just curious if anybody encountered this issue with Arc*: it's a common problem. It does seem to be. short answer: TINs stink. Raster maps created from TINs stink more. My colleague used SRTM DEM converted to grid format of Arc He mentioned that it is fairly common when you convert raster grids in arc to vectors (shapefiles). Especially with very small areas like a 2-3 pixels. Will look into this further. To all the others who gave me a list of Arc mailing list thanks! I am reading a few right now. cheers, maning -- |-|--| | __.-._ |Ohhh. Great warrior. Wars not make one great. -Yoda | | '-._7' |Freedom is still the most radical idea of all -N.Branden| | /'.-c |Linux registered user #402901, http://counter.li.org/ | | | /T |http://esambale.wikispaces.com| | _)_/LI |-|--| ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] newbie: processing lidar with r.in.xyz
andrew haywood wrote: i am currently enjoying the flexibility of using r.in.xyz to process lidar data over native forests in Victoria, Australia. Thank you hamish for such a great tool. cheers, it is always nice to hear what it is used for. At this stage I have 'binned' my data into 20m cells and calculated a number of metrics (including 13 vegetation height percentiles p5 p10 p20 .. p90, p95, p99) based on the z values. What I was wondering is how can I calculate the mean intensity for the associated height percentiles. Just pass the intensity as the z= column instead of the height column? That is why I called it Column number of data values in input file instead of elevation values. It doesn't know/care what the data means. To get around this I have written a script in a proprietary stats package to create the intensity metrics. However, I would prefer to use Grass and opensource tools to do this binning what sort of metrics are you after? Any suggestions would be appreciated My data is in the following format x|y|ground|intensity|class|canopy_top|height 410774.45|5820999.93|773.3|23|10|766.0983886719|7.20161 410763.07|5820999.9|802.27|2|10|763.9434814453|38.3265 410765.47|5820999.94|773.11|90|10|764.2877807617|8.8 410758.09|5820999.95|807.99|47|10|762.8372192383|45.1528 410748.12|5820999.93|760.89|8|10|760.5731811523|0.316819 ... I suspect I should start to learn python. time well spent regardless; good return on investment if you like to think like that. a nice place to start: 10 minutes Python tutorial for programmers of other languages http://www.poromenos.org/tutorials/python you could probably do what you want quite simply in awk as well, but I'd only really suggest that if you know UNIX or C already. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] newbie: processing lidar with r.in.xyz
andrew haywood: To get around this I have written a script in a proprietary stats package to create the intensity metrics. However, I would prefer to use Grass and opensource tools to do this binning Any suggestions would be appreciated oh, if you are looking for an open source stats pacakge, R stats is the king. http://www.r-project.org/ Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Triangles polygons with Arc*
maning sambale: My colleague used SRTM DEM converted to grid format of Arc He mentioned that it is fairly common when you convert raster grids in arc to vectors (shapefiles). Especially with very small areas like a 2-3 pixels. Will look into this further. ok, then try and use GRASS's r.contour on the original data to create the vector instead of Arc's version and see if it is any nicer. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user