Re: [GRASS-user] r.stream.order export to shapefile

2012-08-28 Thread Markus Metz
On Tue, Aug 28, 2012 at 9:05 PM, Micha Silver mi...@arava.co.il wrote: A few more details regarding v.out.ogr when layer=2 I see this ticket, possibly similar: http://trac.osgeo.org/grass/ticket/991 I tried to export to PostGIS, and again, all attribute fields are created but all values are

Re: [GRASS-user] Precision of raster/vector float values

2012-08-26 Thread Markus Metz
On Sun, Aug 26, 2012 at 12:19 AM, Glynn Clements gl...@gclements.plus.com wrote: Markus Metz wrote: So I misunderstood something about IEEE precision, sorry about the confusion! Coming back to the original problem, it seems that it can be solved for DCELL by using .15g instead of .10f

Re: [GRASS-user] Problem with Shape import - v.in.ogr crashed

2012-08-26 Thread Markus Metz
On Sun, Aug 26, 2012 at 10:23 PM, Martin Landa landa.mar...@gmail.com wrote: Hi, 2012/8/26 Valentin Wittich vwitt...@gmail.com: Hi all, I need help with a simple function. I tried to import a Shape file with v.in.ogr and as result the function crashed / aborted (core dumped). If possible,

Re: [GRASS-user] Precision of raster/vector float values

2012-08-24 Thread Markus Metz
On Thu, Aug 23, 2012 at 11:52 PM, Glynn Clements gl...@gclements.plus.com wrote: Markus Metz wrote: So if I understand you correctly: 1.0e-15 is the smallest number for values in attribute columns although the values in a DCELL raster can be smaller? A bottleneck is v.what.rast which

Re: [GRASS-user] Precision of raster/vector float values

2012-08-24 Thread Markus Metz
On Fri, Aug 24, 2012 at 10:40 AM, Johannes Radinger johannesradin...@gmail.com wrote: Hi, So I misunderstood something about IEEE precision, sorry about the confusion! Coming back to the original problem, it seems that it can be solved for DCELL by using .15g instead of .10f, i.e. use the

Re: [GRASS-user] Precision of raster/vector float values

2012-08-24 Thread Markus Metz
On Fri, Aug 24, 2012 at 11:27 AM, Markus Metz markus.metz.gisw...@gmail.com wrote: On Fri, Aug 24, 2012 at 10:40 AM, Johannes Radinger johannesradin...@gmail.com wrote: Hi, So I misunderstood something about IEEE precision, sorry about the confusion! Coming back to the original problem

Re: [GRASS-user] Precision of raster/vector float values

2012-08-23 Thread Markus Metz
On Thu, Aug 23, 2012 at 10:35 AM, Johannes Radinger johannesradin...@gmail.com wrote: Hi, I have a GRASS raster map which contains small values (e.g 2.5e-15) among others. This is (dangerously) close to the IEEE limit of double precision floating point values. If GRASS tests for fp precision,

Re: [GRASS-user] Precision of raster/vector float values

2012-08-23 Thread Markus Metz
On Thu, Aug 23, 2012 at 2:40 PM, Johannes Radinger johannesradin...@gmail.com wrote: I have a GRASS raster map which contains small values (e.g 2.5e-15) among others. This is (dangerously) close to the IEEE limit of double precision floating point values. If GRASS tests for fp precision, it

Re: [GRASS-user] Precision of raster/vector float values

2012-08-23 Thread Markus Metz
Radinger johannesradin...@gmail.com wrote: On Thu, Aug 23, 2012 at 3:23 PM, Markus Metz markus.metz.gisw...@gmail.com wrote: On Thu, Aug 23, 2012 at 2:40 PM, Johannes Radinger johannesradin...@gmail.com wrote: I have a GRASS raster map which contains small values (e.g 2.5e-15) among others

Re: [GRASS-user] update cat column in vector data set

2012-08-16 Thread Markus Metz
On Thu, Aug 16, 2012 at 1:10 PM, Stefan Luedtke slued...@gfz-potsdam.de wrote: Dear friends, I am struggling with the following point. I have created a shape of subbasins based on a raster data set. Since this yields some tiny (pixel size) areas, I run v.clean with a certain threshold to

Re: [GRASS-user] [GRASS-dev] please help review osgeo live dvd summary docs

2012-07-25 Thread Markus Metz
On Wed, Jul 25, 2012 at 5:59 AM, Hamish hamis...@yahoo.com wrote: Markus Metz wrote: An example for vector processing using the basins raster created in the previous step: 1) convert the basins to vector areas with r.to.vect -v feature=area 2) load average elevation for each basin

Re: [GRASS-user] please help review osgeo live dvd summary docs

2012-07-24 Thread Markus Metz
On Mon, Jul 23, 2012 at 2:32 AM, Hamish hamis...@yahoo.com wrote: There's a copy of the quickstart in the GRASS wiki. Note that when I wrote it there was a GRASS workshop at the targeted FOSS4G conference which was going to focus on vector processing, so I didn't bother to put much of that in

Re: [GRASS-user] sqlite-db broken?

2012-07-17 Thread Markus Metz
On Mon, Jul 16, 2012 at 2:14 PM, Johannes Radinger jradin...@gmx.at wrote: Hi again, I did another test now and found out following: It seems that the problem with the sqlite database might be related to the directory of my grass location. When I try my to save the location on my local

