[GRASS-user] I do not understand how region/projection/coordinate system are intertwined

2012-03-14 Thread Fridtjof Schiefenhövel
Hello there! it seems I can't solve this problem: ->create a location with a mercator projection (i want to have degrees as the coordinate system) ->import landsat geotiffs ->have degrees as units showing (e.g. when moving the mouse) this is what I did: -created a location using the projecti

Re: [GRASS-user] using v.to.rast

2012-03-14 Thread Johannes Radinger
Hi Janet, if you choose "use=attr" then the value of the attribute you set with "column=xxx" is used as the value for the raster pixel. You can check your attribute table to find out which attribute values are assign to each road section. If you just want a "fix value" for each raster pixel of t

[GRASS-user] postgreSQL import

2012-03-14 Thread Patrick S.
I try to import various files from PostgreSQL to GRASS in GRASS70. / db.connect driver=pg database='host=,dbname=,schema=' db.login user="sustaincity" pass="+++" db.connect -p db.tables -p |less for i in {1..3}; do v.in.ogr dsn='PG:host= dbname= user=' layer=schema.layer output=lay

[GRASS-user] i.vi.mpi

2012-03-14 Thread Ripsime Khachatryan
Hi all This is my Makefile MODULE_TOPDIR = /usr/local/grass-6.4.2svn/ CC=mpicc --static PGM = i.vi.mpi LIBES = $(GISLIB) $(GMATHLIB) DEPENDENCIES = $(GISDEP) $(GMATHDEP) include $(MODULE_TOPDIR)/include/Make/Module.make ifneq ($(USE_LARGEFILES),) EXTRA_CFLAGS = -D_FILE_OFFSET_BITS=64 endif def

Re: [GRASS-user] I do not understand how region/projection/coordinate system are intertwined

2012-03-14 Thread Stephen Sefick
Is the geotiff in degrees? If you are on unix use gdalinfo to find out the projection of the geotiff. You will have to import in its native projection and then reporject into whatever coordinate system that you would like. I am just drinking my coffee this morning, so I may not be awake yet.

[GRASS-user] Re: postgreSQL import (Patrick S.)

2012-03-14 Thread Patrick S.
Just found the problem with PostgreSQL-import: If the postGIS-layer has multiple geometries (in my case, centroid and polygon), it will be viewed with "db.tables -p", but not imported through "v.in.ogr". Error is misleading as it says "Layer not available". Workaround, was a copy of the layer

Re: [GRASS-user] I do not understand how region/projection/coordinate system are intertwined

2012-03-14 Thread Stephen Sefick
I am not an expert with coordinate projections. I think that UTM coordinates are in meters from a false easting and northing. This is a very good projection for a "small" part of an ellipsoid. All other comments are inline. I have reposted this to the list because maybe somebody with a bett

[GRASS-user] grass prompt when changing locations

2012-03-14 Thread Tyler Smith
Hi, When I change locations within a grass session with g.mapset, the prompt is not updated. If I start grass with: grass -text ~/grassdata/dem/PERMANENT my prompt looks like this: GRASS 6.4.1 (dem):~ > The reference to the location '(dem)' stays, even after changing to another location. How d

Re: [GRASS-user] grass prompt when changing locations

2012-03-14 Thread Micha Silver
On 03/14/2012 03:55 PM, Tyler Smith wrote: Hi, When I change locations within a grass session with g.mapset, the prompt is not updated. If I start grass with: grass -text ~/grassdata/dem/PERMANENT my prompt looks like this: GRASS 6.4.1 (dem):~ > The reference to the locat

Re: [GRASS-user] using v.to.rast

2012-03-14 Thread Janet Choate
Johannes, thank you very much, that was all i needed. so simple, i appreciate the explanation! thanx, Janet On Wed, Mar 14, 2012 at 1:12 AM, Johannes Radinger wrote: > Hi Janet, > > if you choose "use=attr" then the value of the > attribute you set with "column=xxx" is used > as the value for t

