[GRASS-user] r.coin behaving differently

2010-04-14 Thread maning sambale
Hi,

I am using r.coin for two raster maps

map1's range is
Range of data:min = 111  max = 140

map2' range
Range of data:min =   max = 17590


I used to get this result with r.coin
|   cat# |   111 |   112 |   120 |   130 |   140 |
  w cat 0 | w/o cat 0 |
||
|p  |  0.00 |  0.00 |  0.00 |  0.00 |
103.33 |103.33 |103.33 |
|h 1280 |   8376.26 |  51848.35 |395.22 |   5594.52 |
36.80 |  66251.15 |  66251.15 |
|i 1290 |  0.00 |  10521.99 |   3791.04 |636.42 |
162.64 |  15112.09 |  15112.09 |
|l 1330 |  0.00 | 66.13 |577.76 |  0.00 |
141.35 |785.24 |785.24 |
|p 1550 |   2079.69 |  12411.01 |440.71 |   1594.05 |
696.53 |  17221.99 |  17221.99 |
|r 2090 | 51.59 |   3862.45 |  0.00 |  0.00 |
0.00 |   3914.03 |   3914.03 |
|o 2150 | 112366.40 | 197525.87 |  0.00 |  17049.36 |
5081.42 | 332023.05 | 332023.05 |

Now I get a different cat# as row headers
|   cat# |   111 |   112 |   120 |   130 |   140 |
  w cat 0 | w/o cat 0 |
||
|p 1 |  0.00 |  0.00 |  0.00 |  0.00 |103.33 |
   103.33 |103.33 |
|h 3 |   8376.26 |  62370.34 |   4186.26 |   6230.93 |199.44 |
 81363.24 |  81363.24 |
|i 4 |  0.00 | 66.13 |577.76 |  0.00 |141.35 |
   785.24 |785.24 |
|l 7 |   2079.69 |  12411.01 |440.71 |   1594.05 |696.53 |
 17221.99 |  17221.99 |

Any idea why?  I'm running GRASS 6.4.0RC5 from kingchaos mac osx

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: r.coin behaving differently

2010-04-14 Thread maning sambale
Follow-up

On Wed, Apr 14, 2010 at 5:41 PM, maning sambale
emmanuel.samb...@gmail.com wrote:
 map2' range
 Range of data:    min =   max = 17590

The range above is not are ordinal categories, they are simply polygon
codes of administrative regions.

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Questions about r.mfilter

2010-04-14 Thread Glynn Clements

Casey Vandenberg wrote:

 Does anyone know if the r.mfilter will support negative values in the 
 matrix filter file. I was hoping to do a high pass filter on a dataset, 
 so I created a filter for input in the following format, but it will not 
 work.
 
 TITLE   HighPass3x3
 MATRIX3
 -1 -1 -1
 -1 8 -1
 -1 -1 -1
 DIVISOR   9
 TYPES
 
 Input command is:
 r.mfilter input=mag_1...@laindia output=HighPass3x3_Mag12.5 
 filter=/home/cvandenberg/Documents/Filters/HighPass3x3.txt*
 **ERROR: Filter file format error*

This error message occurs if the DIVISOR, TYPE or START keywords are
encountered before a matrix has been defined.

Also, a sequential filter (TYPE S) is incorrect for most cases
(including this one).

There aren't any issues with negative numbers.

Note that you can also use r.neighbors weight=... method=sum to
apply a matrix filter.

-- 
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] Running v.krige in GRASS6.5

2010-04-14 Thread António Rocha

Greetings

I have just installed GRASS6.5svn and I have the following files
v.krige.py in sr/local/grass-6.5.svn/etc/wxpython/scripts
v.krige.pyc in /usr/local/grass-6.5.svn/etc/wxpython/scripts

But in Wxpython command console, everytime I put v.krige I only get:
(Wed Apr 14 11:48:10 
2010) 
v.krige 


(Wed Apr 14 11:48:10 2010) Command finished (0 sec)

Can anyone tell me what might be wrong? Or what am I doing wrong? :)

Thanks
Antonio


__ Information from ESET NOD32 Antivirus, version of virus signature 
database 5027 (20100414) __

The message was checked by ESET NOD32 Antivirus.

http://www.eset.com


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Difficulties with i.pr

2010-04-14 Thread Daniel McInerney

