[GRASS-user] Limits on computation
Greetings Is there any defined limits on computation? I mean, for instance I'm creating a mosaic. Is there any limits on size of images or size of data that is known in GRASS? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Limits on computation
Hi Glynn Thanks Just one question aboht that limit of number of cells: is that limit applicable to CELL or DCELL? Thanks Monica On Mon, May 30, 2011 at 2:17 PM, Glynn Clements gl...@gclements.plus.comwrote: Monica Buescu wrote: Is there any defined limits on computation? I mean, for instance I'm creating a mosaic. Is there any limits on size of images or size of data that is known in GRASS? You may encounter problems with maps with over 2^31 (2147483648) cells or files over 2GiB. The 2GiB file limit shouldn't be an issue for GRASS 7 on Unix with large file support (LFS) enabled, but LFS doesn't work on Windows and is incomplete in 6.x. Some modules have problems with large amounts of data, either requiring too much memory or taking forever to run. In these cases, there isn't any clear limit; it depends upon how much RAM is available, whether the system is 32-bit or 64-bit, and how long is too long. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Error while exporting
I'm expericing a problem while trying to export a reclasses (r.reclass) raster map: I'm getting this error: Exporting to GDAL data type: Byte Checking GDAl data type and nodata value Selected GDAL datatype does not cover data range. GDAL datatype: Byte, range: 0 - 255 Raster map full_LULC@Regional range: -2.14748e+009 - 20 Raster export results in data loss. Input raster map contains cells with NULL-value (no-data). The value 255 will be used to represent no-data values in the input map. You can specify a nodata value with the nodata option. ERROR: Raster export aborted. But If I try to export my original file (before reclass) It exports just fine. My reclass file is like this: 190 191 = 1 14 15 16 20 21 22 = 2 11 12 13 = 3 32 40 41 42 50 60 70 90 91 92 100 101 102 = 4 30 31 110 130 131 132 133 134 135 136 = 5 120 140 141 142 143 144 145 = 6 150 151 152 153 = 7 230 = 8 200 201 202 203 = 9 160 161 162 170 180 181 182 183 184 185 186 187 188 = 10 210 = 11 211 245= 20 * = NULL What might be happening? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Very Quick Question about Geometric Correction module (WiNGRASS6.4.0stable)
Greetings Very Quick Question: my students are running Geometric Correction in Windows (6.4.0 stable release from Sep2010). It is running but the interfaces seems to be freezed (at the top of the window it's saying (Not Responding). But it has processed the data. Question 1: can it be a problem due to the fact that Target and Source images are in different mapsets? (guess not) Question 2: Why isn't launching the final panel saying all Geometric Correction concluded (or something silimar) Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Very Quick Question about Geometric Correction module (WiNGRASS6.4.0stable)
Hi markus The thing is for smaller images is OK but for bigger it runs but that panel is not displayed that bug is related with 6.5 and 7 The georectifier is broken in the most recent svn versions of GRASS 6.5 and 7 (maybe earlier ones too, but I don't know). Starting up the wizard works fine all the way to the end when it is supposed to launch the interface for setting and managing GCP's. This generates an error: I have tried with 6.4.1 (in Windows) and I get the same problem... On Fri, Apr 29, 2011 at 10:43 AM, Markus Neteler nete...@osgeo.org wrote: On Fri, Apr 29, 2011 at 11:30 AM, Monica Buescu monicabuescu1...@gmail.com wrote: Greetings Very Quick Question: my students are running Geometric Correction in Windows (6.4.0 stable release from Sep2010). It is running but the interfaces seems to be freezed (at the top of the window it's saying (Not Responding). But it has processed the data. Question 1: can it be a problem due to the fact that Target and Source images are in different mapsets? (guess not) Question 2: Why isn't launching the final panel saying all Geometric Correction concluded (or something silimar) While I cannot say how it works under Windows, I would like to signal this #1257 Georectifier is broken: fixed in 6.4.1 If you can, please update to GRASS 6.4.1 which has received a series of wxGUI bugfixes. See also http://trac.osgeo.org/grass/wiki/Release/6.4.1-News Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Geometric Correction- Freeze
Greetings Markus (I changed the subject since it's happening in 6.4.1) There is no error. It just freezes. It creates the bands and georrectigfies but the panels, windows are freezed. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Warning- Not enough room in history file for command line (truncated).
Greetings I'm running a function (r.series) and I'm getting this warning at the end: Not enough room in history file for command line (truncated). How can I avoid this? as far as I can see there is no problem but... I would like to avoid this Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Adding labels to Classes
Greetings I have a Land Cover map and I want to substitute Land Cover codes by numbers. Besides r.reclass (1 = 1 Agriculture ; 2 = 2 Forest 3 = 3 Water) what other method does exists in GRASS? Thanks ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Retrieve pixel value for a group of images
Greetings I have created a group of 10 images (i.group) and I wanted to know how can I retrieve a pixel value for all images instead of using query raster map tool wherre I can only retrieve a pixel per image? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Manual page on how to create/edit/delete Mapsets and Locations
Greetings Where can I find, if available, a webpage or manual page on how to define/create/edit/delete Locations and Mapsets? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Add a color to a vector
Greetings I have produced in winGRASS 6.4.0 stable a vector (by converting from raster to vector) and I want to add a color. how can I do that inside a Python Script? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Calling GRASS scripts from GRASS C modules
Greetings I'm building a segmentation and statistical image analysis module (in C) and I want to call a grass script. Is it possible inside a GRASS C module? If so, How can I do that? (or can anyone give me an online example= Thanks Best regards, Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Question on using r.resamp.stats
Greetings a few messages ago I saw that someone advised to use r.resamp.stats to rescale for a coarser scale. Now I have one question: In mapsetX I have a raster (named rasterA) with Spatial resolution 10 and 10. Now in mapset Y my region spatial resolution is 50 and 50. Question: Having mapset Y active, If I apply r.resamp.stats to a raster map in a different mapset (mapset X) will it apply some kind of processing before it applies the r.resamp.stats? (like nearest neighbor or something)? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Using temp files in Grass Script
Greetings I'm doing a GRASS script and I have some intermediate processes (e.gr.mapcalc) where I produce raster maps that are not the final product (a kind of an intermediate product) and that i delete at the end. Since I have to produce a lot my question is: Shall I use tempfile function to create and produce them or I shall create these rasters as normal raster files? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] calculating average value over a base map
Hi Dmistrij I think it was worked. I just need to figure out something with r.statistics. Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] calculating average value over a base map
Greetings Christian Unfortunely it only works with CELL values. For instance, r.median uses more than just CELL. Is there any other alternative? Beucase for mode and median there is specific functions. Thamls Best regards, Monica On Mon, Feb 28, 2011 at 3:55 PM, Christian Schwartze christian.schwar...@uni-jena.de wrote: Monica, try to use r.statistics with the average function/method. Pass you maps through the cover= and base= parameters... Christian. *From:* Monica Buescu monicabuescu1...@gmail.com *Sent:* Monday, February 28, 2011 4:27 PM *To:* GRASS user list grass-user@lists.osgeo.org *Subject:* [GRASS-user] calculating average value over a base map Greetings I have a map with values and a base map created with r.clump. How can I calculate the average of values that falls inside each cluster defined by r.clump? Thanks Best regards, Monica -- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] calculating average value over a base map
Hello Christian. I believe I didn't quite understand what you have said but I think that is not the way. My basemap is a r.clump output in CELL format that ranges from 1 til 1700 or something my cover map is NDVI that ranges from -1 (near) til 1 (near) and is a FCELL. My Objective. calculate for each clump, the average of NDVI pixels If I do: r.statistics -c base=new_clump cover=2007_ndvi method=average I get: ERROR: This module currently only works for integer (CELL) maps Is there any alternative to do this? Thank you Christian and the rest of grass mailing list Monica On Mon, Feb 28, 2011 at 5:08 PM, Christian Schwartze christian.schwar...@uni-jena.de wrote: If you want to deal with raster value types other than CELL, you can use the –c flag in r.statistics in order to process the labels. Therfore, you have to copy your cat values to the category labels first... Regards, Christian. *From:* Monica Buescu monicabuescu1...@gmail.com *Sent:* Monday, February 28, 2011 5:16 PM *To:* Christian Schwartze christian.schwar...@uni-jena.de *Cc:* GRASS user list grass-user@lists.osgeo.org *Subject:* Re: [GRASS-user] calculating average value over a base map Greetings Christian Unfortunely it only works with CELL values. For instance, r.median uses more than just CELL. Is there any other alternative? Beucase for mode and median there is specific functions. Thamls Best regards, Monica On Mon, Feb 28, 2011 at 3:55 PM, Christian Schwartze christian.schwar...@uni-jena.de wrote: Monica, try to use r.statistics with the average function/method. Pass you maps through the cover= and base= parameters... Christian. *From:* Monica Buescu monicabuescu1...@gmail.com *Sent:* Monday, February 28, 2011 4:27 PM *To:* GRASS user list grass-user@lists.osgeo.org *Subject:* [GRASS-user] calculating average value over a base map Greetings I have a map with values and a base map created with r.clump. How can I calculate the average of values that falls inside each cluster defined by r.clump? Thanks Best regards, Monica -- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Error while using variance in r.series function
Greetings I'm using r.series with variance option and I'm getting values that are not according to var function in MATLAB or var function in EXCEL. For average and count the result is exacly the expected but for variance (e.g. bellow is for a pixel [ 73 116 72 81 75 136 101 112 94 67 ] I get 532.011 in EXCEL and in MATLAB and in GRASS i get 478.81 What is happening? Was this already reported? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] renaming a raster map in a Python Script
Greetings I have a GRASS Python Script where I create a raster map and I want to rename it. How can I rename a raster map in a GRASS Python Script? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Check if a mapset exists- in Python Script
Hi Glynn I tried grass.script.mapsets but I get this error: line 59, in main var=grass.mapsets(False) AttributeError: 'module' object has no attribute 'mapsets' previously I had this: import grass.script as grass Am I doing something wrong? Thank you Glynn Monica On Mon, Jan 31, 2011 at 11:38 AM, Glynn Clements gl...@gclements.plus.comwrote: Monica Buescu wrote: How can I check if a certain mapset exists in my active location in a Python Script? grass.script.mapsets(False) returns a list of mapsets in the current location. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Check if a mapset exists- in Python Script
Greetings How can I check if a certain mapset exists in my active location in a Python Script? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Produce KML/KMZ files in GRASS
Greetings I would like to know if it's possible to produce KML/KMZ files from a raster file in GRASS? If not, what tool do you advise? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Calculating slope for raster maps on pixel-based methodology
Greetings I have a bunch of raster maps (5 as a matter of fact) and need to calculate the slope for each pixel. I mean, I need to pick-up all 5 values for each pixel and calculate Slope for that pixel. Now, What I want is to do this for all pixels and obtain a map that represents the slope for each pixel based on 5 raster maps. Is this possible in GRASS? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Creating patches in a raster file
Greetings (It seems that I'm mail bombing to this list but I have a few questions about GRASS usage). I have a classification map and I want to create segments/patches aggregated by its classigfication and assign an patch ID. Let me demonstrate: INPUT: 1 1 2 2 4 4 5 5 5 1 2 2 2 4 5 5 5 1 3 3 2 2 2 5 OUTPUT 1 1 3 3 5 5 6 6 2 1 3 3 3 5 6 6 2 1 4 4 3 3 3 6 I mean, instead of having the classification values I have an Pactch ID label. Is this possible in GRASS' Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Not possible to display vector data
Greetings I have a vector map with the following characteristics (v.info) v.info test01_LU ++ | Layer: test01_LU | | Mapset: Regional | | Location:Localize | | Database:D:\data | | Title: | | Map scale: 1:1 | | Map format: native | | Name of creator: monica | | Organization: | | Source date: Tue Nov 23 11:19:55 2010 | || | Type of Map: vector (level: 2) | | | | Number of points: 0 Number of areas: 76832 | | Number of lines:0 Number of islands:7949 | | Number of boundaries: 172475 Number of faces: 0 | | Number of centroids:76816 Number of kernels:0 | | | | Map is 3D: No | | Number of dblinks: 1 | | | | Projection: Transverse Mercator | | N: 276083.7674S: -300404.804 | | E: 162129.0811W: -119191.4075 | | | | Digitization threshold: 0 | | Comments: | | | ++ (Tue Nov 23 11:41:42 2010) Command finished (0 sec) When I try to display it in a Map Window my wxGUI freezes and it gets non responsive. What might be happen? DO you suggest to wait? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] help, manual and instructions on how to use legend
Greetings I'm trying to find an help page/wiki on use of legend but i'm not able to find anything Can anyone indicate me one? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error when importing data in Windows with Python Script
Hello martin and Grass list as I replied, I'm not running scripts like that. For instance, I have a python script named r.in.landsat I run it as r.inlandsat input=C:\Lsat5_10.tif But I'm not able to run this because of \ it's not being recognized as a valid character. I have just installed WinGRASS nightly version from last week and It's still not available. Is this expected? Thanks On Mon, Nov 8, 2010 at 11:04 AM, Martin Landa landa.mar...@gmail.comwrote: Hi, 2010/11/8 Monica Buescu monicabuescu1...@gmail.com: I have just downloaded and installed WinGRASS r44203 from here (http://josef.fsv.cvut.cz/wingrass/grass64/) and I tried to run a Python Script that uses a file path in Windows and still, it's not using the Windows \ file-paths bars. Can anyone confirm me that this is still not working on this version? just to be sure, you are running script from wxGUI command prompt, e.g. c:\test.py c:\path\to\file ? Martin -- Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error when importing data in Windows with Python Script
Hello Martin I have just downloaded and installed WinGRASS r44203 from here ( http://josef.fsv.cvut.cz/wingrass/grass64/) and I tried to run a Python Script that uses a file path in Windows and still, it's not using the Windows \ file-paths bars. Can anyone confirm me that this is still not working on this version? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error when importing data in Windows with Python Script
Hello Martin I'm running my Python script from its Modules GUI. I insert the path, by-hand... On Mon, Nov 8, 2010 at 11:04 AM, Martin Landa landa.mar...@gmail.comwrote: Hi, 2010/11/8 Monica Buescu monicabuescu1...@gmail.com: I have just downloaded and installed WinGRASS r44203 from here (http://josef.fsv.cvut.cz/wingrass/grass64/) and I tried to run a Python Script that uses a file path in Windows and still, it's not using the Windows \ file-paths bars. Can anyone confirm me that this is still not working on this version? just to be sure, you are running script from wxGUI command prompt, e.g. c:\test.py c:\path\to\file ? Martin -- Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Use of r.mode
Hi there, I have created a raster file from r.mode. It's a big image (is true) and it might me heavy but I'm not being able to get r.info information... In a 4GB RAM memory machine I have r.info raster and after 9 minutes It didn't printed anything. So I'm not being able to retrieve information from this raster. Any suggestion? Shall I create a copy of this file? and is this normal? Thanks Monica On Thu, Oct 21, 2010 at 11:32 AM, Monica Buescu monicabuescu1...@gmail.comwrote: Well, that's not good. Yes, maybe a problem with disk space or memory? l'll have to leave to someone more knowledgeable to answer here. Thanks Micha Probaby in this case is Memory space. Because for Smaller images I don't get any problem and bigger images I'm getting this problem. As anyone experienced such a problem? (GRASS takes too long to display an image or apply r.mapcalc? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] using r.mode- Can different resolutions be a problem?
Thanks Markus. By the way, just one last question: If I have pixel values that I don't want to apply r.mode (for instance 0), can I transform them to null (using r.mapcalc) and then apply r.mode? Isn't there any issue regarding null values? Thanks On Sat, Oct 30, 2010 at 5:34 PM, Markus Neteler nete...@osgeo.org wrote: It will work of course but not necessarily be statistically significant (think for example at a high resolution base map and a just one-pixel cover map...). So some careful reasoning is recommended. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] using r.mode- Can different resolutions be a problem?
Greetings A couple of days ago I sent a few ermails about r.mode. Now I have one questions: is there any issue associated with the fact that, ocasionally, base and cover map may have different resolution? Is this a problem? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Compiling in Windows OSGEO4W- Error with gcc
Greetings I tried to compile a module with this command make -C i.segment MODULE_TOPDIR=/C/OSGeo4W/usr/src/grass-6.4.svn_src_snapshot_2010_09_04 And I got this: make: Entering directory `/C/osgeo4w/usr/src/i.segment' make OBJ.i686-pc-mingw32 make[1]: Entering directory `/C/osgeo4w/usr/src/i.segment' make[1]: `OBJ.i686-pc-mingw32' is up to date. make[1]: Leaving directory `/C/osgeo4w/usr/src/i.segment' gcc -I/c/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_09_04/dist.i686-pc-mingw32/include -I/c/OSGeo4W/apps/gdal-16/include -I/c/OSGeo4W/include -g -O2 -I/c/OSGeo4W/apps/gdal-16/include -I/c/OSGeo4W/include -DPACKAGE=\grassmods\ -I/c/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_09_04/dist.i686-pc-mingw32/include -o OBJ.i686-pc-mingw32/hist.o -c hist.c /bin/sh: gcc: command not found make: *** [OBJ.i686-pc-mingw32/hist.o] Error 127 make: Leaving directory `/C/osgeo4w/usr/src/i.segment' I have in osgeo4w/bin a gcc.exe file so it shopuld be available. How come it's not found? Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Use of r.mode
Well, that's not good. Yes, maybe a problem with disk space or memory? l'll have to leave to someone more knowledgeable to answer here. Thanks Micha Probaby in this case is Memory space. Because for Smaller images I don't get any problem and bigger images I'm getting this problem. As anyone experienced such a problem? (GRASS takes too long to display an image or apply r.mapcalc? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Use of r.mode
Hello About 20,000,000 cells. That should be OK. Do you know how many different values are actually in the data? Yes... In this case is 6608. In other examples it can be 3000 or 2000. Why? I'd try two things here. First can you change to a smaller region and rerun r.mode and try to display the results? Yes. It worked with a much smaller image. Second: the result of r.mode is a reclass of the original values. So I suppose it has to be recalculated each time you want to display. You can make that reclass map into a permanent raster by running: r.mapcalc output2 = output1 I tried this and after 15 minutes running I gave up. Can this be related with memory restrictions? Thanks Micha for your help Micha ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Use of r.mode
Greetings This new topic is related with I need to calculate mode for a base map that is composed by float values. r.mode base=b...@permanent cover=1107...@permanent output=saida And I get: ERROR: reading r.stats output ERROR: No rules specified This means that I'm not being able to apply this function just because my base map is a Float. Any ideas of how to fix this? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Use of r.mode
Hi Mich Thanks for the reply. I think it has worked but everytime I try to open, my winGRASS goes bananas. Just to sum what I have done 1- Import a xyz ASCII file that is a float with a range of [166297.00 - 20855944.0] (not all values are used) and Rows: 4655 ;Columns: 4520 . 2- I did r.mapcalc int() to this file that gets the range of [166297 20855944]. QUESTION: can I apply int to a FLoat with values as big as this? 3- I did r.mode base=int_base coverage= 1107...@permanent (it's a high-resolution image) output=output1 Ok it runs (it takes some time as expected) Thern if I try to display it freezes... Even If I just try r.info output1 it also freezes. Any thoughts aboutt this? (Thanks Micha) On Tue, Oct 19, 2010 at 2:14 PM, Micha Silver mi...@arava.co.il wrote: Monica Buescu wrote: Greetings This new topic is related with I need to calculate mode for a base map that is composed by float values. r.mode base=b...@permanent cover=1107...@permanent output=saida And I get: ERROR: reading r.stats output ERROR: No rules specified This means that I'm not being able to apply this function just because my base map is a Float. Any ideas of how to fix this? How about r.mapcalc base_integer=int(base) ? The base map should be integer (category ) values, where the mode is calculated for all cells with the same cat value, so it indeed should be an integer raster. Thanks Monica This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Calculating mode value for each polygon
Greetings I have an image that is divided in several (300) polygons and, using a data_map raster, I want to calculate the mode value from this image for each polygon. Is it possible? can anyone give me a tip on this? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Calculating mode value for each polygon
I don't think so.. What I want is to calculate the mode (most frequent value) of a raster file for each polygon that this raster covers. What I will have is the most frequent raster value that falls inside each polygon. Is kind of a univariate statistic of a raster for each polygon. Any ideas? (thanks Micha for your suggestion) On Mon, Oct 18, 2010 at 10:19 AM, Micha Silver mi...@arava.co.il wrote: On 18/10/2010 10:11, Monica Buescu wrote: Greetings I have an image that is divided in several (300) polygons and, using a data_map raster, I want to calculate the mode value from this image for each polygon. Is it possible? can anyone give me a tip on this? Have a look at v.what.rast. That might be what you're looking for. Thanks Monica This mail was received via Mail-SeCure System. ___ grass-user mailing listgrass-u...@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silverhttp://www.surfaces.co.il/ Arava Development Co. +972-52-3665918 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Calculating mode value for each polygon
Uhm That is a pittty. Is anyone planning on developing this? On Mon, Oct 18, 2010 at 12:33 PM, Daniel Victoria daniel.victo...@gmail.com wrote: This should be done by v.rast.stats but, from the manual I don't think it includes mode... On Mon, Oct 18, 2010 at 8:32 AM, Monica Buescu monicabuescu1...@gmail.com wrote: I don't think so.. What I want is to calculate the mode (most frequent value) of a raster file for each polygon that this raster covers. What I will have is the most frequent raster value that falls inside each polygon. Is kind of a univariate statistic of a raster for each polygon. Any ideas? (thanks Micha for your suggestion) On Mon, Oct 18, 2010 at 10:19 AM, Micha Silver mi...@arava.co.il wrote: On 18/10/2010 10:11, Monica Buescu wrote: Greetings I have an image that is divided in several (300) polygons and, using a data_map raster, I want to calculate the mode value from this image for each polygon. Is it possible? can anyone give me a tip on this? Have a look at v.what.rast. That might be what you're looking for. Thanks Monica This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver http://www.surfaces.co.il/ Arava Development Co. +972-52-3665918 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Calculating mode value for each polygon
About the sugfgestion from Luisa, of using r.mode it seems exacly what I want... But instead of using polygons I use a raster (coverage map). For me, its ok I have tested and, I'm obtaining a bit weird error on this: r.mode base...@permanent cover=1107...@permanent output=test_outputv0 ERROR: reading r.stats output ERROR: No rules specified My 11073_B is a CELL satellite image (0-1000) my v0 is a FCELL raster with values that range from 1 til 5 (but I don't have that many values only 30 of them but are dispersed between 1 and 5. Any ideas to fix this? THanks Monica On Mon, Oct 18, 2010 at 12:33 PM, Daniel Victoria daniel.victo...@gmail.com wrote: This should be done by v.rast.stats but, from the manual I don't think it includes mode... On Mon, Oct 18, 2010 at 8:32 AM, Monica Buescu monicabuescu1...@gmail.com wrote: I don't think so.. What I want is to calculate the mode (most frequent value) of a raster file for each polygon that this raster covers. What I will have is the most frequent raster value that falls inside each polygon. Is kind of a univariate statistic of a raster for each polygon. Any ideas? (thanks Micha for your suggestion) On Mon, Oct 18, 2010 at 10:19 AM, Micha Silver mi...@arava.co.il wrote: On 18/10/2010 10:11, Monica Buescu wrote: Greetings I have an image that is divided in several (300) polygons and, using a data_map raster, I want to calculate the mode value from this image for each polygon. Is it possible? can anyone give me a tip on this? Have a look at v.what.rast. That might be what you're looking for. Thanks Monica This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver http://www.surfaces.co.il/ Arava Development Co. +972-52-3665918 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Checking if a raster is in current MAPSET
Greetings I have a Py script where I do some processing. But, I want to limit it to the files that are available in my current mapset. How can I do this? I mean, I just to evaluate if the input is in the current/active MAPSET? Thanks Best regards Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Checking if a raster is in current MAPSET
Hello Glynn Hum it sounds interesting. How can I parameterize it? I mean, what is the order of the parameters? About the current mapset I'm already using it thanks :) PS- Where can I find a list of avaliable functions in grass.python? Monica On Fri, Oct 15, 2010 at 2:59 PM, Glynn Clements gl...@gclements.plus.comwrote: Monica Buescu wrote: I have a Py script where I do some processing. But, I want to limit it to the files that are available in my current mapset. How can I do this? I mean, I just to evaluate if the input is in the current/active MAPSET? grass.script.find_file() is a wrapper for g.findfile. Either check the returned mapset or use the mapset= argument. You can use grass.script.gisenv() to determine the current mapset. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Merging with Overlap
Hello, I've been using GRASS spatial merging tool (r.patch) with success but now I have one small problem: My tiles have a small overlap area. Is there any module to merge tiles but with methodologies to overlap these tiles (mean, maximum, etc)? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Stack Dumpfile in WinGRASS installation with OSGEO4w
Hello Mr. Neteler Thanks for the reply. Regarding that: 1- I don't have a virus/trohjan horse ( I runned my anti-virus yesterday) 2- I don't have that Web cam but I might have some software that is messing with Msys. Any idea how can I identify that? Thanks Monica On Tue, Jul 27, 2010 at 6:24 AM, Markus Neteler nete...@osgeo.org wrote: Monica, I searched in the internet for such messages and came across this http://readlist.com/lists/lists.sourceforge.net/mingw-msys/0/341.html Perhaps giving you some idea? Markus On Mon, Jul 26, 2010 at 6:26 PM, Monica Buescu monicabuescu1...@gmail.com wrote: Greetings A new update from myu precvious message (from which I didn't receive any answer): I'm installing weekly snapshot from last Friday and I'm gettingh a lot of Stackdumpfile printed in my mingw window and I'm getting a Lot of errors while building dlls which means that a lot of binaries are not being created. On Fri, Jul 23, 2010 at 11:53 AM, Monica Buescu monicabuescu1...@gmail.com wrote: Hi Since I saw today a message from Helmut, stating that 19-07 Snapshot is to be updated, I decided to install an earlier version (from June). While I was installing GRASS i got a lot of : 0 [main] sh 3880 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 960 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 2256 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 3724 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 772 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 2632 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 3464 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 2056 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump But I got a text stating that it was built without any errors. What is happening? Is this a problem? (because in my first attempt nothing was written in C:\osgeo4w\apps\grass and in second, without any change, it wrote there) Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Stack Dumpfile in WinGRASS installation with OSGEO4w
Thanks Helmut I will give a look on those links I have just tried to compile latest week snapshot and I got, besides my usual stackdumps, a lot of errors while creating libs and dbs: can this be associated with stackdumps? Because with earlier versions I ge tthe stackdumps but not errors in libs. Errors in: /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/datetime /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/gis /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/gmath /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/linkm /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/driver /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/pngdriver /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/psdriver /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/raster /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/bitmap /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/display /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/db/dbmi_base /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/db/dbmi_client /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/db/stubs /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/db/dbmi_driver /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/vask /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/edit /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/form /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/imagery /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/cluster /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/proj /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/segment /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/vector/diglib /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/vector/Vlib /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/vector/vedit /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/sites /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/rst/interp_float /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/g3d /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/gpde /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/symbol /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/init /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/stats /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/arraystats /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/ogsf /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/lib/nviz /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/db/drivers/dbf /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/db/drivers/odbc /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/db/drivers/postgres /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/db/drivers/sqlite /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/db/drivers/ogr /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/db/base /C/osgeo4w/usr/src/grass-6.4.svn_src_snapshot_2010_07_24/db/db.login On Tue, Jul 27, 2010 at 10:04 AM, Helmut Kudrnovsky hel...@web.de wrote: Hi, Hello Mr. Neteler Thanks for the reply. Regarding that: 1- I don't have a virus/trohjan horse ( I runned my anti-virus yesterday) 2- I don't have that Web cam but I might have some software that is messing with Msys. Any idea how can I identify that? Thanks Monica On Tue, Jul 27, 2010 at 6:24 AM, Markus Neteler [neteler at osgeo.org] wrote: Monica, I searched in the internet for such messages and came across this [http://readlist.com/lists/lists.sourceforge.net/mingw-msys/0/341.html] Perhaps giving you some idea? Markus http://old.nabble.com/stackdump-Problem-with-configure-or-msysinfo-td22180797.html http://betterlogic.com/roger/?p=1017 http://lists-archives.org/mingw-users/07238-i-broke-my-mingw.html http://lists-archives.org/mingw-users/12776-msys-sh-stackdump-during-make.html http://trac.osgeo.org/grass/wiki/CompileOnWindows#InstalltheenvironmentforcompilationMingW above links indicate that an update of the msys-files maybe could help. Helmut ___ Neu: WEB.DE De-Mail - Einfach wie E-Mail, sicher wie ein Brief! Jetzt De-Mail-Adresse reservieren: https://produkte.web.de/go/demail02 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Stack Dumpfile in WinGRASS installation with OSGEO4w
Greetings A new update from myu precvious message (from which I didn't receive any answer): I'm installing weekly snapshot from last Friday and I'm gettingh a lot of Stackdumpfile printed in my mingw window and I'm getting a Lot of errors while building dlls which means that a lot of binaries are not being created. On Fri, Jul 23, 2010 at 11:53 AM, Monica Buescu monicabuescu1...@gmail.comwrote: Hi Since I saw today a message from Helmut, stating that 19-07 Snapshot is to be updated, I decided to install an earlier version (from June). While I was installing GRASS i got a lot of : 0 [main] sh 3880 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 960 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 2256 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 3724 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 772 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 2632 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 3464 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 2056 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump But I got a text stating that it was built without any errors. What is happening? Is this a problem? (because in my first attempt nothing was written in C:\osgeo4w\apps\grass and in second, without any change, it wrote there) Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Stack Dumpfile in WinGRASS installation with OSGEO4w
Hi Since I saw today a message from Helmut, stating that 19-07 Snapshot is to be updated, I decided to install an earlier version (from June). While I was installing GRASS i got a lot of : 0 [main] sh 3880 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 960 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 2256 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 3724 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 772 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 2632 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 3464 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 2056 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump But I got a text stating that it was built without any errors. What is happening? Is this a problem? (because in my first attempt nothing was written in C:\osgeo4w\apps\grass and in second, without any change, it wrote there) Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to convert multi-bands clusters in a vectorial
Hi there I have used an external tool to produce clusters from Landsat images. For each input bands, I obtain clusters with a desginated value (mean of the values in that cluster). So, if I use 4 bands, I obtained 4 outputs with clusters. in each image I obtain the same clusters but with different values (depends on pixel values). How did you do this? (which commands)? I have used an external tool based on Segmentation paper. (coded in C++). What I want to obtain is a single map (raster or vectorial) with all clusters. in GRASS, how can I merge this 4 bands in order to obtain a labeled product? it can be a single raster or vectorial but it has to be a single band/layer... Yes please :) Eheheeh here goes :) I have 4 bands and in each one pixel values are calculated based on band cluster mean. I mean, pixels included in a certain cluster, are all with the same value that correspond with average of pixels included in the same cluster. This happens in all 4 bands. What I want is to merge 4 bands in a single one in which I will have labeled clusters with an integer (1,2,3,4,5...) isntead of pixel mean. How can I do this merge? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] How to convert multi-bands clusters in a vectorial
Greetings (This may sound as a complicated question but it's not :) ) I have used an external tool to produce clusters from Landsat images. For each input bands, I obtain clusters with a desginated value (mean of the values in that cluster). So, if I use 4 bands, I obtained 4 outputs with clusters. in each image I obtain the same clusters but with different values (depends on pixel values). What I want to obtain is a single map (raster or vectorial) with all clusters. in GRASS, how can I merge this 4 bands in order to obtain a labeled product? it can be a single raster or vectorial but it has to be a single band/layer... Thanks (if needed, I can better explain this with an example) Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Merge images
Greetings I have 2 landsat scenes acquired one after the other and I wanted to merge them. How can I do that? Thanks Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] managing mapsets in Python Scripts
Greetings I'm doing to image processing and I have 2 questions regarding mapset management 1- How can I define an image destination mapset? I mean should I first define mapset (g.mapset) and run my functions or should I define in output the im...@target_mapset? 2- How can I retrieve the current mapset? What is the grass modules to perform this? Thanks for your help Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Creating a winGRASS distribution
Greetings My name is Monica and I'm a grad students of Environmental Systems analysis. To a final project I have been developing all my work on GRASS (6.4.0) (Linux) which includes not only using GRASS' modules but adding a few ones (simples). My idea was, at the end of the project, give an Linux Installer of GRASS with all my added modules and edited Menubars. But, yesterday I realized that Linux is no good for my department so I need to do exacly this but with winGRASS. I gave a look around and: - imposible doing with winGRASS binary since it does not support to add more files than the ones already included - OSGEO4W not possible because it's installing a GRASS version available at OSGEO server. I would have to substitute that version :) - Any suggestions of doing a stand-alone installer of GRASS plus a few modules of my own? Thanks congrats for the excellent free-open-source software Monica ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
