Re: [GRASS-user] grass env

[maven apache]

Hi:
I found that the grassBatchJob is so excited that I like it,however I wonder
that each time a user who want to run grass in my web application have to
write a .sh file and then call the grass with the BATCH_JOB may cause low
effectivity? after all, this is realated the IO operation, isn't it?
2010/3/3 Hamish 

> ... and if *any* user editable inputs will be visible from the web side of
> the app make 100% sure that you have bounds checked and sanitized every
> single one of them. Stripping all punctuation and limiting the string
> length
> before passing as a module option is a good first step.
>
> I've no idea about java but with unix power tools pipe it through
>  `cut -b 255 | sed -e 's/[^a-zA-Z0-9_]//g'`
>
> to only keep the first 255 chars, and only keep a-z, A-Z, 0-9, and the
> underscore "_".
>
>
> there are chances for buffer overflows and unquoted shell script variables
> all over the place.
>
>
> Hamish
>
>
>
>
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] grass env

[Hamish]

... and if *any* user editable inputs will be visible from the web side of
the app make 100% sure that you have bounds checked and sanitized every
single one of them. Stripping all punctuation and limiting the string length
before passing as a module option is a good first step.

I've no idea about java but with unix power tools pipe it through
  `cut -b 255 | sed -e 's/[^a-zA-Z0-9_]//g'`

to only keep the first 255 chars, and only keep a-z, A-Z, 0-9, and the
underscore "_".


there are chances for buffer overflows and unquoted shell script variables
all over the place.


Hamish



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

Re: [GRASS-user] run grass in linux command -- r.fillnulls errors

[maven apache]

2010/3/3 Markus Neteler 

> On Tue, Mar 2, 2010 at 2:55 PM, maven apache 
> wrote:
> > Hi:
> > When I run grass in command line model, the r.fillnulls can not be
> > completed, I post the input and output in the terminal , I hope some one
> can
> > check it :
> > 
> > GRASS 6.4.0RC5 (spearfish60):~/grass/GrassDataBase > g.gisenv
>
> please post the output of
> g.region -p
>

GRASS 6.4.0RC5 (spearfish60):~ > g.gisenv
GISDBASE=/home/kk/grass/GrassDataBase
LOCATION_NAME=spearfish60
MAPSET=PERMANENT
GRASS_GUI=text

GRASS 6.4.0RC5 (spearfish60):~ > g.region -p
projection: 1 (UTM)
zone:   13
datum:  nad27
ellipsoid:  clark66
north:  4928000
south:  4914000
west:   59
east:   609000
nsres:  100
ewres:  100
rows:   140
cols:   190
cells:  26600

GRASS 6.4.0RC5 (spearfish60):~ > r.in.gdal
input=/home/kk/grass/data/4404.tif output=aa location=al
Location  created
 100%
r.in.gdal complete. Raster map  created.

GRASS 6.4.0RC5 (spearfish60):~ > g.mapset location=al mapset=PERMANENT
Erasing monitors...
Cleaning up temporary files...
WARNING: Your shell continues to use the history for the old mapset
You can switch the history by commands:
history -w; history -r
/home/kk/grass/GrassDataBase/al/PERMANENT/.bash_history;
HISTFILE=/home/kk/grass/GrassDataBase/al/PERMANENT/.bash_history

GRASS 6.4.0RC5 (spearfish60):~ > g.gisenv
GISDBASE=/home/kk/grass/GrassDataBase
LOCATION_NAME=al
MAPSET=PERMANENT
GRASS_GUI=text
MONITOR=cairo7

GRASS 6.4.0RC5 (spearfish60):~ > g.region -p
projection: 3 (Latitude-Longitude)
zone:   0
datum:  wgs84
ellipsoid:  wgs84
north:  90N
south:  90S
west:   180W
east:   180E
nsres:  0:15
ewres:  0:15
rows:   720
cols:   1440
cells:  1036800
GRASS 6.4.0RC5 (spearfish60):~ >




> ...
> > GRASS 6.4.0RC5 (spearfish60):~/grass/GrassDataBase > r.fillnulls input=aa
> > output=abc tension=40 smooth=0.1
> > Locating and isolating NULL areas...
> >  100%
> ...
> > Number of points: 704828
> ...
> > Removing raster 
> > Using segmentation for interpolation...
> > Percent complete:
> ...
> > will require 4874896 bytes of disk space for temp files
> > Killed
>
> I wonder if too much RAM was used and swap exceeded?
>
I run linux in the vitual machine within vmware in the win platform, is it
who cause the problem?

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

[GRASS-user] sistemi di coordinate& proiezioni

[Matteo Toro]

Salve,

ho dei file tiff scannerizzati  ai quali sono associati dei .tfw. questi sono 
del tipo:

 2.11975401429604
 0.00
 0.00
-2.11975401429604
906737.1811447544
   4718928.9772432996

e queste sono le coordinate dell'upper left del foglio igm 141-III (per capirci 
la zona di chieti-abruzzo).

Queste coordinate dovrebbero essere UTM-ED50, in quanto esse dovrebbero essere 
la distanza in metri da greenwich (la 4a riga), e la distanza in metri 
dall'equatore(la 5a riga).

Se questo è giusto, dovrei procedere al fine di  trasformare queste coordinate 
da proiettate (e chilometriche) a geografiche (ed espresse in gradi). Per fare 
questo, dovrei creare una location/mapset in utm-ed50 per importare tutti 
questi .tif. In secondo luogo dovrei creare una target location/mapset con 
epsg:4326, al fine di trasformare le coordinate di queste carte in lat/long 
WGS84. Infine dovrei procedere con l'importazione delle carte dalla location 
originale alla target location: è possibile digitare il comando da terminale 
una sola volta per tutte le carte, oppure anche se non utilizzo il gui, il 
comando nel terminal va ripetuto per ogni singolo .tif?

Siccome penso di essermi perso, vi prego di darmi una mano :D





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

Re: [GRASS-user] Custom Projection Problems

[Paul Kelly]

On Tue, 2 Mar 2010, Casey Vandenberg wrote:


Thanks Paul,

That was a typo on my part when specifying the ellipsoid. By typing in 
clrk66, it does maintain the correct ellipsoid and datum transform 
parameters. It is unnerving though how when I specify the Nad27 datum with an 
associated datum transformation shift, it will reject it and use only the 
nad27 default transformation parameters.


I agree - it could arguably be considered a bug. I've submitted a patch to 
GDAL for consideration, which would fix this - see 
http://trac.osgeo.org/gdal/ticket/3450


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

Re: [GRASS-user] Connecting Vectors to Attribute Tables

[Micha Silver]

Andrew Lewin wrote:


Dear All,

I am having some difficulty viewing an attribute table from a vector 
file that is connected to a sqlite database table.


Here is the terminal snippet of my commands with errors:

GRASS 6.4.0RC5 (SantaMonicaBayMM):~ > v.in.db driver=sqlite 
database=/USers/andrewlewin/MM_SMB.db table=MMSightings1 x=LONG y=LAT 
key=cat out=MMSightings_layer2Writing features...

 100%
Copying attributes...
DBMI-SQLite driver error:
Cannot scan date:"Date"



Is there a column named "Date". Maybe because it's an sqlite reserved 
word, importing the table fails.

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

[GRASS-user] Connecting Vectors to Attribute Tables

[Andrew Lewin]

Dear All,

I am having some difficulty viewing an attribute table from a vector  
file that is connected to a sqlite database table.


Here is the terminal snippet of my commands with errors:

GRASS 6.4.0RC5 (SantaMonicaBayMM):~ > v.in.db driver=sqlite database=/ 
USers/andrewlewin/MM_SMB.db table=MMSightings1 x=LONG y=LAT key=cat  
out=MMSightings_layer2Writing features...

 100%
Copying attributes...
DBMI-SQLite driver error:
Cannot scan date:"Date"

WARNING: Unable to fetch data from table 
WARNING: Unable to copy table
Building topology for vector map ...
Registering primitives...
1897 primitives registered
1897 vertices registered
Building areas...
 100%
0 areas built
0 isles built
Attaching islands...
Attaching centroids...
 100%
Number of nodes: 1803
Number of primitives: 1897
Number of points: 1897
Number of lines: 0
Number of boundaries: 0
Number of centroids: 0
Number of areas: 0
Number of isles: 0
v.in.db complete. 1897 points written to vector map.
GRASS 6.4.0RC5 (SantaMonicaBayMM):~ > v.info MMSightings_layer2
  
+ 
+
 | Layer:
MMSightings_layer2|
 | Mapset:   
PERMANENT |
 | Location: 
SantaMonicaBayMM  |
 | Database:/Users/andrewlewin/ 
grass  |
 |  
Title 
: |
 | Map scale:
1:1   |
 | Map format:   
native|
 | Name of creator:  
andrewlewin   |
 |  
Organization 
:  |
 | Source date: Tue Mar  2 15:03:07  
2010  |
  
| 
|
 |   Type of Map:  vector (level:  
2)  |
  
| 
|
 |   Number of points:   1897Number of areas:   
0  |
 |   Number of lines:0   Number of islands: 
0  |
 |   Number of boundaries:   0   Number of faces:   
0  |
 |   Number of centroids:0   Number of kernels: 
0  |
  
| 
|
 |   Map is 3D:   
No   |
 |   Number of dblinks:   
0|
  
| 
|
 | Projection:  
x,y|
 |   N: 34:58:27.012NS:  
0 |
 |   E: 0W:  
119:18:54.72W |
  
| 
|
 |   Digitization threshold:  
0|
 |
Comments 
:|
  
| 
|
  
+ 
+


GRASS 6.4.0RC5 (SantaMonicaBayMM):~ > v.info -c MMSightings_layer2
ERROR: Database connection for map  is not defined  
in

   DB file
GRASS 6.4.0RC5 (SantaMonicaBayMM):~ > v.db.connect  
map=MMSightings_layer2 table=MMSightings driver=sqlite database=/Users/ 
andrewlewin/MM_SMB.db
The table  is now part of vector map   
and

may be deleted or overwritten by GRASS modules
DBMI-SQLite driver error:
Cannot create index:
create unique index MMSightings_cat on MMSightings ( cat )
table MMSightings has no column named cat

WARNING: Cannot create index
Select privileges were granted on the table
GRASS 6.4.0RC5 (SantaMonicaBayMM):~ > v.db.connect -p  
map=MMSightings_layer2

Vector map  is connected by:
layer <1> table  in database MM_SMB.db> through driver  with key 


There are errors from the point at which I try to make a point vector  
from the table.  Any help would be appreciated.


Andrew

Andrew Lewin, M.Sc., Adv. Dipl. GIS
Email: andrew.le...@sympatico.ca
Tel: (905) 639-4646
Mobile: (905) 630-8441
Skype: andrew.lewin99

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

Re: [GRASS-user] run grass in linux command -- r.fillnulls errors

[Markus Neteler]

On Tue, Mar 2, 2010 at 2:55 PM, maven apache  wrote:
> Hi:
> When I run grass in command line model, the r.fillnulls can not be
> completed, I post the input and output in the terminal , I hope some one can
> check it :
> 
> GRASS 6.4.0RC5 (spearfish60):~/grass/GrassDataBase > g.gisenv

please post the output of
g.region -p

...
> GRASS 6.4.0RC5 (spearfish60):~/grass/GrassDataBase > r.fillnulls input=aa
> output=abc tension=40 smooth=0.1
> Locating and isolating NULL areas...
>  100%
...
> Number of points: 704828
...
> Removing raster 
> Using segmentation for interpolation...
> Percent complete:
...
> will require 4874896 bytes of disk space for temp files
> Killed

I wonder if too much RAM was used and swap exceeded?

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

Re: [GRASS-user] Custom Projection Problems

[Paul Kelly]

On Tue, 2 Mar 2010, Casey Vandenberg wrote:

I am trying to create a new project using a custom proj4 string. The 
parameters I specify are '+proj=utm +zone=6 +ellps=clark66 
+towgs84=-5,135,172'. After creating the project however, and verifying the 
projection settings using *g.proj -j -f -p* I get the following info:


-PROJ_INFO-
name   : Universal Transverse Mercator
proj   : utm
ellps  : wgs84
zone   : 6
towgs84: -5,135,172,0,0,0,0
no_defs: defined
-PROJ_UNITS
unit   : Meter
units  : Meters
meters : 1
+proj=utm +no_defs +zone=6 +a=6378137 +rf=298.257223563 
+towgs84=-5,135,172,0,0,0,0 +to_meter=1


Why is it not maintaining the Clarke 1866 Elliposid?


clark66 is the GRASS abbreviation for that ellipsoid. When you define a 
co-ordinate system using a PROJ.4 string you need to use the PROJ.4 
abbreviation, which is clrk66. So using +ellps=clrk66 should do what you 
want.


Alternatively, if I specify a proj 4 string as follows: '+proj=utm +zone=6 
+datum=nad27 +towgs84=-5,135,172' and verify the created projection using 
g.proj, I get the following output:


-PROJ_INFO-
name   : Universal Transverse Mercator
proj   : utm
datum  : nad27
ellps  : clark66
zone   : 6
no_defs: defined
-PROJ_UNITS
unit   : meter
units  : meters
meters : 1
+proj=utm +no_defs +zone=6 +a=6378206.4 +rf=294.9786982 
+towgs84=-22.000,157.000,176.000 +to_meter=1


Here it is maintaining the correct ellipsoid, but it does not maintain the 
towgs84 datum transformation parameters.


I'm not 100% sure of this (since the string is being interpreted by GDAL 
and not GRASS nor PROJ) but specifying +datum overrides both +ellps and 
+towgs84. You are correctly getting Clarke 66 since it is associated with 
NAD27, but specifying NAD27 seems to be overriding your towgs84 
parameters. I suppose that isn't ideal - but seeing you can work around it 
by not specifying the datum, I guess it's OK.


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

Re: [GRASS-user] Univariate Statistics

[Jarek Jasiewicz]

Luis Lisboa pisze:

Thank you
I realized that there is a r.univar.zone and an r.univar. What is the 
main difference between them? (they seems alike)

Luis

On Tue, Mar 2, 2010 at 9:33 AM, Jose A. Ruiz-Arias > wrote:



You probably should consider R. Take a look at
http://grass.itc.it/statsgrass/grass_geostats.html

Jose

Luis Lisboa escribió:
> Greetings
>
> I need to do some univariate statitstics that goes beyond the
> max/min/median/average/standard deviation. I need to do some
outliers
> detection, distribution analysis and some more basic statistics.
> Since GRASS doesn't perform this I would like to know what are
people
> using? Some statistical Python module (e.g. Statpy) or other
software?
>
> Thanks
>
> Best regards,
> Luis
>

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


--
José A. Ruiz-Arias
Solar Radiation and Atmosphere Modelling Group
http://matras.ujaen.es
Department of Physics
University of Jaén
Campus Lagunillas, Building A3, Dep 066
23071, Jaén (Spain)
Tlf: +34 953 212474
Fax: +34 953 212838






___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
  
r.univar.zonal calculates statistics for more than one area at once. R 
univar only for whole map

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

[GRASS-user] Custom Projection Problems

[Casey Vandenberg]

I am trying to create a new project using a custom proj4 string. The 
parameters I specify are '+proj=utm +zone=6 +ellps=clark66 
+towgs84=-5,135,172'. After creating the project however, and verifying 
the projection settings using *g.proj -j -f -p* I get the following info:


-PROJ_INFO-
name   : Universal Transverse Mercator
proj   : utm
ellps  : wgs84
zone   : 6
towgs84: -5,135,172,0,0,0,0
no_defs: defined
-PROJ_UNITS
unit   : Meter
units  : Meters
meters : 1
+proj=utm +no_defs +zone=6 +a=6378137 +rf=298.257223563 
+towgs84=-5,135,172,0,0,0,0 +to_meter=1


Why is it not maintaining the Clarke 1866 Elliposid?



Alternatively, if I specify a proj 4 string as follows: '+proj=utm 
+zone=6 +datum=nad27 +towgs84=-5,135,172' and verify the created 
projection using g.proj, I get the following output:


-PROJ_INFO-
name   : Universal Transverse Mercator
proj   : utm
datum  : nad27
ellps  : clark66
zone   : 6
no_defs: defined
-PROJ_UNITS
unit   : meter
units  : meters
meters : 1
+proj=utm +no_defs +zone=6 +a=6378206.4 +rf=294.9786982 
+towgs84=-22.000,157.000,176.000 +to_meter=1


Here it is maintaining the correct ellipsoid, but it does not maintain 
the towgs84 datum transformation parameters.




Is this a bug or is my interpretation of creating custom projections 
completely wrong?


Thanks,

Casey

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

Re: [GRASS-user] Univariate Statistics

[Luis Lisboa]

Thank you
I realized that there is a r.univar.zone and an r.univar. What is the main
difference between them? (they seems alike)
Luis

On Tue, Mar 2, 2010 at 9:33 AM, Jose A. Ruiz-Arias wrote:

>
> You probably should consider R. Take a look at
> http://grass.itc.it/statsgrass/grass_geostats.html
>
> Jose
>
> Luis Lisboa escribió:
> > Greetings
> >
> > I need to do some univariate statitstics that goes beyond the
> > max/min/median/average/standard deviation. I need to do some outliers
> > detection, distribution analysis and some more basic statistics.
> > Since GRASS doesn't perform this I would like to know what are people
> > using? Some statistical Python module (e.g. Statpy) or other software?
> >
> > Thanks
> >
> > Best regards,
> > Luis
> > 
> >
> > ___
> > grass-user mailing list
> > grass-user@lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/grass-user
> >
>
>
> --
> José A. Ruiz-Arias
> Solar Radiation and Atmosphere Modelling Group
> http://matras.ujaen.es
> Department of Physics
> University of Jaén
> Campus Lagunillas, Building A3, Dep 066
> 23071, Jaén (Spain)
> Tlf: +34 953 212474
> Fax: +34 953 212838
>
>
>
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] GRASS6.4.0 Ubuntu updates

[Luisa Peña]

Hi

Besides SVN and Linux/GNU generic weekly snapshots, is there any other way
to easily update GRASS for ubuntu? like a special Ubuntu distribvution
 GRASS 6.4.0 Weekly snapshot ?

Thanks

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

Re: [GRASS-user] R-GDAL incompatible libraries

[Kim Besson]

EHEHE this is true inglorious. Me and Gilbert were trying to compile and put
this working and a simple typo was our final error.
Thanks. Now it's working

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

Re: [GRASS-user] R-GDAL incompatible libraries

[Gilbert Ferrara]

EEHEHEHEEH You are absolutely right :)



On Tue, Mar 2, 2010 at 5:58 PM, Leo Hardtke  wrote:

>  Kim
> you have to plot the R object *precip,* not the grass vector *precip30n*
> plot(precip,lwd=2, col="brown") should work!
>
> Leo
>
> Kim Besson escribió:
>
> Following Leo's suggestion I installed GDAL 1.6.3 and I tried again:
>
> > library(spgrass6)
> Loading required package: sp
> Loading required package: rgdal
> Geospatial Data Abstraction Library extensions to R successfully loaded
> Loaded GDAL runtime: GDAL 1.6.3, released 2009/11/19
> Path to GDAL shared files: /usr/share/gdal16
> Loaded PROJ.4 runtime: Rel. 4.6.1, 21 August 2008
> Path to PROJ.4 shared files: (autodetected)
> Loading required package: XML
> GRASS GIS interface loaded with GRASS version: 6.4.0svn
> and location: nc_spm_08
>
> (seems OK?)
>
> > precip <- readVECT6("precip_30ynormals", ignore.stderr)
> OGR data source with driver: GRASS
> Source: "/mnt/GIS/nc_spm_08/PERMANENT/vector/precip_30ynormals/head",
> layer: "1"
> with 136 features and 18 fields
> Feature type: wkbPoint with 3 dimensions
> (seems OK?)
>
> When I did
> > plot(nc_state, axes=TRUE)
> > plot(precip30n, add=TRUE, lwd=2, col="brown")
> Error: object 'precip30n' not found
> Error in plot(precip30n, add = TRUE, lwd = 2, col = "brown") :
>   error in evaluating the argument 'x' in selecting a method for function
> 'plot'
> THe last plot gave that error
>
> But would you say that it seems to be working?
>
> Thank you all for your help and support
> Thanks
> Kim
>
>
> Question: Is there any known problem with GRASS and GDAL1.6.3?
>
> 2010/3/2 Kim Besson 
>
>>
>>  Ok I reinstalled rgdal
>> R CMD INSTALL rgdal
>>
>>  And then:
>> R
>> library(spgrass6)
>>
>>  Output:
>>  Loading required package: sp
>> Loading required package: rgdal
>> Geospatial Data Abstraction Library extensions to R successfully loaded
>> Loaded GDAL runtime: GDAL 1.5.4, released 2009/01/07
>> Path to GDAL shared files: /usr/share/gdal15
>> Loaded PROJ.4 runtime: Rel. 4.6.1, 21 August 2008
>> Path to PROJ.4 shared files: (autodetected)
>> Loading required package: XML
>> GRASS GIS interface loaded with GRASS version: 6.4.0svn
>>  and location: North-Carolina
>>
>>  precip30 <- readVECT6("precip_30ynormals", ignore.stderr=TRUE)
>>  ERROR: Incompatible library version for module. You need to rebuild
>> GRASS
>>or untangle multiple installations.
>>
>>   Same error.
>> I'm running GRASSsvn6.4.0, from latest week, in a Linux Machine with:
>> rgdal_0.6-10
>> spgrass6_0.6-15.
>> R version 2.9.2 (2009-08-24)
>>
>>  Any idea of what might not be compatible?
>>
>>  Thank you
>> Kim
>>
>>  2010/3/2 Markus Neteler 
>>
>>> 2010/3/2 Kim Besson :
>>>
>>> > Ok I have just installed latest GRASS snapshot and I got the following
>>> error
>>> >> precip30n <- readVECT6("precip_30ynormals", ignore.stderr=TRUE)
>>> > ERROR: Incompatible library version for module. You need to rebuild
>>> GRASS
>>> >or untangle multiple installations.
>>> > What I would like to know is the versions of R and spgrass6 that I
>>> should
>>> > install on my system in order to not get this errors?
>>>
>>>  You also need to update the GRASS-GDAL plugin.
>>>
>>> Markus
>>>
>>
>>
> --
>
> ___
> grass-user mailing 
> listgrass-u...@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user
>
>
>
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] R-GDAL incompatible libraries

[Leo Hardtke]

Kim
you have to plot the R object precip, not the grass vector precip30n
plot(precip,lwd=2, col="brown") should work!

Leo

Kim Besson escribió:
Following Leo's suggestion I installed GDAL 1.6.3 and I
tried again:
  
> library(spgrass6)
Loading required package: sp
Loading required package: rgdal
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.6.3, released 2009/11/19
Path to GDAL shared files: /usr/share/gdal16
Loaded PROJ.4 runtime: Rel. 4.6.1, 21 August 2008
Path to PROJ.4 shared files: (autodetected)
Loading required package: XML
GRASS GIS interface loaded with GRASS version: 6.4.0svn
and location: nc_spm_08
  
(seems OK?)
  
> precip <- readVECT6("precip_30ynormals", ignore.stderr)
OGR data source with driver: GRASS 
Source: "/mnt/GIS/nc_spm_08/PERMANENT/vector/precip_30ynormals/head",
layer: "1"
with 136 features and 18 fields
Feature type: wkbPoint with 3 dimensions
(seems OK?)
  
When I did
> plot(nc_state, axes=TRUE)
> plot(precip30n, add=TRUE, lwd=2, col="brown")
Error: object 'precip30n' not found
Error in plot(precip30n, add = TRUE, lwd = 2, col = "brown") : 
  error in evaluating the argument 'x' in selecting a method for
function 'plot'
THe last plot gave that error
  
But would you say that it seems to be working?
  
Thank you all for your help and support
Thanks
Kim
  
  
Question: Is there any known problem with GRASS and GDAL1.6.3?
  
  2010/3/2 Kim Besson 
  


Ok I reinstalled rgdal
R CMD INSTALL rgdal


And then:
R
library(spgrass6)


Output:

Loading required package: sp
Loading required package: rgdal
Geospatial Data Abstraction Library extensions to R
successfully loaded
Loaded GDAL runtime: GDAL 1.5.4, released 2009/01/07
Path to GDAL shared files: /usr/share/gdal15
Loaded PROJ.4 runtime: Rel. 4.6.1, 21 August 2008
Path to PROJ.4 shared files: (autodetected)
Loading required package: XML
GRASS GIS interface loaded with GRASS version: 6.4.0svn

and location: North-Carolina



precip30 <- readVECT6("precip_30ynormals",
ignore.stderr=TRUE)

ERROR: Incompatible library version for module. You need to
rebuild GRASS
       or untangle multiple installations.





Same error. 
I'm running GRASSsvn6.4.0, from latest week, in a Linux
Machine with:
rgdal_0.6-10
spgrass6_0.6-15.
R version 2.9.2 (2009-08-24)


Any idea of what might not be compatible?


Thank you
Kim


2010/3/2 Markus Neteler 

2010/3/2
Kim Besson :
  
  
  > Ok I have just installed latest GRASS snapshot and I
got the following error
>> precip30n <- readVECT6("precip_30ynormals",
ignore.stderr=TRUE)
> ERROR: Incompatible library version for module. You need to
rebuild GRASS
>        or untangle multiple installations.
> What I would like to know is the versions of R and spgrass6 that I
should
> install on my system in order to not get this errors?
  
  
You also need to update the GRASS-GDAL plugin.
  