Hi Joel -
Sorry for the late reply. I only managed
to successfully run i.pr.training and i.pr.features.
I tried the image classification, but it
threw too many errors.

I haven't looked at i.pr for a while, but
may revisit it when time permits.

Best Regards,
Daniel.

joel dinis wrote:

Hello Daniel,

Thank you for the thread.
I'm working with polygons.

Did you classify any image with
the i.pr?

JD

On Mon, Apr 12, 2010 at 4:36 PM, Daniel McInerney
daniel.mciner...@ucd.iewrote:


Hi Joel,
I did some work with i.pr a while back,
but was unsuccessful in getting it to
fully work.

Is your vector file a point or polygon?
I tested i.pr.training with points.

You might get some ideas from this thread:
http://www.mail-archive.com/grass-user@lists.osgeo.org/msg00153.html

Regards,
Daniel.

joel dinis wrote:


Hello to all,

Does anyone knows any working
example of the i.pr add-on?

I have install it, and I'm having
some problems with the i.pr.training.

I do:
i.pr.training map=lis...@p9,lis...@p9,lis...@p9,lis...@p9training=iprteste9
vector=area1 rows=3 cols=3 class=0

And the iprteste9 only has this:
Data type:
GrassTraining
Number of layers:
4
Label:
/home/joel/iprteste9
Data:
Layer_1Layer_2Layer_3Layer_4ClassEastNorth
RowsColsEW-resNS-res

And nothing else.
I suppose that the file would have something more.

Can anyone help me?





___
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] Running v.krige in GRASS6.5

2010-04-14 Thread Martin Landa
Hi,

2010/4/14 António Rocha antonio.ro...@deimos.com.pt:
 I have just installed GRASS6.5svn and I have the following files
 v.krige.py     in sr/local/grass-6.5.svn/etc/wxpython/scripts
 v.krige.pyc     in /usr/local/grass-6.5.svn/etc/wxpython/scripts

 But in Wxpython command console, everytime I put v.krige I only get:
 (Wed Apr 14 11:48:10 2010)
   v.krige

 (Wed Apr 14 11:48:10 2010) Command finished (0 sec)

 Can anyone tell me what might be wrong? Or what am I doing wrong? :)

should be fixed r41859.

Martin

-- 
Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Link to Open office Database

2010-04-14 Thread Francesco Barreca
I want link a point in a map with a Mask of OpenOffice database. How am
I doing?
thanks.

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Running v.krige in GRASS6.5

2010-04-14 Thread Martin Landa
Hi,

2010/4/14 António Rocha antonio.ro...@deimos.com.pt:
 1- Is it possible to back-port it to GRASS6.4?

v.krige's command line interface should work in GRASS 6.4, the GUI
part requires GRASS 6.5+.

 2- How can I know with release I'm using (rX)?

Depends how your GRASS is installed (binary snapshots, distro packages, svn...).

Martin

-- 
Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Re: GRASS - Network Analysis

2010-04-14 Thread G. Allegri
Great Martin, I've just tested it!
What do you think about adding the option to parametrize inputs and outputs?

giovanni

2010/4/10 Martin Landa landa.mar...@gmail.com

 Hi,

 2010/4/10 Martin Landa landa.mar...@gmail.com:
  the modules developed by Daniel Bundala within Google Summer of Code
  2009 is now available in GRASS 7.0 and GRASS 6.5 (r41793).

 I forgot to mention that in GRASS 6.5. you need to run ./configure
 before compilation to generate new Grass.make file.

 Martin

 --
 Martin Landa landa.martin gmail.com * 
 http://gama.fsv.cvut.cz/~landahttp://gama.fsv.cvut.cz/%7Elanda
 ___
 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] Running v.krige in GRASS6.5

2010-04-14 Thread Martin Landa
(please keep discussion in ML)

Hi,

2010/4/14 António Rocha antonio.ro...@deimos.com.pt:
 Regarding first question
 This means that it won't be available or adapted to 6.4 right?

could be in 6.4.1, unfortunately adopting v.krige to wxGUI available
in GRASS 6.4 would be too much work. Anyway v.krige's command line
interface should work after little changes.

 I'm using weekly snapshots. So I have to search for that rX

in this case wait for the next snapshot.

Martin

-- 
Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Re: GRASS - Network Analysis

2010-04-14 Thread G. Allegri
Ops, Martin, I've replied using the wrong thread! :)

thanks,
giovanni

