[GRASS-user] multiprocessing in python

[Leonardo Hardtke]

Dear all,
I am working on a module to extract the phenological parameters (like
timesat) from a time series implemented in python/cython and making use of
gscript and other grass stuff.
It works great on a 256x256 and as the plan is applying it over Australia
at 250m over 17 years, I need to split the process in small tiles. The idea
is to run this processes in parallel and I am having issues implementing
it.

This would be the first part of the process that runs on each tile:

def tile_process(tile_index):
'''
Function for every worker:
Applies any function to the sub_region corresponding to the tile_index.
'''
global Rows
global Cols
global RowBlockSize
global ColBlockSize
global full_region
global dates
global years
global indices
global data_serie
global yr_limits_extra
global yr_limits
global dbif

sub_name='block'
TileRow, TileCol, sr =
sub_region(tile_index,full_region,RowBlockSize,ColBlockSize)
# # Define a temporary region based on the parameters caluculated with
the
start_row = TileRow * RowBlockSize
start_col = TileCol * ColBlockSize
n_rows = sr['rows']
n_cols = sr['cols']

strds = tgis.SpaceTimeRasterDataset(data_serie)
strds.select(dbif=dbif)
maps = strds.get_registered_maps_as_objects(dbif=dbif)

# Numer of time steps
steps = len(maps)
# Make an empty array
#print(steps)
EVI = np.empty([steps,n_rows,n_cols])
# fill the array
for step, map in enumerate(maps):
 map.select(dbif=dbif)
 image_name = map.get_name()+'@'+data_serie.split('@')[1]
 #print("reading: {}".format(image_name))
 EVI[step,:] =
raster2numpy_sub(image_name,start_row,n_rows,start_col,n_cols)
mean = EVI.mean()
print(mean)





and this is how I start the multiprocess pool.

pool.map(tile_process, xrange(RowBlockN*ColBlockN))
pool.close()
pool.join()

and it gives me:

AssertionError: can only test a child process


of course if I do: tile_process(0) or tile_process(1) etc ,the right result
comes out.

Does any of you have experience with this? Any suggestion would be welcome!
Sorry for the messy code. Is still in early stage.


-- 
Dr. Leonardo A. Hardtke
C3 UTS, Scientific Officer
CB04.06.315.06
Email:leonardoandres.hard...@uts.edu.au or leohard...@gmail.com
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] How to get a list of maps from a group in python script

[Leonardo Hardtke]

Oh!
Thanks "g" flag combined with .split() worked

Leo

2016-04-06 11:59 GMT+10:00 Anna Petrášová <kratocha...@gmail.com>:

> On Tue, Apr 5, 2016 at 9:15 PM, Leonardo Hardtke <leohard...@gmail.com>
> wrote:
> > Thanks Laurent for your help,
> > it worked to split the string to a listbut now I have to remove the 1
> > line(group  references the following raster
> > maps), line 2 ('-'), the last line ('-') and the
> > "<>" around the names...
> > I know how to do that but is this the best way to deal with outputs
> of
> > i.group -l in python?
> >
> >
> >
>
>
> You should use i.group -g (Print in shell script style)
>
>
> Anna
>
>
> > 2016-04-06 11:02 GMT+10:00 Laurent C. <lrn...@gmail.com>:
> >>
> >> Hello,
> >>
> >> Try the .splitlines() method on the string.
> >>
> >> Cheers,
> >> Laurent
> >>
> >> El abr. 5, 2016 19:52, "Leonardo Hardtke" <leohard...@gmail.com>
> escribió:
> >>>
> >>> Dear list members,
> >>> I know this is more a python question than a GRASS question, bot for
> sure
> >>> some of you knows the ansewer and is able to help me
> >>>
> >>> I need to get a python list with the maps from a group, like
> >>> ["map1","map2","map"]
> >>>
> >>> But by using gscript.read_command i get string, not a list.
> >>>
> >>> gscript.read_command("i.group", group=groupname, flags="l" )
> >>> 'group  references the following raster
> >>> maps\n-\n<LE72270902002150EDC00_toar.1@landsat>
> >>> \n<LE72270902002150EDC00_toar.2@landsat>
> >>> \n<LE72270902002150EDC00_toar.3@landsat>
> >>> \n<LE72270902002150EDC00_toar.4@landsat>
> >>> \n<LE72270902002150EDC00_toar.5@landsat>
> >>> \n<LE72270902002150EDC00_toar.7@landsat>\n-\n'
> >>>
> >>> Your help is appreciated!
> >>> Leo
> >>> --
> >>> Dr. Leonardo A. Hardtke
> >>> Laboratorio de Teledetección y S.I.G.
> >>> Centro Nacional Patagónico (CONICET)
> >>> Bvd. Brown 2825, 9120
> >>> Puerto Madryn, Chubut, Argentina
> >>>
> >>> ___
> >>> grass-user mailing list
> >>> grass-user@lists.osgeo.org
> >>> http://lists.osgeo.org/mailman/listinfo/grass-user
> >
> >
> >
> >
> > --
> > Dr. Leonardo A. Hardtke
> > Laboratorio de Teledetección y S.I.G.
> > Centro Nacional Patagónico (CONICET)
> > Bvd. Brown 2825, 9120
> > Puerto Madryn, Chubut, Argentina
> >
> > ___
> > grass-user mailing list
> > grass-user@lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/grass-user
>



-- 
Dr. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120
Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] How to get a list of maps from a group in python script

[Leonardo Hardtke]

Thanks Laurent for your help,
it worked to split the string to a listbut now I have to remove the 1
line(group  references the following raster
maps), line 2 ('-'), the last line ('-') and the
"<>" around the names...
I know how to do that but is this the best way to deal with outputs of
i.group -l in python?



2016-04-06 11:02 GMT+10:00 Laurent C. <lrn...@gmail.com>:

> Hello,
>
> Try the .splitlines() method on the string.
>
> Cheers,
> Laurent
> El abr. 5, 2016 19:52, "Leonardo Hardtke" <leohard...@gmail.com> escribió:
>
>> Dear list members,
>> I know this is more a python question than a GRASS question, bot for sure
>> some of you knows the ansewer and is able to help me
>>
>> I need to get a python list with the maps from a group, like
>> ["map1","map2","map"]
>>
>> But by using gscript.read_command i get string, not a list.
>>
>> gscript.read_command("i.group", group=groupname, flags="l" )
>> 'group  references the following raster
>> maps\n-\n<LE72270902002150EDC00_toar.1@landsat>
>>  \n<LE72270902002150EDC00_toar.2@landsat>
>>  \n<LE72270902002150EDC00_toar.3@landsat>
>>  \n<LE72270902002150EDC00_toar.4@landsat>
>>  \n<LE72270902002150EDC00_toar.5@landsat>
>>  \n<LE72270902002150EDC00_toar.7@landsat>\n-\n'
>>
>> Your help is appreciated!
>> Leo
>> --
>> Dr. Leonardo A. Hardtke
>> Laboratorio de Teledetección y S.I.G.
>> Centro Nacional Patagónico (CONICET)
>> Bvd. Brown 2825, 9120
>> Puerto Madryn, Chubut, Argentina
>>
>> ___
>> grass-user mailing list
>> grass-user@lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/grass-user
>>
>