Markus
  
  




  
  
  
  

___
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] Display of Terrain Perameter Maps

[Dylan Beaudette]

On Tuesday 02 March 2010, Rich Shepard wrote:
> On Tue, 2 Mar 2010, Dylan Beaudette wrote:
> > I have found that with a certain amount of smoothing (r.param.scale ,
> > r.neighbors) and the right color scheme (r.colors [-e]) one can generate
> > pretty nice looking maps that display overall trends.
>
> Dylan,
>
>I used r.param.scale extensively; will learn r.neighbors.
>
> > One thing to remember-- most GRASS modules are geared towards analysis--
> > making pretty maps comes second.
>
>That's very true. However, when the analytical results need to be
> presented to non-technical decision-makers (lawyers, regulators,
> politicians) it's vital that the presentations be visually attractive. It
> adds credibility to the analytical results and, like typeset docuements
> rather than word processed ones, makes them more easily accepted.
>

I agree 100% -- All I was getting at was that the nice map-making requires 
some creativity. If you have some time check out GMT-- it can generate some 
very nice looking maps in PS/EPS format. There are some scripts on the GRASS 
add-ons wiki page with some ideas. 

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

Re: [GRASS-user] R-GDAL incompatible libraries

[Alex Mandel]

Kim Besson wrote:
> Following Leo's suggestion I installed GDAL 1.6.3 and I tried again:
> 
>> library(spgrass6)
> Loading required package: sp
> Loading required package: rgdal
> Geospatial Data Abstraction Library extensions to R successfully loaded
> Loaded GDAL runtime: GDAL 1.6.3, released 2009/11/19
> Path to GDAL shared files: /usr/share/gdal16
> Loaded PROJ.4 runtime: Rel. 4.6.1, 21 August 2008
> Path to PROJ.4 shared files: (autodetected)
> Loading required package: XML
> GRASS GIS interface loaded with GRASS version: 6.4.0svn
> and location: nc_spm_08
> 
> (seems OK?)
> 
>> precip <- readVECT6("precip_30ynormals", ignore.stderr)
> OGR data source with driver: GRASS
> Source: "/mnt/GIS/nc_spm_08/PERMANENT/vector/precip_30ynormals/head", layer:
> "1"
> with 136 features and 18 fields
> Feature type: wkbPoint with 3 dimensions
> (seems OK?)
> 
> When I did
>> plot(nc_state, axes=TRUE)
>> plot(precip30n, add=TRUE, lwd=2, col="brown")
> Error: object 'precip30n' not found
> Error in plot(precip30n, add = TRUE, lwd = 2, col = "brown") :
>   error in evaluating the argument 'x' in selecting a method for function
> 'plot'
> THe last plot gave that error
> 
> But would you say that it seems to be working?
> 
> Thank you all for your help and support
> Thanks
> Kim

Looks like a typo this time you put your value into precip and then
tried to call precip30n.

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

Re: [GRASS-user] Display of Terrain Perameter Maps

[Rich Shepard]

On Tue, 2 Mar 2010, Dylan Beaudette wrote:


I have found that with a certain amount of smoothing (r.param.scale ,
r.neighbors) and the right color scheme (r.colors [-e]) one can generate
pretty nice looking maps that display overall trends.


Dylan,

  I used r.param.scale extensively; will learn r.neighbors.


One thing to remember-- most GRASS modules are geared towards analysis--
making pretty maps comes second.


  That's very true. However, when the analytical results need to be
presented to non-technical decision-makers (lawyers, regulators,
politicians) it's vital that the presentations be visually attractive. It
adds credibility to the analytical results and, like typeset docuements
rather than word processed ones, makes them more easily accepted.

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

Re: [GRASS-user] R-GDAL incompatible libraries

[Kim Besson]

Following Leo's suggestion I installed GDAL 1.6.3 and I tried again:

> library(spgrass6)
Loading required package: sp
Loading required package: rgdal
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.6.3, released 2009/11/19
Path to GDAL shared files: /usr/share/gdal16
Loaded PROJ.4 runtime: Rel. 4.6.1, 21 August 2008
Path to PROJ.4 shared files: (autodetected)
Loading required package: XML
GRASS GIS interface loaded with GRASS version: 6.4.0svn
and location: nc_spm_08

(seems OK?)

> precip <- readVECT6("precip_30ynormals", ignore.stderr)
OGR data source with driver: GRASS
Source: "/mnt/GIS/nc_spm_08/PERMANENT/vector/precip_30ynormals/head", layer:
"1"
with 136 features and 18 fields
Feature type: wkbPoint with 3 dimensions
(seems OK?)

When I did
> plot(nc_state, axes=TRUE)
> plot(precip30n, add=TRUE, lwd=2, col="brown")
Error: object 'precip30n' not found
Error in plot(precip30n, add = TRUE, lwd = 2, col = "brown") :
  error in evaluating the argument 'x' in selecting a method for function
'plot'
THe last plot gave that error

But would you say that it seems to be working?

Thank you all for your help and support
Thanks
Kim


Question: Is there any known problem with GRASS and GDAL1.6.3?

2010/3/2 Kim Besson 

>
> Ok I reinstalled rgdal
> R CMD INSTALL rgdal
>
> And then:
> R
> library(spgrass6)
>
> Output:
> Loading required package: sp
> Loading required package: rgdal
> Geospatial Data Abstraction Library extensions to R successfully loaded
> Loaded GDAL runtime: GDAL 1.5.4, released 2009/01/07
> Path to GDAL shared files: /usr/share/gdal15
> Loaded PROJ.4 runtime: Rel. 4.6.1, 21 August 2008
> Path to PROJ.4 shared files: (autodetected)
> Loading required package: XML
> GRASS GIS interface loaded with GRASS version: 6.4.0svn
> and location: North-Carolina
>
> precip30 <- readVECT6("precip_30ynormals", ignore.stderr=TRUE)
> ERROR: Incompatible library version for module. You need to rebuild GRASS
>or untangle multiple installations.
>
> Same error.
> I'm running GRASSsvn6.4.0, from latest week, in a Linux Machine with:
> rgdal_0.6-10
> spgrass6_0.6-15.
> R version 2.9.2 (2009-08-24)
>
> Any idea of what might not be compatible?
>
> Thank you
> Kim
>
> 2010/3/2 Markus Neteler 
>
>> 2010/3/2 Kim Besson :
>>
>> > Ok I have just installed latest GRASS snapshot and I got the following
>> error
>> >> precip30n <- readVECT6("precip_30ynormals", ignore.stderr=TRUE)
>> > ERROR: Incompatible library version for module. You need to rebuild
>> GRASS
>> >or untangle multiple installations.
>> > What I would like to know is the versions of R and spgrass6 that I
>> should
>> > install on my system in order to not get this errors?
>>
>> You also need to update the GRASS-GDAL plugin.
>>
>> Markus
>>
>
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Display of Terrain Perameter Maps

[Dylan Beaudette]

On Monday 01 March 2010, Rich Shepard wrote:
>The output maps from r.slope.aspect. r.resamp.rst, and even r.topidx the
> output maps have values for each cell rather than trends. This is very
> difficult to interpret by non-technical viewers.
>
>I've seen maps of slopes, profile and planar curvature, and topographic
> wetness indices that are more smoothly colored. If one wants to determine
> the curvature of slopes adjacent to a specific cell at its base it would be
> easier to see this if the values do not change with each cell. When we
> stand at the foot of a hill we tend to not focus on each detail of the
> terrain but the overall slope, profile curvature, and planar curvature.
> It's this type of overview that I'm trying to display.
>
>Is there a way to get smoother output maps that are visually more
> informative to non-specialists?
>
> Rich
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user

Hi Rich,

I have found that with a certain amount of smoothing (r.param.scale , 
r.neighbors) and the right color scheme (r.colors [-e]) one can generate 
pretty nice looking maps that display overall trends. 

One thing to remember-- most GRASS modules are geared towards analysis-- 
making pretty maps comes second.

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

Re: [GRASS-user] Display of Terrain Perameter Maps

[Rich Shepard]

On Tue, 2 Mar 2010, Glynn Clements wrote:


Filter the input data with r.neighbors[1], r.mfilter, etc.

[1] I wouldn't use method=mode if you're going to calculate
derivatives. I'm not so sure about method=median either.


  I appreciate learning of options that did not strike me as appropriate
when I looked through the module list.

Thanks, Glynn,

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

Re: [GRASS-user] grass env

[Glynn Clements]

maven apache wrote:

