Re: [GRASS-user] opening .gri anf .bin files in grass gis.

On 2016-05-26 24:23, Markus Neteler wrote: > On Tue, May 24, 2016 at 3:27 PM, Rengifo Ortega wrote: >> Dear Grass Community., >> first of all a greeting to the community. I have a question and hope some of >> you can help me with it. >> >> I have got a .gri and a .bin file

Re: [GRASS-user] Reprojection of raster maps in various platforms: projected coordinates are different

On 2015-10-10 22:56, Markus Neteler wrote: ... is this the upper box? In EU LAEA (EPSG 3035), my SIN point becomes -674447.24192338|9862187.84799817 As a small side-remark, related to the use of EPSG:3035 rather than to the OP's question: EPSG:3035 is related datum ETRS89, which in

Re: [GRASS-user] Run grass console from a cronjob

On 2015-08-26 16:27, Andrea Peri wrote: Hi, I'm try to run a grass session from a cronjob The cron execut a bash script where I try to execute Do you know GRASS_BATCH_JOBs? They run fine, AFAICT. Hermann https://grasswiki.osgeo.org/wiki/GRASS_and_Shell#GRASS_Batch_jobs

Re: [GRASS-user] Getting an error using r.mask in a GRASS shell script

On 2015-01-26 23:57, Thomas Adams wrote: All: I have a GRASS bash shell script where all my GRASS commands (7.0.0beta3) execute just fine except for my call to set a MASK: r.mask --overwrite input=$basin_mask@$MAPSET I get this error... Traceback (most recent call last): File

Re: [GRASS-user] Compile GRASS GIS 7

On 2014-10-11 17:42, stephen sefick wrote: Hello All: I am using the instructions here: http://grasswiki.osgeo.org/wiki/Compile_and_Install To compile GRASS 7 latest source code with svn up for SL Linux 6.5. I have errors in what seem to be all directories. I will provide the transcript of

Re: [GRASS-user] GRASS7: problems with WMS

On 2014-09-19 16:02, Anna Petrášová wrote: I have no idea, the input probably refers to m.proj but m.proj has input parameter. Could you try to run only m.proj (take some example from manual), if it works? Either I am doing something terribly stupid, or m.proj is *really* broken in the

Re: [GRASS-user] GRASS7: problems with WMS

On 2014-09-19 19:29, Anna Petrášová wrote: yes, that's pretty funny. Do you have any local changes (svn diff)? I looked in the recent changes but there is nothing suggesting this error. No diff at all. And until yesterday, I did not even know that m.proj existed. Hermann

Re: [GRASS-user] GRASS7: problems with WMS

On 2014-09-19 19:44, Hermann Peifer wrote: On 2014-09-19 19:29, Anna Petrášová wrote: yes, that's pretty funny. Do you have any local changes (svn diff)? I looked in the recent changes but there is nothing suggesting this error. No diff at all. And until yesterday, I did not even know

Re: [GRASS-user] GRASS7: problems with WMS

On 2014-09-19 23:35, Vaclav Petras wrote: Afterwards, you can also try to recompile GRASS (make distclean ./configure ... make). I actually did that, without any further investigations. After make distclean, I also used: svn up, which brought me from r61824 to r62033. When comparing

[GRASS-user] GRASS7: problems with WMS

Hi All, I am trying to load a WMS using GRASS7, which I compiled from svn/release branch. Connecting to the WMServer, selecting and adding a layer works fine, but I do not get any image into the map display. The GUI's Map Console reports: Unable to determine region, m.proj failed The

Re: [GRASS-user] GRASS7: problems with WMS

mtr 574 Sep 18 19:10 17189.5 -rw-r--r-- 1 peifer mtr9125 Sep 18 19:10 17189.6 -rw-r--r-- 1 peifer mtr8903 Sep 18 19:10 17189.7 -rw-r--r-- 1 peifer mtr 574 Sep 18 19:10 17189.8 -rw-r--r-- 1 peifer mtr 574 Sep 18 19:10 17189.9 On 2014-09-18 8:41, Hermann Peifer wrote: Hi All, I

Re: [GRASS-user] GRASS7: problems with WMS

