[GRASS-user] Surfaces with faults

GRASS GRASS user
1

Hello
I would like to use the procedure described in "Open Source GIS: A
GRASS GIS Approach" second edition page 191 "7.3.5 Surfaces with
faults".
I have a survey of points. Some points are two ditches. I have run a
v.surf.rst with the survey points without the ditch.
Attached is the ascii file of the ditch points (EPSG:32632)
I have therefore operated as follows:

v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
separator=comma
g.region vector=fosso_linea res=0.2
v.to.rast -d --overwrite input=fosso_linea type=line output=fosso /
use=z
v.surf.rst --overwrite segmax=60 npmin=100 dmin=1 dmax=6 /
input=fosso_linea elevation=fosso_dem mask=fosso

the result is with no_data

Some suggestions?
--
--
Perito agrario Enrico Gabrielli
progetto F.A.R.M. www.farm-agroecologia.it
Tessera n. 633 Collegio Periti agrari prov. Di Modena
Biblioteca agricoltura: https://www.zotero.org/groups/aplomb/
https://www.inaturalist.org/observations/bonushenricus

(attachments)

fosso.dat (1.62 KB)

2

Hello
I've just tried with the solution of this thread
https://www.mail-archive.com/grass-user@lists.osgeo.org/msg40640.html
so

v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
separator=comma
g.region vector=fosso_linea res=0.2
v.split input=fosso_linea output=fosso_linea_split length=0.2
v.to.points -i input=fosso_linea_split output=fosso_punti
v.to.rast -d input=fosso_punti type=point output=fosso_dem use=z

but the result is

Command ended with non-zero return code -11 (0 sec)

I believe there was something wrong in v.in.lines or in my ASCII file:
the vector line hasn't attribute

Il giorno lun, 24/07/2023 alle 16.44 +0200, bonushenricus ha scritto:

Hello
I would like to use the procedure described in "Open Source GIS: A
GRASS GIS Approach" second edition page 191 "7.3.5 Surfaces with
faults".
I have a survey of points. Some points are two ditches. I have run a
v.surf.rst with the survey points without the ditch.
Attached is the ascii file of the ditch points (EPSG:32632)
I have therefore operated as follows:
> v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
> separator=comma
> g.region vector=fosso_linea res=0.2
> v.to.rast -d --overwrite input=fosso_linea type=line output=fosso /
> use=z
> v.surf.rst --overwrite segmax=60 npmin=100 dmin=1 dmax=6 /
> input=fosso_linea elevation=fosso_dem mask=fosso
the result is with no_data

Some suggestions?
--

3

Ok
v.to.points creates 2 layer.
The good layer is the 2.
So with
v.to.rast -d input=fosso_punti type=point layer=2 output=fosso_dem
use=z
is correct.
I prefered the way with v.surf.rst from my original survey points, but
the way with v.to.points can help me.

Il giorno lun, 24/07/2023 alle 17.59 +0200, bonushenricus ha scritto:

Hello
I've just tried with the solution of this thread
https://www.mail-archive.com/grass-user@lists.osgeo.org/msg40640.html
so
> v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
> separator=comma
> g.region vector=fosso_linea res=0.2
> v.split input=fosso_linea output=fosso_linea_split length=0.2
> v.to.points -i input=fosso_linea_split output=fosso_punti
> v.to.rast -d input=fosso_punti type=point output=fosso_dem use=z
but the result is
> Command ended with non-zero return code -11 (0 sec)

I believe there was something wrong in v.in.lines or in my ASCII
file:
the vector line hasn't attribute

Il giorno lun, 24/07/2023 alle 16.44 +0200, bonushenricus ha scritto:
> Hello
> I would like to use the procedure described in "Open Source GIS: A
> GRASS GIS Approach" second edition page 191 "7.3.5 Surfaces with
> faults".
> I have a survey of points. Some points are two ditches. I have run
> a
> v.surf.rst with the survey points without the ditch.
> Attached is the ascii file of the ditch points (EPSG:32632)
> I have therefore operated as follows:
> > v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
> > separator=comma
> > g.region vector=fosso_linea res=0.2
> > v.to.rast -d --overwrite input=fosso_linea type=line output=fosso
> > /
> > use=z
> > v.surf.rst --overwrite segmax=60 npmin=100 dmin=1 dmax=6 /
> > input=fosso_linea elevation=fosso_dem mask=fosso
> the result is with no_data
>
> Some suggestions?
> --

4

Hello
the problem isn't in my ASCII file or in v.in.lines
but a bug in v.to.rast, with the -d parameter.
I used gdb

