Re: [GRASS-user] Permission denied to access a mapset in Linux
Dear Glynn, When you said: GRASS won't allow you to select a mapset as the current mapset unless you are the owner of the mapset directory. How can I change the owner of a mapset directory? Regards, Antonio On 12-11-2013 14:46, Glynn Clements wrote: António Rocha wrote: I have a Location and Mapset in a folder in Linux Server that was created by me and I runned some scripts and created some maps etc. I have a colleague that is trying enter in the mapset and it gets Permission Denied. How can I allow her to access the mapset without any problems? You can access any mapset provided that you have the relevant permissions (read permission for files, read and execute permission for directories). GRASS won't allow you to select a mapset as the current mapset unless you are the owner of the mapset directory. If you just need to access maps from a particular mapset, you can add the mapset to the search path with g.mapset, or reference maps using fully-qualified names (map@mapset format). ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Permission denied to access a mapset in Linux
Dear all, I have a Location and Mapset in a folder in Linux Server that was created by me and I runned some scripts and created some maps etc. I have a colleague that is trying enter in the mapset and it gets Permission Denied. How can I allow her to access the mapset without any problems? Regards, Antonio ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Exporting a map into JPEG or PNG- get black image
Dear all, I have a Raster layer with colors and values 1 till 8 and null values. I want to export as JPEG using r.out.gdal but I get a totally black image. Any tips on how to solve this? Regards, Antonio ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] List rasters with valid data inside a defined computational region
Dear all, I would like to know if there is any direct method (of a suggestion of an indirect method) to list raster layers with valid pixels (not null) in a defined computational region? Thank you Regards, Antonio Rocha ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] List rasters with valid data inside a defined computational region
Hello Nikos I would rephrase and say list of raster maps with at least a single not-NULL cell inside a computational region. Regards, Antonio On 25-10-2013 11:26, Nikos Alexandris wrote: On Friday 25 of October 2013 10:07:27 António Rocha wrote: I would like to know if there is any direct method (of a suggestion of an indirect method) to list raster layers with valid pixels (not null) in a defined computational region? You mean only raster maps that do not contain NULL cells at all, inside the computational region? Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Usage of r.random with Land cover maps
dear all, I have a land cover map (values from 1 till 15) and I want to extract randomly 100 pixels from each class. I tried to use the cover option but it only produced an additional column with values from the cover map. My question is: how can I extract exactly 100 pixels from each class without having to create a specific cover map for each class? Thank you Regards, Antonio ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] access denied to MAPSET in linux
Dear all, I'm tryging to access a Location/mapset located in a folder (where I have write/execute/read permission) but the mapset and location were not created by me. When I run grass I get permission denied. How can I get access or how can the original user of the mapset can give access permission to other users? This is all in a Linux server Thanks Antonio ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Get output from r.stats in Python script
Dear all I'm doing some Python scripts and I want to retrieve r.stats output information in a variable. How can I do this? variable= grass.run_command(r.series, .. ) Thanks Antonio ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Produce very-high quality images to include in brochures and flyers
Hi Markus Thanks, Does that work on Windows? Regards, A.Rocha On 03-12-2012 14:57, Markus Neteler wrote: On Fri, Nov 30, 2012 at 4:22 PM, António Rocha antonio.ro...@deimos.com.pt wrote: Hi I have a few satellite images and products and I want to generate very high quality screenshots/images based on them to place on brochures/etc. Any tips on how to do this because if I use r.out.gdal its images quality is not good. You can use d.out.png to turn the monitor content into a high res PNG. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Produce very-high quality images to include in brochures and flyers
Hi I have a few satellite images and products and I want to generate very high quality screenshots/images based on them to place on brochures/etc. Any tips on how to do this because if I use r.out.gdal its images quality is not good. Thanks Antonio ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Road and river extraction using GRASS
Greetings, Is there any GRASS module that can be used to river and road extraction from satellite images? Some feature extraction algorithm Thank you Regards, Antonio Rocha ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Error while saving display on Map Display- wxpython tool
Hi When I try to use Save display to graphic file in Windows I get this error: Traceback (most recent call last): File C:\Osgeo\apps\grass\grass-6.4.2\etc\wxpython\gui_mod ules\mapdisp.py, line 1159, in SaveToFile width, height) File C:\Osgeo\apps\grass\grass-6.4.2\etc\wxpython\gui_mod ules\mapdisp_window.py, line 575, in SaveToFile dc = wx.BufferedPaintDC(self, ibuffer) File C:\Osgeo\apps\Python27\lib\site-packages\wx-2.8-msw- unicode\wx\_gdi.py, line 4990, in __init__ _gdi_.BufferedPaintDC_swiginit(self,_gdi_.new_BufferedPaintD C(*args, **kwargs)) wx._core . PyAssertionError : C++ assertion wxAssertFailure failed at ..\..\src\msw\dcclient.cpp(219) in wxPaintDC::wxPaintDC(): wxPaintDC may be created only in EVT_PAINT handler! Traceback (most recent call last): File C:\Osgeo\apps\grass\grass-6.4.2\etc\wxpython\gui_mod ules\mapdisp.py, line 1159, in SaveToFile width, height) File C:\Osgeo\apps\grass\grass-6.4.2\etc\wxpython\gui_mod ules\mapdisp_window.py, line 575, in SaveToFile dc = wx.BufferedPaintDC(self, ibuffer) File C:\Osgeo\apps\Python27\lib\site-packages\wx-2.8-msw- unicode\wx\_gdi.py, line 4990, in __init__ _gdi_.BufferedPaintDC_swiginit(self,_gdi_.new_BufferedPaintD C(*args, **kwargs)) wx._core . PyAssertionError : C++ assertion wxAssertFailure failed at ..\..\src\msw\dcclient.cpp(219) in wxPaintDC::wxPaintDC(): wxPaintDC may be created only in EVT_PAINT handler! My GRASS version is 6.4.2 (installed with OSGEO4w). Any tips on how to fix this? Thanks Antonio -- António Rocha Earth Observation Applications and Services Av. D. João II, Lt 1.17.01 - 10º, Ed. Torre Zen 1998-023 Lisboa, Portugal Ph: +351 21 893 30 23 Fax: +351 21 896 90 99 http://www.deimos.pt ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Using r.out.kml and scripts in Windows Binary- error
Hello Martin I tried from the command line in the GUI. Ok next time I will give a try with that extension. Thanks On 13-09-2012 10:40, Martin Landa wrote: Hi, 2012/9/13 António Rocha antonio.ro...@deimos.com.pt: I tried to install r.out.kml with g.extension but I got this error: which: svn: unknown command ERROR: svn client required. Please install Subversion first. (maybe it is relevant to indicate that in the help file and a few links to Subversion clients) have you tried to install extension from command line or wxGUI [1]? wxGUI uses another g.extension (Python based) which should download on Windows directly prepared binaries. Martin [1] http://grass.osgeo.org/wiki/WxGUI#Extension_Manager -- António Rocha Earth Observation Applications and Services Av. D. João II, Lt 1.17.01 - 10º, Ed. Torre Zen 1998-023 Lisboa, Portugal Ph: +351 21 893 30 23 Fax: +351 21 896 90 99 http://www.deimos.pt ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Using r.out.kml and scripts in Windows Binary- error
Greetings I have used in Linux the r.out.kml script. Now, I hava installed Grass GIS (WinGRASS-6.4.2-2-Setup.exe) in Windows and I have placed my r.out.kml file in scripts folder and created a r.out.kml.bat file (with the path to the file) in bin folder. When I write r.out.kml in command console I get: Traceback (most recent call last): File C:\GRASS642\etc\wxpython\gui_modules\prompt.py, line 1105, in OnKeyPressed self.parent.RunCmd(cmd) File C:\GRASS642\etc\wxpython\gui_modules\goutput.py, line 523, in RunCmd task = menuform.GUI(show = None).ParseCommand(command) File C:\GRASS642\etc\wxpython\gui_modules\menuform.py, line 1732, in ParseCommand blackList = _blackList) File C:\GRASS642\etc\python\grass\script\task.py, line 460, in parse_interface tree = etree.fromstring(get_interface_description(name)) File C:\GRASS642\Python27\lib\xml\etree\ElementTree.py, line 1283, in XML return parser.close() File C:\GRASS642\Python27\lib\xml\etree\ElementTree.py, line 1636, in close self._raiseerror(v) File C:\GRASS642\Python27\lib\xml\etree\ElementTree.py, line 1488, in _raiseerror raise err xml.etree.ElementTree . ParseError : no element found: line 1, column 0 When I try to run other scripts (e.g. v.report) I get also this error: Traceback (most recent call last): File C:\GRASS642\etc\wxpython\gui_modules\prompt.py, line 1105, in OnKeyPressed self.parent.RunCmd(cmd) File C:\GRASS642\etc\wxpython\gui_modules\goutput.py, line 523, in RunCmd task = menuform.GUI(show = None).ParseCommand(command) File C:\GRASS642\etc\wxpython\gui_modules\menuform.py, line 1732, in ParseCommand blackList = _blackList) File C:\GRASS642\etc\python\grass\script\task.py, line 460, in parse_interface tree = etree.fromstring(get_interface_description(name)) File C:\GRASS642\Python27\lib\xml\etree\ElementTree.py, line 1283, in XML return parser.close() File C:\GRASS642\Python27\lib\xml\etree\ElementTree.py, line 1636, in close self._raiseerror(v) File C:\GRASS642\Python27\lib\xml\etree\ElementTree.py, line 1488, in _raiseerror raise err xml.etree.ElementTree . ParseError : no element found: line 1, column 0 What shall be done? Shall I downgrade to a previous vertsion? Thanks Antonio ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Possibility of defining a Map display windows size when launching?
Greetings Is it possible to predefine a map Display windows size when launching? If yes, where is it defined? Thanks Best regards Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6670 (2029) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Problems while running d.out.file
Hi Sorry to get back to this question again but I was not able to solve this. Ok 1- I have plotted my raster 2- I have added a command layer with this d.barscale bcolor=none at=0,90 3- If I use the Save display to graphic file functionality it works just fine. But If i try to set a smaller Width and Height: a) if I place barscale at 0,90 it gets cut b) If I place barscale at 0,0 it does not get cut but barscale gets over the image Is there any other suggestion on resizing PNG output without loosing image and barscale? Another question: when I launch a new Map Display window: Is it possible to define its size? Thanks Antonio Hamish wrote: Moritz wrote: I can confirm the issue (António, have you posted a bug report ?). The problem is not that the scalebar disappears, but that its location is not adapted to the changed resolution: see the http://geog-pc40.ulb.ac.be/grass/export_png/ for a series of examples. I'm not sure I fully understand the problem, but d.out.file has always worked well for me. the effect of making the width and height of the image is the same as dragging the xmonitor window bigger- the barscale and text is still the same size when viewed at 100% scale. when shrunk to fit on the screen in an image viewer it becomes tiny. You can make text and line widths bigger to compensate but really if you want things to scale well e.g. for a hi-res poster image, you should probably be using ps.map, not screen dumps of the Xmonitors, as the d.* modules were designed for. barscale and legend are placed as a percentage of the display frame, so at least their relative position should stay intact even if the frame gets bigger. if y'all think there's a bug please file a detailed ticket with the steps to reproduce using one of the standard datasets, and I'll look into it. Hamish __ Information from ESET NOD32 Antivirus, version of virus signature database 6639 (2017) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database (2028) __ 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/mailman/listinfo/grass-user
[GRASS-user] Generate a box that matches the exact processing region
Greetings I would like to know how can I generate a polygon for the exact computational region (a square). Is this possible? IF so, how can I do that. Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6654 (2023) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Problems while running Export to PNG file
Hi (I have eddited the email subject) but when you export the map composition the overlays don't get the memo that they should be propor- tionately place themselves not simply as a fixed offset from the bottom-left corner. (? I'm not exactly sure how the overlay positioning works, but assume that's what's happening) I confirm. this is what it is happening. The quick work-around is for the user to use d.legend and d.scalebar as command layers in the main map layer wx gui (somewhat hidden under the add grid or vector labels icon) with manual command line controls. Then the legend or barscale gets rendered directly as part of map rendering, not overplotted as a movable overlay. Ok I was able to add a scale bar as command layer. When I export it the recoomended size it is plotted just fine and with a different size it addapted thje scale bar (it remained a bit bigger). Is it possilve to define Scalebar size (not just position)? Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6648 (2021) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Problems while running d.out.file
Hi Hamish THe problem is with the shrink. You are right it's tiny but in the case of the scale bar it is not shrinked. I mean it keeps the original size. Example: This is the original http://s1183.photobucket.com/albums/x462/toyzerocha/?action=viewcurrent=original.png This is half of the size http://s1183.photobucket.com/albums/x462/toyzerocha/?action=viewcurrent=Half_size.png In this is case it does not disappear it's just cut. I would report this as bug. Or not really? Antonio Hamish wrote: Moritz wrote: I can confirm the issue (António, have you posted a bug report ?). The problem is not that the scalebar disappears, but that its location is not adapted to the changed resolution: see the http://geog-pc40.ulb.ac.be/grass/export_png/ for a series of examples. I'm not sure I fully understand the problem, but d.out.file has always worked well for me. the effect of making the width and height of the image is the same as dragging the xmonitor window bigger- the barscale and text is still the same size when viewed at 100% scale. when shrunk to fit on the screen in an image viewer it becomes tiny. You can make text and line widths bigger to compensate but really if you want things to scale well e.g. for a hi-res poster image, you should probably be using ps.map, not screen dumps of the Xmonitors, as the d.* modules were designed for. barscale and legend are placed as a percentage of the display frame, so at least their relative position should stay intact even if the frame gets bigger. if y'all think there's a bug please file a detailed ticket with the steps to reproduce using one of the standard datasets, and I'll look into it. Hamish __ Information from ESET NOD32 Antivirus, version of virus signature database 6639 (2017) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 6639 (2017) __ 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/mailman/listinfo/grass-user
[GRASS-user] Problems while running d.out.file
Greetings I'm trying to run d.out.file to produce a PNG file with a defined size (and others) but I get the following error: ERROR: You must select a display monitor. SInce I'm running from WxPython interface how can I do that? THanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6633 (2015) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Problems while running d.out.file
Hi Martin About that functionality. d.out.file allows me to define size of output image which is something that I need. WHen I try to do that with the functionality that you mention (e.g. with smaller H and W) I get the error I reported here (http://lists.osgeo.org/pipermail/grass-user/2011-November/062577.html). Michael Barton suggested me to use this function. I really need to produce output maps with a specific size and with image+scale bar. and when I set HxW I loose Scalebar in my output image. Antonio Martin Landa wrote: Hi, 2011/11/16 António Rocha antonio.ro...@deimos.com.pt: I'm trying to run d.out.file to produce a PNG file with a defined size (and others) but I get the following error: ERROR: You must select a display monitor. SInce I'm running from WxPython interface how can I do that? you need to start some monitor, eg. `d.mon start=PNG` and then stop it. Why are you not using functionality available from wxGUI Map Display [1]? Martin [1] http://www.youtube.com/watch?v=PY8aytvfQqAfeature=player_embedded __ Information from ESET NOD32 Antivirus, version of virus signature database 6634 (2016) __ 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/mailman/listinfo/grass-user
[GRASS-user] Running a Python script in GRASS6.4.1-RC2
Greetings I have built a python script and it runs just fine in Windows GRASS. when I try to compile it in Linux it places the script file in script folder but when I try to run It I get the following error in Command Outputwindow: Traceback (most recent call last): File /usr/local/grass-6.4.2RC1/etc/wxpython/gui_modules/p rompt.py, line 1098, in OnKeyPressed self.parent.RunCmd(cmd) File /usr/local/grass-6.4.2RC1/etc/wxpython/gui_modules/g output.py, line 534, in RunCmd task = gtask.parse_interface(command[0]) File /root/Downloads/grass/grass-6.4.2RC1/dist.x86_64 -unknown-linux-gnu/etc/python/grass/script/task.py, line 451, in parse_interface File /usr/lib/python2.4/site- packages/elementtree/ElementTree.py, line 961, in XML return parser.close() File /usr/lib/python2.4/site- packages/elementtree/ElementTree.py, line 1251, in close self._parser.Parse(, 1) # end of data xml.parsers.expat . ExpatError : no element found: line 1, column 0 and it does not launch the modules Interface Any idea of what is happening and how to fix this? Thanks Regards Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6601 (2004) __ 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/mailman/listinfo/grass-user
[GRASS-user] Define statistical bins in GRASS
Greetings I have a raster image that I want to analyze its frequency inside defined bins. How can I define bins size and bins-limits in GRASS? THanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6578 (20111027) __ 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/mailman/listinfo/grass-user
[GRASS-user] Running Python Script in GRASS6.4 in Linux
Greetings I have been able to run a GRASS Python script in 6.4. version in Windows by creating a .bat file in bin and adding this at the top: #!/usr/bin/env python But now, in Linux I'm not able to run this. I habve placed my script file in script folder but I'm getting this error: Traceback (most recent call last): File /usr/lib64/grass-6.4.1/etc/wxpython/gui_modules/prompt.py, line 1070, in OnKeyPressed self.parent.RunCmd(cmd) File /usr/lib64/grass-6.4.1/etc/wxpython/gui_modules/gout put.py, line 516, in RunCmd task = menuform.GUI().ParseInterface(command) File /usr/lib64/grass-6.4.1/etc/wxpython/gui_modules/menu form.py, line 2020, in ParseInterface tree = etree.fromstring(getInterfaceDescription(cmd[0])) File /usr/lib64/grass-6.4.1/etc/wxpython/gui_modules/menu form.py, line 1969, in getInterfaceDescription raise gcmd.GException, _(Unable to fetch interface description for command '%s'. AttributeError : OSError instance has no attribute 'value' What is missing? THanks Antonio -- António Rocha Aerospace Defence Portugal DEIMOS Engenharia, SA Av. D. Joao II, Lt 1.17.01 - 10º, Ed. Torre Zen 1998-023 Lisboa, Portugal Ph: +351 21 893 3010 Fax: +351 21 896 9099 http://www.deimos.pt __ Information from ESET NOD32 Antivirus, version of virus signature database 6562 (20111021) __ 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/mailman/listinfo/grass-user
[GRASS-user] Using r.reclass.area with an NDVI map
Greetings I have an NDVI map but it's not continuous: i mean, I have a few patches with values and areas in my raster map without values. I want to apply r.reclass.area in order to eliminate patches with less than 2ha. My question is: shall I convert NDVI pixels to 1 (creating a mask of valid pixels) or can I just use as it is (with values ranging from -1 to 1)? Because if r.reclass.area uses r.clump 2 pixels with different NDVI values are not categorized toghetter. Right? THanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6561 (20111020) __ 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/mailman/listinfo/grass-user
[GRASS-user] Error in r.regression.line manual- R output
Greetings In r.regression.line R is being defined as residuals. Well, after doing some calculation in other tools I realized that R corresponds to Pearson Correlation (also known as R). Can anyone else confirm me this? THanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6512 (20111003) __ 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/mailman/listinfo/grass-user
[GRASS-user] Linear regression between 2 rasters (old and new)
I I have been using r.series to calculate a linear regression slope for an old and new raster map using r.series. But now I need to also obtain the error. It seems not to be possible to use r.series. Is this true? if yes, what alternative do I have? THanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6506 (20110930) __ 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/mailman/listinfo/grass-user
[GRASS-user] Where GRASS_DOC_BASE variable is stored?
Greetings I would like to know where GRASS_DOC_BASE is stored and/or how is builted? THanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6482 (20110921) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Using non-ASCII characters in Legend display
Hello Moritz Thanks. ISO-8859-1 works. Just one question: I'm trying to find where this parameter was stored and I'm not able tyo find. How can I know where these preferences are stored? (I'm using Windows) Thanks Antonio Moritz Lennert wrote: On 09/09/11 09:58, António Rocha wrote: Hi Martin I know where to set the font and the encoding. What I was asking is a combination of Font and Encoding that allows me to Display those characters in the legend and in Map Display window? iso-8859-1 or UTF8 should do the trick in terms of encodings. Most of the fonts, except for the original basic grass fonts can handle accents, e.g. Vera, Bitstream, DejaVu, ... Moritz __ Information from ESET NOD32 Antivirus, version of virus signature database 6455 (20110911) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 6455 (20110911) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Using non-ASCII characters in Legend display
Hi Martin I know where to set the font and the encoding. What I was asking is a combination of Font and Encoding that allows me to Display those characters in the legend and in Map Display window? Thanks Antonio Martin Landa wrote: 2011/9/5 António Rocha antonio.ro...@deimos.com.pt: Does anyone knows a GRASS_FONT and GRASS_ENCODING configuration tyhat allows to display non-ascii characters in Map Display Windows legend (e.g. ááéã) GUI Settings - Map Display - Font settings - Default font for GRASS displays Martin __ Information from ESET NOD32 Antivirus, version of virus signature database 6448 (20110908) __ 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/mailman/listinfo/grass-user
[GRASS-user] Re: Using non-ASCII characters in Legend display
Hi (back on this again!) Does anyone knows a GRASS_FONT and GRASS_ENCODING configuration tyhat allows to display non-ascii characters in Map Display Windows legend (e.g. ááéã) Thanks Antonio António Rocha wrote: Greetings I'm not being able to display non-ascii characters in Legend (e.g. á à ã â ç). Instead of the chatracter a blank space is displayed. (but with query tool, I can display labels full name in Command Output window) I know that I shall select encoding and Font type that supports these characters. I have tried several of them and no LUCK. Has anyone tried this before? Or anyone has any suggestion? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6436 (20110904) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Where are the configuration and preferences saved?
Hi Well in Windows I'm not being able to figure out where that information is stored Michael Barton wrote: For Mac and Linux at least, they are saved in ~/.grass6 and ~/.grass7 directories. These are hidden in the Mac finder (unless you use something like Invisibility Toggler.app) but accessible through the terminal. I'm not sure where these directories live in Windows. Michael C. Michael Barton Director, Center for Social Dynamics Complexity Professor of Anthropology, School of Human Evolution Social Change Arizona State University voice: 480-965-6262 (SHESC), 480-727-9746 (CSDC) fax: 480-965-7671 (SHESC), 480-727-0709 (CSDC) www: http://www.public.asu.edu/~cmbarton, http://csdc.asu.edu On Aug 12, 2011, at 9:00 AM, grass-user-requ...@lists.osgeo.org grass-user-requ...@lists.osgeo.org wrote: Date: Thu, 11 Aug 2011 15:55:18 +0100 From: Ant?nio Rocha antonio.ro...@deimos.com.pt Subject: [GRASS-user] Where are the configuration and preferences saved? To: GRASS user list grass-user@lists.osgeo.org Message-ID: 4e43ed56.4000...@deimos.com.pt Content-Type: text/plain; charset=ISO-8859-1; format=flowed Greetings I Would like to know in which GRASS file are saved all configuration and definitions? Thanks __ Information from ESET NOD32 Antivirus, version of virus signature database 6370 (20110811) __ 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/mailman/listinfo/grass-user
[GRASS-user] Where are the configuration and preferences saved?
Greetings I Would like to know in which GRASS file are saved all configuration and definitions? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6369 (20110811) __ 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/mailman/listinfo/grass-user
[GRASS-user] Using i.rectify and data type
Greetings I have one question regarding i.rectify's data type output: If I have my input is a CELL map, will my output will be also CELL? I mean, does it keep my data type? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6313 (20110721) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Using i.rectify and data type
Hello Markus Are you sure? Because I have tried i.rectify (1st order, 2nd order and 3rd Order) with a CELL image and my output was also CELL. Markus Metz wrote: 2011/7/21 António Rocha antonio.ro...@deimos.com.pt: Greetings I have one question regarding i.rectify's data type output: If I have my input is a CELL map, will my output will be also CELL? I mean, does it keep my data type? That depends on the resampling method. Only nearest neighbor preserves the datatype, with all other resampling methods the output will be DCELL. This is similar to r.proj, only that r.proj uses FCELL for resampling methods other than nearest neighbor. __ Information from ESET NOD32 Antivirus, version of virus signature database 6314 (20110721) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 6314 (20110721) __ 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/mailman/listinfo/grass-user
[GRASS-user] Re: Wxpython module's interface- multiple entries
Greetings Clearly it needs to be fixed. never the less, I have been able to delete 1 or two entries. By selecting all entries and deleting 1 by selecting the faction of text. but it's not very user friendly Michael Barton wrote: Yes. This has been around for some time. Looks like it may need to have a bug report filed as I cannot find one. It definitely needs to be fixed. The only workaround is to close the module and reopen it. Michael __ Information from ESET NOD32 Antivirus, version of virus signature database 6271 (20110706) __ 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/mailman/listinfo/grass-user
[GRASS-user] GRASS version on yum of CEntOS
Greetings I have just installed GRASS in a CentOS machine, using yum repository but the version is GRASS6.4.0-2. How old is this version? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6246 (20110628) __ 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/mailman/listinfo/grass-user
[GRASS-user] Grass-Windows - Error using tcltk's messages
Hi I have exacly the same problem. Is this some limitation on building these installers? Or it may be some error that both of us did Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6188 (20110607) __ 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/mailman/listinfo/grass-user
[GRASS-user] Pan-sharpening with GRASS
Greetings I would like to know if there is any method or function to perform pan-sharpening (image fusion) in GRASS? THanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6171 (20110601) __ 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/mailman/listinfo/grass-user
[GRASS-user] Limits on size of command line
Greetings I would like to know if there is any limit on command line size? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6164 (20110530) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] using not so common characters in Python Script
Olá Ricardo :) Thanks for the feedback. Just one last question: is there any known limitation for using this encoding in Grass Python Scripts? ( # -*- coding: utf-8 -*- Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6133 (20110518) __ 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/mailman/listinfo/grass-user
[GRASS-user] using not so common characters in Python Script
Greetings I would like to include some sentences with characters like á ã àó in my grass.message strings. It seems it's not possible Because I get this grass.fatal(_('á')) I get: File C:/local/scripts/x.teste, line 41 SyntaxError: Non-ASCII character '\xe1' in file C:/local/scripts/x.teste on line 41, but no encoding declared; see http://www.python.org/peps/pep-0263.html for details Is it possible? Thanks Antonio Rocha __ Information from ESET NOD32 Antivirus, version of virus signature database 6132 (20110518) __ 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/mailman/listinfo/grass-user
[GRASS-user] Re:Translation of buttons
Greetings I just saw this message (bellow) in GRASS-translation mailing list and I decided to duplicate this topic here: Besides the buttons that I mentioned yesterday, I nottest that also when I open a script (e.g. i.tasscap) Only the run button is translated. Close, Stop, Copy and help are not translated. where is this managed? because if Run is translated we will be able to translate the others right? Thanks Monica Some buttons in WxPYthon interface cannot be translated (Help and Quit in front panel; next, cancel and back in Location Wizard and Stop, Cancel Copy in Scripts interface. Is this true? Is it possible to, at least, manually change this buttons labels? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6129 (20110517) __ 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/mailman/listinfo/grass-user
[GRASS-user] Error using grass.Popen
Greetings I'm using grass.Popen like this: externalProcess = grass.Popen(cmd, stdout=grass.PIPE, stderr=grass.PIPE, cwd=tempfile) It seems that I'm getting some error while running this function but the only thing that is printed in my command Output window is: grass.script.core.Popen object at 0x00D6ECD0 How can I get the exact error that I'mgetting? Thanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6003 (20110331) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Error while removing file
Hello Glynn I'm running in WinGRASS in WinXP and I have total access to my folders (as far as I can see) but I get Permission denied for a specific folder in C:\grassdata\Area1\Country\cell_misc and ACCESS DENIED. I have built a very long and complex Script with a lot of files production, deleting, opening external configuration files bla bla and now I'm getting this sort of errors when I try to overwrite over existing data or when I try to delete some files. What might be happening? Thanks Antonio Glynn Clements wrote: António Rocha wrote: Ok it has removed but when I try to write over that file I get errors like Permission DENIED. WHat can I do to fix this? Ensure that all files and directories within the mapset directory are owned by you, and have the necessary permissions (read and write for files, read, write and execute for directories). __ Information from ESET NOD32 Antivirus, version of virus signature database 5956 (20110315) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Error while removing file
Hi again :) Ok the main problem is that I'm using r.statistics with base=new_clump10 and cover= myscript.tmpV1 If I try to remove new_clump I get: g.remove -f rast=new_clump10@area [new_clump10@area@area] is a base map for [myscript.tmpV1@area]. Remove forced. Removing raster LULC_new_clump10@National misc: couldn't be removed So it wasn't removed. My question is: 1- Why it's repeating the Mapset information (new_clump10@area@area) That is odd right? Maybe that is what is messing with the system Thanks Antonio Glynn Clements wrote: António Rocha wrote: Ok it has removed but when I try to write over that file I get errors like Permission DENIED. WHat can I do to fix this? Ensure that all files and directories within the mapset directory are owned by you, and have the necessary permissions (read and write for files, read, write and execute for directories). __ Information from ESET NOD32 Antivirus, version of virus signature database 5957 (20110316) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Error while removing file
Just one small note: If I try to remove tmpV I get this Removing information about reclassed map from [new_clump10@area] failed Removing raster myscript.tmpV1@National And it's removed but then If I try to remove the new_clump10 I get: [new_clump11@area@area] is a base map for [myscript.tmpV1@area]. Remove forced. Removing raster new_clump11@area misc: couldn't be removed António Rocha wrote: Hi again :) Ok the main problem is that I'm using r.statistics with base=new_clump10 and cover= myscript.tmpV1 If I try to remove new_clump I get: g.remove -f rast=new_clump10@area [new_clump10@area@area] is a base map for [myscript.tmpV1@area]. Remove forced. Removing raster LULC_new_clump10@National misc: couldn't be removed So it wasn't removed. My question is: 1- Why it's repeating the Mapset information (new_clump10@area@area) That is odd right? Maybe that is what is messing with the system Thanks Antonio Glynn Clements wrote: António Rocha wrote: Ok it has removed but when I try to write over that file I get errors like Permission DENIED. WHat can I do to fix this? Ensure that all files and directories within the mapset directory are owned by you, and have the necessary permissions (read and write for files, read, write and execute for directories). __ Information from ESET NOD32 Antivirus, version of virus signature database 5957 (20110316) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Error while removing file
Ok it has removed but when I try to write over that file I get errors like Permission DENIED. WHat can I do to fix this? Glynn Clements wrote: António Rocha wrote: I'm trying to remove a file named new_clump@Area and I get this error: [new_clump@Area] is a base map for [myscript.tmpV@Area]. Remove forced. is a base map for is just a warning; it doesn't prevent removal so long as the -f flag is given. misc: couldn't be removed This indicates that there was some problem with removing the cell_misc directory. __ Information from ESET NOD32 Antivirus, version of virus signature database 5955 (20110315) __ 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/mailman/listinfo/grass-user
[GRASS-user] Error while removing file
Greetings I'm trying to remove a file named new_clump@Area and I get this error: [new_clump@Area] is a base map for [myscript.tmpV@Area]. Remove forced. misc: couldn't be removed I'm already using -f flag from g.remove function. myscript.tmpV was a file that I generated from this new_clump. But, now I don't have that file in Area MAPSET ( I already checked in the Area folder) so I wasn't suppose to get this error. How can I remove this raster without getting this nonsense error Thanks Best regards, Antonio R. __ Information from ESET NOD32 Antivirus, version of virus signature database 5952 (20110314) __ 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/mailman/listinfo/grass-user
[GRASS-user] Running an external executable file from a Python Script
Greetings all, I have a python script from where I need to run an external binary (in Windows). This external binary (not developed by me) uses an input_parameter file to configure binary and get other parameters. The problem is that this binary requires that in my active folder I have my parameter file. Example: If I want to run this binary (outside GRASS Python Scripts) i open cmd, go to a folder where I have a parameter file (e.g. c:\delete (cd c:\delete)) Then I can run this binary as long as I have my parameter file in my active folder like this (c:\tool\training.exe). What I mean is that in my active folder I do not need to have my binary only my parameter file. My question is, when I'm running a GRASS python Script what is my active folder in order to place there my Parameter file? Or, is there any way to change my active folder while I'm running GRASS python Script? Hope I was clear enough in this email since this was unexpected for me. Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5941 (20110310) __ 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/mailman/listinfo/grass-user
[GRASS-user] Using output from r.statistics
Greetings I have used r.statistics to calculate the sum of a cover map based on a base map. When I use the Query tool I get these: 598666.544966|9159452.159237||3996|7108 The first two corresponds to Geographic Location the third Base map clump ID and the last one it corresponds to the value of the sum (r.statistics). When I try to use this output, in mapcalc, it does not use the values obtained in r.statistics but the clump value Which does not make much sense. How can I use the r.average valeus of this raster? Thanks Antonio Rocha __ Information from ESET NOD32 Antivirus, version of virus signature database 5918 (20110301) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] using r.resamp.stats- When I reach a tie
Hi Glynn. Thanks for your answer. Is there any other alternative, to resample data, that uses another approach to solve tie events? THanks Antonio Glynn Clements wrote: António Rocha wrote: I'm planning on using r.resamp.stats in order to resample a grid to a coarser resolution. Since I'm using this function with classes I will use mode method (more frequentclass to fall into a destination cell is the chosen). My question is: when I have a tie between two classes, how does GRASS decides which class shall be chosen? In the event of a tie, the lowest-numbered category is chosen. __ Information from ESET NOD32 Antivirus, version of virus signature database 5912 (20110227) __ 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/mailman/listinfo/grass-user
[GRASS-user] using r.resamp.stats- When I reach a tie
Greetings I'm planning on using r.resamp.stats in order to resample a grid to a coarser resolution. Since I'm using this function with classes I will use mode method (more frequentclass to fall into a destination cell is the chosen). My question is: when I have a tie between two classes, how does GRASS decides which class shall be chosen? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5906 (20110225) __ 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/mailman/listinfo/grass-user
[GRASS-user] i.gensig- Error because I'm obtaining negative values in covariance diagonal
Greetings I'm using i.gensig and I'm having some problmes while generating signature files. Why? bcause my covariance calculus is placing negative values in the diagonal which is impossible. Has anyone detected any problem with i.gensig? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5877 (20110215) __ 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/mailman/listinfo/grass-user
[GRASS-user] List files and count
Greetings I need to find, in a certain mapset, a list of files with a name that match a certain string and count them. How can I do this? I'm implementing this in a Py Script Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5774 (20110110) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Generalization or interpolation using the majority
Hi Markus It seems that it does exacly what I want. I just have one question: I need to first set the region spatial resolution to 200 before I run r.resamp.stats right? THanks Antonio Markus Metz wrote: 2010/12/21 António Rocha antonio.ro...@deimos.com.pt: Greetings I have a raster classification map with a 100 meters resolution and I want to interpolate to 200m resolution using GRASS and, the output value must be the most frequent value (majority). r.resamp.stats method=mode would do that Markus M Here goes an examples. Intput: 1 1 1 1 1 2 3 4 1 2 2 3 2 2 3 3 2 2 5 2 2 2 2 2 3 2 5 5 1 3 4 5 Output must be: 1 1 2 3 2 5 2 2 My question is what method do you advise to be used? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5720 (20101221) __ 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/mailman/listinfo/grass-user __ Information from ESET NOD32 Antivirus, version of virus signature database 5720 (20101221) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5720 (20101221) __ 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/mailman/listinfo/grass-user
[GRASS-user] WinGRASS: Error when importing GRIB files
Hi I'm using WinGRASS (6.4.1) in Osgeo4w and I'm having a problem in the import of GRIB files. I'm using r.in.gdal, as I have used in Linux, and I'm obtaining the following error: ERROR 4: `C:\Data\GISDataBase\GRIB_codes\new03.grib' not recognised as a supported file format. What might be missing? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5657 (20101129) __ 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/mailman/listinfo/grass-user
[GRASS-user] i.gensig comment- Signature 1 not positive definite
Greetings For a few training areas I'm obtaining the following comment (which means that my signature files are created but empty): Signature 1 not positive definite What does this means? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5647 (20101125) __ 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/mailman/listinfo/grass-user
[GRASS-user] Re: [GRASS-dev] ERROR: G_realloc: during r.to.vect call
Hello Markus I believe it's ok: projection: 1 (UTM) zone: -24 datum: wgs84 ellipsoid: wgs84 north: 9305500.0001 south: 9095770.0001 west: 444000 east: 689130 nsres: 30 ewres: 30 rows: 6991 cols: 8171 cells: 57123461 And my raster is: | Type of Map: raster Number of Categories: 3 | | Data Type: CELL | | Rows: 6991 | | Columns: 8171 | | Total Cells: 57123461 | |Projection: UTM (zone -24) | |N: 9305500.0001S: 9095770.0001 Res: 30 | |E: 689130W: 444000 Res: 30 | | Range of data:min = 0 max = 3 | So it matches. What might be causing this error? Antonio Markus Neteler wrote: 2010/11/18 António Rocha antonio.ro...@deimos.com.pt: Greetings I'm running GRASS (in Windows) and, for a specific raster, when I try to run v.to.rast I get this error: ERROR: G_realloc: unable to allocate 3376 bytes at areas.c:678 1- Is this expectable? 2- is there any fix for this? If not, does anyone has any suggestion to avoid this error? Did you check the computation region with g.region -p ? Just to be sure that extent and resolution are right... Markus __ Information from ESET NOD32 Antivirus, version of virus signature database 5631 (20101118) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5632 (20101119) __ 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/mailman/listinfo/grass-user
[GRASS-user] Using i.gensig- How to keep original values
Greetings I have a training raster map with values such as 3; 5; 7 ; 12 and, when I run i.gensig those values are not included in the signature file. Which means that, if I run i.maxlik, I need to do some reclassification. Is it possible to run i.gensig and keep class values in the signature? THanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5625 (20101116) __ 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/mailman/listinfo/grass-user
[GRASS-user] Re: v.colors it seems to be working but...
Luisa I believe it is pretty simple now . When in d.vect check -a flag (Get colors from map table column) I believe it works just fine (that final part.. about applying v.colors I had the same problems) Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5626 (20101117) __ 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/mailman/listinfo/grass-user
[GRASS-user] Check if a raster exists in a new region
Greetings I've created a raster based on region (let's call it A). Now, I have changed my region (B for example) and I need to check (in a Py Script) if my raster has any valid value in this new region. How can I do that? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5622 (20101115) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Error launching IE from GRASS for manuals
Helmut is absolutely right. Not only the WinBinaries but also my own WinGRASS installers are all ok when it comes to HTML Browser. only my OSGEO4W is not being able to launch IE. I've fixed ticket 980 in all cases in the nsis-installer, so the available standalone WinGrass-installer should have incorporated all fixes. AFAIK this is an osgeo4w issue. in the grass-starting routine there are following variables set in C:\OSGeo4W\etc\ini\grass set GISBASE=%OSGEO4W_ROOT%\apps\grass\grass-6.4.0svn set WINGISBASE=%OSGEO4W_ROOT%\apps\grass\grass-6.4.0svn set GRASS_SH=%OSGEO4W_ROOT%\apps\msys\bin\sh.exe set GRASS_WISH=%OSGEO4W_ROOT%\bin\wish.exe set GRASS_PYTHON=%OSGEO4W_ROOT%\bin\python.exe set GRASS_PROJSHARE=%OSGEO4W_ROOT%\share\proj set GRASS_HTML_BROWSER=%PROGRAMFILES%\Internet Explorer\iexplore.exe so this should be fixed there. please open a ticket there. best regards Helmut __ Information from ESET NOD32 Antivirus, version of virus signature database 5602 (20101108) __ 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/mailman/listinfo/grass-user
[GRASS-user] Error launching IE from GRASS for manuals
Greetings I have a wiNGRASS installation, based on OSGEO4W GRASS compilation, but I'm having a difficulty on launching manuals. When i run g.manual wxGUI I get this error: Starting browser C:\Program Files\Internet Explorer\iexplore.exe for module wxGUI... C:/OSGeo4W/apps/grass/grass-6.4.1svn/scripts/g.manual: C:\Program Files\Internet Explorer\iexplore.exe: command not found But My IE is installed in C:\Program Files\Internet Explorer\iexplore.exe so It's suppose to work. What can be done to fix this? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5600 (20101108) __ 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/mailman/listinfo/grass-user
[GRASS-user] Error with r.mapcalc
I have a Python script with the following: # Conversion to at-sensor brightness temperature grass.mapcalc($raster_file= $K2/log(($K1/$L)+1), raster_file=raster_file, K2=K2, K1=K1,L= tempL) grass.mapcalc($raster_file= if($raster_file 0, null(), $raster_file), raster_file=raster_file) When it runs the second grass.mapcalc I get the following error: ERROR: Unable to close raster map ERROR: An error occurred while running r.mapcalc Why am I obtaining this error? and How can I avoid it? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5563 (20101026) __ 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/mailman/listinfo/grass-user
[GRASS-user] Re: self compile grass- msys
Hi there I have sent a similar message to this http://lists.osgeo.org/pipermail/grass-user/2010-September/057903.html but to the wingrass mailing list. Since there is no follow up on this subject in this mailing list I wanted to know if anyone found any answer to this? Because I'm having the same problem in winWXP 32. In the past I had this Windows GRASS compilation and now, I decided to have it again but I'm having this problem: checking whether to use Curses... yes checking for curses.h... yes checking curses.h WINDOW structure component... _maxy checking for initscr in -lncurses... no checking for initscr in -lcurses... configure: error: *** Unable to locate curses library. no Regarding -lncurses, my config.log has the following: int main() { WINDOW w; w.maxy = 0; ; return 0; } configure:5691: gcc -c -g -O2-I/c/OSGeo4W/apps/gdal-16/include -I/c/OSGeo4W/include conftest.c 15 configure:5707: checking for initscr in -lncurses configure:5724: gcc -o conftest.exe -g -O2 -I/c/OSGeo4W/apps/gdal-16/include -I/c/OSGeo4W/include -Wl,--export-dynamic,--enable-runtime-pseudo-reloc -L/c/OSGeo4W/apps/gdal-16/lib -L/c/OSGeo4W/lib conftest.c -lncurses 15 c:/mingw/bin/../lib/gcc/mingw32/4.5.0/../../../../mingw32/bin/ld.exe: warning: --export-dynamic is not supported for PE targets, did you mean --export-all-symbols? c:/mingw/bin/../lib/gcc/mingw32/4.5.0/../../../../mingw32/bin/ld.exe: cannot find -lncurses collect2: ld returned 1 exit status configure: failed program was: #line 5713 configure #include confdefs.h /* Override any gcc2 internal prototype to avoid an error. */ /* We use char because int might match the return type of a gcc2 builtin and then its argument prototype would still apply. */ char initscr(); And regarding lcurses, it's the following: int main() { initscr() ; return 0; } configure:5743: checking for initscr in -lcurses configure:5760: gcc -o conftest.exe -g -O2 -I/c/OSGeo4W/apps/gdal-16/include -I/c/OSGeo4W/include -Wl,--export-dynamic,--enable-runtime-pseudo-reloc -L/c/OSGeo4W/apps/gdal-16/lib -L/c/OSGeo4W/lib conftest.c -lcurses 15 c:/mingw/bin/../lib/gcc/mingw32/4.5.0/../../../../mingw32/bin/ld.exe: warning: --export-dynamic is not supported for PE targets, did you mean --export-all-symbols? c:/mingw/bin/../lib/gcc/mingw32/4.5.0/../../../crt2.o:crt1.c:(.text+0x28): undefined reference to `__dyn_tls_init_callback' collect2: ld returned 1 exit status configure: failed program was: #line 5749 configure #include confdefs.h /* Override any gcc2 internal prototype to avoid an error. */ /* We use char because int might match the return type of a gcc2 builtin and then its argument prototype would still apply. */ char initscr(); What is missing in my installation? Because I have checked, and my osgeo4W includes pdcurses-devel Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5509 (20101006) __ 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/mailman/listinfo/grass-user
[GRASS-user] [DWE]- Landsat Cloud detection with GRASS
Greetings I'm trying to find if GRASS includes some-sort of way to detect/eliminatd cloud cover and shadows from LANDSAT Images. The only thing I found was this: http://svn.osgeo.org/grass/grass-addons/imagery/i.landsat.acca/description.html But: 1- has anyone tried/used this? Is this useful and working? 2- can I use it to detect-and-eliminate clouds? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5486 (20100928) __ 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/mailman/listinfo/grass-user
[GRASS-user] temporary files in a Python Script with mapcalc
Greetings I'm doing some image processing python script with requires the calculation of many intermediate rasters. for example: x=a+b+c alpha=x-(b^gama) y=c*alpha+b One possibility is integrating all of it in one expression but I'm afraid of having some minor error in expression so I prefer to separate the calculation. My question is: how can I use temporary files to store x and alpha value? Or I should create normal rasters and, at the end, delete them? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5472 (20100923) __ 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/mailman/listinfo/grass-user
[GRASS-user] Compiling WINGRASS with OSGEO4w
Greetings I have successfully compiled svn6.4.0 versions with OSGEO4W compilation tutorial/wiki. Now, I have also successfully compiled svn6.4.1 version with OSGEO4W following the same steps but, when I select the GRASS -wxpython icon on my Desktop, instead of initiating GRASS6.4.1 is calls Init.bat from GRASSsvn6.4.0. How can I change/manage this? Beucase It's correctly installed without any errors. Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5454 (20100916) __ 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/mailman/listinfo/grass-user
[GRASS-user] Merging two vectors
Hi I have 2 vectors (each one is a square) that are not over-imposed and I need to merge them in a single file/vector. How can i do that? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5454 (20100916) __ 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/mailman/listinfo/grass-user
[GRASS-user] Georrectification- Display an RGB instead of just one band
Hi, I'm using Georrectification module (WxGUI) and I have one question: Before the launch of Georrectification panel, Select image/map to display for GCP creation is it possible to display an RGB composite instead of a single Grayscale image? Thanks Antonio R __ Information from ESET NOD32 Antivirus, version of virus signature database 5372 (20100817) __ 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/mailman/listinfo/grass-user
[GRASS-user] Python script- drop-down menu list of multiple groups
Greetings I'm implementing a Python Script and I need to know one thing related with Modules Interfaces: can I have a drop-down menu list of multiple groups? Instead of just selecting one group (e.g. i.maxlik) select more than one (multiple)? Thanks Best regards Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5371 (20100816) __ 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/mailman/listinfo/grass-user
[GRASS-user] Using r.univar.zonal- calculate most frequent value
Hi there, I have tested r.univar.zonal (pretty useful!!) but now I have one question: - Is it possible to calculate the most frequent value inside a zone? (I don't think so If not, is there any other way to calculate this in GRASS? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5345 (20100805) __ 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/mailman/listinfo/grass-user
[GRASS-user] Retrieve pixel position
Greetings I have retriecve a Pixel coordinates but instead of lat lon I want row/col. As far as I can see, r.what only retrieves Lat/lon not row/col. How can I do this? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5337 (20100803) __ 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/mailman/listinfo/grass-user
[GRASS-user] Re: v.to.rast for a vectorial without Values
Hello Luisa Regarding Hamish suggestion I would say to do this: 1- Define boundaries (without centroids) in Tcl/Tk digitizer 2- v.centroids 3- v.t.rast with attributes as cat. Regarding point 1, I don't know if you should add centroids. Maybe, v.centroids won't work. __ Information from ESET NOD32 Antivirus, version of virus signature database 5319 (20100728) __ 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/mailman/listinfo/grass-user
[GRASS-user] osgeo4w winGRASSinstallation: error with PROJ4
Greetings I've installed GRASS weekly snapshots in Windows (osgeo4w) by using CompileOnWindows Wiki. Today I tried to install latest snapshot (24_07_2010) and I got this error: checking External PROJ.4 version... 470 checking for location of External PROJ.4 library... checking for pj_get_def in -lproj... ./configure: line 7176: 5488 Segmentation fault (core dumped) ( eval $ac_link ) 25 configure: error: *** Unable to locate External PROJ.4 library. no Which is strange because last week I have installed a version from June 2010 and it all worked fine. Any idea of what might be wrong? Thanks Antonio Rocha __ Information from ESET NOD32 Antivirus, version of virus signature database 5314 (20100726) __ 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/mailman/listinfo/grass-user
[GRASS-user] Supervised Remote Sensing
Hello Etienne DELAY YOu have skipped the part where the signatures are produced. Take a look at: http://grass.osgeo.org/gdp/html_grass64/imageryintro.html in Image classification Hope it helps Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5297 (20100721) __ 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/mailman/listinfo/grass-user
[GRASS-user] Error while installing winGRASS in Osgeo4w
Greetings I have successfully installed winGRASS following this instructions: http://trac.osgeo.org/grass/wiki/CompileOnWindows I had a few problems with my current grass installation so I decided to Download grass-6.4.svn_src_snapshot_2010_07_18 version (from here: http://grass.osgeo.org/grass64/source/snapshot/) I have done this: cd /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_18 then: ./mswindows/osgeo4w/package.sh But I get this error: Wed Jul 21 15:57:40 GMTDT 2010: STARTING make distclean makefile:22: include/Make/Platform.make: No such file or directory makefile:23: include/Make/Grass.make: No such file or directory make: *** No rule to make target `include/Make/Grass.make'. Stop. But at C:\osgeo4w\usr\src\grass-6.4.svn_src_snapshot_2010_07_18\include\Make I realized that the files were with a different name Platform.make.in Grass.make.in So I renamed them Then I rerun package.sh and I obtained this: Wed Jul 21 16:05:42 GMTDT 2010: STARTING configure 0 [main] sh 2948 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 3140 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump configure: error: can not find sources in or .. Any idea of what might be wrong? or be missing? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5298 (20100721) __ 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/mailman/listinfo/grass-user
[GRASS-user] Compiling a Python-Script in Windows
Greetings I have developed a Python script for GRASS in Linux, that calls a external module, and now I have tried to compile it in Windows using this command: make -C callexternalmeanshift.py MODULE_TOPDIR=/c/osgeo4w/apps/grass/grass-6.4.0svn but I get this error: Traceback (most recent call last): File C:/osgeo4w/apps/grass/grass-6.4.0svn/etc/wxpython/wx gui.py, line 473, in OnRunCmd self.goutput.RunCmd(cmd, switchPage=False) File C:\osgeo4w\apps\grass\grass-6.4.0svn\etc\wxpython\gu i_modules\goutput.py, line 354, in RunCmd menuform.GUI().ParseCommand(cmdlist, parentframe=self) File C:\osgeo4w\apps\grass\grass-6.4.0svn\etc\wxpython\gu i_modules\menuform.py, line 1863, in ParseCommand xml.sax.parseString(getInterfaceDescription(cmd[0]).decode(e nc).split('\n',1)[1].replace('', '?xml version=1.0 encoding=utf-8?\n', 1).encode(utf-8), File C:\osgeo4w\apps\grass\grass-6.4.0svn\etc\wxpython\gu i_modules\menuform.py, line 1802, in getInterfaceDescription raise IOError, _(Unable to fetch interface description for command '%s'.) % cmd IOError : Unable to fetch interface description for command 'callexternalmeanshift.py'. 'callexternalmeanshift' is not recognized as an internal or external command, operable program or batch file. In Linux it works...Can Anyone help me? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5278 (20100714) __ 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/mailman/listinfo/grass-user
[GRASS-user] object-oriented segmentation
Greetings I believe few persons have already sent a message similar to this one but here goes: I'm looking for a object-opriented segmentator in GRASS: has anyone ever programmed or addpated some module/processing chain in order to get a object-oriented segmentator? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5275 (20100713) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Error while using r.in.gdal for USGS landsat7 images
r.info displays normal information such as Extent, resolution etc... So it seems all ok Hamish wrote: António wrote: I'm using r.in.gdal to import Landsat7 images downloaded from USGS website. These landsat images are distributed in UTM-WGS84 for the specific zone. in this case I'm working with images from Brazil (WGS84 UTM ZONE 24S). So first I created an Location based on EPSG:32724. Here goes the gdalinfo command for the image: /Driver: GTiff/GeoTIFF Files: L5216065_06520070405_B10.TIF Size is 8171, 6991 Coordinate System is: PROJCS[WGS 84 / UTM zone 24N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,-39], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32624]] Origin = (443999.9941792,-694500.000) Pixel Size = (30.000,-30.000)/ As can be seen it's North instead of South. Nevertheless, r.in.gdal imports it and it even says: Projection of input dataset and current location appear to match But when I try to display I get an error window stating: Execution failed: 'd.rast -o map=imp...@permanent --q' Details: Error: Unable to open raster map imp...@permanent And nothing else is stated what does g.list rast say? how about r.info? Hamish __ Information from ESET NOD32 Antivirus, version of virus signature database 5261 (20100708) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5261 (20100708) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Error while using r.in.gdal for USGS landsat7 images
Hello Markus I have just tested that and it seems to be the best solution. But in this case I didn't know that it was all set to north. Markus Neteler wrote: 2010/7/7 António Rocha antonio.ro...@deimos.com.pt: Hi there, I'm using r.in.gdal to import Landsat7 images downloaded from USGS website. These landsat images are distributed in UTM-WGS84 for the specific zone. in this case I'm working with images from Brazil (WGS84 UTM ZONE 24S). So first I created an Location based on EPSG:32724. Here goes the gdalinfo command for the image: /Driver: GTiff/GeoTIFF Files: L5216065_06520070405_B10.TIF Size is 8171, 6991 Coordinate System is: PROJCS[WGS 84 / UTM zone 24N, ... As can be seen it's North instead of South. Nevertheless, r.in.gdal imports it and it even says: Projection of input dataset and current location appear to match but only appear I think - it is a garbage in/garbage out problem. I would assign the correct EPSG code with gdal_translate -a_srs ... badfile.tif goodfile.tif then reimport. To catch all possible (rare) nonsense case in r.in.gdal I consider to be overkill. Markus __ Information from ESET NOD32 Antivirus, version of virus signature database 5261 (20100708) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5261 (20100708) __ 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/mailman/listinfo/grass-user
[GRASS-user] Error while using r.in.gdal for USGS landsat7 images
Hi there, I'm using r.in.gdal to import Landsat7 images downloaded from USGS website. These landsat images are distributed in UTM-WGS84 for the specific zone. in this case I'm working with images from Brazil (WGS84 UTM ZONE 24S). So first I created an Location based on EPSG:32724. Here goes the gdalinfo command for the image: /Driver: GTiff/GeoTIFF Files: L5216065_06520070405_B10.TIF Size is 8171, 6991 Coordinate System is: PROJCS[WGS 84 / UTM zone 24N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,-39], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32624]] Origin = (443999.9941792,-694500.000) Pixel Size = (30.000,-30.000)/ As can be seen it's North instead of South. Nevertheless, r.in.gdal imports it and it even says: Projection of input dataset and current location appear to match But when I try to display I get an error window stating: Execution failed: 'd.rast -o map=imp...@permanent --q' Details: Error: Unable to open raster map imp...@permanent And nothing else is stated Shall I add this issue to the Error/debug trac list? Or how can I fix this? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5258 (20100707) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Retrieve pixel resolution using Python, Script
Thanks Eloi but this only gives me the number of pixels (number of lines and columns) not the pixel resolution. Is there a retrieve this information in a Python script (e.g. using python GDAL)? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5236 (20100629) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Retrieve pixel resolution using Python Script
Hello Mark I know that I can retrieve that information with gdalinfo. But I need to have it on a variable and get/use it in a Python Script. Mark Seibel wrote: If you use gdalinfo on an image, it will output the Pixel Size. Hope that helps. Mark 2010/6/28 António Rocha antonio.ro...@deimos.com.pt: Greetings I'm trying to retrieve an image (not yet imported to GRASS) pixel resolution through a Python Script but so far i was unable to figure out a method to do this. Is this possible? I have used g.proj to retrieve its projection and datum but not the Pixel resolution. THanks for your help Best regards, Antonio R. __ Information from ESET NOD32 Antivirus, version of virus signature database 5234 (20100628) __ 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/mailman/listinfo/grass-user __ Information from ESET NOD32 Antivirus, version of virus signature database 5234 (20100628) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5235 (20100628) __ 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/mailman/listinfo/grass-user
[GRASS-user] Re: Compiling winGRASS with OSGEO4w- #2
Hi Luisa and Hamish and the rest While I was reading your log file I realized two things: 1- You only have a problem with libgrass_display.dll 2- GRASS cannot find -lgrass_raster This is the creation of libgrass_display.dll: gcc -shared -o /c/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/dist.i686-pc-mingw32/lib/libgrass_display.6.4.0svn.dll -L/c/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/dist.i686-pc-mingw32/lib -Wl,--export-dynamic,--enable-runtime-pseudo-reloc -L/c/OSGeo4W/apps/gdal-16/lib -L/c/OSGeo4W/lib OBJ.i686-pc-mingw32/cnversions.o OBJ.i686-pc-mingw32/color_list.o OBJ.i686-pc-mingw32/draw.o OBJ.i686-pc-mingw32/draw2.o OBJ.i686-pc-mingw32/get_win.o OBJ.i686-pc-mingw32/ident_win.o OBJ.i686-pc-mingw32/list.o OBJ.i686-pc-mingw32/popup.o OBJ.i686-pc-mingw32/raster.o OBJ.i686-pc-mingw32/raster2.o OBJ.i686-pc-mingw32/setup.o OBJ.i686-pc-mingw32/symbol.o OBJ.i686-pc-mingw32/tran_colr.o OBJ.i686-pc-mingw32/window.o -lgrass_gis -lgrass_datetime -lxdr -liberty -lws2_32-lz -lintl -lgrass_raster -lgrass_pngdriver -lgrass_driver -lgrass_gis -lgrass_datetime -lxdr -liberty -lws2_32-lz -lintl -lfreetype -liconv -lgrass_gis -lgrass_datetime -lxdr -liberty -lws2_32-lz -lintl -lpng -lz -lgrass_psdriver -lgrass_driver -lgrass_gis -lgrass_datetime -lxdr -liberty -lws2_32-lz -lintl -lfreetype -liconv -lgrass_gis -lgrass_datetime -lxdr -liberty -lws2_32-lz -lintl -lgrass_driver -lgrass_gis -lgrass_datetime -lxdr -liberty -lws2_32 -lz -lintl -lfreetype -liconv -lgrass_gis -lgrass_datetime -lxdr -liberty -lws2_32-lz -lintl ../../lib/gis/OBJ.i686-pc-mingw32/dllmain.o \ (cd /c/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/dist.i686-pc-mingw32/lib; ln -f -s libgrass_display.6.4.0svn.dll /c/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/dist.i686-pc-mingw32/lib/libgrass_display.dll) c:\OSGeo4W\bin\..\lib\gcc\mingw32\3.4.5\..\..\..\..\mingw32\bin\ld.exe: cannot find -lgrass_raster collect2: ld returned 1 exit status make[3]: *** can it be some old library or package that was installed? __ Information from ESET NOD32 Antivirus, version of virus signature database 5196 (20100614) __ 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/mailman/listinfo/grass-user
[GRASS-user] Compiling winGRASS with OSGEO4w- #3
Greetings Since Luisa and few others are posting regarding this topic I decided to perform this compilation of GRASS on MS-WINDOWs and I followed the exact steps but made a few updates in Libraries. MsysCORE- 1.0.14-1 bin-utils 2.20.51-1 mingw-make 3.81-20090910 (newer version) mingw-utils 0.4-1 w32api 3.14 and Bison 2.4.1 Instead of compiling a weekly snapshot I choose RC6 version. But I got some errors like thiese oens: 0 [main] sh 2008 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 1248 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump checking host system type... 0 [main] sh 1264 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump checking for source directory... 0 [main] sh 1848 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump checking for termios.h... ./configure: line 4110: 224 Segmentation fault (core dumped) ( eval $ac_try ) 25 grep: conftest.out: No such file or directory checking for f2c.h... 0 [main] sh 2020 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump checking for sys/timeb.h... yes 0 [main] sh 3868 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump checking for socket... no checking for socket in -lsocket... configure: warning: *** Unable to locate socket() function. no *And the final error was this (which terminated the process)* checking for location of PostgreSQL library... checking for PQsetdbLogin in -lpq... ./configure: line 9206: 3856 Segmentation fault (core dumped) ( eval $ac_link ) 25 no checking for PQsetdbLogin in -lpq... no checking for PQsetdbLogin in -lpq... no checking for PQsetdbLogin in -lpq... configure: error: *** Unable to locate PostgreSQL library. no *Then I reinstalled postgres from osgeo4w and FFTW-devel (in a previous run I had an error with missing fftw-devel)* After that I got an error in checking for Tcl_Init in -ltclUNKNOWN... configure: error: *** Unable to locate Tcl library. no *I had to reinstall tcltk usiong osgeo4w.* Then: checking for location of External PROJ.4 includes... checking for proj_api.h... ./configure: line 7069: 2916 Segmentation fault (core dumped) ( eval $ac_try ) 25 grep: conftest.out: No such file or directory yes checking External PROJ.4 version... ./configure: line 7120: 3712 Segmentation fault (core dumped) ( ./conftest; exit ) 2/dev/null configure: error: *** Could not determine External PROJ.4 version. *I had to reinstall proj4 using osgeo4w* After that it compiled with the following errors (installation log at: http://pastebin.com/qvgwE7iP): Errors in: /osgeo4w/usr/src/grass-6.4.0RC6/lib/db/sqlp /osgeo4w/usr/src/grass-6.4.0RC6/lib/gpde /osgeo4w/usr/src/grass-6.4.0RC6/lib/ogsf /osgeo4w/usr/src/grass-6.4.0RC6/lib/nviz /osgeo4w/usr/src/grass-6.4.0RC6/db/drivers/dbf /osgeo4w/usr/src/grass-6.4.0RC6/db/drivers/postgres /osgeo4w/usr/src/grass-6.4.0RC6/db/drivers/ogr /osgeo4w/usr/src/grass-6.4.0RC6/display/d.grid /osgeo4w/usr/src/grass-6.4.0RC6/display/d.where /osgeo4w/usr/src/grass-6.4.0RC6/general/g.mremove /osgeo4w/usr/src/grass-6.4.0RC6/general/g.proj /osgeo4w/usr/src/grass-6.4.0RC6/general/g.region /osgeo4w/usr/src/grass-6.4.0RC6/general/g.setproj /osgeo4w/usr/src/grass-6.4.0RC6/general/manage/cmd /osgeo4w/usr/src/grass-6.4.0RC6/ps/ps.map /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.external /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.gwflow /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.horizon /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.los /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.mapcalc /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.proj.seg /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.sun2 /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.sunmask /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.to.rast3 /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.to.rast3elev /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.univar2 /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.in.gdal /osgeo4w/usr/src/grass-6.4.0RC6/raster/r.out.gdal /osgeo4w/usr/src/grass-6.4.0RC6/raster3d/r3.cross.rast /osgeo4w/usr/src/grass-6.4.0RC6/raster3d/r3.in.ascii /osgeo4w/usr/src/grass-6.4.0RC6/raster3d/r3.in.v5d /osgeo4w/usr/src/grass-6.4.0RC6/raster3d/r3.gwflow /osgeo4w/usr/src/grass-6.4.0RC6/raster3d/r3.mkdspf /osgeo4w/usr/src/grass-6.4.0RC6/raster3d/r3.out.ascii /osgeo4w/usr/src/grass-6.4.0RC6/raster3d/r3.out.v5d /osgeo4w/usr/src/grass-6.4.0RC6/raster3d/r3.out.vtk /osgeo4w/usr/src/grass-6.4.0RC6/raster3d/r3.stats /osgeo4w/usr/src/grass-6.4.0RC6/raster3d/r3.to.rast /osgeo4w/usr/src/grass-6.4.0RC6/raster3d/base /osgeo4w/usr/src/grass-6.4.0RC6/vector/v.proj /osgeo4w/usr/src/grass-6.4.0RC6/vector/v.to.rast3 /osgeo4w/usr/src/grass-6.4.0RC6/vector/v.vol.rst /osgeo4w/usr/src/grass-6.4.0RC6/vector/v.out.ogr /osgeo4w/usr/src/grass-6.4.0RC6/vector/v.in.ogr /osgeo4w/usr/src/grass-6.4.0RC6/visualization/nviz /osgeo4w/usr/src/grass-6.4.0RC6/visualization/nviz2/cmd To check the errors I went to /raster/r.in.gdal and run make and I got: make: *** No rule to make
[GRASS-user] Re: Compiling winGRASS with OSGEO4w- #2
Greetings In my previous email I used GRASS6.4.0RC6 to compile with OSGEO4w now, following Luisa email, I tried with latest snapshot and I got errors in the following modules: Errors in: /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/lib/db/sqlp /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/db/drivers/dbf /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/db/drivers/postgres /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/db/drivers/ogr /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.barscale /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.colorlist /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.colors /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.colortable /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.erase /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.extend /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.extract /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.font /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.frame /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.geodesic /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.graph /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.grid /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.his /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.histogram /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.info /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.legend /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.linegraph /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.mapgraph /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.measure /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.menu /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.mon/cmd /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.mon/pgms /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.nviz /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.paint.labels /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.path /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.profile /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.rast /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.rast.arrow /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.rast.num /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.rgb /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.rhumbline /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.save /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.text.new /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.thematic.area /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.title /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.vect /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.vect.chart /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.what.rast /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.what.vect /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.where /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/display/d.zoom /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/general/g.mapset /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/imagery/i.ask /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/raster/r.le/r.le.setup /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/raster/r.le/r.le.trace /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/raster/r.mapcalc /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/raster/r.profile /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/raster/wildfire/r.spread /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/vector/v.label /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/vector/v.digit /osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_06_12/vector/v.label.sa It seems that I'm missing the following DLL: - libgrass_display.dll Just like you Luisa. Anyone has an idea of what might be happening here? Thank you Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5197 (20100615) __ 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/mailman/listinfo/grass-user
[GRASS-user] Re: Compiling GRAS on Windows- Log and error
Hi Luisa I'm no expert in OSGEO4W but I believe you should 1- shut down Msys and reboot the PC 2- Try again I believe that you might have 1 problem: You don't have an include folder. I don't know how can you add it to your GRASS installation by using OSGEO4W. But that is not one of the steps in that compiling for MS-WIndows. __ Information from ESET NOD32 Antivirus, version of virus signature database 5195 (20100614) __ 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/mailman/listinfo/grass-user
[GRASS-user] Re: Compiling GRAS on Windows- Log and error
Luisa, I believe Glynn and Martin can help you out on this but, I believe you will have to change grass/include/make (Platform.make ) in order to be adequate with your installation. Is this right Glynn? If so, shouldn't this be added to that Wiki? (if it makes sense, I volunteer for that) By the way, shoudl Luisa change Platform.make from include/make/Platform.make in source code or in the GRASS-desktop that she installed with OSGEO4W? __ Information from ESET NOD32 Antivirus, version of virus signature database 5196 (20100614) __ 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/mailman/listinfo/grass-user
[GRASS-user] Re: Import GRIB files to GRASS
Hi all Well, since I'm only evaluating (for now) how to import GRIB's into GRASS, I wanted to test the r.region solution. Can anyone describe me what to do regarding that suggestion? Hamish, about the patch, I'm using GDAL1.7 so I cannot apply your patch. I'm going to reply on this in GDAL bug tracker. Thanks It is so interesting to see vastly different approaches to solve a similar problem. I'm downloading NOAA forecast data which is in GRIB format. I run the files through degrib and then parse the shape files with a C program using the shapelib libraries and create a new shape file which I import into GRASS with v.in.ogr. I plot the values with d.vect.thematic. The resulting maps can be viewed at: http://www.ars.usda.gov/npa/marc/heatstress That reminds me, I originally set up a location for each region with an eye forward on multithreading the map generation phases. I forgot the reason in the middle of the project and created several regions in one location. I had to revert back to a location for each region when I implemented multithreading but reduced the total map generation time from several hours to 15 minutes. :-) Randy __ Information from ESET NOD32 Antivirus, version of virus signature database 5189 (20100611) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Import GRIB files to GRASS
Hello Hamish. Yest, I used a NOAA's tool to verify it and it's misplaced. I will take a look at those tickets that you mention. About this: /There is a cell-registration bug in GDAL: (grid vs cell-center convention) cells are all offset to the east and south by half a cell. You can fix this by hand with the r.region http://grass.osgeo.org/grass64/manuals/html64_user/r.region.html module. /How can I use r.region to fix this? Hamish wrote: António wrote: I have one question regarding Hamish's link. /There is a cell-registration bug in GDAL: (grid vs cell-center convention) cells are all offset to the east and south by half a cell. You can fix this by hand with the r.region http://grass.osgeo.org/grass64/manuals/html64_user/r.region.html module. /Is this still happening, this error?/ well, the linked GDAL bug tickets are still open: https://trac.osgeo.org/gdal/ticket/2637 https://trac.osgeo.org/gdal/ticket/2550 so I guess so. I'd encourage you to verify for yourself if your data is offset or not, and post that verification (or not) to the appropriate bug report with as much detail as you can; and if the patch attached to the ticket fixes it for you (or not). Hamish __ Information from ESET NOD32 Antivirus, version of virus signature database 5181 (20100608) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5185 (20100609) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Import GRIB files to GRASS
Hi Thomas and GRASS-people THanks for your scripts Thomas I have one question regarding Hamish's link. /There is a cell-registration bug in GDAL: (grid vs cell-center convention) cells are all offset to the east and south by half a cell. You can fix this by hand with the r.region http://grass.osgeo.org/grass64/manuals/html64_user/r.region.html module. /Is this still happening, this error?/ / Thomas Adams wrote: António I have does this; I think it's helpful to first extract the fields you need and group them together, otherwise you will probably end up with a lot of unwanted fields being imported. I have written a Perl script to do this if you want it. I'm sure you would need to edit it to meet your needs. Regards, Tom Hamish wrote: António wrote: I'm planning on using some global meteo datasets that are distributed in GRIB format. has anyone ever tried to parse GRIB NCEP files to GRASS? maybe this helps- http://grass.osgeo.org/wiki/GRIB Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user __ Information from ESET NOD32 Antivirus, version of virus signature database 5181 (20100608) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Select random cells from map
Thanks Hamish I thought that random cells are retrieved from valid map cells. Or not? Antonio Hamish wrote: António wrote: I need to randomly select non-null cells from a given map. I tried r.random.cells but I cannot define a maximum number random cells and some selected random cells are not valid values in my map. is there any other function to randomly select cells from a raster? Hamish: r.to.vect to get a vector area coverage of your raster (do not use the -s flag), then v.random.cover from wiki addons then v.to.rast, and finally r.mapcalc to transfer values from source map to random cell map? (probably we/I should make a new wrapper script to do that automatically) or you might skip the last steps and just use v.random.cover's raster= option to sample the raster map directly. probably the best idea is to fix r.random's cover= option to loop until it has enough points. also if you will be doing a lot of geostats it is probably worth your while to learn about the GRASS-R interface: http://grass.osgeo.org/wiki/R_statistics ... then the world is open to you. António wrote: Thanks for your answer. I tried r.random with my input data and the number of cells and it worked but I do have one question: I didn't understand the objective of cover_map. I used the map from I want to extract random cells as input but what is the objective of cover_map? r.random creates random points all over the current rectangular region block. Adding a cover map makes those random points only occur within the irregular area covered by the cover map. (but r.random cover= is a bit buggy right now; see new ticket in the bug trac'er) Hamish __ Information from ESET NOD32 Antivirus, version of virus signature database 5177 (20100606) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5178 (20100607) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Select random cells from map
Hi. I missed that r.random :) It seems to do what I need. I just need to understand a bit more what is the COVER parameter Thanks nikos Nikos Alexandris wrote: António Rocha: Hi there I need to randomly select non-null cells from a given map. I tried r.random.cells but I cannot define a maximum number random cells and some selected random cells are not valid values in my map. is there any other function to randomly select cells from a raster? Doesn't r.random (not r.random.cells) work? I think it works on/for non- NULL cells. In case r.random.cells and its capabilities are required, couldn't you just create a non-NULL MASK and limit its operation as you wish? Nikos __ Information from ESET NOD32 Antivirus, version of virus signature database 5163 (20100601) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5163 (20100601) __ 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/mailman/listinfo/grass-user
[GRASS-user] Using a polygon as a Mask
Greetings I have polygon that I want to use to define an for which I want to calculate some indicators (almost as a mask). Is it possible to use it directly in mapcalculator or do I have to convert it to RASTER and only then I can use it in mapcalculator as a mask? Thanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5149 (20100527) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Using a polygon as a Mask
Thank you Christian What I am affraid is that when I convert Vector to RASTER I will obtain a mask that is not exacly the polygon in the boundaries of the polygon. For instance, if a pixel is cuttted by the polygons boundaries is it considered part odf the polygon or not? Thanks Antonio Christian Schwartze wrote: There is no way to use a vector geometry in mapcalculator (or mapcalc) directly, hence, you have to convert it to raster before. Christian. -- From: António Rocha antonio.ro...@deimos.com.pt Sent: Thursday, May 27, 2010 2:08 PM To: GRASS user list grass-user@lists.osgeo.org Subject: [GRASS-user] Using a polygon as a Mask Greetings I have polygon that I want to use to define an for which I want to calculate some indicators (almost as a mask). Is it possible to use it directly in mapcalculator or do I have to convert it to RASTER and only then I can use it in mapcalculator as a mask? Thanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5149 (20100527) __ 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/mailman/listinfo/grass-user __ Information from ESET NOD32 Antivirus, version of virus signature database 5149 (20100527) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5150 (20100527) __ 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/mailman/listinfo/grass-user
Re: [GRASS-user] Landsat Stack import and export
Greetings Regarding the possibility of exporting an GROUp, using r.in.gldal, the groups are not listed in that menu. So the user cannot see the available groups. It's still processing but IO have one question: The processing bar restarted at least 3 times, what is happening? Because it file, now, it seems that it's not growing. Is it doing some file-check or something? Thanks Antonio Markus Metz wrote: 2010/5/7 Markus Neteler nete...@osgeo.org: 2010/5/7 António Rocha antonio.ro...@deimos.com.pt: Greetings I've a Landsat image with 6 bands in Tif format. I need to import into grass, do some processing and export the same band in a Stack. The import is easy using r.in.gdal but my difficulty is to export 6 bands in one single file. Is it possible using GRASS? To my knowledge no. Sure it is. Create a group in GRASS, export the group (r.out.gdal input=group_name) to a file format that supports multiple bands, e.g. GeoTIFF. Markus M __ Information from ESET NOD32 Antivirus, version of virus signature database 5094 (20100507) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5100 (20100510) __ 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/mailman/listinfo/grass-user
[GRASS-user] v.out.ogr and region
Greetings I'm exporting data from GRASS to Gtif and Shapefiles. While r.out.gdal obays my Region definition v.out.org no. It exports all extent not the specified region. Wasn't it suppose to export only for defined region? Thanks __ Information from ESET NOD32 Antivirus, version of virus signature database 5102 (20100510) __ 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/mailman/listinfo/grass-user
[GRASS-user] Landsat Stack import and export
Greetings I've a Landsat image with 6 bands in Tif format. I need to import into grass, do some processing and export the same band in a Stack. The import is easy using r.in.gdal but my difficulty is to export 6 bands in one single file. Is it possible using GRASS? Thanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5093 (20100506) __ 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/mailman/listinfo/grass-user