> I have post this problem before, I want to run grass in command line model
> so I can call them in java at that time,however it often cause some
> unexpeced errors.
> So today I change to Ubuntu, and I install grass use sudo apt-get install
> grass.
> 
> Now I also want to call grass in my java code but in the Ubuntu platform, I
> want to know which environment variables are needed, so I  did some tests.
> First I run grass in text model in the ternimal by " grass -text",
> and I chose the location and mapset,it works, then I use command "set >
> /home/kk/var" to get all the environment variablse using by grass ,then I
> open a new ternimal to get the normal env vars by " set > /home/kk/varn" .
> Then  I compare the two files and tried to find the differences. I found
> that there are more variables in the varn rather than var (shown below):
> ---
> GISBASE=/usr/lib/grass64
> GISRC=/tmp/grass6-kk-9581/gisrc
> GIS_LOCK=9581
> GRASS_GNUPLOT='gnuplot -persist'
> GRASS_HTML_BROWSER=x-www-browser
> GRASS_LD_LIBRARY_PATH=/usr/lib/grass64/lib
> GRASS_PAGER=pager
> GRASS_PERL=/usr/bin/perl
> GRASS_PROJSHARE=/usr/share/proj
> GRASS_PYTHON=python
> GRASS_SH=/bin/sh
> GRASS_TCLSH=tclsh
> GRASS_VERSION=6.4.0RC5
> GRASS_WISH=wish
> ---
> So I want to know if these vars are enough?

You also need:

PATH=$GISBASE/bin:$GISBASE/bin:$PATH
LD_LIBRARY_PATH=$GISBASE/lib
PYTHONPATH=$GISBASE/etc/python

[PYTHONPATH is only relevant if running Python scripts which use the
grass.script modules.]

Many of the above are unnecessary.

GRASS_TCLSH and GRASS_PERL aren't used. GRASS_SH is only used on
Windows. GRASS_PROJSHARE defaults to /usr/share/proj if not set. 
GRASS_WISH, GRASS_PYTHON, GRASS_HTML_BROWSER and GRASS_LD_LIBRARY_PATH
are only relevant for interactive use. GRASS_GNUPLOT is only used by
i.spectral.

Although it shouldn't matter, I suggest GRASS_PAGER=cat and
GRASS_MESSAGE_FORMAT=silent.

If you need to use any display commands, you also need to set
GRASS_RENDER_IMMEDIATE=TRUE, and may want to use some of the variables
documented in the "pngdriver" manual page to control the output.

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

Re: [GRASS-user] Display of Terrain Perameter Maps

[Glynn Clements]

Rich Shepard wrote:

>The output maps from r.slope.aspect. r.resamp.rst, and even r.topidx the
> output maps have values for each cell rather than trends. This is very
> difficult to interpret by non-technical viewers.
> 
>I've seen maps of slopes, profile and planar curvature, and topographic
> wetness indices that are more smoothly colored. If one wants to determine
> the curvature of slopes adjacent to a specific cell at its base it would be
> easier to see this if the values do not change with each cell. When we stand
> at the foot of a hill we tend to not focus on each detail of the terrain but
> the overall slope, profile curvature, and planar curvature. It's this type
> of overview that I'm trying to display.
> 
>Is there a way to get smoother output maps that are visually more
> informative to non-specialists?

Filter the input data with r.neighbors[1], r.mfilter, etc.

[1] I wouldn't use method=mode if you're going to calculate
derivatives. I'm not so sure about method=median either.

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

Re: [GRASS-user] R-GDAL incompatible libraries

[Leo Hardtke]

Kim
I think you have to update gdal (1.6.3-2 in ubuntugis-unstable repo)
and maybe use grass 6.4rc5 (ubuntugis)

leo

> library(spgrass6)
Loading required package: sp
Loading required package: rgdal
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.6.3, released
2009/11/19
Path to GDAL shared files: /usr/share/gdal16
Loaded PROJ.4 runtime: Rel. 4.7.1, 23 September 2009
Path to PROJ.4 shared files: (autodetected)
Loading required package: XML
GRASS GIS interface loaded with GRASS version: 6.5.svn
and location: hardtke

Kim Besson escribió:

  
  
  Ok I reinstalled rgdal
  R CMD INSTALL rgdal
  
  
  And then:
  R
  library(spgrass6)
  
  
  Output:
  
  Loading required package: sp
  Loading required package: rgdal
  Geospatial Data Abstraction Library extensions to R successfully
loaded
  Loaded GDAL runtime: GDAL 1.5.4,
released 2009/01/07
  Path to GDAL shared files: /usr/share/gdal15
  Loaded PROJ.4 runtime: Rel. 4.6.1, 21 August 2008
  Path to PROJ.4 shared files: (autodetected)
  Loading required package: XML
  GRASS GIS interface loaded with GRASS version: 6.4.0svn
  and location: North-Carolina
  
  
  
  precip30 <- readVECT6("precip_30ynormals", ignore.stderr=TRUE)
  ERROR: Incompatible library version for module. You need to
rebuild GRASS
         or untangle multiple installations.
  
  
  
  
Same error. 
  I'm running GRASSsvn6.4.0, from latest week, in a Linux Machine
with:
  rgdal_0.6-10
  spgrass6_0.6-15.
  R version 2.9.2 (2009-08-24)
  
  
  Any idea of what might not be compatible?
  
  
  Thank you
  Kim
  
  2010/3/2 Markus Neteler 
  2010/3/2
Kim Besson :
> Ok I have just installed latest GRASS snapshot
and I got the following error
>> precip30n <- readVECT6("precip_30ynormals",
ignore.stderr=TRUE)
> ERROR: Incompatible library version for module. You need to
rebuild GRASS
>        or untangle multiple installations.
> What I would like to know is the versions of R and spgrass6 that I
should
> install on my system in order to not get this errors?


You also need to update the GRASS-GDAL plugin.

Markus

  
  
  
  

___
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] R-GDAL incompatible libraries

[Kim Besson]

Ok I reinstalled rgdal
R CMD INSTALL rgdal

And then:
R
library(spgrass6)

Output:
Loading required package: sp
Loading required package: rgdal
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.5.4, released 2009/01/07
Path to GDAL shared files: /usr/share/gdal15
Loaded PROJ.4 runtime: Rel. 4.6.1, 21 August 2008
Path to PROJ.4 shared files: (autodetected)
Loading required package: XML
GRASS GIS interface loaded with GRASS version: 6.4.0svn
and location: North-Carolina

precip30 <- readVECT6("precip_30ynormals", ignore.stderr=TRUE)
ERROR: Incompatible library version for module. You need to rebuild GRASS
   or untangle multiple installations.

Same error.
I'm running GRASSsvn6.4.0, from latest week, in a Linux Machine with:
rgdal_0.6-10
spgrass6_0.6-15.
R version 2.9.2 (2009-08-24)

Any idea of what might not be compatible?

Thank you
Kim

2010/3/2 Markus Neteler 

> 2010/3/2 Kim Besson :
> > Ok I have just installed latest GRASS snapshot and I got the following
> error
> >> precip30n <- readVECT6("precip_30ynormals", ignore.stderr=TRUE)
> > ERROR: Incompatible library version for module. You need to rebuild GRASS
> >or untangle multiple installations.
> > What I would like to know is the versions of R and spgrass6 that I should
> > install on my system in order to not get this errors?
>
> You also need to update the GRASS-GDAL plugin.
>
> Markus
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] R-GDAL incompatible libraries

[Leo Hardtke]

Hi Kim
I am using R 9.2 (ubuntu repo),
grass6.5 (source), and grass6.4(ubuntugis) works too,
rgdal (source, not the lates 
thttp://cran.r-project.org/src/contrib/Archive/rgdal/rgdal_0.6-21.tar.gz 
), and

spgrass6 (install.packages("spgras6") in R)
Hope it works

Leo Hardtke


Kim Besson escribió:
Ok I have just installed latest GRASS snapshot and I got the following 
error

> precip30n <- readVECT6("precip_30ynormals", ignore.stderr=TRUE)
ERROR: Incompatible library version for module. You need to rebuild GRASS
   or untangle multiple installations.

What I would like to know is the versions of R and spgrass6 that I 
should install on my system in order to not get this errors?


Thanks
Kim

2010/3/2 Kim Besson >


Hi 


I have installed R from Ubuntu repositories. My grass installation
is a Weekly snapshot. 
I had to install rgdal since it doesn't exist in UBUNTU

repositories neither CRAN and the same with spgrass6.

Ok I will rebuild grass and reply again to mailing list

Thanks
Kim

Did you install grass from source? Did you install R, spgrass6 from

source?

- If you have installed grass from source, then try to rebuild
grass
(after having installed R and the rest).

- If you have grass from the ubuntu repositories, the you should
probably stick with R and the rest from the repositories as
well. If the
latter is the case and it does not work, then there is a problem.

Tell us a bit more...

Nikos





