[GRASS-user] r.horizon not working with single step

[Lorenzo Bottaccioli]

Hi Everyone,

I'm running r.horizon with this command in grass 7.0 on unbuntu but doesn't
work

*(Tue Aug  4 12:29:35 2015)
 *
*r.horizon elevation=elevation direction=270 maxdistance=1000
output=horizon distance=0.5*
*(Tue Aug  4 12:29:36 2015) Comando terminato (0 sec) *

If I try with multistep it works

*(Tue Aug  4 12:30:10 2015)
 *
*r.horizon elevation=elevation step=20 maxdistance=1000 output=horizon
distance=0.5*
*Calculating map 1 of 18 (angle 0.00, raster map horizon_000)*

Any suggestion?

Tnx

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

Re: [GRASS-user] External raster from PostGIS

[Lorenzo Bottaccioli]

Hi Markus,

I don't see any copy button in the GUI of r.external.

Tnx

2015-10-23 17:20 GMT+02:00 Markus Neteler <nete...@osgeo.org>:

>
> On Oct 23, 2015 4:40 PM, "Lorenzo Bottaccioli"  <lorenzo.bottacci...@gmail.com>@ <lorenzo.bottacci...@gmail.com>gmail.com
> <lorenzo.bottacci...@gmail.com>> wrote:
> >
> > Hi Markus,
> >
> > I have tried to read the manual of r.external but is not clear how it
> works with DB
>
> as mentioned I would like to improve the manual. I am aware that your task
> is not explained yet.
>
> > , is very clear how it works with files. Using the GUI of r.external and
> connecting to the db i cannot see the raster map.
>
> And for this I asked you to tell us what you tried. There is a copy button
> in the GUI for that purpose.
>
> Best
> Markus
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] External raster from PostGIS

[Lorenzo Bottaccioli]

Hi Markus,

I did not try with any command in the command line, I'have just tried with
the GUI of r.in.gdal selecting database. Can you give a suggestion for the
use of "source=" option in the command line of r.external in order to
select data from a PostGIS db?

Tnx

Lorenzo

2015-10-22 18:34 GMT+02:00 Markus Neteler <nete...@osgeo.org>:

> Hi Lorenzo,
>
> On Thu, Oct 22, 2015 at 6:03 PM, Lorenzo Bottaccioli
> <lorenzo.bottacci...@gmail.com> wrote:
> > Hi Markus,
> >
> > How do I implement such command in GRASS GIS? I have tried with
> r.external
> > because I have done the same thing with v.external. I'have tried as well
> > with r.in.gdal but notthing. If i use v.in.ogr I can import vectors from
> > PostGIS.
>
> I would just like to the command to use. And the error message.
> Perhaps we see why it fails.
>
> Did you use "source=..." of r.external?
>
> Markus
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] External raster from PostGIS

[Lorenzo Bottaccioli]

Hello,

I'm tryng to import a raster file sotred in PostGIS with r.external but i
cannot select the db. If I try with r.in.gdal and i connect to the db I
can't see the table.

Any suggestion?

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

Re: [GRASS-user] External raster from PostGIS

[Lorenzo Bottaccioli]

Hi Markus,

How do I implement such command in GRASS GIS? I have tried with r.external
because I have done the same thing with v.external. I'have tried as well
with r.in.gdal but notthing. If i use v.in.ogr I can import vectors from
PostGIS.

Lorenzo

2015-10-22 17:55 GMT+02:00 Markus Neteler <nete...@osgeo.org>:

> On Thu, Oct 22, 2015 at 12:49 PM, Lorenzo Bottaccioli
> <lorenzo.bottacci...@gmail.com> wrote:
> > Hello,
> >
> > I'm tryng to import a raster file sotred in PostGIS with r.external but i
> > cannot select the db. If I try with r.in.gdal and i connect to the db I
> > can't see the table.
>
> Can you please post the full command? (omitting the password if given
> there)
> Is it similar to what's given in
>
> https://trac.osgeo.org/gdal/wiki/frmts_wtkraster.html#a3.2-Readingrasterdatafromthedatabase
> ?
>
> I would like to add an example to the manual page but do not have a PG
> installation containing a raster map.
>
> Markus
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] External raster from PostGIS

[Lorenzo Bottaccioli]

Hi Markus,

I have done what you have asked me. But the GUi that you obtain with
r.external --ui, asu suggested by Helmut, does not have any possibility to
select the DB so I have tried taking inspiration from v.external and I have
wrtitetn this comand that obviuosly does not work "r.external
source="PG:dbname=test_gis2 user=lorenzo" output=elev" where can I specify
the layer that I have to select? can you give me an example of a wright use
of the command?

best

Lorenzo

2015-10-23 18:40 GMT+02:00 Markus Neteler <nete...@osgeo.org>:

> On Fri, Oct 23, 2015 at 5:43 PM, Lorenzo Bottaccioli
> <lorenzo.bottacci...@gmail.com> wrote:
> > Hi Markus,
> >
> > I don't see any copy button in the GUI of r.external.
>
> If you launch it from the command line, without parameters, you get
> the "native" GUI of r.external which comes with the copy button.
>
> best
> Markus
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] External Database and projections

[Lorenzo Bottaccioli]

Hello,

Thanks for your suggestion. The point is that I need do develop a
distributed infrastructure hence I need to store the raster in to a DB
after the creation. Once I have created the raster immages I'll need to do
some simple multiplication and I thought to do them with r.mapcalc. Then
I'll use a vector map in order to extract the values of the raster.

So i think I'll have to use inevitably a DBMS to store the raster images
produced with grass. I'm I wright?

Tnx again regards

Lorenzo

2015-10-08 8:40 GMT+02:00 Blumentrath, Stefan :

> Hi Dylan,
>
>
>
> Yes, r.external and r.external.out works pretty well with TGIS. Here is a
> discussion in this regards I had mainly with Sören…
>
> https://lists.osgeo.org/pipermail/grass-user/2013-December/069377.html
>
>
>
> I have used r.external/r.external.out in TGIS for a while (with compressed
> GeoTiffs, btw you should not pre-define the predictor for LZW-compression
> as long as you are not sure that all raster maps are of the same type
> (CELL/FCELL/DCELL)), but went back to using internal GRASS format, because
> external links are a bit tricky to use on NFS and from different OS (or if
> data is moved):
>
> https://trac.osgeo.org/grass/ticket/2660
>
>
>
> Another minor (and often neglectable) drawback of using r.external is that
> file names and map names get out of sync if you rename maps. g.rename
> renames only the GRASS internal map name (modifies the link) the underlying
> data file is not touched and stays with the old name…
>
>
>
> Cheers
>
> Stefan
>
>
>
>
>
> *From:* Dylan Beaudette [mailto:dylan.beaude...@gmail.com]
> *Sent:* 8. oktober 2015 00:17
> *To:* Markus Neteler 
> *Cc:* Blumentrath, Stefan ; GRASS user list <
> grass-user@lists.osgeo.org>
> *Subject:* Re: [GRASS-user] External Database and projections
>
>
>
> The use of r.external makes a lot of sense when dealing with very large
> files. Does the use of "external" files work as expected in all of the new
> t.* modules?
>
>
>
> Thanks,
>
> Dylan
>
>
>
> On Wed, Oct 7, 2015 at 12:09 PM, Markus Neteler  wrote:
>
> On Wed, Oct 7, 2015 at 12:55 PM, Blumentrath, Stefan
>  wrote:
> ...
> > My suggestion is to not use PostGIS for big rasters, unless you have to,
> > because you want to use the data in a specific application for example.
>
> Note r.external and r.external.out of GRASS GIS 7 for avoiding data
> duplication:
>
> https://grass.osgeo.org/grass70/manuals/r.external.html
> https://grass.osgeo.org/grass70/manuals/r.external.out.html
>
> Here an example:
>
> # register (rather than import) a GeoTIFF file in GRASS GIS:
> r.external input=terra_lst1km20030314.LST_Day.tif output=modis_celsius
>
> # define output directory for files resulting from subsequent calculations:
> r.external.out directory=$HOME/gisoutput/ format="GTiff"
>
> # perform calculations (here: extract pixels > 20 deg C)
> # store output directly as GeoTIFF file, hence add the .tif extension:
> r.mapcalc "warm.tif = if(modis_celsius > 20.0, modis_celsius, null() )"
>
> # cease GDAL output connection and turn back to write standard GRASS
> raster files:
> r.external.out -r
>
> # use the result elsewhere
> qgis $HOME/gisoutput/warm.tif
>
> Best
> 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
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] External Database and projections

[Lorenzo Bottaccioli]

Hi Stefan,

With a distributed infrastructure I mean that the the computation and the
data are deployed in a distributed platform. The data then will be
available trough WWW but the aim of the infrastructure is not on the
visualization of the maps but on providing the information in JSON or
similar format regarding the information of the raster maps.

Cheers

Lorenzo