-- 
Dr. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120
Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] How to get a list of maps from a group in python script

[Leonardo Hardtke]

Dear list members,
I know this is more a python question than a GRASS question, bot for sure
some of you knows the ansewer and is able to help me

I need to get a python list with the maps from a group, like
["map1","map2","map"]

But by using gscript.read_command i get string, not a list.

gscript.read_command("i.group", group=groupname, flags="l" )
'group  references the following raster
maps\n-\n
 \n
 \n
 \n
 \n
 \n\n-\n'

Your help is appreciated!
Leo
-- 
Dr. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120
Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] problem with python scripting

[Leonardo Hardtke]

Thanks Paulo and Sören for helping!
overwrite=True made the trick! Is this documented somewhere and I just
missed it??? Or should I fill a ticket?

Leo



2016-04-05 17:06 GMT+10:00 Sören Gebbert <soerengebb...@googlemail.com>:

> Please try gscript.mapcalc(expr... , overwrite=True)
> Am 05.04.2016 09:02 schrieb "Leonardo Hardtke" <leohard...@gmail.com>:
>
>> Hi Paulo and thanks for the answer... but It gives me a very similar
>> error if I write the full flag (The only difference is TypeError... str
>> if I use 'o' and int if I use 'overwrite' )
>>
>>  gscript.mapcalc('endmember_rstr.1 = null()', flags="overwrite")
>>
>> ERROR: output map  exists. To overwrite, use the
>>--overwrite flag
>> Traceback (most recent call last):
>>   File "", line 1, in 
>>   File "/usr/lib/grass70/etc/python/grass/script/raster.py", line 103, in
>> mapcalc
>> quiet=quiet, verbose=verbose, overwrite=overwrite)
>>   File "/usr/lib/grass70/etc/python/grass/script/core.py", line 516, in
>> write_command
>> return handle_errors(returncode, returncode, args, kwargs)
>>   File "/usr/lib/grass70/etc/python/grass/script/core.py", line 312, in
>> handle_errors
>> returncode=returncode)
>>   File "/usr/lib/grass70/etc/python/grass/exceptions/__init__.py", line
>> 68, in __init__
>> msg = _("Module run %s %s ended with error") % (module, code)
>> TypeError: 'int' object is not callable
>>
>> Thanks for your help!
>>
>> 2016-04-05 16:56 GMT+10:00 Paulo van Breugel <p.vanbreu...@gmail.com>:
>>
>>>
>>>
>>> On 05-04-16 08:53, Leonardo Hardtke wrote:
>>>
>>> Dear users,
>>> I am writing a script and I'm having problems with the flags for mapcalc
>>> gscript.mapcalc('endmember_rstr.1 = null()', flags="o")
>>>
>>>
>>> If  you are using GRASS 7.+, the overwrite flag should be written in
>>> full, i.e.,
>>>
>>> gscript.mapcalc('endmember_rstr.1 = null()', flags="overwrite")
>>>
>>>
>>>
>>> ERROR: output map  exists. To overwrite, use the
>>>--overwrite flag
>>> Traceback (most recent call last):
>>>   File "", line 1, in 
>>>   File "/usr/lib/grass70/etc/python/grass/script/raster.py", line 103,
>>> in mapcalc
>>> quiet=quiet, verbose=verbose, overwrite=overwrite)
>>>   File "/usr/lib/grass70/etc/python/grass/script/core.py", line 516, in
>>> write_command
>>> return handle_errors(returncode, returncode, args, kwargs)
>>>   File "/usr/lib/grass70/etc/python/grass/script/core.py", line 312, in
>>> handle_errors
>>> returncode=returncode)
>>>   File "/usr/lib/grass70/etc/python/grass/exceptions/__init__.py", line
>>> 68, in __init__
>>> msg = _("Module run %s %s ended with error") % (module, code)
>>> TypeError: 'str' object is not callable
>>>
>>> with other modules flags are working as I expect, like
>>>
>>> gscript.run_command("g.region", flags="p")
>>> projection: 1 (UTM)
>>> zone:   -20
>>> datum:  wgs84
>>> ellipsoid:  wgs84
>>> north:  5326556.02809
>>> south:  5228355.32592
>>> west:   290546.512535
>>> east:   444008.446334
>>> nsres:  30.00326984
>>> ewres:  30.0081998
>>> rows:   3273
>>> cols:   5114
>>> cells:  16738122
>>>
>>> Am I missing something?
>>>
>>> Thanks!
>>> --
>>> Dr. Leonardo A. Hardtke
>>> Laboratorio de Teledetección y S.I.G.
>>> Centro Nacional Patagónico (CONICET)
>>> Bvd. Brown 2825, 9120
>>> Puerto Madryn, Chubut, Argentina
>>>
>>>
>>> ___
>>> grass-user mailing 
>>> listgrass-user@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user
>>>
>>>
>>>
>>
>>
>> --
>> Dr. Leonardo A. Hardtke
>> Laboratorio de Teledetección y S.I.G.
>> Centro Nacional Patagónico (CONICET)
>> Bvd. Brown 2825, 9120
>> Puerto Madryn, Chubut, Argentina
>>
>> ___
>> grass-user mailing list
>> grass-user@lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/grass-user
>>
>


-- 
Dr. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120
Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] problem with python scripting

[Leonardo Hardtke]

Hi Paulo and thanks for the answer... but It gives me a very similar error
if I write the full flag (The only difference is TypeError... str if I
use 'o' and int if I use 'overwrite' )

 gscript.mapcalc('endmember_rstr.1 = null()', flags="overwrite")

ERROR: output map  exists. To overwrite, use the
   --overwrite flag
Traceback (most recent call last):
  File "", line 1, in 
  File "/usr/lib/grass70/etc/python/grass/script/raster.py", line 103, in