___
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] Raster attributes

[Ned Horning]

Hamish,

Thanks for the tips. Doing the raster to vector and then vector to 
raster conversion probably won't work so well for this application. The 
image is quite large and I don't want to risk modifying the data in the 
conversion process.


I'll see what I can do with regard to writing a script. I have a decent 
understanding of r.reclass and r.category but the problem is that I'm 
not sure how to join the database to the raster. I'll see if the 
r.rast4d module can provide some hints.


Ned

On 3/2/2010 2:38 AM, Hamish wrote:

Ned Horning wrote:
   

Is it possible to link a raster image
to a set of attributes stored in a database using GRASS?
 

categorical (CELL) raster maps have very basic support via the r.category
module; the category value of the cell can hold the key column ID.

for something more advanced see:
  
http://grass.osgeo.org/wiki/GRASS_SoC_Ideas_2008#Time_series_management_in_GRASS
  http://grass.osgeo.org/wiki/Time_series_in_GRASS
  http://trac.osgeo.org/grass/browser/grass-addons/raster/r.rast4d


   

I am working with the Harmonized World Soil Database which
consist of an image and a database that holds a set of soil
attributes. These attributes can be linked/joined to the
image since the pixel values are IDs that appear in the
database. I would like to output a series of images, each
representing a different attribute such as soil unit,
texture, soil depth. Is this possible in GRASS?
 

I think it would be quite possible to write a nice script that does this with 
r.reclass + r.category. It would be a nice addition to wiki addons :)

I think it could be done without the complication of an extra SQL layer,
but whatever works.


probably the easiest idea though is to simply run 'r.to.vect -v
feature=area' and then v.db.connect. You will be limited by the number
of cells though, probably if the raster is bigger than 1600x1600 you
might run into memory problems. (depends on how many repeated cell values/
clumps are in the map)


good luck,
Hamish





   


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

Re: [GRASS-user] Terrain Curve Values

[Rich Shepard]

On Tue, 2 Mar 2010, Markus Neteler wrote:


Here the answer from Helena (not on this list):



yes - convex is positive and concave is negative. It is the same as in 
v.surf.rst.


  Thank you both.

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

[GRASS-user] run grass in linux command -- r.fillnulls errors

[maven apache]

Hi:
When I run grass in command line model, the r.fillnulls can not be
completed, I post the input and output in the terminal , I hope some one can
check it :

GRASS 6.4.0RC5 (spearfish60):~/grass/GrassDataBase > g.gisenv
GISDBASE=/home/kk/grass/GrassDataBase
LOCATION_NAME=spearfish60
MAPSET=PERMANENT
GRASS_GUI=text
MONITOR=cairo7

GRASS 6.4.0RC5 (spearfish60):~/grass/GrassDataBase > r.in.gdal
input=/home/kk/grass/data/4404.tif output=aa location=al
Location  created
 100%
r.in.gdal complete. Raster map  created.

GRASS 6.4.0RC5 (spearfish60):~/grass/GrassDataBase > g.mapset location=al
mapset=PERMANENT
Erasing monitors...
Cleaning up temporary files...
WARNING: Your shell continues to use the history for the old mapset
You can switch the history by commands:
history -w; history -r
/home/kk/grass/GrassDataBase/al/PERMANENT/.bash_history;
HISTFILE=/home/kk/grass/GrassDataBase/al/PERMANENT/.bash_history

GRASS 6.4.0RC5 (spearfish60):~/grass/GrassDataBase > g.gisenv
GISDBASE=/home/kk/grass/GrassDataBase
LOCATION_NAME=al
MAPSET=PERMANENT
GRASS_GUI=text
MONITOR=cairo7

GRASS 6.4.0RC5 (spearfish60):~/grass/GrassDataBase > r.null map=aa
setnull=-1.267e+30--1.268e+30
 100%

GRASS 6.4.0RC5 (spearfish60):~/grass/GrassDataBase > r.fillnulls input=aa
output=abc tension=40 smooth=0.1
Locating and isolating NULL areas...
 100%
Reading input raster map ...
 100%
Finding buffer zones...
 100%
Writing output raster map ...
 100%
 100%
Creating interpolation points...
WARNING: Default driver / database set to:
 driver: dbf
 database: $GISDBASE/$LOCATION_NAME/$MAPSET/dbf/
Extracting points...
 100%
Building topology for vector map ...
Registering primitives...
704828 primitives registered
704828 vertices registered
Building areas...
 100%
0 areas built
0 isles built
Attaching islands...
Attaching centroids...
 100%
Number of nodes: 704828
Number of primitives: 704828
Number of points: 704828
Number of lines: 0
Number of boundaries: 0
Number of centroids: 0
Number of areas: 0
Number of isles: 0
r.to.vect complete.
Interpolating 704828 points
Note: The following warnings may be ignored.
Removing raster 
Using segmentation for interpolation...
Percent complete:
Loading data from attribute table ...
Reading lines from vector map ...
Reading nodes from vector map ...
 100%
 100%
WARNING: strip exists with insufficient data


The number of points from vector map is 704828
The number of points outside of 2D/3D region 0
The number of points being used is 704828
Processing all selected output files
will require 4874896 bytes of disk space for temp files
Percent complete:
  76%
WARNING: taking too long to find points for interpolation--please change
 the warnnings .
  78%
WARNING: taking too long to find points for interpolation--please change
WARNING: taking too long to find points for interpolation--please change
  80%
 the warnnings ..
  94%
Killed
Note: Above warnings may be ignored.
Patching fill data into NULL areas...
WARNING: Raster map  not found
ERROR: One or more input raster maps not found
Removing raster 
Removing raster 
Removing raster 
WARNING: Raster map  not found
WARNING:  nothing removed
Removing vector 
Filled raster map is: abc
ERROR: Raster map  not found in current mapset
Done.
--
It seems that the r.fillnulls does not work, and when I enter the location
al and mapset PERMANENT , there is just the rast map aa, not abc. I woder
why?
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] grass env

[maven apache]

It seems that the batch job is easy,however I still have some questions:
1) during my process the location and mapset maybe changed, can it be set in
the script file?
2) My application is to support multiple users, so each user should have a
gisrc file, can it be set before run the batch job?

2010/3/2 Markus Neteler 

> On Tue, Mar 2, 2010 at 12:05 PM, maven apache 
> wrote:
> > Hi:
> > I have post this problem before, I want to run grass in command line
> model
> > so I can call them in java at that time,however it often cause some
> > unexpeced errors.
> > So today I change to Ubuntu, and I install grass use sudo apt-get install
> > grass.
> >
> > Now I also want to call grass in my java code but in the Ubuntu platform,
> I
> > want to know which environment variables are needed, so I  did some
> tests.
> > First I run grass in text model in the ternimal by " grass -text",
> > and I chose the location and mapset,it works, then I use command "set >
> > /home/kk/var" to get all the environment variablse using by grass ,then I
> > open a new ternimal to get the normal env vars by " set > /home/kk/varn"
> .
> > Then  I compare the two files and tried to find the differences. I found
> > that there are more variables in the varn rather than var (shown below):
> > ---
> > GISBASE=/usr/lib/grass64
> > GISRC=/tmp/grass6-kk-9581/gisrc
> > GIS_LOCK=9581
> > GRASS_GNUPLOT='gnuplot -persist'
> > GRASS_HTML_BROWSER=x-www-browser
> > GRASS_LD_LIBRARY_PATH=/usr/lib/grass64/lib
> > GRASS_PAGER=pager
> > GRASS_PERL=/usr/bin/perl
> > GRASS_PROJSHARE=/usr/share/proj
> > GRASS_PYTHON=python
> > GRASS_SH=/bin/sh
> > GRASS_TCLSH=tclsh
> > GRASS_VERSION=6.4.0RC5
> > GRASS_WISH=wish
> > ---
> > So I want to know if these vars are enough?
>
> I have no idea. But it is way easier to use GRASS in batch mode:
>
> http://grass.osgeo.org/wiki/GRASS_and_Shell#GRASS_Batch_jobs
>
> Markus
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] R-GDAL incompatible libraries

[Markus Neteler]

2010/3/2 Kim Besson :
> Ok I have just installed latest GRASS snapshot and I got the following error
>> precip30n <- readVECT6("precip_30ynormals", ignore.stderr=TRUE)
> ERROR: Incompatible library version for module. You need to rebuild GRASS
>        or untangle multiple installations.
> What I would like to know is the versions of R and spgrass6 that I should
> install on my system in order to not get this errors?

You also need to update the GRASS-GDAL plugin.

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

Re: [GRASS-user] grass env

[Markus Neteler]

On Tue, Mar 2, 2010 at 12:05 PM, maven apache  wrote:
> Hi:
> I have post this problem before, I want to run grass in command line model
> so I can call them in java at that time,however it often cause some
> unexpeced errors.
> So today I change to Ubuntu, and I install grass use sudo apt-get install
> grass.
>
> Now I also want to call grass in my java code but in the Ubuntu platform, I
> want to know which environment variables are needed, so I  did some tests.
> First I run grass in text model in the ternimal by " grass -text",
> and I chose the location and mapset,it works, then I use command "set >
> /home/kk/var" to get all the environment variablse using by grass ,then I
> open a new ternimal to get the normal env vars by " set > /home/kk/varn" .
> Then  I compare the two files and tried to find the differences. I found
> that there are more variables in the varn rather than var (shown below):
> ---
> GISBASE=/usr/lib/grass64
> GISRC=/tmp/grass6-kk-9581/gisrc
> GIS_LOCK=9581
> GRASS_GNUPLOT='gnuplot -persist'
> GRASS_HTML_BROWSER=x-www-browser
> GRASS_LD_LIBRARY_PATH=/usr/lib/grass64/lib
> GRASS_PAGER=pager
> GRASS_PERL=/usr/bin/perl
> GRASS_PROJSHARE=/usr/share/proj
> GRASS_PYTHON=python
> GRASS_SH=/bin/sh
> GRASS_TCLSH=tclsh
> GRASS_VERSION=6.4.0RC5
> GRASS_WISH=wish
> ---
> So I want to know if these vars are enough?

I have no idea. But it is way easier to use GRASS in batch mode:

http://grass.osgeo.org/wiki/GRASS_and_Shell#GRASS_Batch_jobs

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

Re: [GRASS-user] R-GDAL incompatible libraries

[Νίκος Αλεξανδρής]

On Tue, 2010-03-02 at 12:05 +, Kim Besson wrote:
> Ok I have just installed latest GRASS snapshot

Latest snapshot of grass65 or grass70?

>  and I got the following error
> > precip30n <- readVECT6("precip_30ynormals", ignore.stderr=TRUE)
> ERROR: Incompatible library version for module. You need to rebuild
> GRASS
>or untangle multiple installations.
> 
> 
> What I would like to know is the versions of R and spgrass6 that I
> should install on my system in order to not get this errors?

Check the last posts in this thread:

-http://lists.osgeo.org/pipermail/grass-stats/2009-January/000899.html
and afterwards,

and perhaps:
-http://n2.nabble.com/error-Incompatible-library-version-for-module-You-need-to-rebuild-GRASS-or-untangle-multiple-install-td3347115.html

