[GRASSLIST:6490] Re: Sites and GSLIB
samuel verstraete
samuel.verstraete at gmail.com
Fri Apr 15 10:25:33 EDT 2005
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... :)
> > > >(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.
> > > > > > >
> > > > > > >
> > > > > > >
> > > > >
> > >
> >
>
More information about the grass-user
mailing list