mapcalc
quiet=quiet, verbose=verbose, overwrite=overwrite)
  File "/usr/lib/grass70/etc/python/grass/script/core.py", line 516, in
write_command
return handle_errors(returncode, returncode, args, kwargs)
  File "/usr/lib/grass70/etc/python/grass/script/core.py", line 312, in
handle_errors
returncode=returncode)
  File "/usr/lib/grass70/etc/python/grass/exceptions/__init__.py", line 68,
in __init__
msg = _("Module run %s %s ended with error") % (module, code)
TypeError: 'int' object is not callable

Thanks for your help!

2016-04-05 16:56 GMT+10:00 Paulo van Breugel <p.vanbreu...@gmail.com>:

>
>
> On 05-04-16 08:53, Leonardo Hardtke wrote:
>
> Dear users,
> I am writing a script and I'm having problems with the flags for mapcalc
> gscript.mapcalc('endmember_rstr.1 = null()', flags="o")
>
>
> If  you are using GRASS 7.+, the overwrite flag should be written in full,
> i.e.,
>
> gscript.mapcalc('endmember_rstr.1 = null()', flags="overwrite")
>
>
>
> ERROR: output map  exists. To overwrite, use the
>--overwrite flag
> Traceback (most recent call last):
>   File "", line 1, in 
>   File "/usr/lib/grass70/etc/python/grass/script/raster.py", line 103, in
> mapcalc
> quiet=quiet, verbose=verbose, overwrite=overwrite)
>   File "/usr/lib/grass70/etc/python/grass/script/core.py", line 516, in
> write_command
> return handle_errors(returncode, returncode, args, kwargs)
>   File "/usr/lib/grass70/etc/python/grass/script/core.py", line 312, in
> handle_errors
> returncode=returncode)
>   File "/usr/lib/grass70/etc/python/grass/exceptions/__init__.py", line
> 68, in __init__
> msg = _("Module run %s %s ended with error") % (module, code)
> TypeError: 'str' object is not callable
>
> with other modules flags are working as I expect, like
>
> gscript.run_command("g.region", flags="p")
> projection: 1 (UTM)
> zone:   -20
> datum:  wgs84
> ellipsoid:  wgs84
> north:  5326556.02809
> south:  5228355.32592
> west:   290546.512535
> east:   444008.446334
> nsres:  30.00326984
> ewres:  30.0081998
> rows:   3273
> cols:   5114
> cells:  16738122
>
> Am I missing something?
>
> Thanks!
> --
> Dr. Leonardo A. Hardtke
> Laboratorio de Teledetección y S.I.G.
> Centro Nacional Patagónico (CONICET)
> Bvd. Brown 2825, 9120
> Puerto Madryn, Chubut, Argentina
>
>
> ___
> grass-user mailing 
> listgrass-user@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user
>
>
>


-- 
Dr. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120
Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] problem with python scripting

[Leonardo Hardtke]

Dear users,
I am writing a script and I'm having problems with the flags for mapcalc
gscript.mapcalc('endmember_rstr.1 = null()', flags="o")

ERROR: output map  exists. To overwrite, use the
   --overwrite flag
Traceback (most recent call last):
  File "", line 1, in 
  File "/usr/lib/grass70/etc/python/grass/script/raster.py", line 103, in
mapcalc
quiet=quiet, verbose=verbose, overwrite=overwrite)
  File "/usr/lib/grass70/etc/python/grass/script/core.py", line 516, in
write_command
return handle_errors(returncode, returncode, args, kwargs)
  File "/usr/lib/grass70/etc/python/grass/script/core.py", line 312, in
handle_errors
returncode=returncode)
  File "/usr/lib/grass70/etc/python/grass/exceptions/__init__.py", line 68,
in __init__
msg = _("Module run %s %s ended with error") % (module, code)
TypeError: 'str' object is not callable

with other modules flags are working as I expect, like

gscript.run_command("g.region", flags="p")
projection: 1 (UTM)
zone:   -20
datum:  wgs84
ellipsoid:  wgs84
north:  5326556.02809
south:  5228355.32592
west:   290546.512535
east:   444008.446334
nsres:  30.00326984
ewres:  30.0081998
rows:   3273
cols:   5114
cells:  16738122

Am I missing something?

Thanks!
-- 
Dr. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120
Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Is Watershed image segmentation program in GRASS GIS?

[Leonardo Hardtke]

Hi there...

I used r.grow to segment dNBR images to map wildfires. See:


https://www.researchgate.net/publication/269692754_Semi-automated_mapping_of_burned_areas_in_semi-arid_ecosystems_using_MODIS_time-series_imagery

Hope it helps...


2015-11-11 3:47 GMT-03:00 Uttam Kumar :

>
>
> Hi All,
>
> Is watershed image segmentation program or a similar module available in
> GRASS GIS?
>
> I understand that we have r.watershed, i.segment and r.grow in GRASS GIS.
> But I do not know which one best suits my problem.
>
> The problem is:
>
> I have two layers of raster data -
> 1.) a gradient image (obtained by applying a 3 x 3 sobel filter on a water
> index using r.mapcalc), and
> 2.) a seed image of 2 classes labelled as "1" and "2" and remaining pixels
> as "0".
>
> I want to use the seed image and the gradient image (which shows direction
> of change in intensity) information to perform segmentation using the
> watershed algorithm and label all "0" pixels as either 1 or 2.
>
> r.watershed does not use watershed segmentation algorithm. i.segment is
> something close but I am not sure.
>
> Please suggest.
>
> Any comment is appreciated.
>
> Uttam.
>
>
>
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>



-- 
Dr. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120
Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] problem with t.rast.mapcalc (grass71 daily ubuntu 14.04)

[Leonardo Hardtke]

Thaks
If it is a bug, I should fill a bug report so that some developer can fix
it!

Cheers

2015-09-09 15:17 GMT-03:00 Veronica Andreo <veroand...@gmail.com>:

> Hi Leo,
>
> Yes, I had the same problem some time ago and the solution was to change
> the name, as you mentioned...
>
> Maybe a bug? Dunno :)
>
> Cheers,
> Vero
>
> El sep 9, 2015 3:05 PM, "Leonardo Hardtke" <leohard...@gmail.com>
> escribió:
>
>> Dear GRASS users (and developers?!),
>>
>> I'm facing a problem with an expression in the t.rast.mapcalc module...
>>
>> t.rast.mapcalc inputs=QA_mask,ndvi,ndvi_smooth_spacetime
>> expression="if(QA_mask==0,ndvi,ndvi_smooth_spacetime)" output=ndvi_fixed
>> basename=ndvi_fixed
>> It throws the following error:
>> Starting temporal sampling...
>> Starting mapcalc computation...
>> Invalid map
>> 
>> Parse error
>> ERROR: parse error
>> ERROR: Error while mapcalc computation
>>
>> The QA_mask maps are:
>> QA_mask_01|PERMANENT|2009-01-01 00:00:00|2009-01-17 00:00:00
>> QA_mask_02|PERMANENT|2009-01-17 00:00:00|2009-02-02 00:00:00
>> ..
>>
>> The ndvi_smooth_spacetime maps are:
>> ndvi_smooth_spacetime_0|PERMANENT|2009-01-17 00:00:00|2009-02-02 00:00:00
>> ndvi_smooth_spacetime_1|PERMANENT|2009-02-02 00:00:00|2009-02-18 00:00:00
>> ...
>>
>> The ndvi maps are
>> MOD13Q1.A2009001.h12v12.005.2009019145228_01_4326|PERMANENT|2009-01-01
>> 00:00:00|2009-01-17 00:00:00
>> MOD13Q1.A2009017.h12v12.005.2009035110230_01_4326|PERMANENT|2009-01-17
>> 00:00:00|2009-02-02 00:00:00
>> M
>>
>> As you see, there is no such a map like..
>> MOD13Q1.A2009017.h12v12.005.2009035110230_01_4326_smooth_spacetime_0
>>
>> I guess the module is mixing ndvi (a temporal dataset) with
>> _smooth_spacetime... the second part of the name
>>
>> I know it work ok by renaming the ndvi_smooth_spacetime by something
>> else... but... is this a bug or shoul I use more '`" to protect names !?
>>
>> Thanks!
>> --
>> Dr. Leonardo A. Hardtke
>> Laboratorio de Teledetección y S.I.G.
>> Centro Nacional Patagónico (CONICET)
>> Bvd. Brown 2825, 9120
>> Puerto Madryn, Chubut, Argentina
>>
>> ___
>> grass-user mailing list
>> grass-user@lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/grass-user
>>
>


-- 
Dr. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120
Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] problem with t.rast.mapcalc (grass71 daily ubuntu 14.04)

[Leonardo Hardtke]

Dear GRASS users (and developers?!),

I'm facing a problem with an expression in the t.rast.mapcalc module...

t.rast.mapcalc inputs=QA_mask,ndvi,ndvi_smooth_spacetime
expression="if(QA_mask==0,ndvi,ndvi_smooth_spacetime)" output=ndvi_fixed
basename=ndvi_fixed
It throws the following error:
Starting temporal sampling...
Starting mapcalc computation...
Invalid map

Parse error
ERROR: parse error
ERROR: Error while mapcalc computation

The QA_mask maps are:
QA_mask_01|PERMANENT|2009-01-01 00:00:00|2009-01-17 00:00:00
QA_mask_02|PERMANENT|2009-01-17 00:00:00|2009-02-02 00:00:00
..

The ndvi_smooth_spacetime maps are:
ndvi_smooth_spacetime_0|PERMANENT|2009-01-17 00:00:00|2009-02-02 00:00:00
ndvi_smooth_spacetime_1|PERMANENT|2009-02-02 00:00:00|2009-02-18 00:00:00
...

The ndvi maps are
MOD13Q1.A2009001.h12v12.005.2009019145228_01_4326|PERMANENT|2009-01-01
00:00:00|2009-01-17 00:00:00
MOD13Q1.A2009017.h12v12.005.2009035110230_01_4326|PERMANENT|2009-01-17
00:00:00|2009-02-02 00:00:00
M

As you see, there is no such a map like..
MOD13Q1.A2009017.h12v12.005.2009035110230_01_4326_smooth_spacetime_0

I guess the module is mixing ndvi (a temporal dataset) with
_smooth_spacetime... the second part of the name

I know it work ok by renaming the ndvi_smooth_spacetime by something
else... but... is this a bug or shoul I use more '`" to protect names !?

Thanks!
-- 
Dr. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120
Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] t.rast.mapcalc question

[Leonardo Hardtke]

Oh! Thanks!
I was using GRASS 7.0.1, which has no t.rast.algebra module.
I didn't noticed there where a t.rast.mapcalc and a t.rast.algebra in newer
versions...
I installed grass71 and as you said,

t.rast.algebra expression=ndvi_diff = ndvi[0] - ndvi[1] basename=ndvi_diff
worked perfect! Much easy than the old mapcalc and for loop version!
 I'm really enjoying the grass 7 temporal framework!
 Thanks for the fast (and accurate) answer!

