Re: [GRASS-user] detrending time series maps

[Ivan Marchesini via grass-user]

Hi veronica



I see what you did with creating the days time series. In that way you 
acknowledge irregular gaps, right?

yes this is the reason


However, why do you multiply by days strds? From my understanding, 
detrending by subtracting the results of a model obbeys this rule: 
value(t) = observed(t) - predicted(t). Then, this 
mystrds-(regression_offset+regression_slope*days) should be 
mystrds-(regression_offset+regression_slope*mystrds).


hmm

may be I'm wrong but I used r.regression.series to assess a relationship 
between mystrds values and time (days)


As a consequence the offset and slope maps I obtain are b and a in the 
following linear equation


y=ax+b

i.e.

predicted_mystrds=slope*days+offset

This is why I suppose that for detrending I need to do:

mystrds-(predicted_mystrds)

i.e.

mystrds-(slope*days+offset)


I'm I wrong?

thank you

Ivan





Best,
Vero

El jue, 28 dic 2023 a las 8:14, Ivan Marchesini 
() escribió:


Dear Veronica

I think I found a simple solution using temporal raster modules.
Here is an example:

#evaluating info of the strds
eval `t.info <http://t.info> mystrds -g`

#getting the starting day (of the year, 0-365) of my strds
startday=$(date -d "$start_time" "+%j")

#Creating a new strds where each pixel has the value of the count
of the days starting from the start_day of my strds (the start day
in my dataset is in 2016)
t.rast.mapcalc inputs=mystrds
expression="(start_year()-2016)*365-${startday} +start_doy()"
output=days basename=days nprocs=xxx --o