Maybe there is (still?) a problem with grass70... ? Not sure (since I
haven't been messing around for a while).

Nikos


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

Re: [GRASS-user] problem with spgrass6

[Leo Hardtke]

Hi, thanks for the fast replay

//
Daniel McInerney escribió:

Hi Leo,
these type of questions are best sent to the
grass-st...@list.osgeo.org [1] list (cc'd here),
which discusses GRASS-GIS, R & spgrass6 issues.

What does sessionInfo() give?


sesioninfo() gives
/R version 2.9.2 (2009-08-24)
x86_64-pc-linux-gnu

locale:
LC_CTYPE=es_AR.UTF-8;LC_NUMERIC=C;LC_TIME=es_AR.UTF-8;LC_COLLATE=es_AR.UTF-8;LC_MONETARY=es_AR.UTF-8;LC_MESSAGES=es_AR.UTF-8;LC_PAPER=es_AR.UTF-8;LC_NAME=es_AR.UTF-8;LC_ADDRESS=es_AR.UTF-8;LC_TELEPHONE=es_AR.UTF-8;LC_MEASUREMENT=es_AR.UTF-8;LC_IDENTIFICATION=es_AR.UTF-8

attached base packages:
[1] stats graphics  grDevices utils datasets  methods   base


other attached packages:
[1] tree_1.0-27 spgrass6_0.6-15 XML_2.6-0   rgdal_0.6-21  
[5] sp_0.9-59   rkward_0.5.1  


loaded via a namespace (and not attached):
[1] grid_2.9.2  lattice_0.17-25 tools_2.9.2   /

to Hamish
`/media/COMP/grassdb//hardtke/Bisigato/.tmp/rs7/ndvi_sum'

the "//" in the middle of the path makes me suspect that there is 
a$VARIABLE there which is supposed to be set but is empty for some reason.


   ?

   (does that file exist?)
 


The file dosn't exist...
But when I import another raster (that works), no temporary file is 
created at this location... so I don't understand.


x=readRAST6(c("distancia_costa"))
/media/COMP/grassdb//hardtke/Bisigato/cellhd/distancia_costa has GDAL 
driver GRASS

and has 956 rows and 786 columns

> ls /media/COMP/grassdb//hardtke/Bisigato/.tmp/rs7/
> 29.0  29.1  29.2  29.3  29.4  57.0  57.1  57.2  57.3  671.0  68.0

When I work in grass6.4 , there is no problem and everyting is imported.

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

Re: [GRASS-user] R-GDAL incompatible libraries

[Kim Besson]

Ok I have just installed latest GRASS snapshot and I got the following error
> precip30n <- readVECT6("precip_30ynormals", ignore.stderr=TRUE)
ERROR: Incompatible library version for module. You need to rebuild GRASS
   or untangle multiple installations.

What I would like to know is the versions of R and spgrass6 that I should
install on my system in order to not get this errors?

Thanks
Kim

2010/3/2 Kim Besson 

> Hi
>
> I have installed R from Ubuntu repositories. My grass installation is a
> Weekly snapshot.
> I had to install rgdal since it doesn't exist in UBUNTU repositories
> neither CRAN and the same with spgrass6.
>
> Ok I will rebuild grass and reply again to mailing list
>
> Thanks
> Kim
>
> Did you install grass from source? Did you install R, spgrass6 from
>
>> source?
>>
>> - If you have installed grass from source, then try to rebuild grass
>> (after having installed R and the rest).
>>
>> - If you have grass from the ubuntu repositories, the you should
>> probably stick with R and the rest from the repositories as well. If the
>> latter is the case and it does not work, then there is a problem.
>>
>> Tell us a bit more...
>>
>> Nikos
>>
>>
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Creating mosaics of images

[Νίκος Αλεξανδρής]

Gilbert Ferrara wrote: 
> Is there any function or application to create image mosaics in GRASS?

gdalwarp outside grass, r.patch inside grass. Check the archive [1] (and
maybe the wiki [2]) for examples (like [3]).

Regards, Nikos

---
[1] http://n2.nabble.com/forum/Search.jtp?node=1837860&query=mosaic
[2] http://grass.osgeo.org/wiki/Main_Page
[3] http://grass.osgeo.org/wiki/ASTER_topography


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

Re: [GRASS-user] Creating mosaics of images

[Alberto Pettazzi]

Hi,

check with d.frame command


Gilbert Ferrara escribió:

Greetings

Is there any function or application to create image mosaics in GRASS?

Thanks

Gilber


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



--

Alberto Pettazzi


MeteoGalicia - Departamento de Climatología y Observación

Consellería de Medio Ambiente, Territorio e Infraestruturas

Rúa de Roma, 6

15707 Santiago de Compostela. A Coruña


Teléfono: +34-881-999646


e-mail: alberto.petta...@meteogalicia.es 




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

[GRASS-user] Re: [r.resample] null values for min and max range of data values

[baobazz]

Hi,

Thanks for your answer!

I didn't launch r.in.gdal after r.resample. It's just part of the raster
history info given by the r.info command.

Here is the complete workflow without the info commands. The goal is to
produce an image with a less precise resolution.
1)r.in.gdal input=C:/tests/GRASS/rasterToVect/BDP_08_1234_6134_010_L93.tif
output=raster6133_2 
2)g.region res=2
3)r.resample input=raster6133_2 output=raster.resampled

-raster6133_2 with no null range of values
-raster.resampled with null range of values (:-(()

If you see why there are null range of values..

Thanks!



-- 
View this message in context: 
http://n2.nabble.com/r-resample-null-values-for-min-and-max-range-of-data-values-tp4631171p4660500.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] Creating mosaics of images

[Gilbert Ferrara]

Greetings

Is there any function or application to create image mosaics in GRASS?

Thanks

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

[GRASS-user] grass env

[maven apache]

Hi:
I have post this problem before, I want to run grass in command line model
so I can call them in java at that time,however it often cause some
unexpeced errors.
So today I change to Ubuntu, and I install grass use sudo apt-get install
grass.

Now I also want to call grass in my java code but in the Ubuntu platform, I
want to know which environment variables are needed, so I  did some tests.
First I run grass in text model in the ternimal by " grass -text",
and I chose the location and mapset,it works, then I use command "set >
/home/kk/var" to get all the environment variablse using by grass ,then I
open a new ternimal to get the normal env vars by " set > /home/kk/varn" .
Then  I compare the two files and tried to find the differences. I found
that there are more variables in the varn rather than var (shown below):
---
GISBASE=/usr/lib/grass64
GISRC=/tmp/grass6-kk-9581/gisrc
GIS_LOCK=9581
GRASS_GNUPLOT='gnuplot -persist'
GRASS_HTML_BROWSER=x-www-browser
GRASS_LD_LIBRARY_PATH=/usr/lib/grass64/lib
GRASS_PAGER=pager
GRASS_PERL=/usr/bin/perl
GRASS_PROJSHARE=/usr/share/proj
GRASS_PYTHON=python
GRASS_SH=/bin/sh
GRASS_TCLSH=tclsh
GRASS_VERSION=6.4.0RC5
GRASS_WISH=wish
---
So I want to know if these vars are enough?
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] R-GDAL incompatible libraries

[Kim Besson]

Hi

I have installed R from Ubuntu repositories. My grass installation is a
Weekly snapshot.
I had to install rgdal since it doesn't exist in UBUNTU repositories neither
CRAN and the same with spgrass6.

Ok I will rebuild grass and reply again to mailing list

Thanks
Kim

Did you install grass from source? Did you install R, spgrass6 from

> source?
>
> - If you have installed grass from source, then try to rebuild grass
> (after having installed R and the rest).
>
> - If you have grass from the ubuntu repositories, the you should
> probably stick with R and the rest from the repositories as well. If the
> latter is the case and it does not work, then there is a problem.
>
> Tell us a bit more...
>
> Nikos
>
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Apply a segmentation to derive objects

[Nikos Dumakis]

Hi

Regarsding segmentation, has anyone implemented a method to segment an image
in a unsupervised way? (i.smap is supervised).

Or anyone developed some method to do image segmentation?


Thanks
Nikos

On Mon, Jan 18, 2010 at 11:22 PM, Markus Neteler  wrote:

> On Mon, Jan 18, 2010 at 6:30 PM, Nikos Dumakis 
> wrote:
> > Greetings all
> > I've been following GRASS-user discussion related with
> pos-classigfication,
> > classification, aggregation and so on. But, while I was readingh those
> > emails I figured that I will need to perform, before classification, an
> > image segmentation to derive objects.
>
> You don't need to perform a segmentation before doing classification but
> i.smap is doing it in a combined way.
>
> > Can anyone give me a clue abvout that because in GRASS book, I only
> founded
> > Image segmentation related with processing large quantities of vectorial
> > data (DEM's).
>
> (DEMs are raster data)
>
> What you need to do is to prepare the statistics with i.gensigset. See
> "radiometric & geometric  supervised" in
>
> http://grass.osgeo.org/wiki/Image_classification
>
> (The book is explaining that as well)
>
> Hope this helps,
> Markus
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] POW() in SQLite?

[Hamish]

> Hamish wrote:
> > I am trying to use POW(x,n) with SQLite but it fails:
...
> > no such function: POW

Glynn:
> Nope; not in SQLite:
> http://www.sqlite.org/lang_corefunc.html
> 
..
> > Similar problem if I try x^n:
> > unrecognized token: "^"
> 
> SQLite doesn't have a power operator:
> http://www.sqlite.org/lang_expr.html#binaryops


 
> > this random webhit suggests that it should be possible:
> >   http://osdir.com/ml/sqlite-users/2009-03/msg00248.html
> 
> That says:
..
> Note: "in addition to the SQLite default". IOW, those
> functions are added by SQLiteSpy, and are not part of SQLite.

ah, skimming error.


Rich:
> I'd manipulate those values first, then pass them to SQLite.


solution:

#no good:
#v.db.update basemap column=volume value="POW(length,3)" --verbose

#good:
v.db.select basemap column=cat,length | awk -F'|' \
  '{printf("UPDATE basemap SET volume=%f WHERE cat=%d;\n", $2^3, $1)}' \
  | db.execute



thanks,
Hamish



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

Re: [GRASS-user] Features Extraction

[Hamish]

ab der wrote:
> i'd hope to post those scripts but i don't know how or
> where to post them

see
http://grass.osgeo.org/wiki/GRASS_AddOns
https://trac.osgeo.org/grass/wiki/HowToContribute#WriteaccesstotheGRASS-Addons-SVNrepository



regards,
Hamish



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

Re: [GRASS-user] Export/define region with Map Scale

[Franz Schiller]

Greetings Hamish



> yes. the ps.map module knows its output plots will be on paper,
> and it knows how big that piece of paper is. therefore it can
> correctly calculate and use the scale for PostScript (and so
> PDF format too) hardcopies.
>
>
Ok. So I just need to define the scale-output . Seems reasonably and easy :)


> > 2- Is it possible to set region using this map scale
> > definition in order to have by default this setting?
>
> not directly, but you could write a script wrapping g.region
> and d.info which after calibration could keep you somewhat
> close to your desired scale. That's probably more trouble than
> it's worth though.


yes. It seems to. It seems easier to define Scale only in the output

Thanks

Franz


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

Re: [GRASS-user] ps.map issue

[Hamish]

Rich wrote:
>   This is for my information
> only now as I don't need an answer right away.
> 
>   I have a ps.map script that uses a raster map, two
> boundary vline maps,
> and a vpoints map. There are also clauses for text,
> colortable, and vlegend.
> It make for a pretty map.
> 
>   However, if I add one more vline map (the stream
> channels), the colortable
> does not display and ghostscript complains about a problem.
> I cannot figure
> out why. I changed the sequence in which elements were
> processed and drawn, but this made no difference.
> 
>   Has anyone seen anything like this before?

not that I know of.
does it happen if the vlines map is the only instruction given?
ie can you try to isolate it and submit a test case to the bug trac'er?


thanks
Hamish




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

Re: [GRASS-user] Univariate Statistics

[Jose A. Ruiz-Arias]

You probably should consider R. Take a look at
http://grass.itc.it/statsgrass/grass_geostats.html

Jose

Luis Lisboa escribió:
> Greetings
>
> I need to do some univariate statitstics that goes beyond the
> max/min/median/average/standard deviation. I need to do some outliers
> detection, distribution analysis and some more basic statistics.
> Since GRASS doesn't perform this I would like to know what are people
> using? Some statistical Python module (e.g. Statpy) or other software?
>
> Thanks
>
> Best regards,
> Luis
> 
>
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>   


-- 
José A. Ruiz-Arias
Solar Radiation and Atmosphere Modelling Group
http://matras.ujaen.es
Department of Physics
University of Jaén
Campus Lagunillas, Building A3, Dep 066
23071, Jaén (Spain)
Tlf: +34 953 212474
Fax: +34 953 212838

 

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

Re: [GRASS-user] Univariate Statistics

[Jarosław Jasiewicz]

Luis Lisboa pisze:

Greetings

I need to do some univariate statitstics that goes beyond the 
max/min/median/average/standard deviation. I need to do some outliers 
detection, distribution analysis and some more basic statistics.
Since GRASS doesn't perform this I would like to know what are people 
using? Some statistical Python module (e.g. Statpy) or other software?


Thanks

Best regards,
Luis


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

The easiest way is grass + R + spgrass6,
Jarek
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Univariate Statistics

[Luis Lisboa]

Greetings

I need to do some univariate statitstics that goes beyond the
max/min/median/average/standard deviation. I need to do some outliers
detection, distribution analysis and some more basic statistics.
Since GRASS doesn't perform this I would like to know what are people using?
Some statistical Python module (e.g. Statpy) or other software?

Thanks

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

Re: [GRASS-user] Terrain Curve Values

[Hamish]

Rich wrote:
>   Just to confirm: the curve
> output from r.slope.aspect is positive when the
> curve is convex and negative when the curve is concave.
> Correct?
> 
>   There's nothing on the man page about this. Perhaps
> a sentence could be
> added.


please run an experiment after creating a pseudo surface with
r.surf.volcano (wiki addons) using the gaussian method, and report
back to us. Then we can add the confirmed answer to the help page.

suggestions for pseudo DEM:  (after r.sun wiki page)
  r.surf.volcano out=gauss method=gaussian kurtosis=1


regards,
Hamish



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

Re: [GRASS-user] Display of Terrain Perameter Maps

[Hamish]

Rich:
>   Is there a way to get smoother output maps

run eg 'r.neighbors method=mode' to smooth/denoise maps

'r.param.scale param=feature' also offers a simplified version.

both modules allow you to change the search window size. the bigger the
window, the more smoothing.


Hamish



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