Thread 1 "v.to.rast" received signal SIGSEGV, Segmentation fault.
(gdb) bt
#0 0x0000000000000000 in ?? ()
#1 0x00007ffff7ef7091 in G_plot_point ()

   from /usr/local/grass83/lib/libgrass_gis.8.3.so

#2 0x0000555555558d46 in do_lines ()
#3 0x000055555555abc6 in vect_to_rast ()
#4 0x0000555555557ad6 in main ()
(gdb)

Il giorno lun, 24/07/2023 alle 17.59 +0200, bonushenricus ha scritto:

Hello
I've just tried with the solution of this thread
https://www.mail-archive.com/grass-user@lists.osgeo.org/msg40640.html
so
> v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
> separator=comma
> g.region vector=fosso_linea res=0.2
> v.split input=fosso_linea output=fosso_linea_split length=0.2
> v.to.points -i input=fosso_linea_split output=fosso_punti
> v.to.rast -d input=fosso_punti type=point output=fosso_dem use=z
but the result is
> Command ended with non-zero return code -11 (0 sec)

I believe there was something wrong in v.in.lines or in my ASCII
file:
the vector line hasn't attribute

Il giorno lun, 24/07/2023 alle 16.44 +0200, bonushenricus ha scritto:
> Hello
> I would like to use the procedure described in "Open Source GIS: A
> GRASS GIS Approach" second edition page 191 "7.3.5 Surfaces with
> faults".
> I have a survey of points. Some points are two ditches. I have run
> a
> v.surf.rst with the survey points without the ditch.
> Attached is the ascii file of the ditch points (EPSG:32632)
> I have therefore operated as follows:
> > v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
> > separator=comma
> > g.region vector=fosso_linea res=0.2
> > v.to.rast -d --overwrite input=fosso_linea type=line output=fosso
> > /
> > use=z
> > v.surf.rst --overwrite segmax=60 npmin=100 dmin=1 dmax=6 /
> > input=fosso_linea elevation=fosso_dem mask=fosso
> the result is with no_data
>
> Some suggestions?
> --

5

Please provide your input data and computational region settings. I
was not able to reproduce the issue with my data.

Māris.

otrd., 2023. g. 25. jūl., plkst. 07:17 — lietotājs bonushenricus
(<bonushenricus76@gmail.com>) rakstīja:

Hello
the problem isn't in my ASCII file or in v.in.lines
but a bug in v.to.rast, with the -d parameter.
I used gdb
> Thread 1 "v.to.rast" received signal SIGSEGV, Segmentation fault.
> (gdb) bt
> #0 0x0000000000000000 in ?? ()
> #1 0x00007ffff7ef7091 in G_plot_point ()
   from /usr/local/grass83/lib/libgrass_gis.8.3.so
> #2 0x0000555555558d46 in do_lines ()
> #3 0x000055555555abc6 in vect_to_rast ()
> #4 0x0000555555557ad6 in main ()
> (gdb)

Il giorno lun, 24/07/2023 alle 17.59 +0200, bonushenricus ha scritto:
> Hello
> I've just tried with the solution of this thread
> https://www.mail-archive.com/grass-user@lists.osgeo.org/msg40640.html
> so
> > v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
> > separator=comma
> > g.region vector=fosso_linea res=0.2
> > v.split input=fosso_linea output=fosso_linea_split length=0.2
> > v.to.points -i input=fosso_linea_split output=fosso_punti
> > v.to.rast -d input=fosso_punti type=point output=fosso_dem use=z
> but the result is
> > Command ended with non-zero return code -11 (0 sec)
>
> I believe there was something wrong in v.in.lines or in my ASCII
> file:
> the vector line hasn't attribute
>
> Il giorno lun, 24/07/2023 alle 16.44 +0200, bonushenricus ha scritto:
> > Hello
> > I would like to use the procedure described in "Open Source GIS: A
> > GRASS GIS Approach" second edition page 191 "7.3.5 Surfaces with
> > faults".
> > I have a survey of points. Some points are two ditches. I have run
> > a
> > v.surf.rst with the survey points without the ditch.
> > Attached is the ascii file of the ditch points (EPSG:32632)
> > I have therefore operated as follows:
> > > v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
> > > separator=comma
> > > g.region vector=fosso_linea res=0.2
> > > v.to.rast -d --overwrite input=fosso_linea type=line output=fosso
> > > /
> > > use=z
> > > v.surf.rst --overwrite segmax=60 npmin=100 dmin=1 dmax=6 /
> > > input=fosso_linea elevation=fosso_dem mask=fosso
> > the result is with no_data
> >
> > Some suggestions?
> > --
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

