Re: [GRASS-user] r.profile grass64_rc5 - grass7
Hi Hamish, many thank you for your messege, Hamish wrote: Francesco Mirabella wrote: I am trying to use the r.profile module on a DEM. I have tried on both grass64_rc5 and grass7 with wxpython gui. In 64_rc5 there still exist the old option which allowed the user to interactively select the end-points by mouse, but the function does not appear to work as I get: r.profile -i input=DEM output=profile Using resolution [90.0161] No socket to connect to for monitor cairo7. - I can use it only with the oldtcltk Right, the new GUIs do not use Xmon windows (started with d.mon) and in general you can not use d.* and interactive xmon modules with them. You have to use the built-in GUI tools instead. Ok, this is what I rechoned In this case use the GUI profile tool, which is in the map display window just to the right of the magnifying glass buttons. The button looks like a little line graph. In grass 7 the possibility to choose the end-points with mouse is not there any more as the interactive Xmonitors have been completely removed there. (no more 'd.mon x0') By the way, the display profiler in the monitor works good in 64_rc5 - not i grass7 (I guess it's the old d.profile) The wxGUI profile tool from the map display window button? Yes, it is really a nice tool with beatifull graphical rendering but I cannot see any way to export the data in ascii format (distance-elevation) to load them into a spreadsheet if you use the query raster map button in the map display window GUI you will see the easting,northing,value text printed to the layer manager output window. you can then Save output from there or highlight+right-click Copy the end points to a text file and then pass that to r.profile. Not not very smooth- file a wish in the trac bug/wish system to save the profile values to a text file from the gui profile tool if you like. John Tate wrote: I used d.what.rast as I had ground control points, from GPS, for the profiles with which to validate the DEM. Add columns and update them with d.what.rast and export (.csv), then I put that together with the distance info from the GPS data. v.what.rast might help, or r.what. (I prefer v.rast.stats with a small buffer around each point to smooth any noise; some scripts in addons for aiding with that) Francesco: by browsing the man pages I have found that r.transect should do the work. Howevere at the moment it works by telling the module east, north, azimuth and distance, r.transect is just a wrapper around the r.profile module. All it does is convert the heading,distance to x2,y2 for you. while maybe it would be good to have also another couple of options like: 1 choose a vector make with sql query to tell the module where to sample elevations v.what.rast 2 be able to draw a cross sections and get and ascii file distance-elevations (this could be achieved by d.profile if there was such an option) isn't this what the r.profile output gives? Yes, however, maybe the possibility to make a straight ascii export of distance-elevation data from the wxGUI profile tool would be good. (d.profile did this before by giving plot.A, plot.B etc..). At present the user can export images (png, svg ect.) but not the data from the wxGUI profile tool. I think this would be good if one needs to plot more data along a topographic profile (in the same graph). Of course at present the r.profile command-line workaround works good, but maybe a more intuitive tools would be good (in the wxGUI mood :-)) best wishes Francesco At present the only way I can get a topographic profile to load into a spreadsheet is to use the oldtcltk, query the raster (d.where) and then use the d.where output as the input data for r.profile (still in oldtcltk) which gives me an ascii distance-elevation file. Indetail: cat d.where-output | r.profile input=dem output=elevation_points fwiw you can simplify: d.where | r.profile input=rastermap output=- Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.out.gdal and INTERLEAVE=PIXEL
Hello everybody, I need to export some raster data from GRASS to ArcView; as suggested on the wiki [1] I used the option INTERLEAVE=PIXEL with r.out.gdal: r.out.gdal in=elevation.10m format=GTiff type=Float32 createopt=INTERLEAVE=PIXEL out=output.tif The problem is that, using gdalinfo, INTERLEAVE is still set to BAND. Any idea? Using GRASS 6.4.0svn and GDAL 1.7.1 Thanks Luca -- View this message in context: http://n2.nabble.com/r-out-gdal-and-INTERLEAVE-PIXEL-tp4611081p4611081.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] open source alternative to ecognition?
Hello everybody, I am currently looking for an open source alternative to eCognition. I know the issue about how to emulate eCognitions object-based approach to image classification into GRASS has been dicussed before within this list, but I was wondering if there are any news considering this. As I am new to GRASS and to this mailing list, I would be especially interested in your opinion about the feasibility of such an emulation and would be thankful for any remarks or information about ongoing projects and already existing tools. Cheers, Marc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
SV: [GRASS-user] Help installing from http://trac.osgeo.org/grass/browser/grass-addons/
Thank you for responding. You are right. I installed grass with sudo and now tried as a normal user. It seems now that the permission and owner should be right. However, despite this, I have got a new error message today: U ps.output Checked out revision 41157. Compiling ps.output... /usr/local/grass-6.5.svn/include/Make/Module.make:25: warning: overriding commands for target `install' /usr/local/grass-6.5.svn/include/Make/Rules.make:82: warning: ignoring old commands for target `install' make: *** No rule to make target `default', needed by `first'. Stop. ERROR: Compilation failed, sorry. Please check above error messages. Any ideas ? Martin Fra: neteler.os...@gmail.com [neteler.os...@gmail.com] p#229; vegne av Markus Neteler [nete...@osgeo.org] Sendt: 19. februar 2010 16:10 Til: Martin Album Ytre-Eide Kopi: grass-user@lists.osgeo.org Emne: Re: [GRASS-user] Help installing from http://trac.osgeo.org/grass/browser/grass-addons/ On Fri, Feb 19, 2010 at 11:48 AM, Martin Album Ytre-Eide martin.album.ytre-e...@nrpa.no wrote: Hello. I updated my GRASS 6.5 and now I get a different error: mkdir -p /usr/local/src/grass6_devel/dist.i686-pc-linux-gnu/etc/paint/patterns/ /usr/bin/install -c -m 644 patterns/*.eps /usr/local/src/grass6_devel/dist.i686-pc-linux-gnu/etc/paint/patterns/ Installing ps.output... touch: cannot touch `/usr/local/grass-6.5.svn/g.extension.12845': Permission denied ERROR: Cannot write to /usr/local/grass-6.5.svn, installation failed The permission should be alright, but the catalog /usr/local/grass-6.5.svn/g.extension.12845 does not seem to exist. Please check the file permissions in /usr/local/grass-6.5.svn I suspect that it was installed with root permissions/sudo before and now you try as normal user or so. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Rectify an image
How can r.region put my image in place? Because it's a geostationary data so the image corners don't have the same coordinates (e.g. Top Left and Top Right don't have the same latitude). But I do have some GCP's for the image. So r.region doesn't help me (as far as I can see). Nikos On Sun, Feb 21, 2010 at 1:35 PM, Maciej Sieczka msiec...@sieczka.orgwrote: Nikos Dumakis pisze: Greetings I have an image without any GCP's but, teorethically it's in WGS84. I have all pixel positions outside that file. My question is, what should I do in GRASS to rectify/project this image? I thought in two solutions If the image doesn't require warping, but is just missing a georereference, use r.region to put in place. Maciek -- Maciej Sieczka http://www.sieczka.org ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Error while installing latest Weekly Snapshot
With the newer snapshot, it's ok THanks Markus Antonio António Rocha wrote: Hi I have just installed latest Generic GNU/Linux Weekly Snapshot for Linux in my Ubuntu OS (http://grass.itc.it/grass64/binary/linux/snapshot/grass-6.4.0svn-i686-pc-linux-gnu-13_02_2010.tar.gz). I configured with: ./configure --with-wxwidgets --with-blas --with-lapack --wils --with-cxx --with-python=/usr/bin/python2.6-config --with-gdal --enable-largefile --with-wxpython --with-sqlite --with-fftw --with-motif --with-tcltk-includes=/usr/include/tk --with-freetype-includes=/usr/include/freetype2 --enable-largefile When It starts, while it loads the first panel, an Error Window appears With a python error stating: Image file is not of type 3. After I select Location and Mapset the following error appears in terminal window: GRASS 6.4.0svn (nc_spm_08):~ Traceback (most recent call last): File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 1808, in module sys.exit(main()) File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 1801, in main app = GMApp(workspaceFile) File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 1728, in __init__ wx.App.__init__(self, False) File /usr/lib/python2.6/dist-packages/wx-2.8-gtk2-unicode/wx/_core.py, line 7978, in __init__ self._BootstrapApp() File /usr/lib/python2.6/dist-packages/wx-2.8-gtk2-unicode/wx/_core.py, line 7552, in _BootstrapApp return _core_.PyApp__BootstrapApp(*args, **kwargs) File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 1744, in OnInit workspace = self.workspaceFile) File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 136, in __init__ self.toolbar = self.__createToolBar() File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 345, in __createToolBar self.AddToolbarButton(self.toolbar, *each) File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 1076, in AddToolbarButton tool = toolbar.AddLabelTool(id=wx.ID_ANY, label=label, bitmap=icon, shortHelp=help) File /usr/lib/python2.6/dist-packages/wx-2.8-gtk2-unicode/wx/_controls.py, line 3764, in AddLabelTool shortHelp, longHelp, clientData) File /usr/lib/python2.6/dist-packages/wx-2.8-gtk2-unicode/wx/_controls.py, line 3686, in DoAddTool return _controls_.ToolBarBase_DoAddTool(*args, **kwargs) ValueError: invalid null reference in method 'ToolBarBase_DoAddTool', expected argument 4 of type 'wxBitmap const ' The commands seems to be working but the general GUI is not. If I insert a command in terminal windows, such as r.in.gdal, the r.in.gdal frame appears. If I do: g.gui wxpython I get: Traceback (most recent call last): File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 1808, in module sys.exit(main()) File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 1801, in main app = GMApp(workspaceFile) File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 1728, in __init__ wx.App.__init__(self, False) File /usr/lib/python2.6/dist-packages/wx-2.8-gtk2-unicode/wx/_core.py, line 7978, in __init__ self._BootstrapApp() File /usr/lib/python2.6/dist-packages/wx-2.8-gtk2-unicode/wx/_core.py, line 7552, in _BootstrapApp return _core_.PyApp__BootstrapApp(*args, **kwargs) File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 1744, in OnInit workspace = self.workspaceFile) File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 136, in __init__ self.toolbar = self.__createToolBar() File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 345, in __createToolBar self.AddToolbarButton(self.toolbar, *each) File /usr/local/grass-6.4.0svn/etc/wxpython/wxgui.py, line 1076, in AddToolbarButton tool = toolbar.AddLabelTool(id=wx.ID_ANY, label=label, bitmap=icon, shortHelp=help) File /usr/lib/python2.6/dist-packages/wx-2.8-gtk2-unicode/wx/_controls.py, line 3764, in AddLabelTool shortHelp, longHelp, clientData) File /usr/lib/python2.6/dist-packages/wx-2.8-gtk2-unicode/wx/_controls.py, line 3686, in DoAddTool return _controls_.ToolBarBase_DoAddTool(*args, **kwargs) ValueError: invalid null reference in method 'ToolBarBase_DoAddTool', expected argument 4 of type 'wxBitmap const ' Should I report to the error/bug list? Best regards Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 4879 (20100219) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 4886 (20100222) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 4886 (20100222) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org
Re: [GRASS-user] Rectify an image
Nikos, as Maciek said, your image seems to need warping, so I fear you have no other solution than georectify it. Vincent Le lundi 22 février 2010 à 11:03 +, Nikos Dumakis a écrit : How can r.region put my image in place? Because it's a geostationary data so the image corners don't have the same coordinates (e.g. Top Left and Top Right don't have the same latitude). But I do have some GCP's for the image. So r.region doesn't help me (as far as I can see). Nikos On Sun, Feb 21, 2010 at 1:35 PM, Maciej Sieczka msiec...@sieczka.org wrote: Nikos Dumakis pisze: Greetings I have an image without any GCP's but, teorethically it's in WGS84. I have all pixel positions outside that file. My question is, what should I do in GRASS to rectify/project this image? I thought in two solutions If the image doesn't require warping, but is just missing a georereference, use r.region to put in place. Maciek -- Maciej Sieczka http://www.sieczka.org ___ 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] Location wizard error
I have just installed latest snapshot and it's stil not listing EPSG's codes and descriptions. Anyway to debug this in order to understand what might not be ok? (GRASS version svn6.4.0 in Linux) On Fri, Feb 19, 2010 at 3:08 PM, Markus Neteler nete...@osgeo.org wrote: Hi, then please ask via list or package maintainer, I have Mandriva here and not much ideas about Karmic... Markus On Fri, Feb 19, 2010 at 1:25 PM, Gilbert Ferrara gilbertferrara1...@gmail.com wrote: Hello No. I installed Karmic GRASS binary and then svn6.4.0... nothing else On Fri, Feb 19, 2010 at 9:14 AM, Markus Neteler nete...@osgeo.org wrote: On Fri, Feb 19, 2010 at 10:06 AM, Gilbert Ferrara gilbertferrara1...@gmail.com wrote: According to download file, my SVN is from 16-01-2010 so it should be working. Is this problem it probably due to my GRASSbinary installation? Sorry, I have no idea. Perhaps you have some leftover old files there from a previous installation? By the way, in Location wizard what are the units for Spatial Resolution? meters, Km? The units are map units which depend on the projection. The selected project defines them. m On Thu, Feb 18, 2010 at 7:28 PM, Markus Neteler nete...@osgeo.org wrote: Hi Gilbert, On Thu, Feb 18, 2010 at 6:24 PM, Gilbert Ferrara gilbertferrara1...@gmail.com wrote: Hi Markus Ok I have first installed (Ubuntu 9.10) GRASS6.4.0~rc5.2 from KARMIC and then I installed GRASSsvn6.4.0, from GRASS website, at with a weekly snapshot from January 2010. I suppose that the karmic might be old but, since I have installed GRASSsvn6.4.0 from Jan2010, it would be ok, right? I think that it was fixed on 2010-jan-09. So it depends if yours is older or newer. Markus On Thu, Feb 18, 2010 at 4:09 PM, Markus Neteler nete...@osgeo.org wrote: On Thu, Feb 18, 2010 at 3:52 PM, Gilbert Ferrara gilbertferrara1...@gmail.com wrote: Hi I'm using GRASS6.4.0 in a Linux machine and, everytime I define a new Location, with Location wizard, The select EPSG code table is empty, I mean no EPSG is listed there... Any idea of what migh be happening or what file am I missing? To my knowledge, this bug has been fixed - how old is your installation? Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.profile grass64_rc5 - grass7
On Monday 22 February 2010 09:22:55 Francesco Mirabella wrote: Hi Hamish, many thank you for your messege, Hamish wrote: Francesco Mirabella wrote: I am trying to use the r.profile module on a DEM. I have tried on both grass64_rc5 and grass7 with wxpython gui. In 64_rc5 there still exist the old option which allowed the user to interactively select the end-points by mouse, but the function does not appear to work as I get: r.profile -i input=DEM output=profile Using resolution [90.0161] No socket to connect to for monitor cairo7. - I can use it only with the oldtcltk Right, the new GUIs do not use Xmon windows (started with d.mon) and in general you can not use d.* and interactive xmon modules with them. You have to use the built-in GUI tools instead. Ok, this is what I rechoned In this case use the GUI profile tool, which is in the map display window just to the right of the magnifying glass buttons. The button looks like a little line graph. In grass 7 the possibility to choose the end-points with mouse is not there any more as the interactive Xmonitors have been completely removed there. (no more 'd.mon x0') By the way, the display profiler in the monitor works good in 64_rc5 - not i grass7 (I guess it's the old d.profile) The wxGUI profile tool from the map display window button? Yes, it is really a nice tool with beatifull graphical rendering but I cannot see any way to export the data in ascii format (distance-elevation) to load them into a spreadsheet if you use the query raster map button in the map display window GUI you will see the easting,northing,value text printed to the layer manager output window. you can then Save output from there or highlight+right-click Copy the end points to a text file and then pass that to r.profile. Not not very smooth- file a wish in the trac bug/wish system to save the profile values to a text file from the gui profile tool if you like. John Tate wrote: I used d.what.rast as I had ground control points, from GPS, for the profiles with which to validate the DEM. Add columns and update them with d.what.rast and export (.csv), then I put that together with the distance info from the GPS data. v.what.rast might help, or r.what. (I prefer v.rast.stats with a small buffer around each point to smooth any noise; some scripts in addons for aiding with that) Cheers for pointing out my obvious mistake, v.what.rast is what I meant. I looked at v.rast.stats. Nice that it updates the table but for this purpose, it produces a lot of extra unneeded columns (shame you can't select which stats to add to the table - or perhaps you can?). Francesco: by browsing the man pages I have found that r.transect should do the work. Howevere at the moment it works by telling the module east, north, azimuth and distance, r.transect is just a wrapper around the r.profile module. All it does is convert the heading,distance to x2,y2 for you. while maybe it would be good to have also another couple of options like: 1 choose a vector make with sql query to tell the module where to sample elevations v.what.rast 2 be able to draw a cross sections and get and ascii file distance-elevations (this could be achieved by d.profile if there was such an option) isn't this what the r.profile output gives? Yes, however, maybe the possibility to make a straight ascii export of distance-elevation data from the wxGUI profile tool would be good. (d.profile did this before by giving plot.A, plot.B etc..). At present the user can export images (png, svg ect.) but not the data from the wxGUI profile tool. I think this would be good if one needs to plot more data along a topographic profile (in the same graph). Of course at present the r.profile command-line workaround works good, but maybe a more intuitive tools would be good (in the wxGUI mood :-)) best wishes Francesco At present the only way I can get a topographic profile to load into a spreadsheet is to use the oldtcltk, query the raster (d.where) and then use the d.where output as the input data for r.profile (still in oldtcltk) which gives me an ascii distance-elevation file. Indetail: cat d.where-output | r.profile input=dem output=elevation_points fwiw you can simplify: d.where | r.profile input=rastermap output=- Hamish ___ 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] r.out.gdal and INTERLEAVE=PIXEL
thedok78 wrote: I need to export some raster data from GRASS to ArcView; as suggested on the wiki [1] I used the option INTERLEAVE=PIXEL with r.out.gdal: r.out.gdal in=elevation.10m format=GTiff type=Float32 createopt=INTERLEAVE=PIXEL out=output.tif The problem is that, using gdalinfo, INTERLEAVE is still set to BAND. Any idea? Using GRASS 6.4.0svn and GDAL 1.7.1 yes, I can reproduce this: (grass 6.5 + gdal 1.5.2) #spearfish G6.5 r.out.gdal elevation.dem out=elev_pix.tif crea=INTERLEAVE=PIXEL G6.5 gdalinfo elev_pix.tif | grep INTERLEAVE INTERLEAVE=BAND Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: calculating maximum and minimum diameters with, GRASS or R and GRASS?
Is there no simpler way? e.g using the centroid as a center spot and then find the part of the polygon which is furthest away and use that distance for max radius? (and min dist for min radius) or make grass create a circle increasing in size untill it is fully outside the polygon then use that radius? Is this beyond GRASS? There must be a way using R - although I am still learning so I would be greatful of any instructions of how to do this in R (including which libraries to use) I may have to go to a R forum, any good ones? -- View this message in context: http://n2.nabble.com/Re-calculating-maximum-and-minimum-diameters-with-GRASS-or-R-and-GRASS-tp4577146p4611842.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] Edge file for r.le.pixel
Is there someone who have defined an edge file for r.le.pixel (e2 -analysis)? I have three attributes classes in my data, but my all efforts for the edge file end up with an announcement: "The edge file (r.le.para/edge) is incorrect since more/less than the 3 attributes found in the map are listed in the edge file)." What would be the right format for the file? I am running GRASS 6.4.0 under Linux. Anne Laita from Finland Anne LaitaDepartment of Biological and Environmental ScienceUniversity of JyväskyläP.O. Box 35,FIN-40014 University of JyväskyläFinland phone: +358 14 260 4206 or +358 40 541 9299e-mail: anne.la...@bytl.jyu.fiHomepage:http://users.jyu.fi/~annholo/___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Location wizard error
Gilbert Ferrara wrote: I have just installed latest snapshot and it's stil not listing EPSG's codes and descriptions. Anyway to debug this in order to understand what might not be ok? (GRASS version svn6.4.0 in Linux) In general don't be surprised if a generic linux build doesn't work properly on Ubuntu (or other distro) as location of system support files move about in an uncontrollable way. On Ubuntu 9.10/Karmic with the proj package installed the path to the EPSG file is /usr/share/proj/epsg. I've got the UbuntuGIS grass 6.4.0rc5 on ubuntu 9.10 running here, and the location wizard looks fine. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: Fwd: [GRASS-user] I.class empty Xmonitor
Greetings I have installed latest weekly snapshot update and it's fully working. Thank you Hamish for all your help Luisa 2010/2/18 Luisa Peña luisapena1...@gmail.com ah, ok. For the beta of the OSGeo live-dvd I have xUbuntu 9.10 on a bootable USB stick so I can test. Yes, I get the same problem. Uooou :) I thought i was going nuts regarding this :D This is with version 6.4.0~rc5-3~karmic1 from UbuntuGIS's launchpad.ppa repository. Are you using the same? Ok I have installed from Synaptic package manager, GRASS6.4.0~rc5.2 from KARMIC at /usr/lib/grass64. But aftter that I installed GRASSsvn6.4.0, from GRASS website, at /usr/local/grass-6.4.0svn with a weekly snapshot from January 2010 Since I'm using grasssvn6.4.0 does it make a difference??? How can I solve this? unistall Karmic GRASS? Thank you basically 6.4.0rc5 is about 6 months old now and this is almost certainly already fixed since then. the fix is going to be us releasing 6.4.0rc6, and then the various distro teams packaging that. I don't know when that will be, but I can only say that we are making good progress right now and I hope it will happen very soon. Regarding « d.mon x0; i.class » command, is it possible to have a wxpython GUI menu with that command? that's essentially what the existing menu item does. (or is supposed to be doing) (I will have non-experience users using i.class). fair enough. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] winGRASS
anyone developing winGRASS without using MSYS/Cygwin. (ie using mingw). If any of you working on this please let me know. I am happy to join your team. if you are porting GRASS GIS to windows without cygwin support please post svn source of winGRASS. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: [GRASS-dev] winGRASS
Hi, 2010/2/22 Mohammed Rashad rasha...@yahoo.com: anyone developing winGRASS without using MSYS/Cygwin. (ie using mingw). If any of you working on this please let me know. I am happy to join your team. if you are porting GRASS GIS to windows without cygwin support please post svn source of winGRASS. http://trac.osgeo.org/grass/wiki/CompileOnWindows 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] r.out.gdal and INTERLEAVE=PIXEL
On Mon, Feb 22, 2010 at 1:40 PM, Hamish hamis...@yahoo.com wrote: thedok78 wrote: I need to export some raster data from GRASS to ArcView; as suggested on the wiki [1] I used the option INTERLEAVE=PIXEL with r.out.gdal: r.out.gdal in=elevation.10m format=GTiff type=Float32 createopt=INTERLEAVE=PIXEL out=output.tif The problem is that, using gdalinfo, INTERLEAVE is still set to BAND. Any idea? Using GRASS 6.4.0svn and GDAL 1.7.1 yes, I can reproduce this: (grass 6.5 + gdal 1.5.2) #spearfish G6.5 r.out.gdal elevation.dem out=elev_pix.tif crea=INTERLEAVE=PIXEL G6.5 gdalinfo elev_pix.tif | grep INTERLEAVE INTERLEAVE=BAND IMHO best will be to open a GDAL ticket for this. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Location wizard error
In general don't be surprised if a generic linux build doesn't work properly on Ubuntu (or other distro) as location of system support files move about in an uncontrollable way. To avoid this kind of problems, what is your suggestion in order to update GRASS, using weekly snapshots? On Ubuntu 9.10/Karmic with the proj package installed the path to the EPSG file is /usr/share/proj/epsg. I've got the UbuntuGIS grass 6.4.0rc5 on ubuntu 9.10 running here, and the location wizard looks fine. Mine is pointing to: /usr/local/share/proj/epsg. I changed to your path and it worked. Question: Where do I define this path in order to change it automatically when I install/update GRASS? Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] g.region
Hi, with GRASS 6.4.0RC5 (Debian). I try to set the region with 250x250 meters of resolution but I have this result: GRASS 6.4.0RC5 (WGS84_33n):~ g.proj -p -PROJ_INFO- name : Universe Transverse Mercator proj : utm datum : wgs84 ellps : wgs84 zone : 33 no_defs: defined -PROJ_UNITS unit : metre units : metres meters : 1 GRASS 6.4.0RC5 (WGS84_33n):~ g.region nsres=250 ewres=250 GRASS 6.4.0RC5 (WGS84_33n):~ g.region -p projection: 1 (UTM) zone: 33 datum: wgs84 ellipsoid: wgs84 north: 4697295.976235 south: 4694979.412653 west: 265478.186612 east: 268489.71927 nsres: 257.39595356 ewres: 250.96105483 rows: 9 cols: 12 cells: 108 I try to set 250*250 but the result is 257.3*250.9. And if I try to convert a raster to a vector, the vector (point or line) is shifted from the raster. http://i49.tinypic.com/b7x7id.jpg Thanks, Alfredo ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] g.region
On Mon, 22 Feb 2010, Alfredo Alessandrini wrote: nsres: 257.39595356 ewres: 250.96105483 I try to set 250*250 but the result is 257.3*250.9. Alfredo, Use the g.region align switch, '-a'. Write: g.region res=250 -ap Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.to.vect Produces Too Many Areas
I've done something incorrectly but cannot find what that is. Ran r.water.outlet and produced the output basin map. Then I applied r.null to that map setting nulls to zeros. Next I applied r.to.vect to that map (with the '-s' option) setting feature to area. When I look at the aspect map overlaid by the sub-basin map I see the 5 extraneous cells that should not be there. You can see them on the attached screen shot. This did not happen the last time I used r.water.outlet, and the drainage map does not suggest to me why these 5 cells are included. Help would be appreciated in learning what happened and how to get this tail off the map. Thanks, Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.to.vect Produces Too Many Areas
I've done something incorrectly but cannot find what that is. Ran r.water.outlet and produced the output basin map. Then I applied r.null to that map setting nulls to zeros. Next I applied r.to.vect to that map (with the '-s' option) setting feature to area. When I look at the aspect map overlaid by the sub-basin map I see the 5 extraneous cells that should not be there. You can see them on the attached screen shot. This did not happen the last time I used r.water.outlet, and the drainage map does not suggest to me why these 5 cells are included. Help would be appreciated in learning what happened and how to get this tail off the map. Thanks, Richattachment: screenshot.jpg___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.to.vect Produces Too Many Areas
Since now a vector, can you select and delete those extraneous blocks? Mark On Mon, Feb 22, 2010 at 3:14 PM, Rich Shepard rshep...@appl-ecosys.com wrote: I've done something incorrectly but cannot find what that is. Ran r.water.outlet and produced the output basin map. Then I applied r.null to that map setting nulls to zeros. Next I applied r.to.vect to that map (with the '-s' option) setting feature to area. When I look at the aspect map overlaid by the sub-basin map I see the 5 extraneous cells that should not be there. You can see them on the attached screen shot. This did not happen the last time I used r.water.outlet, and the drainage map does not suggest to me why these 5 cells are included. Help would be appreciated in learning what happened and how to get this tail off the map. Thanks, Rich ___ 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] Question about least cost path
On Mon, 22 Feb 2010, M S wrote: Is the path terminating in local minimums? Mark, I don't know. It is my understanding that r.drain overcomes pits/sinks in a DEM by traversing lowest locations until goes off map. I did not use r.drain, but r.watershed. The last time I ran this sequence -- before correcting the reprojection and resolution -- there were no spurious cells. Thanks, Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.to.vect Produces Too Many Areas
On Mon, 22 Feb 2010, M S wrote: Since now a vector, can you select and delete those extraneous blocks? Mark, I'll try. Would still like to learn why the occur. Thanks, Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: grass-user Digest, Vol 46, Issue 63
Gary, Message: 2 Date: Mon, 22 Feb 2010 05:00:09 -0800 (PST) From: Gary Nobles garynob...@yahoo.com Subject: [GRASS-user] Re: calculating maximum and minimum diameters with, GRASS or R and GRASS? To: grass-user@lists.osgeo.org Message-ID: 1266843609307-4611842.p...@n2.nabble.com Content-Type: text/plain; charset=us-ascii Is there no simpler way? e.g using the centroid as a center spot and then find the part of the polygon which is furthest away and use that distance for max radius? (and min dist for min radius) or make grass create a circle increasing in size untill it is fully outside the polygon then use that radius? Is this beyond GRASS? There must be a way using R - although I am still learning so I would be greatful of any instructions of how to do this in R (including which libraries to use) I may have to go to a R forum, any good ones? There are a few approaches that occur to me. One thing I think you'll need is a map showing just the centroid locations (as points). You can do this by using v.extract, eg: v.extract input map output=output map type=centroid Then you could use v.distance measure distance from the centroids to all the boundaries in the original map, capture this in a database table, and select the min and max associated with each centroid. A question for others in the list: if you use the -a flag in v.distance, from centroid to line, I presume the measurement from the centroid is to all nodes on the line? Richard ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Bug or my fault?
Thanks Michael but I don't think this is the problem in my case, I am using the default DBF database, no DB connections required. Also remember that it all works well when I use TclTk as the GUI, and all works well when the number of DBF files in the mapset is less than approx 188. I haven't been able to find any other cause so far, but it could be Python related rather than GRASS, when I kill the GUI that hangs I get an error message to the effect that Python.exe is not responding. But why should the number of DBF files be relevant??? -- View this message in context: http://n2.nabble.com/Bug-or-my-fault-tp4609390p4614994.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] r.to.vect Produces Too Many Areas
On Mon, 22 Feb 2010, M S wrote: Since now a vector, can you select and delete those extraneous blocks? Mark, I've tried, by selecting both 'delete features' and 'delete nodes', but when I redraw the display they're still there. The wxPython GUI help web page doesn't tell me just how to do this. I exit and confirm the deletions, but they're still there. Probably user error. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: r.to.vect Produces Too Many Areas
Hi Rich, Message: 4 Date: Mon, 22 Feb 2010 12:14:03 -0800 (PST) From: Rich Shepard rshep...@appl-ecosys.com Subject: [GRASS-user] r.to.vect Produces Too Many Areas To: grass-us...@lists.osgeo.org Message-ID: alpine.lnx.2.00.1002221213470.20...@salmo.appl-ecosys.com Content-Type: text/plain; charset=us-ascii I've done something incorrectly but cannot find what that is. Ran r.water.outlet and produced the output basin map. Then I applied r.null to that map setting nulls to zeros. r,water.outlet produces a map with zeros instead of nulls, so you better must apply r.null to set zeros _to_ nulls Next I applied r.to.vect to that map (with the '-s' option) setting feature to area. Also, try without '-s' HTH Margherita -- Eng. Margherita Di Leo Ph.D. Candidate Methods and Technologies for Environmental Monitoring Department of Environmental Engineering and Physics (DIFA) University of Basilicata Campus Macchia Romana 85100 - Potenza Italy Office: +39-0971205363 Fax: +39-0971205160 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: r.to.vect Produces Too Many Areas
On Mon, 22 Feb 2010, Margherita Di Leo wrote: r,water.outlet produces a map with zeros instead of nulls, so you better must apply r.null to set zeros _to_ nulls Margherita, I stated it backwards. I understand that the setnull command does just that. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.to.vect Produces Too Many Areas
I haven't used the python GUI much. How about v.digit to graphically delete the blocks? Mark On Feb 22, 2010, at 4:33 PM, Rich Shepard rshep...@appl-ecosys.com wrote: On Mon, 22 Feb 2010, M S wrote: Since now a vector, can you select and delete those extraneous blocks? Mark, I've tried, by selecting both 'delete features' and 'delete nodes', but when I redraw the display they're still there. The wxPython GUI help web page doesn't tell me just how to do this. I exit and confirm the deletions, but they're still there. Probably user error. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.watershed: Interpreting Visual Output Map
When I read the man page for the visual output map of r.watershed I see that it's supposed to help me display results. I assume that it does this by coloring the accumulation of surface runoff with sub-basin sizes between 0 and the threshold value. When I display the map with d.rast.leg The color scale of the legend is blank. There are values from zero on the top to the threshold value at the bottom, but no colors are displayed. How am I to interpret this map? How does 'surface runoff accumulation with values modified for easy display' relate to minimum basin size? What are the units of surface runoff accumulation? Thanks in advance, Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Location wizard error
Hamish: In general don't be surprised if a generic linux build doesn't work properly on Ubuntu (or other distro) as location of system support files move about in an uncontrollable way. Gilbert Ferrara wrote: To avoid this kind of problems, what is your suggestion in order to update GRASS, using weekly snapshots? If you must have ubuntu-specific weekly snapshots, recompile the package from source. see the debian/ dir in the GRASS source code for details. Once you have the DebianGIS control files downloaded, often it is as easy as checking out the SVN for the 6.4 releasebranch, then in a weekly cron job do something like make distclean \ svn up \ debuild clean \ debuild binary \ dpkg -i ../grass*.deb (use pbuilder or cowbuilder instead of debuild if you like) On Ubuntu 9.10/Karmic with the proj package installed the path to the EPSG file is /usr/share/proj/epsg. Mine is pointing to: /usr/local/share/proj/epsg. I changed to your path and it worked. Question: Where do I define this path in order to change it automatically when I install/update GRASS? the default is set at build time, but you can override it. tcl/tk GUI: # add this to ~/.bashrc (or set before starting GRASS) GRASS_PROJSHARE=/usr/share/proj export GRASS_PROJSHARE wxGUI: Config - Preferences - Projection tab - EPSG file /usr/share/proj/epsg (sorry no file picker for that yet) or: edit $GISBASE/etc/Init.sh and look for GRASS_PROJSHARE Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user