Lindsay Mico | 1 Mar 2007 01:13
Picon
Favicon

[GRASS-user] using LLID numbers with the linear referencing system

Fellow GRASS users,

I would like to build a new point data layer from a large database of stream 
surveys which have 2 pieces of spatial info:

1) an LLID number for that stream
2) distance from the start of that stream

Based on what I have seen, I would use the linear referencing system 
modules.  I am not sure if I need to build a new reference layer as I have a 
stream layer with the LLID info already in it.  Additionally, I could not 
find an example of how v.lrs.segment works in practice.  I have attempted to 
summarize things as I understand it to help make it easier to answer my 
questions.  As I understand it the input parameters for v.lrs.segment would 
be as follows:

1) input= the stream layer with the associated LLID info
2) output = the name of the output layer
3) llayer = whether it is a point or line (is it 1 for line and 0 for 
point?)
4) rsdriver = I don't know what this means
5) rsdatabase = the location of the file with the LLID # and distance
6) rstable = the name of the file which has the LLID number and distance?  
will it automatically add in the other data in the database?  Or do I do 
that later using a spreadsheet program on the .dbf

As I am still learning how to do this, I expect I will have more questions 
in the future.  Hopefully this will be of use to others as well.  Many 
thanks in advance.  Cheers.

(Continue reading)

Hamish | 1 Mar 2007 01:14
X-Face
Picon
Favicon

Re: [GRASS-user] 64 bit or parallel processing

Glynn Clements wrote:
> It would be more feasible to modify individual modules; e.g. a
> multi-threaded version of r.mapcalc might be feasible. But unless the
> raster I/O code was made thread-safe, you would still be limited by
> the rate at which map data could be read and written by a single
> thread, so it would only be worthwhile for cases where the bulk of the
> overhead was in the actual calculations.

Right, so working on individual modules which can typically take 10+
hours to run (*.rst, watersheds, viewsheds (incl. r.sun), indices from
RS data, etc) seems like a much better focus of effort for the greatest
gain:work ratio. And as seen in the MPI add-ons this is exactly where
the previous ad-hoc effort has gone.

IIRC Helena said that the RST modules have changed substantially sice
GRASS 5 and the s.surf.rst.mpi work would have to be largely redone.

Any thoughts on gains that could be made by MPI'ing the segmentation
lib? Do do modules using that usually do so for memory needs not
processing speed?

Hamish

_______________________________________________
grassuser mailing list
grassuser <at> grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

Charles Ehlschlaeger | 1 Mar 2007 01:52
Favicon

RE: [GRASS-user] 64 bit or parallel processing

Most GRASS raster functions, r.mapcalc and similar, won't speed up by more
than a couple of percent from parallelization. 90% of the time spent by
these commands comes from moving the data off the hard drive and getting it
to the CPU and visa versa (based on parallel processing research I did in
the mid '90s). For "normal" GRASS raster commands, having four hard drives
on RAID 0 will dramatically increase GRASS performance more than anything
one thing you could do to a computer (or GRASS itself). The reason ORNL
parallelized my code was that the computational requirements was over c^2,
where c is the number of grid cells on a map, for some of the analyses they
wanted to do.

As Hamish wrote in recent email, good candidates for parallelization would
be viewshed analysis or other functions requiring many cells to be read to
get a value for each output grid cell. r.watershed won't benefit from
parallelization.

Chuck Ehlschlaeger, Associate Professor & GIS Center Director
Department  of  Geography,      Western  Illinois  University
215 Tillman Hall,   1 University Circle,  Macomb,  IL   61455
cre111 <at> wiu.edu,    phone: 309-298-1841,     fax: 309-298-3003

-----Original Message-----
From: grassuser-bounces <at> grass.itc.it [mailto:grassuser-bounces <at> grass.itc.it]
On Behalf Of Glynn Clements
Sent: Wednesday, February 28, 2007 3:33 PM
To: Hamish
Cc: grassuser <at> grass.itc.it
Subject: Re: [GRASS-user] 64 bit or parallel processing

Hamish wrote:
(Continue reading)

Jonathan Greenberg | 1 Mar 2007 02:06
Picon

RE: [GRASS-user] 64 bit or parallel processing

I'm now working at UCD on a project with > 65 (!) hyperspectral images, and
a LOT of the slow downs stem from I/O intensive processes -- with that said,
there are plenty of "standard" raster processing routines that start
becoming processor intensive (e.g. linear spectral unmixing) that, with a
decent i/o setup (we have a lot of cheap-o 2-drive RAID0s that we back up to
a permanent location, this nearly doubles our I/O speed).  Not having any
idea of how grass's raster access backbone functions, it does seem you could
have a generic tiling routine that could be called from all of the raster
commands that would be highly parallizable (I note that RSI ENVI does this
out of the box now, it detects the # of processors and simply starts forking
a single image across all of the processors, I assume one line/tile per
processor).  Mapcalc and related processes would be a good start, and keep
in mind *spatial* analyses (e.g. texture windows) as you design it.  

My suggestion to the GRASS programmers: if you decide this is useful, I'd
HIGHLY recommend trying to do this in an out-of-the-box approach, e.g.
having an install-time "mp" capability -- if it requires being a) experts in
raster processing, b) experts in GRASS and c) experts in parallel coding,
you'll end up with 3 people doing SMP :)  