6

Thank you.
the CRS is EPSG:32632
fosso.dat I attached

<v.in.lines -z --overwrite input=fosso.dat output=fosso_linea

separator=comma

g.region vector=fosso_linea res=0.2
v.split input=fosso_linea output=fosso_linea_split length=0.1
v.to.points -i input=fosso_linea_split output=fosso_punti
v.to.rast -d input=fosso_punti type=point output=fosso_dem

use=z

Il giorno mar, 25/07/2023 alle 08.33 +0300, Maris Nartiss ha scritto:

Please provide your input data and computational region settings. I
was not able to reproduce the issue with my data.

Māris.

otrd., 2023. g. 25. jūl., plkst. 07:17 — lietotājs bonushenricus
(<bonushenricus76@gmail.com>) rakstīja:
>
> Hello
> the problem isn't in my ASCII file or in v.in.lines
> but a bug in v.to.rast, with the -d parameter.
> I used gdb
> > Thread 1 "v.to.rast" received signal SIGSEGV, Segmentation fault.
> > (gdb) bt
> > #0 0x0000000000000000 in ?? ()
> > #1 0x00007ffff7ef7091 in G_plot_point ()
> from /usr/local/grass83/lib/libgrass_gis.8.3.so
> > #2 0x0000555555558d46 in do_lines ()
> > #3 0x000055555555abc6 in vect_to_rast ()
> > #4 0x0000555555557ad6 in main ()
> > (gdb)
>
>
>
> Il giorno lun, 24/07/2023 alle 17.59 +0200, bonushenricus ha
> scritto:
> > Hello
> > I've just tried with the solution of this thread
> > https://www.mail-archive.com/grass-user@lists.osgeo.org/msg40640.html
> > so
> > > v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
> > > separator=comma
> > > g.region vector=fosso_linea res=0.2
> > > v.split input=fosso_linea output=fosso_linea_split length=0.2
> > > v.to.points -i input=fosso_linea_split output=fosso_punti
> > > v.to.rast -d input=fosso_punti type=point output=fosso_dem
> > > use=z
> > but the result is
> > > Command ended with non-zero return code -11 (0 sec)
> >
> > I believe there was something wrong in v.in.lines or in my ASCII
> > file:
> > the vector line hasn't attribute
> >
> > Il giorno lun, 24/07/2023 alle 16.44 +0200, bonushenricus ha
> > scritto:
> > > Hello
> > > I would like to use the procedure described in "Open Source
> > > GIS: A
> > > GRASS GIS Approach" second edition page 191 "7.3.5 Surfaces
> > > with
> > > faults".
> > > I have a survey of points. Some points are two ditches. I have
> > > run
> > > a
> > > v.surf.rst with the survey points without the ditch.
> > > Attached is the ascii file of the ditch points (EPSG:32632)
> > > I have therefore operated as follows:
> > > > v.in.lines -z --overwrite input=fosso.dat output=fosso_linea
> > > > /
> > > > separator=comma
> > > > g.region vector=fosso_linea res=0.2
> > > > v.to.rast -d --overwrite input=fosso_linea type=line
> > > > output=fosso
> > > > /
> > > > use=z
> > > > v.surf.rst --overwrite segmax=60 npmin=100 dmin=1 dmax=6 /
> > > > input=fosso_linea elevation=fosso_dem mask=fosso
> > > the result is with no_data
> > >
> > > Some suggestions?
> > > --
> _______________________________________________
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user

(attachments)

fosso.dat (1.62 KB)

7

Almost eureka
spline with tension only for point of the ditch
is possible with r.resample.rst

r.resamp.rst --overwrite input=fosso_punti_rast@granello ew_res=0.1

ns_res=0.1 elevation=fosso_dem_rst@granello overlap=2 tension=40.
but the result in width of the line is too big
and also the profile is with "terraces"

Il giorno mar, 25/07/2023 alle 08.33 +0300, Maris Nartiss ha scritto:

Please provide your input data and computational region settings. I
was not able to reproduce the issue with my data.

Māris.

