[GRASS-user] Problem in v.surf.rst

2013-01-17 Thread Brian Sanjeewa Rupasinghe
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

2013-01-16 Thread Brian Sanjeewa Rupasinghe
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

2013-01-16 Thread Brian Sanjeewa Rupasinghe
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

2013-01-16 Thread Brian Sanjeewa Rupasinghe
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

2013-01-15 Thread Brian Sanjeewa Rupasinghe
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

2013-01-14 Thread Brian Sanjeewa Rupasinghe
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

2013-01-13 Thread Brian Sanjeewa Rupasinghe
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

2013-01-13 Thread 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


[GRASS-user] Raster DEM generation

2013-01-12 Thread 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


Re: [GRASS-user] Exporting GRASS vector to ShapeFile

2013-01-12 Thread Brian Sanjeewa Rupasinghe
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

2013-01-11 Thread 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


[GRASS-user] Exporting GRASS vector to ShapeFile

2013-01-11 Thread 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] Reclassification RASTER

2013-01-10 Thread Brian Sanjeewa Rupasinghe
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

2013-01-10 Thread Brian Sanjeewa Rupasinghe
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

2013-01-10 Thread Brian Sanjeewa Rupasinghe
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

2013-01-08 Thread Brian Sanjeewa Rupasinghe
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

2013-01-08 Thread Brian Sanjeewa Rupasinghe
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

2013-01-07 Thread Brian Sanjeewa Rupasinghe
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

2012-12-20 Thread Brian Sanjeewa Rupasinghe
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

2012-12-19 Thread Brian Sanjeewa Rupasinghe
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