Re: [GRASS-user] sqlite-db broken?

2012-07-17 Thread Markus Metz
Hi, On Mon, Jul 16, 2012 at 4:02 PM, Johannes Radinger jradin...@gmx.at wrote: Hi, after some search in the internet, it seems the sqlite doesn't work very well with samba-network drives (mounted via cifs). There is somehow an issue with the locking of the db. I think I found the reason

Re: [GRASS-user] dpi output from ps.map

2012-07-14 Thread Markus Metz
On Sat, Jul 14, 2012 at 3:14 AM, Gerald Nelson nelson.ger...@gmail.com wrote: We need to have a raster eps or pdf that has a dpi of 300. Reading through the old list entries it appears that ps.map had a hard coded limit of 75 dpi. The spatial detail represented in the output of ps.map is

Re: [GRASS-user] Invalid region and coordinates when trying to reproject.

2012-07-02 Thread Markus Metz
On Mon, Jul 2, 2012 at 9:03 AM, Marcello Gorini gor...@gmail.com wrote: Makus: What are the extents of the vector in latlon before projecting into the mercator location? Note that regions spanning 180 degrees longitude or more can not be reprojected to mercator locations due to mathematical

Re: [GRASS-user] How to implement this kind of operation?

2012-07-02 Thread Markus Metz
On Mon, Jul 2, 2012 at 10:56 AM, DavidRA theboss...@gmail.com wrote: I guess I could use the code of r.series, modifying it a bit to obtain a new function. The problem is that... man, that code is hard to understand. I do have some notions of C, but here I'm finding a lot of trouble to see what

Re: [GRASS-user] calculate mode value on moving window (with mapalgebra) with null value cells

2012-06-20 Thread Markus Metz
, if the border problem arouse from reprojecting with a method other than nearest, you can use in trunk as resampling method one of bilinear_f,cubic_f,lanczos_f with r.proj to avoid the border effect. Markus M giovanni 2012/6/19 Markus Metz markus.metz.gisw...@googlemail.com On Tue, Jun 19

Re: [GRASS-user] calculate mode value on moving window (with mapalgebra) with null value cells

2012-06-19 Thread Markus Metz
On Tue, Jun 19, 2012 at 12:39 PM, G. Allegri gioha...@gmail.com wrote: Hi Marcello, thanks, but my problem is not finding the border. I already have it (with its own category) but assigning it the value from the surrounding (not null) value, e.g. max/min/etc. The problem with kernel

Re: [GRASS-user] calculate mode value on moving window (with mapalgebra) with null value cells

2012-06-19 Thread Markus Metz
On Tue, Jun 19, 2012 at 6:15 PM, G. Allegri gioha...@gmail.com wrote: AFAICT, r.neighbors ignores NULL values and assigns the new value from the surrounding non-NULL values. You could then patch your original map with the output of r.neighbors. I have considered r.neighbors but I need to

Re: [GRASS-user] Feet-to-Meter Conversion

2012-06-18 Thread Markus Metz
On Mon, Jun 18, 2012 at 5:46 PM, Rich Shepard rshep...@appl-ecosys.com wrote:  The standard conversion factor from feet to meters is 1 ft = 0.3048 m. However, the conversion factor for US Survey Feet (frequently the distance unit in the US State Plane Coordinate System) is 1 US ft = 0.3600 m.

Re: [GRASS-user] r.viewshed add-on is taking too long

2012-06-13 Thread Markus Metz
On Wed, Jun 13, 2012 at 12:21 AM, Chaulio Ferreira chau...@hotmail.com wrote: Markus M, Thanks for answering. I'll try to use the default 500mb. But I still think these differences in the times are too big. I was testing a bit and can confirm the differences. Still, r.viewshed is much faster

Re: [GRASS-user] 回覆: Re: ridge extraction from DEM

2012-06-12 Thread Markus Metz
On Tue, Jun 12, 2012 at 11:28 AM, Margherita Di Leo direg...@gmail.com wrote: Hi, there is also r.convergence add-on: http://grass.osgeo.org/wiki/GRASS_AddOns#r.convergence BTW, the topographic convergence index is also available in r.watershed in GRASS 7. Markus M Regards, madi On

Re: [GRASS-user] Import grd/gri-Raster from R-package 'dismo'

2012-06-12 Thread Markus Metz
On Tue, Jun 12, 2012 at 11:41 AM, Johannes Radinger jradin...@gmx.at wrote: Hi, I doing species distribution modeling in R using the package 'dismo'[1]. The predict() function produces a output raster file with the prediction of the species distribution. The file which is created is a grd

Re: [GRASS-user] NVIZ can't allocate enough memory

2012-06-12 Thread Markus Metz
On Tue, Jun 12, 2012 at 5:23 AM, Glynn Clements gl...@gclements.plus.com wrote: Daniel Lee wrote: ERROR: G_malloc: unable to allocate 18446744071974792324 bytes at gsds.c:575 18446744071974792324 = 0x9899ac84 Does anybody know what could be the problem? If I'm interpreting the

Re: [GRASS-user] r.viewshed add-on is taking too long

2012-06-12 Thread Markus Metz
On Tue, Jun 12, 2012 at 9:52 PM, Chaulio Ferreira chau...@hotmail.com wrote: Hello everyone, I've been trying to compute some viewsheds using the r.viewshed module (available in http://grass.osgeo.org/wiki/GRASS_AddOns), but it's taking to long to finish. I ran this module on a 5000x5000