otrd., 2023. g. 25. jūl., plkst. 07:17 — lietotājs bonushenricus
(<bonushenricus76@gmail.com>) rakstīja:
>
> Hello
> the problem isn't in my ASCII file or in v.in.lines
> but a bug in v.to.rast, with the -d parameter.
> I used gdb
> > Thread 1 "v.to.rast" received signal SIGSEGV, Segmentation fault.
> > (gdb) bt
> > #0 0x0000000000000000 in ?? ()
> > #1 0x00007ffff7ef7091 in G_plot_point ()
> from /usr/local/grass83/lib/libgrass_gis.8.3.so
> > #2 0x0000555555558d46 in do_lines ()
> > #3 0x000055555555abc6 in vect_to_rast ()
> > #4 0x0000555555557ad6 in main ()
> > (gdb)
>
>
>
> Il giorno lun, 24/07/2023 alle 17.59 +0200, bonushenricus ha
> scritto:
> > Hello
> > I've just tried with the solution of this thread
> > https://www.mail-archive.com/grass-user@lists.osgeo.org/msg40640.html
> > so
> > > v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
> > > separator=comma
> > > g.region vector=fosso_linea res=0.2
> > > v.split input=fosso_linea output=fosso_linea_split length=0.2
> > > v.to.points -i input=fosso_linea_split output=fosso_punti
> > > v.to.rast -d input=fosso_punti type=point output=fosso_dem
> > > use=z
> > but the result is
> > > Command ended with non-zero return code -11 (0 sec)
> >
> > I believe there was something wrong in v.in.lines or in my ASCII
> > file:
> > the vector line hasn't attribute
> >
> > Il giorno lun, 24/07/2023 alle 16.44 +0200, bonushenricus ha
> > scritto:
> > > Hello
> > > I would like to use the procedure described in "Open Source
> > > GIS: A
> > > GRASS GIS Approach" second edition page 191 "7.3.5 Surfaces
> > > with
> > > faults".
> > > I have a survey of points. Some points are two ditches. I have
> > > run
> > > a
> > > v.surf.rst with the survey points without the ditch.
> > > Attached is the ascii file of the ditch points (EPSG:32632)
> > > I have therefore operated as follows:
> > > > v.in.lines -z --overwrite input=fosso.dat output=fosso_linea
> > > > /
> > > > separator=comma
> > > > g.region vector=fosso_linea res=0.2
> > > > v.to.rast -d --overwrite input=fosso_linea type=line
> > > > output=fosso
> > > > /
> > > > use=z
> > > > v.surf.rst --overwrite segmax=60 npmin=100 dmin=1 dmax=6 /
> > > > input=fosso_linea elevation=fosso_dem mask=fosso
> > > the result is with no_data
> > >
> > > Some suggestions?
> > > --
> _______________________________________________
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user

8

And at the end I found the solutions:
after a v.to.rast at the desired resolution
I can choose the line width with r.buffer
and at the end I can use v.surf.rst, with raster line mask

v.surf.rst --overwrite input=fosso_punti zcolumn=z

elevation=fosso_dem_rst mask=line3d_rast smooth_column=vdop segmax=60
dmin=0.2 dmax=5
!!!

Il giorno mar, 25/07/2023 alle 06.17 +0200, bonushenricus ha scritto:

Hello
the problem isn't in my ASCII file or in v.in.lines
but a bug in v.to.rast, with the -d parameter.
I used gdb
> Thread 1 "v.to.rast" received signal SIGSEGV, Segmentation fault.
> (gdb) bt
> #0 0x0000000000000000 in ?? ()
> #1 0x00007ffff7ef7091 in G_plot_point ()
from /usr/local/grass83/lib/libgrass_gis.8.3.so
> #2 0x0000555555558d46 in do_lines ()
> #3 0x000055555555abc6 in vect_to_rast ()
> #4 0x0000555555557ad6 in main ()
> (gdb)

Il giorno lun, 24/07/2023 alle 17.59 +0200, bonushenricus ha scritto:
> Hello
> I've just tried with the solution of this thread
> https://www.mail-archive.com/grass-user@lists.osgeo.org/msg40640.html
> so
> > v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
> > separator=comma
> > g.region vector=fosso_linea res=0.2
> > v.split input=fosso_linea output=fosso_linea_split length=0.2
> > v.to.points -i input=fosso_linea_split output=fosso_punti
> > v.to.rast -d input=fosso_punti type=point output=fosso_dem use=z
> but the result is
> > Command ended with non-zero return code -11 (0 sec)
>
> I believe there was something wrong in v.in.lines or in my ASCII
> file:
> the vector line hasn't attribute
>
> Il giorno lun, 24/07/2023 alle 16.44 +0200, bonushenricus ha
> scritto:
> > Hello
> > I would like to use the procedure described in "Open Source GIS:
> > A
> > GRASS GIS Approach" second edition page 191 "7.3.5 Surfaces with
> > faults".
> > I have a survey of points. Some points are two ditches. I have
> > run
> > a
> > v.surf.rst with the survey points without the ditch.
> > Attached is the ascii file of the ditch points (EPSG:32632)
> > I have therefore operated as follows:
> > > v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
> > > separator=comma
> > > g.region vector=fosso_linea res=0.2
> > > v.to.rast -d --overwrite input=fosso_linea type=line
> > > output=fosso
> > > /
> > > use=z
> > > v.surf.rst --overwrite segmax=60 npmin=100 dmin=1 dmax=6 /
> > > input=fosso_linea elevation=fosso_dem mask=fosso
> > the result is with no_data
> >
> > Some suggestions?
> > --

