[GRASS-user] Fwd: 1st Call for Papers: Geomorphometry 2009, 29 August - 2nd September, Zurich, Switzerland
Apologies for cross postings FIRST CALL FOR PAPERS Geomorphometry 2009 29 August - 2 September 2009 Zurich, Switzerland http://2009.GEOMORPHOMETRY.ORG e-mail: [EMAIL PROTECTED] PROGRAM CHAIRS Ross Purves University of Zurich Stephan Gruber University of Zurich Tomislav Hengl University of Amsterdam KEY DATES Workshop proposals due 14 January 2009 Extended abstracts due 1 March 2009 Notification of acceptance 1 April 2009 Final camera-ready digital manuscripts due 1 May 2009 Author registration deadline 15 May 2009 Early registration deadline 15 May 2009 Geomorphometry 2009 Workshops 29 August 30 August 2009 Geomorphometry 2009 31 August - 2 September 2009 AIMS AND SCOPE The aim of Geomorphometry 2009 is to bring together researchers to present and discuss developments in the field of quantitative modelling and analysis of elevation data. Geomorphometry is the science of quantitative land-surface analysis and description at diverse spatial scales. It draws upon mathematical, statistical and image-processing techniques and interfaces with many disciplines including hydrology, geology, computational geometry, geomorphology, remote sensing, geographic information science and geography. The conference aims to attract leading researchers in geomorphometry presenting methodological advances in the field and to provide young researchers with an opportunity to present new results. The Geomorphometry 2009 conference will continue a series initiated by the Terrain Analysis and Digital Terrain Modelling conference hosted by Nanjing Normal University in November 2006. Possible topics include, but are not limited to: - Extraction of land-surface parameters from DEMs - Implications of novel data sources - Identification and classification of land-surface objects - Uncertainty in geomorphometry - Semantics of land-surface description - Visualisation in geomorphometry - Implications of scale and resolution - Flow and hydrological modelling using DEMs - Efficient methods for application to large data sets - Novel applications of geomorphometry - Planetary geomorphometry We specifically aim at papers with new methodological insights and thus papers which simply describe the application of GIS are discouraged. CONFERENCE PROGRAMME The conference programme will be based around a single track of papers, all of which will be subject to review in the form of extended abstracts by members of the programme committee. Criteria for paper acceptance will include relevance to the conference, novelty, scientific significance, relation to previous work in the domain and the quality of presentation. The proceedings will be made available both digitally and as printed working materials to attendees at the time of the conference and archived online. A special issue of a journal, to which authors will be invited to submit full papers after the conference is also planned. WORKSHOPS Geomorphometry will host up to three workshops, each with 15-30 attendees on the 29th and 30th August. We invite applications to host a workshop on a theme related to the main conference. Workshops should primarily take the form of either tutorials in a particular method or technique, or provide the opportunity for detailed discussion of upcoming topics. They should not simply be mini-conferences. If you are interested in organising a workshop, please mail [EMAIL PROTECTED] with WORKSHOP as the subject line, as well as a 1-page description of the workshop with the following headings: intended aims and scope, intended audience, outline workshop programme and technical requirements. SUBMISSIONS Prospective authors will be invited to submit extended abstracts of up to 2000 words by the above deadline through the EasyChair system. Formatting instructions and detailed information on using the submission system will be available in due course. Extended abstracts must be original works by the authors, not be currently under review in the same form by another outlet and not submitted elsewhere prior to the notification date. SCIENTIFIC COMMITTEE Alexander Brenning University of Waterloo, Canada Ian Evans Durham University, UK Peter Fisher University of Leicester, UK John Gallant CSIRO, Australia Paul Gessler University of Idaho, USA Stephan Gruber University of Zurich, Switzerland Tomislav Hengl University of Amsterdam, Netherlands Oliver Korup WSL, Switzerland Helena Mitasova North Carolina State University, USA Scott Peckham Rivix, USA Hannes Reuter Joint Research Centre, Italy Robert Weibel University of Zurich, Switzerland John Wilson University of Southern California, USA Jo Wood City University, UK Ralph Straumann University of Zurich, Switzerland Ross Purves University of Zurich, Switzerland Qiming Zhou Hong Kong Baptist University, Hong Kong -- +-+ Carlos Henrique Grohmann - Geologist D.Sc. a.k.a Guano - Linux User #89721 carlos dot grohmann at gmail dot com
[GRASS-user] running grass on a cluster
Dears, I am thinking about making my own linux cluster. Any good advice on running GRASS in a machine like this? cheers Carlos -- +-+ Carlos Henrique Grohmann - Geologist D.Sc. a.k.a Guano - 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
Re: [GRASS-user] nnbathy
I have not seen the original code / license, but would it be possible to directly integrate the code into GRASS? This would remove the dependency of an external nnbathy library. This would be great. cheers 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
Re: [GRASS-user] create 3d plane from vector lines
Thanks Soren, I got the effect I wanted in Paraview, using the linear extrusion filter. cheers Carlos On Wed, Sep 10, 2008 at 04:54, Sören Gebbert [EMAIL PROTECTED] wrote: Hi Carlos, how about using v.delaunay? If you have two lines for each fault, you can triangulate the space in between. But im not sure if this works with vertical lines. If you export those lines to VTK, you can triangulate vertical planes with the VTK delaunay triangulator in ParaView. Soeren Original-Nachricht Datum: Tue, 9 Sep 2008 23:32:37 -0300 Von: Carlos \\Guâno\\ Grohmann [EMAIL PROTECTED] An: CC: grass-user grass-user@lists.osgeo.org Betreff: Re: [GRASS-user] create 3d plane from vector lines Hello José. v.surf.rst is not exactly what I want here, although I could duplicate my vector lines with different elevations and then interpolate the values to get _raster_ surfaces. I was thinking more in a v.extrude way, just that it currently extrude areas to form facets or shifts lines, but it doesn't create planes from line (AFAICT). cheers Carlos On Tue, Sep 9, 2008 at 22:59, José María Michia [EMAIL PROTECTED] wrote: Hi Carlos. You can use v.surf.rst. Description (from command help): Spatial approximation and topographic analysis from given ... isoline data in vector format to floating point raster format using regularized spline with tension. Basic usage (remember to adjust region settings): v.surf.rst input=contour layer=0 elev=surface Where: contour : name of the map containing your vector lines (iso lines or level contours). surface : name of the raster map, obtained from interpolation of the vector lines. More info here: http://grass.itc.it/grass62/manuals/html62_user/v.surf.rst.html I hope that this help you. José Saludos (sorry for my bad english) 2008/9/9 Carlos Guâno Grohmann [EMAIL PROTECTED] Hi all, I'm stuck here, trying to create 3d vertical surfaces from 2d vector lines (to represent faults in nviz). hints? TIA 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 -- +---+ 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 -- GMX Kostenlose Spiele: Einfach online spielen und Spaß haben mit Pastry Passion! http://games.entertainment.gmx.net/de/entertainment/games/free/puzzle/6169196 -- +---+ 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
Re: [GRASS-user] create 3d plane from vector lines
Hey Mate, you know, simple dipping surfaces are easy, r.plane will do the job, but _complex_ dipping surfaces (like real faults) are still a work in progress.. cheers G On Wed, Sep 10, 2008 at 21:00, Paulo Marcondes [EMAIL PROTECTED] wrote: 2008/9/10 Carlos Guâno Grohmann [EMAIL PROTECTED]: Thanks Soren, I got the effect I wanted in Paraview, using the linear extrusion filter. Hey pal, does that method allows for representation of dipping surfaces? You know you should think about that too =] -- Paulo Marcondes = PU1/PU2PIX -22.915 -42.224 = GG86jc -- +---+ 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] create 3d plane from vector lines
Hi all, I'm stuck here, trying to create 3d vertical surfaces from 2d vector lines (to represent faults in nviz). hints? TIA 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
Re: [GRASS-user] create 3d plane from vector lines
Hello José. v.surf.rst is not exactly what I want here, although I could duplicate my vector lines with different elevations and then interpolate the values to get _raster_ surfaces. I was thinking more in a v.extrude way, just that it currently extrude areas to form facets or shifts lines, but it doesn't create planes from line (AFAICT). cheers Carlos On Tue, Sep 9, 2008 at 22:59, José María Michia [EMAIL PROTECTED] wrote: Hi Carlos. You can use v.surf.rst. Description (from command help): Spatial approximation and topographic analysis from given ... isoline data in vector format to floating point raster format using regularized spline with tension. Basic usage (remember to adjust region settings): v.surf.rst input=contour layer=0 elev=surface Where: contour : name of the map containing your vector lines (iso lines or level contours). surface : name of the raster map, obtained from interpolation of the vector lines. More info here: http://grass.itc.it/grass62/manuals/html62_user/v.surf.rst.html I hope that this help you. José Saludos (sorry for my bad english) 2008/9/9 Carlos Guâno Grohmann [EMAIL PROTECTED] Hi all, I'm stuck here, trying to create 3d vertical surfaces from 2d vector lines (to represent faults in nviz). hints? TIA 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 -- +---+ 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] Trend-surface analysis with GRASS+R: did something changed in sp package?
Hi all (sorry for cross-posting) I trying to do some trend-surface analysis in R using the spatial package (so I can go up to order 6). It's been a while since the last I did it, and my previously used steps are not working this time: R library(spgrass6) library(spatial) G-gmeta6() G gisdbase/home/guano/grassdata locationamsul_wgs84 mapset brasil rows3120 columns 2640 north 0 south -26 west-56 east-34 nsres 0.0083 ewres 0.0083 projection +proj=longlat +a=6378137 +rf=298.257223563 +no_defs +towgs84=0.000,0.000,0.000 # GridTopology with 0.1 degree resolution grd - GridTopology(cellcentre.offset=c(-55, -25), cellsize=c(0.1,0.1), cells.dim=c(560, 230)); #here is my data: aftaSE - readVECT6(afta_SE) # first-order TSA idade1-surf.ls(1,x=coords[,1], y=coords[,2], z=aftaSE$idade) ~so far, so good. now I want to get things as SpatialGridDataFrame, so I can get them back into GRASS: resid1-residuals(idade1); resid2-as.data.frame(resid1); resid.trend-SpatialPointsDataFrame(coords,resid2); resid.img - SpatialGridDataFrame(grd, resid.trend); Error in validObject(.Object) : invalid class SpatialGridDataFrame object: invalid object for slot data in class SpatialGridDataFrame: got class SpatialPointsDataFrame, should be or extend class data.frame This is where I stopped. I had these steps from emails in this list, and I remember they used to work. Has anything changed I am not aware of? TIA 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
Re: [GRASS-user] errors compiling on Zenwalk
There seems to be a problem with Your PNG library. Specifically, you appear to have multiple, incompatible versions installed with the same name. Thanks for your help, Glinn, but it doesn't look like I have more than one libpng installed: root[guano]# netpkg list |grep installed |grep png [I][l] Found installed libpng-1.2.27-i486-1.tgz on the repository 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
Re: [GRASS-user] Is the SRTM dataset a DEM or a DSM?
IMO its a DEM, since the z value represents elevation, including tree canopy and urban features. I personally don't like the term DSM because it can be related with any kind of Surface, and DEM (and DTM, when the z values represent the altitude) is related with the _topographical_ surface. cheers Carlos On Wed, Jun 25, 2008 at 15:31, Nikos Alexandris [EMAIL PROTECTED] wrote: In http://srtm.csi.cgiar.org/SRTMdataProcessingMethodology.asp it's called only DEM. In some papers it's called DSM. I'd rather say it's a DSM. But I would also like to know if there are objections ;-) Thank you, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- +---+ 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
Re: [GRASS-user] RST - how does it chooses from near points?
This is Helena's answer: Carlos, first point found at a given location is preserved and any consequent points within distance smaller than dmin are ignored (the z-values are not taken into account). I have done some experiments with 3D distance as well if you need to take into account the z - it is very easy to change that in the code but that would require adding a flag and I am not sure that many users would need it. If your data are contours, use v.generalize to thin the points first - it speeds up the processing a lot and the results are much better - see here: http://skagit.meas.ncsu.edu/~helena/grasswork/interpgen.html If you want to use mean value of your close points - the best approach is through r.in.xyz - see the example in its man page Helena cheers On Wed, May 28, 2008 at 5:19 AM, Markus Neteler [EMAIL PROTECTED] wrote: Maybe best ask to Helena who is on the road and not subscribed to this list AFAIK... please post the answer then! markus On Tue, May 27, 2008 at 9:39 PM, 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 -- Open Source Geospatial Foundation http://www.osgeo.org/ http://www.grassbook.org/ -- +---+ 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] RST - how does it chooses from near points?
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] rectangle behind colorscale in ps.map?
I just noticed that, in ps.map, while the scalebar has an option for a filled background, colorscale doesn't. I think it would be a nice option. What do you guys think? 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] geoid data from TOPEX
Hello all. I downloaded some geoid data from ftp://topex.ucsd.edu/pub/global_geoid_2min/ there are three files: geoid.egm96.grd ftp://topex.ucsd.edu/pub/global_geoid_2min/geoid.egm96.grd geoid.img.9.2 ftp://topex.ucsd.edu/pub/global_geoid_2min/geoid.img.9.2 geoid.img.9.2.ersftp://topex.ucsd.edu/pub/global_geoid_2min/geoid.img.9.2.ers the grd is a GMT grid with 7.5min spacing. The other is a ERMapper image file with 2min resolution. I managed to import geoid.img.9.2 with r.in.bin: r.in.bin -s -b input=/home/guano/geodata/TOPEX/geoid.img.9.2 output=geoid.img bytes=2 north=90 south=-90 east=0 west=360 rows=6336 cols=10800 --overwrite but there is a shift in position. When I compare with vector data for continents limits, I see that it is shifted to north. (see attached) Anyone have used this data? cheers 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. attachment: geoid.jpg___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] removal of duplicate points
maybe a script to remove each fifth line of the ascii? carlos On Wed, Mar 5, 2008 at 3:27 PM, Nikos Alexandris [EMAIL PROTECTED] wrote: On Wed, 2008-03-05 at 09:58 -0800, Dylan Beaudette wrote: Hi, Hi Dylan I have some squares which are stored as polygons. When converting these to points (v.to.points) I am left with the 5 points (first and last point are duplicated) that define the corners of the square. Is there a vector module (v.edit / v.clean) that could remove duplicate points? I would like to replace the manual operation with v.digit with something that can be scripted. Note that there are multiple squares, so a simple v.out.ascii | head -4 | v.in.ascii would be a little awkward. Ideas? Not really a solution but maybe... If polygons have first and last points duplicated what happens if they are converted in polylines? Duplicate coordinates still there? If not it could be an easy method to convert boundaries to lines and continue your task. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- +---+ 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] r.resamp.rsf fails
Hello all I just tried to resample a raster from 0:0:30 to 0:30 resolution, and calculating the curvatures at the same time. I have about 15Gb of free disk space, but after a while, the command fails complaining about not enough disk space... here is the output: Processing all selected output files will require 18304 bytes of disk space for temp files. Temporarily changing the region to desired resolution ... Bitmap mask created Changing back to the original region ... Percent complete: Not enough disk space--cannot write files interpolate() failed dnorm in mainc after grid before out1= 0.108012 split_and_interpolate() failed Running 6.3-SVN. cheers 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] a few suggestions and questions
Hello all first I'd like to make a few suggestions to enhance (IMO) our favorite GIS: - normalize the options in all modules. I noticed that while most commands take --overwrite, some (r.neighbors) use -o - some modules require a dummy value when output a list (r.mask, r.proj, v.proj..). Is this really necessary? this is a very confusing option (specially for newbies). - also, I think that _all_ fields (gui speaking now) where a raster/vector map is to be entered should have the select button, no matter if it is a field for new maps, because it makes the whole thing consistent, and it is easier for new users to understant how things work. now some issues: is plan curvature in r.param.scale broken? I can't get a map of it. profc works fine. Also is it possible to make r.param.scale work in latlong regions? (please?) this is just out of curiosity: if there are too many points to be used for interpolation by RST, how does it chooses which points will be used? thats all for now, cheers 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] GRASS tutorial in Portuguese (pt-br)
Hello all. This is just to announce that a new tutorial is available, in Brazilian Portuguese (pt-br). I wrote it over the last year, and tested it in short-course I gave last week. The plans are to keep updating it, with new sections and errors corrections. Some of it is based on the GDF tutorial (Dassau et al), and some is based on my own experience. http://www.igc.usp.br/pessoais/guano/grass.html cheers 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] circular pallete?
Hello all. does anyone know how I can draw a circular colorscale? maybe using some illustration program, saving as eps? The intent id to draw a circular colorscale for aspect maps, so it would look nicer than a rectangular one. cheers 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
Re: [GRASS-user] circular neighborhoods?
Yes, Eric, that is what it does for a _square_ neighborhood. I am a little confuse about the actual shape of the _circular_ neighborhood. using your examples: will a 3x3 circular neighborhood look like this:? ---o--- o-X-o ---o--- and a 5x5? ---o-o-o--- o-o-o-o-o o-o-X-o-o o-o-o-o-o ---o-o-o--- maybe now my question is more clear... thanks Carlos On Thu, Feb 21, 2008 at 3:07 PM, Patton, Eric [EMAIL PROTECTED] wrote: I was thinking about the circular neighborhoods that r.neighbors use. How exactly are they shaped? I mean, in a 3x3 window, does the circular one looks like a cross? and then it start to look more like a circle as the size increases? As far as I understand it, A 3x3 window would look like this: o-o-o o-X-o o-o-o Where the 'X' is the current cell being processed. Similarly, a 5X5 window would look like this: o-o-o-o-o o-o-o-o-o o-o-X-o-o o-o-o-o-o o-o-o-o-o r.neighbors will perform calculations on all the cells I've marked 'o' and and assign the output value to 'X', according to whatever method has been chosen (i.e., average, median, min, max, etc.). ~ Eric. -- +---+ 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
Re: [GRASS-user] circular neighborhoods?
Thanks for the quick answer Tom. Carlos On Thu, Feb 21, 2008 at 3:52 PM, Tom Russo [EMAIL PROTECTED] wrote: On Thu, Feb 21, 2008 at 03:11:43PM -0300, we recorded a bogon-computron collision of the [EMAIL PROTECTED] flavor, containing: Yes, Eric, that is what it does for a _square_ neighborhood. I am a little confuse about the actual shape of the _circular_ neighborhood. using your examples: will a 3x3 circular neighborhood look like this:? ---o--- o-X-o ---o--- and a 5x5? ---o-o-o--- o-o-o-o-o o-o-X-o-o o-o-o-o-o ---o-o-o--- maybe now my question is more clear... Piping in here after a quick glance at the source code... When the -c option to r.neighbors is selected, a neighborhood mask of booleans is created by this loop (in pseudocode): neighborhoodDistance=neighborhoodSize/2 (in integer math) for (i=0; ineigborhoodSize; i++) for (j=0; jneighborhoodSize; j++) mask[i][k]= ( (i-neighborhoodDistance)^2 + (j-neighborhoodDistance)^2 = (neighborhoodDistance)^2 So in a 3x3 example, neighborhoodSize=3, neighborhoodDistance=1, and your first diagram has os where the mask is true. The corners of the square wouldn't be in the mask (as in Eric's 3x3), because those points would have the left hand side of the inequality equal to 2, making the mask false. I haven't actually *tried* the code, but that is what the block of code in gather.c for circular neighborhoods says it should be doing. I was thinking about the circular neighborhoods that r.neighbors use. How exactly are they shaped? I mean, in a 3x3 window, does the circular one looks like a cross? and then it start to look more like a circle as the size increases? As far as I understand it, A 3x3 window would look like this: o-o-o o-X-o o-o-o Where the 'X' is the current cell being processed. Similarly, a 5X5 window would look like this: o-o-o-o-o o-o-o-o-o o-o-X-o-o o-o-o-o-o o-o-o-o-o r.neighbors will perform calculations on all the cells I've marked 'o' and and assign the output value to 'X', according to whatever method has been chosen (i.e., average, median, min, max, etc.). ~ Eric. -- +---+ 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 -- Tom RussoKM5VY SAR502 DM64ux http://www.swcp.com/~russo/ Tijeras, NM QRPL#1592 K2#398 SOC#236 AHTB#1 http://kevan.org/brain.cgi?DDTNM And, isn't sanity really just a one-trick pony anyway? I mean all you get is one trick, rational thinking, but when you're good and crazy, oooh, oooh, oooh, the sky is the limit! --- The Tick -- +---+ 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
Re: [GRASS-user] Interpolation with faultlines?
The attached picture is a sketch I made based on what I remember from a conference, long time ago, about some software that did interpolate with faultlines and how it worked. Carlos On Thu, Feb 21, 2008 at 4:42 PM, Maciej Sieczka [EMAIL PROTECTED] wrote: Moritz Lennert pisze: On 13/02/08 18:40, G. Allegri wrote: It's the same question I was going to ask. It would be an important advance in Grass, as I'm obliged to use commercial sw to do it (Esri) both to interpolate with barriers and to do 3D TINs with constraints. The answer was (and as far as I know still is) to either run v.surf.rst several times, once for each part within fault lines (using mask) or to use r.surf.nnbathy (AddOns). r.surf.nnbathy does not support fault lines either. FWIW, there is surfit [1] (FOSS). It supports fault lines AIUI - linear features which don't have a value assigned, but only depict a curve at which the interpolation should break. They call break/fault lines inequalities. I still haven't used surfit, but from examples [2] it seems interesting. [1]http://surfit.sourceforge.net [2]http://surfit.sourceforge.net/surfit/examples.html Maciek ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- +---+ 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. attachment: interp.png___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.out.gmt -- color palette with gaps?
Hi all. I have this color shaded relief map created with r.his and r.composite. I want to use it as a base map in GMT, sou I went to the AddOns page to get the r.out.gmt script, which creates the color palette (.cpt) for me. But when I try to create the map with GMT, I get this msg: [EMAIL PROTECTED] grdimage america_sul_srtm30plus_shade_color.grd -JM15c -R-56/-34/-26/0 -Camerica_sul_srtm30plus_shade_color.cpt amsul.ps grdimage: GMT Fatal Error: Color palette table america_sul_srtm30plus_shade_color.cpt has gaps - aborts! I remember using this script a while ago, and it worked fine.. thanks for any hint. 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] Interpolation with faultlines?
I remember that some time ago, there was some discussion on interpolation using faultlines. Any advances in that subject? cheers 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
Re: [GRASS-user] wxgui issue ( No module named wx)
Thanks Martin. I always forget about that dependency (fresh installed system).. cheers Carlos On Feb 11, 2008 2:04 PM, Martin Landa [EMAIL PROTECTED] wrote: No, 2008/2/11, Carlos Guâno Grohmann [EMAIL PROTECTED]: Hi, I just updated from svn, to test the new wxgui, but I got some errors: GRASS 6.3.svn (amsul_wgs84):~ Traceback (most recent call last): File /usr/local/grass-6.3.svn/etc/wxpython/wxgui.py, line 53, in module import gui_modules.utils as utils File /usr/local/grass-6.3.svn/etc/wxpython/gui_modules/utils.py, line 20, in module import gcmd File /usr/local/grass-6.3.svn/etc/wxpython/gui_modules/gcmd.py, line 29, in module import wx ImportError: No module named wx Should I just rename some file? Or is it more serious? no, you just need to install wxPython = 2.8.1, see http://svn.osgeo.org/grass/grass/trunk/gui/wxpython/README Martin -- Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa * -- +---+ 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] wxgui issue ( No module named wx)
Hi, I just updated from svn, to test the new wxgui, but I got some errors: GRASS 6.3.svn (amsul_wgs84):~ Traceback (most recent call last): File /usr/local/grass-6.3.svn/etc/wxpython/wxgui.py, line 53, in module import gui_modules.utils as utils File /usr/local/grass-6.3.svn/etc/wxpython/gui_modules/utils.py, line 20, in module import gcmd File /usr/local/grass-6.3.svn/etc/wxpython/gui_modules/gcmd.py, line 29, in module import wx ImportError: No module named wx Should I just rename some file? Or is it more serious? thanks 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
Re: [GRASS-user] How to calculate the area of a rastermap
r.surf.area? r.stats? Carlos On Jan 14, 2008 5:47 AM, Philipp Steigenberger [EMAIL PROTECTED] wrote: Hi list, I know there is a posibility to calculate the area of a rastermap, but I cannot find the command. I have a water.outlet map. The area data range is 1 to 1 Thankyou best wishes Philipp ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- +---+ 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