Re: [GRASS-user] cut off floating map

2012-06-10 Thread Markus Metz
On Sat, Jun 9, 2012 at 4:50 PM, Salvatore Mellino salvatore.mell...@gmail.com wrote: Hi, is it possible to cut off a floating point map to the 4th decimal digit? Which is the procedure? Try r.mapcalc rounded = int(fp_map * 1) / 1.0 Markus M

Re: [GRASS-user] Why r.watershed Yields Narrow Basins

2012-06-10 Thread Markus Metz
On Sat, Jun 9, 2012 at 6:38 PM, Rich Shepard rshep...@appl-ecosys.com wrote: On Sat, 9 Jun 2012, Margherita Di Leo wrote: Try also multiple flow direction (MFD) algorithm ( -f  flag in r.watershed) madi,  That was much worse. The basins were even less realistic. I also tried the -4

Re: [GRASS-user] patch boundaries

2012-06-08 Thread Markus Metz
On Fri, Jun 8, 2012 at 5:58 AM, Daniel Victoria daniel.victo...@gmail.com wrote: Hi all, I'm having a hard time trying to get one boundary polygon file from four smaller shapes. I'm using wingrass7, daily build from today. This is how it goes. I have four boundaries which I imported and

Re: [GRASS-user] Import R point dataframe as points in GRASS

2012-06-08 Thread Markus Metz
On Fri, Jun 8, 2012 at 2:44 PM, Johannes Radinger jradin...@gmx.at wrote: Hi, I have a R data.frame containing around 60 columns, two of them are X and Y respectively. What is the recommended way to get these points imported into GRASS? Maybe the simplest way is to save the data.frame as

Re: [GRASS-user] Attribute transfer from one table to another based on position of vector items

2012-06-08 Thread Markus Metz
On Fri, Jun 8, 2012 at 2:25 PM, Jaromír Kamler kam...@centrum.cz wrote: Hi, I have two vectors imported from SHP. One contains lines (areas), which represents sites and in other are texts with cost of sites. I imported them to GRASS. Texts with costs of ground are now points with cost in

Re: [GRASS-user] Report from ongoing GRASS GIS Community Sprint in Prague

2012-05-30 Thread Markus Metz
cavall...@faunalia.it wrote: Il 28/05/2012 21:17, Markus Metz ha scritto: not to java SEXTANTE. While on it (the blaming), I think it is time for Paolo to apologize for his rude [1] comments [2][3][4][5][6][7] earlier in the thread. I'm more than happy to apologize if I have hurt

Re: [GRASS-user] v.dissolve renames the attribute column

2012-05-29 Thread Markus Metz
Moritz Lennert wrote: On 28/05/12 12:30, Markus Metz wrote: On Mon, May 28, 2012 at 11:52 AM, spiderplant0spiderpla...@gmail.com  wrote: Hi. It doesn't add any new information but I wanted the .dbf file to be there so I can open it up with other tools and check the contents for debug

Re: [GRASS-user] v.dissolve renames the attribute column

2012-05-28 Thread Markus Metz
On Mon, May 28, 2012 at 11:52 AM, spiderplant0 spiderpla...@gmail.com wrote: Hi. It doesn't add any new information but I wanted the .dbf file to be there so I can open it up with other tools and check the contents for debug purposes etc. And I wanted to change the column to the correct name

Re: [GRASS-user] Report from ongoing GRASS GIS Community Sprint in Prague

2012-05-28 Thread Markus Metz
On Mon, May 28, 2012 at 2:47 PM, Paolo Cavallini cavall...@faunalia.it wrote: Il 28/05/2012 14:33, Benjamin Ducke ha scritto: SEXTANTE and gvSIG CE share a bug tracker at http://bugs.gvsigce.org, a decision that was made together with Victor Olaya a while Hi all. Sorry for the

Re: [GRASS-user] Report from ongoing GRASS GIS Community Sprint in Prague - i.atcoor

2012-05-28 Thread Markus Metz
Daniel Victoria wrote: Crap... I looked out the windows and realized the atcoor vs. flaash comparison I just did had a major flaw. I also have a good linear regression between uncorrected radiance and reflectance from bothe flaash and atcoor. So a linear regression between flaash and atcoor

Re: [GRASS-user] Report from ongoing GRASS GIS Community Sprint in Prague

2012-05-28 Thread Markus Metz
On Mon, May 28, 2012 at 6:02 PM, Paolo Cavallini cavall...@faunalia.it wrote: Il 28/05/2012 17:56, José Antonio Canalejo Alonso ha scritto: we have put also a lot effort in SEXTANTE and in gvSIG/gvSIG CE. I like a lot these projects and also QGIS and I think it's positiv to see both  

Re: [GRASS-user] v.reclass and sqlite

2012-05-27 Thread Markus Metz
On Sun, May 27, 2012 at 5:01 PM, Micha Silver mi...@arava.co.il wrote: Looking at v.dissolve ( which uses v.reclass ) I see that it fails when the vector has an sqlite based attribute table: This sqlite issue should already be fixed in GRASS 7, and v.dissolve waits for testing before the fix

Re: [GRASS-user] v.dissolve renames the attribute column

