[GRASS-user] Re: Bug or my fault?
Maybe we should all switch to OSX The problem so far appears to be either specific to my own machine setup, or specific to Windows and Grass 6.4.0/6.5. The details are all in this thread. -- View this message in context: http://n2.nabble.com/Bug-or-my-fault-tp4609390p4635889.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
[GRASS-user] Re: Bug or my fault?
Thanks Marcus, I think our messages crossed, here is the result (see my previou message): I have just installed 6.5.svn-r41167-1 from the daily Win releases, into my Win 7. It crashes at start-up with: C:\programs\GRASS-65>grass65svn.bat Traceback (most recent call last): File "C:/programs/GRASS-65/etc/wxpython/gis_set.py", line 37, in from gui_modules import help File "c:/osgeo4w/usr/src/grass6_devel/dist.i686-pc-mingw32/etc/wxpython/gui_mo dules/help.py", line 25, in ImportError: No module named agw.customtreectrl Traceback (most recent call last): File "C:/programs/GRASS-65/etc/wxpython/wxgui.py", line 63, in import wx.lib.agw.customtreectrl as CT ImportError: No module named agw.customtreectrl C:\programs\GRASS-65> -- View this message in context: http://n2.nabble.com/Bug-or-my-fault-tp4609390p4634804.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
[GRASS-user] Re: Bug or my fault?
Some more info. My main platform is Win7 Pro, GRASS 6.4.0 standalone, Python 2.6.4. But I have tried other platforms with identical setup, with the following results: 1. All works well under Ubuntu 9.10. 2. But problem exists also in WinXP so it is not Win 7 specific, but it is Windows specific. You don't need a lot of vector files. You can have a single vector, then go to the dbf directory and copy the single .dbf file in that same directory (Windows will automatically give it a new name). Then copy all again in the same directory and so on, after only 8 iterations you will have 256 files. Open the original vector file in Grass and try to list the attribute table. Unless there is something very specifically wrong with my setup, you should experience the same problem I have. And yes it is slow, very slow in fact, compared to Ubuntu which is almost instantaneous in its response, and Win 7 is even slower than XP running as a virtual machine! I don't have grass 6.5 installed, so I cannot try it. If ther are Win binaries available for it I will give it a go, but I need some more time. Thanks again for your help -- View this message in context: http://n2.nabble.com/Bug-or-my-fault-tp4609390p4629131.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
[GRASS-user] Re: Bug or my fault?
Thanks for your help Michael. Here is the status. Issued v.db.connect and all prints as expected by you and all vectors have a corresponding .dbf file. This is what I do. I have a mapset with 270 vectors. I open GRASS64 with wxpython GUI and select one vector in the middle, call it vect135. I right click to list the attributes. The system hangs with error "Python.exe not responding" I restart GRASS64 with TclTk GUI and load the same vector, press the attributes button and all attributes show up perfectly. I then delete the top 100 .dbf files from the dbf directory only (not from the vector directory). I restart GRASS64 with the python GUI again, load the same vector (vect135), right click and the attribute table shows up perfectly. I then restore the 100 .dbf files I had previously deleted and now delete the bottom 100. Restart GRASS64 with python GUI, load vect135, right click and the attribute table shows up perfectly. As soon as the number of dbf files increases to approx 188 (the actual number seems to depend on which vectors are involved) the system hangs when opening the attribute table, but only with the python GUI. I have now restructured my location so I have a number of mapsets and have distributed the vectors among the new mapsets, so they all have less than 100 files each, and all works perfectly. But I am really puzzled and would dearly like to know why this happens only to me! By the way, the problem occurs running GRASS with or without MSys. -- View this message in context: http://n2.nabble.com/Bug-or-my-fault-tp4609390p4623610.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
[GRASS-user] Re: Bug or my fault?
Thanks Michael but I don't think this is the problem in my case, I am using the default DBF database, no DB connections required. Also remember that it all works well when I use TclTk as the GUI, and all works well when the number of DBF files in the mapset is less than approx 188. I haven't been able to find any other cause so far, but it could be Python related rather than GRASS, when I kill the GUI that hangs I get an error message to the effect that Python.exe is not responding. But why should the number of DBF files be relevant??? -- View this message in context: http://n2.nabble.com/Bug-or-my-fault-tp4609390p4614994.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
[GRASS-user] Re: Bug or my fault?
The only culprit I could find so far is the number of vector or dbf files in the Location. If the number is greater than 188 then the wxPython GUI goes into an infinite loop. Reduce it to 188 and it works as it should. It does not matter which ones are the 188 files, as long as there are no more than 188. I am reluctant to suggest this is the cause of the problem, but it certainly looks that way so far. Any more plausible suggestions anybody? -- View this message in context: http://n2.nabble.com/Bug-or-my-fault-tp4609390p4610112.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
[GRASS-user] Bug or my fault?
Under certain conditions, which I have not yet been able to pinpoint, the python gui hangs forever when right clicking on vector map to list the attribute table. Change the gui to tcltk and the problem disappears, the attributes for the same map are listed correctly. I am using WINGRASS 6.4.0 standalone, stable release, under Win 7 Pro. and Python 2.6.4 I will persevere to try and pinpoint the problem a bit better, but I would really like to know if anybody else has had a similar experience. Ownership and permissions of all the files involved appear to be ok, although one never knows for sure with Win 7! thanks -- View this message in context: http://n2.nabble.com/Bug-or-my-fault-tp4609390p4609390.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
[GRASS-user] New Wiki entry in the FAQs
For those of you who are interested, I have created a new FAQ entry in the Wiki describing how to use GRASS from Python using the Eclipse IDE. Can be found here: http://grass.osgeo.org/wiki/Using_Eclipse_to_develop_GRASS_Python_programs Cheers George -- View this message in context: http://n2.nabble.com/New-Wiki-entry-in-the-FAQs-tp2879085p2879085.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
[GRASS-user] GRASS, Eclipse, Python behavioural problems in Linux9.04 vs WinXP
Hi, I am suffering from few inconsistencies when running the same software levels of Grass (6.4.rc4), Eclipse (3.4.2) and Pydev 1.4.5., Python 2.6 on a WinXP platform vs Ubuntu 9.04, as follows. Using Eclipse to run or debug this short Python program: import os import grass vectorname="map2000vect" type="vector" mapset="PERMANENT" env = grass.gisenv() print env res = grass.find_file(vectorname,type,mapset) print res rem = grass.read_command('g.remove', vect=vectorname) print rem behaves faultlessly under WinXP but when run under Ubuntu although it runs without errors, only the first two commands are actually executed and produce output, the g.remove appears to be executed but results in no removal and no output. The same source program runs correctly however when executed in Ubuntu as a python script from the Grass terminal. I was hoping that some Grass user out there that has tried a similar software combination might be able to throw some light even though I am not sure this really is a Grass problem. Many thanks -- View this message in context: http://n2.nabble.com/GRASS%2C-Eclipse%2C-Python-behavioural-problems-in-Linux9.04-vs-WinXP-tp2737691p2737691.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] Splitting shapefile
You could use v.mkgrid to generate a grid. Then loop over all categories in this grid (v.categories) and use v.overlay to intersect. Thank you Markus, I had no problem in creating the grid with: v.mkgrid --overwrite map=grid_vect grid=20,20 box=1,1 (only struck a small bug in the GUI which does not let you enter two values for the grid (eg 20,20) so I had to enter the command in the command box). But had no success at all using v.category and v.overlay to create individual map tiles for each of the 400 tiles in the grid. I suspect I am not using v.category properly but could not find an example that could help me. My understanding is that for v.overlay to work, the b-input must be a single tile, then the AND intersection will extract the a-input features for the area under that tile. But I could not find a way to extract a single tile out of the grid. What am I missing? Do I need to reduce the current region to the single tile before v.overlay and then loop through all the tiles that way? George -- View this message in context: http://n2.nabble.com/Splitting-shapefile-tp2536771p2551084.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
[GRASS-user] Splitting shapefile
Hi , I have a large vector shapefile I would like to split into many smaller tiles of equal size. Can somebody suggest a way to do this in GRASS 6.4 (either Windows or Linux). Also is it possible to use scripts (for batch processing) under the Windows version of GRASS? Many thanks George -- View this message in context: http://n2.nabble.com/Splitting-shapefile-tp2536771p2536771.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] d.rast.edit does not save edited map although Save button is pressed
thanks Glynn, here is what I am doing. After I type the command two new windows open up: one a grid with a faint outline of the input raster maps, the other what I consider to be a navigation window, to help displaying on the grid the area to be edited. I want to extend the track21_raster map line, a walking track. So I click on the cell nearest to the last colored cell in the grid and enter the two category values I have, in the bottom grid display boxes: one called value (I enter a 1, the track segment value), the other aspect (I enter the elevation, 640) then press enter. (There is a third box called also Value, which I leave blank as I don't know what to place there!). A new whitish cell is created. I repeat the process for a few cells then save, but nothing is saved. The output map is identical to the input. The following is displayed in the terminal: WARNING: Raster map not found WARNING: nothing removed Is it possible that the new cells are painted in white therefore invisible when seen in the display? I have not found a way to change their color. Any further help you can give will be most welcome. By the way I am on GRASS 6.3 and Ubuntu 8.04. Glynn Clements wrote: > > > georgew wrote: > >> Hi, I must be doing something wrong with the following: >> >> d.rast.edit input=track21_r...@permanent output=track21_edited_rast >> aspect=contours21_r...@permanent width=640 height=480 size=12 rows=200 >> cols=200 --overwrite >> >> The edited data is not saved in the output map which remains identical to >> the input map, although I pressed the Save button in the grid display. >> All the data is in my home directory (and anyway using "sudo >> d.rast.edit" >> results in "d.rast.edit not found" error.) >> Any clues what I might be doing wrong. ? > > What is the nature of the editing? If the only changes are setting > specific cells to null, that's one thing which d.rast.edit can't > handle. You would need to set the cells to some unused value then > subsequently change them to nulls with r.null. > > If you're changing cells to non-null values, then I don't know what > the problem is. > > -- > Glynn Clements > ___ > grass-user mailing list > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user > > -- View this message in context: http://www.nabble.com/d.rast.edit-does-not-save-edited-map-although-Save-button-is-pressed-tp21011399p21020243.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
[GRASS-user] d.rast.edit does not save edited map although Save button is pressed
Hi, I must be doing something wrong with the following: d.rast.edit input=track21_r...@permanent output=track21_edited_rast aspect=contours21_r...@permanent width=640 height=480 size=12 rows=200 cols=200 --overwrite The edited data is not saved in the output map which remains identical to the input map, although I pressed the Save button in the grid display. All the data is in my home directory (and anyway using "sudo d.rast.edit" results in "d.rast.edit not found" error.) Any clues what I might be doing wrong. ? Many thanks George -- View this message in context: http://www.nabble.com/d.rast.edit-does-not-save-edited-map-although-Save-button-is-pressed-tp21011399p21011399.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
[GRASS-user] Automatic calculation of coordinates, elevation and distance from context region.
Hi, some background to my problem: I have downloaded a shapefile containing a walking track with coordinates and no elevation, and downloaded another contours shapefile with coordinates and elevation, for the same region as the track. Thanks to this forum I can now extract the walking track elevation (and distance between points) from the contour lines at the point where the track intersects with the contours. I now would like to extend the walking track, preferably through a GUI, by creating new points or lines, whose coordinates and elevation should be calculated by the GUI from the context in which the track is placed, and be able to process the resulting track, i.e. extract elevation (and distance between points), as though they were part of the original shape file. Is this possible with GRASS and any suggestions how best to achieve this? I hope I have explained this clearly enough, if not please let me know. Many thanks, -- View this message in context: http://www.nabble.com/Automatic-calculation-of-coordinates%2C-elevation-and-distance-from-context-region.-tp20926649p20926649.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] Problem trying to create new wiki FAQ
The new FAQ is now in place: http://grass.osgeo.org/wiki/How_to_create_an_elevation_profile#A: Thanks Martin. Martin Landa-2 wrote: > > Hi, > > 2008/12/6 georgew <[EMAIL PROTECTED]>: >> Hi, I tried several times to obtain edit permission to the FAQs so as to >> insert a new one but without success. The error is always the same, my >> email >> is not validated. I tried for a couple of days to have it validated but I >> never receive the validation email. (The email address I supplied on >> registration is correct). What else can I do?... > > when you create wiki account you should receive authenticate e-mail. > Anyway now your account should be activated. Let me know in the case > of problems. > > Martin > > -- > Martin Landa * http://gama.fsv.cvut.cz/~landa * > ___ > grass-user mailing list > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user > > -- View this message in context: http://www.nabble.com/Problem-trying-to-create-new-wiki-FAQ-tp20867557p20878725.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
[GRASS-user] Problem trying to create new wiki FAQ
Hi, I tried several times to obtain edit permission to the FAQs so as to insert a new one but without success. The error is always the same, my email is not validated. I tried for a couple of days to have it validated but I never receive the validation email. (The email address I supplied on registration is correct). What else can I do?... -- View this message in context: http://www.nabble.com/Problem-trying-to-create-new-wiki-FAQ-tp20867557p20867557.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] Elevation profi le from intersecting shapefiles
Much progress! Micha Silver wrote: > > Set resolution to 10m. X 10m. with > g.region -s -p res=10 > Hamish_b wrote: > > Multiple displays would need multiple > WIND files but the processing module (eg v.to.rast) wouldn't know which > was the appropriate one to use. Thus computational region refers to what > is stored in the mapset's WIND file. > Moritz Lennert-2 wrote: > > You are aware that this means that the Contour.shp did not contain any > projection info and that GRASS thus consideres it as unprojected ? > Setting the resolution BEFORE converting to raster was the answer, thank you. I now see the contours in the raster map as well as in the original vector file. As for the null projection, thanks for pointing it out (I was not aware!) but since the purpose of the exercise was to create an elevation profile from tabular data perhaps it is not critical. So I went ahead with the original purpose, i.e. to create an elevation profile for the track. Micha Silver wrote: > > r.surf.contour to create an elevation surface from the rasterized contours > r.profile to get the elevations along the tracks. > So I entered: r.surf.contour [EMAIL PROTECTED] output=elevmap_raster --overwrite and few seconds later it completed its job without errors. Then went ahead with: r.profile -i [EMAIL PROTECTED] output=- null=* I followed the track with the left mouse button and a two column list of values was printed on stdout, the distance and the elevation. Success! My only remaining question is: Is there a way to use the track raster file to input the values instead of using the mouse, so that the process can be more or less automated? Thank you all for your great help. I'll write a summary and post it in the Wiki as an FAQ or short tutorial as suggested by Markus Neteler. grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- View this message in context: http://www.nabble.com/Elevation-profile-from-intersecting-shapefiles%E2%80%8F-tp20744104p20844940.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] Elevation profi le from intersecting shapefiles
Hamish, please forgive my ignorance but I am not sure what you mean. As far as I know g.region is the means to change the resolution, whether through the command line or the GUI. Can you please explain how to access and change the computational region settings. hamish_b wrote: > > [please don't post in html, it makes it hard to read] > > > georgew wrote: >> PS I tried to play around with resolution in the Region dialog but the >> display never changed from the same chunky coloured squares. > > > you need to change the computational region, not just the display region > settings. the computational region is what the processing modules will > use. > > > Hamish > > > > > > ___ > grass-user mailing list > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user > > -- View this message in context: http://www.nabble.com/Elevation-profile-from-intersecting-shapefiles%E2%80%8F-tp20744104p20828846.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] Elevation profi le from intersecting shapefiles
Thanks Micha, here we go: Micha Silver wrote: > > How many shapefiles did you have in the GRASSDATA directory? > I had the two shapefiles, Contour and Track, but I got your message loud and clear about creating separate vectors, which I did and displayed as follows: v.in.ogr -o dsn=/home/george/GRASSDATA/Contour.shp output=contours min_area=0.0001 type=line snap=-1 --overwrite v.in.ogr -o dsn=/home/george/GRASSDATA/Track_cl.shp output=tracks min_area=0.0001 type=line snap=-1 --overwrite d.vect [EMAIL PROTECTED] color=0:0:0 lcolor=0:0:0 fcolor=170:170:170 display=shape type=point,line,boundary,area icon=basic/circle size=5 layer=1 lsize=8 xref=left yref=center llayer=1 d.vect [EMAIL PROTECTED] color=255:0:0 lcolor=0:0:0 fcolor=170:170:170 display=shape type=point,line,boundary,area icon=basic/circle size=5 width=3 layer=1 lsize=8 xref=left yref=center llayer=1 No problems so far, all displayed correctly. Micha Silver wrote: > > At this point I'd check to make sure the elevation data came thru: > v.info -c [EMAIL PROTECTED] Displaying column types/names for database connection of layer 1: INTEGER|cat DOUBLE PRECISION|ID CHARACTER|CRT_DATE CHARACTER|MOD_DATE DOUBLE PRECISION|ELEVATION CHARACTER|DESIGNATED CHARACTER|NAT_FORM CHARACTER|DEFINITION v.info -c [EMAIL PROTECTED] Displaying column types/names for database connection of layer 1: INTEGER|cat DOUBLE PRECISION|ID CHARACTER|NAME CHARACTER|CRT_DATE CHARACTER|MOD_DATE CHARACTER|TRACK_USE CHARACTER|TRACK_TYPE CHARACTER|TRK_STATUS DOUBLE PRECISION|area DOUBLE PRECISION|len v.db.select [EMAIL PROTECTED] layer=1 column=ELEVATION fs=| ELEVATION 1180 1220 1200 etc... etc... So far all is as it should be (I think..) Convert vector to raster: (no problems reported) v.to.rast [EMAIL PROTECTED] output=contours_rast use=attr type=line layer=1 column=ELEVATION value=1 rows=10 --overwrite (... Converted points/lines: 630 of 630) Display raster map (PROBLEM:only displays large, brightly coloured rectangles, no contours d.rast [EMAIL PROTECTED] -o How did you set up the Location/Mapset ? You might have a problem with the resolution. What do you get from: r.info testrastmap g.region -p /quote> I set the Location in GRASS's opening welcome screen, using the Contour.shp file to provide the boundaries. Here is the result: r.info contours_rast ++ | Layer:contours_rast Date: Thu Dec 4 08:50:19 2008 | | Mapset: test02 Login of Creator: george | | Location: test02 | | DataBase: /home/george/GRASSDATA | | Title:Labels ( contours_rast ) | | Timestamp: none | || | | | Type of Map: raster Number of Categories: 0 | | Data Type:DCELL | | Rows: 20 | | Columns: 20 | | Total Cells: 400 | |Projection: x,y | |N: 6020724.5S: 6013982.586 Res: 337.0957 | |E: 2481099.136W: 2467034.878 Res: 703.2129 | | Range of data:min = 520.00 max = 1840.00 | | | | Data Source: | |Vector Map: [EMAIL PROTECTED] in mapset test02 | |Original scale from vector map: 1:1 | | | | Data Description: | |generated by v.to.rast | | | | Comments: | |v.to.rast input="[EMAIL PROTECTED]" output="contours_rast" use="attr" \ | |type="line" layer=1 column="ELEVATION" value=1 rows=10 | | | ++ GRASS 6.3.0 (test02):~ > g.region -p projection: 0 (x,y) zone: 0 north: 6020724.5 south: 6013982.586 west: 2467034.878 east:
Re: [GRASS-user] Elevation profi le from intersecting shapefiles
Dylan, using v.surf.rst gives me "Segmentation fault", using: v.surf.rst [EMAIL PROTECTED] layer=1 zcolumn=ELEVATION elev=kahu-surf-rast tension=40. segmax=40 npmin=300 dmin=168.547850 dmax=842.739250 zmult=1.0 --overwrite I am on GRASS 6.3, Ubuntu 8.04, AMD64 Any suggestions please? Maris, I cannot find v.profile on my GRASS6.3! Regards Dylan Beaudette-2 wrote: > > Dude- It is saturday night. Your request is for a fairly simple > operation, but there are several steps. > > Some hints: > > 1. convert your contours into a DEM: > v.surf.rst > -- View this message in context: http://www.nabble.com/Elevation-profile-from-intersecting-shapefiles%E2%80%8F-tp20744104p20808111.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] Elevation profi le from intersecting shapefiles
Micha, I followed your instructions with the following results: Import data from 2 shapefiles, Contour (layer 1)and Track (layer 2): (no problems reported) v.in.ogr dsn=/home/george/GRASSDATA/ output=testvectmap min_area=0.0001 snap=-1 Display vector map: (no problems, map displayes correctly) d.vect [EMAIL PROTECTED] color=0:0:0 lcolor=0:0:0 fcolor=170:170:170 display=shape type=point,line,boundary,area icon=basic/circle size=5 layer=-1 lsize=8 xref=left yref=center llayer=1 Convert vector to raster and thin: (no problems reported) v.to.rast [EMAIL PROTECTED] output=testrastmap use=attr type=point,line,area layer=1 column=ELEVATION value=1 rows=10 r.thin [EMAIL PROTECTED] output=thinrastmap iterations=200 Display raster map (PROBLEM:only displays large, brightly coloured rectangles, no contours d.rast [EMAIL PROTECTED] -o At this point I stopped as I got the impression something was already wrong but could not work out what. Also, should I not rasterize the track data as well? Any help would be appreciated. If you require the source shapefiles I have no problems sending them to you. Micha Silver wrote: > > > Yes. > I think you should be looking at converting the shapefile contour lines > to a raster > > v.in.org to import the contour lines and tracks into grass > v.to.rast and r.thin to convert the contour lines to a raster map. > r.surf.contour to create an elevation surface from the rasterized > contours > r.profile to get the elevations along the tracks. > > Post back with specific questions if you hit a snag. > > Regards, > Micha > > ___ > grass-user mailing list > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user > > -- View this message in context: http://www.nabble.com/Elevation-profile-from-intersecting-shapefiles%E2%80%8F-tp20744104p20803902.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] Elevation profi le from intersecting shapefiles
I am indebted to all of you for taking the time to answer. Will follow up on all your suggestions. Season's greetings to all.... georgew wrote: > > Hi, I am totally new to GRASS and overwhelmed by its scope and > -- View this message in context: http://www.nabble.com/Elevation-profile-from-intersecting-shapefiles%E2%80%8F-tp20744104p20761847.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] Elevation profi le from intersecting shapefiles
What, no takers? Just a Yes or No will be already of great help to me, anything more, like some hint on a possible solution would be cream on top. Thanks all the same for reading this. George georgew wrote: > > Hi, I am totally new to GRASS and overwhelmed by its scope and not sure > whether the following is possible: > I have a set of topographic maps in shapefile format containing contours > and walking tracks. The track data does not include elevation (only lat. > and long.), but the contour lines do include elevation. Is it possible > (and how?) to create an elevation profile for the length of the walking > track using the intersecting points between the track and the contour > lines to obtain the elevation values? > Many thanks for any assistance you can give me. > George > -- View this message in context: http://www.nabble.com/Elevation-profile-from-intersecting-shapefiles%E2%80%8F-tp20744104p20754964.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