2015-08-27 10:04 GMT-03:00 Veronica Andreo veroand...@gmail.com:

 Hola Leo!

 You should use t.rast.algebra (
 https://grass.osgeo.org/grass71/manuals/t.rast.algebra.html) from GRASS
 GIS 7.1 for that kind of temporal algebra operation. The command should be
 like:

 t.rast.algebra expression=ndvi_a = ndvi[1] - ndvi[0] basename=ndvi_a

 For more examples on temporal modules, you can see:
 https://grasswiki.osgeo.org/wiki/Temporal_data_processing

 HTH,
 Vero

 2015-08-27 9:52 GMT-03:00 Leonardo Hardtke leohard...@gmail.com:

 Dear GRASS users:
 I'm experienced grass user...  I use GRASS since I started my PhD
 research in landscape ecology back in 2009...
 I'm trying to calculate the difference of NDVI between consecutive dates.
 I used to do it like:

 i=0
 for map in $(g.list rast pattern=MOD13Q1*_01_*); do
 NDVI[i]=$map
 echo Adding:  $map
 let i=i+1
 done

 for id in $(seq 0 $(echo ${#NDVI[@]})); do
 echo $id
  r.mapcalc ndvi_dif2_$id = ${NDVI[$id+1]} - ${NDVI[id]}
 done

 But i thing there must be an easier way using the grass 7 temporal
 framework. So i registred my dataset and tried the following

 t.rast.list input=ndvi@ts
 MOD13Q1.A2009001.h12v12.005.2009019145228_01_4326|ts|2009-01-01
 00:00:00|2009-01-17 00:00:00
 MOD13Q1.A2009017.h12v12.005.2009035110230_01_4326|ts|2009-01-17
 00:00:00|2009-02-02 00:00:00

 
 MOD13Q1.A2010337.h12v12.005.2010355100650_01_4326|ts|2010-12-03
 00:00:00|2010-12-19 00:00:00
 MOD13Q1.A2010353.h12v12.005.2011005202724_01_4326|ts|2010-12-19
 00:00:00|2011-01-04 00:00:00


 t.rast.mapcalc input=ndvi expression=ndvi[1] - ndvi[0] output=ndvi_a
 basename=ndvi_a

 But the output is not good (see attached file: left is OK, rigth is wrong)
 Am I doing something wrong?
 Thanks!
 Dr. Leonardo A. Hardtke
 Laboratorio de Teledetección y S.I.G.
 Centro Nacional Patagónico (CONICET)
 Bvd. Brown 2825, 9120
 Puerto Madryn, Chubut, Argentina

 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user





-- 
Dr. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120
Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] read raster by block...

[Leonardo Hardtke]

Hi list,
I made a gdal-phyton script that reads big images (a 230 images
timeseries 4000x4000 px) by blocks, something like
data = band.ReadAsArray(j, i, xBlSize, yBlSize)
someone knows how to do that in a Grass-phyton script (or has a better idea! )?
Thaks
Leonardo Hardtke

-- 
Lic. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120
Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] wx.class

[Leonardo Hardtke]

Hi,
I tried to use it, but

   GRASS 6.4.0RC6 (Spearfish):~  ./wx.class.py
   ./wx.class.py: line 19: !
   @package wx.class.py

   @brief Interface implementation of i.class module without Xterm.
   Many thanks to Markus Neteler for his help which allowed me to
   know about i.class and its usage.
   Classes:
   - IClass
   - BufferedWindow2
   - IClassApp


   (C) 2006-2010 by the GRASS Development Team
   This program is free software under the GNU General Public
   License (=v2). Read the file COPYING that comes with GRASS
   for details.

   @author Mohammed Rashad K.M rashadkm [at] gmail [dot] com
   : command not found
   import: unable to grab mouse `': Resource temporarily unavailable @
   xwindow.c/XSelectWindow/8993.
   import: unable to grab mouse `': Resource temporarily unavailable @
   xwindow.c/XSelectWindow/8993.
   import: unable to grab mouse `': Resource temporarily unavailable @
   xwindow.c/XSelectWindow/8993.
   import: unable to grab mouse `': Resource temporarily unavailable @
   xwindow.c/XSelectWindow/8993.
   ./wx.class.py: line 28: syntax error near unexpected token `('
   ./wx.class.py: line 28: `gbase = os.getenv(GISBASE) '

any idea?
Leonardo A. Hardtke

On 09/08/2010 06:34 PM, Mohammed Rashad wrote:


Hi all,

 I had just created a wx version of i.class module which generates 
spectral signatures for an image by allowing the user to outline 
regions of interest. The resulting signature file can be used as input 
for i.maxlik or as a seed signature file for i.cluster.


Actual i.class module requires Xterm and work only in systems having 
Xserver installed. I created python version of i.class which works on 
all platform. Now windows users can use i.class module.


Currently development is in beta stage. wx.class is tested and works 
only for NC dataset available from grassbook.org 
http://grassbook.org site.
It uses wxpython paintDC to draw histograms which must be changed to 
use python-matplotlib for drawing smooth histograms as the output 
needed to be.


I dont want to add a new dependency to grass as other developers but I 
am forced to do :(


current wx.class doesn't need any other dependencies. The needs some 
cleanup but as I am busy with some other works related to GRASS I 
don't have much time to spend on cleaning up the code. But within a 
month or to you can have the fully functional wx.class.


Now I am hosting to GRASS-Addons repo because the code needs 
testing...(although needs some cleanup)


A README is included with source which shows how to use the module. 
But wx.class is made to be user-friendly.



Any comments or suggestions are really welcomed :)


--
Rashad


___
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] Cluster whit area?

[Leonardo Hardtke]

Hi, you are missing something!

library(spgrass6) #start grassR interface
Vine1- readVECT6(vineyardheroiq...@etienne) #read vector as SpatialDataFrame
Vine1.df=as.data.frame(Vine1) # Transform to dataframe

library(clustTool) #start GUI

On 07/29/2010 05:17 AM, Etienne DELAY (deust 2010) wrote:

Hi Leonardo
this morning I try your trick :-p
But can you explain in more detail the handling ..

I run Grass
Run R in GARSS
library(spgrass6)
Vine1.df=as.data.frame(Vine1)
He make an erreur
Erreur dans as.data.frame(Vine1)

So if you can get me more help... :-)
Thank tou


Le 28/07/10 20:21, Leonardo Hardtke a écrit :
   

Hi!
I never used clustTool before.
  I made a try and i think I found the trick :-).
  You have to transform your (SpatialPointDataFrame) data to a dataframe...
 Vine1.df=as.data.frame(Vine1)
  and then you can Set acctive data in the GUI!
hope it helps!

Leonardo
On 07/27/2010 09:06 PM, Etienne DELAY (deust 2010) wrote:
 

Le 27/07/10 20:35, Hamish a écrit :

   

Leonardo wrote:

 

you can run R in a grass terminal, and load data via
spgrass6 package,
readRAST6 or readVECT6.

   

see  http://grass.osgeo.org/wiki/R_statistics


 

Thank you Leonardo, but I read it, without understanding much more about
the use of clustTool... (it's ok whit readVECT ;-) )
I dont understand...

library(spgrass6)

Vine1- readVECT6(vineyardheroiq...@etienne)
library(clustTool)

it's run a graphical windows...
an I can't choose my data ;-(

so I need an other help

thank you


   

___
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] Cluster whit area?

[Leonardo Hardtke]

Hi!
I never used clustTool before.
 I made a try and i think I found the trick :-).
 You have to transform your (SpatialPointDataFrame) data to a dataframe...
Vine1.df=as.data.frame(Vine1)
 and then you can Set acctive data in the GUI!
hope it helps!

Leonardo
On 07/27/2010 09:06 PM, Etienne DELAY (deust 2010) wrote:

Le 27/07/10 20:35, Hamish a écrit :
   

Leonardo wrote:
 

you can run R in a grass terminal, and load data via
spgrass6 package,
readRAST6 or readVECT6.
   

see  http://grass.osgeo.org/wiki/R_statistics

 