--j

-----Original Message-----
From: grassuser-bounces <at> grass.itc.it [mailto:grassuser-bounces <at> grass.itc.it]
On Behalf Of Hamish
Sent: Wednesday, February 28, 2007 4:14 PM
To: Glynn Clements
Cc: grassuser <at> grass.itc.it
Subject: Re: [GRASS-user] 64 bit or parallel processing

(Continue reading)

Dana Nibby | 1 Mar 2007 03:23
Picon
Favicon

[GRASS-user] Open Source GIS Gathering San Diego (OSG-SD) - June 18th, 2007

Folks, save the datum.

We're tentatively planning a 2nd Annual Open Source GIS Gathering to be held on Monday, June 18th - probably starting early afternoon, ending late in the evening. The same Monday of the Plenary Session(s) of another famous GIS conference held in San Diego.

Don't be lost in space... San Diego State University's Visualization Lab is easy to find via trolley from the San Diego Convention Center area.

We filled the Lab in 2006, and arrangements are being made to accommodate more. Suggested ideas for this year include technical workshops and roundtable discussions.
 
If you, or anyone you know, would be interested in presenting on an open source GIS topic, or conducting a technical workshop, please contact Dana Nibby here: danaspatial <at> yahoo.com.

There's now a mailing list for OSG-SD 2007. If you think you may be interested in atten ding, please sign-up below:

http://lists.telascience.org/mailman/listinfo/osg-sd

Spread the word to the extent possible. Geotag, you're "it."


Spatial Regards,

    Dana

Check out the all-new Yahoo! Mail beta - Fire up a more powerful email and get things done faster.
_______________________________________________
grassuser mailing list
grassuser <at> grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser
Glynn Clements | 1 Mar 2007 03:27

Re: [GRASS-user] 64 bit or parallel processing


Hamish wrote:

> Any thoughts on gains that could be made by MPI'ing the segmentation
> lib? Do do modules using that usually do so for memory needs not
> processing speed?

Modules which use the segment library normally do so because they
access the raster data in an order other than top-to-bottom.

If an algorithm reads the data in a linear fashion, there's no point
in using the segment library.

Ultimately, we need to look into replacing the core raster I/O code. 
However, this needs to be done incrementally; we can't afford to put
everything on hold while we re-write everything which uses
G_{get,put}_raster_row() etc (i.e. most of GRASS).

--

-- 
Glynn Clements <glynn <at> gclements.plus.com>

_______________________________________________
grassuser mailing list
grassuser <at> grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

Hamish | 1 Mar 2007 05:56
X-Face
Picon
Favicon

Re: [GRASS-user] help with ogr2ogr and .shp file

antonio rodriguez wrote:
> I need some help and advice to transform a .shp file to a format 
> readable by GRASS.

from within your lat lon GRASS location:

v.in.ogr location=

to import the shape file and create a new location for it based on its
.prj projection info. (assuming it has a .prj file. if not you'll have
to create the needed location manually and import using v.in.ogr from
within the new location)

> Extent: (100528.742188, 3991489.250000) - (610083.250000, 4288506.900918)
> Layer SRS WKT: (unknown)

Looks like you'll have to create a new location and set the projection
settings manually...

then (in the lat/lon location) suck it into the lat/lon projection with
v.proj.

Hamish

_______________________________________________
grassuser mailing list
grassuser <at> grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

Hamish | 1 Mar 2007 06:07
X-Face
Picon
Favicon

Re: [GRASS-user] r.reclass clarification - interactive mode

Roy Sanderson wrote:
> 
> I'm having problems using r.reclass interactively, in that it doesn't
> seem to be possible to display the full-screen text menu giving the
> category labels and category numbers of the old raster values, and the
> fields for the new values alongside.  Instead it is prompting me for
> the old and new values to be entered line by line.  Could you clarify
> whether this is a change from Grass 5.0 to 6.0, or have I not
> installed 6.0 correctly (using 6.2.2cvs, but had the same problem with
> 6.0).

The full screen interactive mode of r.reclass was removed for GRASS 6.
Try running r.report on the raster as a guide to the old cat values.

Hopefully one day soon we will have a fancy new wxPython GUI editor for
this.

Hamish

_______________________________________________
grassuser mailing list
grassuser <at> grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

Hamish | 1 Mar 2007 08:13
X-Face
Picon
Favicon

Re: [GRASS-user] cannot display vector file

