Re: [GRASS-user] GRASS Location covering multiple UTM zones and SpatiaLite
Hello Stefan, As Markus mentioned, although SpatiaLite uses SQLite, it is a different entity. In the past, I used native dbms of GRASS. For the new project, I am thinking of switching to SpatiaLite, if it would be practical. The aim here is to record settlements (points), weather stations (points), stream network (lines), and basins (polygons) in one db to make some calculations easier; rather than extracting data from GRASS, exporting them into a separate SpatiaLite file. analyzing, and putting data back into GRASS' attribute data. It seemed to me that doing it with SpatiaLite would be easier but I am open to suggestions. I just started so I want to make sure that I build this right. As for QGIS (2.18.2) - GRASS (v. 7.0) issue –running on Mac OS 10.10.5 Yosemite–, when I install GRASS plugin, the default path for GRAS installation files is set to " /Users/Shared/unix/grass-qgis-mtlion/grass-7.0.4 ". I do not have GRASS installation files there and this creates an initiate error "GRASS init error: GRASS cannot found..." The only folder related to GRASS in my Users folder is "grass_dev" (sub-folders for v.6 and v.7) and some read-me files. I saw on the internet that some people suggested to move (,app) icons out of the folders but I am not certain because it was an old email thread. When I provide the correct path: AFAIK it is " /Applications/GRASS/GRASS-7. 0.app/Contents/MacOS " that does not seem to resolve the issue. All GRASS versions I have are kept in one folder titled GRASS under Applications. So, I am not sure how to solve this one. BÜLENT On Tue, Mar 28, 2017 at 11:32 AM, Blumentrath, Stefan < stefan.blumentr...@nina.no> wrote: > Hi Bülent, > > > > GRASS 7 uses SQLite for attribute management by default, so except for the > spatial part the DB backend it is identical. And a couple of modules can > handle points without topology. > > If you give us more information on what exactly you are planning to do > with your vector data it would be more easy to suggest a solution… E.g. are > you using polygons as input or output? What are you going to do with them… > > BTW. it should be fully possible to load GRASS data in QGIS 2.18. What is > your OS and GRASS version and what is the “path” issue you mentioned? > > > > Cheers > > Stefan > > > > > > > > *From:* Bulent Arikan [mailto:bulent.ari...@gmail.com] > *Sent:* tirsdag 28. mars 2017 09.37 > *To:* Blumentrath, Stefan <stefan.blumentr...@nina.no> > *Subject:* Re: [GRASS-user] GRASS Location covering multiple UTM zones > and SpatiaLite > > > > Hi Stefan, > > > > Thank you for your response! I guess that makes sense and I will define a > location for UTM 36 since it covers the majority of the region as well as > it is between 35 and 37. That would probably help to keep any distortion at > minimum. > > > > I understand that there is still not full compatibility with SpatiaLite, > so I will probably think something else. > > > > Cheers, > > > > BÜLENT ARIKAN > > Assistant Professor > > Department of Ecology and Evolution > > Eurasia Institute of Earth Sciences > Istanbul Technical University > Maslak-Sarıyer > Istanbul > 34469 > > Turkey > > http://web.itu.edu.tr/~bulentarikan/index.html > [image: cid:image001.jpg@01D2A7AA.458E7A90] > > > > On Mar 28, 2017, at 10:07 AM, Blumentrath, Stefan < > stefan.blumentr...@nina.no> wrote: > > > > Hi Bülent, > > > > Here in Norway we face the same situation. The country covers UTM zones 31 > to 36. However, the standard the public authorities suggest is to use the > average (UTM 33) when running analysis on the whole country. > > Of course this leads to imprecision and distortion, but if that is > significant or not depends on your application (esp. scale of the problem > you are going to solve)… > > > > Since GRASS has a topological vector format, linking external vector data > can cause issues in your analysis, esp. with polygons. However, points can > be fine in many cases… Please have a look at the mauals of v.external [1] > and v.external.out [2]. > > > > Cheers, > > Stefan > > > > 1: https://grass.osgeo.org/grass72/manuals/v.external.html > > 2: https://grass.osgeo.org/grass72/manuals/v.external.out.html > > > -- Asst. Prof. BÜLENT ARIKAN Department of Ecology and Evolution Eurasia Institute of Earth Sciences Istanbul Technical University Maslak-Sariyer, Istanbul TURKEY 34469 ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] GRASS Location covering multiple UTM zones and SpatiaLite
Dear List, I have two separate questions for Grass 7.0, which has been running on Mac OS 10.10.5 Yosemite: 1) My project encompasses three UTM zones (35-36-37 in northern hemisphere). About 80% of the area falls within UTM 36, 15% in UTM 37, and 5% in UTM 35. I know that in principle, each GRASS location should be defined to cover one UTM zone for accuracy. I will be doing a lot of terrain analyses, including hydrology. I simply cannot work in three separate zones/locations since this is one continuous piece of land. I can work in LatLong location but I am not sure if this will cause any problems in the accuracy of raster analyses. I will appreciate any suggestions. 2) For the same project, I have been thinking of using SpatiaLite rather than PostgreSQL or other DBF managers. I will have a lot of point data and some polygons. I could not find much about Grass GIS-SpatiaLite compatibility (except for a book in Italian). Is there a resource where I can get informed about using SpatiaLite in Grass? I understand that QGIS has SpatiaLite functionality but I noticed that I could not load my Grass data in QGIS (2.18.2 Las Palmas) because of a “path” issue. If possible, my priority would be to go with Grass-Spatialite combo. Thank you for your time! BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] issues with addon scripts
Dear List, I am having issues with adding on scripts to GRASS 7.0 stable (running Mac OS 10.10.2 Yosemite). I am the Administrator. I have the “Modules” folder that contains bin, scripts, etc. in /Users/bulentarikan/Library/GRASS/7.0 The add-on base is defined as: bash: /Users/bulentarikan/Library/GRASS/7.0/Modules: is a directory The path is: bash: /Applications/GRASS_7.0/GRASS-7.0.app/Contents/MacOS/bin:/Applications/GRASS_7.0/GRASS-7.0.app/Contents/MacOS/scripts:/Users/bulentarikan/Library/GRASS/7.0/Modules/bin:/Users/bulentarikan/Library/GRASS/7.0/Modules/scripts:-framework: The script I am trying to add on is not available in the add-on repository as far as I know. I used BatChmod to make sure that I have the permissions of the (.py) script folder. Then I moved the script into the “scripts” folder located under (right click on the GRASS 7.0 app ContentsMacOS). In the scripts folder, all scripts have black icons with the word ‘exec’ on the upper left corner. The script I dropped into the scripts folder did not have that, so I used: chmod +x /Users/bulentarikan/Desktop/r.landscape.evol7.py However, that did not change anything! The file is still not executable. So, my first question is, how can I make the script file an executable file? In the Terminal, I call the script by typing its name. I get the error below: Does this mean that I also need to have the executable file in /Users/bulentarikan/Library/GRASS/7.0/Modules? Thank you! BÜLENT ARIKAN Assistant Professor Department of Ecology and Evolution Assistant Director Eurasia Institute of Earth Sciences Istanbul Technical University Maslak-Sarıyer Istanbul 34469 Turkey http://web.itu.edu.tr/~bulentarikan/index.html http://web.itu.edu.tr/~bulentarikan/index.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] R.stream.* add ons for GRASS 6.4 or 6.5 on Mac OS 10.8.5 (Mountain Lion)
Dear List, As far as I understand, for some of us R.stream.* add ons are not accessible right now due to bugs in these scripts that cause g.extension to fail. Until that is resolved, I am wondering if there is any user that can send me a zip file containing the 'exec' files of the following modules. I mostly use GRASS 6.4.3 and 6.5 versions on Mac OS Mountain Lion (10.8.5). If you have it for GRASS 7, it is OK as well. r.ipso r.stream.basins r.stream.distance r.stream.extract r.stream.order r.stream.preview r.stream.stats Thank you! BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.stream.* addons issue
Hi Ralf, If you do not mind, can you, please tell me which version of Grass you are running? I would also very much appreciate if you could send me the details of your g.extension parameters (esp, the url and path information). Maybe a screenshot! Thank you so much for your time. Bulent On Oct 19, 2013, at 10:51, Ralf Schäfer sena...@ecotoxicology.de wrote: Dear all I can confirm that g.extension now works for me on a mac (10.8.5) in the default path. Tried r.stream.order and r.stream.distance and both were successfully installed. Makes things much easier, great! Nice weekend Ralf Am 19.10.2013 um 05:23 schrieb William Kyngesburye wokl...@kyngchaos.com: On Oct 18, 2013, at 9:50 PM, Bulent Arikan wrote: Unfortunately, it still does not work. I downloaded the stable version of Grass from your website and I am having the same issue. Like I said, I haven't tried it from the GUI. It's different code (python) than the CLI version (shell script), so maybe something isn't working right in python. Try running g.extension in the GRASS Terminal, not from the GUI. Works fine for me in the default GRASS_ADDON_PATH prefix (not global). It is not recommended to install globally on OS X, that's the GRASS.app. Mainly, it's just not the Mac way to add something to an existing application. Global installation on OS X should be the /Library/GRASS addons folder, but I don't think there is a default behavior in g.extension for that yet. You woud have to manually specify the prefix and make sure you are logged in as an administrator user (both the -s and -u flags force the prefix to be the GRASS app). Something like: g.extension ext=r.stream.basins prefix=/Library/GRASS/6.4/Modules Meanwhile, I am trying to remember an old trick I learned from a friend. I saved the .html file as a plain text file and dropped it into the Scripts folder in ContentsMacOS. I used to type a command starting with 'chmod' to make that plain text file a Unix executable file. It worked in the past for several addons. That sounds dangerous - making html executable, anything could happen if there was some valid script code in the html. BÜLENT On Oct 19, 2013, at 1:39 AM, William Kyngesburye wokl...@kyngchaos.com wrote: Actually, I think g.extension should work now, I forgot to update the readme. I haven't tried it from the GUI, or installing system-wide. On Oct 18, 2013, at 3:36 PM, Bulent Arikan wrote: Thank you Ralf! BÜLENT On Oct 18, 2013, at 11:21 PM, Ralf Schäfer sena...@ecotoxicology.de wrote: Dear Bulent g.extension does not work on Mac to my knowledge. There is a detailed description of an alternative how to build modules on Mac using a few lines in the Shell. If you install GRASS from kyngchaos, he provides a description in a readme file. http://www.kyngchaos.com/software/grass I do not know whether this description is compatible with the installation from Michael Barton. Best regards Ralf ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user - William Kyngesburye kyngchaos*at*kyngchaos*dot*com http://www.kyngchaos.com/ Theory of the Universe There is a theory which states that if ever anyone discovers exactly what the universe is for and why it is here, it will instantly disappear and be replaced by something even more bizarrely inexplicable. There is another theory which states that this has already happened. -Hitchhiker's Guide to the Galaxy 2nd season intro ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user - William Kyngesburye kyngchaos*at*kyngchaos*dot*com http://www.kyngchaos.com/ Mon Dieu! but they are all alike. Cheating, murdering, lying, fighting, and all for things that the beasts of the jungle would not deign to possess - money to purchase the effeminate pleasures of weaklings. And yet withal bound down by silly customs that make them slaves to their unhappy lot while firm in the belief that they be the lords of creation enjoying the only real pleasures of existence - the wisdom of Tarzan ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.stream.* addons issue
Thanks Ralf, I think the second half of your email is reminding me that I should do the same and check my add on settings. I will give it a try! Thank you! Asst. Prof. BÜLENT ARIKAN Department of Ecology and Evolution Eurasia Institute of Earth Sciences Istanbul Technical University Maslak-Sarıyer Istanbul-TURKEY 34469 On Oct 19, 2013, at 11:36, Ralf Schäfer sena...@ecotoxicology.de wrote: Hi Bulent I am running version 6.4.3 from Williams page: http://www.kyngchaos.com/software/grass I do not change any parameters in g.extension. In the shell window I type g.extension and then enter the extension in the respective window. However, this translates to g.extension extension=r.stream.order svnurl=http://svn.osgeo.org/grass/grass-addons/grass6 Checking my gis environment with set | grep GRASS gives for the Addon Path: GRASS_ADDON_PATH=/Users/ralfs/Library/GRASS/6.4/Modules/bin:/Users/ralfs/Library/GRASS/6.4/Modules/scripts:/Library/GRASS/6.4/Modules/bin:/Library/GRASS/6.4/Modules/scripts Hope this helps Ralf Am 19.10.2013 um 09:59 schrieb Bulent Arikan bulent.ari...@gmail.com: Hi Ralf, If you do not mind, can you, please tell me which version of Grass you are running? I would also very much appreciate if you could send me the details of your g.extension parameters (esp, the url and path information). Maybe a screenshot! Thank you so much for your time. Bulent On Oct 19, 2013, at 10:51, Ralf Schäfer sena...@ecotoxicology.de wrote: Dear all I can confirm that g.extension now works for me on a mac (10.8.5) in the default path. Tried r.stream.order and r.stream.distance and both were successfully installed. Makes things much easier, great! Nice weekend Ralf Am 19.10.2013 um 05:23 schrieb William Kyngesburye wokl...@kyngchaos.com: On Oct 18, 2013, at 9:50 PM, Bulent Arikan wrote: Unfortunately, it still does not work. I downloaded the stable version of Grass from your website and I am having the same issue. Like I said, I haven't tried it from the GUI. It's different code (python) than the CLI version (shell script), so maybe something isn't working right in python. Try running g.extension in the GRASS Terminal, not from the GUI. Works fine for me in the default GRASS_ADDON_PATH prefix (not global). It is not recommended to install globally on OS X, that's the GRASS.app. Mainly, it's just not the Mac way to add something to an existing application. Global installation on OS X should be the /Library/GRASS addons folder, but I don't think there is a default behavior in g.extension for that yet. You woud have to manually specify the prefix and make sure you are logged in as an administrator user (both the -s and -u flags force the prefix to be the GRASS app). Something like: g.extension ext=r.stream.basins prefix=/Library/GRASS/6.4/Modules Meanwhile, I am trying to remember an old trick I learned from a friend. I saved the .html file as a plain text file and dropped it into the Scripts folder in ContentsMacOS. I used to type a command starting with 'chmod' to make that plain text file a Unix executable file. It worked in the past for several addons. That sounds dangerous - making html executable, anything could happen if there was some valid script code in the html. BÜLENT On Oct 19, 2013, at 1:39 AM, William Kyngesburye wokl...@kyngchaos.com wrote: Actually, I think g.extension should work now, I forgot to update the readme. I haven't tried it from the GUI, or installing system-wide. On Oct 18, 2013, at 3:36 PM, Bulent Arikan wrote: Thank you Ralf! BÜLENT On Oct 18, 2013, at 11:21 PM, Ralf Schäfer sena...@ecotoxicology.de wrote: Dear Bulent g.extension does not work on Mac to my knowledge. There is a detailed description of an alternative how to build modules on Mac using a few lines in the Shell. If you install GRASS from kyngchaos, he provides a description in a readme file. http://www.kyngchaos.com/software/grass I do not know whether this description is compatible with the installation from Michael Barton. Best regards Ralf ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user - William Kyngesburye kyngchaos*at*kyngchaos*dot*com http://www.kyngchaos.com/ Theory of the Universe There is a theory which states that if ever anyone discovers exactly what the universe is for and why it is here, it will instantly disappear and be replaced by something even more bizarrely inexplicable. There is another theory which states that this has already happened. -Hitchhiker's Guide to the Galaxy 2nd season intro ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user - William Kyngesburye kyngchaos
[GRASS-user] r.stream.* addons issue
Dear List, I am trying to add on several (r.stream.*) modules from (http://grasswiki.osgeo.org/wiki/AddOns/GRASS_6). I am running Mac OS 10.8.5 and I tried this both in GRASS 6.4.3 and 6.5 (latest builds by Michael Barton). I tried both typing the command into GUI and using g.extension (both command dialog and fetchinstall methods). In g.extension, I type the name of the extension (r.stream.basins), check -s flag (install system-wide, provide the SVN URL (https://svn.osgeo.org/grass/grass-addons/grass6/raster/r.stream.basins), I keep the prefix path at default ($GRASS_ADDON_PATH) (GRASS should ignore this anyway because of the -s flag) I am getting the error below: Traceback (most recent call last): File /Applications/GRASS_6.4.3/GRASS-6.4.3.app/Contents/M acOS/etc/gui/scripts/g.extension.py, line 828, in module sys.exit(main()) File /Applications/GRASS_6.4.3/GRASS-6.4.3.app/Contents/M acOS/etc/gui/scripts/g.extension.py, line 817, in main install_extension() File /Applications/GRASS_6.4.3/GRASS-6.4.3.app/Contents/M acOS/etc/gui/scripts/g.extension.py, line 499, in install_extension ret = install_extension_other() File /Applications/GRASS_6.4.3/GRASS-6.4.3.app/Contents/M acOS/etc/gui/scripts/g.extension.py, line 520, in install_extension_other gui_list = list_wxgui_extensions(print_module = False) File /Applications/GRASS_6.4.3/GRASS-6.4.3.app/Contents/M acOS/etc/gui/scripts/g.extension.py, line 273, in list_wxgui_extensions f = urlopen(url) File /System/Library/Frameworks/Python.framework/Versions /2.6/lib/python2.6/urllib2.py, line 126, in urlopen return _opener.open(url, data, timeout) File /System/Library/Frameworks/Python.framework/Versions /2.6/lib/python2.6/urllib2.py, line 397, in open response = meth(req, response) File /System/Library/Frameworks/Python.framework/Versions /2.6/lib/python2.6/urllib2.py, line 510, in http_response 'http', request, response, code, msg, hdrs) File /System/Library/Frameworks/Python.framework/Versions /2.6/lib/python2.6/urllib2.py, line 435, in error return self._call_chain(*args) File /System/Library/Frameworks/Python.framework/Versions /2.6/lib/python2.6/urllib2.py, line 369, in _call_chain result = func(*args) File /System/Library/Frameworks/Python.framework/Versions /2.6/lib/python2.6/urllib2.py, line 518, in http_error_default raise HTTPError(req.get_full_url(), code, msg, hdrs, fp) urllib2.HTTPError: HTTP Error 404: Not Found I tried adding on r.stream.distance as well but I am getting a similar error. Any suggestions? Thank you, BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.stream.* addons issue
Thank you Ralf! BÜLENT On Oct 18, 2013, at 11:21 PM, Ralf Schäfer sena...@ecotoxicology.de wrote: Dear Bulent g.extension does not work on Mac to my knowledge. There is a detailed description of an alternative how to build modules on Mac using a few lines in the Shell. If you install GRASS from kyngchaos, he provides a description in a readme file. http://www.kyngchaos.com/software/grass I do not know whether this description is compatible with the installation from Michael Barton. Best regards Ralf ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.stream.* addons issue
Unfortunately, it still does not work. I downloaded the stable version of Grass from your website and I am having the same issue. Meanwhile, I am trying to remember an old trick I learned from a friend. I saved the .html file as a plain text file and dropped it into the Scripts folder in ContentsMacOS. I used to type a command starting with 'chmod' to make that plain text file a Unix executable file. It worked in the past for several addons. BÜLENT On Oct 19, 2013, at 1:39 AM, William Kyngesburye wokl...@kyngchaos.com wrote: Actually, I think g.extension should work now, I forgot to update the readme. I haven't tried it from the GUI, or installing system-wide. On Oct 18, 2013, at 3:36 PM, Bulent Arikan wrote: Thank you Ralf! BÜLENT On Oct 18, 2013, at 11:21 PM, Ralf Schäfer sena...@ecotoxicology.de wrote: Dear Bulent g.extension does not work on Mac to my knowledge. There is a detailed description of an alternative how to build modules on Mac using a few lines in the Shell. If you install GRASS from kyngchaos, he provides a description in a readme file. http://www.kyngchaos.com/software/grass I do not know whether this description is compatible with the installation from Michael Barton. Best regards Ralf ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user - William Kyngesburye kyngchaos*at*kyngchaos*dot*com http://www.kyngchaos.com/ Theory of the Universe There is a theory which states that if ever anyone discovers exactly what the universe is for and why it is here, it will instantly disappear and be replaced by something even more bizarrely inexplicable. There is another theory which states that this has already happened. -Hitchhiker's Guide to the Galaxy 2nd season intro ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.stream.* addons issue
Thanks William! BÜLENT On Oct 19, 2013, at 6:23 AM, William Kyngesburye wokl...@kyngchaos.com wrote: On Oct 18, 2013, at 9:50 PM, Bulent Arikan wrote: Unfortunately, it still does not work. I downloaded the stable version of Grass from your website and I am having the same issue. Like I said, I haven't tried it from the GUI. It's different code (python) than the CLI version (shell script), so maybe something isn't working right in python. Try running g.extension in the GRASS Terminal, not from the GUI. Works fine for me in the default GRASS_ADDON_PATH prefix (not global). It is not recommended to install globally on OS X, that's the GRASS.app. Mainly, it's just not the Mac way to add something to an existing application. Global installation on OS X should be the /Library/GRASS addons folder, but I don't think there is a default behavior in g.extension for that yet. You woud have to manually specify the prefix and make sure you are logged in as an administrator user (both the -s and -u flags force the prefix to be the GRASS app). Something like: g.extension ext=r.stream.basins prefix=/Library/GRASS/6.4/Modules Meanwhile, I am trying to remember an old trick I learned from a friend. I saved the .html file as a plain text file and dropped it into the Scripts folder in ContentsMacOS. I used to type a command starting with 'chmod' to make that plain text file a Unix executable file. It worked in the past for several addons. That sounds dangerous - making html executable, anything could happen if there was some valid script code in the html. BÜLENT On Oct 19, 2013, at 1:39 AM, William Kyngesburye wokl...@kyngchaos.com wrote: Actually, I think g.extension should work now, I forgot to update the readme. I haven't tried it from the GUI, or installing system-wide. On Oct 18, 2013, at 3:36 PM, Bulent Arikan wrote: Thank you Ralf! BÜLENT On Oct 18, 2013, at 11:21 PM, Ralf Schäfer sena...@ecotoxicology.de wrote: Dear Bulent g.extension does not work on Mac to my knowledge. There is a detailed description of an alternative how to build modules on Mac using a few lines in the Shell. If you install GRASS from kyngchaos, he provides a description in a readme file. http://www.kyngchaos.com/software/grass I do not know whether this description is compatible with the installation from Michael Barton. Best regards Ralf ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user - William Kyngesburye kyngchaos*at*kyngchaos*dot*com http://www.kyngchaos.com/ Theory of the Universe There is a theory which states that if ever anyone discovers exactly what the universe is for and why it is here, it will instantly disappear and be replaced by something even more bizarrely inexplicable. There is another theory which states that this has already happened. -Hitchhiker's Guide to the Galaxy 2nd season intro ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user - William Kyngesburye kyngchaos*at*kyngchaos*dot*com http://www.kyngchaos.com/ Mon Dieu! but they are all alike. Cheating, murdering, lying, fighting, and all for things that the beasts of the jungle would not deign to possess - money to purchase the effeminate pleasures of weaklings. And yet withal bound down by silly customs that make them slaves to their unhappy lot while firm in the belief that they be the lords of creation enjoying the only real pleasures of existence - the wisdom of Tarzan ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Wingrass Raster Layer Manager issue
Dear List, I am having trouble with the Layer Manager in winGRASS 6.4.3.1 downloaded from (http://grass.osgeo.org/grass64/binary/mswindows/native/). This is a PC running Windows 7. The problem is that when I want to add a raster or vector layer, that separate window that needs to pop up is not coming. Rather, the Layer Manager creates a layer with no name. I tried to use the rename option by right clicking. I types the name of the map but it dis not work. I also could not remove the layer. I downloaded and installed the same application package on a different PC with the same OS and I did not have the issue. This suggests to me that this is an issue about the computer I am using. I have not seen anything like that. Am I missing an additional file? Should I download and install .NET Framework (although, GRASS launches fine). Thank you, BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Wingrass Raster Layer Manager issue
Hi, Thank you for your response and suggestions. I figured it out soon after I posted the question. Although GRASS launched fine, it was missing .NET framework. I installed v 4.5 of that and it works now! Thank you again! BÜLENT On Oct 10, 2013, at 7:45 PM, Vaclav Petras wenzesl...@gmail.com wrote: Hi Bülent, can you check the Command console tab in the Layer Manager which is the second item from the left at the bottom of Layer Manager window? There might be an error message in there. You may also check the (black) windows command line. The layer with the label raster (double click to change properties) is created immediately when you press the Add raster map layer button. If there is no text, it is a problem. This seems that this will be hard to reproduce, so it would be great help if you double check all consoles, all opened windows, try this with nightly build of GRASS 7, do some screenshots (1)... it's hard to say what is important. For now it could be enough if you just check the Command console. Vaclav (1) You may want to send the screenshots only to me, to avoid file size limit check in the mailing list. On Thu, Oct 10, 2013 at 4:04 AM, Bulent Arikan bulent.ari...@gmail.com wrote: Dear List, I am having trouble with the Layer Manager in winGRASS 6.4.3.1 downloaded from (http://grass.osgeo.org/grass64/binary/mswindows/native/). This is a PC running Windows 7. The problem is that when I want to add a raster or vector layer, that separate window that needs to pop up is not coming. Rather, the Layer Manager creates a layer with no name. I tried to use the rename option by right clicking. I types the name of the map but it dis not work. I also could not remove the layer. I downloaded and installed the same application package on a different PC with the same OS and I did not have the issue. This suggests to me that this is an issue about the computer I am using. I have not seen anything like that. Am I missing an additional file? Should I download and install .NET Framework (although, GRASS launches fine). Thank you, BÜLENT ___ 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] v.transform issue
Dear List, I digitized a vegetation map and now I want to transform a centroid using 'v.transform' only using (x,y) values. I am getting this failure message on the Command Console of the Layer Manager: Command 'd.vect map=SSOOF_map_vtransform@Vegetation_30m' failed Details: Coor files of vector map SSOOF_map_vtransform@Vegetation_30m is larger than it should be (240 bytes excess) I am working in a UTM location at 30m res. This is a tiny centroid and the module runs without error message. I tried to transform the vector with and without -s flag (use transformation parameters instead of points). What is the problem? Thank you! BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] GIS Lock issue
Dear List, I am running GRASS 6.5 (2012) on Snow Leopard. I want to access a Mapset, have GRASS run analyses, and save new maps on an external HD. I use this external HD to back up the same computer I am running the GRASS. When I direct GRASS to browse the 'grassdata' folder on the external HD and start GUI from a Mapset (Permanent or others) I get this error: /Volumes/MacBook Pro Backup/Backups.backupdb/Bulent Arikan’s MacBook Pro/2011-10-25-113623/Ari Bub/Users/bulentarikan/grassdata/TR_WGS84_37N/PERMANENT/.gislock: Permission denied ERROR: /Applications/GRASS_6.5/GRASS-6.5.app/Contents/MacOS/etc/lock: Unable to properly access /Volumes/MacBook Pro Backup/Backups.backupdb/Bulent Arikan’s MacBook Pro/2011-10-25-113623/Ari Bub/Users/bulentarikan/grassdata/TR_WGS84_37N/PERMANENT/.gislock Please notify system personel. logout [Process completed] I am not running concurrent sessions and I know that GIS Lock is a feature but is there a way to get around this, so that I can save time by directly having some very large-size maps stored in external HD? Thank you. BÜLENT ARIKAN, PhD Senior Research Fellow Research Center for Anatolian Civilizations Koç University İstiklal Caddesi No: 181 Merkez Han Beyoğlu - ISTANBUL TURKEY 34433 (+ 90) 212-393-6036 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Grass 6.4.2 on Ubuntu 11.10
Thank you for your input. Yes, I misspoke: I should have said that I was looking for relatively bug-free version of GRASS for Linux. I must also mention that the version of Linux installed on computers is beyond my control since it is a lab maintained by others. The version the computers are running is Ubuntu 10.10 32 bit ( https://launchpad.net/+help-soyuz/ppa-sources-list.html ). The archaic looking GRASS 6.4 came from ( https://launchpad.net/~ubuntugis/+archive/ubuntugis-unstable ): they reached this link from the following address ( http://grass.fbk.eu/download/software.php#g64x ). I realize that this is a mirror site in Europe and I think that this last link has some problematic versions of GRASS. I had Windows version (6.4) installed from this address for my class last semester and there were a lot of problems with mapcalc, gregion, and v.in.ogr. I will try finding the link for GRASS 7, which had problems so that it can be fixed. As for Linux GRASS 6.4 issue: I only typed GRASS in Terminal and GUI came up however, I will try what you suggested (grass64 -gui). Meanwhile, can anyone suggest a link for GRASS 6.4 for the version of Ubuntu I have to work with? Thanks, Bulent On Mar 5, 2012, at 2:32 PM, Markus Neteler wrote: On Mon, Mar 5, 2012 at 12:59 PM, barikan bulent.ari...@gmail.com wrote: Hi, As a follow-up on this issue, and as a person who never installed GRASS on Linux, I wanted to ask the Forum if there is a version of GRASS 6.x that will run on 32-bit system, which has a similar look to the most recent Unix interface. If you install a Linux distro package, it should come in either 64 or 32 bit depending on your system. Unfortunately some distros are quite lagging behind with the packaging. I installed 6.4.2 from the US Tier#2 web site of GRASS (the link gave the date of 2009) Can you please indicate the link? Anything which should be fixed? but I am sure that the interface looked archaic. Probably you did not start with -gui parameter and GRASS used an old (tcl) setting from a previous run of GRASS on that machine. Try (later it will remember): grass64 -gui It is important that the version is free of bugs since it will be used for a workshop designed for beginners. I will appreciate your suggestions. Say minimized number of bugs - bugfree software hardly exists. GRASS 6.4.2 should be very good for teaching, I used in happily in a course recently. * I also installed GRASS 7 but there are issues with vector import and basic analyses. Please indicate which issues you found (otherwise we cannot fix them). Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to import as vector?
Thank you for the suggestions! Yes, it will be a lot of work and time consuming to digitize it. I was also suggested to use Inkscape to autotrace it and save it in (svg). However, saving it as (svg) does not help because 'v.in.ogr' does not seem to import that format. I even saved it as (dxf) but the import failed –nothing showed up after importing it as vector layer– with 'v.in.ogr' and 'v'in.dxf'. Is there a separate module for importing (svg) files or should I use a different file extension? Thank you, Bulent On Thu, Jan 12, 2012 at 1:49 AM, Markus Neteler nete...@osgeo.org wrote: On Wed, Jan 11, 2012 at 1:25 PM, Bulent Arikan bulent.ari...@gmail.com wrote: Dear List, I found a plan / architectural drawing online. I am not sure about the resolution but JPG version is 25 Mb and TIFF version is 410 Kb. I want to import the drawing as vector into 'xy location' and then rectify it. I already did rectify the image as raster but I will need a vector layer. Looking at 'v.in.ogr' I am not sure how to bring images of such formats into GRASS as vector. Raster to vector transformation does not work. Since it is a raster map you have (at least) two options: - use it as background for a digitizer (likely a lot of work) - try image processing to get out raster lines, then r.thin and subsequently r.to.vect with lines. This much depends on the JPG quality. - (there are some autotracer tools, too) Markus -- BÜLENT ARIKAN, PhD Senior Research Fellow Research Center for Anatolian Civilizations Koç University İstiklal Caddesi No: 181 Merkez Han Beyoğlu - ISTANBUL TURKEY 34433 (+ 90) 212-393-6036 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Connecting table
OK! That makes sense now. I have one last standing issue about this. It seems that the second layer I created using 'v.db.addtable' (after setting the dbf to SQL and importing a centroid –vector– map) has a column for CAT but it does not have a number. I checked the book to make sure and it says it should have a number and it does not have to be continuous from the first layer. I tried assigning a CAT# using 'v.category' and 'v.db.connect' but I still do not see a number in CAT column. I think I should have a number there so that my ...WHERE cat = argument makes sense. Any ideas how to correct this? Thank you, Bulent On Thu, Jan 5, 2012 at 12:27 PM, Micha Silver mi...@arava.co.il wrote: ** On 01/05/2012 10:49 AM, Bulent Arikan wrote: Hi, I really would like to learn more about the SQL command-based operations in GRASS. To summarize: I have created a new mapset where the driver is set to SQL. I have created second layers for all of my centroids (the tables in the second layers have a different name like Db_runivar ). I used 'v.db.addtable' to have GRASS create columns to upload the results from extended statistics of r.univar text files (e.g., Db_runivar.txt) for each centroid. Now, I have been trying to get 'db.execute' to work. One suggestion was that I use a loop for automated updating (I am assuming that I have to type this in the Terminal): while read l; do \ col=`echo $l | cut -d= -f1`; val=`echo $l | cut -d= -f2`; \ sqlite3 sqlite.db UPDATE centroid SET ${col}=${val} ;; \ done univar.txt I came up with a statement (below) by looking at the manual in 'r.univar' to create an SQL command file: sed -e '1d' Dogubayazit_runivar.txt | awk -F'l' '{print UPDATE Dogubayazit_runivar SET non_null_cells = $2, null_cells = $3, min = $4, max = $5, range = $6, mean = $7, mean_of_abs = $8, stddev = $9, variance = $10, coeff_var = $11, sum = $12, sum_abs = $13, first_quart = $14, median = $15, third_quart = $16, perc_90 = $17 WHERE cat = $1;}' The above will create a list of update statements, one for each row in the runivar.txt file. You can't feed db.execute a list of statements, only one at a time. So you'll have to dump the output of the above command into an intermediary file, then use the db.execute input parameter to specify that file of SQL commands. It might be worth noting that if you have many UPDATE's, then doing them one by one with an individual call to the database each time might take a bit of time. In this case, going back to sqlite, you could wrap all the UPDATE statements between a BEGIN...COMMIT clause, then feed that file to db.execute (or directly to sqlite), and it will open one connection, and push all the updates at once. For many 1000's of rows this would probably be much faster. that can be used in 'db.execute' but I am getting a syntax error: DBMI-SQLite driver error: Error in sqlite3_prepare(): near sed: syntax error ERROR: Error while executing: 'sed -e '1d' What seems to be the problem? Thank you for your time. Bulent On Wed, Jan 4, 2012 at 9:16 AM, Micha Silver mi...@arava.co.il wrote: On 01/03/2012 10:44 PM, Bulent Arikan wrote: Thank you! This seems to be an interesting solution to the problem especially because it involves changing the driver from DBF to SQLITE. I ended up using 'v.rast.stats', which adds columns for extended statistics and uploads values all at once. The only thing is: data are recorded as part of the same layer (i.e., Layer 1). So, I cheated! Yes, v.rast.stats is definitely the way to go if you want raster univariate statistics pushed into a polygon vector. However, I started trying your method. I created a new mapset to use sqlite and copied some vector maps from a mapset where default driver is DBF. I defined SQLITE as the new driver in the new mapset using 'db.connect' (driver=sqlite, database='$GISDBASE/$LOCATION_NAME/$MAPSET/sqlite.db', no flags checked). Then I wanted to use 'v.db.connect' for a vector map but I cannot select table name etc. I also tried just using 'v.db.connect' but I still cannot see anything under Table. So. I am confused about how to define a new driver in a mapset. I do as follows: eval `g.gisenv` (This creates the environment variables for GISDBASE, etc.) Now: db.connect driv=sqlite database=$GISDBASE/$LOCATION_NAME/$MAPSET/sqlite.db Next: g.copy vect=old,new to create a new copy of the original vector. The original will still have its attrib table as dbf. THe new copy will have an sqlite based attribute table. Now you can do: v.db.addcol new_vector col=..., ... . BTW, I think that v.rast.stats will automatically create the needed columns . Cheers, Micha Thank you again, Bulent On Tue, Jan 3, 2012 at 9:13 PM, Micha Silver mi...@arava.co.il wrote: On 01/03/2012 04:50 PM, Bulent Arikan wrote: Dear List, I rasterized a centroid and ran 'r.univar', whose extended
Re: [GRASS-user] Connecting table
Hi, I really would like to learn more about the SQL command-based operations in GRASS. To summarize: I have created a new mapset where the driver is set to SQL. I have created second layers for all of my centroids (the tables in the second layers have a different name like Db_runivar ). I used 'v.db.addtable' to have GRASS create columns to upload the results from extended statistics of r.univar text files (e.g., Db_runivar.txt) for each centroid. Now, I have been trying to get 'db.execute' to work. One suggestion was that I use a loop for automated updating (I am assuming that I have to type this in the Terminal): while read l; do \ col=`echo $l | cut -d= -f1`; val=`echo $l | cut -d= -f2`; \ sqlite3 sqlite.db UPDATE centroid SET ${col}=${val} ;; \ done univar.txt I came up with a statement (below) by looking at the manual in 'r.univar' to create an SQL command file: sed -e '1d' Dogubayazit_runivar.txt | awk -F'l' '{print UPDATE Dogubayazit_runivar SET non_null_cells = $2, null_cells = $3, min = $4, max = $5, range = $6, mean = $7, mean_of_abs = $8, stddev = $9, variance = $10, coeff_var = $11, sum = $12, sum_abs = $13, first_quart = $14, median = $15, third_quart = $16, perc_90 = $17 WHERE cat = $1;}' that can be used in 'db.execute' but I am getting a syntax error: DBMI-SQLite driver error: Error in sqlite3_prepare(): near sed: syntax error ERROR: Error while executing: 'sed -e '1d' What seems to be the problem? Thank you for your time. Bulent On Wed, Jan 4, 2012 at 9:16 AM, Micha Silver mi...@arava.co.il wrote: ** On 01/03/2012 10:44 PM, Bulent Arikan wrote: Thank you! This seems to be an interesting solution to the problem especially because it involves changing the driver from DBF to SQLITE. I ended up using 'v.rast.stats', which adds columns for extended statistics and uploads values all at once. The only thing is: data are recorded as part of the same layer (i.e., Layer 1). So, I cheated! Yes, v.rast.stats is definitely the way to go if you want raster univariate statistics pushed into a polygon vector. However, I started trying your method. I created a new mapset to use sqlite and copied some vector maps from a mapset where default driver is DBF. I defined SQLITE as the new driver in the new mapset using 'db.connect' (driver=sqlite, database='$GISDBASE/$LOCATION_NAME/$MAPSET/sqlite.db', no flags checked). Then I wanted to use 'v.db.connect' for a vector map but I cannot select table name etc. I also tried just using 'v.db.connect' but I still cannot see anything under Table. So. I am confused about how to define a new driver in a mapset. I do as follows: eval `g.gisenv` (This creates the environment variables for GISDBASE, etc.) Now: db.connect driv=sqlite database=$GISDBASE/$LOCATION_NAME/$MAPSET/sqlite.db Next: g.copy vect=old,new to create a new copy of the original vector. The original will still have its attrib table as dbf. THe new copy will have an sqlite based attribute table. Now you can do: v.db.addcol new_vector col=..., ... . BTW, I think that v.rast.stats will automatically create the needed columns . Cheers, Micha Thank you again, Bulent On Tue, Jan 3, 2012 at 9:13 PM, Micha Silver mi...@arava.co.il wrote: On 01/03/2012 04:50 PM, Bulent Arikan wrote: Dear List, I rasterized a centroid and ran 'r.univar', whose extended stats are saved as a text file. Then, I ran 'v.db.addtable' to create a second layer in the attribute table of the centroid and I had columns added in this second layer using 'v.db.addcol'. I want GRASS to upload values from the text file but I am not sure how to perform this task. Is there a shorter way of creating a second layer and uploading values from the text file or what should I do next so that the columns in the second layer will be populated using the text file I have? I'm not sure about a shorter way, but here's an option: First setup your mapset to save attributes to sqlite. v.db.connect centroid driver=sqlite database= Now add the columns which will accept r.univar values to the sqlite table: v.db.addcol centroid col=n double, null_cells double, cells double, min double, max double, range double, mean double, mean_of_abs double, stddev double, variance double, coeff_var double, sum double Run r.univar on your raster, putting results into a text file r.univar your_rast -g univar.txt Now do this loop to update values for the centroid: while read l; do \ col=`echo $l | cut -d= -f1`; val=`echo $l | cut -d= -f2`; \ sqlite3 sqlite.db UPDATE centroid SET ${col}=${val} ;; \ done univar.txt Assuming you want to run this for several centroids/areas, you'll probably want to add a WHERE clause to the UPDATE statement so as to put values for only one certain row (centroid). HTH, Micha Thank you for your time. GRASS 6.5 svn on Snow Leopard -- BÜLENT This mail was received via Mail-SeCure System
[GRASS-user] Connecting table
Dear List, I rasterized a centroid and ran 'r.univar', whose extended stats are saved as a text file. Then, I ran 'v.db.addtable' to create a second layer in the attribute table of the centroid and I had columns added in this second layer using 'v.db.addcol'. I want GRASS to upload values from the text file but I am not sure how to perform this task. Is there a shorter way of creating a second layer and uploading values from the text file or what should I do next so that the columns in the second layer will be populated using the text file I have? Thank you for your time. GRASS 6.5 svn on Snow Leopard -- BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Connecting table
Thank you! This seems to be an interesting solution to the problem especially because it involves changing the driver from DBF to SQLITE. I ended up using 'v.rast.stats', which adds columns for extended statistics and uploads values all at once. The only thing is: data are recorded as part of the same layer (i.e., Layer 1). So, I cheated! However, I started trying your method. I created a new mapset to use sqlite and copied some vector maps from a mapset where default driver is DBF. I defined SQLITE as the new driver in the new mapset using 'db.connect' (driver=sqlite, database='$GISDBASE/$LOCATION_NAME/$MAPSET/sqlite.db', no flags checked). Then I wanted to use 'v.db.connect' for a vector map but I cannot select table name etc. I also tried just using 'v.db.connect' but I still cannot see anything under Table. So. I am confused about how to define a new driver in a mapset. Thank you again, Bulent On Tue, Jan 3, 2012 at 9:13 PM, Micha Silver mi...@arava.co.il wrote: ** On 01/03/2012 04:50 PM, Bulent Arikan wrote: Dear List, I rasterized a centroid and ran 'r.univar', whose extended stats are saved as a text file. Then, I ran 'v.db.addtable' to create a second layer in the attribute table of the centroid and I had columns added in this second layer using 'v.db.addcol'. I want GRASS to upload values from the text file but I am not sure how to perform this task. Is there a shorter way of creating a second layer and uploading values from the text file or what should I do next so that the columns in the second layer will be populated using the text file I have? I'm not sure about a shorter way, but here's an option: First setup your mapset to save attributes to sqlite. v.db.connect centroid driver=sqlite database= Now add the columns which will accept r.univar values to the sqlite table: v.db.addcol centroid col=n double, null_cells double, cells double, min double, max double, range double, mean double, mean_of_abs double, stddev double, variance double, coeff_var double, sum double Run r.univar on your raster, putting results into a text file r.univar your_rast -g univar.txt Now do this loop to update values for the centroid: while read l; do \ col=`echo $l | cut -d= -f1`; val=`echo $l | cut -d= -f2`; \ sqlite3 sqlite.db UPDATE centroid SET ${col}=${val} ;; \ done univar.txt Assuming you want to run this for several centroids/areas, you'll probably want to add a WHERE clause to the UPDATE statement so as to put values for only one certain row (centroid). HTH, Micha Thank you for your time. GRASS 6.5 svn on Snow Leopard -- BÜLENT This mail was received via Mail-SeCure System. ___ grass-user mailing listgrass-user@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver GIS Consultant, Arava Development Co.http://www.surfaces.co.il -- BÜLENT ARIKAN, PhD Senior Research Fellow Research Center for Anatolian Civilizations Koç University İstiklal Caddesi No: 181 Merkez Han Beyoğlu - ISTANBUL TURKEY 34433 (+ 90) 212-393-6036 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] 'v.distance'
Dear List, I have a points vector file, which I use as input in 'v.distance' to create a new vector file of lines (showing the minimum distance to the coast). The output is naturally a vector file that contains lines between each point and the sea. However, I realize that the output map cannot be queried. Even if I have GRASS connect the new (lines) vector file with the points as the second layer, querying of the lines does not return results. What I want to do is to assign CAT numbers to each line (it can be same with the CAT# of its point of origin) so that I can display lines individually in the new lines vector file for further analyses. Thank you for your help! Bulent ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Installing GRASS on Windows
Dear List, I have a student who needs to install GRASS preferably on his laptop (running Vista) or desktop (running XP). He has troubles with installation. Are there additional programs that need to be installed? What is the best resource to find the (.exe) files? Thanks, -- BÜLENT ARIKAN, PhD Senior Research Fellow Research Center for Anatolian Civilizations Koç University İstiklal Caddesi No: 181 Merkez Han Beyoğlu - ISTANBUL TURKEY 34433 (+ 90) 212-393-6036 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] vector database connection
Dear List, I prepared a vector map using v.distance, showing the minimum distance from a series of vector points (i.e., cities) to a coastline (another vector map of area type). Naturally, the minimum distances are shown as lines in this new vector map. When I want to see the attribute table, I get the following warning Database connection for vector map is not defined in DB file. Define new connection in manage layers tab. However, I am not sure how to define such connection since there is not a (.dbf) for this new vector map in the folder. Do I connect this new vector map of lines to the original vector map of points (cities) as a second layer? If so, how does one do that? Eventually, I want to get an elevation profile using these minimum distance lines. I think the most appropriate module for this is r.profile (which means that I have to rasterize the vector lines) but if you have other suggestions I will be grateful. Thank you, -- BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] making fcell raster maps using climate data
Dear List, I have an Excel file that contains precipitation data in 13 columns (mean values for each month and annual average) from a region between certain years (in rows): eventually showing temporal changes in precipitation values. The data are taken at a weather station whose coordinates are known. I want to use this file in order to make an fcell raster map in GRASS. After some research I found few leads and mostly on interpolation of data. However, I am confused about the actual process turning the data in Excel into a raster map. I will greatly appreciate any suggestions as to which module to use and how to arrange the data. Thank you, -- BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.surf.rst question
Dear List, I have 30m ASTER DEM from a region (UTM-37) which I am trying to interpolate using ' v.surf.rst ' to obtain 10m resolution DEMs. I prefer this module since I will also need additional maps such as slope, aspect, prof curv, etc. I will try describing my problem below and I will appreciate any suggestions/ideas. Since the original 30m res. image covers a large area, I cropped it into several tiles using 'r.mapcalc'. Each tile has about 12 million cells in 30m res. and 115 million cells in 10m res. and takes an average of 4 hrs to complete one run of v.surf.rst. I sample each tile using 'r.random' at 50%. Then I change the region settings in 'g.region': I align region to resolution and set the resolution to 10m before running v.surf.rst http://v.surf.rs. I ran the module at different tension settings (55., 40., 25. and 10.) while keeping smoothing parameter (.1) and other parameters at default settings. When I add legend to my 10m DEM from the toolbar in Map Display, the elevation range is usually a couple meters above or below the original map. After each run, I set my region resolution back to 30m, and using r.mapcalc, I prepare a difference map [ (30m DEM) - (10m DEM) ]. In this map, I get extremely high and low values (i.e., the range in the difference map can be -70 to 183 suggesting to me that there is both over- and undershooting), which seem to be limited to mountain tops and some canyons. When I run 'r.report', I can see that actually these extreme values are limited to very few cells (about 3-5% of 12 million cells) and the main variation for 95% of cells is within more reasonable range such as +/- 3 meters. I recently became aware of some accuracy issues related to ASTER DEMs especially in mountainous areas however, I have the same issue for flat parts of this landscape as well. My question is: How can I fix the issue of having these few cells with extreme values? I am thinking that: - I may be having this problem because I am trying to interpolate images that are too large (crop tiles of smaller areas?), - I may be able to fix the negative cells by making them NULL and then using ' r.fill.nulls' module but even if this is advisable, I do not know of a solution for higher than normal value cells. Again, I appreciate your ideas and suggestions. Thank you, -- BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] interpolating images with coastal areas
Dear List, I am working on interpolating 30m ASTER DEMs to 10m resolution. I am using GRASS 6.4.1 svn on Mac OS 10.6.6. Since the imagery covers a large area, I crop the 30m DEM into individual tiles in order to make the project more manageable. The interpolation process (using v.surf.rst) has been working fine and the values in resultant 10m DEMs have been quite accurate, especially for inlands. However, I have some problems with interpolating tiles where coastline is present. In these tiles, the range of values in 10m DEMs go below zero. In some cases I have values such as -21 meters whereas the original –30m– tiles have only zero. Obviously, I am concerned about the accuracy of the end product. Do I get high negative values because of the cells with zero value? Should I zoom further in and avoid having zero-value cells? (which is a problem with coastlines where they are not straight) Does any one have a suggestion about how to deal with this problem? * I sample the original imagery at 50% in r.random, * I use the default parameter values in v.surf.rst (tension= 40.) and use the 30m DEM as MASKMAP, * I align region to resolution, using g.region. Thanks! -- BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] UTM location
Dear List, I have several ASTER GDEM (originally in Latlong, 30m res) images, which are reprojected into corresponding UTM locations. My research area covers two zones: UTM 37 and 38. Is there a way of creating a UTM location that can accept data from both zones? I think I saw a thread in recent list serve emails but I could not find it. I apologize if my question has recently been answered. I -- BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] patching
Dear List, I have several ASTER GDEM (originally Latlong) images that I want to reinterpolate from 30 to 10 m res. I also need to patch these images. I am not sure whether: - I should patch them in Latlong using r.series, reproject into UTM and interpolate them using r.resamp.rst or - reproject images individually into UTM, interpolate individually using r.resamp.rst and then patch them using r.series. So far, I tried the second option. After cubic reprojection and interpolation to 10m the tiles have little borders of NULL areas when I display them together. They don't go away when patched. I will appreciate any advice Thanks. -- BÜLENT ARIKAN, PhD School of Human Evolution and Social Change Arizona State University Tempe - AZ 85287-2402 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.in.ascii problem
Dear List, I have a tab-delimited text (.txt) file and I am experimenting with importing points into GRASS (6.5 svn) in Snow Leopard. Here is the content of the file: 675166.97 4450516.81 246 Orhan 675135.87 4450516.98 250 Osman 675156.68 4450501.48 250 Sehitlik 675155.74 4450560.93 243 SaatKulesi 675192.48 4450504.39 245 Kuledenilerisi The first two columns are UTM (zone 36) east and north readings, the third is the elevation and the last column is the explanation for each point. In 'v.in.ascii' module, I specified point as the format and typed tab (I also tried \t). I left x and y column numbers at default (1 for x, 2 for y), I did not specify z or column definition. Here is the error message i am getting: Scanning input for column types... Maximum input row length: 9 Maximum number of columns: 1 Minimum number of columns: 1 ERROR: y column number minimum last column number (incorrect field separator?) * I tried importing this file without decimals at the end of each coordinate, and used ' | ' instead of tab as delimiter. The result is the same. What am I doing wrong? Thanks, -- BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] EPSG code for KML
Dear List, I occasionally use Google Earth –GE– (kml files) and I am aware of the fact that there is a bit of distortion (i.e., metric offset) between GE and other projections for a variety of reasons. I am wondering if anyone knows whether there is an EPSG code for KML format so that I can create a LOCATION in Grass using this code and avoid/minimize such distortions when I import vector files. Thanks, -- BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.in.gpsbabel
Hi, I am running 6.5 svn on Snow Leopard. I am trying to import points and tracks from a Garmin device (connected via USB). Using GPS Babel (v 1.4.1 for Mac), I prepared several files (.txt, .gdb, and .gpi) and tried to import these using v.in.gpsbabel (UTM wgs84 location) –k flag checked–. I am getting this error message: ERROR: The gpsbabel program was not found, please install it first. http://gpsbabel.sourceforge.net As silly as it sounds, I am not sure if GPS Babel needs installing. After the download, a disk image appeared and double clicking launched the program (i.e., the window that retrieves data from the device and writes files). I did not install anything because it launched and it seems to be working (i.e., it writes files). I may be skipping some steps, so I will appreciate any suggestions. Thanks, -- BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.in.gpsbabel-2
Hi, Hamish, thanks for your suggestion about fixing the GPS Babel application issue. It worked and now I can run the module. I am back with another issue, again ' v.in.gpsbabel' module! My GPS unit records data in UTM projection with WGS84 datum. I just attempted to import the waypoints into GRASS using a Garmin MapSource text (tab delimited) file. I did not check -k flag (do not transform from WGS84). I kept the format box in the Required tab in default and provided PROJ.4 parameter for UTM [ +proj=utm +zone=36 +ellps=WGS84 +datum=WGS84 +units=m +no_defs ] in the Proj box of the Required tab, since this module assumes lat/long projection with WGS84 datum. Here are the error messages I receive: Loading Waypoints from /Users/barikan2/Grmn.txt... [ERROR] SERIAL: tcgetattr error: Inappropriate ioctl for device [ERROR] Cannot open serial port '/Users/barikan2/Grmn.txt' GARMIN:Can't init /Users/barikan2/Grmn.txt /Users/barikan2/Desktop/GRASS-6.5.app/Contents/MacOS/scripts /v.in.gpsbabel: line 263: /Users/barikan2/grassdata/TR_WGS84_36N/PERMANENT/.tmp /Bulent-Arikans-iMac.local/1470.0.xcsv: No such file or directory /Users/barikan2/Desktop/GRASS-6.5.app/Contents/MacOS/scripts /v.in.gpsbabel: line 263: [: : integer expression expected ERROR: Error loading data from gpsbabel I used few other Garmin file formats such as Garmin Points of Interest (.gpi) and Garmin MapSource (.mps) and the result does not change. I assume I am doing something wrong with the Projection detail! I will appreciate suggestions/advice. Thanks, -- BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.in.gpsbabel-3
Thanks for all the replies. It seems that I forgot to specify the format ( _txt) in the Required tab. So, it fixed the problem and now I can import all the waypoints and tracks. However, I am puzzled with one detail. As I mentioned: - I set my GPS unit to record measurements in UTM, - the text file prepared by GPS Babel has header showing UTM as projection, - even if I specify the projection parameter according to PROJ.4 in the Required tab of the module, the attribute tables for the resultant waypoint and track vector files are in lat/long format. How can this happen? I never work with lat/long. BÜLENT ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.proj or r.resamp.rst issues
Thanks Glynn. That makes sense now. I did not know about the inherent problem of overshooting in cubic interpolation. ASTER DEMs I have seem to be pretty consistent in terms of data ranges and representation of values (based on r.report), so I am not sure about the noise effect in my case. However, does this mean that I am sort of out of corrective options if I want to use cubic reprojection of imagery from Latlong to UTM? I mean I either use nearest method for error-free reprojection or I have to be careful about the results of cubic method. Finally, when you say use higher values for tension parameter, I guess you mean values lower than default 40. such as 30. or 20. Bulent On Fri, May 7, 2010 at 5:08 AM, Glynn Clements gl...@gclements.plus.comwrote: Bulent Arikan wrote: I am running GRASS 6.5 svn (Snow Leopard). I have several ASTER GDEMs (Latlong, 30m res.), which I reprojected into UTM using both 'nearest' and 'cubic' methods ('r.proj'). Only in some imagery that are reprojected in cubic, I ended up having couple of cells (literally, 1-2 cells out of 8 million in average) with minus (-) values. For example, in a DEM where the elevation values are between 800-2600 meters, I have cell values between -150 and -85 meters. This does not seem to be an issue in reprojected imagery with the nearest method. I am not sure how these minus values are introduced at the first place. Cubic interpolation can introduce overshoot, as can other forms of spline interpolation. Linear and nearest-neighbor interpolation don't have this issue. With r.resamp.rst, the problem can be alleviated to a degree by using higher values for the tension= parameter. Also, if your data is noisy, this will tend to exaggerate the gradients, making overshoot more likely. Filtering the data first will reduce the errors. -- Glynn Clements gl...@gclements.plus.com -- BÜLENT ARIKAN School of Human Evolution and Social Change Arizona State University Tempe - AZ 85287-2402 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.proj or r.resamp.rst issues
Hi, I am running GRASS 6.5 svn (Snow Leopard). I have several ASTER GDEMs (Latlong, 30m res.), which I reprojected into UTM using both 'nearest' and 'cubic' methods ('r.proj'). Only in some imagery that are reprojected in cubic, I ended up having couple of cells (literally, 1-2 cells out of 8 million in average) with minus (-) values. For example, in a DEM where the elevation values are between 800-2600 meters, I have cell values between -150 and -85 meters. This does not seem to be an issue in reprojected imagery with the nearest method. I am not sure how these minus values are introduced at the first place. Since I interpolate GDEMs from 30m to 10m using 'r.resamp.rst' (keeping tension at 40.), I need to fix this problem. I thought of making negative value cells null but this may create more problems during interpolation (actually, I tried and r.resamp.rst could not complete interpolation). If I take absolute values of cells to make negative cells positive, then I introduce cells with unrealistic elevation data. It seems like after interpolation –without any correction of negative values–, the number of negative value cells do not increase, they are still less than 1% of total cells. I will appreciate any suggestions why cubic method produces negative values in some imagery and how to fix this before or after interpolation. * I did not try changing default r.resamp.rst values, such as tension, because most of the imagery I interpolated have very accurate results; the elevation difference between interpolated (10m res) and original (30m res) is around 1 meter. Thank you for your comments. -- BÜLENT ARIKAN School of Human Evolution and Social Change Arizona State University Tempe - AZ 85287-2402 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.proj and r.resamp.rst issues -2
Hi, I might have found a solution to my problem, which I shared with you earlier today. I am writing about it in case someone else has a similar issue and this may be a way of solving the problem. However, I am not certain if this is the easiest or the most practical way of getting the results. The problem was the presence of few negative value cells in 'cubic' reprojected DEMs (from Latlong to UTM). I ran r.mapcalc to make these negative cells NULLs. Then, I ran 'r.fillnulls' in that part of the DEM and patched the 'filled' map with the original map (the one that has NULLs). Finally, I interpolated the resultant (i.e., patched) imagery from 30m to 10m. The values in 10m DEM look accurate and free of negative value cells. If you know another way of doing it and share with me, I will greatly appreciate. -- BÜLENT ARIKAN School of Human Evolution and Social Change Arizona State University Tempe - AZ 85287-2402 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] 'v.surf.rst' VS. 'r.resamp.rst'
Hi, I'd like to thank everyone for suggestions on my previous question about the interpolation of DEMs from 30 to 10 m res. and problems I had with 'v.surf.rst' module. I am still experimenting interpolation with it and i hope to get a good interpolated DEM finally. Hamish suggested I use 'r.resamp.rst' instead and I did (with the default values). The results look really impressive; it took equally long time for a single DEM (about 90 minutes). I used the resultant DEM to make slope and aspect maps; nothing seems weird. I also used the mapcalc module to check the difference between 30m and 10m DEMs and it is within -1 and 1 m range. However, I used 'r.resamp.rst' on one imagery only. So, I am writing to ask for suggestions about which module to use for 30m to 10m interpolation. I understand that 'v.surf.rst' and 'r.resamp.rst' use the same approach (i.e., spline tension) but in different ways (which is not surprising since GRASS offers a lot of alternative ways of manipulating data) and given my limited experience with GIS, I will appreciate suggestions from the more experienced users. Thanks, * I am running a slightly old version (late 2009) of GRASS 6.5 on Mac OS 10.5.8 -- BÜLENT ARIKAN School of Human Evolution and Social Change Arizona State University Tempe - AZ 85287-2402 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] advice on interpolation
Hi, I have an ASTER GDEM (Latlong, 30 m resolution, from NASA web site) imagery. I want to produce 10m res DEM from it. Here is what I have done and what my problem is. I really appreciate suggestions. - I reprojected the imagery into a UTM location (using the 'nearest' method), - Using 'g.region', I set my region to match the original DEM - I set my computational extents from the display, - then, I sampled 50% of the original imagery (in r.random) and created a vector map (it took long time to do this and it has around 5 million points in it), - I set the resolution to 10 meters using 'g.region', - then, using 'v.surf.rst', I tried to make a new DEM -10m resol- from the sampled vector map and masking the operation with the original -30m- DEM. I do not change the default settings and I use a smoothing parameter of '.1' and i use the value column from the vector map for approximation in the process. I get a lot of error messages during interpolation, like; G_ludcmp() failed! n=0 taking too long to find points for interpolation--please change the region to area where your points are. Continuing calculations... I am not sure but may be my region settings is wrong (i.e., too wide) although I used the imagery to set them. I never waited for interpolation to run completely. Thank you for suggestions! -- BÜLENT ARIKAN School of Human Evolution and Social Change Arizona State University Tempe - AZ 85287-2402 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] help
Hi, I have been using r.param.scale and inverting DEM method for identifying peaks. I am specifically interested in finding the high spots on a landscape; not just the highest single cell. It has been suggested that r.prominence may be of help. I realise that this comes in a file that needs to be compiled (.c extension). Before compiling, I will appreciate any insight on what it does. Thank you -- BÜLENT ARIKAN School of Human Evolution and Social Change Arizona State University Tempe - AZ 85287-2402 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.prominence
Hi, I have been using r.param.scale and inverting DEM method for identifying peaks. I am specifically interested in finding the high spots on a landscape; not just the highest single cell. It has been suggested that r.prominence may be of help. I realise that this comes in a file that needs to be compiled (.c extension). Before compiling, I will appreciate any insight on what it does. Thank you -- BÜLENT ARIKAN School of Human Evolution and Social Change Arizona State University Tempe - AZ 85287-2402 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] inverting DEMs
Hi, I am trying to invert a DEM (i.e., make peaks basins and vice versa) but there does not seem to be a module for this. I am assuming this can be done through ' r.mapcalc ' which I have very limited experience. I also want to get ' r.prominence ' from the Add-ons site but the text file is not downloadable. Should I copy and paste to create a text file that contains the script? Thank you for all the suggestions, -- BÜLENT ARIKAN School of Human Evolution and Social Change Arizona State University Tempe - AZ 85287-2402 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] ' v.in.db ' question
Hi, I am trying to import an OpenOffice dbf into a wgs84 PERMANENT Location. I did this several times in the past with older versions of GRASS. I enter all the necessary information (input table, column names for x, y, z and CAT#). The Connection tab shows dbf for the driver name and database name reads $GISDBASE/$LOCATION_NAME/$MAPSET/ dbf/ I get this message: DBMI-DBF driver error: Table 'Iron_II_KNOWN_Sites' doesn't exist Missing column CAT in table Iron_II_KNOWN_Sites I will greatly appreciate any suggestions. BÜLENT ARIKAN SHESC-ASU___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user