2015-10-14 16:09 GMT+02:00 Blumentrath, Stefan <stefan.blumentr...@nina.no>:

> Hi Lorenzo,
>
>
>
> What do you mean with a «distributed infrastructure”?
>
>
>
> Do you want to make the rasters available via WWW or in your
> company/institution?
>
> There are other options to share data than PostGIS… Within an organization
> your GRASS GIS database on NFS would be a much, much more efficient one.
>
> For WWW-publishing WCS is an option…
>
>
>
> Anyway, writing copies of your raster maps to DB (PostGIS) – if that is a
> fix requirement - is less of a problem (if they are not too heavy), than
> writing to DB, reading from DB for processing in GRASS and writing result
> back to DB…
>
>
>
> If your aim is to add “zonal statistics” to a vector layer in PostGIS I
> would do the raster operation in GRASS and then join the result to PG…
>
>
>
> Cheers
>
> Stefan
>
>
>
>
>
> *From:* Lorenzo Bottaccioli [mailto:lorenzo.bottacci...@gmail.com]
> *Sent:* 14. oktober 2015 15:44
> *To:* Blumentrath, Stefan <stefan.blumentr...@nina.no>
> *Cc:* Dylan Beaudette <dylan.beaude...@gmail.com>; Markus Neteler <
> nete...@osgeo.org>; GRASS user list <grass-user@lists.osgeo.org>
>
> *Subject:* Re: [GRASS-user] External Database and projections
>
>
>
> Hello,
>
> Thanks for your suggestion. The point is that I need do develop a
> distributed infrastructure hence I need to store the raster in to a DB
> after the creation. Once I have created the raster immages I'll need to do
> some simple multiplication and I thought to do them with r.mapcalc. Then
> I'll use a vector map in order to extract the values of the raster.
>
> So i think I'll have to use inevitably a DBMS to store the raster images
> produced with grass. I'm I wright?
>
> Tnx again regards
>
> Lorenzo
>
>
>
> 2015-10-08 8:40 GMT+02:00 Blumentrath, Stefan <stefan.blumentr...@nina.no
> >:
>
> Hi Dylan,
>
>
>
> Yes, r.external and r.external.out works pretty well with TGIS. Here is a
> discussion in this regards I had mainly with Sören…
>
> https://lists.osgeo.org/pipermail/grass-user/2013-December/069377.html
>
>
>
> I have used r.external/r.external.out in TGIS for a while (with compressed
> GeoTiffs, btw you should not pre-define the predictor for LZW-compression
> as long as you are not sure that all raster maps are of the same type
> (CELL/FCELL/DCELL)), but went back to using internal GRASS format, because
> external links are a bit tricky to use on NFS and from different OS (or if
> data is moved):
>
> https://trac.osgeo.org/grass/ticket/2660
>
>
>
> Another minor (and often neglectable) drawback of using r.external is that
> file names and map names get out of sync if you rename maps. g.rename
> renames only the GRASS internal map name (modifies the link) the underlying
> data file is not touched and stays with the old name…
>
>
>
> Cheers
>
> Stefan
>
>
>
>
>
> *From:* Dylan Beaudette [mailto:dylan.beaude...@gmail.com]
> *Sent:* 8. oktober 2015 00:17
> *To:* Markus Neteler <nete...@osgeo.org>
> *Cc:* Blumentrath, Stefan <stefan.blumentr...@nina.no>; GRASS user list <
> grass-user@lists.osgeo.org>
> *Subject:* Re: [GRASS-user] External Database and projections
>
>
>
> The use of r.external makes a lot of sense when dealing with very large
> files. Does the use of "external" files work as expected in all of the new
> t.* modules?
>
>
>
> Thanks,
>
> Dylan
>
>
>
> On Wed, Oct 7, 2015 at 12:09 PM, Markus Neteler <nete...@osgeo.org> wrote:
>
> On Wed, Oct 7, 2015 at 12:55 PM, Blumentrath, Stefan
> <stefan.blumentr...@nina.no> wrote:
> ...
> > My suggestion is to not use PostGIS for big rasters, unless you have to,
> > because you want to use the data in a specific application for example.
>
> Note r.external and r.external.out of GRASS GIS 7 for avoiding data
> duplication:
>
> https://grass.osgeo.org/grass70/manuals/r.external.html
> https://grass.osgeo.org/grass70/manuals/r.external.out.html
>
> Here an example:
>
> # register (rather than import) a GeoTIFF file in GRASS GIS:
> r.external input=terra_lst1km20030314.LST_Day.tif output=modis_celsius
>
> # define out