2012-05-27 Thread Markus Metz
On Sun, May 27, 2012 at 2:38 PM, spiderplant0 spiderpla...@gmail.com wrote: When I run v.dissolve on my vector... *v.dissolve --o input=pass2@mapset1 layer=1 column=fttcp output=pass3* The dissolving appears to work but the column name in the attribute table has changed from 'fttcp' to

Re: [GRASS-user] XY trace plots

2012-05-20 Thread Markus Metz
Hamish wrote: [...] I've just added an experimental module in GRASS addons called r.to.vect.lines.   http://grass.osgeo.org/wiki/GRASS_AddOns#r.to.vect.lines ps- if any python+ctypes experts want to have a look over it please do... especially to fix the G_is_null_value() call which is

[GRASS-user] Re: Help - Output difference between ARCGIS enviroment settings and GRASS g.region function

2012-05-09 Thread Markus Metz
On Tue, May 8, 2012 at 9:41 PM, Alejandro Coca Castro aco...@gmail.com wrote: Hi Mark, i followed your instructions, so i put in command-line This g.region rast=MASK@default nsres=0:00:07.5 ewres=0:00:07.5 will adjust, i.e. modify the resolution in order to fit into the extents. The extents

Re: [GRASS-user] r.watershed MFD strange basin

2012-05-09 Thread Markus Metz
On Wed, May 9, 2012 at 9:19 AM, ivan marchesini ivan.marches...@gmail.com wrote: Dear Markus, thank you for your answer. Please report if the stream segments are as expected, because the determination of half-basins and basins is based on the stream segments. this is the map of streams:

Re: [GRASS-user] Re: Help - Output difference between ARCGIS enviroment settings and GRASS g.region function

2012-05-09 Thread Markus Metz
of Arcmap ASCII... I think instead of r.out.arc, you can also use r.out.gdal format=AAIGrid and see if that gives a better match. Markus M Mark, i set the region as you tried... Thanks in advance, 2012/5/9 Markus Metz markus.metz.gisw...@googlemail.com On Tue, May 8, 2012 at 9:41 PM

Re: [GRASS-user] r.watershed MFD strange basin

2012-05-08 Thread Markus Metz
On Mon, May 7, 2012 at 1:08 PM, ivan marchesini ivan.marches...@gmail.com wrote: Dear Grass users, I have noted a strange result from r.watershed using the MDF option (with all levels of convergence). This is a part of my DEM: http://imageshack.us/photo/my-images/690/demv.png/ This is

Re: [GRASS-user] Help - Output difference between ARCGIS enviroment settings and GRASS g.region function

2012-05-08 Thread Markus Metz
On Tue, May 8, 2012 at 6:09 PM, Alejandro Coca Castro aco...@gmail.com wrote: Hi, i tried to convert the following enviroments settings to GRASS g.region and r.mask for the next ARCGIS output (ASCII file):     fileTIF = rx:\test\fileTIF.tif     MASK = rx:\test\mask_file     env.extent =

Re: [GRASS-user] Re: Help - Output difference between ARCGIS enviroment settings and GRASS g.region function

2012-05-08 Thread Markus Metz
On Tue, May 8, 2012 at 6:27 PM, Alejandro Coca Castro aco...@gmail.com wrote: Hi, i forgot... below I have included the gdalinfo from both ouputs (ArcMap and GRASS). Any ideas as to what else I should check would be appreciated. ArcMap output gdalinfo The extents and resolution seem

Re: [GRASS-user] Relationships of Raster Vector Attribute Data

2012-05-04 Thread Markus Metz
On Thu, May 3, 2012 at 12:49 AM, Rich Shepard rshep...@appl-ecosys.com wrote: On Thu, 3 May 2012, Markus Neteler wrote: It is in the cats/ subdirectory. Markus,  Interesting: # 166 categories Watershed basins 0.00 0.00 0.00 0.00 dem.basin lines 1-4/4 (END) They are created if the

Re: [GRASS-user] Re: Pixels values ​​are changed when the layer is exported to vector (column value in attribute tables)

2012-05-04 Thread Markus Metz
On Thu, May 3, 2012 at 10:07 PM, Marcello Benigno benigno.marce...@gmail.com wrote: Any solution? I also tried using the r.thin but keep having the same problem If you are using GRASS 6.4.1, this is a bug in v.out.ogr that has been fixed in 6.4.2. Markus M Thanks in advance, 2012/5/2

Re: [GRASS-user] Problem with v.out.ogr?

2012-04-30 Thread Markus Metz
If this happens with GRASS 6.x, I think the expected output can be obtained by adding the -c flag: Export features with category (labeled) only. Otherwise all features are exported v.out.ogr type=line input=esw_g26a format=KML dsn=main.kml layer=7 -c The meaning of the -c flag has been reverted

Re: [GRASS-user] internal buffer of polygon

2012-04-29 Thread Markus Metz
On Thu, Apr 19, 2012 at 6:17 PM, Markus Neteler nete...@osgeo.org wrote: Hi all, does anyone have a new(er) trick for negative vector buffers (in GRASS 7)? Internal buffers of areas are now (r51573) supported in trunk with negative distances if GEOS support is enabled, i.e. v.buffer

Re: [GRASS-user] WinGrass 7