Thank you Leonardo, but I read it, without understanding much more about
the use of clustTool... (it's ok whit readVECT ;-) )
I dont understand...

library(spgrass6)

Vine1- readVECT6(vineyardheroiq...@etienne)
library(clustTool)

it's run a graphical windows...
an I can't choose my data ;-(

so I need an other help

thank you

   


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Cluster whit area?

[Leonardo Hardtke]

On 07/28/2010 08:50 PM, Etienne DELAY (deust 2010) wrote:

I made my excuses to you, it works perfectly ... thank you very much
now I can go to bed happy :D

   

Glad to hear it worked!
Leonardo
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Cluster whit area?

[Leonardo Hardtke]

Hi,
you can run R in a grass terminal, and load data via spgrass6 package, 
readRAST6 or readVECT6.

Or export raster and read with rgdal.

Hope it helps.

Leonardo Hardtke
On 07/27/2010 01:30 PM, Etienne DELAY (deust 2010) wrote:

I'm come back :-)
I want to try clustTool in GRASS,
so I reset the region setting ,
I run R in grass,
I run library(clustTool)
That give mi a graphical environement, but I can't choose the data where
I want to work!!! I just can choose exemple data...

may be a basic question how to work on my data?
sorry for this noob question...:- S
   


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] r.what error

[Leonardo Hardtke]

Self updating
by doing
g.region rast=cleanNDVI_..
r.what . on the terminal gives the correct result...
but... wyh g.gui does it good without g.region???
Thanks
Leonardo Hardtke

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Compiling r.seg

[Leonardo Hardtke]

Hi,
I compiled r.seg in grass 6.5 (Ubuntu 10.04) but it must be the same in 
grass 6.4. Just copy r.seg directory to the raster subdirectory of the 
grass source code tree.

cd /usr/local/grass-6.4.0svn/ratster/r.seg
make
when done...
cd /usr/local/grass-6.4.0svn
make
sudo make install.

Hope it helps

On 06/17/2010 08:22 AM, katrin eggert wrote:

Hi there

I have tried to compile r.seg, a Segmentation developed by Trento 
University. I Know this is not part of GRASS6.4 and not even an addon 
but can you help me on this? (where I got the module 
http://www.ing.unitn.it/~grass/ http://www.ing.unitn.it/%7Egrass/)


This is the command:sudo make -C r.seg 
MODULE_TOPDIR=/usr/local/grass-6.4.0svn

The output:
make: Entering directory `/home/knrr/Desktop/GRASS_dvp/r.seg'
aseglib
make[1]: Entering directory `/home/knrr/Desktop/GRASS_dvp/r.seg/aseglib'
make[1]: Leaving directory `/home/knrr/Desktop/GRASS_dvp/r.seg/aseglib'
src
make[1]: Entering directory `/home/knrr/Desktop/GRASS_dvp/r.seg/src'
gcc -L/usr/local/grass-6.4.0svn/lib -Wl,--export-dynamic 
-Wl,-rpath-link,/usr/local/grass-6.4.0svn/lib-o 
/usr/local/grass-6.4.0svn/bin/r.seg OBJ.i686-pc-linux-gnu/main.o 
 -laseg -lgrass_gis -lgrass_datetime -lz -lm  -lz

make htmlcmd
make[2]: Entering directory `/home/knrr/Desktop/GRASS_dvp/r.seg/src'
make /usr/local/grass-6.4.0svn/docs/html/r.seg.html 
HTMLSRC=/usr/local/grass-6.4.0svn/bin/r.seg

make[3]: Entering directory `/home/knrr/Desktop/GRASS_dvp/r.seg/src'
if [ /usr/local/grass-6.4.0svn/bin/r.seg !=  ] ; then 
GISRC=/usr/local/grass-6.4.0svn/demolocation/grassrc64 
GISBASE=/usr/local/grass-6.4.0svn 
PATH=/usr/local/grass-6.4.0svn/bin:$PATH 
LD_LIBRARY_PATH=/usr/local/grass-6.4.0svn/bin:/usr/local/grass-6.4.0svn/lib: 
LC_ALL=C /usr/local/grass-6.4.0svn/bin/r.seg --html-description  
/dev/null | grep -v '/body\|/html'  r.seg.tmp.html ; true ; fi
/usr/local/grass-6.4.0svn/tools/mkhtml.sh r.seg ; mkdir -p 
/usr/local/grass-6.4.0svn/docs/html ; /usr/bin/install -c  -m 644 
r.seg.tmp.html /usr/local/grass-6.4.0svn/docs/html/r.seg.html ; for 
file in  *.png *.jpg ; do head -n 1 $file | grep '^#!'  /dev/null ; 
if [ $? -ne 0 ] ; then /usr/bin/install -c  -m 644 $file 
/usr/local/grass-6.4.0svn/docs/html ; fi done 2 /dev/null ; true

make[3]: Leaving directory `/home/knrr/Desktop/GRASS_dvp/r.seg/src'
make[2]: Leaving directory `/home/knrr/Desktop/GRASS_dvp/r.seg/src'
make mancmd
make[2]: Entering directory `/home/knrr/Desktop/GRASS_dvp/r.seg/src'
make /usr/local/grass-6.4.0svn/man/man1/r.seg.1 
MANSRC=/usr/local/grass-6.4.0svn/docs/html/r.seg.html

make[3]: Entering directory `/home/knrr/Desktop/GRASS_dvp/r.seg/src'
mkdir -p /usr/local/grass-6.4.0svn/man/man1
GRASS_PERL=/usr/bin/perl VERSION_NUMBER=6.4.0svn sh 
/usr/local/grass-6.4.0svn/tools/g.html2man/g.html2man 
/usr/local/grass-6.4.0svn/docs/html/r.seg.html 
/usr/local/grass-6.4.0svn/man/man1/r.seg.1 1

sh: Can't open /usr/local/grass-6.4.0svn/tools/g.html2man/g.html2man
make[3]: *** [/usr/local/grass-6.4.0svn/man/man1/r.seg.1] Error 2
make[3]: Leaving directory `/home/knrr/Desktop/GRASS_dvp/r.seg/src'
make[2]: *** [mancmd] Error 2
make[2]: Leaving directory `/home/knrr/Desktop/GRASS_dvp/r.seg/src'
make[1]: *** [cmd] Error 2
make[1]: Leaving directory `/home/knrr/Desktop/GRASS_dvp/r.seg/src'
make: Leaving directory `/home/knrr/Desktop/GRASS_dvp/r.seg'

Does it seems to all ok?
At my GRASS binary folder I don't have an executable 
(application/x-executable) but a KSeg document (application/x-kseg) 
that can be runned from GRASS. Is this normal/expected?


Thanks
Katrin


___
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] random subset in R

[Leonardo Hardtke]

Hi list,
I'm tryng to do a simple random subset of spatial points in R and don't 
know how

I tried:

library(spgrass6)
hs=readVECT6(hs...@hot_spot) (3500 point) #SpatialPointDataFrame
hs_smpl=spsample(hs, 200, random)  #-SpatialPoints