[GRASS-user] External Database and projections

[Lorenzo Bottaccioli]

Hello everyone,

I not very expert in this field and searching on the web didn't clarify my
dubs.
I wanted to store a raster DEM and some vectors map in to an external DB,
GRASS will use the DEM in order to calculate various raster map of solar
radiation and need to store the again in the external DB. Than I'll use
r.mapcal to do some calculation of on the raster maps and i'll use the
vector map to extract this information. I don't need to store the
calculation I just need to perform the calculation and use the results
immediately.

My question is which could be a feasible solution for this? how do I manage
the fact that the input raster and vector map can have different
projection? What about if I want to define a default projection of the DB
and convert the maps when I load them in the DB?

The DB will probably be Postgis, do you have any other suggestion?

Thanks a lot

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

Re: [GRASS-user] v.rast.stats error

[Lorenzo Bottaccioli]

I'have created the vector with r.to.vect, before I have liked
v.exeternel.out to the PostGIS DB. If I run v.what.rast there is no problem
but I believe that the results are different.

2015-11-20 17:25 GMT+01:00 Martin Landa <landa.mar...@gmail.com>:

> Hi,
>
> 2015-11-20 17:24 GMT+01:00 Lorenzo Bottaccioli <
> lorenzo.bottacci...@gmail.com>:
> > v.rast.stats -c map=areasuit@lorenzo raster=opr@lorenzo
> column_prefix=opr_
> > method=average
> > WARNING: No areas selected from vector map 
> > ERROR: No categories found in raster map
> >
> > The vector is contained in a PostGIS db.
>
> how did you connect vector map from PostGIS DB, using v.external? Ma
>
> --
> Martin Landa
> http://geo.fsv.cvut.cz/gwiki/Landa
> http://gismentors.cz/mentors/landa
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] v.rast.stats error

[Lorenzo Bottaccioli]

Hi,

I'm reciving this error if I v.rast.stats :

*v.rast.stats -c map=areasuit@lorenzo raster=opr@lorenzo column_prefix=opr_
method=average*
*WARNING: No areas selected from vector map *
*ERROR: No categories found in raster map*

The vector is contained in a PostGIS db.

Any suggestion?

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

Re: [GRASS-user] v.rast.stats error

[Lorenzo Bottaccioli]

Hi Martin,

Here is the output of v.info. I need to use PostGIS connection because I
have to expose the results via web-service. The fact is that I multiple
raster file for which I have to perform some simple mapalgebra, and then
querry the results from them and expose them trough web services. Any
suggestion for this?

 ++
 | Name:areasuit
   |
 | Mapset:  lorenzo
  |
 | Location:prova
  |
 | Database:/home/lorenzo/grassData
  |
 | Title:
  |
 | Map scale:   1:1
  |
 | Name of creator: lorenzo
  |
 | Organization:
   |
 | Source date: Fri Nov 20 17:15:02 2015
   |
 | Timestamp (first layer): none
   |
 ||
 | Map format:  PostGIS (PostgreSQL)
   |
 | DB table:public.areasuit
  |
 | DB name: gis_1
  |
 | Geometry column: geom
   |
 | Feature type:polygon
  |
 | Topology:pseudo (simple features)
   |
 ||
 |   Type of map: vettoriale (level: 2)
  |
 |
   |
 |   Number of points:   0   Number of centroids:  4002
  |
 |   Number of lines:0   Number of boundaries: 4022
  |
 |   Number of areas:4022Number of islands:4022
  |
 |
   |
 |   Map is 3D:  No
  |
 |   Number of dblinks:  1
   |
 |
   |
 |   Projection: Latitudine-Longitudine
  |
 |
   |
 |   N:  45:04:24.507374NS:  45:03:03.227962N
  |
 |   E:   7:40:23.868382EW:   7:38:42.814094E
  |
 |
   |
 |   Digitization threshold: 0
   |
 |   Comment:
  |
 |
   |
 ++

2015-11-20 17:33 GMT+01:00 Martin Landa <landa.mar...@gmail.com>:

> Hi,
>
> 2015-11-20 17:28 GMT+01:00 Lorenzo Bottaccioli <
> lorenzo.bottacci...@gmail.com>:
> > I'have created the vector with r.to.vect, before I have liked
> > v.exeternel.out to the PostGIS DB. If I run v.what.rast there is no
> problem
> > but I believe that the results are different.
>
> OK, now I understand, v.what.rast operates on points/centoids,
> v.rast.stats on areas. There is something wrong with areas. Generally
> speaking I cannot suggest to use PostGIS links in GRASS for production
> work. Use native format instead.
>
> Ma
>
> --
> Martin Landa
> http://geo.fsv.cvut.cz/gwiki/Landa
> http://gismentors.cz/mentors/landa
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] r.mapcalc multiprocess

