[GRASS-user] v.surf.nnbathy

GRASS GRASS user
1

Dear list,

Is there any working version of v.surf.nnbathy around? From this recent post
http://www.nabble.com/do-not-export-categories-with-v.out.ascii-t4142590.html

it did not sound like an impossible thing. I am looking to replace v.surf.idw with a more robust interpolation method (natural neighbour?) in my mapping routine -- my input is scattered 3D vector points. I am just checking to see if somebody else already did it before trying to modify r.surf.nnbathy myself (I am not really a programmer).

Thanks and regards,

Luigi

2

I haven't played with this yet, but I might do it. Give me a few days, ok?

cheers

Carlos

On 8/22/07, Luigi Ponti <lponti@infinito.it> wrote:

Dear list,

Is there any working version of v.surf.nnbathy around? From this recent post
http://www.nabble.com/do-not-export-categories-with-v.out.ascii-t4142590.html

it did not sound like an impossible thing. I am looking to replace
v.surf.idw with a more robust interpolation method (natural neighbour?)
in my mapping routine -- my input is scattered 3D vector points. I am
just checking to see if somebody else already did it before trying to
modify r.surf.nnbathy myself (I am not really a programmer).

Thanks and regards,

Luigi

_______________________________________________
grassuser mailing list
grassuser@grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

--
+-----------------------------------------------------------+
              Carlos Henrique Grohmann - Guano
  Visiting Researcher at Kingston University London - UK
  Geologist M.Sc - Doctorate Student at IGc-USP - Brazil
Linux User #89721 - carlos dot grohmann at gmail dot com
+-----------------------------------------------------------+
_________________
"Good morning, doctors. I have taken the liberty of removing Windows
95 from my hard drive."
--The winning entry in a "What were HAL's first words" contest judged
by 2001: A SPACE ODYSSEY creator Arthur C. Clarke

Can't stop the signal.

3

Luigi, I just went to look at r.surf.nnbathy, and we _can_ use it to
interpolate 3D vector points.
just don't enter any raster, and use the second option, for a 3D vector.

works well, and it is fast.

cheers

Carlos

On 8/22/07, Luigi Ponti <lponti@infinito.it> wrote:

Dear list,

Is there any working version of v.surf.nnbathy around? From this recent post
http://www.nabble.com/do-not-export-categories-with-v.out.ascii-t4142590.html

it did not sound like an impossible thing. I am looking to replace
v.surf.idw with a more robust interpolation method (natural neighbour?)
in my mapping routine -- my input is scattered 3D vector points. I am
just checking to see if somebody else already did it before trying to
modify r.surf.nnbathy myself (I am not really a programmer).

Thanks and regards,

Luigi

_______________________________________________
grassuser mailing list
grassuser@grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

--
+-----------------------------------------------------------+
              Carlos Henrique Grohmann - Guano
  Visiting Researcher at Kingston University London - UK
  Geologist M.Sc - Doctorate Student at IGc-USP - Brazil
Linux User #89721 - carlos dot grohmann at gmail dot com
+-----------------------------------------------------------+
_________________
"Good morning, doctors. I have taken the liberty of removing Windows
95 from my hard drive."
--The winning entry in a "What were HAL's first words" contest judged
by 2001: A SPACE ODYSSEY creator Arthur C. Clarke

Can't stop the signal.

4

Carlos: Thanks for your prompt reply.

I am using r.surf.nnbathy 1.9, 2007.01.02, and if I don’t enter any raster name, it complains:

ERROR: Raster map not found

How can the script work if I don’t give it an input name? I did not understand what is the “second option” you mentioned in your message, so I am asking whether the following command line corresponds to what you mean:

r.surf.nnbathy input= output=3D_vector_name alg=nn

Sorry for my bad intuition and thanks again,

Luigi

Carlos “Guâno” Grohmann wrote:

5

Luigi (and everyone else)

can you try this version?

I must be going crazy or something. Apparently, I did the changes in
r.surf.nnbathy back when I started that thread you cited, about a
month ago.. but I forgot to test it and submit it..

I did a fast test here and it worked fine, with a 3D vector points layer.

hope it works!

cheers

Carlos

On 8/22/07, Luigi Ponti <lponti@infinito.it> wrote:

Carlos: Thanks for your prompt reply.

I am using r.surf.nnbathy 1.9, 2007.01.02, and if I don't enter any raster
name, it complains:

ERROR: Raster map not found
How can the script work if I don't give it an input name? I did not
understand what is the "second option" you mentioned in your message, so I
am asking whether the following command line corresponds to what you mean:

r.surf.nnbathy input= output=3D_vector_name alg=nn
Sorry for my bad intuition and thanks again,

Luigi

Carlos "Guâno" Grohmann wrote:
Luigi, I just went to look at r.surf.nnbathy, and we _can_ use it to
interpolate 3D vector points.
just don't enter any raster, and use the second option, for a 3D vector.

works well, and it is fast.

cheers