#fitting the trend equation
r.regression.series xseries="`t.rast.list in=days columns=name
sep=, format=line`" yseries="` t.rast.list in=mystrds columns=name
sep=, format=line`"
out=regression_offset,regression_slope,regression_rsq,regression_t
meth=offset,slope,rsq,t

#detrending
t.rast.mapcalc input=mystrds,days expression="mystrds_detrend
= mystrds-(regression_offset+regression_slope*days)"
output=mystrds_detrend basename=mystrds_detrend nprocs=xxx
    method=start --o


Best

Ivan




On 23/12/23 14:53, Ivan Marchesini wrote:


Hi Veronica

Thank you. It goes in the direction of my idea evn if my problem
is exactly trying to take into account the correct gaps between
that data

I have another idea.

if it works I will come back here to explain how I did

thank you again

Ivan


On 22/12/23 13:45, Veronica Andreo wrote:

Hello Ivan,

AFAIU you could use the slope and offset maps from t.rast.series
within t.rast.algebra to detrend the values of the maps within
the strds, something like "detrended_strds = trend_strds -
(trend_strds*map(slope) + map(offset))". Others suggest, to
detrend by subtracting the previous value, i.e. that would imply
using the temporal algebra with the temporal index, something
like "detrended_strds = trend_strds[1] - trend_strds[0]".

I haven't tested any of these, just a couple of ideas ;-)
However, I do not know how this might interact with seasonality
within data, or irregular gaps.

hth somehow
Vero

El vie, 22 dic 2023 a las 5:10, Ivan Marchesini via grass-user
() escribió:

Dear colleagues

I would like to the advantage of the t.* modules for
detrending a strd.

In the strd I have earth observation data irregularly
sampled (2 or 3
times per month), in the period November-February, for 7
years. They are
not equally spaced (i.e gaps have different duration)

A simple t.rast.series analysis (opion=slope,offset)
highlights that
probably there is a descending trend when considering the
maps ordered
by id.

I would like to fit a proper time depending fitting curve
for each pixel
and then subtract the function from the real data.

any hints on how I can do this task exploiting the GRASS GIS
modules or
some simple bash/python scripting?

thank you

Ivan




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

Re: [GRASS-user] detrending time series maps

[Ivan Marchesini via grass-user]

Dear Veronica

I think I found a simple solution using temporal raster modules. Here is 
an example:


#evaluating info of the strds
eval `t.info mystrds -g`

#getting the starting day (of the year, 0-365) of my strds
startday=$(date -d "$start_time" "+%j")

#Creating a new strds where each pixel has the value of the count of the 
days starting from the start_day of my strds (the start day in my 
dataset is in 2016)
t.rast.mapcalc inputs=mystrds 
expression="(start_year()-2016)*365-${startday} +start_doy()" 
output=days basename=days nprocs=xxx --o


#fitting the trend equation
r.regression.series xseries="`t.rast.list in=days columns=name sep=, 
format=line`" yseries="` t.rast.list in=mystrds columns=name sep=, 
format=line`" 
out=regression_offset,regression_slope,regression_rsq,regression_t 
meth=offset,slope,rsq,t


#detrending
t.rast.mapcalc input=mystrds,days expression="mystrds_detrend 
= mystrds-(regression_offset+regression_slope*days)" 
output=mystrds_detrend basename=mystrds_detrend nprocs=xxx method=start --o



Best

Ivan




On 23/12/23 14:53, Ivan Marchesini wrote:


Hi Veronica

Thank you. It goes in the direction of my idea evn if  my problem is 
exactly trying to take into account the correct gaps between that data


I have another idea.

if it works I will come back here to explain how I did

thank you again

Ivan


On 22/12/23 13:45, Veronica Andreo wrote:

Hello Ivan,

AFAIU you could use the slope and offset maps from t.rast.series 
within t.rast.algebra to detrend the values of the maps within the 
strds, something like "detrended_strds = trend_strds - 
(trend_strds*map(slope) + map(offset))". Others suggest, to detrend 
by subtracting the previous value, i.e. that would imply using the 
temporal algebra with the temporal index, something like 
"detrended_strds = trend_strds[1] - trend_strds[0]".


I haven't tested any of these, just a couple of ideas ;-) However, I 
do not know how this might interact with seasonality within data, or 
irregular gaps.


hth somehow
Vero

El vie, 22 dic 2023 a las 5:10, Ivan Marchesini via grass-user 
() escribió:


Dear colleagues

I would like to the advantage of the t.* modules for detrending a
strd.

In the strd I have earth observation data irregularly sampled (2
or 3
times per month), in the period November-February, for 7 years.
They are
not equally spaced (i.e gaps have different duration)

A simple t.rast.series analysis (opion=slope,offset) highlights that
probably there is a descending trend when considering the maps
ordered
by id.

I would like to fit a proper time depending fitting curve for
each pixel
and then subtract the function from the real data.

any hints on how I can do this task exploiting the GRASS GIS
modules or
some simple bash/python scripting?

thank you

Ivan




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

Re: [GRASS-user] detrending time series maps

[Ivan Marchesini via grass-user]

Hi Veronica

Thank you. It goes in the direction of my idea evn if  my problem is 
exactly trying to take into account the correct gaps between that data


I have another idea.

if it works I will come back here to explain how I did

thank you again

Ivan


On 22/12/23 13:45, Veronica Andreo wrote:

Hello Ivan,

AFAIU you could use the slope and offset maps from t.rast.series 
within t.rast.algebra to detrend the values of the maps within the 
strds, something like "detrended_strds = trend_strds - 
(trend_strds*map(slope) + map(offset))". Others suggest, to detrend by 
subtracting the previous value, i.e. that would imply using the 
temporal algebra with the temporal index, something like 
"detrended_strds = trend_strds[1] - trend_strds[0]".


I haven't tested any of these, just a couple of ideas ;-) However, I 
do not know how this might interact with seasonality within data, or 
irregular gaps.


hth somehow
Vero

El vie, 22 dic 2023 a las 5:10, Ivan Marchesini via grass-user 
() escribió:


Dear colleagues

I would like to the advantage of the t.* modules for detrending a
strd.

In the strd I have earth observation data irregularly sampled (2 or 3
times per month), in the period November-February, for 7 years.
They are
not equally spaced (i.e gaps have different duration)

A simple t.rast.series analysis (opion=slope,offset) highlights that
probably there is a descending trend when considering the maps
ordered
by id.

I would like to fit a proper time depending fitting curve for each
pixel
and then subtract the function from the real data.

any hints on how I can do this task exploiting the GRASS GIS
modules or
some simple bash/python scripting?

thank you

Ivan




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

[GRASS-user] detrending time series maps

[Ivan Marchesini via grass-user]

Dear colleagues

I would like to the advantage of the t.* modules for detrending a strd.

In the strd I have earth observation data irregularly sampled (2 or 3 
times per month), in the period November-February, for 7 years. They are 
not equally spaced (i.e gaps have different duration)


A simple t.rast.series analysis (opion=slope,offset) highlights that 
probably there is a descending trend when considering the maps ordered 
by id.


I would like to fit a proper time depending fitting curve for each pixel 
and then subtract the function from the real data.


any hints on how I can do this task exploiting the GRASS GIS modules or 
some simple bash/python scripting?


thank you

Ivan




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

Re: [GRASS-user] largest number of consecutive days where rainfall >1 mm

[Ivan Marchesini]

Dear Veronica

thank you for your answer

I understand your point but I still believe that the code is not working 
properly for my aims


I try to explain myself with an example.

According to the example in the tutorial:

the first command creates a weekly mask based on my data (daily 
precipitation in 2021):


t.rast.aggregate input=test2021 output=weekly_mask basename=mask_week 
granularity="1 week" method=count nprocs=10


The second command add to the weekly value (7) the value of the 
consecutive days in the week period


t.rast.algebra base=test_rr_maggiore_1mm 
expression="test_consecutive_days_rr_maggiore_1mm = weekly_mask 
{+,contains,l} if(test2021 < 1 && test2021[-1] < 1 || test2021[1] < 1 && 
test2021 < 1, 1, 0)" nprocs=10


In the attached file

https://drive.google.com/file/d/1eusufL8JLysac0wLfWC6oub0IM79e-KO/view?usp=sharing

I show an extract of the the values of a given pixel in the study area.

In the first week the value of the map 
"test_consecutive_days_rr_maggiore_1mm is 10  (7 + 3). 7 is the weekly 
base. 3 is the number of times in which a rainy day is preceded [-1] or 
followed [+1] by another rainy day. This happen in 2021-01-05, 
2021-01-09, 2021-01-10 (in fact 3 days).


What is not fitting my needs here is the fact that between the 
2021-01-05 and 2021-01-09  there are other days without rainfall. 
Actually we have not 3 consecutive days in this week!! We have 2 
consecutive days: 2021-01-09, 2021-01-10 and another rainfall day at the 
start of the week.


I would like to have, in this first week, the value of 9 (7+2) not 10 
(7+3), because I'm interested in the annual maximum length of wet spell, 
in days.


I hope I was able to explain myself

Thank you again

Ivan









On 18/08/22 15:36, Veronica Andreo wrote:

Hi Ivan,

It does indeed count consecutive days meeting a condition per week. 
That's what the [-1] means in the t.rast.algebra command, i.e., it is 
the temporal neighbourhood modifier. If you then want to get the 
largest weekly consecutive days meeting the condition within a month 
you can aggregate the "weekly_consecutive" time series with 
t.rast.aggregate granularity="1 month" method="maximum"


HTH,
Vero

El jue, 18 ago 2022 a las 14:42, Ivan Marchesini 
() escribió:


Hi Markus

thank you

Yes I saw that algorithm but, if I'm not wrong, it does not fit
with my
needs. The calculation made returns the number of days, in a given
interval (here week), that meets a certain condition. Not sure
they are
consecutive. As an example: if negative temperature occur in
monday and
thuesday and then in friday, saturday and sunday then the code
count 5
for that week. In my opinion what we need is 3.

I hope I was enough clear

In any case it would be good to find a solution for this type of
problem
because these type of indexes are  used for climate data analysis
(e.g.:
https://www.ecad.eu/download/millennium/millennium.php
https://www.climdex.org/learn/indices/ )

thank you

Ivan



On 17/08/22 20:04, Markus Neteler wrote:
> Hi Ivan,
>
> On Wed, Aug 17, 2022 at 9:07 AM Ivan Marchesini
>  wrote:
>> Dear grass user
>>
>> I have daily rainfall strds
>>
>> I would like to obtain a layer of the largest number of
consecutive days
>> where rainfall >1 mm
> May this script code help?
>

https://grasswiki.osgeo.org/wiki/Temporal_data_processing#How_to_count_consecutive_days_that_meet_a_certain_condition?
>
> Markus
>
>> Do you have any suggestions on how to calculate this climate
index using
>> grass's time modules?
>>
>> thank you very much
>>
>> Ivan
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] largest number of consecutive days where rainfall >1 mm

[Ivan Marchesini]

Hi Markus

thank you

Yes I saw that algorithm but, if I'm not wrong, it does not fit with my 
needs. The calculation made returns the number of days, in a given 
interval (here week), that meets a certain condition. Not sure they are 
consecutive. As an example: if negative temperature occur in monday and 
thuesday and then in friday, saturday and sunday then the code count 5 
for that week. In my opinion what we need is 3.


I hope I was enough clear

In any case it would be good to find a solution for this type of problem 
because these type of indexes are  used for climate data analysis (e.g.: 
https://www.ecad.eu/download/millennium/millennium.php 
https://www.climdex.org/learn/indices/ )


thank you

Ivan



On 17/08/22 20:04, Markus Neteler wrote:

Hi Ivan,

On Wed, Aug 17, 2022 at 9:07 AM Ivan Marchesini
 wrote:

Dear grass user

I have daily rainfall strds

I would like to obtain a layer of the largest number of consecutive days
where rainfall >1 mm

May this script code help?
https://grasswiki.osgeo.org/wiki/Temporal_data_processing#How_to_count_consecutive_days_that_meet_a_certain_condition?

Markus


Do you have any suggestions on how to calculate this climate index using
grass's time modules?

thank you very much

Ivan

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

[GRASS-user] largest number of consecutive days where rainfall >1 mm

[Ivan Marchesini]

Dear grass user

I have daily rainfall strds

I would like to obtain a layer of the largest number of consecutive days 
where rainfall >1 mm


Do you have any suggestions on how to calculate this climate index using 
grass's time modules?


thank you very much

Ivan

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

Re: [GRASS-user] temporal raster accumulation

[Ivan Marchesini]

Thank you Stefan and Markus

I'll check out which one of the different solutions works better for me

thank you again!!!

On 09/08/22 13:43, Markus Neteler wrote:

Hi Ivan,

In addition, see also

https://grasswiki.osgeo.org/wiki/Temporal_data_processing#Aggregation_of_seasonal_data_using_time_ranges

best,
Markus

On Tue, Aug 9, 2022 at 11:40 AM Stefan Blumentrath
 wrote:

Hi Ivan,

A simple, quick-and-dirty approach could be:
t.rast.algebra basename=rain5 expression="rain5=rain + rain[-1] + rain[-2] + 
rain[-3] + rain[-4] + rain[-5]"

Temporal selection (start_date(rain)...) is probably more elegant...

See:
https://grass.osgeo.org/grass78/manuals/t.rast.algebra.html

Cheers
Stefan

-Original Message-
From: grass-user  On Behalf Of Ivan 
Marchesini
Sent: tirsdag 9. august 2022 09:45
To: GRASS user list 
Subject: [GRASS-user] temporal raster accumulation

Dear GRASS users

I'm a little bit confused about which can be the best module/tool to perform 
the following task I have a strd having a daily granularity. It contains total 
daily cumulated rainfall.

I would like to create a new strd where, for each day, I have the cumulated 
rainfall in the antecedent 5 days

can you kindly suggest which is the best tool (or combination of tools) to 
reach this aim?

thank you very much

Ivan

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

[GRASS-user] temporal raster accumulation

[Ivan Marchesini]

Dear GRASS users

I'm a little bit confused about which can be the best module/tool to 
perform the following task
I have a strd having a daily granularity. It contains total daily 
cumulated rainfall.


I would like to create a new strd where, for each day, I have the 
cumulated rainfall in the antecedent 5 days


can you kindly suggest which is the best tool (or combination of tools) 
to reach this aim?


thank you very much

Ivan





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

Re: [GRASS-user] stvds

[Ivan Marchesini]

Thank you for your answer Veronica,

in the meantime I found an addons that does exactly what I need: sample 
a strds at specific point in space and time defined by means of a point 
vector layer having an attribute column with dates (tool prepared by 
@Luca Delucchi: thank you)


I put the link here for anyone that in the future may have a similar 
problem:


https://grass.osgeo.org/grass82/manuals/addons/t.rast.what.aggr.html

thank you again

Ivan



On 24/07/22 13:07, Veronica Andreo wrote:

Hello Ivan,

Space-time datasets in GRASS are map based, so, you either need to 
create one vector map per date and then create the time series and 
register the maps there, or in the case of STVDS, you could also add a 
different layer per date to the same vector map and register those as 
in: 
https://grass.osgeo.org/grass80/manuals/t.register.html#north-carolina-dataset 
(see the "Register a vector map with layers representing time steps" 
example)


HTH,
Vero

El vie, 22 jul 2022 a las 14:34, Ivan Marchesini 
() escribió:


Dear GRASS users,

I have a question about the creation of a spatio-temporal vector
datasets

I have a vector points map representing sample sites and it
includes a
field of the attribute table describing the sampling date (a
different
day, not sequential, for each point).

I would like to transform the vector point map into a stvds in
order to
use it with t.vect.what.strds.

My question is: is it possible to register the vector point map as
a stvds?

thank you

Ivan




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

[GRASS-user] stvds

[Ivan Marchesini]

Dear GRASS users,

I have a question about the creation of a spatio-temporal vector datasets

I have a vector points map representing sample sites and it includes a 
field of the attribute table describing the sampling date (a different 
day, not sequential, for each point).


I would like to transform the vector point map into a stvds in order to 
use it with t.vect.what.strds.


My question is: is it possible to register the vector point map as a stvds?

thank you

Ivan




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

[GRASS-user] running a simple python script on windows

[Ivan Marchesini]

Dear Grass friends,

I prepared a simple python3 script that runs correctly using Grass 7.8 
on GNU/Linux.


I put the module in the script folder of /usr/lib/grass78.  The parser 
is able to open the interface and the module runs.


Someone asked me to run the same script using Grass on Windows (the 
standalone version 7.8).


I put the script on the  "C:\Program Files\GRASS GIS 7.8\scripts" folder.

When I run a the script,  even only with the --help option 
(mySimpleGrassScript --help) from the Layer Manager Console, i get the 
error listed at the end of this e-mail.


May I kindly ask your help to understand what I have to check or to 
modify in order to run the code on window ?


thank you very much

Ivan


The error I get:

__

Traceback (most recent call last):
  File "C:\Program Files\GRASS GIS
7.8\etc\python\grass\script\task.py", line 474, in
get_interface_description

stderr=PIPE)
  File "C:\Program Files\GRASS GIS
7.8\etc\python\grass\script\core.py", line 57, in __init__

.format(args[0]))
OSError
:
Cannot find the executable mySimpleGrassScript
During handling of the above exception, another exception
occurred:
Traceback (most recent call last):
  File "C:\Program Files\GRASS GIS
7.8\gui\wxpython\gui_core\forms.py", line 2866, in
ParseCommand

blackList=_blackList)
  File "C:\Program Files\GRASS GIS
7.8\etc\python\grass\script\task.py", line 522, in
parse_interface

tree = etree.fromstring(get_interface_description(name))
  File "C:\Program Files\GRASS GIS
7.8\etc\python\grass\script\task.py", line 502, in
get_interface_description

"\n\nDetails: <{det}>").format(cmd=cmd, det=e))
grass.exceptions
.
ScriptError
:
Unable to fetch interface description for command
''.
Details: 
During handling of the above exception, another exception
occurred:
Traceback (most recent call last):
  File "C:\Program Files\GRASS GIS
7.8\gui\wxpython\core\gconsole.py", line 493, in RunCmd

task = GUI(show=None).ParseCommand(command)
  File "C:\Program Files\GRASS GIS
7.8\gui\wxpython\gui_core\forms.py", line 2868, in
ParseCommand

raise gcmd.GException(e.value)
core.gcmd
.
GException
:
Unable to fetch interface description for command
''.
Details: 
During handling of the above exception, another exception
occurred:
Traceback (most recent call last):
  File "C:\Program Files\GRASS GIS
7.8\gui\wxpython\gui_core\prompt.py", line 413, in
OnKeyPressed

self._runCmd(self.GetCurLine()[0].strip())
  File "C:\Program Files\GRASS GIS
7.8\gui\wxpython\gui_core\prompt.py", line 120, in _runCmd

self.promptRunCmd.emit(cmd=cmd)
  File "C:\Program Files\GRASS GIS
7.8\etc\python\grass\pydispatch\signal.py", line 229, in
emit

dispatcher.send(signal=self, *args, **kwargs)
  File "C:\Program Files\GRASS GIS
7.8\etc\python\grass\pydispatch\dispatcher.py", line 349, in
send

**named
  File "C:\Program Files\GRASS GIS
7.8\etc\python\grass\pydispatch\robustapply.py", line 60, in
robustApply

return receiver(*arguments, **named)
  File "C:\Program Files\GRASS GIS
7.8\gui\wxpython\gui_core\goutput.py", line 120, in 

self._gconsole.RunCmd(command=cmd))
  File "C:\Program Files\GRASS GIS
7.8\gui\wxpython\core\gconsole.py", line 496, in RunCmd

message=unicode(e),
NameError
:
name 'unicode' is not defined



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

Re: [GRASS-user] In-river (hydraulic) sediment transport dynamics

[Ivan Marchesini]

Hi Rich,

I think that this project can be interesting for you

https://sourceforge.net/projects/caesar-lisflood/

Ivan


On 26/07/19 02:47, Rich Shepard wrote:

On Thu, 25 Jul 2019, stephen sefick wrote:


The army core of engineers hec-ras might be useful.


Windoze only. They don't know anything but Microsoft.

Thanks,

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

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

Re: [GRASS-user] Fwd: Global overview of GRASS HPC resources/installations ?

[Ivan Marchesini]

Hi.. Alex

I think that your experience with CRESCO can be interesting..

thank for sharing it..

Ivan


On 25/05/2018 15:56, Alex Peloso wrote:


If you think it can be useful, i could provide you with information 
about CRESCO to be included in the "Supercomputing" wiki page.




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

Re: [GRASS-user] g.gui.animation: strds and stvds in teh same frame

[Ivan Marchesini]

Thank you Anna and Nikos
so you gave me a reason to try that way.
Thank you
Ivan

Il Dom 28 Gen 2018, 01:31 Nikos Alexandris <n...@nikosalexandris.net> ha
scritto:

> Ivan Marchesini:
>
> [..]
>
> >> I've done the same thing a couple of years ago and it seems to me that
> at
> >> that time it wasn't possible to show in the same frame of
> g.gui.animation
> >> the evolution of strds and stvds..
> >>
> >> do I remember well? is it still that case?
>
> Anna:
>
> >I think you should be able to show both, but they must have the same
> >temporal type (relative, absolute).
>
> And the same size, otherwise it won't work (from what I've quickly
> tested).
>
> Nikos
>


-- 

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

[GRASS-user] g.gui.animation: strds and stvds in teh same frame

[Ivan Marchesini]

Hi guys

I'm starting to work on an animation exploiting raster and vector 
spatial temporal dataset.


I've done the same thing a couple of years ago and it seems to me that 
at that time it wasn't possible to show in the same frame of 
g.gui.animation the evolution of strds and stvds..


do I remember well? is it still that case?

best regards

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

Re: [GRASS-user] Grass: Finding depression from contour

[Ivan Marchesini]

There is a useful addon for this: r.geomorphon .
An alternative can be also
r.param.scale
Best
Ivan

Il Dom 11 Dic 2016, 13:22 Markus Neteler  ha scritto:

> On Sun, Dec 11, 2016 at 5:01 AM, Faisal Kabir 
> wrote:
> > Hi,
> > I've labeled contour (DEM) map using d.vect. Now I would like to find the
> > depression [contour with low elev value are inside the contour of higher
> > elev] automatically. I can view it in the Map display window which is
> > tedious as I've to zoom in and out for several times.
> >
> > Is there any method available in GRASS?
>
> Yes: in r.fill.dir for example.
>
> I have taken your question as motivation to add an example to the
> manual page and updated it on the GRASS GIS Web server right away:
>
> https://grass.osgeo.org/grass72/manuals/r.fill.dir.html#examples
>
> It shows how to run sink-filling and how to identify the filled areas.
>
> Please note that filling depressions is not needed nor desired when
> calculating watersheds with r.watershed in case that this is your
> scope.
>
> Best,
> Markus
>
> --
> Markus Neteler
> http://www.mundialis.de - free data with free software
> http://grass.osgeo.org
> http://courses.neteler.org/blog
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user

-- 

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

[GRASS-user] grass and numpy (garray.array). problem when saving array to grass layer

[Ivan Marchesini]

Hi all,
I'm using grass gis 7.0 on Ubuntu
I've used the python library to import a grass raster layer in numpy 
(from grass.script import array as garray)
then I would like to export a derived array:

myarray.write(mapname="tmp",overwrite=True)

but I obtain this error:

ERROR: Required parameter  not set:
(Binary raster file to be imported)
that is triggered by r.in.bin

has someone encountered the same problem and solved it?

many thanks

Ivan




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

Re: [GRASS-user] r.reclass problem

[Ivan Marchesini]

Hi Eugenio..
use r.recode instead..

ciao


On Wed, 2015-06-03 at 12:27 +0200, Eugenio Trumpy wrote:
 
 hello everybody, 
 
 I have a problem with r.reclass.
 I would like to use such command by using a text file to reclassify a
 raster grid. 
 The text file is:
 
 -5 thru -3.5 = 1
 -3.49 thru -2.5 = 2
 -2.49 thru -1.0 = 3
 -0.99 thru 0 = 4
 0.01 thru 10 = 5
 
 However I got problem, because pixel should be in class 3 follow in
 class 4, I don't understand why? am I wronging the text file??
 
 Can somebody help me?
 
 best
 
 Eugenio 
 ___
 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] Create a raster starting from a coordinate pair

[Ivan Marchesini]

Hi Andrea,
may be you need something like that:
http://grass.osgeo.org/grass65/manuals/r.in.xyz.html

I'm not sure if that was the problem.

all the best

Ivan


On Tue, 2013-12-10 at 05:01 -0800, aborruso wrote:
 Hi Johannes,
 thank you.
 
 Johannes Radinger wrote
  What do you mean exactly with creating a raster from a coordinate pair?
 
 You are rigth, I have to be more understandable.
 
 I have defined region, with a defined extent and a defined resolution.
 Starting from a coordinate pair I would like to create a boolean raster
 (with the same extent and resolution of my region) that register my pair
 (in example 1 for my pixel and 0 for all the rest).
 
 Best regards,
 
 Andrea
 
 
 
 
 --
 View this message in context: 
 http://osgeo-org.1560.x6.nabble.com/Create-a-raster-starting-from-a-coordinate-pair-tp5093513p5093523.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

Re: [GRASS-user] segmented operations in r.drain and r.los

[Ivan Marchesini]

Hi 

 - the same question applies to r.los - can I calculate several visible
 areas with a single command?

r.los allows to define the area that is visible from one point. 
if you need to analyze more than one point you have to redo the analysis
(may be using a simple for cicle in bash)

Regards

Ivan
 


 
 -- 
 ---
 SOCIEDADE DE HISTÓRIA NATURAL - Departamento de Informação Geográfica
 
 Apartado 25 2564-909 Torres Vedras Portugal
 Sede e Biblioteca: rua Cavaleiros da Espora Dourada, 27A 2560 Torres
 Vedras
 
 Laboratório de Paleontologia e Paleoecologia: Polígono Industrial do
 Alto do Ameal 2565-641 Ramalhal
 http://alt-shn.blogspot.com
 www.alt-shn.org
 ___
 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.neighbors velocity

[Ivan Marchesini]

Dear Soeren, Hamish, Markus M., Markus N,.
sorry for the delay in the answer but I didn't believe that my question
could have determined so many answers and I was on holiday for 10 days
with no way to test your code.
First of all a quick answer to Markus N.
My collegue is working on high resolution images and DEMs. In this case,
in particular, he is working on a 1 meter resolution raster map
concerning landslides. He said me that, due to the landslides
dimensions, he need to use a kernel of 501*501 in order to catch the
signatures of the phenomena (I don't know exactly the details.. I'm
sorry).

I'm not a C developer but reading your e-mails it seems that the
performances of the C codes (and r.neigbors is written in C) strongly
depend on the compiler.
does it mean that compiling in a different way (I don't know how) the
r.neigbours module we can obtain better results?

We have tested the last code of Soeren in the same machine where the
proprietary software (it is ENVI 5.x) showed those good performances.  

These are the results:

gcc -Wall -fopenmp -lgomp -Ofast main.c -o neighbor
export OMP_NUM_THREADS=1
time ./neighbor 5000 5000 23
real0m16.598s
user0m16.477s
sys0m0.080s

export OMP_NUM_THREADS=2
time ./neighbor 5000 5000 23
real0m8.977s
user0m17.573s
sys0m0.080s

export OMP_NUM_THREADS=4
time ./neighbor 5000 5000 23
real0m5.993s
user0m20.277s
sys0m0.088s

export OMP_NUM_THREADS=6
time ./neighbor 5000 5000 23
real0m4.784s
user0m25.770s
sys0m0.096s


Many thanks for your answers
 
 
 


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

Re: [GRASS-user] r.neighbors velocity

[Ivan Marchesini]

Hi Markus
you are perfectly right

the region is 4312*5576
the moving window 501
GRASS is the stable version on a machine with 8 core and 32 gb RAM.
Ubuntu 12.04

it seems that the proprietary software is able to perform the analysis
in 2/3 seconds

:-|

ciao


On Fri, 2013-06-28 at 00:02 +0200, Markus Neteler wrote:
 On Thu, Jun 27, 2013 at 9:01 AM, Ivan Marchesini
 ivan.marches...@gmail.com wrote:
  Hi all,
  A friend of mine (having skill also in some proprietary remote sensing
  softwares) is testing GRASS for some task.
  He is quite happy about the results but he was really surprised by the
  time interval that r.neighbor take for doing the analysis. The time is
  really large in his opinion
 
 Please post some indications: computational region size and
 moving window size, also which hardware/operating system.
 Otherwise it is hard to say anything...
 
 Markus


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

[GRASS-user] r.neighbors velocity

[Ivan Marchesini]

Hi all,
A friend of mine (having skill also in some proprietary remote sensing
softwares) is testing GRASS for some task.
He is quite happy about the results but he was really surprised by the
time interval that r.neighbor take for doing the analysis. The time is
really large in his opinion. Do you feel the same in your experience? is
r.mfilter faster? do you have other ideas to perform this kind of kernel
based analyses using GRASS or other OS software?

many thanks 

Ivan   



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

Re: [GRASS-user] r.in.wms and PCN webservices

[Ivan Marchesini]

Many thanks stephan.. I will try with the new release.. Ciao

inviato da smartphone
Il giorno 31/mag/2013 18:52, Štěpán Turek stepan.tu...@seznam.cz ha
scritto:

 Hi Ivan and pierluigi,

 Hy pierluigi
 a web GET request is made of a url and, after a ?, of a set of params.
 I think that r.in.wms try to put a ? at the and of the url.
 but the wms address you use already contain a ?
 so you have to use, for the url, only the http addrees and, for the
 urlparams, the remaining part of the address.


 it worked as you described, however GetCapabilities request ignored
 urlparams parameter.  After r56519 it should be able also to use url with
 ?. GetCapabilities request should also work both ways now (r56520) (using
 urlparams or including parameters directly into url).


 Best

 Stepan



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

Re: [GRASS-user] r.in.wms and PCN webservices

[Ivan Marchesini]

Hi pierluigi

this should work in grass7 (don't know in grass 6.4) :

r.in.wms --overwrite output=tile1 url=http://wms.pcn.minambiente.it/ogc
urlparams=map=/ms_ogc/WMS_v1.3/raster/ortofoto_colore_06.map
layers=OI.ORTOIMMAGINI.2006. srs=32633 format=jpeg maxcols=2000
maxrows=2000 method=near wms_version=1.3.0

please use urlparams 

ciao

ivan

On Fri, 2013-05-31 at 14:11 +0200, Ing. Pierluigi De Rosa wrote:
 Dear GRASS User
 did anyone tried to use r.in.wms to get data from pcn[1] webservices?
 
 I want to download aerial photo from this URL:
 http://wms.pcn.minambiente.it/ogc?map=/ms_ogc/WMS_v1.3/raster/ortofoto_colore_06.map
 
 if I run
 r.in.wms -c
 url=http://wms.pcn.minambiente.it/ogc?map=/ms_ogc/WMS_v1.3/raster/ortofoto_colore_06.map
  wms_version=1.3.0
 
 it retries me an error.
 below details.
 
 bWarning/b:  ms_newMapObj(): [MapServer Error]: msLoadMap():
 MS_DEFAULT_MAPFILE_PATTERN validation failed.
  in bC:\ms4w\Apache\htdocs\wms_capabilities_v130.php/b on line
 b72/bbr /
 br /
 bWarning/b:  ms_newMapObj(): [MapServer Error]: msEvalRegex():
 String failed expression test.
 
  in bC:\ms4w\Apache\htdocs\wms_capabilities_v130.php/b on line
 b72/bbr /
 br /
 bFatal error/b:  Uncaught exception 'MapScriptException' with
 message 'Failed to open map file
 quot;/ms_ogc/WMS_v1.3/raster/ortofoto_colore_06.map?SERVICE=WMSquot;'
 in C:\ms4w\Apache\htdocs\wms_capabilities_v130.php:72
 
 Stack trace:
 #0 C:\ms4w\Apache\htdocs\wms_capabilities_v130.php(72):
 ms_newMapObj('/ms_ogc/WMS_v1')
 #1 {main}
   thrown in bC:\ms4w\Apache\htdocs\wms_capabilities_v130.php/b on
 line b72/bbr /
 
 
 [1] http://www.pcn.minambiente.it/PCNDYN/catalogowms.jsp


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

Re: [GRASS-user] Parallel GRASS jobs .....

[Ivan Marchesini]

  hi. Even If we do not use it at the end, you could try r.tileset to
create region tiles.

Ivan

inviato da smartphone
Il giorno 28/nov/2012 10:47, Andranik Hayrapetyan andranik@gmail.com
ha scritto:

 Good day.

 In grass wiki page about Parallel GRASS jobs (
 http://grasswiki.osgeo.org/wiki/Parallel_GRASS_jobs), they suggest
 2 approaches. One of them (which I want to implement and try) is

- *split map into spatial chunks (possibly with overlap to gain smooth
results)*

 But they say nothing about how you can split the map into chunks .
 Is there anybody who tried this kind of parallelism?
 How I can split the map into chunks to do the calculation in every
 chunk separately?

 Thanks in advance.

 ___
 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] [Fwd: Fwd: OGRS - GRASS list]

[ivan marchesini]

I apologize for cross-posting:

___
Subject: OGRS 2012 :: last days to register ...

There's still a little time to register. Think of it !

Discover the final program, including the workshop
descriptions and the discussion topics details.

Best regards,
OGRS 2012 program committee

# # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
# # # # #
Open Source Geospatial Research  Education Symposium
Yverdon-les-Bains, Switzerland, 24 to 26 October 2012
Go to www.ogrs2012.org and follow us on twitter.com/OGRS2012
# # # # # # # # # # # # # # # # # # # # # # # # # # # # # #
# # # # # 

_

Dr. Marj Tonini

Responsable de recherche

Bureau 3550
Tél +41 216923537

Centre de Recherche en Environnement Terrestre - CRET
Faculté des Géosciences et de l’Environnement
Quartier UNIL – Mouline, Batiment Géopolis
1015 Lausanne - Switzerland 

__


-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: ivan.marches...@irpi.cnr.it
ivan.marches...@gmail.com
fax (mailfax): +39 1782092534
jabber: geoiva...@jabber.org
skype: geoivan73


signature.asc
Description: This is a digitally signed message part
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] importing or reading a small part of a tiff file

[ivan marchesini]

Hi Hamish

Il giorno ven, 15/06/2012 alle 16.09 -0700, Hamish ha scritto:
 which leads me to an alternate solution to your problem-

:-)
many thanks...
it was a egg of Columbus  
http://en.wikipedia.org/wiki/Egg_of_Columbus

many thanks

Ivan


-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: ivan.marches...@irpi.cnr.it
ivan.marches...@gmail.com
fax (mailfax): +39 1782092534
jabber: geoiva...@jabber.org
skype: geoivan73


signature.asc
Description: This is a digitally signed message part
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] importing or reading a small part of a tiff file

[ivan marchesini]

Hi,
is there a simple way to import only a small part of a tiff file (by
mean, for example, of a bbox)?
some gdal parameter for r.in.gdal? 
many thanks..

Ivan


-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: ivan.marches...@irpi.cnr.it
ivan.marches...@gmail.com
fax (mailfax): +39 1782092534
jabber: geoiva...@jabber.org
skype: geoivan73



signature.asc
Description: This is a digitally signed message part
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

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

[ivan marchesini]

Dear Markus,
thank you for your answer.

 Please report if the stream segments are as expected, because the
 determination of half-basins and basins is based on the stream
 segments.

this is the map of streams:
http://imageshack.us/photo/my-images/232/streams.png/

it seems fine..


 If the stream segments seem out of place, something is wrong with the
 identification of stream segments and you will need to look at flow
 accumulation and check if that makes sense.

this is the accumulation map:
http://imageshack.us/photo/my-images/37/accumulation.png/

the same as above...


 If the stream segments seem to be ok, something is wrong with
 determining the left and right side of a stream.

this is the map of basins:
http://imageshack.us/photo/my-images/11/basins.png/

these are the half_basins:
http://imageshack.us/photo/my-images/4/halfbasins.png/


As you can see the problem arises for the half_basins and not only in
one place. There are a lot of problems all around.

I cannot explain which is the reason of this behaviour..

if you like I can send you the DEM


ciao

Ivan











-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: ivan.marches...@irpi.cnr.it
ivan.marches...@gmail.com
fax (mailfax): +39 1782092534
jabber: geoiva...@jabber.org
skype: geoivan73


signature.asc
Description: This is a digitally signed message part
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

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

[ivan marchesini]

Hi Markus

 Should be fixed now in all branches.

Thank you very much..   :-)

I will test it in few minuts..

I'm using grass 7 trunk..

Ivan


-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: ivan.marches...@irpi.cnr.it
ivan.marches...@gmail.com
fax (mailfax): +39 1782092534
jabber: geoiva...@jabber.org
skype: geoivan73


signature.asc
Description: This is a digitally signed message part
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

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

[ivan marchesini]

It works perfectly.. !!  :-)
thank you very much

Ivan



Il giorno mer, 09/05/2012 alle 12.23 +0200, Markus Metz ha scritto:
 On Wed, May 9, 2012 at 9:30 AM, Markus Metz
 markus.metz.gisw...@googlemail.com wrote:
  On Wed, May 9, 2012 at 9:19 AM, ivan marchesini
  ivan.marches...@gmail.com wrote:
  Dear Markus,
  thank you for your answer.
 
  Please report if the stream segments are as expected, because the
  determination of half-basins and basins is based on the stream
  segments.
 
  this is the map of streams:
  http://imageshack.us/photo/my-images/232/streams.png/
 
  it seems fine..
 
 
  If the stream segments seem out of place, something is wrong with the
  identification of stream segments and you will need to look at flow
  accumulation and check if that makes sense.
 
  this is the accumulation map:
  http://imageshack.us/photo/my-images/37/accumulation.png/
 
  the same as above...
 
 
  If the stream segments seem to be ok, something is wrong with
  determining the left and right side of a stream.
 
  this is the map of basins:
  http://imageshack.us/photo/my-images/11/basins.png/
 
  these are the half_basins:
  http://imageshack.us/photo/my-images/4/halfbasins.png/
 
 
  As you can see the problem arises for the half_basins and not only in
  one place. There are a lot of problems all around.
 
  OK, now I have an idea where to look for the problem.
 
 Should be fixed now in all branches.
 
 Markus M

-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: ivan.marches...@irpi.cnr.it
ivan.marches...@gmail.com
fax (mailfax): +39 1782092534
jabber: geoiva...@jabber.org
skype: geoivan73


signature.asc
Description: This is a digitally signed message part
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] r.watershed MFD strange basin

[ivan marchesini]

Dear Grass users,
I have noted a strange result from r.watershed using the MDF option
(with all levels of convergence).

This is a part of my DEM:

http://imageshack.us/photo/my-images/690/demv.png/


This is the result obtained using SDF:

http://imageshack.us/photo/my-images/843/27420800.png/


This is the result obtained using MDF:

http://imageshack.us/photo/my-images/190/36629290.png/



Please have a look to the last picture.
In the central part of the image (near to the pointer) there is a little
basin connected to the larger one that is on the other side of the
valley.

It seems an error.

did someone else encounter a similar problem?

many thanks

Ivan


-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: ivan.marches...@irpi.cnr.it
ivan.marches...@gmail.com
fax (mailfax): +39 1782092534
jabber: geoiva...@jabber.org
skype: geoivan73


signature.asc
Description: This is a digitally signed message part
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] How one can find the starting and end point of a line in a vector file or how one can connect two line end to end using C code

[ivan marchesini]

don't know if can help but v.to.db give you coordinates of starting and
ending points of vector lines

to connect them I use v.edit

this is my pythin code

startc=grass.read_command(v.to.db, flags=p, map=myline,
type=line, option=start, units=meters , quiet=True)
startx=float(startc.split(|)[1])
starty=float(startc.split(|)[2])
startz=float(startc.split(|)[3])
endc=grass.read_command(v.to.db, flags=p, map=myline,
type=line, option=end, units=meters)  
endx=float(endc.split(|)[1])
endy=float(endc.split(|)[2])
endz=float(endc.split(|)[3])

tmpfile = grass.tempfile()
pfile = open(tmpfile, 'w')
print pfile, L, 2, 1
print pfile, startx, starty, startz
print pfile, endx, endy, endz
print pfile, 1, i
pfile.close()
grass.run_command(v.edit, input=tmpfile, flags=n, tool=add,
map=mynewline, snap=node, thresh=1)


hope this help

Ivan


Il giorno sab, 03/03/2012 alle 11.24 +0530, SWAPAN GHOSH ha scritto:
 Hello,
  
 Any one can solve me the problem like- How one can find the starting
 and end point of a line in a vector file or how one can connect two
 line end to end using C code. If there is any change needed in code
 like C or  Python Script please suggest me. How I can do it.
  
 Regards,
  
 Swapan
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user

-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: ivan.marches...@irpi.cnr.it
ivan.marches...@gmail.com
fax (mailfax): +39 1782092534
jabber: geoiva...@jabber.org
skype: geoivan73


signature.asc
Description: This is a digitally signed message part
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] g.parser inside g.parser

[ivan marchesini]

Hi all,
We just created a script that creates another script.
I'll try to explain:


suppose we have a big text file (init.file) organized like that
___
parameter1  parameter2
0.3 0.9

parameter3  parameter4
10001

.
___

but really more complex 


we have just created a script (script1.sh) like this
(pseudocode)




#!/bin/bash
#%Module
#% description: description
#% keywords: blah
#%End
#%option
#% key: initfile
#% gisprompt: string
#% type: string
#% description: initfile
#% required : yes
#%end
if [ $1 != @ARGS_PARSED@ ] ; then
exec g.parser $0 $@
fi
- read the passed init.file
- parse the existing init.file and assign each parameter value to a
different variable
- echo 
#!/bin/bash
#%Module
#% description: description
#% keywords: raster
#%End
#%option
#% key: param1
#% gisprompt: string
#% type: string
#% description: DEM
#% required : yes
#% answer: first_variable_parsed_from_init.file
#%end


if [ $1 != @ARGS_PARSED@ ] ; then
exec g.parser $0 $@
fi

./another_fortran_code param1=something param2=something .
  script2.sh
./script2.sh
_


if we execute the script1 from command line

./script1.sh initfile=init.file

a new gui windows is opened from the execution of the script2.sh 


but when we simply run 

./script1.sh

the first gui window starts, we fill the form concerning the initfile
key, and then we click run, the second gui window (that produced by
script2.sh) doesn't start.

It seems a problem related to g.parser but we can't understand how to
solve it...

any help is really apreciated...

Thanks...

Ivan

-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: ivan.marches...@irpi.cnr.it
ivan.marches...@gmail.com
fax (mailfax): +39 1782092534
jabber: geoiva...@jabber.org
skype: geoivan73


signature.asc
Description: This is a digitally signed message part
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] strange result of v.overlay in 6.4

[ivan marchesini]

Hi Markus

 which GRASS version on which OS?
 Latest and greatest stable (i.e. 6.4?)

the release was 6.4 stable


 please consider to put a screenshot somewhere. Even better, the
 data or a reproducible example.

yes...
you are right..
I can't show my data but I can try to generate some synthetic data..

ciao

Ivan

 
 Markus

-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: ivan.marches...@irpi.cnr.it
ivan.marches...@gmail.com
fax (mailfax): +39 1782092534
jabber: geoiva...@jabber.org
skype: geoivan73


signature.asc
Description: This is a digitally signed message part
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] fuzzy

[ivan marchesini]

Dear Grass Gis Users,

I'm looking for applications, papers, modules, concerning the fuzzy
logic applied by means of grass gis (or Grass+R)

can someone give me some address, link, resource ???

many thanks...

Ivan



-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: marches...@unipg.it
ivan.marches...@gmail.com
fax (home): +39(0)5782830887
jabber: geoiva...@jabber.org
skype: geoivan73


signature.asc
Description: Questa è una parte del messaggio	firmata digitalmente
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] GRASS 7.0 svn location/map creation using console

[ivan marchesini]

Hi,
please try starting grass with:

grass -text

cheers

Ivan



Il giorno gio, 30/07/2009 alle 08.39 +0200, Enrico Gallo ha scritto:
 2009/7/30 Jachym Cepicky jachym.cepi...@gmail.com:
  Start GRASS somehow and run
 
  g.mapset -c mapset=new
 
  j
 
 hi Jachym , hi list
 
 AFAIK, g.mapset -c is OK for new mapset in existing location, but
 what about new location?
 
 We have to use ssh -X and gui?
 
 Enrico  Gallo
 
  On Wed, Jul 29, 2009 at 09:25:18PM +0200, Enrico Gallo wrote:
  Hi list,
 
  I am running remotely GRASS 7.0 from svn on a workstation using ssh.
 
  As classical location/mapset prompt is not displayed in 7.0,
  how can I create new locations/mapset using console?
 
  Thanks in advance,
 
  regards,
 
  Enrico Gallo
  ___
  grass-user mailing list
  grass-user@lists.osgeo.org
  http://lists.osgeo.org/mailman/listinfo/grass-user
 
  --
  Jachym Cepicky
  e-mail: jachym.cepicky gmail com
  URL: http://les-ejk.cz
  GPG: http://www.les-ejk.cz/pgp/JachymCepicky.pgp
  Key fingerprint: 0C6D 0EAE 76BD 506C F299  ED8A C8AB 74B8 08D4 E08F
 
  -BEGIN PGP SIGNATURE-
  Version: GnuPG v1.4.9 (GNU/Linux)
 
  iEYEARECAAYFAkpxJW0ACgkQyKt0uAjU4I9zswCfUBvjDl3SMLyF8YDwA2bcVb9/
  z/oAni+vns+9bUoJQEd4iqYe0RypgWlQ
  =zSM8
  -END PGP SIGNATURE-
 
 
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user
 
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: marches...@unipg.it
ivan.marches...@gmail.com
fax (home): +39(0)5782830887
jabber: geoiva...@jabber.org
skype: geoivan73


signature.asc
Description: Questa è una parte del messaggio	firmata digitalmente
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Building topologically correct stream network

[ivan marchesini]

Hi..

If I well remember, the problems of 0 values in the category data were
related to braiding or similar artifacts deriving from rasterization
process

r.strahler used with threshold not too small should clean network and
solve the task...

 
here [1] (in italian.. sorry...) you can find some pictures about this
topic...

sorry but it was more than one year ago and I don't remember well all
these problems..



[1]
http://www.grassmeeting2008.unipg.it/interventi/Cencetti_et_al_licensed.pdf







Il giorno mer, 22/07/2009 alle 19.35 -0700, Hamish ha scritto:
 stephen wrote:
  Well, I am not sure what I am going to do.  I would look at the
  code, but I am not adept at C.
 
 none the less, it doesn't hurt to look, depending on how it is written
 it might be understandable / tweakable.
 
 
 MMetz:
   What exactly is wrong with topology? The existence of these dangles?
 
 note the v.strahler help page makes special mention that if you used
 r.watershed you should check the topology. Maybe it makes a bitter deal
 out of that than is really the case?
 
 S:
  Maybe there is nothing wrong.  On further reflection the area that
  I am trying to use this procedure on is very flat and there are
  intermittent streams and wetlands (site visit) all over the
  study area.  These non-connected lines maybe showing
  intermittentancy?
 
 do you see a ladder effect in the r.to.vect river network?
 (like a braided river with all joins at right angles)
 maybe that is the problematic topology the module is concerned about?
 
 
 Hamish
 
 
 
   
 
 ___
 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] v.trees3d

[ivan marchesini]

Hi all, hi Jachym,

I'm looking for v.trees3d but it seems that the link which points to it
is brocken...

do someone can help me sending the code or showing another working link?

many thanks

Ivan


-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: marches...@unipg.it
ivan.marches...@gmail.com
tel: +39(0)755853760
fax (home): +39(0)5782830887
jabber: geoiva...@jabber.org
skype: geoivan73

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

Re: [GRASS-user] Plot 3D Drillholes in Grass GIS

[ivan marchesini]

Hi...
we have worked on doing volumetric interpolation of dip and direction
obtained from survey of geological strata (but if you have a lot of data
coming from more drill holes the problem is exactly the same:

1) we consider the normal to the plane representing in each point the
badding strata

2) we interpolate (3d interpolation with v.vol.rst) the three components
(xyz) of this normal, obtaining a voxel map for each component

3) we combine back this three components to obtain a volumetric
interpolated map (voxel) of geological dip and direction

4) now we can obtaine slices of this two maps  

but if your problem is simply visualization of the data coming from a
single hole then a possible solution is to create inclined planes
(small) that pass from the point where bedding has been collected.
Then represents this planes into NVIZ

cheers

ivan



Il giorno mar, 10/03/2009 alle 12.30 -0700, Casey Vandenberg ha scritto:
 I was wondering if anybody is familiar with ways of plotting drill hole 
 data in Grass (i.e. to visualize in NVIZ or to manage attributes, query 
 etc...) based on point features with known dip, strike and azimuth angles.
 
 I have come across a paper (Development of Visualization Tools for 
 Geologic Information on Grass GIS) which seems to do just that, however 
 as far as I can tell it has not been released or provided to the Grass 
 Community.
 
 If anybody is familiar with anything else, it would be greatly appreciated.
 
 Thanks,
 
 Casey
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user
 
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: marches...@unipg.it
ivan.marches...@gmail.com
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: geoiva...@jabber.org

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

Re: [GRASS-user] v.generalize for area boundaries?

[ivan marchesini]

Hi..
please try r.neighbors

cheers..
Ivan


Il giorno mar, 17/02/2009 alle 18.00 +0100, FAROUX STEPHANIE ha scritto:
 Hi,
 Is it possible with grass to replace missing values in a map by 
 surrounding values by a sort of interpolation, knowing that sometimes 
 many missing values are numerous one near another? If yes, how?
 Thank you very much
 Stéphanie
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user
 
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: marches...@unipg.it
ivan.marches...@gmail.com
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: geoiva...@jabber.org

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

Re: [GRASS-user] raster map of the azimuth to a point

[ivan marchesini]

thanks Hamish..
I'll take a look to your suggestions...
:-)

Ivan


Il giorno ven, 06/02/2009 alle 17.04 -0800, Hamish ha scritto:
 ivan marchesini wrote:
  I have a single point vector map.
  
  I would like to obtain a raster map where each cell represents the value
  of the azimuth to the point raster map.
  I mean I would like to know which direction someone have to look to see
  the point in the vector map.
  
  any suggestion?
 
 many years ago I locally modified r.los to output azimuth angle instead of
 elevation angle. If you search the mailing list archives you might find a
 link to the code.
 
 screenshot:
  http://bambi.otago.ac.nz/hamish/pics/milf_we_los-exp2.png
 
 r.cva or r.viewshed addon modules might help too.
 
 
 Hamish
 
 
 
   
 
 
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: marches...@unipg.it
ivan.marches...@gmail.com
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: geoiva...@jabber.org

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

Re: [GRASS-user] raster map of the azimuth to a point

[ivan marchesini]

This is my solution at the moment:


given a point of coordinates x,y

_
r.mapcalc px=(x()-x)
r.mapcalc py=(y()-y)

r.mapcalc azimuth=if(px=0  py=0,90-atan(py/px),if(px=0 
py0,90-atan(py/px),if(px0  py0,270-atan(py/px),if(px0 
py=0,270-atan(py/px)   
_


cheers

Ivan



Il giorno ven, 06/02/2009 alle 14.27 +0100, ivan marchesini ha scritto:
 Dear Grass Users,
 I have a single point vector map.
 
 I would like to obtain a raster map where each cell represents the value
 of the azimuth to the point raster map.
 I mean I would like to know which direction someone have to look to see
 the point in the vector map.
 
 any suggestion?
 
 many thanks 
 
 Ivan
 
 
 
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: marches...@unipg.it
ivan.marches...@gmail.com
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: geoiva...@jabber.org

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

[GRASS-user] raster map of the azimuth to a point

[ivan marchesini]

Dear Grass Users,
I have a single point vector map.

I would like to obtain a raster map where each cell represents the value
of the azimuth to the point raster map.
I mean I would like to know which direction someone have to look to see
the point in the vector map.

any suggestion?

many thanks 

Ivan



-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: marches...@unipg.it
ivan.marches...@gmail.com
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: geoiva...@jabber.org

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

[GRASS-user] v.vol.rst problem

[ivan marchesini]

Dear grass users..
I'm working with grass6 devel.
Trying to use v.vol.rst with option cellout I have obtained a strange
result:

The elev 3D map is correctly created (I have seen it by means of nviz)
but the cellout map is created as a small map (like a miniature of the
map that I would expect) that is placed at the up-left corner of the
region.. (I'm using monitors and I haven't MASK in my mapset). The
remaining part of the cellout map has only values = 0.

I think I can obtain the same thing using r3.cross.rast starting from
the elev map created from v.vol.rst, but: could the problem I have
found be a little bug? or it is a problem o f zmult or wmult (but the
resolution is the same along the x,y,z directions...) 

thanks...

Ivan




Il giorno ven, 09/01/2009 alle 12.41 -0800, Dylan Beaudette ha scritto:
 On Wednesday 07 January 2009, Helena Mitasova wrote:
  On Jan 7, 2009, at 11:19 AM, Dylan Beaudette wrote:
   On Tue, Jan 6, 2009 at 2:33 PM, Markus Neteler nete...@osgeo.org
  
   wrote:
   Dylan,
  
   On Tue, Jan 6, 2009 at 10:01 PM, Dylan Beaudette
  
   debeaude...@ucdavis.edu wrote:
   Hi,
  
   For some crazy reason I was under the impression that it is
   possible to do
   interpolation with a covariate with v.vol.rst. Are there any
   examples on how
   to parameterize this module, when a 2D surface is requested,
   rather than a 3D
   volume. I noticed the 'cellinp' argument for a cross-section, but
   this is not
   quite what I am after. I am looking to do something very similar to
   interpolation of rainfall data, taking into account the
   orographic effect of
   terrain.
  
   This was my main business (say, of our cluster) over the last
   months :)
   You can do that. I am using the elevation model as auxiliary
   variable:
  
   # something like this:
   v.vol.rst in=vectpoints cellinp=dem wcolumn=pointval cellout=rst2d
  
   cellout delivers the 2D map, extracted from the volume along the
   dem map.
  
   Hope this helps
   Markus
  
   Thanks Markus. One more question: have you found a good compromise in
   the 3D region settings- i.e. some ratio of horizontal:vertical
   resolution that gives good results and doesn't take too long to
   compute?
 
  you need to set your vertical resolution based on the spatial
  variability you want to capture
  and then don't forget to use appropriate zmult that will stretch the
  vertical distances so that
  they are about the same magnitude as horizontal distances (e.g. if
  your horizontal res. is 10m
  and vertical is 0.1m you need zmult around 100, if you want
  anisotropic effects it may be less -
  see more on how to find the parameters here:
  http://skagit.meas.ncsu.edu/~helena/gmslab/papers/
  TGIS2002_Hofierka_et_al.pdf
 
  I hope this helps, Helena
 
   Cheers,
  
   Dylan
 
 Thanks for the tips Helena. I'll check on that paper, and report back with 
 the 
 results.
 
 Dylan
 
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: marches...@unipg.it
ivan.marches...@gmail.com
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: geoiva...@jabber.org

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

Re: [GRASS-user] shell script problem

[ivan marchesini]

hmmm 
are you starting the script inside grass?

probably not... or I'm wrong?

Ivan



Il giorno gio, 04/12/2008 alle 17.06 +, Mario Giacomello ha scritto:
 Dear all,
 
 Apologies for the total beginner's question, but despite my efforts
 and searches through the GRASS mailing list archives I have not been
 able to solve my problem. I am trying to write a simple grass script
 in order to automated a process, but I have not very far yet. One
 thing that I would like to do is to be able to change the region, but
 GRASS seems to come up with the whenever I , for example, try to use
 the following simple lines:
 
 #!/bin/sh
 
 g.region rast=map1
 r.stats -ln map1  file1
 
 The message I get is: command not found. What am I doing wrong? Is
 this problem related to the setting of environment variables?
 
 I am using GRASS 6.3.0  installed on Ubuntu 8.10
 
 thanks in advance,
 
 Mario Giacomello
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user
 
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

Re: [GRASS-user] r.terraflow

[ivan marchesini]

Il giorno mar, 29/07/2008 alle 02.22 +0200, G. Allegri ha scritto:
 Thanks Luca.
 The coding mapping seems correct. So, zeroes are for unsolved pixels...
 I'll do some tests to compare the results with r.watershed. At now
 they seem quite different. I'm not sure if I should trust onflow
 directions from r.terraflow to use them in r.water.outlet...
 

h in my opinion the r.terraflow output is not usefull for
r.water.outlet but if you find a solution please let we know   ;-)

 bye

Ivan




 http://lists.osgeo.org/mailman/listinfo/grass-user
 
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

[GRASS-user] r.watershed speed-up

[ivan marchesini]

First results of one test I did: excellent!!!

1392776 cells dem with steep slopes but also plain areas

r.watershed.fast: 7 seconds
r.watershed:  4 min and 14 seconds

no difference at all between the two generated network paths...

markus: many thanks..
 

Ivan


Il giorno lun, 28/07/2008 alle 23.14 +0200, Markus Metz ha scritto:
 Hello list,
 
 I'm not sure if this list is the right place or rather the developer list.
 For the A * Search algorithm in r.watershed (memory version without -m 
 flag set, r.watershed.ram), I implemented a binary min-heap instead of a 
 linear array to store points in ascending order relative to elevation. 
 This improves the speed of SECTION 2: A * Search considerably, the 
 larger the map, the larger the speed improvement. A region with about 
 1,800,000 cells is processed about 37 times faster than with the 
 original routine. A region with about 11,000,000 cells is processed 
 about 170 times faster than with the original routine (46.7sec vs. 
 2h12m9sec on my system). SECTION 4: Watershed Determination takes now 
 the longest.
 
 Several tests (different region sizes, different DEMs) showed that the 
 results of the binary heap version are very similar but not 100% 
 identical to the original r.watershed, because of a slightly different 
 treatment of cells with equal elevation in a binary min-heap compared to 
 a linear array. Differences are found in difficult areas (equal 
 elevation of several neighbouring cells) where there are several 
 possible solutions for how water could flow. Still, compared with other 
 hydrology analysis methods, e.g. r.terraflow, the results can be 
 regarded as identical.
 
 The original code was only altered where necessary, only the sorting 
 method is new, everything else is unchanged.
 Memory usage increases a bit, because a binary heap needs its own index 
 for each analysed point (in addition to the other indices already needed 
 by r.watershed.ram).
 
 My question is if there is interest in this faster version of 
 r.watershed and if someone wants to test it.
 The source code is available on http://markus.metz.giswork.googlepages.com/
 
 Regards,
 
 Markus
 
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user
 
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

Re: [GRASS-user] v.net option connect

[ivan marchesini]

oh
yes...

I supposed it was a stupid question  :-)

Thanks... Martin...

ciao



Il giorno mer, 23/07/2008 alle 16.05 +0200, Martin Landa ha scritto:
 Ciao Ivan,
 
 2008/7/23 ivan marchesini [EMAIL PROTECTED]:
  I have a network of 6379 roads
  I have a map with a single point
 
  I run:
 
  v.net input=road points=point operation=connect output=connected_point
  thresh=(very high thresh) alayer=1 nlayer=2
 
  then (to check the output):
 
  v.build connected_point
 
  and I see that the number of lines is 6381
 
  i.e. the module v.net creates two lines for one point???
 
 the network arc(edge) need to be splitted, e.g.
 
 --
 
   +
 
 -
 
 +-
|
|
   +
 
 In this case you get two lines more.
 
 Cheers, Martin
 
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

[GRASS-user] v.net option connect

[ivan marchesini]

may be it is a stupid question but I can't understand this problem:

I have a network of 6379 roads
I have a map with a single point

I run:

v.net input=road points=point operation=connect output=connected_point
thresh=(very high thresh) alayer=1 nlayer=2

then (to check the output):

v.build connected_point

and I see that the number of lines is 6381

i.e. the module v.net creates two lines for one point???

many thanks...

Ivan






-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

Re: [GRASS-user] new user setting up a new mapset/project

[ivan marchesini]

hi

 So when I set up the project, which coord should I use?

You must create two location: the first one UTM WGS84 and the second one
long/lat WGS84

you import your long/lat data into the proper location and then you can
use v.proj or r.proj to reproject them into the UTM WGS84 location

another solution is to use gdalwarp or ogr2ogr but could be hard if
you are just starting to learn foss4g softwares..

 
 When I went to use UTM, I did not understand how to get the north,  
 south, east, west info. Grass did not like the tiff, and associated  
 tfw for the base map.

grass imports tiff files...
 you must use the --o option (override projection) to import data with
r.in.gdall



 using Ubuntu, Grass 6.3
 

hope this helps...


 
 Thanks,
 Tom


Ivan


 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user
 
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

Re: [GRASS-user] Stereo image tools

[ivan marchesini]

I Terry...
take a look to efoto project
http://www.efoto.eng.uerj.br/efoto

ciao





Il giorno gio, 19/06/2008 alle 15.29 +1000, Terry Duell ha scritto:
 Hullo All,
 I would like to try to generate some 3d data from stereo image pairs, and
 have looked at the stereo-0.2b package at
 http://grass.itc.it/outgoing/grass5/.
 It looks as though it was aimed at problems that are similar to what I have
 in mind, but unfortunately it fails with errors when I try to build it, and I
 don't think I can solve the problems. There is a binary as well, but it also
 has problems.
 As this package dates back quite a few years someone may have done some
 further development/bug fixing, but thus far I have not been able to find any
 other references to it on the web.
 Does anyone know of any later versions of this software?
 
 Cheers,
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

[GRASS-user] anisotropy and 3d interpolation

[ivan marchesini]

Dear all...
I have this kind of data:

two parallel lines of tomographic seismic profiles, that means (for who
is not a geologist), that I have, along these lines, regularly spaced 3d
points (like a vertical grid of points under the earth surface) each one
having as an attribute a value of the P waves velocity (P waves are
longitudinal seismic waves, high velocity means presence of rocks, low
velocity of sand or clay).
The house for which I did this research is just in the middle between
these two lines (see the next map). I would like to obtain a 3d velocity
model of the terrain just under the house. 


^ North axis
|
|
|    Tomograpphic profile1
| _
||_|  house
|    Tomograpphic profile2
|
|
|---  East axis 



It is simple to obtain a 3d voxel map of the area under the house
(v.vol.rst) but it is not simple to obtain an anisotropic
interpolation because v.vol.rst doesn't have possibility to define
an ellipsoid of the anisotropy as it is possible to do with v.surf.rst.

Since, for geological reasons, it is clear that there is velocity
continuity (isotropy) along the est-nord plane but not along the Z
plane, I would like to know if someone can suggest me a way to solve
this problem in grass (3d interpolation with anisotropy), or in case,
using some others free-GPL softwares (R?, or there's something else??). 

many thanks for the time you losed reading this post...  :-)

Ivan






-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

Re: [GRASS-user] selecting upstream river segments

[ivan marchesini]

Hi Markus and Silvia...
many thanks for your answer..

I have tried to use the solution proposed by Markus..
it works but it selects all the connencted stream... instead and I would
like to obtain only those upstream to the point I choose...
removing the others downstream

I have thougth this solution:
suppose I'm interested to the upstream part of a stream network starting
from a choosen segment of the stream network, so I can:
1) use v.to.db to find the downstream node of the selected segment
2) use v.edit to move this node (disconneting it from the downstream
segment)
3) get the new coordinates of the node
4) use the procedure suggested by markus (v.net.iso) starting from the
moved downstream node coordinates

What do you think about this solution?

many thanks 

Ivan






Il giorno mar, 22/04/2008 alle 10.08 +0200, Markus Neteler ha scritto:
 Hi Ivan,
 
 please see
 http://www.nabble.com/Streams-extraction-from-streams-map-td10267212.html#a10306465
 
 Markus
 
 On Mon, Apr 21, 2008 at 8:58 AM, ivan marchesini [EMAIL PROTECTED] wrote:
  Dear Grass users,
 
   I have a vector stream network map of a single basin (so only a stream
   tree)...
 
   each stream segment is correctly oriented downstream (the direction of
   each vector is from higher to lower elevation..)
 
   I would like to select a random segment of the tree and automatically
   obtain only the upstream segments..
   i.e. I would like to obtain a new map containing only the upstream
   stream segments of the selected one... (that is the same of removing
   those downstream)
 
   obviously I must implement this process into a shell script...
 
   Could you suggest a simple way to solve this little problem?
   Can I use the direction of the vectors to solve this problem?
 
   Many thanks
 
   Ivan
 
 
 
 
   --
   Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
   Preferisco formati liberi.
   Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
   I prefer free formats.
   http://it.wikipedia.org/wiki/Formato_aperto
   http://en.wikipedia.org/wiki/Open_format
 
   Ivan Marchesini
   Department of Civil and Environmental Engineering
   University of Perugia
   Via G. Duranti 93/a
   06125
   Perugia (Italy)
   Socio fondatore GFOSS Geospatial Free and Open Source Software 
  http://www.gfoss.it
   e-mail: [EMAIL PROTECTED]
  [EMAIL PROTECTED]
   tel: +39(0)755853760
   fax (university): +39(0)755853756
   fax (home): +39(0)5782830887
   jabber: [EMAIL PROTECTED]
 
   ___
   grass-user mailing list
   grass-user@lists.osgeo.org
   http://lists.osgeo.org/mailman/listinfo/grass-user
 
 
 
 
-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

[GRASS-user] selecting upstream river segments

[ivan marchesini]

Dear Grass users,

I have a vector stream network map of a single basin (so only a stream
tree)...

each stream segment is correctly oriented downstream (the direction of
each vector is from higher to lower elevation..)

I would like to select a random segment of the tree and automatically
obtain only the upstream segments.. 
i.e. I would like to obtain a new map containing only the upstream
stream segments of the selected one... (that is the same of removing
those downstream)

obviously I must implement this process into a shell script...

Could you suggest a simple way to solve this little problem? 
Can I use the direction of the vectors to solve this problem?

Many thanks 

Ivan




-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

[GRASS-user] r.watershed massive dataset

[ivan marchesini]

Dear friends,
I'm trying to use r.watershed to obtain dainage and stream map of a very
large dem..
The region containing the dem is of around 260.000.000 of cells but
masking out null cells (with a raster MASK) the number of residuals dem
cells is of around 140.000.000.
I have tried without the -m flag but the process stops after few minutes
due to memory allocation problems (the computer has only 2 GB of RAM). 
Now I have started the process with the -m flag

Do you have any suggestion to optimize the process?

thanks 

Ivan



-- 
Ti prego di cercare di non inviarmi files .dwg, .doc, .xls, .ppt.
Preferisco formati liberi.
Please try to avoid to send me .dwg, .doc, .xls, .ppt files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

[GRASS-user] v.kernel

[ivan marchesini]

Dear users...
I'm trying to use v.kernel...
but I didn't well understand the meaning of the raster density map
created...
do you know what does the values of the raster map mean?
are they an estimation of the number of points per unit width? (in my
case the number of point per square meter?)

many thanks...

Ivan



-- 
Ti prego di cercare di non inviarmi files .doc, .xls, .ppt, .dwg.
Preferisco formati liberi.
Please try to avoid to send me  .doc, .xls, .ppt, .dwg files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

Re: [GRASS-user] Snap across layers

[ivan marchesini]

Hi Martin and Luca...
I have just tried with v.out.ascii as you suggested...

I have two test vector layer (pippo and pippo1)
each one contains only one area (without centroid)
they are really similar and vertex are quite near each other (~30
meters)

I have tried these commands:

v.edit a3 tool=create
v.out.ascii pippo1 format=standard | v.edit a3 tool=add bgmap=pippo
thre=200 snap=vertex

it works but, at the end I obtain only a copy of the first map (i.e.
a3 is exactly equal to pippo1)
It seems that only add works... but not the snap and the threshold...

any suggestion?

many thanks

Ivan






Il giorno mar, 29/01/2008 alle 15.52 +0100, Martin Landa ha scritto:
 Hi,
 
 2008/1/28, Luca Casagrande [EMAIL PROTECTED]:
  Markus Neteler ha scritto:
   - v.edit does it:
 bgmap   Name of background vector map(s) snap to
   - the new digitizer in wxgrass does it (not yet enabled by default)
 
 Luca:
  Before moving to my data, I am doing some tests with 2 vector layers:
 
  http://img337.imageshack.us/my.php?image=schermatamapdisplay1vp5.png
 
  Basically the small one, should snap to the bigger one.
  The command line I use is this:
 
  v.edit map=origine bgmap=destinazione cat=1-999 tool=snap thresh=-1
 
  But nothing change.
 
  In the log I have got a lot of line like this, both if i set thresh -1
  or 500
  [..]
  D1/1: Vect_set_thresh(): thresh = 0.00
  [..]
 
 
 tool 'snap' in v.edit uses library fn Vect_snap_lines() which doesn't
 support snapping to the vector features from background map layer.
 
 Snapping mode (no, node, vertex) and background map option works
 currently only for tool='add'. For now, try to create an empty vector
 map, export to ascii format the map you want to modify and update the
 empty map, e.g.
 
 v.edit a3 tool=create
 v.out.ascii a2 format=standard | v.edit a3 tool=add bgmap=a1 id=1-999
 thre=500 snap=node
 
 Martin
 
-- 
Ti prego di cercare di non inviarmi files .doc, .xls, .ppt, .dwg.
Preferisco formati liberi.
Please try to avoid to send me  .doc, .xls, .ppt, .dwg files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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

Re: [GRASS-user] raster module for travel time calculations (updated)

[ivan marchesini]

Thank you very much for your work Kristian...
really...
Cheers 
Ivan

Il giorno dom, 30/12/2007 alle 09.44 +0100, [EMAIL PROTECTED] ha
scritto:
 Dear GRASS developers and users,
 
 I am using GRASS GIS for some hydrological applications. For some questions 
 like flood prediction it is useful to estimate the travel time of surface 
 water in a basin to the outlet. Therefor I've developed a raster program 
 r.traveltime which calculates the travel time for each cell in a basin from 
 outputs of r.watershed and r.fill.dir. There is much literature about this 
 procedure available but there's a lag of available software (GRASS and even 
 ArcGIS do not come with this function). The personal need an the lag of 
 available programs were the reason for me to develop a new program for GRASS 
 GIS. I am not a programming expert, so the program is very simple and far 
 away from being perfect.
 I submitted this tool to the GRASS Wiki addon page.
 
 In November I have already posted r.traveltime to the GRASS weblist mailing 
 list. Recently I've found some bugs in the program. I hope it works correctly 
 now.
 
 r.traveltime is available here:
 http://home.arcor.de/jesbergwetter/main_sw.html
 
 Best regards
 
 Kristian Förster
 
 
 Viel oder wenig? Schnell oder langsam? Unbegrenzt surfen + telefonieren
 ohne Zeit- und Volumenbegrenzung? DAS TOP ANGEBOT FÜR ALLE NEUEINSTEIGER
 Jetzt bei Arcor: günstig und schnell mit DSL - das All-Inclusive-Paket
 für clevere Doppel-Sparer, nur  29,95 Euro  inkl. DSL- und ISDN-Grundgebühr!
 http://www.arcor.de/rd/emf-dsl-2
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user
 
-- 
Ti prego di cercare di non inviarmi files .doc, .xls, .ppt, .dwg.
Preferisco formati liberi.
Please try to avoid to send me  .doc, .xls, .ppt, .dwg files.
I prefer free formats.
http://it.wikipedia.org/wiki/Formato_aperto
http://en.wikipedia.org/wiki/Open_format

Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
Socio fondatore GFOSS Geospatial Free and Open Source Software 
http://www.gfoss.it
e-mail: [EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: +39(0)755853760
fax (university): +39(0)755853756
fax (home): +39(0)5782830887
jabber: [EMAIL PROTECTED]

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