[Lorenzo Bottaccioli]

Hi every one,

I have just compiled grass with phtreds but when I run r.mapcalc it runs in
only with one theard. How ca I fix this? Is it possible to check if is
correclty compiled with phtreads?

Tnx

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

[GRASS-user] Pygrass reading raster issue

[Lorenzo Bottaccioli]

Hi all,

I'm tring to read a raster map from my grassData.

I'm using this code

import os
import grass.script.setup as gsetup
from grass.pygrass import raster

gisbase = os.environ['GISBASE']
gisdbase = os.path.abspath("/home/myuser/grassData")
location = 'mylocation'
mapset = 'mymapset'
gsetup.init(gisbase, gisdbase, location, mapset)

elev=raster.RasterRow('elevation')
elev.open()
for row in elev:
print row

But I get all values with nan, where is the error?

If I do `elev.info()` I get this result

elevation@
rows: 4502
cols: 5554
north: 45.0743663458 south: 45.0500262069 nsres:5.40651685671e-06
east:  7.67457811722 west: 7.6445503225 ewres:5.40651687473e-06
range: 233.050003052, 310.049987793
proj: 3

Tnx

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

Re: [GRASS-user] Pygrass reading raster issue

[Lorenzo Bottaccioli]

 g.region('-p')
Out[40]:
{'b': 0.0,
 'cells': 25004108,
 'cells3': 25004108,
 'cols': 5554,
 'cols3': 5554,
 'depths': 1,
 'e': 7.6745781172,
 'ewres': 5.40651687472501e-06,
 'ewres3': 5.40651687472501e-06,
 'n': 45.0743663458333,
 'nsres': 5.40651685670714e-06,
 'nsres3': 5.40651685670714e-06,
 'rows': 4502,
 'rows3': 4502,
 's': 45.0500262069444,
 't': 1.0,
 'tbres': 1.0,
 'w': 7.6445503225}


before I have set g.region('raster=elevation')

So the region should be OK!

2016-03-14 18:10 GMT+01:00 Helmut Kudrnovsky :

>  >> [...] nsres:5.40651685671e-06
>  >> [...] ewres:5.40651687473e-06
> > Region is Ok!
>
> g.region -p ?
>
>
>
>
> -
> best regards
> Helmut
> --
> View this message in context:
> http://osgeo-org.1560.x6.nabble.com/Pygrass-reading-raster-issue-tp5256213p5256279.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 mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Compile error

[Lorenzo Bottaccioli]

Hi,

I can not compile grass-gis from source. I have download grass-gis 7.0.3
source and I have run configure with ./configure
--with-gdal=/usr/local/bin/gdal-config --with-pthread --with-opengl
--with-openmp --with-opencl --with-freetype-includes=/usr/include/freetype2
--with-geos.

When i do make i get this error:

Started compilation: mer  9 mar 2016, 11.16.21, CET
--
Errors in:
/home/lorenzo/tools/grass-7.0.3/display/d.path
/home/lorenzo/tools/grass-7.0.3/display/d.vect
/home/lorenzo/tools/grass-7.0.3/display/d.vect.chart
/home/lorenzo/tools/grass-7.0.3/display/d.vect.thematic
/home/lorenzo/tools/grass-7.0.3/general/g.copy
/home/lorenzo/tools/grass-7.0.3/general/g.findfile
/home/lorenzo/tools/grass-7.0.3/general/g.list
/home/lorenzo/tools/grass-7.0.3/general/g.region
/home/lorenzo/tools/grass-7.0.3/general/g.remove
/home/lorenzo/tools/grass-7.0.3/general/g.rename
/home/lorenzo/tools/grass-7.0.3/general/manage/lister
/home/lorenzo/tools/grass-7.0.3/raster/r.carve
/home/lorenzo/tools/grass-7.0.3/raster/r.contour
/home/lorenzo/tools/grass-7.0.3/raster/r.cost
/home/lorenzo/tools/grass-7.0.3/raster/r.drain
/home/lorenzo/tools/grass-7.0.3/raster/r.external
/home/lorenzo/tools/grass-7.0.3/raster/r.flow
/home/lorenzo/tools/grass-7.0.3/raster/r.in.gdal
/home/lorenzo/tools/grass-7.0.3/raster/r.out.gdal
/home/lorenzo/tools/grass-7.0.3/raster/r.random
/home/lorenzo/tools/grass-7.0.3/raster/r.reclass
/home/lorenzo/tools/grass-7.0.3/raster/r.region
/home/lorenzo/tools/grass-7.0.3/raster/r.resamp.bspline
/home/lorenzo/tools/grass-7.0.3/raster/r.resamp.rst
/home/lorenzo/tools/grass-7.0.3/raster/r.sim/r.sim.water
/home/lorenzo/tools/grass-7.0.3/raster/r.sim/r.sim.sediment
/home/lorenzo/tools/grass-7.0.3/raster/r.stream.extract
/home/lorenzo/tools/grass-7.0.3/raster/r.to.vect
/home/lorenzo/tools/grass-7.0.3/raster/r.volume
/home/lorenzo/tools/grass-7.0.3/raster/r.walk
/home/lorenzo/tools/grass-7.0.3/raster/r.what
/home/lorenzo/tools/grass-7.0.3/vector/v.buffer
/home/lorenzo/tools/grass-7.0.3/vector/v.build
/home/lorenzo/tools/grass-7.0.3/vector/v.build.polylines
/home/lorenzo/tools/grass-7.0.3/vector/v.category
/home/lorenzo/tools/grass-7.0.3/vector/v.class
/home/lorenzo/tools/grass-7.0.3/vector/v.clean
/home/lorenzo/tools/grass-7.0.3/vector/v.colors
/home/lorenzo/tools/grass-7.0.3/vector/v.colors.out
/home/lorenzo/tools/grass-7.0.3/vector/v.db.connect
/home/lorenzo/tools/grass-7.0.3/vector/v.db.select
/home/lorenzo/tools/grass-7.0.3/vector/v.delaunay
/home/lorenzo/tools/grass-7.0.3/vector/v.distance
/home/lorenzo/tools/grass-7.0.3/vector/v.drape
/home/lorenzo/tools/grass-7.0.3/vector/v.edit
/home/lorenzo/tools/grass-7.0.3/vector/v.extract
/home/lorenzo/tools/grass-7.0.3/vector/v.extrude
/home/lorenzo/tools/grass-7.0.3/vector/v.generalize
/home/lorenzo/tools/grass-7.0.3/vector/v.hull
/home/lorenzo/tools/grass-7.0.3/vector/v.info
/home/lorenzo/tools/grass-7.0.3/vector/v.in.ascii
/home/lorenzo/tools/grass-7.0.3/vector/v.in.db
/home/lorenzo/tools/grass-7.0.3/vector/v.in.dxf
/home/lorenzo/tools/grass-7.0.3/vector/v.in.region
/home/lorenzo/tools/grass-7.0.3/vector/v.kcv
/home/lorenzo/tools/grass-7.0.3/vector/v.kernel
/home/lorenzo/tools/grass-7.0.3/vector/v.label
/home/lorenzo/tools/grass-7.0.3/vector/v.lidar.correction
/home/lorenzo/tools/grass-7.0.3/vector/v.lidar.edgedetection
/home/lorenzo/tools/grass-7.0.3/vector/v.lidar.growing
/home/lorenzo/tools/grass-7.0.3/vector/v.lrs/v.lrs.create
/home/lorenzo/tools/grass-7.0.3/vector/v.lrs/v.lrs.segment
/home/lorenzo/tools/grass-7.0.3/vector/v.lrs/v.lrs.label
/home/lorenzo/tools/grass-7.0.3/vector/v.lrs/v.lrs.where
/home/lorenzo/tools/grass-7.0.3/vector/v.proj
/home/lorenzo/tools/grass-7.0.3/vector/v.mkgrid
/home/lorenzo/tools/grass-7.0.3/vector/v.neighbors
/home/lorenzo/tools/grass-7.0.3/vector/v.net
/home/lorenzo/tools/grass-7.0.3/vector/v.net.alloc
/home/lorenzo/tools/grass-7.0.3/vector/v.net.allpairs
/home/lorenzo/tools/grass-7.0.3/vector/v.net.bridge
/home/lorenzo/tools/grass-7.0.3/vector/v.net.centrality
/home/lorenzo/tools/grass-7.0.3/vector/v.net.components
/home/lorenzo/tools/grass-7.0.3/vector/v.net.connectivity
/home/lorenzo/tools/grass-7.0.3/vector/v.net.distance
/home/lorenzo/tools/grass-7.0.3/vector/v.net.flow
/home/lorenzo/tools/grass-7.0.3/vector/v.net.iso
/home/lorenzo/tools/grass-7.0.3/vector/v.net.path
/home/lorenzo/tools/grass-7.0.3/vector/v.net.salesman
/home/lorenzo/tools/grass-7.0.3/vector/v.net.spanningtree
/home/lorenzo/tools/grass-7.0.3/vector/v.net.steiner
/home/lorenzo/tools/grass-7.0.3/vector/v.net.timetable
/home/lorenzo/tools/grass-7.0.3/vector/v.net.visibility
/home/lorenzo/tools/grass-7.0.3/vector/v.normal
/home/lorenzo/tools/grass-7.0.3/vector/v.out.ascii
/home/lorenzo/tools/grass-7.0.3/vector/v.out.dxf
/home/lorenzo/tools/grass-7.0.3/vector/v.out.pov
/home/lorenzo/tools/grass-7.0.3/vector/v.out.svg
/home/lorenzo/tools/grass-7.0.3/vector/v.out.vtk
/home/lorenzo/tools/grass-7.0.3/vector/v.outlier
/home/lorenzo/tools/grass-7.0.3/vector/v.overlay