2012-04-24 Thread Markus Metz
On Tue, Apr 24, 2012 at 4:13 PM, Daniel Victoria daniel.victo...@gmail.com wrote: Hi all, Just installed the latest Wingrass 7, build 51507-212, from http://wingrass.fsv.cvut.cz/grass70/ and when I try to run Grass 7 I get the following error at the windows comand prompt:

Re: [GRASS-user] r.stream.extract: map mis-match error

2012-04-23 Thread Markus Metz
On Mon, Apr 23, 2012 at 12:17 AM, Rich Shepard rshep...@appl-ecosys.com wrote: On Sun, 22 Apr 2012, Markus Metz wrote: I guess so because you said that r.stream.extract works with the unweighed accumulation map. How exactly did you calculate the weight? Markus, r.param.scale in=dem_10m

Re: [GRASS-user] r.stream.extract: map mis-match error

2012-04-22 Thread Markus Metz
On Fri, Apr 20, 2012 at 7:01 PM, Rich Shepard rshep...@appl-ecosys.com wrote: On Fri, 20 Apr 2012, Markus Metz wrote: Try to set the region and MASK before the analysis, then don't modify these settings during the analysis. Something like g.region vect=basins align=dem_proj g.region res=10m

Re: [GRASS-user] r.stream.extract: map mis-match error

2012-04-20 Thread Markus Metz
On Thu, Apr 19, 2012 at 10:57 PM, Rich Shepard rshep...@appl-ecosys.com wrote:  On one project the r.stream.extract error, Accumulation map does not match elevation map!, was caused by a mask on the elevation map. With this project the mask was removed, and no longer limits the DEM region. It

Re: [GRASS-user] Vector not getting imported in GRASS v7.0 latest svn.

2012-04-18 Thread Markus Metz
On Wed, Apr 18, 2012 at 10:26 AM, Pankaj Kr Sharma pks...@gmail.com wrote: Vector not getting imported in GRASS v7.0 latest svn. After completion of the import process, following error message is displayed : Command 'd.vect map=%@%' failed Details: The bounding box of the map is outside the

Re: [GRASS-user] v.db.* commands in GRASS-7

2012-04-12 Thread Markus Metz
Andres Kuusk wrote: Hi! It seems to me that some v.db.* modules do not do what is described in the manual. I am using GRASS-7.0 r51370 and Postgres data base. In the manual I read: It is possible to link the geographic objects in a vector map to one or more tables. All vector objects

[GRASS-user] Re: [GRASS-dev] how to fix coor file size error

2012-04-11 Thread Markus Metz
On Wed, Apr 11, 2012 at 12:19 AM, Michael Barton michael.bar...@asu.edu wrote: In GRASS 7, how can you fix the following error? Coor file of vector map coast3@SAA2012_survey_analysis is larger than it should be (154662962 bytes excess) That can happen if a module exists with an error before

Re: [GRASS-user] r.stream.extract Syntax Problem

2012-04-10 Thread Markus Metz
[Answering to your latest post in a new thread, but I prefer to keep it all in one thread.] On Tue, Apr 10, 2012 at 1:54 AM, Rich Shepard rshep...@appl-ecosys.com wrote: On Mon, 9 Apr 2012, Markus Metz wrote: The computational region should be aligned to the DEM in order to get accurate

Re: [GRASS-user] v.to.rast in grass70

2012-04-10 Thread Markus Metz
On Mon, Apr 9, 2012 at 12:08 PM, Paulo van Breugel p.vanbreu...@gmail.com wrote: I can confirm there are issues in how GRASS 7 handles region settings. In my case, when using v.to.rast, the region settings are reset to the default, i.e., the raster is converted using the default region settings

Re: [GRASS-user] Re: r.stream.extract Syntax Problem

2012-04-09 Thread Markus Metz
Considering the length of this thread and the various issues mentioned, it seems that I need to update the manual of r.stream.extract a bit. I will try to address the issues: r.stream.extract is designed to work with a DEM only or with a DEM and a flow accumulation map as input. This flow

Re: [GRASS-user] trouble with r.terraflow

2012-04-09 Thread Markus Metz
tomte wrote: Dear folks, I am using GRASS 6.4.1 (2011) on Mac (mac os x 10.7.3). Generally I am running GRASS without any trouble but calling r.terraflow from the command line generates the following error comments: MM error: limit =0B. allocating 16B. limit exceeded by 24B. Assertion

Re: [GRASS-user] Re: trouble with r.terraflow

2012-04-09 Thread Markus Metz
On Mon, Apr 9, 2012 at 9:09 PM, tomte tomte.tume...@gmx.ch wrote: below the output of g.region -p GRASS 6.4.1 (kana_cat0206):~/GIS/kanada/modelling/nka_05m/catch_0206_SA g.region -p projection: 1 (UTM) zone:       11 datum:      nad83 ellipsoid:  grs80 north:      5652250 south:      

Re: [GRASS-user] Unable to get rid of duplicate polygons

2012-04-04 Thread Markus Metz
are not messing up your polygon layer. On 04/02/2012 09:09 PM, David J. Bakeman wrote: Markus Metz wrote: On Sun, Apr 1, 2012 at 11:45 PM, David J. Bakemandbake...@comcast.net wrote: David J. Bakeman wrote: Markus Neteler wrote: On Sun, Apr 1, 2012 at 9:45 PM, David J. Bakemandbake