2010/4/14 Martin Landa landa.mar...@gmail.com

 Hi,

 2010/4/14 G. Allegri gioha...@gmail.com:
  Great Martin, I've just tested it!
  What do you think about adding the option to parametrize inputs and
 outputs?

 seems that you are speaking about the modeler. Parametrization of
 inputs/outputs/variables is in TODO.

 Martin

 --
 Martin Landa landa.martin gmail.com * 
 http://gama.fsv.cvut.cz/~landahttp://gama.fsv.cvut.cz/%7Elanda

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] r.mapcalc and nested if

2010-04-14 Thread thedok78

Hello everybody,
I need to run a complex r.mapcalc expression, but I can't understand how to
use nested if.
This is the expression:
r.mapcalc
'map2=if(mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1])==1,0,mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1]))'

The idea is this:
if the cell value is 1, do the mode else skip that cell; 
if the mode is equal to 1 assign 0 to that cell else assign the mode value.

The r.mapcalc expression I just posted produce wrong result..
Any idea?

Thanks
Luca
-- 
View this message in context: 
http://n2.nabble.com/r-mapcalc-and-nested-if-tp4901916p4901916.html
Sent from the Grass - Users mailing list archive at Nabble.com.
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: r.mapcalc and nested if

2010-04-14 Thread thedok78

Sorry, this is the expression:
r.mapcalc 'map2=if(map ==1
,if(mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1])==1,0,mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1]),map))'

Thanks
Luca
-- 
View this message in context: 
http://n2.nabble.com/r-mapcalc-and-nested-if-tp4901916p4901925.html
Sent from the Grass - Users mailing list archive at Nabble.com.
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Stratified sampling in GRASS

2010-04-14 Thread joel dinis
Hello to all,

I'm trying to do a stratified
sampling in GRASS, and I found
nothing in GRASS about it.

Does anyone knows how to do it?

Thank in advance
JD
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Stratified sampling in GRASS

2010-04-14 Thread Dylan Beaudette
On Wednesday 14 April 2010, joel dinis wrote:
 Hello to all,

 I'm trying to do a stratified
 sampling in GRASS, and I found
 nothing in GRASS about it.

 Does anyone knows how to do it?

 Thank in advance
 JD

Here is one approach with R:

http://casoilresource.lawr.ucdavis.edu/drupal/node/644

This example can be modified to read directly from a GRASS database with the 
spgrass6 package for R.

Cheers,
Dylan

-- 
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Error in GRASS6.4 while creating image group

2010-04-14 Thread Nikos Dumakis
Hi
I'm using GRASS6.4.svn from latest week snapshot in Ubuntu and I was testing
a few functions while I got an error when creating a image group with
NorthCarolina Dataset.

I did:
i.group group=gr04 subgroup=sub input=lsat7_2000...@landsat
,lsat7_2000...@landsat,lsat7_2000...@landsat,lsat7_2000...@landsat
,lsat7_2000...@landsat,lsat7_2000...@landsat

And I got:

Adding raster map lsat7_2000...@landsat to group
Adding raster map lsat7_2000...@landsat to group
Adding raster map lsat7_2000...@landsat to group
Adding raster map lsat7_2000...@landsat to group
Adding raster map lsat7_2000...@landsat to group
Adding raster map lsat7_2000...@landsat to group
Adding raster map lsat7_2000...@landsat to subgroup
Adding raster map lsat7_2000...@landsat to subgroup
Adding raster map lsat7_2000...@landsat to subgroup
Adding raster map lsat7_2000...@landsat to subgroup
Adding raster map lsat7_2000...@landsat to subgroup
Adding raster map lsat7_2000...@landsat to subgroup

Then when I open i.cluster GUI I selected gr04 group and, sub doesn't appear
in subgroup. It only appears Not a selectable element. If I leave it blank
it doesn't run (sub group is a required option). But, if I write  sub it
works...

Is this suppose to happen?

Nikos
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] temporary maps/space during r.fill.dir execution

2010-04-14 Thread Christian Schwartze
Hi grass users,

does anybody know how to avoid the temporary use of disk space during execution 
of r.fill.dir? For large input files it seems that large temporary files are 
written as well since I got messages like no disk space left (running in a 
minimized environment where not so much disk space is available).

Thanks a lot for any hints...

Regards,
Christian.___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Error with sub-group while using i.gensigset