[GRASS-user] Why Does Re-projection Increase Disk Space Used?

2012-03-14 Thread Rich Shepard
I have a source 28M DEM in a directory with geographic coordinates in Lat/Long. When I re-project that to another location and dataset with Lambert Conformal Conic projection the disk space consumed bloats to 106G. That's more than a 20x increase. Why? In the destination location/directory th

Re: [GRASS-user] Why Does Re-projection Increase Disk Space Used? [UPDATE]

2012-03-14 Thread Rich Shepard
On Wed, 14 Mar 2012, Rich Shepard wrote: In the destination location/directory there's no .tmp/ subdirectory with a huge file. Actually, there are two such files. Am I correct in assuming that once the re-projection process is completed I can delete those two files (called -1 and -2) wit

Re: [GRASS-user] grass prompt when changing locations

2012-03-14 Thread Tyler Smith
On Wed, Mar 14, 2012 at 10:33 AM, Micha Silver wrote: > On 03/14/2012 03:55 PM, Tyler Smith wrote: > > > my prompt looks like this: > > GRASS 6.4.1 (dem):~ > > > The reference to the location '(dem)' stays, even after changing to > another location. How do I make Grass update this? > > > The new e

Re: [GRASS-user] v.overlay missing polygon

2012-03-14 Thread Falko Engel
Am 13.03.2012 07:53, schrieb Vincent Bain: Falko, if you try to make things match, looks like v.overlay returns a correct result. Try the following command -- threshold 5 m is arbitrary, of course : v.edit map=plot type=boundary tool=move move=0,0 thresh=-1,5,0 bgmap=veg snap=vertex bbox=3607160

Re: [GRASS-user] Why Does Re-projection Increase Disk Space Used? [UPDATE]

2012-03-14 Thread Mark Seibel
Hi. On Wed, Mar 14, 2012 at 11:33 AM, Rich Shepard wrote: > On Wed, 14 Mar 2012, Rich Shepard wrote: > > Related question: if they're no longer needed, why doesn't grass remove > them from the hard drive? > If you exit GRASS does it purge them? Mark

Re: [GRASS-user] Why Does Re-projection Increase Disk Space Used? [UPDATE]

2012-03-14 Thread Rich Shepard
On Wed, 14 Mar 2012, Mark Seibel wrote: If you exit GRASS does it purge them? Mark, Not on my laptop on which I ran the r.proj module. At the bash prompt I used rm to remove them. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http:

[GRASS-user] nviz: Multiple Vector Maps Draped on Raster DEM

