[GRASSLIST:6478] Sites and GSLIB

Dear list member: a couple of questions:

1- Migrating from Grass5.4 to Grass 6 I see that "sites" data object no longer exists. Personally, I found this a little drawback (clearly this could be a drawback only for me and maybe I was using the incorrect procedure: tell me !!!!!), especially when trying to importing huge amount of point cloud data (i.e. for example when using data collected by means of airborne laser scanning technology: millions of points). I found really difficult to import elevation data and execute interpolation (compared when using grass5.4 with sites): it seems that I need to link the imported vector points to a database table in order to define the column category. Then I`m wondering if using the "vector" format there is more memory consumption than when using "sites" format. Have you any suggestions?

2- I`m interested in visualizing 3D grid file generated by geostatistical tools (mainly GSLIB library: SAsim, Sisim, 3dKT....) in NViz...have you any example file of a 3d grid text file in order to import it in Grass? I would like to write a small code to convert the GSlib format to Grass format

Thank you in advance for your time

Sebastiano T.

I must say that i'm not "100%" sure of the conventions either...
but if you check
r3.out.ascii
(the manual page on the inet,
http://www.grass-kr.org/html/r3.out.ascii.html)
you will find more information on the grid specifications...

I can also send you the source code of my "gslib2g3d" program (which is
a big name for that). It's written fortran so if you're familiar with
the gslib code it should be fairly easy to understand...

(btw i've put grasslist back in the "to" field)

gr,
Samuel

On Fri, 2005-04-15 at 14:42 +0200, seba wrote:

Hi
yes, they are continuos......
Well, I installed grass on fedora...and seems to work more or less properly
(NVIZ works fine).
Could you send me an example ascii grid file to put in grass....
The guide isn`t clear about grid convenctions.
Bye
Sebastiano

At 13.57 15/04/2005, you wrote:
>Hi,
>
>I'm trying to visualize concentration so i guess that are continuous
>attributes.
>I never posted this on the list because i'm basically pretty new in
>Grass...
>I've been able to set up the region and to import my data but
>vizualisation seems to require extra effort... :slight_smile:
>(And in anyway i think there is still something lacking on my grass
>installation...
>Imo there's something wrong with the ebuild file for gentoo so i'm
>trying now a manual recompile)
>
>gr,
>
>Sam
>
>
>
>
>
>On Fri, 2005-04-15 at 11:05 +0200, seba wrote:
> > Hello Samuel
> > Good to know....
> > Are you tryng to visualizing continuos attributes (i.e. concentration,
> > elevation, ....) or discrete (categories) ?
> > Did you get any answer about this issue???
> > Bye
> > Sebastiano
> >
> > At 10.06 15/04/2005, you wrote:
> > >Hello seba...
> > >
> > >I'm actually trying *exactly* the same...
> > >I was writing a fortran program and actually it is working... I can
> > >import my 3d Gslib rasters in grass6... But i just fail to visualize
> > >them...
> > >
> > >Grass is actually reading them and he recognizes the coordinates correct
> > >but the visualization itself seems to be problematic...
> > >
> > >Samuel V.
> > >
> > >
> > >On Fri, 2005-04-15 at 09:53 +0200, seba wrote:
> > > > Dear list member: a couple of questions:
> > > >
> > > > 1- Migrating from Grass5.4 to Grass 6 I see that "sites" data object no
> > > > longer exists. Personally, I found this a little drawback (clearly this
> > > > could be a drawback only for me and maybe I was using the incorrect
> > > > procedure: tell me !!!!!), especially when trying to importing huge
> amount
> > > > of point cloud data (i.e. for example when using data collected by
> > > means of
> > > > airborne laser scanning technology: millions of points). I found really
> > > > difficult to import elevation data and execute interpolation (compared
> > > when
> > > > using grass5.4 with sites): it seems that I need to link the imported
> > > > vector points to a database table in order to define the column
> category.
> > > > Then I`m wondering if using the "vector" format there is more memory
> > > > consumption than when using "sites" format. Have you any suggestions?
> > > >
> > > > 2- I`m interested in visualizing 3D grid file generated by
> geostatistical
> > > > tools (mainly GSLIB library: SAsim, Sisim, 3dKT....) in NViz...have
> you
> > > any
> > > > example file of a 3d grid text file in order to import it in Grass? I
> > > would
> > > > like to write a small code to convert the GSlib format to Grass format
> > > >
> > > >
> > > > Thank you in advance for your time
> > > >
> > > > Sebastiano T.
> > > >
> > > >
> > > >
> >

What I wonder.... if i understand the manual correct... I should have a
command... g3.region to initialise my 3d environment... this command
seems to be absent...

On Fri, 2005-04-15 at 14:55 +0200, samuel verstraete wrote:

I must say that i'm not "100%" sure of the conventions either...
but if you check
r3.out.ascii
(the manual page on the inet,
http://www.grass-kr.org/html/r3.out.ascii.html)
you will find more information on the grid specifications...

I can also send you the source code of my "gslib2g3d" program (which is
a big name for that). It's written fortran so if you're familiar with
the gslib code it should be fairly easy to understand...

(btw i've put grasslist back in the "to" field)

gr,
Samuel

On Fri, 2005-04-15 at 14:42 +0200, seba wrote:
> Hi
> yes, they are continuos......
> Well, I installed grass on fedora...and seems to work more or less properly
> (NVIZ works fine).
> Could you send me an example ascii grid file to put in grass....
> The guide isn`t clear about grid convenctions.
> Bye
> Sebastiano
>
> At 13.57 15/04/2005, you wrote:
> >Hi,
> >
> >I'm trying to visualize concentration so i guess that are continuous
> >attributes.
> >I never posted this on the list because i'm basically pretty new in
> >Grass...
> >I've been able to set up the region and to import my data but
> >vizualisation seems to require extra effort... :slight_smile:
> >(And in anyway i think there is still something lacking on my grass
> >installation...
> >Imo there's something wrong with the ebuild file for gentoo so i'm
> >trying now a manual recompile)
> >
> >gr,
> >
> >Sam
> >
> >
> >
> >
> >
> >On Fri, 2005-04-15 at 11:05 +0200, seba wrote:
> > > Hello Samuel
> > > Good to know....
> > > Are you tryng to visualizing continuos attributes (i.e. concentration,
> > > elevation, ....) or discrete (categories) ?
> > > Did you get any answer about this issue???
> > > Bye
> > > Sebastiano
> > >
> > > At 10.06 15/04/2005, you wrote:
> > > >Hello seba...
> > > >
> > > >I'm actually trying *exactly* the same...
> > > >I was writing a fortran program and actually it is working... I can
> > > >import my 3d Gslib rasters in grass6... But i just fail to visualize
> > > >them...
> > > >
> > > >Grass is actually reading them and he recognizes the coordinates correct
> > > >but the visualization itself seems to be problematic...
> > > >
> > > >Samuel V.
> > > >
> > > >
> > > >On Fri, 2005-04-15 at 09:53 +0200, seba wrote:
> > > > > Dear list member: a couple of questions:
> > > > >
> > > > > 1- Migrating from Grass5.4 to Grass 6 I see that "sites" data object no
> > > > > longer exists. Personally, I found this a little drawback (clearly this
> > > > > could be a drawback only for me and maybe I was using the incorrect
> > > > > procedure: tell me !!!!!), especially when trying to importing huge
> > amount
> > > > > of point cloud data (i.e. for example when using data collected by
> > > > means of
> > > > > airborne laser scanning technology: millions of points). I found really
> > > > > difficult to import elevation data and execute interpolation (compared
> > > > when
> > > > > using grass5.4 with sites): it seems that I need to link the imported
> > > > > vector points to a database table in order to define the column
> > category.
> > > > > Then I`m wondering if using the "vector" format there is more memory
> > > > > consumption than when using "sites" format. Have you any suggestions?
> > > > >
> > > > > 2- I`m interested in visualizing 3D grid file generated by
> > geostatistical
> > > > > tools (mainly GSLIB library: SAsim, Sisim, 3dKT....) in NViz...have
> > you
> > > > any
> > > > > example file of a 3d grid text file in order to import it in Grass? I
> > > > would
> > > > > like to write a small code to convert the GSlib format to Grass format
> > > > >
> > > > >
> > > > > Thank you in advance for your time
> > > > >
> > > > > Sebastiano T.
> > > > >
> > > > >
> > > > >
> > >
>

Hello

Under you will find the file on which I'm experimenting on...
So it would be 2 levels and a 10*10 grid x and y

north: 10.0
south: 0.0
east: 10.0
west: 0.0
top: 4.0
bottom:0.0
rows:10
cols:10
levels:2
2.0 2.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
2.0 2.0 1.0 1.0 3.0 1.0 2.0 2.0 1.5 5.0
2.0 2.0 1.0 1.0 3.0 1.0 2.0 2.0 2.5 7.5
2.0 2.0 1.0 1.0 2.0 1.0 1.5 3.0 4.0 10.0
1.0 1.0 1.0 1.0 1.0 1.0 1.0 2.0 3.0 5.0
2.0 2.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
2.0 2.0 1.0 1.0 3.0 1.0 2.0 2.0 1.5 5.0
2.0 2.0 1.0 1.0 3.0 1.0 2.0 2.0 2.5 7.5
2.0 2.0 1.0 1.0 2.0 1.0 1.5 3.0 4.0 10.0
1.0 1.0 1.0 1.0 1.0 1.0 1.0 2.0 3.0 5.0
2.0 2.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
2.0 2.0 1.0 1.0 3.0 1.0 2.0 2.0 1.5 5.0
2.0 2.0 1.0 1.0 3.0 1.0 2.0 2.0 2.5 7.5
2.0 2.0 4.0 1.0 2.0 1.0 1.5 3.0 4.0 10.0
1.0 1.0 1.0 5.0 1.0 1.0 1.0 2.0 3.0 5.0
2.0 2.0 1.0 1.0 6.0 1.0 1.0 1.0 1.0 1.0
2.0 2.0 1.0 1.0 3.0 1.0 2.0 2.0 1.5 5.0
2.0 2.0 1.0 1.0 3.0 6.0 2.0 2.0 2.5 7.5
2.0 2.0 1.0 3.0 4.0 1.0 1.5 3.0 4.0 10.0
1.0 1.0 1.0 4.0 1.0 1.0 1.0 2.0 3.0 5.0

Import with r3.in.ascii

and run r3.info on it I get :

r3.info map=inside

+----------------------------------------------------------------------+
|Layer : inside Date: ?? |
|Mapset: samuel Login Of Creator:?? |
|Location: grass2 |
|Database: /home/samuel/grass2 |
|Title: ( ??) |
|----------------------------------------------------------------------|
| |
| Type of Map: 3d cell Number of Categories:0 |
| Rows: 10 |
| Columns: 10 |
| Depths : 2 |
| Total Cells: 200 |
| Projection: x,y (zone 0) |
| N: 10 S: 0 Res:1 |
| E: 10 W: 0 Res:1 |
| T: 4 B: 0 Res:2 |
| |
+----------------------------------------------------------------------+

So this seems pretty ok... although I'm not entirely sure about the
Number of Categories thing... Pretty weird that it would be zero...

But if i run "nviz" on it... And I go browsing for a raster... It
doesn't pop up...

gr,

Samuel

On Fri, 2005-04-15 at 15:56 +0200, samuel verstraete wrote:

What I wonder.... if i understand the manual correct... I should have a
command... g3.region to initialise my 3d environment... this command
seems to be absent...

On Fri, 2005-04-15 at 14:55 +0200, samuel verstraete wrote:
> I must say that i'm not "100%" sure of the conventions either...
> but if you check
> r3.out.ascii
> (the manual page on the inet,
> http://www.grass-kr.org/html/r3.out.ascii.html)
> you will find more information on the grid specifications...
>
> I can also send you the source code of my "gslib2g3d" program (which is
> a big name for that). It's written fortran so if you're familiar with
> the gslib code it should be fairly easy to understand...
>
> (btw i've put grasslist back in the "to" field)
>
> gr,
> Samuel
>
>
> On Fri, 2005-04-15 at 14:42 +0200, seba wrote:
> > Hi
> > yes, they are continuos......
> > Well, I installed grass on fedora...and seems to work more or less properly
> > (NVIZ works fine).
> > Could you send me an example ascii grid file to put in grass....
> > The guide isn`t clear about grid convenctions.
> > Bye
> > Sebastiano
> >
> > At 13.57 15/04/2005, you wrote:
> > >Hi,
> > >
> > >I'm trying to visualize concentration so i guess that are continuous
> > >attributes.
> > >I never posted this on the list because i'm basically pretty new in
> > >Grass...
> > >I've been able to set up the region and to import my data but
> > >vizualisation seems to require extra effort... :slight_smile:
> > >(And in anyway i think there is still something lacking on my grass
> > >installation...
> > >Imo there's something wrong with the ebuild file for gentoo so i'm
> > >trying now a manual recompile)
> > >
> > >gr,
> > >
> > >Sam
> > >
> > >
> > >
> > >
> > >
> > >On Fri, 2005-04-15 at 11:05 +0200, seba wrote:
> > > > Hello Samuel
> > > > Good to know....
> > > > Are you tryng to visualizing continuos attributes (i.e. concentration,
> > > > elevation, ....) or discrete (categories) ?
> > > > Did you get any answer about this issue???
> > > > Bye
> > > > Sebastiano
> > > >
> > > > At 10.06 15/04/2005, you wrote:
> > > > >Hello seba...
> > > > >
> > > > >I'm actually trying *exactly* the same...
> > > > >I was writing a fortran program and actually it is working... I can
> > > > >import my 3d Gslib rasters in grass6... But i just fail to visualize
> > > > >them...
> > > > >
> > > > >Grass is actually reading them and he recognizes the coordinates correct
> > > > >but the visualization itself seems to be problematic...
> > > > >
> > > > >Samuel V.
> > > > >
> > > > >
> > > > >On Fri, 2005-04-15 at 09:53 +0200, seba wrote:
> > > > > > Dear list member: a couple of questions:
> > > > > >
> > > > > > 1- Migrating from Grass5.4 to Grass 6 I see that "sites" data object no
> > > > > > longer exists. Personally, I found this a little drawback (clearly this
> > > > > > could be a drawback only for me and maybe I was using the incorrect
> > > > > > procedure: tell me !!!!!), especially when trying to importing huge
> > > amount
> > > > > > of point cloud data (i.e. for example when using data collected by
> > > > > means of
> > > > > > airborne laser scanning technology: millions of points). I found really
> > > > > > difficult to import elevation data and execute interpolation (compared
> > > > > when
> > > > > > using grass5.4 with sites): it seems that I need to link the imported
> > > > > > vector points to a database table in order to define the column
> > > category.
> > > > > > Then I`m wondering if using the "vector" format there is more memory
> > > > > > consumption than when using "sites" format. Have you any suggestions?
> > > > > >
> > > > > > 2- I`m interested in visualizing 3D grid file generated by
> > > geostatistical
> > > > > > tools (mainly GSLIB library: SAsim, Sisim, 3dKT....) in NViz...have
> > > you
> > > > > any
> > > > > > example file of a 3d grid text file in order to import it in Grass? I
> > > > > would
> > > > > > like to write a small code to convert the GSlib format to Grass format
> > > > > >
> > > > > >
> > > > > > Thank you in advance for your time
> > > > > >
> > > > > > Sebastiano T.
> > > > > >
> > > > > >
> > > > > >
> > > >
> >
>

I've found a little error in the help of v.in.ascii in the:
file:///Applications/Grass/grass61cvs.app/Contents/Resources/grass-6.1.cvs/docs/html/v.in.ascii.html

It is written:

cat points.dat | v.in.ascii out=mypoints xcol=2 ycol=3 catcol=1 columns='cat int, x double, y double, label varchar(20)'

but it is correct:

cat points.dat | v.in.ascii output=mypoints x=2 y=3 cat=1 columns='cat int, x double, y double, label varchar(20)'

I hope this is useful
Giacomo

I've found a little error in the help of v.in.ascii in the:
file: [...]/v.in.ascii.html

It is written:

cat points.dat | v.in.ascii out=mypoints xcol=2 ycol=3 catcol=1
columns='cat int, x double, y double, label varchar(20)'

but it is correct:

cat points.dat | v.in.ascii output=mypoints x=2 y=3 cat=1 columns='cat
int, x double, y double, label varchar(20)'

I hope this is useful
Giacomo

fixed in CVS, thanks for noticing.

Hamish

1- Migrating from Grass5.4 to Grass 6 I see that "sites" data object
no longer exists. Personally, I found this a little drawback (clearly
this could be a drawback only for me and maybe I was using the
incorrect procedure: tell me !!!!!),

"Sites" are now replaced with "vector points". Use the v.* modules now.
Most functionality, plus a lot more, still exists.

especially when trying to importing huge amount of point cloud data
(i.e. for example when using data collected by means of airborne
laser scanning technology: millions of points). I found really
difficult to import elevation data and execute interpolation (compared
when using grass5.4 with sites): it seems that I need to link the
imported vector points to a database table in order to define the
column category. Then I`m wondering if using the "vector" format
there is more memory consumption than when using "sites" format. Have
you any suggestions?

there was a memory leak bug in the 6.0.0 release. It has been fixed now
and I have loaded 1.15 million point LIDAR data into GRASS without
problems. You will need to compile from the latest CVS snapshot to get
this fix.

2- I`m interested in visualizing 3D grid file generated by
geostatistical tools (mainly GSLIB library: SAsim, Sisim, 3dKT....)
in NViz...have you any example file of a 3d grid text file in order
to import it in Grass? I would like to write a small code to convert
the GSlib format to Grass format

see
http://grass.ibiblio.org/grid3d/index.html
http://grass.ibiblio.org/grass60/manuals/html60_user/raster3D.html
http://grass.ibiblio.org/grass60/manuals/html60_user/r3.in.ascii.html

Hamish

What I wonder.... if i understand the manual correct... I should have
a command... g3.region to initialise my 3d environment... this command
seems to be absent...

merged into g.region

Hamish

On Fri, Apr 15, 2005 at 02:55:14PM +0200, samuel verstraete wrote:

I must say that i'm not "100%" sure of the conventions either...
but if you check
r3.out.ascii
(the manual page on the inet,
http://www.grass-kr.org/html/r3.out.ascii.html)
you will find more information on the grid specifications...

I'm not sure if above site is updated regularly. Better to
go here:
http://grass.itc.it/gdp/manuals.php
or to an official mirror site.

The specifications are described in the Programmer's Manual
http://grass.itc.it/devel/index.php#prog

Just to avoid confusion

Markus

On Fri, Apr 15, 2005 at 03:56:24PM +0200, samuel verstraete wrote:

What I wonder.... if i understand the manual correct... I should have a
command... g3.region to initialise my 3d environment... this command
seems to be absent...

It was integrated into g.region. See

http://grass.gdf-hannover.de/twiki/bin/view/GRASS/ModulePortingList

Markus