2010-04-14 Thread Luisa Peña
Greetings
I'm freely following this tutorial (
http://grass.osgeo.org/wiki/Image_classification) and I'm getting a odd
error.
I'm sending this command:
i.gensigset trainingma...@landsat group=g...@landsat subgroup=sub
signaturefile=sig maxsig=5

and I get this:
Finding training classes...
5 classes found
Reading raster maps...
Clustering class 1 (1753 pixels)...
Number of subclasses is 5
Clustering class 2 (2810 pixels)...
Number of subclasses is 5
Clustering class 3 (2772 pixels)...
Number of subclasses is 5
Clustering class 4 (2177 pixels)...
Number of subclasses is 5
Clustering class 5 (6060 pixels)...
Number of subclasses is 5
unable to create signature file sig for subgroup sub of group g...@landsat
ERROR: Unable to create signature file sig

I think it might have somethibg to do with i.group because everytime I open
i.group the followiong is printed to the Command output:
Invalid element 'g...@landsat@landsat'
Invalid element 'g...@landsat@landsat'
Invalid element 'g...@landsat@landsat'
Invalid element 'g...@landsat@landsat'
All these are groups that I have made

If I do  group -l group=g...@landsat
I get:

group gr06 references the following raster maps
-
lsat7_2000...@landsatlsat7_2000...@landsat
lsat7_2000...@landsatlsat7_2000...@landsat
lsat7_2000...@landsat
-
i.group complete.

What might be wrong?
Thanks
Luisa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Re: r.mapcalc and nested if

2010-04-14 Thread Jarek Jasiewicz

thedok78 pisze:

Sorry, this is the expression:
r.mapcalc 'map2=if(map ==1
,if(mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1])==1,0,mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1]),map))'

Thanks
Luca
  

seem that much easier will be:

r.neighbors input=your_map method=mode size=3 output=mode_output
r.mapcalc 'first_map=if(your_map==1,output_mode,your_map)'
r.mapclac 'final=if(first_map=1,0,first_map'

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Pango ERROR

2010-04-14 Thread António Rocha

Hi