[GRASS-user] r.external.out store data in memory

[Lorenzo Bottaccioli]

Hi list,

I wanted to know if was possible and how to store produced rasters in
memory (as an rasterio object or similar) and do not store/save the in the
disk, because I just need to do some statistics and then delete them. This
if possibele will speed up my workflow because I'll avoid I/O from the disk.

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

Re: [GRASS-user] r.mapcalc multiprocess

[Lorenzo Bottaccioli]

Hi Stefan,

I have just updated my grass-gis version to 7.0.3, berfore I was using
7.0.0 version and mapcalc was running in multitheard mode and was much
faster. My question is how can i make it run again in multitheading?

Regards,

Lorenzo

2016-03-10 22:39 GMT+01:00 Blumentrath, Stefan <stefan.blumentr...@nina.no>:

> Hi,
>
>
>
> As far as I know it is not exactly recommended to compile with pthread,
> and if I recall the latest discussion correctly, the number of default
> workers was deliberately set to 1 due to issues with multithreaded
> r.mapcalc.
>
>
>
> Don`t know which was the latest discussion on this topic, but it has been
> discussed / explained in detail on the developer list, e.g.:
>
> https://lists.osgeo.org/pipermail/grass-dev/2012-June/058749.html
>
> https://lists.osgeo.org/pipermail/grass-dev/2008-November/041047.html
>
>
>
> Hope that helps a bit.
>
>
>
> Kind regards,
>
> Stefan
>
>
>
> *From:* grass-user [mailto:grass-user-boun...@lists.osgeo.org] *On Behalf
> Of *Lorenzo Bottaccioli
> *Sent:* 10. mars 2016 17:14
> *To:* GRASS user list <grass-user@lists.osgeo.org>
> *Subject:* [GRASS-user] r.mapcalc multiprocess
>
>
>
> Hi every one,
>
> I have just compiled grass with phtreds but when I run r.mapcalc it runs
> in only with one theard. How ca I fix this? Is it possible to check if is
> correclty compiled with phtreads?
>
> Tnx
>
> Lorenzo
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Pygrass reading raster issue

[Lorenzo Bottaccioli]

Region is Ok!

2016-03-14 14:45 GMT+01:00 Moritz Lennert <mlenn...@club.worldonline.be>:

> Le Mon, 14 Mar 2016 13:45:42 +0100,
> Lorenzo Bottaccioli <lorenzo.bottacci...@gmail.com> a écrit :
>
> > Hi all,
> >
> > I'm tring to read a raster map from my grassData.
> >
> > I'm using this code
> >
> > import os
> > import grass.script.setup as gsetup
> > from grass.pygrass import raster
> >
> > gisbase = os.environ['GISBASE']
> > gisdbase = os.path.abspath("/home/myuser/grassData")
> > location = 'mylocation'
> > mapset = 'mymapset'
> > gsetup.init(gisbase, gisdbase, location, mapset)
> >
> > elev=raster.RasterRow('elevation')
> > elev.open()
> > for row in elev:
> > print row
> >
> > But I get all values with nan, where is the error?
> >
> > If I do `elev.info()` I get this result
> >
> > elevation@
> > rows: 4502
> > cols: 5554
> > north: 45.0743663458 south: 45.0500262069 nsres:5.40651685671e-06
> > east:  7.67457811722 west: 7.6445503225 ewres:5.40651687473e-06
> > range: 233.050003052, 310.049987793
> > proj: 3
> >
>
>
> Check your computational region settings.
>
> Moritz
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Parallel proccess with pygrass and mapcalc

[Lorenzo Bottaccioli]

Hi I'm writing a code in python able to perform several mapcalc operation
on various rasters. I'm using GRASS 7.0.1 and Python 2.7. I'have my files
stored in /home/myuser/rasters/ and they are linked to GRASS with
r.external. The out put of the operation are going out of GRASS in
/home/myuser/tmp/ with r.external.out.

I have written this code from the example provided in (
https://grass.osgeo.org/grass70/manuals/libpython/pygrass.modules.interface.html?highlight=parallelmodulequeue#pygrass.modules.interface.module.ParallelModuleQueue)
 but there is no difference in terms of time between setting nproces=1 or
nproces=8. How is possible? How can I achieve improvements?

The task that the code need to do is to read data from a file containing
datetime and 2 coefficients k1 and k2. Than need to multiply and sum the
input raster as out = input1 * k1 + input2 * k2.

The code that i have written is:


import pandas as pd
import os
import grass.script.setup as gsetup
import grass.script as g
import time
import copy
from grass.pygrass.modules import Module, ParallelModuleQueue

def main():

   #set GRASS LOCATION AND MAPSET
gisbase = os.environ['GISBASE']   # Grass 7.0svn
gisdbase = os.path.abspath("/home/myuser/grassData")
location = 'mylocation' # Grass Location.
mapset = 'mymapset'
gsetup.init(gisbase, gisdbase, location, mapset)

   #READ INPUT DATA FOR RASTER CALC
df = pd.read_csv('input.csv', sep=";", index_col='Date Time',
decimal=',')
df.index = pd.to_datetime(df.index, unit='s')


month={1:'17',2:'47',3:'75',4:'105',5:'135',6:'162',7:'198',8:'228',9:'258',10:'288',11:'318',12:'344'}

hour={4:'04',5:'05',6:'06',7:'07',8:'08',9:'09',10:'10',11:'11',12:'12',13:'13',14:'14',15:'15',16:'16',17:'17',18:'18',19:'19',20:'20',21:'21',22:'22'}
minute={0:'00',15:'15',30:'30',45:'45'}
directory='/home/myuser/raster/'
tmp='/home/myuser/tmp/'
g.run_command('r.external.out', directory=tmp, format="GTiff")
 #MAPCALC
start_time = time.time()
mapcalc_list = []
mapcalc = Module("r.mapcalc", overwrite=True, run_=False)
queue = ParallelModuleQueue(nprocs=8)
for dfix in df.index:
if 5<=dfix.hour<20:
input1 = 'input1_' + month[dfix.month] + '_' + hour[dfix.hour]
+ minute[dfix.minute]
input2 = 'input2_' + month[dfix.month] + '_' + hour[dfix.hour]
+ minute[dfix.minute]
out = ' " ' + str(dfix.date()) + '_' + str(dfix.time()) + ' " '
new_mapcalc = copy.deepcopy(mapcalc)
mapcalc_list.append(new_mapcalc)
m = new_mapcalc(expression="%s =
%s*%i+%s*%i"%(out,input1,df.ix[dfix,'k1'],input2,df.ix[dfix,'k2']))
queue.put(m)
queue.wait()


print("--- %s seconds ---" % (time.time() - start_time))



Best,

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

[GRASS-user] Workflow for DSM from LiDAR

[Lorenzo Bottaccioli]

High list,

I have a dense LiDAR data ~8pt/m2, how can I process this data in order to
obtain a dense digital surface model for example 50cm of relsolution? I'll
like to have good detail of the rooftops for solar analysis , in particolar
I'll like to have not to much noise in the slope and aspect maps obtained
from the DSM.

Thanks,

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

[GRASS-user] r.sun and OpenCL not working

[Lorenzo Bottaccioli]

Hi list,

I have configure grass7.04 -with-opencl
--with-opencl-includes=/path/include --with-opencl-libs=/path/lib and I try
to run it on a Intel Xeon dual socket 32-thread. The compilation say that
OpenCl support is activated but when i run r.sun I don't see any
multi-theard I uses only one CPU.

Any suggestion to solve the issue?

Best,

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

[GRASS-user] Interpolate daily rainfall data over DEM

[Lorenzo Bottaccioli]

Hi all,

Which is the best tool or tool-chain to interpolate the daily rainfall in
(mm) over a DEM? What if I want to evaluate the overall presence of water
in the terrain taking in to account evaporation?

Ciao,

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