the question is... how do I attach the right data?
Thanks!

Leonardo Hardtke

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Python problem

[Leonardo Hardtke]

Hi!
I'm trying to develop some pyton code without luck(Grass64 on a Ubuntu
9.04 python 2.6)
I start using the example from the wiki
http://grass.osgeo.org/wiki/GRASS_and_Python, but I got error messages...

Traceback (most recent call last):
  File /usr/lib/grass64/etc/wxpython/gui_modules/menuform.py, line 1931,
in module
GrassGUIApp( grassTask( sys.argv[1] ) ).MainLoop()
  File /usr/lib/grass64/etc/wxpython/gui_modules/menuform.py, line 335, in
__init__
xml.sax.parseString( getInterfaceDescription( grassModule ) ,
processTask( self ) )
  File /usr/lib/grass64/etc/wxpython/gui_modules/menuform.py, line 1748,
in getInterfaceDescription
raise IOError, _(Unable to fetch interface description for command
'%s'.) % cmd
IOError: Unable to fetch interface description for command 'test.py'.

Any hint?
Thanks





-- 
Lic. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120 Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Problem with Grass 6.4 and R 2.9.2 on Ubuntu 9.04

[Leonardo Hardtke]

Hi,
I'm doing my first steps with the R-Grass interface but...

library(spgrass6)
Loading required package: sp
Loading required package: rgdal
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.6.1, released 2009/05/11
Path to GDAL shared files: /usr/share/gdal16
Loaded PROJ.4 runtime: Rel. 4.6.1, 21 August 2008
Path to PROJ.4 shared files: (autodetected)
Loading required package: XML
GRASS GIS interface loaded with GRASS version: 6.4.0RC5
and location: monte
 a=readVECT6(area_estudio)
ERROR: Incompatible library version for module. You need to rebuild GRASS
   or untangle multiple installations.

I installed grass65 from source, is this the problem? How do I fix it?
Thanks!


-- 
Lic. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120 Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] problem in Grass 6.4 rc3 Ubuntu 9.04

[Leonardo Hardtke]

Hi,
I'm having problems running Grass 6.4 rc3 on Ubuntu 9.04
When trying to open an existing grass database i see the following error.

*** glibc detected *** g.proj: double free or corruption (fasttop):
0x09b8fee0 ***
=== Backtrace: =
/lib/tls/i686/cmov/libc.so.6[0xb78b7604]
/lib/tls/i686/cmov/libc.so.6(cfree+0x96)[0xb78b95b6]
/usr/lib/libstdc++.so.6(_ZdlPv+0x21)[0xb6b4dad1]
/usr/lib/libstdc++.so.6(_ZNSs4_Rep10_M_destroyERKSaIcE+0x1d)[0xb6b2b6fd]
/usr/lib/libstdc++.so.6(_ZNSsD1Ev+0x4c)[0xb6b2d10c]
/lib/tls/i686/cmov/libc.so.6(__cxa_finalize+0xc1)[0xb7876f51]
/usr/local/lib/libgeos-3.1.1.so[0xb6984e64]
/usr/local/lib/libgeos-3.1.1.so[0xb6a2dde0]
/lib/ld-linux.so.2[0xb80e2a73]
/lib/tls/i686/cmov/libc.so.6(exit+0x89)[0xb7876bb9]
g.proj[0x804ab47]
/lib/tls/i686/cmov/libc.so.6(__libc_start_main+0xe5)[0xb785e775]
g.proj[0x8049ad1]
=== Memory map: 
08048000-0804d000 r-xp  08:1a 125207 /usr/lib/grass64/bin/g.proj
0804d000-0804e000 r--p 4000 08:1a 125207 /usr/lib/grass64/bin/g.proj
0804e000-0804f000 rw-p 5000 08:1a 125207 /usr/lib/grass64/bin/g.proj
09b8f000-09bb rw-p 09b8f000 00:00 0  [heap]
b590-b5921000 rw-p b590 00:00 0
b5921000-b5a0 ---p b5921000 00:00 0
b5a5-b5a8f000 r--p  08:1a 16817
/usr/lib/locale/es_AR.utf8/LC_CTYPE
b5a8f000-b5a93000 rw-p b5a8f000 00:00 0
b5a93000-b5a96000 r-xp  08:1a 2666   /lib/libgpg-error.so.0.3.0
b5a96000-b5a97000 rw-p 2000 08:1a 2666   /lib/libgpg-error.so.0.3.0
b5a97000-b5a99000 r-xp  08:1a 2672   /lib/libkeyutils-1.2.so
b5a99000-b5a9a000 r--p 1000 08:1a 2672   /lib/libkeyutils-1.2.so
b5a9a000-b5a9b000 rw-p 2000 08:1a 2672   /lib/libkeyutils-1.2.so
b5a9b000-b5a9c000 rw-p b5a9b000 00:00 0
b5a9c000-b5aa3000 r-xp  08:1a 10231
/usr/lib/libkrb5support.so.0.1
b5aa3000-b5aa4000 r--p 6000 08:1a 10231
/usr/lib/libkrb5support.so.0.1
b5aa4000-b5aa5000 rw-p 7000 08:1a 10231
/usr/lib/libkrb5support.so.0.1
b5aa5000-b5ac7000 r-xp  08:1a 10223  /usr/lib/libk5crypto.so.3.1
b5ac7000-b5ac8000 r--p 00022000 08:1a 10223  /usr/lib/libk5crypto.so.3.1
b5ac8000-b5ac9000 rw-p 00023000 08:1a 10223  /usr/lib/libk5crypto.so.3.1
b5ac9000-b5adf000 r-xp  08:1a 37009  /usr/lib/libsasl2.so.2.0.22
b5adf000-b5ae r--p 00015000 08:1a 37009  /usr/lib/libsasl2.so.2.0.22
b5ae-b5ae1000 rw-p 00016000 08:1a 37009  /usr/lib/libsasl2.so.2.0.22
b5ae1000-b5af3000 r-xp  08:1a 6269   /lib/tls/i686/cmov/
libresolv-2.9.so
b5af3000-b5af4000 r--p 00011000 08:1a 6269   /lib/tls/i686/cmov/
libresolv-2.9.so
b5af4000-b5af5000 rw-p 00012000 08:1a 6269   /lib/tls/i686/cmov/
libresolv-2.9.so
b5af5000-b5af7000 rw-p b5af5000 00:00 0
b5af7000-b5af9000 r-xp  08:1a 2647   /lib/libcom_err.so.2.1
b5af9000-b5afa000 r--p 1000 08:1a 2647   /lib/libcom_err.so.2.1
b5afa000-b5afb000 rw-p 2000 08:1a 2647   /lib/libcom_err.so.2.1
b5afb000-b5afc000 rw-p b5afb000 00:00 0
b5afc000-b5b8b000 r-xp  08:1a 10229  /usr/lib/libkrb5.so.3.3
b5b8b000-b5b8d000 r--p 0008e000 08:1a 10229  /usr/lib/libkrb5.so.3.3
b5b8d000-b5b8e000 rw-p 0009 08:1a 10229  /usr/lib/libkrb5.so.3.3
b5b8e000-b5cc1000 r-xp  08:1a 91452
/lib/i686/cmov/libcrypto.so.0.9.8
b5cc1000-b5cc9000 r--p 00132000 08:1a 91452
/lib/i686/cmov/libcrypto.so.0.9.8
b5cc9000-b5cd6000 rw-p 0013a000 08:1a 91452
/lib/i686/cmov/libcrypto.so.0.9.8
b5cd6000-b5cda000 rw-p b5cd6000 00:00 0
b5cda000-b5d1c000 r-xp  08:1a 91454
/lib/i686/cmov/libssl.so.0.9.8
b5d1c000-b5d1d000 ---p 00042000 08:1a 91454
/lib/i686/cmov/libssl.so.0.9.8
b5d1d000-b5d1e000 r--p 00042000 08:1a 91454
/lib/i686/cmov/libssl.so.0.9.8
b5d1e000-b5d21000 rw-p 00043000 08:1a 91454
/lib/i686/cmov/libssl.so.0.9.8
b5d21000-b67f7000 r--p  08:1a 10179  /usr/lib/libicudata.so.38.1
b67f7000-b67f8000 r--p 00ad5000 08:1a 10179  /usr/lib/libicudata.so.38.1
b67f8000-b6915000 r-xp  08:1a 10191  /usr/lib/libicuuc.so.38.1
b6915000-b6916000 ---p 0011d000 08:1a 10191  /usr/lib/libicuuc.so.38.1
b6916000-b691e000 r--p 0011d000 08:1a 10191  /usr/lib/libicuuc.so.38.1
b691e000-b691f000 rw-p 00125000 08:1a 10191  /usr/lib/libicuuc.so.38.1
b691f000-b6922000 rw-p b691f000 00:00 0
b6922000-b6a65000 r-xp  08:1a 242302 /usr/local/lib/
libgeos-3.1.1.so
b6a65000-b6a69000 r--p 00142000 08:1a 242302 /usr/local/lib/
libgeos-3.1.1.so
b6a69000-b6a6b000 rw-p 00146000 08:1a 242302 /usr/local/lib/
libgeos-3.1.1.so
b6a6b000-b6a95000 r-xp  08:1a 135513 /lib/libgcc_s.so.1
b6a95000-b6a96000 r--p 00029000 08:1a 135513 /lib/libgcc_s.so.1
b6a96000-b6a97000 rw-p 0002a000 08:1a 135513 /lib/libgcc_s.so.1
b6a97000-b6b7c000 r-xp  08:1a 135509
/usr/lib/libstdc++.so.6.0.12
b6b7c000-b6b8 r--p 000e4000 08:1a 135509
/usr/lib/libstdc++.so.6.0.12
b6b8-b6b81000 rw-p 000e8000 08:1a 135509
/usr/lib/libstdc++.so.6.0.12
b6b81000-b6b88000 