I'm using GRASS6.4svn in an Ubuntu Machine and, when I start GRASS I get 
this warning at command line:
(python:3305): Pango-CRITICAL **: pango_layout_get_cursor_pos: assertion 
`index = 0  index = layout-length' failed

What might be causing this? and what are the consequences of this?

Thanks

Best regards
Antonio



__ Information from ESET NOD32 Antivirus, version of virus signature 
database 5029 (20100414) __

The message was checked by ESET NOD32 Antivirus.

http://www.eset.com


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Re: r.mapcalc and nested if

2010-04-14 Thread Milton Cezar Ribeiro
Hi there,

may be the last command need double ==

r.mapclac 'final=if(first_map==1,0,first_map'

cheers

milton

2010/4/14 Jarek Jasiewicz jar...@amu.edu.pl

 thedok78 pisze:

 Sorry, this is the expression:
 r.mapcalc 'map2=if(map ==1

 ,if(mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1])==1,0,mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1]),map))'

 Thanks
 Luca


 seem that much easier will be:

 r.neighbors input=your_map method=mode size=3 output=mode_output
 r.mapcalc 'first_map=if(your_map==1,output_mode,your_map)'
 r.mapclac 'final=if(first_map=1,0,first_map'


 ___
 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] Error with sub-group while using i.gensigset

2010-04-14 Thread Markus Neteler
On Wed, Apr 14, 2010 at 7:41 PM, Luisa Peña luisapena1...@gmail.com wrote:
 Greetings
 I'm freely following this tutorial
 (http://grass.osgeo.org/wiki/Image_classification) and I'm getting a odd
 error.
 I'm sending this command:
 i.gensigset trainingma...@landsat group=g...@landsat subgroup=sub
 signaturefile=sig maxsig=5

Do you really need all @landsat indications?

 and I get this:
 Finding training classes...
 5 classes found
 Reading raster maps...
 Clustering class 1 (1753 pixels)...
 Number of subclasses is 5
 Clustering class 2 (2810 pixels)...
 Number of subclasses is 5
 Clustering class 3 (2772 pixels)...
 Number of subclasses is 5
 Clustering class 4 (2177 pixels)...
 Number of subclasses is 5
 Clustering class 5 (6060 pixels)...
 Number of subclasses is 5
 unable to create signature file sig for subgroup sub of group g...@landsat
 ERROR: Unable to create signature file sig

 I think it might have somethibg to do with i.group because everytime I open
 i.group the followiong is printed to the Command output:
 Invalid element 'g...@landsat@landsat'
 Invalid element 'g...@landsat@landsat'
 Invalid element 'g...@landsat@landsat'
 Invalid element 'g...@landsat@landsat'
 All these are groups that I have made

 If I do  group -l
 group=g...@landsat   I get:

 group gr06 references the following raster maps
 -
 lsat7_2000...@landsat    lsat7_2000...@landsat
 lsat7_2000...@landsat    lsat7_2000...@landsat
 lsat7_2000...@landsat
 -
 i.group complete.

 What might be wrong?

I dunno but perhaps the @landsat isn't properly stripped off in some cases.
Could you retry to make the group etc in the landsat mapset itself in
order to understand if @landsat makes the difference?

Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Error in GRASS6.4 while creating image group

2010-04-14 Thread Markus Neteler
On Wed, Apr 14, 2010 at 7:05 PM, Nikos Dumakis nikosdu1...@gmail.com wrote:
 Hi
 I'm using GRASS6.4.svn from latest week snapshot in Ubuntu and I was testing
 a few functions while I got an error when creating a image group with
 NorthCarolina Dataset.

 I did:
 i.group group=gr04 subgroup=sub
 input=lsat7_2000...@landsat,lsat7_2000...@landsat,lsat7_2000...@landsat,lsat7_2000...@landsat,lsat7_2000...@landsat,lsat7_2000...@landsat

 And I got:

 Adding raster map lsat7_2000...@landsat to group
 Adding raster map lsat7_2000...@landsat to group
 Adding raster map lsat7_2000...@landsat to group
 Adding raster map lsat7_2000...@landsat to group
 Adding raster map lsat7_2000...@landsat to group
 Adding raster map lsat7_2000...@landsat to group
 Adding raster map lsat7_2000...@landsat to subgroup
 Adding raster map lsat7_2000...@landsat to subgroup
 Adding raster map lsat7_2000...@landsat to subgroup
 Adding raster map lsat7_2000...@landsat to subgroup
 Adding raster map lsat7_2000...@landsat to subgroup
 Adding raster map lsat7_2000...@landsat to subgroup

 Then when I open i.cluster GUI I selected gr04 group and, sub doesn't appear
 in subgroup. It only appears Not a selectable element. If I leave it blank
 it doesn't run (sub group is a required option). But, if I write  sub it
 works...

 Is this suppose to happen?

I would say no. It reminds me of
http://trac.osgeo.org/grass/ticket/780

Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] temporary maps/space during r.fill.dir execution

2010-04-14 Thread Christian Schwartze

g.region -p says:

   projection: 0 (x,y)
   zone: 0
   north: 6257261.73938969
   south: 6124108.23119424
   west: 650218.74629604
   east: 775529.77427479
   nsres: 84.32774427
   ewres: 69.6946763
   rows: 1579
   cols: 1798
   cells: 2839042

Hence, working on nearly 3 mill. cells...

Regards,
Christian.
--
From: Markus Neteler nete...@osgeo.org
Sent: Wednesday, April 14, 2010 10:07 PM
To: Christian Schwartze christian.schwar...@uni-jena.de
Cc: GRASS user list grass-user@lists.osgeo.org
Subject: Re: [GRASS-user] temporary maps/space during r.fill.dir execution


On Wed, Apr 14, 2010 at 11:01 AM, Christian Schwartze
christian.schwar...@uni-jena.de wrote:

Hi grass users,

does anybody know how to avoid the temporary use of disk space during
execution of r.fill.dir? For large input files it seems that large 
temporary

files are written as well since I got messages like no disk space left
(running in a minimized environment where not so much disk space is
available).


The command is actually memory demanding. What's the output of
g.region -p
? On how many cells do you work on?

Markus 


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] temporary maps/space during r.fill.dir execution

2010-04-14 Thread Markus Neteler
On Wed, Apr 14, 2010 at 11:01 AM, Christian Schwartze
christian.schwar...@uni-jena.de wrote:
 Hi grass users,

 does anybody know how to avoid the temporary use of disk space during
 execution of r.fill.dir? For large input files it seems that large temporary
 files are written as well since I got messages like no disk space left
 (running in a minimized environment where not so much disk space is
 available).

The command is actually memory demanding. What's the output of
g.region -p
? On how many cells do you work on?

Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] r.average and r.info

2010-04-14 Thread William Hudspeth
Hello,

I am using r.average to calculate the mean pixel value in a base map
aggregating pixels in a cover map with the same value. I know that I
have 54 different categories (e.g. parcels) in the cover map, but when I
run r.info, it only outputs 47 of them, with a  on the final line.
Where are the rest of the results?

Thanks, Bill

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Error in GRASS6.4 while creating image group

2010-04-14 Thread Nikos Dumakis
Does anyone reports the same? because that happened with Win32 and I'm using
a snapshot from 10/4/2010 in Ubuntu9.10.
It happens the same in i.cluster and in all the others menus that required
groups...
On Wed, Apr 14, 2010 at 9:11 PM, Markus Neteler nete...@osgeo.org wrote:

  On Wed, Apr 14, 2010 at 7:05 PM, Nikos Dumakis nikosdu1...@gmail.com
 wrote:
  Hi
  I'm using GRASS6.4.svn from latest week snapshot in Ubuntu and I was
 testing
  a few functions while I got an error when creating a image group with
  NorthCarolina Dataset.
 
  I did:
  i.group group=gr04 subgroup=sub
  input=lsat7_2000...@landsat,lsat7_2000...@landsat,lsat7_2000...@landsat
 ,lsat7_2000...@landsat,lsat7_2000...@landsat,lsat7_2000...@landsat
 
  And I got:
 
  Adding raster map lsat7_2000...@landsat to group
  Adding raster map lsat7_2000...@landsat to group
  Adding raster map lsat7_2000...@landsat to group
  Adding raster map lsat7_2000...@landsat to group
  Adding raster map lsat7_2000...@landsat to group
  Adding raster map lsat7_2000...@landsat to group
  Adding raster map lsat7_2000...@landsat to subgroup
  Adding raster map lsat7_2000...@landsat to subgroup
  Adding raster map lsat7_2000...@landsat to subgroup
  Adding raster map lsat7_2000...@landsat to subgroup
  Adding raster map lsat7_2000...@landsat to subgroup
  Adding raster map lsat7_2000...@landsat to subgroup
 
  Then when I open i.cluster GUI I selected gr04 group and, sub doesn't
 appear
  in subgroup. It only appears Not a selectable element. If I leave it
 blank
  it doesn't run (sub group is a required option). But, if I write  sub
 it
  works...
 
  Is this suppose to happen?

 I would say no. It reminds me of
 http://trac.osgeo.org/grass/ticket/780

 Markus

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] temporary maps/space during r.fill.dir execution

2010-04-14 Thread Martin Landa
Hi,

2010/4/14 Markus Neteler nete...@osgeo.org:
 does anybody know how to avoid the temporary use of disk space during
 execution of r.fill.dir? For large input files it seems that large temporary
 files are written as well since I got messages like no disk space left
 (running in a minimized environment where not so much disk space is
 available).

 The command is actually memory demanding. What's the output of

btw, is there any particular reason why data are kept in temp files
and then written to raster maps?

Martin

-- 
Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] isolines and kml files interpolation

2010-04-14 Thread Ricardo Rodríguez
hello all, a question, points are interpolated by any method in kml and get
a new file with interpolation or isolines in kml,
thanks for your answers

Ricardo Rodríguez
Ing. Topográfico
Univalle
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Error with sub-group while using i.gensigset

2010-04-14 Thread Luisa Peña
I have re-downloaded NC_spm_08 from GRASS website and the already available
group,  lsat7_2...@landsat doesn't display sub-group in i.cluster or other
functions..
Nevertheless, I tried what Markus said:
i.group group=t02 subgroup=sub input=lsat5_1987_10,lsat5_1987_20 (altough
this is not performed by GUI but by me on the command output)
Adding raster map lsat5_1987...@landsat to group
Adding raster map lsat5_1987...@landsat to group
Adding raster map lsat5_1987...@landsat to subgroup
Adding raster map lsat5_1987...@landsat to subgroup
i.group complete.

And then at i.cluster when I select t02, in subgroup there is no sub-group
only NOT SELECTABLE ELEMENT.

What am I missing/or might be wrong?
Luisa

2010/4/14 Markus Neteler nete...@osgeo.org

 On Wed, Apr 14, 2010 at 7:41 PM, Luisa Peña luisapena1...@gmail.com
 wrote:
  Greetings
  I'm freely following this tutorial
  (http://grass.osgeo.org/wiki/Image_classification) and I'm getting a odd
  error.
  I'm sending this command:
  i.gensigset trainingma...@landsat group=g...@landsat subgroup=sub
  signaturefile=sig maxsig=5

 Do you really need all @landsat indications?

  and I get this:
  Finding training classes...
  5 classes found
  Reading raster maps...
  Clustering class 1 (1753 pixels)...
  Number of subclasses is 5
  Clustering class 2 (2810 pixels)...
  Number of subclasses is 5
  Clustering class 3 (2772 pixels)...
  Number of subclasses is 5
  Clustering class 4 (2177 pixels)...
  Number of subclasses is 5
  Clustering class 5 (6060 pixels)...
  Number of subclasses is 5
  unable to create signature file sig for subgroup sub of group
 g...@landsat
  ERROR: Unable to create signature file sig
 
  I think it might have somethibg to do with i.group because everytime I
 open
  i.group the followiong is printed to the Command output:
  Invalid element 'g...@landsat@landsat'
  Invalid element 'g...@landsat@landsat'
  Invalid element 'g...@landsat@landsat'
  Invalid element 'g...@landsat@landsat'
  All these are groups that I have made
 
  If I do  group -l
  group=g...@landsat   I
 get:
 
  group gr06 references the following raster maps
  -
  lsat7_2000...@landsatlsat7_2000...@landsat
  lsat7_2000...@landsatlsat7_2000...@landsat
  lsat7_2000...@landsat
  -
  i.group complete.
 
  What might be wrong?

 I dunno but perhaps the @landsat isn't properly stripped off in some cases.
 Could you retry to make the group etc in the landsat mapset itself in
 order to understand if @landsat makes the difference?

 Markus

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Can't initialize wxpython GRASS6.4

2010-04-14 Thread Jenny Turner
Greetings
I have just installed in my old laptop a Ubuntu Karmic version of GRASS
(6.4). I know it's not new but I'm getting an ODD error that it's not
suppose to happen.
I have placed nc_spm_08 data in /home/jenny/data folder and everytime I do
sudo grass64 -wxpython
When it's open GRASS GIS layer Manager I get an error window  stating:
Unable to get current geographic extent. Force quiting wxGUI. Please run
manually g.region to fix the problem.

This is weird because it's suppose to work right? Any ideas?
Jenny
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Can't initialize wxpython GRASS6.4

2010-04-14 Thread Martin Landa
Hi,

2010/4/15 Jenny Turner jennyturner1...@gmail.com:

 I have just installed in my old laptop a Ubuntu Karmic version of GRASS
 (6.4). I know it's not new but I'm getting an ODD error that it's not
 suppose to happen.
 I have placed nc_spm_08 data in /home/jenny/data folder and everytime I do
 sudo grass64 -wxpython
 When it's open GRASS GIS layer Manager I get an error window  stating:
 Unable to get current geographic extent. Force quiting wxGUI. Please run
 manually g.region to fix the problem.

 This is weird because it's suppose to work right? Any ideas?

try to run `g.region -p` from command line. What it says?

Martin

-- 
Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Can't initialize wxpython GRASS6.4

2010-04-14 Thread Martin Landa
Hi,

2010/4/15 Jenny Turner jennyturner1...@gmail.com:
 It says
 g.region: error while loading shared libraries: libgdal.so.1: cannot open
 shared object file: No such file or directory

 But I have gdal installed and libgdal and libgdal1-dev

try `sudo ldconfig`. Martin

-- 
Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Question with Python-SWIG example

2010-04-14 Thread AhmadKhaled

Hi,

I was trying to run the example given on the Grass and Python wiki for
interface between Python-SWIG-GRASS. I ran the raster example, and I got the
following error: 

TypeError: 'int' object is unsubscriptable

Which kinds of make sense, because presumably the G_get_raster_row function
returns an integer, to show if the reading was successful or not. So we can
not iterate on that to get the row.
I think that the function returns the result in the intrast variable here,
but I can not access it, because it says that it is a Swig Object of type
'void *'. Any idea on how I can read the rows of the map?

Thanks in advance,
Ahmad 
-- 
View this message in context: 
http://n2.nabble.com/Question-with-Python-SWIG-example-tp4904954p4904954.html
Sent from the Grass - Users mailing list archive at Nabble.com.
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user