be possible to send me the WMS service URL you are trying to connect? Best Stepan -- Původní zpráva -- Od: Hermann Peifer pei...@gmx.eu Komu: GRASS user list grass-user@lists.osgeo.org Datum: 18. 9. 2014 19:24:44 Předmět: Re: [GRASS-user] GRASS7: problems with WMS Just to add

Re: [GRASS-user] Add-on compilation errors in GRASS 7

On 2014-09-09 8:57, Martin Landa wrote: it's related to the problem on the server where modules.xml is generated. I am working on solving it. Hi Martin, I assume that at some point, there will be a full-size modules.xml file for grass7, i.e. comparable to

Re: [GRASS-user] r.info output after r.reclass

On 2014-08-24 5:44, Hamish wrote: Hermann wrote: I am using [1] as rules for r.reclass to shift a map's value range of -100..100 to 0..200. fyi, see also r.recode, r.rescale, and r.mapcalc. I am aware. I don't want to create new map with 200 000 x 200 000 px (all of Europe in 20m

Re: [GRASS-user] Grass Installation on mac - not working

On 2013-02-07 9:16, Erika Zarate Torres wrote: Dear GRASS community, I am new to GRASS and have a Mac OS X Lion 10.7.5. I followed the instructions in http://www.kyngchaos.com/software/grass and installed the following: 1-) GDAL complete 1.9 2-) Free Type 3-) cairo 4-) PIL 5-) GRASS.app 6.4.2-5

Re: [GRASS-user] Grib file: strange error

On 03/07/2012 14:22, Gabriele N. wrote: Hi everyone. Using the software zygrib (Also seen here http://grass.osgeo.org/wiki/GRIB) I download the file .grb of the area of interest. When importing the file in GRASS I get many raster ... ok but each raster is stretched along the y It looks

Re: [GRASS-user] Error when I start GRASS

On 03/06/2012 19:27, Prof. Ricardo Zúñiga wrote: ValueError: unknown locale: UTF-8 Error in GUI startup. If necessary, please report this error to the GRASS developers. Switching to text mode now. What is the output of locale on your system? Hermann

Re: [GRASS-user] GRASS region corruption

On 26/05/2012 20:26, spiderplant0 wrote: Does anyone know where what file the region is set in? I.e. is there some file I can delete to reset it. Look for a file named WIND, in the given mapset. I was about to point you to [1], then I noted that this page tells you the name of the default

[GRASS-user] Re: r.grow: Manual Page Needs Additional Information

On 26/04/2012 10:05, Markus Neteler wrote: On Tue, Apr 24, 2012 at 6:40 PM, Rich Shepardrshep...@appl-ecosys.com wrote: I've just been made aware of the r.grow module as part of the solution to the map mis-match issue with r.stream.extract. The r.grow manual page is quite clear on how to use

[GRASS-user] Re: r.mapcalc: Syntax Error?

On 05/04/2012 14:27, Rich Shepard wrote: On Thu, 5 Apr 2012, Hermann Peifer wrote: GRASS 6.5.svn (Washington-Kinross):~/grassdata r.mapcalc weight = if(tan_curv_5.lm 0, -100 * tan_curv_5.lm, \ if(tan_curv_7.lm 0, -100 * tan_curv_7.lm, \ if(tan_curv_11.lm 0, -100 * tan_curv_11.lm, 0.01

[GRASS-user] Re: Converting a XYZ positions in WGS84 to UTM29N

On 27/03/2012 19:41, katrin eggert wrote: Greetings I have an excel file with XY and point ID and I want to: 1- Create an Shapefile with this points (they were obtaiend with a GPS so it's WGS84) 2- Reproject to UTM29N can anyone just give me a few tips of the steps I need to do? - Save the

[GRASS-user] Re: NetCDF rotated pole

On 19/03/2012 18:50, Eduardo Klein wrote: Hi, Is there any way to read a NetCDF file which has a rotated pole Mercator projection and have it inside GRASS with the North pole at North? I'm trying to use AVISO SSH products on GRASS 6.4.1 and gdal 1.8.0 Thanks! This is possible but works in

[GRASS-user] Re: r.region for intersection with vector map?

On 14/03/2012 19:02, Rich Shepard wrote: I assume the reason is that the region is a square. Indeed. What command could I pass to r.mapcalc that would clip the raster map to the basins themselves Define a MASK before running r.mapcalc

Re: [GRASS-user] Re: r.region for intersection with vector map?

On 12/03/2012 10:13, Moritz Lennert wrote: On 12/03/12 05:57, Hermann Peifer wrote: On 12/03/2012 00:16, Helmut Kudrnovsky wrote: but r.region defines the extent of the whole raster map. Why would one want to do this, compared to leaving the extent as detected during import? Use cases

[GRASS-user] Re: r.region for intersection with vector map?

On 11/03/2012 17:35, Rich Shepard wrote: If I run r.region map=dem_east vect=analytical_boundary Try: g.region vect=analytical_boundary the raster map's region is that of the vector map. But, I don't believe that this removes all cells beyond that region. Advice on how to eliminate the

[GRASS-user] Re: r.region for intersection with vector map?

On 12/03/2012 00:16, Helmut Kudrnovsky wrote: but r.region defines the extent of the whole raster map. Why would one want to do this, compared to leaving the extent as detected during import? Hermann ___ grass-user mailing list

[GRASS-user] Re: Speeding Up Reprojection of Large Raster Maps

On 10/03/2012 20:07, Daniel Victoria wrote: ...The problem is that, if the grass command changes the region settings, things might not work. You can start additional GRASS sessions in other mapsets, change the region there, then work in parallel on various chunks of the big map. Hermann

[GRASS-user] Re: Possible bug in r.sum and in r.stats?

On 09/03/2012 00:16, Glynn Clements wrote: I would expect it to be the other way around. Most raster modules respect the current region and mask, because it happens automatically. Modules which don't want to use the current region have to explicitly set the working region (typically based upon

[GRASS-user] Re: Possible bug in r.sum and in r.stats?

On 08/03/2012 10:43, Rainer M Krug wrote: As some commands do use the region and mask, while others don't, it might be a good idea to make this clear in the manual - I know, it is stated if they do, but not if they don't. An sub-section under the Description might be an option? This

[GRASS-user] Zonal statistics with 3 maps

Hi, I am wondering how to re-produce in GRASS 6.4.2 the following data processing step which I have gotten in some sort of ArcGIS notation: Population_exposed = zonalstatistics( layer1, “value”, layer2 * layer3, sum, # ) layer1 and layer2 are CELL, layer3 is FCELL r.statistics doesn't

Re: [GRASS-user] r.univar: total null and non-null cells: -1580967296

On 11/09/2011 23:26, Hamish wrote: try r48240 in trunk. I downloaded grass7 from trunk and its r.univar now reports for the same map I used earlier: total null and non-null cells: 271400 (rather than -1580967296) Thanks for the fix. Will it also find its way into version 6.4.2?

[GRASS-user] r.univar: total null and non-null cells: -1580967296

Hi, r.univar gives me negative cell counts, whereas r.stats -c produces correct results (I hope ;-): (...) 41 118441 42 193 43 135 44 8189 * 2704177984 I guess this is the same integer overflow issue reported earlier [1]. I am just a bit puzzled that after changing to a 64-bit Linux OS,

Re: [GRASS-user] Check for non-overlapping raster maps

On 12/08/2011 08:52, Glynn Clements wrote: Hermann Peifer wrote: I have a mapset with 44 raster maps. I want to check that they do NOT have overlapping data areas. My thought is to do something like: r.mapcalc checkmap=isnull(map1) + isnull(map2) + ... + isnull(map44) If the min value

[GRASS-user] Re: Comparing 2 large Corine Land Cover maps

On 24/07/2011 00:15, Paulo van Breugel wrote: Maybe r.cross can be of help. It creates a cross product of the category variables of a layer, i.e., a map with separate categories for each unique combination of the two CLC maps (the category labels gives you the combination of input maps for each

[GRASS-user] Re: Check for non-overlapping raster maps

On 13/08/2011 07:12, Hamish wrote: Hermann wrote: Where I would see (perhaps) a small deficit in the communication or documentation is a simple overview about which (raster) modules respect the current region and resolution, and which don't. The issue went over the list [1], but as far as I can

[GRASS-user] Check for non-overlapping raster maps

Hi, I have a mapset with 44 raster maps. I want to check that they do NOT have overlapping data areas. My thought is to do something like: r.mapcalc checkmap=isnull(map1) + isnull(map2) + ... + isnull(map44) If the min value of the resulting checkmap is 43, then I know that the maps are

[GRASS-user] Re: Using non-ASCII chartacters in r.reclass rules file

On 30/07/2011 23:12, Glynn Clements wrote: Hermann Peifer wrote: Question: It seems that it is not possible to use non-ascii characters. Is it true? Is there any way to add labels to the classes and use characters (such as á;à;ó; â) I don't know; what error (or other unexpected behaviour

[GRASS-user] Re: Using non-ASCII chartacters in r.reclass rules file

On 29/07/2011 20:26, Glynn Clements wrote: Luisa Peña wrote: I want to add labels to a Land use/cover map with integer numbers. I have produced a reclass file where I associate a number with a land Cover class name 1 = Water 2 = Forests ... and It works fine Question: It seems that it is not

Re: [GRASS-user] Comparing 2 large Corine Land Cover maps

, 24 Jul 2011 00:15:29 +0200 Von: Paulo van Breugel p.vanbreu...@gmail.com An: Hermann Peifer pei...@gmx.eu CC: grass-user grass-user@lists.osgeo.org Betreff: Re: [GRASS-user] Comparing 2 large Corine Land Cover maps Maybe r.cross can be of help. It creates a cross product of the category

Re: [GRASS-user] Comparing 2 large Corine Land Cover maps

Glynn Clements wrote: Alternatively, the r.cross step could be done with r.mapcalc as e.g.: r.mapcalc 'result = map1 * 100 + map2' Thanks for the hint. I should have thought about it myself. r.mapcalc doesn't do short-circuit evaluation, which is why long if-else chains are

[GRASS-user] Comparing 2 large Corine Land Cover maps

Hi, I have 2 Corine Land Cover (CLC) raster maps, from 2000 and 2006. Both maps have the same extent, the same resolution (100m) and use the same 44 CLC categories. The maps are pretty large: 67 000 cols x 58 000 rows = 3 886 000 000 cells. I'd like to categorise the changes between the 2

[GRASS-user] Re: Eliminating values based on a null

On 16/02/2011 18:12, Marcello Gorini wrote: Jenny: in my base map (called X) I have a few null values and, in Y map I want to eliminate pixels that match with null values in X. How can I do this? besides using mapcalc r.series -n input=X,Y output=Y_without_X_nulls method=average (or any

[GRASS-user] Re: 2 questions about usage of g.region

On 03/02/2011 07:57, Micha Silver wrote: katrin eggert wrote: 3- Just to confirm this, when an image is imported with r.in.gdal it is kept with its original resolution right? it only gets the active region spatial resolution when I process it right (e.g. mapcalc or somethingg) AFAIK, yes.

[GRASS-user] Re: batch conversion of ESRI to MapInfo

ogr2ogr -f MapInfo File mapinfo/$file.tab $file You could cut away the shp suffix like this: ogr2ogr -f MapInfo File mapinfo/${file%shp}tab $file I am not fully sure what real world co-ordinates means, but this additional switch might do what you expect: -t_srs epsg:4326 Hermann On

Re: [GRASS-user] Re: Import question

On 23/09/2010 03:47, Hamish wrote: The fundamental and fantastic thing to realize about r.reclass is that is doesn't actually write a new map, it creates a virtual map based on classes of the original. It acts as a seamless filter between you and the data. With a bit of imagination it can be an

[GRASS-user] Re: Import question

On 22/09/2010 12:19, Glynn Clements wrote: Hermann Peifer wrote: By the way: is there a convenient overview of raster modules that do NOT respect the current region? In general: r.in.* ignore the region, and import rasters cell-for-cell. r.resamp.* and r.proj ignore the region for reading

Re: [GRASS-user] Re: Import question

On 21/09/2010 08:25, Hamish wrote: to sum up: * r.in.gdal differs from most raster modules in that it imports the entire map and ignores the current region. By the way: is there a convenient overview of raster modules that do NOT respect the current region? Hermann

[GRASS-user] Re: Import question

On 20/09/2010 18:09, Adam Dershowitz, Ph.D., P.E. wrote: I have a large ortho image, that is in MrSID format (A full county at high res). The file is ~3.5 gig. It covers a much larger area then I need. So I set my region to the size that I do need, and then imported. The import took

[GRASS-user] Re: which file do I use? .adf files

I guess you want to run r.in.gdal, or something. So you could follow this advice: To open the coverage select the coverage directory, or an .adf file (such as hdr.adf) from within it. http://www.gdal.org/frmt_various.html#AIG Hermann On 17/09/2010 02:33, stephen sefick wrote: I know that

[GRASS-user] Re: Importing multiple files with r.in.xyz

On 11/08/2010 15:07, Hanlie Pretorius wrote: This works nicely, except that my input files have names like: 3B42.000201.0.6.nc.lieb.txt 3B42.000201.3.6.nc.lieb.txt and I want the output rasters to be named 3B42.000201.0.6 3B42.000201.3.6 In fact, it would be best if they could be named:

[GRASS-user] Re: Patch raster maps very slow - Is there a way to speed this up?

I would try to warp patch the data with GDAL command line tools (gdalwarp in this case), then import the result into GRASS. See also http://trac.osgeo.org/gdal/wiki/UserDocs/GdalWarp#WillincreasingRAMincreasethespeedofgdalwarp Hermann Original Message Subject: Patch

[GRASS-user] Re: Unable to open gridded data

Illegal line in header ncols: 180 My guess would be that the colon is the culprit. See also http://en.wikipedia.org/wiki/ESRI_grid Hermann Original Message Subject: Unable to open gridded data From: Muhammad Rahiz muhammad.ra...@ouce.ox.ac.uk To: Date: 11/12/2009

[GRASS-user] Re: Unable to open gridded data

Oxford University Centre for the Environment, University of Oxford South Parks Road, Oxford, OX1 3QY, United Kingdom Tel: +44 (0)1865-285194 Mobile: +44 (0)7854-625974 Email: muhammad.ra...@ouce.ox.ac.uk Hermann Peifer wrote: Illegal line in header ncols: 180 My guess would

[GRASS-user] Re: r.mapcalc changes color table

Glynn Clements wrote: Tim Michelsen wrote: I have a large georeferenced digital map. I successfully imported it into GRASS. Now I would like to cut only a small rectangle -- my investigation area -- out of this large map. I used the following approach: 1) import the TIFF 2) select the part

Re: [GRASS-user] 2 r.neighbors questions

should be 1) from http://grass.itc.it/grass64/manuals/html64_user/r.neighbors.html Thanks, Hermann Glynn Clements wrote: Hermann Peifer wrote: I started using r.neighbors and came across the following issues: I have a CELL type input raster and calculate weighted sums. The weights

[GRASS-user] 2 r.neighbors questions

Hi, I started using r.neighbors and came across the following issues: I have a CELL type input raster and calculate weighted sums. The weights are float, so the calculation result should be float, but r.neighbors output raster is again CELL type with truncated integers. How do I get to

[GRASS-user] Re: v.out.ogr kml does not add name field

Works fine for me with Grass 6.4 from SVN. Which version are you using? Hermann maning sambale wrote: Hi, I have a grass vector of provincial boundaries db.describe -c generalized_adminncols: 3 nrows: 415 Column 1: cat:INTEGER:11 Column 2: region:CHARACTER:15 Column 3:

[GRASS-user] Re: How to convert from point to line ?

sgw00...@nifty.com wrote: Dear all I am a beginner user of grass gis. So am I. But as nobody else seems to have a better idea... here is what I would do: a) reformat yourdata into GMT format (yourdata.gmt) with a small AWK script $ awk -f reformat_to_gmt.awk yourdata yourdata.gmt b)

[GRASS-user] r.statistics method=distribution: no percentages for negative basemap values

Hi, Is there a (good) reason why there are no (or better: wrong) percentages in case of negative basemap values? Here some sampe results: basemap covermap percentage -70 -7 inf -70 63 inf -50 -7 inf -50 45 inf As far as I

[GRASS-user] r.stats: negative cell counts and percentages

Hi, I have the following raster map: Type of Map: raster Number of Categories: 255 Data Type:CELL Rows: 58000 Columns: 67000 Total Cells: 388600 r.stats (-c and -p) tells me: 1 6082321 -0.15% 2 13691259 2 -3.35%

Re: [GRASS-user] r.stats: negative cell counts and percentages

/09, Hermann Peifer pei...@gmx.eu wrote: From: Hermann Peifer pei...@gmx.eu Subject: [GRASS-user] r.stats: negative cell counts and percentages To: grass-user@lists.osgeo.org Received: Thursday, 1 October, 2009, 10:54 PM Hi, I have the following raster map: Type of Map: raster

[GRASS-user] Re: r.stats: negative cell counts and percentages

Hamish wrote: Hermann wrote: svn up https://svn.osgeo.org/grass/grass/branches/releasebranch_6_4 grass64_release not that it matters here, but do you use 32 bit or 64 bit CPU+OS? 32 what is the result of r.univar ? total null and non-null cells: -408967296 total null cells:

[GRASS-user] How to crop away NODATA margins

Hi, I see that d.out.file has a switch: -m (Do NOT crop away margins). So I assume that by default, white margins are auto-cropped. What I am looking for is some sort of auto-crop feature which I could use with r.out.gdal, in order to get rid of abundant NODATA margins in raster layers.

Re: [GRASS-user] r.statistics base=name cover=name method=count ?

Markus wrote: On Fri, Sep 25, 2009 at 8:26 AM, Hermann Peifer pei...@gmx.eu wrote: Hi, I see that there is a method=distribution which gives me the percentage values, but is there an easy way to get to a count of cover layer values ? Maybe r.clump and then r.report with counts? Thanks

Re: [GRASS-user] How to crop away NODATA margins

g.region zoom=myrast align=myrast I overlooked this in the documentation, I have to admit :-( Sorry for the noise and thanks to both of you. Hermann ___ grass-user mailing list grass-user@lists.osgeo.org

Re: [GRASS-user] r.statistics base=name cover=name method=count ?

Markus wrote: Maybe r.univar.zonal from grass-addons can help? It could give you the number of cells of the cover layer (count of cover layer values) per base layer (zoning map in r.univar.zonal) class. Not sure if this is what you are looking for. This is actually what I am looking for,

[GRASS-user] Re: Strange r.mapcalc results

Glynn Clements wrote: Hermann Peifer wrote: Thanks indeed for the enlightening. I am new to the mapcalc business and my obviously wrong (AWK-based) thinking was that mapcalc would do a lazy evaluation of the conditions, i.e. 1 1 || something would always be 1, whatever something was. Now I

[GRASS-user] Re: Strange r.mapcalc results

r.univar based comparison: What if there was a small geometric shift between the 2 maps. Would I notice if both maps have enough NULLs around the edges? r.mapcalc diffmap = map1 - map2 How would I separate NULL,NULL cases from cases where one map has a value, but the other has NULL?

[GRASS-user] Re: Strange r.mapcalc results

Glynn Clements wrote: Hermann Peifer wrote: I am trying to find out if 2 raster layers are exactly identical. I am not an experienced mapcalc user, so I am wondering, where the 3200204 NULL values come from, in Test1. The result in Test2 is what I expected. My obviously wrong understanding

[GRASS-user] Strange r.mapcalc results

Hi, I am trying to find out if 2 raster layers are exactly identical. I am not an experienced mapcalc user, so I am wondering, where the 3200204 NULL values come from, in Test1. The result in Test2 is what I expected. My obviously wrong understanding was that I wouldn't need the special

[GRASS-user] Re: Exporting attributes to GMT ASCII format using v.out.ogr

Hamish wrote: Eric: I'm using v.out.ogr to export GRASS vectors to GMT format. The geometry exports fine, but the attributes don't seem to be in the GMT file after the export finishes. I can see the attribute column names in the GMT header portion of the file, but all I get in the data columns