9

Thanks for the data. I was able to reproduce the segfault. I opened a
bug report so it doesn't get lost.
https://github.com/OSGeo/grass/issues/3105

Māris.

otrd., 2023. g. 25. jūl., plkst. 08:41 — lietotājs bonushenricus
(<bonushenricus76@gmail.com>) rakstīja:

Thank you.
the CRS is EPSG:32632
fosso.dat I attached
><v.in.lines -z --overwrite input=fosso.dat output=fosso_linea
separator=comma
> g.region vector=fosso_linea res=0.2
> v.split input=fosso_linea output=fosso_linea_split length=0.1
> v.to.points -i input=fosso_linea_split output=fosso_punti
> v.to.rast -d input=fosso_punti type=point output=fosso_dem
use=z

Il giorno mar, 25/07/2023 alle 08.33 +0300, Maris Nartiss ha scritto:
> Please provide your input data and computational region settings. I
> was not able to reproduce the issue with my data.
>
> Māris.
>
> otrd., 2023. g. 25. jūl., plkst. 07:17 — lietotājs bonushenricus
> (<bonushenricus76@gmail.com>) rakstīja:
> >
> > Hello
> > the problem isn't in my ASCII file or in v.in.lines
> > but a bug in v.to.rast, with the -d parameter.
> > I used gdb
> > > Thread 1 "v.to.rast" received signal SIGSEGV, Segmentation fault.
> > > (gdb) bt
> > > #0 0x0000000000000000 in ?? ()
> > > #1 0x00007ffff7ef7091 in G_plot_point ()
> > from /usr/local/grass83/lib/libgrass_gis.8.3.so
> > > #2 0x0000555555558d46 in do_lines ()
> > > #3 0x000055555555abc6 in vect_to_rast ()
> > > #4 0x0000555555557ad6 in main ()
> > > (gdb)
> >
> >
> >
> > Il giorno lun, 24/07/2023 alle 17.59 +0200, bonushenricus ha
> > scritto:
> > > Hello
> > > I've just tried with the solution of this thread
> > > https://www.mail-archive.com/grass-user@lists.osgeo.org/msg40640.html
> > > so
> > > > v.in.lines -z --overwrite input=fosso.dat output=fosso_linea /
> > > > separator=comma
> > > > g.region vector=fosso_linea res=0.2
> > > > v.split input=fosso_linea output=fosso_linea_split length=0.2
> > > > v.to.points -i input=fosso_linea_split output=fosso_punti
> > > > v.to.rast -d input=fosso_punti type=point output=fosso_dem
> > > > use=z
> > > but the result is
> > > > Command ended with non-zero return code -11 (0 sec)
> > >
> > > I believe there was something wrong in v.in.lines or in my ASCII
> > > file:
> > > the vector line hasn't attribute
> > >
> > > Il giorno lun, 24/07/2023 alle 16.44 +0200, bonushenricus ha
> > > scritto:
> > > > Hello
> > > > I would like to use the procedure described in "Open Source
> > > > GIS: A
> > > > GRASS GIS Approach" second edition page 191 "7.3.5 Surfaces
> > > > with
> > > > faults".
> > > > I have a survey of points. Some points are two ditches. I have
> > > > run
> > > > a
> > > > v.surf.rst with the survey points without the ditch.
> > > > Attached is the ascii file of the ditch points (EPSG:32632)
> > > > I have therefore operated as follows:
> > > > > v.in.lines -z --overwrite input=fosso.dat output=fosso_linea
> > > > > /
> > > > > separator=comma
> > > > > g.region vector=fosso_linea res=0.2
> > > > > v.to.rast -d --overwrite input=fosso_linea type=line
> > > > > output=fosso
> > > > > /
> > > > > use=z
> > > > > v.surf.rst --overwrite segmax=60 npmin=100 dmin=1 dmax=6 /
> > > > > input=fosso_linea elevation=fosso_dem mask=fosso
> > > > the result is with no_data
> > > >
> > > > Some suggestions?
> > > > --
> > _______________________________________________
> > grass-user mailing list
> > grass-user@lists.osgeo.org
> > https://lists.osgeo.org/mailman/listinfo/grass-user