[GRASS-user] GIPE INSTALL UBUNTU

[Leonardo Hardtke]

Hi!
I'm not linux expert, and i couldn't install GIPE (and other add-ons) on my
Ubuntu 9.04 (GRASS 6.4 RC4 les-ejk repo)...

I downloaded grass source, (apt-get source  grass) and svn add-ons, then I
changed the first lines on the GIPE install script..

GRASSSVNDIR=/home/leomint/Descargas/fuente/grass-6.4.0~rc4
ADDONSVNDIR=/home/leomint/Descargas/fuente/grass-addons

when I ran the script, it doesn't worked

leom...@leo-cenpat ~/Descargas/fuente/grass-addons/gipe $ sudo
./install_yanns_grass.sh

Makefile.gipe:13: ../include/Make/Platform.make: No existe el fichero ó
directorio
Makefile.gipe:17: ../include/Make/Dir.make: No existe el fichero ó
directorio
make: *** No hay ninguna regla para construir el objetivo
`../include/Make/Dir.make'.  Alto.
cp: no se puede efectuar `stat' sobre
«/home/leomint/Descargas/fuente/grass-addons/i.*»: No existe el fichero ó
directorio
Makefile.gipe:3: ..//include/Make/Dir.make: No existe el fichero ó
directorio
make: *** No hay ninguna regla para construir el objetivo
`..//include/Make/Dir.make'.  Alto.
cp: no se puede crear el fichero regular
«/home/leomint/Descargas/fuente/gui/tcltk/gis.m/»: No existe el fichero ó
directorio
cp: no se puede crear el fichero regular
«/home/leomint/Descargas/fuente/gui/wxpython/gui_modules/»: No existe el
fichero ó directorio
./install_yanns_grass.sh: line 75: cd: /home/leomint/Descargas/fuente/gui:
No existe el fichero ó directorio
Makefile:8: ../../include/Make/Module.make: No existe el fichero ó
directorio
make: *** No hay ninguna regla para construir el objetivo
`../../include/Make/Module.make'.  Alto.
make: *** No hay ninguna regla para construir el objetivo `install'.  Alto.

Thanks in advise
Leonardo Hardtke


-- 
Lic. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120 Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] GIPE install ubuntu

[Leonardo Hardtke]

Hi,

i'm new to GRASS and not an linux expert. I need help to install GIPE. I
installed the GRASS 6.4 RC4 on UBUNTU 9.04 from the les-ejk and downloaded
the GIPE svn to my home directory, what shoud I do now?
Thanks!

-- 
Lic. Leonardo A. Hardtke
Laboratorio de Teledetección y S.I.G.
Centro Nacional Patagónico (CONICET)
Bvd. Brown 2825, 9120 Puerto Madryn, Chubut, Argentina
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] (sin asunto)

[Leonardo Hardtke]

Hi,
first, sorry about my english...
I have lot of points (over 60.000), some with the same attribute (Date), and
are very near (clustered). I don't now how to get only 1 point for each
cluster and whit its attributes. I followed the example of the archsites
from v.buffer, but it does not work well for me...


Thanks!
Leonardo
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user