Re: [GRASS-user] Unable to get rid of duplicate polygons

2012-04-02 Thread Markus Metz
On Sun, Apr 1, 2012 at 11:45 PM, David J. Bakeman dbake...@comcast.net wrote: David J. Bakeman wrote: Markus Neteler wrote: On Sun, Apr 1, 2012 at 9:45 PM, David J. Bakemandbake...@comcast.net  wrote: I didn't find an answer in the archives so. I have a shapefile of polygons and some of

Re: [GRASS-user] Import a catalog of ascii *.arc files and merge to a raster

2012-04-01 Thread Markus Metz
Lars Forseth wrote: Thanks Daniel. But no, I was not able to make your suggestion work; it complains that it dos not know the command 'For'... In Linux shell and msys this must be lowercase, i.e. for i in *.arc do r.in.arc $i output name here done or for i in *.arc do r.in.arc in=$i

Re: [GRASS-user] Vector map too big for v.proj

2012-03-28 Thread Markus Metz
On Tue, Mar 27, 2012 at 8:21 PM, Kurt Springs ferret_b...@mac.com wrote: I am having a problem reprojecting vector maps from lat/long to Irish National Grid.  I have been using v.proj: v.proj input=Lakes location=IrelandGeo output=Lakes1 and got: pj_transform() failed: latitude or

Re: [GRASS-user] Vegetation index with OpenMP

2012-03-28 Thread Markus Metz
On Tue, Mar 27, 2012 at 1:17 PM, Ripsime Khachatryan kul...@mail.ru wrote: Hi all I would like to asks if anybody calculated the vegetation index using OpenMP, and can anybody give me detailed information how to do that. I performed the same operation using MPI, so I am familiar with

Re: [GRASS-user] Re: Grid Mismatch while importing with 6.4.1 and 6.4.2

2012-03-28 Thread Markus Metz
On Wed, Mar 28, 2012 at 2:06 PM, antonio.ro...@deimos.com.pt wrote: Greetings Daniel lee suggested me that the problem might be the defined region. But, as far as I know, Defined region has no impact over importing data. This is true right? Yes, r.in.gdal ignores the current region. The

Re: [GRASS-user] Re: Grid Mismatch while importing with 6.4.1 and 6.4.2

2012-03-28 Thread Markus Metz
mentioned only a grid mismatch between the two imported versions. Comparing r.info of these to gdalinfo of the original raster will tell you which version is correct. Markus M On Wed, 28 Mar 2012 14:49:42 +0200, Markus Metz markus.metz.gisw...@googlemail.com wrote: On Wed, Mar 28, 2012 at 2:06 PM

Re: [GRASS-user] Re: problem using r.watershed

2012-03-27 Thread Markus Metz
On Mon, Mar 26, 2012 at 5:32 PM, Filipe Silva Dias filipesd...@gmail.com wrote: Thanks, it took a while but it worked. You are welcome to compare it to alternatives ;) On Mon, Mar 26, 2012 at 3:11 PM, Markus Metz markus.metz.gisw...@googlemail.com wrote: On Sun, Mar 25, 2012 at 2:59 PM

Re: [GRASS-user] Re: problem using r.watershed

2012-03-26 Thread Markus Metz
On Sun, Mar 25, 2012 at 11:49 AM, Markus Metz markus.metz.gisw...@googlemail.com wrote: On Sun, Mar 25, 2012 at 12:30 PM, Filipe Silva Dias filipesd...@gmail.com wrote: I switched the resolution from 30x30 to 100x100 and now it works. Problem solved. On Sun, Mar 25, 2012 at 11

Re: [GRASS-user] Re: problem using r.watershed

2012-03-25 Thread Markus Metz
On Sun, Mar 25, 2012 at 12:30 PM, Filipe Silva Dias filipesd...@gmail.com wrote: I switched the resolution from 30x30 to 100x100 and now it works. Problem solved. On Sun, Mar 25, 2012 at 11:12 AM, Filipe Silva Dias filipesd...@gmail.com wrote: I'm using GRASS 6.4.2 and I'm on Windows 7

Re: [GRASS-user] Re: problem using r.watershed

2012-03-25 Thread Markus Metz
Sensing und dem GIS-Lab Marburg. Am 25. März 2012 12:49 schrieb Markus Metz markus.metz.gisw...@googlemail.com: On Sun, Mar 25, 2012 at 12:30 PM, Filipe Silva Dias filipesd...@gmail.com wrote: I switched the resolution from 30x30 to 100x100 and now it works. Problem solved. On Sun

Re: [GRASS-user] G_realloc error with v.voronoi

2012-03-19 Thread Markus Metz
On Mon, Mar 19, 2012 at 7:03 PM, spiderplant0 spiderpla...@gmail.com wrote: I have a shapefile containing points that I want to make into a voronoi diagram. I import the shapefile into GRASS and then run v.voronoi, but I then get this error... C:\Program Files (x86)\Quantum GIS

Re: [GRASS-user] How to create polygon from two selected contour line