Carlos

On 8/22/07, Luigi Ponti <lponti@infinito.it> wrote:

Dear list,

Is there any working version of v.surf.nnbathy around? From this recent post
http://www.nabble.com/do-not-export-categories-with-v.out.ascii-t4142590.html

it did not sound like an impossible thing. I am looking to replace
v.surf.idw with a more robust interpolation method (natural neighbour?)
in my mapping routine -- my input is scattered 3D vector points. I am
just checking to see if somebody else already did it before trying to
modify r.surf.nnbathy myself (I am not really a programmer).

Thanks and regards,

Luigi

_______________________________________________
grassuser mailing list
grassuser@grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

--
+-----------------------------------------------------------+
              Carlos Henrique Grohmann - Guano
  Visiting Researcher at Kingston University London - UK
  Geologist M.Sc - Doctorate Student at IGc-USP - Brazil
Linux User #89721 - carlos dot grohmann at gmail dot com
+-----------------------------------------------------------+
_________________
"Good morning, doctors. I have taken the liberty of removing Windows
95 from my hard drive."
--The winning entry in a "What were HAL's first words" contest judged
by 2001: A SPACE ODYSSEY creator Arthur C. Clarke

Can't stop the signal.

(attachments)

r.surf.nnbathy (6.45 KB)

6

I’m catching up your discuss in a hurry, but if it can help, for a couple of weeks I’ve been using r.surf.nnbathy to generate TINs :
http://www.mail-archive.com/grassuser@grass.itc.it/msg00891.html
I turn my vector 3D points into a raster, then make the interpolation (with the ‘l’ option) with r.surf.nnbathy, it works fine for me.
And sorry if I flunk !

VB

Le mercredi 22 août 2007 à 23:13 +0100, Carlos “Guâno” Grohmann a écrit :

Luigi (and everyone else)

can you try this version?

I must be going crazy or something. Apparently, I did the changes in
r.surf.nnbathy back when I started that thread you cited, about a
month ago.. but I forgot to test it and submit it..

I did a fast test here and it worked fine, with a 3D vector points layer.

hope it works!

cheers

Carlos

On 8/22/07, Luigi Ponti <[lponti@infinito.it](mailto:lponti@infinito.it)> wrote:
>
>  Carlos: Thanks for your prompt reply.
>
>  I am using r.surf.nnbathy 1.9, 2007.01.02, and if I don't enter any raster
> name, it complains:
>
> ERROR: Raster map not found
>  How can the script work if I don't give it an input name? I did not
> understand what is the "second option" you mentioned in your message, so I
> am asking whether the following command line corresponds to what you mean:
>
> r.surf.nnbathy input= output=3D_vector_name alg=nn
>  Sorry for my bad intuition and thanks again,
>
>  Luigi
>
>
>  Carlos "Guâno" Grohmann wrote:
>  Luigi, I just went to look at r.surf.nnbathy, and we _can_ use it to
> interpolate 3D vector points.
> just don't enter any raster, and use the second option, for a 3D vector.
>
> works well, and it is fast.
>
> cheers
>
> Carlos
>
>
>
> On 8/22/07, Luigi Ponti <[lponti@infinito.it](mailto:lponti@infinito.it)> wrote:
>
>
>  Dear list,
>
> Is there any working version of v.surf.nnbathy around? From this recent post
> [http://www.nabble.com/do-not-export-categories-with-v.out.ascii-t4142590.html](http://www.nabble.com/do-not-export-categories-with-v.out.ascii-t4142590.html)
>
> it did not sound like an impossible thing. I am looking to replace
> v.surf.idw with a more robust interpolation method (natural neighbour?)
> in my mapping routine -- my input is scattered 3D vector points. I am
> just checking to see if somebody else already did it before trying to
> modify r.surf.nnbathy myself (I am not really a programmer).
>
> Thanks and regards,
>
> Luigi
>
> _______________________________________________
> grassuser mailing list
> [grassuser@grass.itc.it](mailto:grassuser@grass.itc.it)
> [http://grass.itc.it/mailman/listinfo/grassuser](http://grass.itc.it/mailman/listinfo/grassuser)
>
>
>
>
>
>

_______________________________________________
grassuser mailing list
[grassuser@grass.itc.it](mailto:grassuser@grass.itc.it)
[http://grass.itc.it/mailman/listinfo/grassuser](http://grass.itc.it/mailman/listinfo/grassuser)
7

Thanks Carlos,

Your version works fine but I had to use the Hamish trick posted at http://www.nabble.com/do-not-export-categories-with-v.out.ascii-t4142590.html

v.out.ascii input=“${VECTOR}” format=point dp=3 | cut -f1-3 -d’|’ | tr ‘|’ ‘,’ > $TMP.${PROG}.input_xyz

because the fs (field separator) option is not available in the GRASS-6.2.2 version of v.out.ascii that I am using.

Kind regards,

Luigi

Carlos “Guâno” Grohmann wrote: