Re: [GRASS-user] 6.4 Release Candidate versions
Roy wrote: Quick query - are there any plans to provide Grass 6.4 binaries in the current Ubuntu LTS format? Unfortunately this is still Ubuntu 8.04 with the next LTS not due until the middle of 2010. you'd have to check with the UbuntuGIS team. http://lists.osgeo.org/mailman/listinfo/ubuntu Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] 6.4 Release Candidate versions
Hamish wrote: Roy wrote: Quick query - are there any plans to provide Grass 6.4 binaries in the current Ubuntu LTS format? Unfortunately this is still Ubuntu 8.04 with the next LTS not due until the middle of 2010. you'd have to check with the UbuntuGIS team. http://lists.osgeo.org/mailman/listinfo/ubuntu Hamish https://launchpad.net/~ubuntugis/+archive/ubuntugis-unstable 6.4 RC5 for Hardy is available and I use it. I don't know of any plans not to keep it up to date until the next LTS comes out. And yes the UbuntuGIS mailing list is the place to keep up on that stuff. Alex ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error while onitializing GRASS6.4 with wxpython
Hi, 2009/12/15 Nikos Dumakis nikosdu1...@gmail.com: I have installed and compiled GRASS6.4 in an Linux server. It's working OK in text mode and it opens the menu with tcltk. But when I run grass64 -wxpython I get this errors: Traceback (most recent call last): File /usr/local/grass-6.4.0RC5/etc/wxpython/gis_set.py, line 33, in ? from gui_modules import globalvar File /usr/local/grass-6.4.0RC5/etc/wxpython/gui_modules/globalvar.py, line 55, in ? CheckForWx() File /usr/local/grass-6.4.0RC5/etc/wxpython/gui_modules/globalvar.py, line 44, in CheckForWx except (ImportError, ValueError, wxversion.VersionError), e: UnboundLocalError: local variable 'wxversion' referenced before assignment Error in GUI startup. If necessary, please report this error to the GRASS developers. Switching to text mode now. Do anyone has a clue of what is not right with WXPYTHON GUI startup? try 6.4 from svn. The bug should be already fixed. Martin -- Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error while onitializing GRASS6.4 with wxpython
2009/12/16 Martin Landa landa.mar...@gmail.com: try 6.4 from svn. The bug should be already fixed. Anyway install wxversion on your system. Martin -- Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] g.region (was v.to.rast conversion)
Hi Markus, I'm not using the GUI because I need to use GRASS as automated tool to convert a shape file to a raster in GTIFF. I tried to do what you have explained in many ways but without success. I know that I have to execute something like: v.in.ogr ... v.to.rast ... r.out.gdal ... But the v.in.ogr returns the g.region error also with the location= option. I tried to do another operation (maybe simpler), like reading a raster image and get information on the region, diong: r.in.gdal -o -e input=$HOME/raster.tif output=grass_rast location=$LOCATION Then setting the new location (as you suggested) and running the g.region g.region -p and this is the output: Location test created 100% r.in.gdal complete. Raster map grass_rast.1 created. Copying 144 GCPS in points file for grass_rast GCPs have the following OpenGIS WKT Coordinate System:: GEOGCS[WGS 84,DATUM[WGS_1984,SPHEROID[WGS 84,6378137,298.2572235629972,AUTHORITY[EPSG,7030]],AUTHORITY[EPSG,6326]],PRIMEM[Greenwich,0],UNIT[degree,0.0174532925199433],AUTHORITY[EPSG,4326]] projection: 0 (x,y) zone: 0 north: 2030 south: 0 west: 0 east: 1354 nsres: 1 ewres: 1 rows: 2030 cols: 1354 cells: 2748620 It seems to work properly... Why I can't run the v.in.ogr? Thank you for help. Markus Metz wrote: Hi, According to the grass book, the v.in.ogr command should include the location option: v.in.ogr -o -e dsn=./test_data/test_data.shp output=grass_map location=$LOCATION then go to the newly created location (set the environment variables accordingly and update .grassrc6), set the region extends and resolution to your demands, run v.to.rast, then r.out.gdal. IMHO, it is easier to use the GUI. Hope that helps, Markus M -- Giacomo ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.out.ascii output modification
Hello everybody, As the title says, I'm trying to use r.profile with the output of r.drain. I know I cannot use directly the output, and instead I need to transform r.drain to a vector line and transform it to a ascii point file (v.out.ascii). So, when the output form r.drain has been transformed to a ascii point file, I cannot use it because there is a third column (the category) preventing me from piping it into r.profile. I tried deleting the category but this doesn't work. I can always use awk to delete characters in the ascii file that I don't need, but I was wondering if there is a quicker way to do it. Basically what I need is to get from this ascii point file: X Y CAT xxx|yyy|1 xxx|yyy|1 to: XY xxx|yyy without using a awk-like method. Is it possible? Thank you, Maxime -- View this message in context: http://n2.nabble.com/v-out-ascii-output-modification-tp4175638p4175638.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] g.region (was v.to.rast conversion)
Hi, Giacomo Piva wrote: Hi Markus, I'm not using the GUI because I need to use GRASS as automated tool to convert a shape file to a raster in GTIFF. The GUI is good to develop and debug the script before you create an automated tool that may produce bogus results. I tried to do what you have explained in many ways but without success. I know that I have to execute something like: v.in.ogr ... v.to.rast ... r.out.gdal ... But the v.in.ogr returns the g.region error also with the location= option. Then try a valid location to start with, as explained by Glynn. I think you can use any valid location, also one of the sample datasets. When using an existing valid location, it should be enough to give the location name and mapset name as in your script. Then do v.in.ogr -oe location=$LOCATION That obviously only works properly if there is a *.prj file accompanying the *.shp file, otherwise you end up with a xy location. In the absence of a *.prj file you must either ask the authors of the shapefile or guess and manually create the location with the correct projection info. Markus M I tried to do another operation (maybe simpler), like reading a raster image and get information on the region, diong: r.in.gdal -o -e input=$HOME/raster.tif output=grass_rast location=$LOCATION Then setting the new location (as you suggested) and running the g.region g.region -p and this is the output: Location test created 100% r.in.gdal complete. Raster map grass_rast.1 created. Copying 144 GCPS in points file for grass_rast GCPs have the following OpenGIS WKT Coordinate System:: GEOGCS[WGS 84,DATUM[WGS_1984,SPHEROID[WGS 84,6378137,298.2572235629972,AUTHORITY[EPSG,7030]],AUTHORITY[EPSG,6326]],PRIMEM[Greenwich,0],UNIT[degree,0.0174532925199433],AUTHORITY[EPSG,4326]] projection: 0 (x,y) zone: 0 north: 2030 south: 0 west: 0 east: 1354 nsres: 1 ewres: 1 rows: 2030 cols: 1354 cells: 2748620 It seems to work properly... Why I can't run the v.in.ogr? Thank you for help. Markus Metz wrote: Hi, According to the grass book, the v.in.ogr command should include the location option: v.in.ogr -o -e dsn=./test_data/test_data.shp output=grass_map location=$LOCATION then go to the newly created location (set the environment variables accordingly and update .grassrc6), set the region extends and resolution to your demands, run v.to.rast, then r.out.gdal. IMHO, it is easier to use the GUI. Hope that helps, Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] V.overlay very slow
By the way, i don't know if what i've done is correct and can be a usefull trick, or it's a rookie mistake, or a casuality. I hope someone knows why this happens, and could help anyone. This is: After trying the v.patch, v.clean(break,rmdupl), v.build thing, it was the moment to create a new column for all geometries, in order to run v.dissolve. Then, v.dissolve doesn't worked. I see that v.info gives different information (number of boundaries, areas, islands...) of the same map, obtained from v.overlay's or v.patch,v.clean... Then, I've tryed this: v.patch input=quad1,quad2,quad3.quadN output=quadS v.overlay ainput=quadS binput=quadS operation=or output=MAP And then, I've obtained the same map that running v.overlay quad1 or quad2 - sum v.overlay quad3 or sum - sum v.overlay quadN or sum - MAP but MUCH MORE FASTER, obviously. After that, I passed that MAP to create the common column and running v.dissolve, getting the right result. -- View this message in context: http://n2.nabble.com/V-overlay-very-slow-tp4169783p4175821.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] V.overlay very slow
incanus wrote: By the way, i don't know if what i've done is correct and can be a usefull trick, or it's a rookie mistake, or a casuality. I hope someone knows why this happens, and could help anyone. This is: After trying the v.patch, v.clean(break,rmdupl), v.build thing, v.build is not needed after v.clean. it was the moment to create a new column for all geometries, in order to run v.dissolve. Then, v.dissolve doesn't worked. I see that v.info gives different information (number of boundaries, areas, islands...) of the same map, obtained from v.overlay's or v.patch,v.clean... Can you post that information, or better the output of v.build for each map (not all the input maps, only the results of v.patch and v.overlay)? Look out for number of incorrect boundaries, number of duplicate centroids and number of centroids outside area, that means that there is still something wrong with vector topology. Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] V.overlay very slow
that's an example with quad3 + quad4 OUTPUT OF V.PATCH Patching vector map qu...@andresfr@andresfr... 42 primitives registered 922 vertices registered Building areas... 26 areas built 24 isles built Attaching islands... Attaching centroids... Number of nodes: 40 Number of primitives: 42 Number of points: 0 Number of lines: 0 Number of boundaries: 26 Number of centroids: 16 Number of areas: 26 Number of isles: 24 Number of areas without centroid: 10 Intersections at borders will have to be snapped Lines common between files will have to be edited The header information also may have to be edited v.patch complete. 2 vector maps patched Patching vector map qu...@andresfr@andresfr... Building topology for vector map test... Registering primitives... #OUTPUT OF V.OVERLAY Copying vector objects from vector map qu...@andresfr... Copying vector objects from vector map qu...@andresfr... Collecting input attributes... Collecting input attributes... Building partial topology... Building topology for vector map test2... Registering primitives... 26 primitives registered 906 vertices registered Number of nodes: 24 Number of primitives: 26 Number of points: 0 Number of lines: 0 Number of boundaries: 26 Number of centroids: 0 Number of areas: - Number of isles: - Breaking lines... Removing duplicates... Cleaning boundaries at nodes... 26 areas built Attaching islands... Number of primitives: 37 Number of lines: 0 Querying vector map qu...@andresfr... Number of nodes: 47 Building topology for vector map test2... 47 primitives registered 22 isles built Attaching centroids... Number of nodes: 43 Number of points: 0 Number of areas: 26 Number of isles: 22 v.overlay complete. Number of areas without centroid: 10 Number of lines: 0 Number of boundaries: 31 Number of centroids: 16 Number of primitives: 47 Attaching islands... 26 areas built 925 vertices registered Building areas... Registering primitives... Rebuilding topology... Number of areas: 26 Number of isles: 22 Number of areas without centroid: 10 Number of primitives: 53 Number of points: 0 Number of lines: 0 Number of boundaries: 37 Number of centroids: 16 Attaching centroids... Building topology for vector map test2... Attaching centroids... Writing centroids... Querying vector map qu...@andresfr... Number of boundaries: 37 Number of centroids: 0 Number of areas: 26 Number of isles: 22 Number of areas without centroid: 26 Number of points: 0 Number of nodes: 31 22 isles built Attaching islands... Building topology for vector map test2... Building areas... Markus Metz-2 wrote: incanus wrote: By the way, i don't know if what i've done is correct and can be a usefull trick, or it's a rookie mistake, or a casuality. I hope someone knows why this happens, and could help anyone. This is: After trying the v.patch, v.clean(break,rmdupl), v.build thing, v.build is not needed after v.clean. it was the moment to create a new column for all geometries, in order to run v.dissolve. Then, v.dissolve doesn't worked. I see that v.info gives different information (number of boundaries, areas, islands...) of the same map, obtained from v.overlay's or v.patch,v.clean... Can you post that information, or better the output of v.build for each map (not all the input maps, only the results of v.patch and v.overlay)? Look out for number of incorrect boundaries, number of duplicate centroids and number of centroids outside area, that means that there is still something wrong with vector topology. Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- View this message in context: http://n2.nabble.com/V-overlay-very-slow-tp4169783p4176118.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] V.overlay very slow
Um, sorry, I was looking for the output of v.build for the cleaned result of v.patch and for the overlay result, I'm not (yet) interested in the output of v.patch and v.overlay. What I'm looking for is the output of v.build map=patched_and_cleaned_map and v.build map=overlay_map incanus wrote: that's an example with quad3 + quad4 OUTPUT OF V.PATCH Patching vector map qu...@andresfr@andresfr... 42 primitives registered 922 vertices registered Building areas... 26 areas built 24 isles built Attaching islands... Attaching centroids... Number of nodes: 40 Number of primitives: 42 Number of points: 0 Number of lines: 0 Number of boundaries: 26 Number of centroids: 16 Number of areas: 26 Number of isles: 24 Number of areas without centroid: 10 Intersections at borders will have to be snapped Lines common between files will have to be edited The header information also may have to be edited v.patch complete. 2 vector maps patched Patching vector map qu...@andresfr@andresfr... Building topology for vector map test... Registering primitives... #OUTPUT OF V.OVERLAY Copying vector objects from vector map qu...@andresfr... Copying vector objects from vector map qu...@andresfr... Collecting input attributes... Collecting input attributes... Building partial topology... Building topology for vector map test2... Registering primitives... 26 primitives registered 906 vertices registered Number of nodes: 24 Number of primitives: 26 Number of points: 0 Number of lines: 0 Number of boundaries: 26 Number of centroids: 0 Number of areas: - Number of isles: - Breaking lines... Removing duplicates... Cleaning boundaries at nodes... 26 areas built Attaching islands... Number of primitives: 37 Number of lines: 0 Querying vector map qu...@andresfr... Number of nodes: 47 Building topology for vector map test2... 47 primitives registered 22 isles built Attaching centroids... Number of nodes: 43 Number of points: 0 Number of areas: 26 Number of isles: 22 v.overlay complete. Number of areas without centroid: 10 Number of lines: 0 Number of boundaries: 31 Number of centroids: 16 Number of primitives: 47 Attaching islands... 26 areas built 925 vertices registered Building areas... Registering primitives... Rebuilding topology... Number of areas: 26 Number of isles: 22 Number of areas without centroid: 10 Number of primitives: 53 Number of points: 0 Number of lines: 0 Number of boundaries: 37 Number of centroids: 16 Attaching centroids... Building topology for vector map test2... Attaching centroids... Writing centroids... Querying vector map qu...@andresfr... Number of boundaries: 37 Number of centroids: 0 Number of areas: 26 Number of isles: 22 Number of areas without centroid: 26 Number of points: 0 Number of nodes: 31 22 isles built Attaching islands... Building topology for vector map test2... Building areas... Markus Metz-2 wrote: incanus wrote: By the way, i don't know if what i've done is correct and can be a usefull trick, or it's a rookie mistake, or a casuality. I hope someone knows why this happens, and could help anyone. This is: After trying the v.patch, v.clean(break,rmdupl), v.build thing, v.build is not needed after v.clean. it was the moment to create a new column for all geometries, in order to run v.dissolve. Then, v.dissolve doesn't worked. I see that v.info gives different information (number of boundaries, areas, islands...) of the same map, obtained from v.overlay's or v.patch,v.clean... Can you post that information, or better the output of v.build for each map (not all the input maps, only the results of v.patch and v.overlay)? Look out for number of incorrect boundaries, number of duplicate centroids and number of centroids outside area, that means that there is still something wrong with vector topology. Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] How to compile GRASS 6.3 with gstat 2.4.4 on Mac OS X?
Hello GRASS-Users, i would like to build gstat 2.4.4 on Mac OS X 10.5.8 with GRASS 6.4 support (gstat is downloaded from http://www.gstat.org/download.html). I've searched for manuals or mailing lists in the web, but i've find nothing. Can someone help me? Please note, I am a greenhorn! Thank you in advance, Andreas___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Python problem
Hi! I'm trying to develop some pyton code without luck(Grass64 on a Ubuntu 9.04 python 2.6) I start using the example from the wiki http://grass.osgeo.org/wiki/GRASS_and_Python, but I got error messages... Traceback (most recent call last): File /usr/lib/grass64/etc/wxpython/gui_modules/menuform.py, line 1931, in module GrassGUIApp( grassTask( sys.argv[1] ) ).MainLoop() File /usr/lib/grass64/etc/wxpython/gui_modules/menuform.py, line 335, in __init__ xml.sax.parseString( getInterfaceDescription( grassModule ) , processTask( self ) ) File /usr/lib/grass64/etc/wxpython/gui_modules/menuform.py, line 1748, in getInterfaceDescription raise IOError, _(Unable to fetch interface description for command '%s'.) % cmd IOError: Unable to fetch interface description for command 'test.py'. Any hint? Thanks -- Lic. Leonardo A. Hardtke Laboratorio de Teledetección y S.I.G. Centro Nacional Patagónico (CONICET) Bvd. Brown 2825, 9120 Puerto Madryn, Chubut, Argentina ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.out.ascii output modification
Maxime Phaneuf wrote: As the title says, I'm trying to use r.profile with the output of r.drain. I know I cannot use directly the output, and instead I need to transform r.drain to a vector line and transform it to a ascii point file (v.out.ascii). So, when the output form r.drain has been transformed to a ascii point file, I cannot use it because there is a third column (the category) preventing me from piping it into r.profile. I tried deleting the category but this doesn't work. I can always use awk to delete characters in the ascii file that I don't need, but I was wondering if there is a quicker way to do it. Basically what I need is to get from this ascii point file: X Y CAT xxx|yyy|1 xxx|yyy|1 to: X Y xxx|yyy without using a awk-like method. Is it possible? cat $xyfile | cut -f1,2 -d'|' Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Python problem
Leonardo Hardtke wrote: I'm trying to develop some pyton code without luck(Grass64 on a Ubuntu 9.04 python 2.6) I start using the example from the wiki http://grass.osgeo.org/wiki/GRASS_and_Python, but I got error messages... IOError: Unable to fetch interface description for command 'test.py'. Any hint? Do the following work: test.py --help test.py --interface-description -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user