2012-02-25 Thread Markus Metz
On Sat, Feb 25, 2012 at 7:52 PM, Michael Barton michael.bar...@asu.edu wrote: The trick for many things in GRASS is to think raster instead of thinking vector. Make a raster DEM from your map (if you don't already have one). Query the raster (e.g., with r.mapcalc or r.reclass) to create a

Re: [GRASS-user] The algorithm r.cost is not working

2012-02-23 Thread Markus Metz
On Thu, Feb 23, 2012 at 11:29 AM, Moritz Lennert mlenn...@club.worldonline.be wrote: On 23/02/12 07:19, Markus Metz wrote: On Wed, Feb 22, 2012 at 11:31 AM, Elisa Solano Villareal villar...@fbk.eu  wrote: Thanks for your answer. The cost surface used as input for r.cost must have only

Re: [GRASS-user] The algorithm r.cost is not working

2012-02-23 Thread Markus Metz
) == 0' failed. Can you try r.cost with percent_memory=30 (default is 100)? Markus M Please, can you give more information. Thanking you in advance, Elisa From: Markus Metz [markus.metz.gisw...@googlemail.com] Sent: Tuesday, February 14, 2012 12:39

Re: [GRASS-user] Buffered categories

2012-02-23 Thread Markus Metz
On Thu, Feb 23, 2012 at 8:14 PM, Hamish hamis...@yahoo.com wrote: Moritz wrote: If you create buffers around points which have category values and you would like to assign those values to the buffers, you can then (after v.category) do something like this: v.db.addtable yourbuffers

Re: [GRASS-user] The algorithm r.cost is not working

2012-02-22 Thread Markus Metz
: Markus Metz [markus.metz.gisw...@googlemail.com] Sent: Tuesday, February 14, 2012 12:39 PM To: Elisa Solano Villareal Cc: grass-user@lists.osgeo.org Subject: Re: [GRASS-user] The algorithm r.cost is not working On Tue, Feb 14, 2012 at 11:56 AM, Elisa Solano Villareal villar...@fbk.eu wrote

Re: [GRASS-user] get values for each point within lines

2012-02-15 Thread Markus Metz
Milton Cezar Ribeiro wrote: Dear Grass gurus I have a vector map with lines, and I need to get raster values each 5 meters. So I used v.split to spring my lines in segments of 5 meters, and after that I used v.to.points to generate points on each node of my line segments. Now I am trying

Re: [GRASS-user] How to remove small vectors in a shapefile in GRASS 6.4.1

2012-02-15 Thread Markus Metz
On Wed, Feb 15, 2012 at 5:46 PM, lucien.blandenier lucien.blanden...@unine.ch wrote: Hello, I want to remove small vector in a shp. I tried with v.generalize but I only managed to remove segment of a line. The image shows  in black the vectors I want to clean. And in green is the result I

Re: [GRASS-user] The algorithm r.cost is not working

2012-02-14 Thread Markus Metz
On Tue, Feb 14, 2012 at 11:56 AM, Elisa Solano Villareal villar...@fbk.eu wrote: Hi there, I have a question, I am working with r.cost algotithm but is not running for my whole spatial data but when a put just a piece of the map it works. The main problem is that the algorithm does not

Re: [GRASS-user] two questions: snap and close lines

2012-02-13 Thread Markus Metz
... Markus M Il 27/01/2012 18:53, Markus Metz ha scritto: The only two options I can think of is giving v.edit a try, otherwise hand-editing. Urgh! Too bad, thanks anyway. Why does breaking not work? Maybe you need remove duplicates and clean small angles at nodes. You could try v.clean

Re: [GRASS-user] How to link or combine adjacent vector lines

2012-02-12 Thread Markus Metz
On Sat, Feb 11, 2012 at 7:49 PM, Stefan Luedtke slued...@gfz-potsdam.de wrote: Sorry for the delay... and thanks for the fast answer. Adding a column and then running a command in the sense of v.dissolve is not an option because there are just to many. I was thinking of a way that combines

Re: [GRASS-user] v.generalize threshold and boyle's algorithm

