Re: [GRASS-user] Loosing color table when reprojecting raster file

2012-02-13 Thread Markus Neteler
On Mon, Feb 13, 2012 at 3:42 AM, Kurt Springs ferret_b...@mac.com wrote: Hi folks, I was preparing a some maps for a presentation and ran into a snag. I am trying to reproject raster maps from a Lat/Long to Irish national grid.  I've done this before.  This time I found that I had no color.

Re: [GRASS-user] How to create a project

2012-02-13 Thread Moritz Lennert
On 11/02/12 17:32, albertsaribek...@rambler.ru wrote: thank You for response, I'm interested in is it possible to create some project for GRASS without using the GRASS? GRASS uses its own internal format (and directory structure) for data. You, therefore, cannot create a project such as

Re: [GRASS-user] How to link or combine adjacent vector lines

2012-02-13 Thread Moritz Lennert
On 12/02/12 22:20, Stefan Luedtke wrote: That did the trick, almost. They join each other -so v.build polylines ... worked, but, the table is exactly the same as it was before. The command output said, that the number of lines was reduced to almost 50%. Furthermore, if I check the cat number

Re: [GRASS-user] Re: How to rescale an aspect map

2012-02-13 Thread Moritz Lennert
On 13/02/12 08:31, leonidas wrote: I need the most suitable areas (with SE aspect) to have the value of 255 and the unsuitable (North aspect) to have the value of 0. Grass returns N=90 and SE=315 and flat=0. So the minimum value of aspect raster is flat areas and the max is areas with 359

[GRASS-user] d.rast.arrows in GRASS6.4.2RC3 wxpython GUI

2012-02-13 Thread Aldo Clerici
GRASS users and developers. In GRASS6.4.2RC3 with wxpython GUI the option 'Add raster cell arrows' seems not to work properly. It seems that the arrows in the Map Display cannot zoomed in.. It can be checked by the following steps: Load and display the aspect (Spearfish) map. Open the

[GRASS-user] Re: How to rescale an aspect map

2012-02-13 Thread leonidas
Moritz, how did you choose 170? -- View this message in context: http://osgeo-org.1560.n6.nabble.com/How-to-rescale-an-aspect-map-tp4377272p4464502.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list

Re: [GRASS-user] v.to.rast in a python loop, error

2012-02-13 Thread Johannes Radinger
Dear Glynn, Dear Moritz, Dear Scott, Have you tried adding flags = '-o' as an argument to grass.run_command? That will overwrite the output raster if it already exists. Actually, it will raise a ScriptError if you're running anything later than 6.4.1. If you want to pass the

Re: [GRASS-user] Re: How to rescale an aspect map

2012-02-13 Thread Moritz Lennert
On 13/02/12 11:16, leonidas wrote: Moritz, how did you choose 170? 315 = 255 90 = 0 Between the two you have 135 degrees (45 from 315 to 360 plus 90 from 1 to 90). 45 / 135 = 1/3, so 360 is at 255 - 255/3 = 170. Obviously this means that you go from 255 to 0 much quicker when you move

Re: [GRASS-user] v.to.rast in a python loop, error

2012-02-13 Thread Moritz Lennert
On 13/02/12 11:32, Johannes Radinger wrote: And my output with although all maps were removed manually with g.remove and cross checked: Launching script '/Users/Johannes Radinger/Desktop/test2.py'... (Mon Feb 13 11:12:27 2012) /Users/Johannes Radinger/Desktop/test2.py

Re: [GRASS-user] d.rast.arrows in GRASS6.4.2RC3 wxpython GUI

2012-02-13 Thread Hamish
Aldo wrote: Something wrong in my procedure or is this a bug in wxpython GUI (with tcltk GUI the procedure works fine)? I think it's a known bug in the wxGUI. It's very similar to the d.rast.num module. you might also try the d.barb module from the addons repository, maybe more luck with

Re: [GRASS-user] How to link or combine adjacent vector lines

2012-02-13 Thread Hamish
Moritz wrote: With options first and multi v.build.polylines just copies the attribute tables as is, i.e. it does not delete any lines for which there are no more category values in the vector file. To create an attribute table that corresponds to the current vector amp with the following

Re: [GRASS-user] How to link or combine adjacent vector lines

2012-02-13 Thread Moritz Lennert
On 13/02/12 12:06, Hamish wrote: Moritz wrote: With options first and multi v.build.polylines just copies the attribute tables as is, i.e. it does not delete any lines for which there are no more category values in the vector file. To create an attribute table that corresponds to the current

[GRASS-user] Re: v.to.rast in a python loop, error

2012-02-13 Thread Margherita Di Leo
Hi, Have you tried adding flags = '-o' as an argument to grass.run_command? That will overwrite the output raster if it already exists. Actually, it will raise a ScriptError if you're running anything later than 6.4.1. If you want to pass the --overwrite flag to the module,

Re: [GRASS-user] v.to.rast in a python loop, error

2012-02-13 Thread Johannes Radinger
Original-Nachricht Datum: Mon, 13 Feb 2012 11:46:21 +0100 Von: Moritz Lennert mlenn...@club.worldonline.be An: Johannes Radinger jradin...@gmx.at CC: grass-user@lists.osgeo.org, gl...@gclements.plus.com Betreff: Re: [GRASS-user] v.to.rast in a python loop, error On

[GRASS-user] Re: v.to.rast in a python loop, error

2012-02-13 Thread Johannes Radinger
Hi, Original-Nachricht Datum: Mon, 13 Feb 2012 12:20:18 +0100 Von: Margherita Di Leo direg...@gmail.com An: grass-user@lists.osgeo.org, jradin...@gmx.at, gl...@gclements.plus.com, scott.sinclair...@gmail.com Betreff: Re: v.to.rast in a python loop, error Hi,

[GRASS-user] Not being able to produce decent-look color-composites

2012-02-13 Thread Luisa Peña
Greetings I'have 3 layers of Landsat (Green, Red and Near-Infrared) and I'm trying to produce a false-color composite (r.composite( with those but... I'm getting very frustrated because I'm producing very bad look composites. My layers are reflectances from 0-1.8 or something. Shall I rescale

Re: [GRASS-user] v.to.rast in a python loop, error

2012-02-13 Thread Moritz Lennert
On 13/02/12 12:26, Johannes Radinger wrote: Original-Nachricht Datum: Mon, 13 Feb 2012 11:46:21 +0100 Von: Moritz Lennertmlenn...@club.worldonline.be An: Johannes Radingerjradin...@gmx.at CC: grass-user@lists.osgeo.org, gl...@gclements.plus.com Betreff: Re: [GRASS-user]

Re: [GRASS-user] v.to.rast in a python loop, error

2012-02-13 Thread Johannes Radinger
Original-Nachricht Datum: Mon, 13 Feb 2012 13:18:58 +0100 Von: Moritz Lennert mlenn...@club.worldonline.be An: Johannes Radinger jradin...@gmx.at CC: grass-user@lists.osgeo.org, gl...@gclements.plus.com Betreff: Re: [GRASS-user] v.to.rast in a python loop, error On

Re: [GRASS-user] two questions: snap and close lines

2012-02-13 Thread Markus Metz
On Sun, Feb 12, 2012 at 8:58 PM, Paolo Cavallini cavall...@faunalia.it wrote: Sample location here: http://int.faunalia.it/~paolo/test.tar.gz Confirmed by other users. Thanks. Hmm. Breaking does indeed not work, even combining it with a bunch of other cleaning tools. No quick solution here...

[GRASS-user] Possible databases for GRASS

2012-02-13 Thread Thorsten
Hello List, how flexible is GRASS with regards to the DB where the GIS-data are stored? I read in the docs that GRASS can now work with many popular relational DB systems. But what about an entirely unconventional data store like an object-oriented database implemented in lisp? Is it just a

Re: [GRASS-user] Loosing color table when reprojecting raster file

2012-02-13 Thread Kurt Springs
The map is a digital elevation model. Interpolation is set to nearest. On Feb 13, 2012, at 3:46 AM, Markus Neteler wrote: On Mon, Feb 13, 2012 at 3:42 AM, Kurt Springs ferret_b...@mac.com wrote: Hi folks, I was preparing a some maps for a presentation and ran into a snag. I am trying to

Re: [GRASS-user] Possible databases for GRASS

2012-02-13 Thread Moritz Lennert
On 13/02/12 15:49, Thorsten wrote: Hello List, how flexible is GRASS with regards to the DB where the GIS-data are stored? I read in the docs that GRASS can now work with many popular relational DB systems. This relates to storing attribute values, not the spatial feature information. GRASS

Re: [GRASS-user] v.to.rast in a python loop, error

2012-02-13 Thread Moritz Lennert
On 13/02/12 15:19, Johannes Radinger wrote: Original-Nachricht Datum: Mon, 13 Feb 2012 13:18:58 +0100 Von: Moritz Lennertmlenn...@club.worldonline.be An: Johannes Radingerjradin...@gmx.at CC: grass-user@lists.osgeo.org, gl...@gclements.plus.com Betreff: Re: [GRASS-user]

Re: [GRASS-user] Using r.buffer on image with geographic coordinates

2012-02-13 Thread Markus Neteler
On Fri, Mar 25, 2011 at 2:42 PM, Ned Horning horn...@amnh.org wrote: Hi - I am trying to create buffers for large (~1GB) line and polygon vector files. I tried v.buffer but after letting it run overnight with no apparent progress I killed the process. I'm not aware of any software capable of

Re: [GRASS-user] v.to.rast in a python loop, error

2012-02-13 Thread Moritz Lennert
On 13/02/12 16:11, Moritz Lennert wrote: On 13/02/12 15:19, Johannes Radinger wrote: #!/usr/bin/env python import sys import time import grass.script as grass import grass.script.setup as gsetup def main(): print start sys.stdout.flush() grass.run_command(g.list, type = rast)

Re: [GRASS-user] Loosing color table when reprojecting raster file

2012-02-13 Thread Markus Neteler
On Mon, Feb 13, 2012 at 3:55 PM, Kurt Springs ferret_b...@mac.com wrote: The map is a digital elevation model.  Interpolation is set to nearest. For a DEM I would avoid nearest to obtain a better (smooth) result. Markus ___ grass-user mailing list

[GRASS-user] d.zoom.keys -- shell script for navigation in GRASS monitors from keyboard

2012-02-13 Thread Alexander Muriy
Dear GRASS'ers :) Let me introduce (a bit hacky) shell script ( https://raw.github.com/amuriy/GRASS-scripts/master/d.zoom.keys) for navigating in classic X-monitors with the typical keys (arrows and - / +). It works in two modes: from terminal or directly from X-monitor (requires xev and xdotool

Re: [GRASS-user] Not being able to produce decent-look color-composites

2012-02-13 Thread Markus Neteler
On Mon, Feb 13, 2012 at 12:50 PM, Luisa Peña luisapena1...@gmail.com wrote: Greetings I'have 3 layers of Landsat (Green, Red and Near-Infrared) and I'm trying to produce a false-color composite (r.composite( with those but... I'm getting very frustrated because I'm producing very bad look

Re: [GRASS-user] Not being able to produce decent-look color-composites

2012-02-13 Thread Michael Barton
Also, it usually helps to do a histogram equalization on each band prior to compositing. i.landsat.rgb does some of this kind of enhancement automatically. Michael On Feb 13, 2012, at 10:00 AM, grass-user-requ...@lists.osgeo.org wrote: Date: Mon, 13 Feb 2012 17:22:50 +0100 From: Markus

Re: [GRASS-user] How to link or combine adjacent vector lines

2012-02-13 Thread Stefan Luedtke
I tried the short one of course , worked perfectly ;-) Thanks a lot for your efforts. Cheers, stefan On Mon, 2012-02-13 at 03:06 -0800, Hamish wrote: Moritz wrote: With options first and multi v.build.polylines just copies the attribute tables as is, i.e. it does not delete any

[GRASS-user] Re: Possible databases for GRASS

2012-02-13 Thread Thorsten
Moritz Lennert mlenn...@club.worldonline.be writes: On 13/02/12 15:49, Thorsten wrote: Hello Moritz, how flexible is GRASS with regards to the DB where the GIS-data are stored? I read in the docs that GRASS can now work with many popular relational DB systems. This relates to storing

Re: [GRASS-user] Not being able to produce decent-look color-composites

2012-02-13 Thread Hamish
Luisa wrote: I'have 3 layers of Landsat (Green, Red and Near-Infrared) and I'm trying to produce a false-color composite (r.composite( with those but... I'm getting very frustrated because I'm producing very bad look composites. My layers are reflectances from 0-1.8 or something. Shall I

[GRASS-user] Rotating/sizing labels based on DB column?

2012-02-13 Thread Asher Hoskins
Hello. I'm using some Ordnance Survey street-level Mastermap data which implements street names, house numbers, etc. as map points linked to a database containing the text to print, its size and its rotation. I can produce basic map labels using v.label: v.label labels=xxx map=g26a

Re: [GRASS-user] v.to.rast in a python loop, error

2012-02-13 Thread Hamish
anyone know if python has a portable version of 'sync' from linux, and/or 'wait' from bash? Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Loosing color table when reprojecting raster file

2012-02-13 Thread Kurt Springs
I tried shifting the interpolation method. Nothing worked. When I try to change the color rules, in 6.4.1 it runs the change, but doesn't add any color. When I try in 7.0 svn, it grays the apply and okay button. Kurt On Feb 13, 2012, at 10:58 AM, Markus Neteler wrote: On Mon, Feb 13, 2012

Re: [GRASS-user] Not being able to produce decent-look color-composites

2012-02-13 Thread Hamish
MarkusN wrote: Did you try with i.landsat.rgb? I didn't check if it works with reflectance, too. BTW: I just added a screenshot example to http://grass.osgeo.org/grass64/manuals/html64_user/i.landsat.rgb.html see also http://grass.osgeo.org/screenshots/rs.php and the commissioning trials

Re: [GRASS-user] Not being able to produce decent-look color-composites

2012-02-13 Thread Glynn Clements
Luisa Peña wrote: I'have 3 layers of Landsat (Green, Red and Near-Infrared) and I'm trying to produce a false-color composite (r.composite( with those but... I'm getting very frustrated because I'm producing very bad look composites. My layers are reflectances from 0-1.8 or something. Shall

Re: [GRASS-user] v.to.rast in a python loop, error

2012-02-13 Thread Glynn Clements
Hamish wrote: anyone know if python has a portable version of 'sync' from linux, and/or 'wait' from bash? There is no reason to use 'sync' or an equivalent. sync causes modified cache blocks to be written to the underlying block device. As processes read from the cache (not the physical

Re: [GRASS-user] v.to.rast in a python loop, error

2012-02-13 Thread Glynn Clements
Johannes Radinger wrote: Looks weird as 2nd g.list over is printed after the ERROR, resp. the ERROR is produced by v.to.rast before the print statement is processed... This is almost certainly an artifact of the GUI. Python's print statement writes to stdout by default, while GRASS error

[GRASS-user] Handling 700 Dtm files in Grass

2012-02-13 Thread Paul Shapley
Hi, I need to merge 700 arcinfo ascii files so they can be imported into Grass then saved as a single geotif, How can this best be done?. It is possible to import them one at a time but this is not practical. -- Paul J. Shapley ___ grass-user mailing

Re: [GRASS-user] Handling 700 Dtm files in Grass

2012-02-13 Thread Hamish
Paul wrote: I need to merge 700 arcinfo ascii files so they can be imported into Grass then saved as a single geotif, How can this best be done?. It is possible to import them one at a time but this is not practical. a simple unix shell loop will do it, (use with MSys on MS-Windows) for file