[GRASS-user] Problem in v.surf.rst
Hi, As i mentioned in my couple of previous mails, I found that v.surf.rst function does not work for either 2D point or iso-lines unless "Use Z coordinates for approximation" flag is ticked together with the selection of the *name* of the attribute column with values to be used for approximation. When both are chosen, it runs and gives the expected output.However I think there is a minor bug in this function and it should be rectified and made work for 2D data sets without selection of Z flag. The version of GRASS i used is WINGRASS 7 installed through OSGEO4W installer. Cheers, Brian. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Interpolation with v.surf.rst
Yes. there is no point in using Z coordinate if the input is 2D. However, in the current GRASS 7, even if the input is 2D, if -z flag is not used, v.surf.rst does not work. On Wed, Jan 16, 2013 at 4:24 PM, Martin Landa wrote: > Hi, > > 2013/1/16 Markus Metz : > > If you want to use elevation values from zcolumn, you should not use > > the -z flag. Try instead > > v.surf.rst input=isoline@TOPOIKDATA zcolumn=height elev=elev_1m > > I was thinking to omit such -z flags in G7. > > * if zcolumn is given, use values from this column > * if zcolumn is *not* given use z coordinates (if vector is 2D, print > error). > > What do you think? Something similar was done for v.out.ogr (-z -> > !-2) some weeks ago. Motto would be: to use z coordinates if the map > is 3D and the user didn't enter an option or flag which forces to use > values from attribute table. > > Martin > > -- > Martin Landa * http://geo.fsv.cvut.cz/~landa > ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Interpolation with v.surf.rst
Hi, When i omit -z flag, v.surf.rst function keeps on running and nothing happens. Here is what i used. v.surf.rst input=isoline@TOPOIKDATA zcolumn=height elev=elev_1m Brian On Wed, Jan 16, 2013 at 3:24 PM, Markus Metz wrote: > On Wed, Jan 16, 2013 at 10:16 AM, Brian Sanjeewa Rupasinghe > wrote: > > Hi, > > > > I have 2D isolines in GRASS vector format with elevation in an attribute > > field called height. > > I used v.surf.rst function to create raster DEM. However, after > creation, it > > outputs a DEM > > with value 0 elevation only. The description of the function is given as > > below. What would > > be the possible cause for this? I am using wingrass 7. > > > > cheers, Brian. > > > > v.surf.rst -z input=isoline@TOPOIKDATA zcolumn=height elev=elev_1m > > If you want to use elevation values from zcolumn, you should not use > the -z flag. Try instead > v.surf.rst input=isoline@TOPOIKDATA zcolumn=height elev=elev_1m > > Also make sure that the computational region extents and resolution > are set correctly. > > HTH, > > Markus M > ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Interpolation with v.surf.rst
Hi, I have 2D isolines in GRASS vector format with elevation in an attribute field called height. I used v.surf.rst function to create raster DEM. However, after creation, it outputs a DEM with value 0 elevation only. The description of the function is given as below. What would be the possible cause for this? I am using wingrass 7. cheers, Brian. v.surf.rst -z input=isoline@TOPOIKDATA zcolumn=height elev=elev_1m Reading features from vector map ... Ignoring 86730 points (too dense) Processing segments... WARNING: Taking too long to find points for interpolation - please change the region to area where your points are. Continuing calculations... WARNING: Taking too long to find points for interpolation - please change the region to area where your points are. Continuing calculations... WARNING: Taking too long to find points for interpolation - please change the region to area where your points are. Continuing calculations... WARNING: Taking too long to find points for interpolation - please change the region to area where your points are. Continuing calculations... v.surf.rst complete. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Problem in function execution - GRASS 7
Hi, I am using WinGRASS 7 in a machine having Windows 7 OS. When i execute some functions such as v.category, v.rename and v.copy with menu, they keep on running without any message. To terminate, i use click on STOP. I tested this behaviour both on WinGRASS installed through OSGEO4W and nightly built binary for windows on 15th Decemeber 2012. However, both have this problem. Cheers, Brian. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Join points and contours
Hi, I have two 2D vector layers; spot heights and contours. Both represent elevation with an attribute field (say: level). Now i need to combine both vectors into one vector so that i can use combined vector to generate elevation raster. What is the best method to do that? Cheers, Brian. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Raster DEM generation
Thanks. I have not mistyped column name. In GRASS 7 that is what you get. On Mon, Jan 14, 2013 at 12:38 PM, Maris Nartiss wrote: > As it says. You have mistyped column name. Check column name and then try > again. > > Still in Your case I would try to interpolate with v.surf.rst first > and then use result for analysis. > > Maris. > > 2013/1/13 Brian Sanjeewa Rupasinghe : > > Hi Maris, > > > > Many thanks for the advice. Further, i am going to use r.surf.contour to > > generate a raster elevation DEM for viewshed analysis. The grass vector > > contour layer that i have needs to be converted to raster contour. For > this > > i tried using v.to.rast function. Elevation field of contour layer is; > > EL_Elevati which is a double precision data type. Following are the > > parameters i used: > > > > v.to.rast input = terrain@topo1K type = line output = rasterCont > use=attr > > attrcolumn = EL_Elevati labelcolumn = cat > > > > However, when i run this it gives following error: > > > > dbmi: protocol error > > ERROR: Column not found > > > > Any suggestion on this is appreciated. > > > > Cheers, Brian. > > > > > > On Sun, Jan 13, 2013 at 2:10 PM, Maris Nartiss > wrote: > >> > >> Hello, > >> there are two options: > >> 1st - join both datasets into single one (GRASS vectors can contain > >> different geometries in single vector dataset) and use v.surf.rst (it > >> will use both line and point data) > >> 2nd - convert lines to points with v.to.points and then join into > >> single point dataset. Interpolate with any tool. > >> > >> Good luck, > >> Maris. > >> > >> 2013/1/12 Brian Sanjeewa Rupasinghe : > >> > Hi all, > >> > > >> > have two shape files containing height information of a certain area. > >> > One > >> > file consists of height information in contours (isolines) and the > >> > othesr > >> > file consists of height information with spot heights. In order to > >> > create a > >> > complete DEM i need heights of both layers. I need to create a raster > >> > DEM > >> > with winGRASS 7 using these information. What would be the best > >> > technique to > >> > generate raster DEM with these two types of data? Do i need to > generate > >> > some > >> > spot heights from contours and then integerate with existing spot > >> > heights > >> > for DEM generation or any other method? > >> > > >> > Cheers, Brian > >> > > >> > ___ > >> > 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
Re: [GRASS-user] Raster DEM generation
Hi Maris, Many thanks for the advice. Further, i am going to use r.surf.contour to generate a raster elevation DEM for viewshed analysis. The grass vector contour layer that i have needs to be converted to raster contour. For this i tried using v.to.rast function. Elevation field of contour layer is; EL_Elevati which is a double precision data type. Following are the parameters i used: v.to.rast input = terrain@topo1K type = line output = rasterCont use=attr attrcolumn = EL_Elevati labelcolumn = cat However, when i run this it gives following error: dbmi: protocol error ERROR: Column not found Any suggestion on this is appreciated. Cheers, Brian. On Sun, Jan 13, 2013 at 2:10 PM, Maris Nartiss wrote: > Hello, > there are two options: > 1st - join both datasets into single one (GRASS vectors can contain > different geometries in single vector dataset) and use v.surf.rst (it > will use both line and point data) > 2nd - convert lines to points with v.to.points and then join into > single point dataset. Interpolate with any tool. > > Good luck, > Maris. > > 2013/1/12 Brian Sanjeewa Rupasinghe : > > Hi all, > > > > have two shape files containing height information of a certain area. One > > file consists of height information in contours (isolines) and the othesr > > file consists of height information with spot heights. In order to > create a > > complete DEM i need heights of both layers. I need to create a raster DEM > > with winGRASS 7 using these information. What would be the best > technique to > > generate raster DEM with these two types of data? Do i need to generate > some > > spot heights from contours and then integerate with existing spot heights > > for DEM generation or any other method? > > > > Cheers, Brian > > > > ___ > > 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] Raster DEM generation
Hi all, have two shape files containing height information of a certain area. One file consists of height information in contours (isolines) and the othesr file consists of height information with spot heights. In order to create a complete DEM i need heights of both layers. I need to create a raster DEM with winGRASS 7 using these information. What would be the best technique to generate raster DEM with these two types of data? Do i need to generate some spot heights from contours and then integerate with existing spot heights for DEM generation or any other method? Cheers, Brian ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Exporting GRASS vector to ShapeFile
Hi Maris, Thanks. Please see following the build info of the Grass Vector file that i need to convert to a shapefile. Cheers, Brian. (Sat Jan 12 14:25:02 2013) v.build map=MaskVector@user1 Building topology for vector map ... Registering primitives... 184 primitives registered 4948 vertices registered Building areas... 110 areas built 91 isles built Attaching islands... Attaching centroids... Number of nodes: 101 Number of primitives: 184 Number of points: 0 Number of lines: 0 Number of boundaries: 120 Number of centroids: 64 Number of areas: 110 Number of isles: 91 (Sat Jan 12 14:25:02 2013) Command finished (0 sec) On Sat, Jan 12, 2013 at 12:52 PM, Maris Nartiss wrote: > Please, post output of v.build > > 2013/1/12 Brian Sanjeewa Rupasinghe : > > Hi, > > > > I used v.category. But it too keeps on running failing to add categories. > > > > > > On Fri, Jan 11, 2013 at 11:28 PM, Maris Nartiss > wrote: > >> > >> Try to run on Your vector dataset: > >> v.category option=add type=centroid > >> Then export with v.out.ogr. > >> > >> Reason - v.out.ogr is working correctly and Your centroids have no > >> CATs and thus are not exported. That would explain "not correctly > >> created" - an empty shapefile. > >> > >> Hope it helps, > >> Maris. > >> > >> 2013/1/11 Brian Sanjeewa Rupasinghe : > >> > Hi, > >> > > >> > I am using GRASS 7svn in windows. Using v.out.ogr function in GRASS i > >> > converted a set of polygons (some of which have holes as well) into a > >> > shape > >> > file. > >> > However, when i do not tick the box in Selection pane: ‘Also export > >> > features > >> > without category (not labelled) otherwise only features with category > >> > are > >> > exported, > >> > the function keeps on executing and shapefile is not correctly > created. > >> > When > >> > i tick this, the shapefile is created but filling all holes as > separate > >> > polygons. Any comment on this behaviour? > >> > > >> > Cheers, Brian. > >> > > >> > ___ > >> > 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
Re: [GRASS-user] Exporting GRASS vector to ShapeFile
Hi, I used v.category. But it too keeps on running failing to add categories. On Fri, Jan 11, 2013 at 11:28 PM, Maris Nartiss wrote: > Try to run on Your vector dataset: > v.category option=add type=centroid > Then export with v.out.ogr. > > Reason - v.out.ogr is working correctly and Your centroids have no > CATs and thus are not exported. That would explain "not correctly > created" - an empty shapefile. > > Hope it helps, > Maris. > > 2013/1/11 Brian Sanjeewa Rupasinghe : > > Hi, > > > > I am using GRASS 7svn in windows. Using v.out.ogr function in GRASS i > > converted a set of polygons (some of which have holes as well) into a > shape > > file. > > However, when i do not tick the box in Selection pane: ‘Also export > features > > without category (not labelled) otherwise only features with category are > > exported, > > the function keeps on executing and shapefile is not correctly created. > When > > i tick this, the shapefile is created but filling all holes as separate > > polygons. Any comment on this behaviour? > > > > Cheers, Brian. > > > > ___ > > 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] Exporting GRASS vector to ShapeFile
Hi, I am using GRASS 7svn in windows. Using v.out.ogr function in GRASS i converted a set of polygons (some of which have holes as well) into a shape file. However, when i do not tick the box in Selection pane: ‘Also export features without category (not labelled) otherwise only features with category are exported, the function keeps on executing and shapefile is not correctly created. When i tick this, the shapefile is created but filling all holes as separate polygons. Any comment on this behaviour? Cheers, Brian. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Reclassification RASTER
Hi, I create each viewshed raster map for every observer location with flag (b) which denotes 0-non-visible and 1-visible values. Then i add all viewshed rasters to get integrated raster. So my integrated raster contains values >=1. In the raster, all values > 1 represents visible area. Therefore, i need to reclassify the image by setting all values > 1 to 1. If i use the following Python code to have output Name: MASK, is it ok? (out = integraed_raster) grass.mapcalc("MASK=if(($ReclassiName > 1),1,1)", ReclassiName = out) ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Use of g.mlist inside r.series
Many Thanks Martin. It worked fine. Cheers, Brian On Thu, Jan 10, 2013 at 6:55 PM, Martin Landa wrote: > Hi, > > 2013/1/10 Brian Sanjeewa Rupasinghe : > > grass.run_command ("r.series", input="`g.mlist pattern = 'viewshed*' sep > =, > > `", output = 'xx', method = 'sum') > > you cannot use bash-like constructions in Python script. Try > > maps = grass.mlist_strings(type = 'rast', pattern = 'viewshed*') > grass.run_command ("r.series", input = maps, output = 'xx', method = 'sum') > > Martin > > -- > Martin Landa * http://geo.fsv.cvut.cz/~landa > ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Use of g.mlist inside r.series
Hi, I am using following Python code to add all raster maps with wildcard name: viewshed*. grass.run_command ("r.series", input="`g.mlist pattern = 'viewshed*' sep =, `", output = 'xx', method = 'sum') But it does not work saying: WARNING: Illegal filename <`g.mlist pattern = 'viewshed*' sep =>. Character < > not allowed. ERROR: Raster map <`g.mlist pattern = 'viewshed*' sep => not found What is the possible error? I am using GRASS in Windows7 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Using a variable inside GRASS commands with Python code/ grass.mapcalc
Hi, Thanks again. Since i use GRASS 7.0, it will not work there. I do not think it can either be installed in GRASS for Windows through Add-on extensions. What about using grass.mapcalc to integrate each viewshed if there is a method of handling variables in it? If not how about using r.series with sum? Your advice is appreciated since i am pretty new to GRASS. On Tue, Jan 8, 2013 at 9:19 PM, Leonidas Liakos wrote: > Do you need cumulative viewshed analysis? > Maybe r.cva is for you: > http://www.ucl.ac.uk/~tcrnmar/GIS/r.cva.html > > -- > *Απο:* Brian Sanjeewa Rupasinghe > *Προς:* Anna Kratochvílová ; > grass-user@lists.osgeo.org > *Στάλθηκε:* 2:45 μ.μ. Τρίτη, 8 Ιανουαρίου 2013 > *Θέμα:* Re: [GRASS-user] Using a variable inside GRASS commands with > Python code/ grass.mapcalc > > Hi, > > Many thanks. It worked for me. Then i have a similar problem in > GRASS.mapcalc. What i need next is to add viewshed of each observer > location cumulatively while looping through the observer locations in order > to have final integrated viewshed analysis map. So i used the following > command in Python (some code ommitted) > > for i in: > grass.run_command('r.viewshed', input = rinput, output = 'viewshed', > coordinate = [x,y], obs_elev = oelv, tgt_elev = th, memory = 4098, flags = > 'b', overwrite = True, quiet = True) > > grass.mapcalc("viewshed_cum = viewshed + viewshed_cum", overwrite = > True, quiet = True) > > It seems that expressions like sum = c + sum does not work inside > grass.mapcalc.Is <http://grass.mapcalc.is/> there any alternative or > other way around to get this done? I am using GRASS 7 in Windows 7. > > Cheers, Brian > > On Tue, Jan 8, 2013 at 1:25 PM, Anna Kratochvílová > wrote: > > Hi Brian, > > On Tue, Jan 8, 2013 at 6:51 AM, Brian Sanjeewa Rupasinghe > wrote: > > Hi all > > > > I am generating viewshed of each observer location by iterating through a > > set of observer coordinate file through > > Python code in Windows 7 . Below is the python command used in grass for > > each iteration. According to this, each time > > loop iterates, output is overwritten with the same name. Now what i need > is > > to keep each output of all the observer locations. For this, i need to > > concatanate output name (i.e. 'viewshed' with loop index i which will be > > 'viewshed' + str(i) ). How is that possible within > > this grass.run_command? > > > > just set a variable within each loop and use it in the run_command: > for i in ... : > out = 'viewshed' + str(i) > grass.run_command('r.viewshed', input = rinput, output = out, ...) > > > Regards, > Anna > > > grass.run_command('r.viewshed', input = rinput, output = 'viewshed', > > coordinate = [x,y], obs_elev = oelv, tgt_elev = th, memory = 4098, flags > = > > 'b', overwrite = True, quiet = True) > > > > Cheers, Brian > > > > ___ > > 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 mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Using a variable inside GRASS commands with Python code/ grass.mapcalc
Hi, Many thanks. It worked for me. Then i have a similar problem in GRASS.mapcalc. What i need next is to add viewshed of each observer location cumulatively while looping through the observer locations in order to have final integrated viewshed analysis map. So i used the following command in Python (some code ommitted) for i in: grass.run_command('r.viewshed', input = rinput, output = 'viewshed', coordinate = [x,y], obs_elev = oelv, tgt_elev = th, memory = 4098, flags = 'b', overwrite = True, quiet = True) grass.mapcalc("viewshed_cum = viewshed + viewshed_cum", overwrite = True, quiet = True) It seems that expressions like sum = c + sum does not work inside grass.mapcalc.Is there any alternative or other way around to get this done? I am using GRASS 7 in Windows 7. Cheers, Brian On Tue, Jan 8, 2013 at 1:25 PM, Anna Kratochvílová wrote: > Hi Brian, > > On Tue, Jan 8, 2013 at 6:51 AM, Brian Sanjeewa Rupasinghe > wrote: > > Hi all > > > > I am generating viewshed of each observer location by iterating through a > > set of observer coordinate file through > > Python code in Windows 7 . Below is the python command used in grass for > > each iteration. According to this, each time > > loop iterates, output is overwritten with the same name. Now what i need > is > > to keep each output of all the observer locations. For this, i need to > > concatanate output name (i.e. 'viewshed' with loop index i which will be > > 'viewshed' + str(i) ). How is that possible within > > this grass.run_command? > > > > just set a variable within each loop and use it in the run_command: > for i in ... : > out = 'viewshed' + str(i) > grass.run_command('r.viewshed', input = rinput, output = out, ...) > > > Regards, > Anna > > > grass.run_command('r.viewshed', input = rinput, output = 'viewshed', > > coordinate = [x,y], obs_elev = oelv, tgt_elev = th, memory = 4098, flags > = > > 'b', overwrite = True, quiet = True) > > > > Cheers, Brian > > > > ___ > > 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] Using a variable inside GRASS commands with Python code
Hi all I am generating viewshed of each observer location by iterating through a set of observer coordinate file through Python code in Windows 7 . Below is the python command used in grass for each iteration. According to this, each time loop iterates, output is overwritten with the same name. Now what i need is to keep each output of all the observer locations. For this, i need to concatanate output name (i.e. 'viewshed' with loop index i which will be 'viewshed' + str(i) ). How is that possible within this grass.run_command? grass.run_command('r.viewshed', input = rinput, output = 'viewshed', coordinate = [x,y], obs_elev = oelv, tgt_elev = th, memory = 4098, flags = 'b', overwrite = True, quiet = True) Cheers, Brian ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Installing r.viewshed add-on in GRASS 6.4.3 in Windows 7
Hi, Markus I could install r.viewshed without any trouble with GRASS 6.4.3svn from SVN nightly built version. It does not work with stable winGRASS 6 probably a bug!. Cheers, Brian. On Wed, Dec 19, 2012 at 10:24 PM, Markus Neteler wrote: > (please keep the list in the loop) > > On Wed, Dec 19, 2012 at 9:44 PM, Markus Neteler wrote: > > On Wed, Dec 19, 2012 at 9:42 PM, Brian Sanjeewa Rupasinghe > > wrote: > >> Thanks. Yes with Graphical user interface. In that I used running with > >> Coomand output, but the same error. > >> Initially i got and error saying following: > >> > >> ERROR: Unable to create 'C:/Program Files/GRASS GIS 6.4.3RC1\man\man1': > >> [Error 5] Access is denied: 'C:/Program Files/GRASS GIS 6.4.3RC1\\man' > >> (Wed Dec 19 20:41:32 2012) Command finished (0 sec) > >> > > I am not on Windows myself, hence no idea. > > Commonly we suggest to not install GRASS in paths with white space. > > Regarding your other offlist question: no compiler needed since the > installer will just fetch the ZIP'ed executable from our server and > unpacks it. > > Markus > ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Installing r.viewshed add-on in GRASS 6.4.3 in Windows 7
Hi all, I am a newbie to GRASS and i recently installed GRASS 6.4.3 for Windows 7. When I try to add r.viewshed add-on, it gives an error as follows: g.extension.py extension=r.viewshed svnurl= http://svn.osgeo.org/grass/grass-addons/grass6 Downloading precompiled GRASS Addons ... ERROR: GRASS Addons not found (Wed Dec 19 20:17:56 2012) Command finished (1 sec) Could somebody help to solve this? Cheers, Brian. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user