Re: [GRASS-user] Moving windows on raster categorical map

2008-11-05 Thread Glynn Clements
Milton Cezar Ribeiro wrote: > I have a categorical raster map with 5 classes. > Now I need to estimate some stats using different-sized windows > like 3x3, 5x5..25x25. The stats the I would like > to compute is to count the number of classes on > the windows, to estimate the percentage of an > cl

Re: [GRASS-user] Moving windows on raster categorical map

2008-11-05 Thread Dylan Beaudette
On Wednesday 05 November 2008, Milton Cezar Ribeiro wrote: > Dear Grass-experts > > I have a categorical raster map with 5 classes. > Now I need to estimate some stats using different-sized windows > like 3x3, 5x5..25x25. The stats the I would like > to compute is to count the number of classes on

[GRASS-user] Moving windows on raster categorical map

2008-11-05 Thread Milton Cezar Ribeiro
Dear Grass-experts I have a categorical raster map with 5 classes. Now I need to estimate some stats using different-sized windows like 3x3, 5x5..25x25. The stats the I would like to compute is to count the number of classes on the windows, to estimate the percentage of an class of interest (like

Re: [GRASS-user] resampling and patching DEM

2008-11-05 Thread Glynn Clements
Silvia Simoni wrote: > I need to patch dems with different resolutions (1x1m) and (2x2 m) and > resample them to a coarser resolution (10x10 m). Which is the correct > procedure? > > a. resampling the 1x1 dem to 2x2 using r.resampl.stats -w (to be more > accurate) > patch them all at 2x2,

Re: [GRASS-user] Bug?: Projecting the same region from location to location and back gives a different result!!!!

2008-11-05 Thread Glynn Clements
Corrado wrote: > a very strange behaviour indeed: > > 1) location "A", one region defined TM, datum osgb36, region > s=0,w=0,e=80,n=130 > > 2) location "B", one region defined lat/long, wgs84 180W,180E,60N,60S > > in location A, I use v.in.region to draw a vector called "outbox", I lo

Re: A solution? (including resolution) - Re: [GRASS-user] Re: Help: Converting a raster map between locations (from wgs84 lat / long to UTM)

2008-11-05 Thread Glynn Clements
Corrado wrote: > > > does the r.proj resample the raster? > > > > > > The documentation seems to suggest you have to use g.region -m, or > > > may be I did not understand it properly. > > > > Like most raster commands, the output from r.proj has bounds and > > resolution which match the curr

Re: [GRASS-user] Rename multiple sqlite columns at once

2008-11-05 Thread Moritz Lennert
On 05/11/08 19:36, Nikos Alexandris wrote: On Wed, 2008-11-05 at 17:33 +0100, Moritz Lennert wrote: On 05/11/08 15:17, Nikos Alexandris wrote: On Tue, 2008-11-04 at 10:35 +0100, Paolo Craveri wrote: Hello Nikos The command v.db.renamecol map=test_rencol column="AREA","area" does not work. I

Re: [GRASS-user] Rename multiple sqlite columns at once

2008-11-05 Thread Nikos Alexandris
On Wed, 2008-11-05 at 19:51 +0100, Nikos Alexandris wrote: > On Wed, 2008-11-05 at 19:36 +0100, Nikos Alexandris wrote: > > On Wed, 2008-11-05 at 17:33 +0100, Moritz Lennert wrote: > > > On 05/11/08 15:17, Nikos Alexandris wrote: > > > > On Tue, 2008-11-04 at 10:35 +0100, Paolo Craveri wrote: > [..

Re: [GRASS-user] Rename multiple sqlite columns at once

2008-11-05 Thread Nikos Alexandris
On Wed, 2008-11-05 at 19:36 +0100, Nikos Alexandris wrote: > On Wed, 2008-11-05 at 17:33 +0100, Moritz Lennert wrote: > > On 05/11/08 15:17, Nikos Alexandris wrote: > > > On Tue, 2008-11-04 at 10:35 +0100, Paolo Craveri wrote: [...] I tried another example using "cut" this time. #!/bin/sh C

Re: [GRASS-user] Rename multiple sqlite columns at once

2008-11-05 Thread Nikos Alexandris
On Wed, 2008-11-05 at 17:33 +0100, Moritz Lennert wrote: > On 05/11/08 15:17, Nikos Alexandris wrote: > > On Tue, 2008-11-04 at 10:35 +0100, Paolo Craveri wrote: > >> Hello Nikos > >> > >>> The command > >>> v.db.renamecol map=test_rencol column="AREA","area" > >>> > >>> does not work. I assume tha

[GRASS-user] resampling and patching DEM

2008-11-05 Thread Silvia Simoni
Hi, I need to patch dems with different resolutions (1x1m) and (2x2 m) and resample them to a coarser resolution (10x10 m). Which is the correct procedure? a. resampling the 1x1 dem to 2x2 using r.resampl.stats -w (to be more accurate) patch them all at 2x2, resample the patched dem to 1

Re: [GRASS-user] Exporting raster map into unusual ascii file

2008-11-05 Thread Dylan Beaudette
On Wednesday 05 November 2008, Corrado wrote: > Dear friends, > > I have some 36 raster maps in one location, in the mapset PERMANENT: 12 for > the average monthly precipitation, 12 for the min temperature, 12 for the > max temperature. They are on the same region at the same resolution, with > nsr

Re: [GRASS-user] Rename multiple sqlite columns at once

2008-11-05 Thread Moritz Lennert
On 05/11/08 15:17, Nikos Alexandris wrote: On Tue, 2008-11-04 at 10:35 +0100, Paolo Craveri wrote: Hello Nikos The command v.db.renamecol map=test_rencol column="AREA","area" does not work. I assume that trying to use 'tr' for example wouldn't help (?). Or am I doing something wrong? Try t

Re: [GRASS-user] inner and outer buffer for a boundary

2008-11-05 Thread Moritz Lennert
On 05/11/08 14:46, Martin Landa wrote: Hi, 2008/11/5 Moritz Lennert <[EMAIL PROTECTED]>: Sounds like the lines do not have category values associated to them (which is done automatically during import), so you probably need to run v.category, before running v.buffer. is there any reason why v

Re: [GRASS-user] i.pca vs. r.covar/m.eigensystem/r.mapcalc

2008-11-05 Thread Dylan Beaudette
On Wed, Nov 5, 2008 at 5:39 AM, Nikos Alexandris <[EMAIL PROTECTED]> wrote: > On Tue, 2008-10-28 at 09:39 -0700, Dylan Beaudette wrote: >> On Monday 27 October 2008, Nikos Alexandris wrote: >> > My apologies for BUMPing :-) >> > >> > On Sun, 2008-10-12 at 00:20 +0200, Nikos Alexandris wrote: >> > >

[GRASS-user] Another off-topic question about sqlite, MySQL, Postgresql

2008-11-05 Thread Nikos Alexandris
Running Amarok (music player) for the first time asks which database to use (to store music info). The "wizard" dialog states: MySQL and Postgresql are faster than sqlite, but require additional setup. I though that sqlite is faster (but with some multi-user access limitations). So, what is faster

Re: [GRASS-user] Exporting raster maps into structured ascii file (X, Y, value, ...., value)

2008-11-05 Thread Nikos Alexandris
On Wed, 2008-11-05 at 10:38 +, Corrado wrote: > Dear friends, > > I have some 36 raster maps in one location, in the mapset PERMANENT: 12 for > the average monthly precipitation, 12 for the min temperature, 12 for the max > temperature. They are on the same region at the same resolution, wit

Re: [GRASS-user] Rename multiple sqlite columns at once

2008-11-05 Thread Nikos Alexandris
On Tue, 2008-11-04 at 10:35 +0100, Paolo Craveri wrote: > Hello Nikos > > > The command > > v.db.renamecol map=test_rencol column="AREA","area" > > > > does not work. I assume that trying to use 'tr' for example wouldn't > > help (?). Or am I doing something wrong? > > > > > Try this: > > > ==

Re: [GRASS-user] inner and outer buffer for a boundary

2008-11-05 Thread Otto Dassau
Hi Moritz, On Wed, 05 Nov 2008 14:25:24 +0100 Moritz Lennert <[EMAIL PROTECTED]> wrote: > On 05/11/08 13:46, Otto Dassau wrote: > > Hi, > > > > I tried to create a bufferzone (inner and outer area) around an area > > boundary > > using following steps (reproduced with spearfish quads data): > >

Re: [GRASS-user] inner and outer buffer for a boundary

2008-11-05 Thread Martin Landa
Hi, 2008/11/5 Moritz Lennert <[EMAIL PROTECTED]>: > Sounds like the lines do not have category values associated to them (which > is done automatically during import), so you probably need to run > v.category, before running v.buffer. is there any reason why v.buffer should ignore features withou

Re: [GRASS-user] i.pca vs. r.covar/m.eigensystem/r.mapcalc

2008-11-05 Thread Nikos Alexandris
On Tue, 2008-10-28 at 09:39 -0700, Dylan Beaudette wrote: > On Monday 27 October 2008, Nikos Alexandris wrote: > > My apologies for BUMPing :-) > > > > On Sun, 2008-10-12 at 00:20 +0200, Nikos Alexandris wrote: > > > I wanted to check/verify some i.pca results so I calculated 2 principal > > > comp

Re: [GRASS-user] inner and outer buffer for a boundary

2008-11-05 Thread Moritz Lennert
On 05/11/08 13:46, Otto Dassau wrote: Hi, I tried to create a bufferzone (inner and outer area) around an area boundary using following steps (reproduced with spearfish quads data): # change boundary to line v.type input=quads output=quadsline type=boundary,line # remove centroids from line map

[GRASS-user] Bug?: Projecting the same region from location to location and back gives a different result!!!!

2008-11-05 Thread Corrado
Dear friends, a very strange behaviour indeed: 1) location "A", one region defined TM, datum osgb36, region s=0,w=0,e=80,n=130 2) location "B", one region defined lat/long, wgs84 180W,180E,60N,60S in location A, I use v.in.region to draw a vector called "outbox", I load it in location

[GRASS-user] inner and outer buffer for a boundary

2008-11-05 Thread Otto Dassau
Hi, I tried to create a bufferzone (inner and outer area) around an area boundary using following steps (reproduced with spearfish quads data): # change boundary to line v.type input=quads output=quadsline type=boundary,line # remove centroids from line map v.edit map=quadsline type=centroid tool

Re: [GRASS-user] Exporting raster maps into structured ascii file (X, Y, value, ...., value)

2008-11-05 Thread Daniel Victoria
Hi, Take a look at the r.what man page... http://grass.osgeo.org/gdp/html_grass63/r.what.html I did something similar but for a few points. What you will have to do is create a list with the x,y coordinates that you want and then pass it to r.what You could create your x,y list by hand or, using

[GRASS-user] Exporting raster map into unusual ascii file

2008-11-05 Thread Corrado
Dear friends, I have some 36 raster maps in one location, in the mapset PERMANENT: 12 for the average monthly precipitation, 12 for the min temperature, 12 for the max temperature. They are on the same region at the same resolution, with nsres = ewres (2000). I would like to ouput (export) the

[GRASS-user] Exporting raster maps into structured ascii file (X, Y, value, ...., value)

2008-11-05 Thread Corrado
Dear friends, I have some 36 raster maps in one location, in the mapset PERMANENT: 12 for the average monthly precipitation, 12 for the min temperature, 12 for the max temperature. They are on the same region at the same resolution, with nsres = ewres (2000). I would like to ouput (export) the

Re: A solution? (including resolution) - Re: [GRASS-user] Re: Help: Converting a raster map between locations (from wgs84 lat / long to UTM)

2008-11-05 Thread Corrado
On Tuesday 04 November 2008 19:17:23 Glynn Clements wrote: > Corrado wrote: > > does the r.proj resample the raster? > > > > The documentation seems to suggest you have to use g.region -m, or > > may be I did not understand it properly. > > Like most raster commands, the output from r.proj has

Re: [GRASS-user] v.select output vector is empty

2008-11-05 Thread Moritz Lennert
On 04/11/08 17:57, Kevin Webb wrote: Thank you for your reply Hamish. 2 questions based on your suggestion... At 11:18 PM 10/29/2008, Hamish wrote: Kevin Webb wrote: > I am attempting to execute a "point-in-polygon" process on 2 vectors using > v.select and the output vector is empty, it co