2012-03-14 Thread Rich Shepard
From a thread a couple of months ago I learned that nviz will allow only a single vector map to be draped over the raster base map. For presentation purposes using ps.map I need multiple vector maps to be displayed on the raster map. Is the best approach to use v.overlay (with the union/OR op

Re: [GRASS-user] nviz: Multiple Vector Maps Draped on Raster DEM

2012-03-14 Thread Juan David Rondón Dí­az
On 14/03/12 12:47, Rich Shepard wrote:   From a thread a couple of months ago I learned that nviz will allow only a single vector map to be draped over the raster base map. For presentation purposes using ps.map I need multiple vector maps to be display

[GRASS-user] v.patch

2012-03-14 Thread Patrick S.
Dear list, I am trying to merge several douzen vectorlayers in a loop. These are representing zones of neighboring municipalities, that have been imported from PostgreSQL and might overlap. Overlapping is not expected to be of relevant size, so it is more important to keep the number of objec

[GRASS-user] Re: r.region for intersection with vector map?

2012-03-14 Thread Hermann Peifer
On 14/03/2012 19:02, Rich Shepard wrote: I assume the reason is that the region is a square. Indeed. What command could I pass to r.mapcalc that would clip the raster map to the basins themselves Define a MASK before running r.mapcalc http://grass.osgeo.org/grass64/manuals/html64_user/r.ma

[GRASS-user] import text file with v.in.ascii and include labels

2012-03-14 Thread Janet Choate
Hi users, i am attempting to import a text file of coordinates and labels with v.in.ascii. i have tried the examples on the v.in.ascii help page, and looked through other posts, but still not getting it to do what i want. my text file looks like this: 1|559347.085|4895479.776597|LookoutCreek

Re: [GRASS-user] Re: r.region for intersection with vector map?

2012-03-14 Thread Rich Shepard
On Wed, 14 Mar 2012, Hermann Peifer wrote: Define a MASK before running r.mapcalc http://grass.osgeo.org/grass64/manuals/html64_user/r.mask.html Hermann, Duh! Of course! Thanks, Rich ___ grass-user mailing list grass-user@lists.osgeo.org http:/

[GRASS-user] Display Vector Areas Without Centroid Symbol

2012-03-14 Thread Rich Shepard
I don't see an option in the GUI properties dialog box to hide or remove the little 'x' marking area centroids. Is there a way to not display those when the vector map overlays a raster map? Rich ___ grass-user mailing list grass-user@lists.osgeo.org

Re: [GRASS-user] Display Vector Areas Without Centroid Symbol

2012-03-14 Thread Rich Shepard
On Wed, 14 Mar 2012, Martin Landa wrote: Selection -> Feature types uncheck 'centroid' Martin, Thank you. I keep overlooking that line Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Display Vector Areas Without Centroid Symbol

2012-03-14 Thread Martin Landa
2012/3/14 Rich Shepard : >  I don't see an option in the GUI properties dialog box to hide or remove > the little 'x' marking area centroids. Is there a way to not display those > when the vector map overlays a raster map? Selection -> Feature types uncheck 'centroid' Martin -- Martin Landa *

Re: [GRASS-user] Why Does Re-projection Increase Disk Space Used?

2012-03-14 Thread Daniel Victoria
Some things come to mind: 1) The resolution has changed 2) the number of bits per pixel has changed 3) the region of the reprojected file is bigger than the original file so you have lots of nulls. Daniel On Wed, Mar 14, 2012 at 12:28 PM, Rich Shepard wrote: >  I have a source 28M DEM in a direc

Re: [GRASS-user] import text file with v.in.ascii and include labels

2012-03-14 Thread Falko Engel
On 14.03.2012 20:50, Janet Choate wrote: cat gages.txt | v.in.ascii out=gages x=2 y=3 cat=1 \ columns='x double precision, y double precision, cat int, name varchar(15)' Hi Janet, I think your command is just mixed up a little. Try: cat gages.txt | v.in.ascii out=gages x=3 y=2 cat=1 columns='

[GRASS-user] Re: i.vi.mpi

2012-03-14 Thread Alexander Muriy
Hello Ripsime. This is my Makefile > > > MODULE_TOPDIR = /usr/local/grass-6.4.2svn/ > CC=mpicc --static > PGM = i.vi.mpi > LIBES = $(GISLIB) $(GMATHLIB) > DEPENDENCIES = $(GISDEP) $(GMATHDEP) > include $(MODULE_TOPDIR)/include/Make/Module.make > ifneq ($(USE_LARGEFILES),) > EXTRA_CFLAGS = -D_FILE

Re: [GRASS-user] v.overlay missing polygon

2012-03-14 Thread Vincent Bain
Never tried this. I am generally very cautious with snapping. Not knowing exactly how grass handles it but one should be aware of possible topological alterations, duplicate features, etc. Vincent > > Dear Vincent, > > I was able to reproduce your code. This is definitely an improvement! > For

RE: [GRASS-user] Re: i.vi.mpi

2012-03-14 Thread Chemin, Yann (IWMI)
Hi all, Added a GRASS7 version of i.vi.mpi in addons/grass7/imagery/ Yann From: grass-user-boun...@lists.osgeo.org [grass-user-boun...@lists.osgeo.org] on behalf of Alexander Muriy [amu...@gmail.com] Sent: Thursday, March 15, 2012 12:12 PM To: grass-user@lists.os