2012-02-07 Thread Markus Metz
Johannes Radinger wrote: Hello, after some month a tried an script again which calls v.generalize. It seems that this module has changed in version 6.4.2 and 6.5SVN. Before I could use following: grass.run_command(v.generalize,                      input = input,                      

Re: [GRASS-user] v.generalize: strange behaviour on win?

2012-01-30 Thread Markus Metz
On Sun, Jan 29, 2012 at 1:12 PM, Paolo Cavallini cavall...@faunalia.it wrote: Il 28/01/2012 11:58, Markus Metz ha scritto: are left have mixed up attributes. That was fixed such that boundaries are not generalized if the generalization would damage topology. For OK, now I see. I think

Re: [GRASS-user] v.generalize: strange behaviour on win?

2012-01-28 Thread Markus Metz
On Sat, Jan 28, 2012 at 9:12 AM, Paolo Cavallini cavall...@faunalia.it wrote: Il 27/01/2012 19:24, Markus Metz ha scritto: v.generalize is broken in 6.4.1 oh, so it's the different version - I didn't know v.generalize was changed recently, thanks for letting me know. strangely enough

Re: [GRASS-user] two questions: snap and close lines

2012-01-27 Thread Markus Metz
On Fri, Jan 27, 2012 at 1:02 PM, Paolo Cavallini cavall...@faunalia.it wrote: Il 26/01/2012 20:41, Markus Metz ha scritto: On Mon, Jan 23, 2012 at 8:00 PM, Paolo Cavallini cavall...@faunalia.it wrote: Hi all. A couple of issues I ran into today: - I have a layer with some good and some bad

Re: [GRASS-user] v.generalize: strange behaviour on win?

2012-01-27 Thread Markus Metz
On Fri, Jan 27, 2012 at 1:32 PM, Paolo Cavallini cavall...@faunalia.it wrote: Il 26/01/2012 20:53, Markus Metz ha scritto: You would need to provide the input, not the output as an example. https://int.faunalia.it/~paolo/province.zip Considering the detail of the output and the threshold

Re: [GRASS-user] v.generalize: strange behaviour on win?

2012-01-27 Thread Markus Metz
On Fri, Jan 27, 2012 at 7:01 PM, Paolo Cavallini cavall...@faunalia.it wrote: Il 27/01/2012 19:02, Markus Metz ha scritto: same output. That is, if Province_general3 has been generated on windows, I do get the same result on Linux, no difference between the OS's. Are you using an older grass

Re: [GRASS-user] two questions: snap and close lines

2012-01-26 Thread Markus Metz
On Mon, Jan 23, 2012 at 8:00 PM, Paolo Cavallini cavall...@faunalia.it wrote: Hi all. A couple of issues I ran into today: - I have a layer with some good and some bad polygons; how do I import it into GRASS letting the bad ones snapping onto the good ones, and not vice versa? You could try

Re: [GRASS-user] v.generalize: strange behaviour on win?

2012-01-26 Thread Markus Metz
On Wed, Jan 25, 2012 at 11:03 AM, Paolo Cavallini cavall...@faunalia.it wrote: Hi all. Commands like: v.generalize input=province@Mapset_stefano type=area layer=1 -c type=boundary method=douglas threshold=10 look_ahead=7 reduction=50 slide=0.5 angle_thresh=3 degree_thresh=0

Re: [GRASS-user] v.edit tool=delete

2012-01-20 Thread Markus Metz
On Wed, Jan 18, 2012 at 8:32 PM, Martin Landa landa.mar...@gmail.com wrote: Hi, currently `v.edit` does not modify linked attribute tables when deleting features, eg. $ v.edit map=a tool=delete cat=1 $ v.db.select a where=cat=1 -c | wc -l 1 I am not sure if this behaviour is preferred

Re: [GRASS-user] How to rescale

2012-01-09 Thread Markus Metz
On Mon, Jan 9, 2012 at 8:42 PM, leonidas leonidas_lia...@yahoo.gr wrote: I created a raster using *r.grow.distance* I need to rescale that file to a scale 0-255 but in reverse order where high values in distance file will be converted in low values to the rescaled file. I'm using the

Re: [GRASS-user] v.water.outlet

2012-01-06 Thread Markus Metz
On Thu, Jan 5, 2012 at 12:16 PM, Paolo Cavallini cavall...@faunalia.it wrote: Il 04/01/2012 20:55, Micha Silver ha scritto: I guess you mean r.water.outlet? Yes, sorry. 1- The input must be the drainage *direction* map from r.watershed. ok 2- And it's crucial that the easting and

Re: [GRASS-user] v.water.outlet

2012-01-06 Thread Markus Metz
On Thu, Jan 5, 2012 at 12:12 PM, Massimiliano Cannata massimiliano.cann...@supsi.ch wrote: Moreover, check that your dem has sinks filled (r.fill.dir) Since the drainage map has been created by r.watershed, sink filling is not only not needed but does more harm than good and should not be used

Re: [GRASS-user] v.water.outlet

2012-01-06 Thread Markus Metz
On Fri, Jan 6, 2012 at 9:42 AM, Moritz Lennert mlenn...@club.worldonline.be wrote: On 06/01/12 09:23, Markus Metz wrote: On Thu, Jan 5, 2012 at 12:12 PM, Massimiliano Cannata massimiliano.cann...@supsi.ch  wrote: Moreover, check that your dem has sinks filled (r.fill.dir) Since

Re: [GRASS-user] using a vector map as target map for raster map to georectify

2012-01-04 Thread Markus Metz
On Mon, Oct 31, 2011 at 1:27 PM, Moritz Lennert mlenn...@club.worldonline.be wrote: On 31/10/11 11:34, Markus Metz wrote: On Mon, Oct 31, 2011 at 11:22 AM, Moritz Lennert mlenn...@club.worldonline.be  wrote: On 31/10/11 10:08, Markus Metz wrote: On Sat, Oct 29, 2011 at 11:25 AM, Moritz

Re: [GRASS-user] memory requirement of v.net.allpairs module

2011-12-27 Thread Markus Metz
: 1. r.grow by one cell  on area of interest. (It's an important step. Thanks to Markus Metz for this tip.) Time taken by the module: approximately 2 hours. 2. r.stream.extract on my area of interest.(42000*42000 cells) with 100 sq.km. as threshold (12345 in number of cells.) Above you said

Re: [GRASS-user] LIDAR tools and topology (unable to open ... on level 2)

2011-12-18 Thread Markus Metz
On Fri, Dec 16, 2011 at 7:03 AM, Micha Silver mi...@arava.co.il wrote: I find that after running each of the v.lidar.* the resulting vector lacks topology, so I can't view the categories: i.e. (after v.lidar.correct): GRASS 6.4.2RC1 (ITM):~/GIS/DEM/LIDAR_EinYahav v.category south_correct

<    4   5   6   7   8   9   10   11   12   13   >