[GRASS-user] Fwd: [GRASS-dev] New WinGrass Release
For your information... Note: This is a separate package from the OSGeo4W installation (which is packaged by someone else). -- Forwarded message -- From: Colin Nielsen colin.niel...@gmail.com Date: Tue, Jun 2, 2009 at 9:14 PM Subject: [GRASS-dev] New WinGrass Release To: grass-wind...@lists.osgeo.org, grass-dev grass-...@lists.osgeo.org I've just put together a new native WinGrass package based on releasebranch_6_4 r37703 (almost a month's worth of updates). Download here: http://grass.osgeo.org/download/index.php#g64x Please test and report any problems. -Colin ___ grass-dev mailing list grass-...@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-dev ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Icelandic / international characters in databases and ps.map
Hi, After initially being put off by seeing complicated scripts, I've started to use psmap. I wanted to use the postscript fill patterns to create a geological map that would reproduce in black and white. It's really good - the commands are pretty easy and the results look clear and excellent. Anyway. I have digitised some of the main peaks in the region, whose names contain Icelandic characters. Their names contain Icelandic characters. Using v.db.select, they appear to have stored correctly. When I use labels created from the vector, the characters appear as gaps with d.labels and as nonsense symbols using ps.map. I read on the psmap g.manual page about character encoding, but don't know how it applies to the labels/database. Can anyone help? Cheers John -- Dr John Stevenson Postdoctoral Research Associate School of Earth, Atmospheric and Environmental Sciences Williamson Building (Room 2.42) University of Manchester Manchester M13 9PL, UK tel. +44(0)161 306 6585; fax. +44(0)161 306 9361; john.steven...@manchester.ac.uk ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.in.ogr problem
Hello Maris, Hello all! Problem solved (I have to put this on topic?) I'm not sure if the file I've tested is so big to GRASS, but I left the v.in.ogr working for a hole afternoon and nothing. So I tried a smaller data and still nothing. Then I did some search on the net and I found a tip that solved my problem. At QGIS forum I found a topic discussing the same problem (http://forum.qgis.org/viewtopic.php?f=2t=4451#p7490). A post there send me to a Jachym Cepicky's post at another forum ( http://n2.nabble.com/grass-gdal-ubuntu-issue-solved-(-)-td2352498.html) with the solution. Well, I reintalled GRASS by Jachim's repository, and the v.in.ogr worked fine. I have also tested this tool in QGIS (grass plugin) and there was no problem. By the way, grass imported that big data in a glimpse! : ) Cheers! Rafael Hello. How many vector features are in that file? Could it be that v.in.ogr is still processing input file? On some files it may take long time to import vector data. Maris. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Icelandic / international characters in databases and ps.map
John A Stevenson wrote: After initially being put off by seeing complicated scripts, I've started to use psmap. I wanted to use the postscript fill patterns to create a geological map that would reproduce in black and white. It's really good - the commands are pretty easy and the results look clear and excellent. Anyway. I have digitised some of the main peaks in the region, whose names contain Icelandic characters. Their names contain Icelandic characters. Using v.db.select, they appear to have stored correctly. When I use labels created from the vector, the characters appear as gaps with d.labels and as nonsense symbols using ps.map. I read on the psmap g.manual page about character encoding, but don't know how it applies to the labels/database. Can anyone help? First, any text passed to ps.map must be in ISO-8859-1; other encodings won't work (although if your PostScript interpreter has suitable fonts, you can fix this by manually editing prolog.ps). If you use FreeType fonts, d.* commands will work with any encoding known to iconv (from the command line, you need to specify the encoding via d.font; I don't recall how the GUI deals with this). -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] wxpython gui and v.digit
Hi, I have installed from source Grass65 devel_branch today. It compiled without problems. When I run g.gui wxpython, it gives me the following warning: WARNING: Vector digitizer is not available (No module named grass6_wxvdigit) How do I enable this module? Thanks, John -- Dr John Stevenson Postdoctoral Research Associate School of Earth, Atmospheric and Environmental Sciences Williamson Building (Room 2.42) University of Manchester Manchester M13 9PL, UK tel. +44(0)161 306 6585; fax. +44(0)161 306 9361; john.steven...@manchester.ac.uk ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Icelandic / international characters in databases and ps.map - SOLVED
Hi, Solved. Firstly, I had to compile grass with freetype (which I didn't have) and freetype includes. ./configure ... --with-freetype --with-freetype-includes=/usr/include/freetype2/ I backed up my vector (as peaks_iso), then, as described in the thread, I opened the its dbf file in Open Office using character encoding ISO-8859-1. I had to retype the Icelandic characters into the database, and saved it. I remade the labels using peaks_iso and it worked fine. Thank you. John Miha Staut wrote: This thread might help (see the ending comment by Glynn). http://n2.nabble.com/Character-encoding-problem-td2185961.html#a2191493 Miha --- On Wed, 3/6/09, John A Stevenson john.steven...@manchester.ac.uk wrote: From: John A Stevenson john.steven...@manchester.ac.uk Subject: [GRASS-user] Icelandic / international characters in databases and ps.map To: GRASS user list grass-user@lists.osgeo.org Date: Wednesday, 3 June, 2009, 10:44 AM Hi, After initially being put off by seeing complicated scripts, I've started to use psmap. I wanted to use the postscript fill patterns to create a geological map that would reproduce in black and white. It's really good - the commands are pretty easy and the results look clear and excellent. Anyway. I have digitised some of the main peaks in the region, whose names contain Icelandic characters. Their names contain Icelandic characters. Using v.db.select, they appear to have stored correctly. When I use labels created from the vector, the characters appear as gaps with d.labels and as nonsense symbols using ps.map. I read on the psmap g.manual page about character encoding, but don't know how it applies to the labels/database. Can anyone help? Cheers John -- Dr John Stevenson Postdoctoral Research Associate School of Earth, Atmospheric and Environmental Sciences Williamson Building (Room 2.42) University of Manchester Manchester M13 9PL, UK tel. +44(0)161 306 6585; fax. +44(0)161 306 9361; john.steven...@manchester.ac.uk ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Dr John Stevenson Postdoctoral Research Associate School of Earth, Atmospheric and Environmental Sciences Williamson Building (Room 2.42) University of Manchester Manchester M13 9PL, UK tel. +44(0)161 306 6585; fax. +44(0)161 306 9361; john.steven...@manchester.ac.uk -- Dr John Stevenson Postdoctoral Research Associate School of Earth, Atmospheric and Environmental Sciences Williamson Building (Room 2.42) University of Manchester Manchester M13 9PL, UK tel. +44(0)161 306 6585; fax. +44(0)161 306 9361; john.steven...@manchester.ac.uk ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Deleting vertices based on a list of points
On Jun 3, 2009, at 4:50 AM, Moritz Lennert wrote: On 02/06/09 20:50, Dwight Needels wrote: Hi all, I have two vectors: Vector A has lines and Vector B has points that share the same coordinates as a subset of the vertices in the lines. I would like to delete all of the vertices/nodes in Vector A that lie directly underneath the points in Vector B (or perhaps lie within a threshold). I suspect that v.edit tool=vertexdel can do this, but I cannot figure out how to pass the coordinates of a set of points from Vector B to v.edit to be used to select vertices in Vector A for deletion. Could somebody give me a hand with this (or tell me it is not possible)? Feed the output of v.to.db -p option=coor on Vector B into the coords= parameter of v.edit on Vector A, i.e. something like this: v.edit A tool=vertexdel coords=`v.to.db -p B option=coor --quiet | awk -F; '{printf%f,%f,,$2,$3} END{printf\n}'` Moritz Thanks! This worked great with two small changes... using | instead of ; for the field terminator and specifying layer=2 because the points were generated by v.to.points, which puts the point data on layer 2. -Dwight ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Forcing polyline nodes through a point
On May 28, 2009, at 10:04 PM, Dwight Needels wrote: I am trying to create a map of hiking trails based on GPS data. If anyone can point me to a tutorial that discusses the techniques and issues involved, I would appreciate it. Meanwhile, a couple of specific questions. I have a series of GPS tracks that have been brought into GRASS as polylines, and a series of time-averaged waypoints that have been brought into GRASS as points. The lines represent segments of color- blazed trails, and the waypoints represent higher-quality measurements of trail termini and intersections. My questions were about how to force a vector line to pass exactly through specified points in such a way that they would continue to pass through those points after smoothing. Thank you Hamish for pointing out several very useful tools for me to concentrate on learning. Also, thank you Moritz for the awk magic required to pipe point coordinates from one vector into the v.edit command operating on a second vector. With your indispensable help I was able to come up with something that does almost exactly what I want. If anybody can see any gotchas or better ways to accomplish these steps, I would appreciate the feedback. Thanks, -Dwight Starting with a line (track) vector and a point (waypoint) vector, modify the line vector so that it passes exactly through each point. This handles three cases; 1) where a vector line terminus undershoots the defined terminus point, 2) where a vector line terminus overshoots the defined terminus point, and 3) where the vector line passes within the threshold of a point somewhere in the middle of the line (a fork, hard curve or switchback). Threshold (distance from points) should be chosen to include at least two vertices for each line segment of interest, without including vertices from other line segments. Line vector: trail Point vector: trail_nodes Resulting line vector: trail_clean Summary: # Split the line vector at closest spot to each point within a threshold, delete superfluous points (v.net op=connect; v.edit tool=delete) # Delete added lines and nodes to leave a gap (v.to.points -n; v.edit tool=vertexdel) # Bridge the gaps, but have each bridge pass through one of the original points (v.distance -a; v.edit tool=copy) # Remove zero length lines, merge bridges across the gaps in the line vector (v.clean tool=rmline; v.edit tool=merge) # Add back a node at each of the original points, remove zero length lines (v.net op=connect; v.clean tool=rmline) # Clean up intermediate vectors (g.remove) Steps: v.net trail points=trail_nodes out=trail_net op=connect thresh=100 --o v.edit tool=delete map=trail_net layer=2 cats=1- v.to.points -n in=trail_net out=trail_net_nodes --o v.edit trail_net tool=vertexdel coords=`v.to.db -p trail_net_nodes layer=2 option=coor --quiet | awk -F| '{printf%f,%f,,$2,$3} END{printf\n}'` v.distance -p -a from=trail_nodes to=trail_net from_type=point to_type=line output=trail_bridges dmax=100 upload=dist column=dist -- overwrite v.edit map=trail_net tool=copy bgmap=trail_bridges ids=1- v.clean input=trail_net output=trail_net_clean type=line tool=rmline thresh=-1 --o v.edit tool=merge map=trail_net_clean ids=1- v.net trail_net_clean points=trail_nodes out=trail_plus_nodes op=connect thresh=100 --o v.clean input=trail_plus_nodes output=trail_clean type=line tool=rmline thresh=-1 --o g.remove vect = trail_net,trail_net_clean,trail_bridges,trail_net_nodes,trail_plus_nodes Nodes were added back (the polyline was split at each of the original points) so that the trail could be smoothed while still guaranteeing that the vector will pass through the original points. # Smooth/straighten (v.generalize, method=snakes) v.generalize input=trail_clean output=trail_smooth method=snakes alpha=1 beta=1 --o ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Averaging multiple vector lines
I have a GRASS vector that originated as multiple GPS tracks from walking a particular trail segment on several different days. Is there a good way to average these lines to get a single line? I want to minimize GPS accuracy errors by averaging across multiple days and also minimize precision errors (random jumping around on a single day) while still maintaining the shape of the trail with all of its twists and turns. I have been able to generate a composite vector by using a combination of v.to.rast, r.grow, r.thin, r.to.vect, v.clean, and v.generalize method=douglas. This method works pretty well when the lines remain close together, but it is very dependent on picking a value for the r.grow radius that fills in all of the gaps between the multiple tracks. If one track is quite different than the others in even a single region of the vector, this requires a relatively large radius value. Moreover, the final vector is located about midway between the two extremes rather than being weighted toward where the majority of tracks fall. It seems like there would be a way to calculate some sort of sliding average of the coordinates that fall within a certain size window, perhaps after using v.to.points with a small dmax (5 ft?) to generate a fairly dense set of points. Ideally, the calculation window could be wider perpendicular to the direction of the line than it is along the direction of the line. From day to day tracks are often within 10 to 20 ft of each other, but it is not uncommon for two tracks to be 30 ft away from each other at some points. Any ideas? -Thanks, -Dwight ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user