Dave Kent wrote:
> 
> I recently installed Mr. Kyngesburye's CVS 6.3 070220 binaries and  
> new frameworks (updated from previous version).
> 
> When I try to display a vector in GIS manager  I get:
> 
> Size of 'coor' file differs from value saved in topo file.
> Please rebuild topology for vector 'EGEWells <at> Hartney'
> Cannot open topo file for vector 'EGEWells <at> Hartney'.
> 
> The vector is simple 2D point vector with a few attributes.
> I get the same error message from  d.vect in the terminal.
> 
> I tried to rebuild topology with v.build but the repaired file would  
> not display and would then no longer display in previous CVS version.
> Is there something I need to do with data when moving to new CVS  
> version?

You only need to rebuild if you had a very old version of GRASS you
might have to rebuild with v.build. Did the data files get damaged?
If migrating from GRASS 5 you might need to use v.convert.

does v.info show any features in the maps?

Hamish

_______________________________________________
grassuser mailing list
grassuser <at> grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser

Roger Bivand | 1 Mar 2007 10:56
Picon

RE: [GRASS-user] 64 bit or parallel processing

By the way, there is a forthcoming paper in Computers and Geosciences at:

http://www.sciencedirect.com/science?_ob=ArticleURL&_udi=B6V7D-4MY0TX9-8&_user=10&_coverDate=01%2F30%2F2007&_alid=543285067&_rdoc=1&_fmt=summary&_orig=search&_cdi=5840&_sort=d&_docanchor=&view=c&_ct=6&_acct=C000050221&_version=1&_urlVersion=0&_userid=10&md5=898ebe84b51937827e5e0b86a348340c

(we only see the abstract without paying):

Implementation of a parallel high-performance visualization technique in 
GRASS GIS by Alexandre Sorokine - the scripts are at:

http://www.ornl.gov/sci/gist/software/grass/

The approach involves tiling a region for visualization on a big wall from 
a head node to rendering nodes.

Roger

On Wed, 28 Feb 2007, Jonathan Greenberg wrote:

> I'm now working at UCD on a project with > 65 (!) hyperspectral images, and
> a LOT of the slow downs stem from I/O intensive processes -- with that said,
> there are plenty of "standard" raster processing routines that start
> becoming processor intensive (e.g. linear spectral unmixing) that, with a
> decent i/o setup (we have a lot of cheap-o 2-drive RAID0s that we back up to
> a permanent location, this nearly doubles our I/O speed).  Not having any
> idea of how grass's raster access backbone functions, it does seem you could
> have a generic tiling routine that could be called from all of the raster
> commands that would be highly parallizable (I note that RSI ENVI does this
> out of the box now, it detects the # of processors and simply starts forking
> a single image across all of the processors, I assume one line/tile per
> processor).  Mapcalc and related processes would be a good start, and keep
> in mind *spatial* analyses (e.g. texture windows) as you design it.  
> 
> My suggestion to the GRASS programmers: if you decide this is useful, I'd
> HIGHLY recommend trying to do this in an out-of-the-box approach, e.g.
> having an install-time "mp" capability -- if it requires being a) experts in
> raster processing, b) experts in GRASS and c) experts in parallel coding,
> you'll end up with 3 people doing SMP :)  
> 
> --j
> 
> -----Original Message-----
> From: grassuser-bounces <at> grass.itc.it [mailto:grassuser-bounces <at> grass.itc.it]
> On Behalf Of Hamish
> Sent: Wednesday, February 28, 2007 4:14 PM
> To: Glynn Clements
> Cc: grassuser <at> grass.itc.it
> Subject: Re: [GRASS-user] 64 bit or parallel processing
> 
> Glynn Clements wrote:
> > It would be more feasible to modify individual modules; e.g. a
> > multi-threaded version of r.mapcalc might be feasible. But unless the
> > raster I/O code was made thread-safe, you would still be limited by
> > the rate at which map data could be read and written by a single
> > thread, so it would only be worthwhile for cases where the bulk of the
> > overhead was in the actual calculations.
> 
> Right, so working on individual modules which can typically take 10+
> hours to run (*.rst, watersheds, viewsheds (incl. r.sun), indices from
> RS data, etc) seems like a much better focus of effort for the greatest
> gain:work ratio. And as seen in the MPI add-ons this is exactly where
> the previous ad-hoc effort has gone.
> 
> IIRC Helena said that the RST modules have changed substantially sice
> GRASS 5 and the s.surf.rst.mpi work would have to be largely redone.
> 
> Any thoughts on gains that could be made by MPI'ing the segmentation
> lib? Do do modules using that usually do so for memory needs not
> processing speed?
> 
> 
> Hamish
> 
> _______________________________________________
> grassuser mailing list
> grassuser <at> grass.itc.it
> http://grass.itc.it/mailman/listinfo/grassuser
> 
> 
> _______________________________________________
> grassuser mailing list
> grassuser <at> grass.itc.it
> http://grass.itc.it/mailman/listinfo/grassuser
> 

--

-- 
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: Roger.Bivand <at> nhh.no

_______________________________________________
grassuser mailing list
grassuser <at> grass.itc.it
http://grass.itc.it/mailman/listinfo/grassuser


Gmane