[GRASS-user] r.basin matplotlib error

2015-04-30 Thread maning sambale
Hi,

I'm running r.basin using GRASS 6.4.4 from KyngChaos.

r.basin map=dem_90_int@PERMANENT prefix=out easting=712026.954
northing=1212651.537 threshold=40 --overwrite

I get this error, but the process continues.

out_dem_90_int_dist2out Done
##
Traceback (most recent call last):
  File /Users/maning/Library/GRASS/6.4/Modules/bin/r.ipso, line 57,
in module
import matplotlib #required by windows
ImportError: No module named matplotlib
##
##
Traceback (most recent call last):
  File /Users/maning/Library/GRASS/6.4/Modules/bin/r.wf, line 48, in module
import matplotlib #required by windows
ImportError: No module named matplotlib
##
Reading maps...

I was able to generate the maps and the Morphometric parameters.
What does this error mean and does it have any implication in the
Morphometric parameters results?



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
https://epsg4253.wordpress.com/
http://twitter.com/maningsambale
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] r.basin matplotlib error

2015-04-30 Thread maning sambale
Dear Margherita,

Thanks!

On Thu, Apr 30, 2015 at 4:35 PM, Margherita Di Leo direg...@gmail.com wrote:
 Matplotlib is a dependency needed for the graphic output, being the
 hypsographic curve and the width function. If you don't need those graphics,
 this error will not affect the result that you get. For what concerns the
 reason why this dependency is missing in windows, I don't know, perhaps
 packaging issue?

This is what I suspect.  I'm sure I have matplotlib installed (might
be some conflicting packages), will investigate further.
Not a biggie for not getting the graphics (would be nice to have of
course) as long as the results are OK.

Thanks again!



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
https://epsg4253.wordpress.com/
http://twitter.com/maningsambale
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] SRTM 1-arc (~30m) available outside USA and Africa

2015-01-13 Thread maning sambale
Just noticed it today and would like to share.  Last year, NASA/JPL
announced that they will finally release SRTM ~30m to areas outside
the USA starting with Africa.

I noticed today that data is now available for my country.  Might be
available for your areas too.

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Issues with resampling to equal pixel size for FCELL rasters

2014-12-22 Thread maning sambale
Hi,

I have a DTM I want to resample. The original data was in lat/lon SRS
which I re-projected to UTM and imported to GRASS GIS.

GRASS 6.4.3 (pampanga_wshed):~  r.info ifsar_dtm
 ++
 | Layer:ifsar_dtm  Date: Mon Dec 22 10:07:06 2014|
 | Mapset:   PERMANENT  Login of Creator: maning  |
 | Location: pampanga_wshed   |
 | DataBase: /Users/maning/GRASSDATA/ |
 | Title: ( ifsar_dtm )   |
 | Timestamp: none|
 ||
 ||
 |   Type of Map:  raster   Number of Categories: 255 |
 |   Data Type:FCELL  |
 |   Rows: 0  |
 |   Columns:  53950  |
 |   Total Cells:  2397538000 |
 |Projection: UTM (zone 51)   |
 |N: 1842577.192205S: 1620332.25097015   Res: 5.00101128  |
 |E: 420942.75364933W: 151243.07137293   Res: 4.99906733  |
 |   Range of data:min = -1  max = 3.402823e+38   |
 ||
 |   Data Description:|
 |generated by r.in.gdal  |
 ||
 |   Comments:|
 |r.in.gdal -e input=ifsar_dtm_utmz51n.vrt output=ifsar_dtm   |
 ||
 ++

Because of the origin SRS, the pixel size is not the same for ew_res and nw_res.
I tried running several resampling techniques to match the pixel size
equal to 5.

First, I tried r.resamp.interp.

g.region rast=ifsar_dtm_subset -pm
projection: 1 (UTM)
zone:   51
datum:  wgs84
ellipsoid:  wgs84
north:  1683792
south:  1645429
west:   220827
east:   268500
nsres:  5.00104289
ewres:  4.99926594
rows:   7671
cols:   9536
cells:  73150656

g.region res=5 -ap
projection: 1 (UTM)
zone:   51
datum:  wgs84
ellipsoid:  wgs84
north:  1683795
south:  1645425
west:   220825
east:   268500
nsres:  5
ewres:  5
rows:   7674
cols:   9535
cells:  73171590

r.resamp.interp input=ifsar_dtm_subset
output=ifsar_dtm_subset_interp_bic method=bicubic --overwrite
r.info ifsar_dtm_subset_interp_bic
 ++
 | Layer:ifsar_dtm_subset_interp_bicDate: Tue Dec 23 13:17:51 2014|
 | Mapset:   test1  Login of Creator: maning  |
 | Location: pampanga_wshed   |
 | DataBase: /Users/maning/GRASSDATA/ |
 | Title:Resample by bicubic interpolation ( ifsar_dtm_subset_interp_bic  |
 | Timestamp: none|
 ||
 ||
 |   Type of Map:  raster   Number of Categories: 255 |
 |   Data Type:DCELL  |
 |   Rows: 7674   |
 |   Columns:  9535   |
 |   Total Cells:  73171590   |
 |Projection: UTM (zone 51)   |
 |N:1683795S:1645425   Res: 5 |
 |E: 268500W: 220825   Res: 5 |
 |   Range of data:min = -2.72093853362067  max = 1157.05376343899|
 ||
 |   Data Source: |
 |ifsar_dtm_subset|
 |Source map NS res: 5.00104289   EW res: 4.99926594  |
 ||
 |   Data Description: 

Re: [GRASS-user] 10day accumulation which t.rast module to use?

2014-12-04 Thread maning sambale

 Oh, I see... then, t.rast.accumulate is what you need. The default method
 mean does the sum.

So mean is actually sum?


 Maybe, something like this (please re-check the cycle and granularity
 parameters, i didn't test)

 t.rast.accumulate input=3hr_rainrate output=10day_acum_precip
 base=10day_acc_prec start=2001-01-01 end=2001-01-20 cycle=10 days
 granularity=1 days method=mean
I tried this, end should be stop

t.rast.accumulate --overwrite input=10day_hourly
output=10day_accumulation base=10day_accumulation \
start=2014-11-04 11:00:00 stop=2014-12-04 11:00:00 cycle=10 days
granularity=1 days method=mean

The process wouldn't stop, my area is very small so I expect it to
finish in less than an hour:
 g.region -p
projection: 3 (Latitude-Longitude)
zone:   0
datum:  wgs84
ellipsoid:  wgs84
north:  22N
south:  4N
west:   115E
east:   155E
nsres:  0:15
ewres:  0:15
rows:   72
cols:   160
cells:  11520

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] 10day accumulation which t.rast module to use?

2014-12-04 Thread maning sambale
Dear Soeren,


 The temporal raster algebra module is very experimental, but i hope
 that this will work.
Thanks will try this.


-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] GRASS71 development packages for Ubuntu 13.10 and 14.04

2014-12-04 Thread maning sambale
GRASS71 binaries at launchpad's grass/grass-devel is only for Ubuntu 12.04.1.
Any chance to get one for 13.10 and 14.04.  Maybe from other repo?

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] 10day accumulation which t.rast module to use?

2014-12-03 Thread maning sambale
Hi,

I have 3hourly rainrate data for a given region in a GRASS temporal
database (strds).

I want to aggregate the data by summing all raster from the previous
10 days to calculate the accumulated rainfall.
The temporal modules have several tools:
t.rast.series
t.rast.accumulate
t.rast.aggregate
t.rast.mapcalc

From the manual, the appropriate module seems to be t.rast.accumulate
but it does not have sum in the method parameters. Currently, I'm
running several r.series commands for each day like below:

r.series input=\
`t.rast.list input=3hr_rainrate where=\
('${DATE} 00:00:00' = datetime(start_time, '3 hours')) \
and ('${DATE} 00:00:00' = datetime(start_time, '10 days'))  \
col=name method=comma` \
output=test_2014-07-20 method=sum

Then, register each layer in another strds:

#Create a new STRDS
t.create output='10day_accumulation' type='strds' temporaltype='absolute' \
title='10day_accumulation' \
description=10 day rainfall accumulation


g.mlist type='rast', pattern='test*', separator=','

#Register the r.series output to the strds
t.register -i type=rast input=10day_accumulation \
start=2014-07-10 00:00:00 increment=1 days separator=, \
maps=`g.mlist type='rast', pattern='test*', separator=','`


Is there another approach using directly any of the t.rast.* modules.
For a time-line visualization, here's what I want to do:
https://dl.dropboxusercontent.com/u/2096185/sample_accumulation.png

Thanks!

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] 10day accumulation which t.rast module to use?

2014-12-03 Thread maning sambale
Dear Vero,

Thanks for the reply.  But t.rast.aggregate does not do what I intend
to do.  My rainrate is from a span of 20 days.
t.rast.aggregate produces (as expected) two rasters, one for each of
the 10day accumulation.  What I want for example is like a moving sum,
where for each day, it calculates the sum for the past 10 days.

On Thu, Dec 4, 2014 at 4:31 AM, Veronica Andreo veroand...@gmail.com wrote:
 Hi Emmanuel,

 Hi,

 I have 3hourly rainrate data for a given region in a GRASS temporal
 database (strds).

 I want to aggregate the data by summing all raster from the previous
 10 days to calculate the accumulated rainfall.
 The temporal modules have several tools:
 t.rast.series
 t.rast.accumulate
 t.rast.aggregate
 t.rast.mapcalc

 From the manual, the appropriate module seems to be t.rast.accumulate
 but it does not have sum in the method parameters. Currently, I'm
 running several r.series commands for each day like below:

 r.series input=\
 `t.rast.list input=3hr_rainrate where=\
 ('${DATE} 00:00:00' = datetime(start_time, '3 hours')) \
 and ('${DATE} 00:00:00' = datetime(start_time, '10 days'))  \
 col=name method=comma` \
 output=test_2014-07-20 method=sum

 Then, register each layer in another strds:

 #Create a new STRDS
 t.create output='10day_accumulation' type='strds' temporaltype='absolute'
 \
 title='10day_accumulation' \
 description=10 day rainfall accumulation


 g.mlist type='rast', pattern='test*', separator=','

 #Register the r.series output to the strds
 t.register -i type=rast input=10day_accumulation \
 start=2014-07-10 00:00:00 increment=1 days separator=, \
 maps=`g.mlist type='rast', pattern='test*', separator=','`


 Is there another approach using directly any of the t.rast.* modules.
 For a time-line visualization, here's what I want to do:
 https://dl.dropboxusercontent.com/u/2096185/sample_accumulation.png

 Thanks!

 If i understood right, something like

 t.rast.aggregate input=3hr_rainrate output=10day_rain basename=10day_rain
 granularity='10 days' method=sum

 should do the what you need.

 HTH,

 Vero



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] temporal mapcalc 3hour precipitation to hourly

2014-12-02 Thread maning sambale
Hi,

A have a temporal raster of 3hourly_rainrate.
 t.rast.list input=3hr_rainrate
name|mapset|start_time|end_time
3hr_rainrate_2014-11-13_0300|PERMANENT|2014-11-13 03:00:00|2014-11-13 06:00:00
3hr_rainrate_2014-11-13_0600|PERMANENT|2014-11-13 06:00:00|2014-11-13 09:00:00

The values are the average hourly rainfall at 3hr snapshot.  I want to
get the total rainfall which is (raster * 3).  Using t.rast.mapcalc, I
get this error:

 t.rast.mapcalc input=3hr_rainrate output=hourly_rainrate 
 basename=hourly_rainrate expression=(3hr_rainrate) * (3.0) nprocs=4
Start temporal sampling
Start mapcalc computation
Invalid map 3hr_rainrate_2014
Invalid map 13_0300
Parse error
ERROR: parse error
...
ERROR: parse error
Invalid map 3hr_rainrate_2014
Invalid map 13_1200
Parse error
ERROR: parse error
ERROR: Error while mapcalc computation

I'm using GRASS 7.0.0.  Thanks!



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] t.rast.accumulate ValueError: could not convert string to float:

2014-11-12 Thread maning sambale
Hi,

I have an 3hr raintrate data in GRASS STRDs. I tested running a 10day
accumulation but I get valuerror.
Details below:

GRASS 7.0.0 (rainfall):~/trmm  t.info 3hr_rainrate
 + Space Time Raster Dataset -+
 ||
 + Basic information -+
 | Id:  3hr_rainrate@PERMANENT
 | Name: .. 3hr_rainrate
 | Mapset:  PERMANENT
 | Creator: ... maning
 | Temporal type: . absolute
 | Creation time: . 2014-11-12 18:44:43.567312
 | Modification time:.. 2014-11-12 19:04:43.527267
 | Semantic type:.. mean
 + Absolute time -+
 | Start time:. 2014-07-09 00:00:00
 | End time:... 2014-07-20 03:00:00
 | Granularity: 3 hours
 | Temporal type of maps:.. interval
 + Spatial extent +
 | North:.. 22.0
 | South:.. 4.0
 | East:..  155.0
 | West:... 115.0
 | Top: 0.0
 | Bottom:. 0.0
 + Metadata information --+
 | Raster register table:..
raster_map_register_491d12cc26d2460d9591bd699f6beb43
 | North-South resolution min:. 0.25
 | North-South resolution max:. 0.25
 | East-west resolution min:... 0.25
 | East-west resolution max:... 0.25
 | Minimum value min:.. 0.0
 | Minimum value max:.. 0.0
 | Maximum value min:.. 9.39
 | Maximum value max:.. 38.32
 | Aggregation type:... None
 | Number of registered maps:.. 89
 
 |
 ++

GRASS 7.0.0 (rainfall):~/trmm  t.rast.accumulate
input=3hr_rainrate@PERMANENT output=10day_accumulation
start=2014-07-09 granularity='1 days' cycle='10 days'
basename=10day_accum
Processing cycle 2014-07-09 00:00:00 - 2014-07-19 00:00:00
Traceback (most recent call last):
  File /usr/lib/grass70/scripts/t.rast.accumulate, line 507, in module
main()
  File /usr/lib/grass70/scripts/t.rast.accumulate, line 393, in main
limits_lower = float(limits_vals[0])
ValueError: could not convert string to float:


-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Importing netCDF where raster names is the time variable

2014-11-11 Thread maning sambale
Hi,

I'm trying to import multidimensional netCDF files slicing the bands
according to time.  I want to use the new temopral module in GRASS70
for further processing.

My process are as follows:

1. Download the netCDF usinf ncks, i.e.:
$ ncks -O -v r -d time,${MIN},${MAX} -d lon,115.00,155.00 -d lat,4.00,22.00 \
http://disc2.nascom.nasa.gov:80/dods/3B42RT_V7_rainrate \
3hourly_trmm_2.nc

This is a 3hourly rainfal  data.

2. Running ncdump, the data are as follows:
$ ncdump -ct 3hourly_trmm_2.nc
netcdf 3hourly_trmm_2 {
dimensions:
lat = 72 ;
lon = 160 ;
time = 89 ;
variables:
double lat(lat) ;
lat:grads_dim = y ;
lat:grads_mapping = linear ;
lat:grads_size = 480 ;
lat:units = degrees_north ;
lat:long_name = latitude ;
lat:minimum = -59.875 ;
lat:maximum = 59.875 ;
lat:resolution = 0.25f ;
double lon(lon) ;
lon:grads_dim = x ;
lon:grads_mapping = linear ;
lon:grads_size = 1440 ;
lon:units = degrees_east ;
lon:long_name = longitude ;
lon:minimum = 0.125 ;
lon:maximum = 359.875 ;
lon:resolution = 0.25f ;
float r(time, lat, lon) ;
r:_FillValue = -9.f ;
r:missing_value = -9.f ;
r:long_name = hourly rain rate (mm/hr)  ;
double time(time) ;
time:grads_dim = t ;
time:grads_mapping = linear ;
time:grads_size = 131078 ;
time:grads_min = 00z01mar2000 ;
time:grads_step = 3hr ;
time:units = days since 1-1-1 00:00:0.0 ;
time:long_name = time ;
time:minimum = 00z01mar2000 ;
time:maximum = 15z08jan2045 ;
time:resolution = 0.125f ;

// global attributes:
:title = Three Hourly TRMM and Other Satellite Rainfall (3B42RT) ;
:Conventions = COARDS\n,
GrADS ;
:dataType = Grid ;
:history = Tue Nov 11 15:30:44 2014: ncks -O -v r -d
time,41945,42033 -d lon,115.00,155.00 -d lat,4.00,22.00
http://disc2.nascom.nasa.gov:80/dods/3B42RT_V7_rainrate
3hourly_trmm_2.nc\n,
Mon Nov 10 23:05:02 GMT-00:00 2014 : imported by GrADS
Data Server 2.0 ;
:NCO = 4.4.2 ;
data:

 lat = 4.125, 4.375, 4.625, 4.875, 5.125, 5.375, 5.625, 5.875, 6.125, 6.375,
6.625, 6.875, 7.125, 7.375, 7.625, 7.875, 8.125, 8.375, 8.625, 8.875,
 ... ;

 lon = 115.125, 115.375, 115.625, 115.875, 116.125, 116.375, 116.625,
116.875, 117.125, 117.375, 117.625, 117.875, 118.125, 118.375, 118.625,
  ... ;

 time = 2014-07-09 03, 2014-07-09 06, 2014-07-09 09, 2014-07-09 12,
2014-07-09 15, 2014-07-09 18, 2014-07-09 21, 2014-07-10,
 ... ;
}

3. I want to import it in a GRASS database where the raster filename
correspond to the date/time.

script.run_command('r.in.gdal', input='3hourly_trmm_2.nc', output='g',
flags='o')
script.list_pairs(type='rast')

[('g.1', 'PERMANENT'),
...
 ('g.85', 'PERMANENT')]

As you can see, the band/raster names were just an autoincrement of
the netCDF layers.
Any advice on doing this in Python?  Thanks!


-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Archiving GPM 3hr rainrate in GRASS

2014-11-06 Thread maning sambale
Thanks Markus.  Already registered and can access ftp.  Unfortunately,
processed (L3) rainfall data will be released by Dec 2014.
Will just wait then. :)

On Thu, Nov 6, 2014 at 4:27 PM, Markus Neteler nete...@osgeo.org wrote:
 On Mon, Nov 3, 2014 at 10:18 AM, maning sambale
 emmanuel.samb...@gmail.com wrote:
 Has anyone here able to archive and load GPM dataset into GRASS?  I
 was able to get TRMM netcdf in my area of interest before, but I can't
 find the tools to automate GPM downloads.

 Thanks!

 [0] http://www.nasa.gov/mission_pages/GPM/main/

 They state

 All data are freely available through the NASA's Precipitation
 Processing System at http://pps.gsfc.nasa.gov;

 -- Register and search for GPM and TRMM data, order custom subsets
 and set up subscriptions using PPS Data Products Ordering Interface
 (STORM).

 At time the server seems to be down?

 Markus



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Archiving GPM 3hr rainrate in GRASS

2014-11-03 Thread maning sambale
Has anyone here able to archive and load GPM dataset into GRASS?  I
was able to get TRMM netcdf in my area of interest before, but I can't
find the tools to automate GPM downloads.

Thanks!

[0] http://www.nasa.gov/mission_pages/GPM/main/
-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Possible to use vector attribute for r.carve's depth parameter

2014-10-13 Thread maning sambale
Hi,

I have a series of vector polygon layer where there is an attribute
for estimated depth of the rivers/streams.  I want to use this in
r.carve [0] to improve the river/stream channel of my DEM.
As I understand it, r.carve [0] only uses a single float value to
carve the stream channel.

Instead of looping a series of r.carve run for each depth attribute,
is it possible to just use the attribute of the vector polygon?

Thanks!

[0] http://grass.osgeo.org/grass70/manuals/r.carve.html

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] r.hazard.flood does not respect region settings?

2014-03-06 Thread maning sambale
Hi,

I'm testing r.hazard.flood and noticed that it computes the flood and
mti layers on the full region of the elevation raster instead of the
pre-defined region settings.

The relvant code I found from the r.hazard.flood is this:
   # Detect cellsize of the DEM
info_region = grass.read_command('g.region', flags = 'p', rast =
'%s' % (r_elevation))
dict_region = grass.parse_key_val(info_region, ':')
resolution = (float(dict_region['nsres']) + float(dict_region['ewres']))/2
grass.message(Cellsize : %s  % resolution)

Would it be possible to either respect the current region settings or
add a flag to choose between the dem region settings or current region
settings?


-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] GRASS7 fftw-3 error using r.sim.water

2013-09-06 Thread maning sambale
I downloaded grass7 from Michael Barton's page [0] and then tried r.sim.water:

GRASS 7.0.svn (dtm):~  r.sim.water elevin=dtm_1m dxin=dtm_1m_dx
dyin=dtm_1m_dy rain=rain manin=manning infil=infilt depth=depth
disch=discharge err=error --overwrite --verbose
dyld: Library not loaded: /Library/Frameworks/FFTW3.framework/Versions/3/FFTW3
  Referenced from:
/Applications/GRASS/GRASS-7.0.app/Contents/MacOS/lib/libgrass_gmath.7.0.svn.dylib
  Reason: image not found
Trace/BPT trap: 5





[0] http://grassmac.wikidot.com/
-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Graphical Modeler Output

2012-11-06 Thread maning sambale
Dear Paul,

Do you have x2 monitor activated?

d.mon start=x2
d.monsize setmonitor=x2 setwidth=601 setheight=391




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] implementing frost filter in grass

2012-08-14 Thread maning sambale
Hi,

I'm trying to implement the FROST [0] SAR despeckling filter in a
GRASS script.  The code according to the referenced link is below:

#
The implementation of this filter consists of defining a circularly
symmetric filter with a set of weighting values M for each pixel:

M = exp(- A * T)

where:

   A = DAMP * (V / I**2)

   T   is the absolute value of the pixel distance
   from the centre pixel to its neighbours in
   the filter window

   DAMPis the exponential damping factor

   V   is the variance of the grey-level in
   the filter window

   I**2is the square of the mean grey level in the
   filter window

The resulting grey-level value R for the smoothed pixel is:

R = (P1*M1 + P2*M2 + ... +  Pn*Mn) / (M1 + M2 + ... + Mn)

where:

P1 .. Pn   are grey levels of each pixel in filter window

M1 .. Mn   are weights (as defined above) for each pixel
#


Computing for A is straightforward using r.mapcalc.  I'm having
difficulty in deriving T which was defined as the absolute value of
the pixel distance
from the centre pixel to its neighbours in the filter window

I investigated r.cost and r.distance but it is not suited for a moving
window operation (perhaps I'm wrong).  It maybe possible  r.mapcalc's
row-column offset
but it will loose the ability to interactively define a window size.

Any tips on how to get the T and final R values?  Advance thanks!


[0] 
http://web.archive.org/web/20071216081728/http://www.pcigeomatics.com/cgi-bin/pcihlp/FFROST

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] creating custom filter with r.mfilter

2012-08-06 Thread maning sambale
Any working code I can start with to create custom filters?  I
basically want to run the speckle filters in GRASS.  Sample filters
are explained here:
http://resources.arcgis.com/en/help/main/10.1/index.html#//009t01z700

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] r.forestfrag AddOn

2012-03-02 Thread maning sambale
Dear Stefan,

If i remember it right, the window size is an input to the r.neighbors command.
In the script there is an option to slect a windows size:

#%option
#% key: window
#% type: integer
#% description: window size default is 5
#% answer : 5
#% required : yes
#%END

You can modify the r.neighbors section like this:

r.neighbors input=A output=C method=sum size=$GIS_OPT_window --o
r.neighbors input=B output=D method=sum size= $GIS_OPT_window --o

On Sun, Feb 26, 2012 at 5:51 PM, Stefan Sylla stefansy...@gmx.de wrote:
 By now, the script can only calculate the forest fraction for a 3x3 window.
 Any ideas on how to automatically generate different window-sizes would be
 welcome.



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Re: grass and palsar processing

2012-02-15 Thread maning sambale
Thanks markus for this.

I added an initial workflow for the ENVI tutorial I attended my goal
is to port this to GRASS.  Most of the process is already available
in the pdf Alexander provided.

However, I'm having difficulty importing images from the KC mosaic
product [0].  In the training I attended the mosaic product was
converted to GTiff using ASF MapReady [1].  According to ASF MapReady
the format is CEOS L1.

I tried using r.in.gdal and r.in.bin (based on the pdf) but I cannot
import the KC mosaic


[0] http://www.eorc.jaxa.jp/ALOS/en/kc_mosaic/kc_mosaic.htm
[1] http://www.asf.alaska.edu/downloads/software_tools/mapready/version_history


On Sat, Feb 11, 2012 at 7:54 AM, Markus Neteler nete...@osgeo.org wrote:
 Hi,

 I am almost there:http://grass.osgeo.org/wiki/Radar_data_processing
 only imagemagick currently fails to make the thumbnails after debian
 upgrade.will check asap and upload the missing files.
 On Fri, Feb 10, 2012 at 7:43 PM, Alexander Muriy amu...@gmail.com wrote:
 Also there is PDF http://gis-lab.info/forum/download/file.php?id=2370,
 copied from http://www.safe.iis.u-tokyo.ac.jp/, about SAR  GRASS.

 thanks, added on the Wiki page.

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



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] grass and palsar processing

2012-02-10 Thread maning sambale
Hi,

Just a heads up if there are GRASS gurus who are using radar
(specifically PALSAR) for land cover classification.  I just got back
from a radar workshop using ENVI 4.7.  I am interested in porting
the methodology (import, pre-processing, classification) to GRASS.

Just looking for specific people I can collaborate with in the future. :)

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] grass and palsar processing

2012-02-10 Thread maning sambale
Sweet! Thanks markus.
On Feb 10, 2012 7:16 PM, Markus Neteler nete...@osgeo.org wrote:

 On Fri, Feb 10, 2012 at 12:15 PM, Martin Landa landa.mar...@gmail.com
 wrote:
  2012/2/10 Markus Neteler nete...@osgeo.org:
  Seems they have done it:
 
 http://www.osgeo.jp/wiki/index.php?title=Introduction_to_Radar_data_processing
 
  could be this page move to http://grass.osgeo.org/wiki? Seems to be
  very useful article...

 Will do that!

 Markus

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


Re: [GRASS-user] measuring area on a 3d polygon

2011-11-16 Thread maning sambale
Thanks, will try this.

On Thu, Nov 17, 2011 at 3:41 AM, Benjamin Ducke bendu...@fastmail.fm wrote:
 I think Michael's reply (see below) to my message is
 actually the answer to your question ;)

 Ben

 On Wednesday, November 16, 2011 10:13 AM, Michael Barton
 michael.bar...@asu.edu wrote:
 r.surf.area

 will calculate a topographically corrected area for rasters. So if you
 convert your vector areas to raster, this will do the trick.

 Michael
 __
 C. Michael Barton
 Director, Center for Social Dynamics  Complexity
 Professor of Anthropology, School of Human Evolution  Social Change
 Arizona State University
 Tempe, AZ  85287-2402
 USA

 voice:  480-965-6262 (SHESC), 480-727-9746 (CSDC)
 fax:          480-965-7671(SHESC), 480-727-0709 (CSDC)
 www:    http://csdc.asu.edu, http://shesc.asu.edu
               http://www.public.asu.edu/~cmbarton


 --
 Benjamin Ducke
 {*} Geospatial Consultant
 {*} GIS Developer

  bendu...@fastmail.fm


 On Wednesday, November 16, 2011 11:01 AM, maning sambale
 emmanuel.samb...@gmail.com wrote:
 Hi,

 I have a series polygons located in a mountainous terrain.  I can
 convert the 2d vector into 2.5d using v.drape and a DEM, however, it
 only converts 2D vector point or line data and not areas.  My task is
 to calculate the area in hectares of the polygon incorporating the
 elevation/surface of the polygon.  I then need to subdivide the main
 polygon into smaller polygons using a preferred area in hectares.

 How to calculate 2.5d polygon areas in grass?

 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --
 ___
 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




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] measuring area on a 3d polygon

2011-11-15 Thread maning sambale
Hi,

I have a series polygons located in a mountainous terrain.  I can
convert the 2d vector into 2.5d using v.drape and a DEM, however, it
only converts 2D vector point or line data and not areas.  My task is
to calculate the area in hectares of the polygon incorporating the
elevation/surface of the polygon.  I then need to subdivide the main
polygon into smaller polygons using a preferred area in hectares.

How to calculate 2.5d polygon areas in grass?

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] osgeo philippines ailing list

2011-10-10 Thread maning sambale
[Apologies for cross-posting]

Just a short announcement for anyone interested.

Public mailinglist for the OSGeo-Philippine local chapter is available at:

http://lists.osgeo.org/mailman/listinfo/philippines

This will be the main discussion list for organizing the OSGeo
Philippine Local Chapter.  See you there!

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] unable to open GRASS mapset under Q-GIS 1.6

2011-06-01 Thread maning sambale
In your file browser, got to your grass database and delete the
.gislock file (you need to show hidden files to see the .gislock or
any file with a preceding dot.


On Wed, Jun 1, 2011 at 6:31 PM, Niels Thevs th...@uni-greifswald.de wrote:
 Dear all,

 I want to digitize shapefiles in Q-GIS based on raster layers in a GRASS
 mapset. But, Q-GIS always complains that the mapset canno be opened, because
 it is used by another programme. I tried on several computers, but I only
 was able to open a GRASS mapset under Q-GIS on a computer, where I never
 started GRASS before. How could I solve this problem. I work under Windows
 XP with Q-GIS 1.6.

 Regards

 Niels

 --



 -
 Dr. Niels Thevs
 Institute of Botany and Landscape Ecology
 Greifswald University
 Grimmer Strasse 88
 17487 Greifswald
 Germany

 Tel.: +49-3834-86-4131
 Fax:  +49-3834-86-4114

 -


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




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Convert Floating Point Raster to Integer Raster

2011-05-18 Thread maning sambale
r.mapcalc integer_map = round(fp_maps)

On Thu, May 19, 2011 at 7:04 AM, Michael Morin mmo...@cassblue.com wrote:
 Techies,

 In ArcInfo I used to write .aml code to covert Floating Point Grids to 
 Integer Grids.  Doing this produced a .vat table in the new Integer grid 
 thereby giving me elevation data which i could then used to develop painted 
 DEMS.

 Question...?  Would r.mapcalc or r.reclass be the tools I would 
 use?..could someone share some example code how I would make this 
 conversion?

 I will sum.

 Bonjour.

 Michael D. Morin
 mmo...@cassblue.com





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




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Re: Produce KML/KMZ files in GRASS

2011-01-14 Thread maning sambale
I think its worth pointing out is that according to the ogr docs,
gdal/ogr automatically re-projects the data in EPSG:4326.  I'm not
sure though if this is correct for all cases (in the past I always
make sure that my data in GRASS is in EPSG:4326 before converting to
KML).

http://www.gdal.org/ogr/drv_kml.html
KML has been accepted as an OGC standard, and is supported in one way
or another on the major GeoBrowsers. Note that KML by specification
uses only a single projection, EPSG:4326. All OGR KML output will be
presented in EPSG:4326. As such OGR will create layers in the correct
coordinate system and transform any geometries.

On Fri, Jan 14, 2011 at 4:33 PM, Stefan fabios...@gmx.de wrote:

 Monica,
 if you're using GRASS 6.4, you can also use v.out.ogr with format=KML.
 It works perfectly for google-earth.
 All the best
 Stefan
 --
 View this message in context: 
 http://osgeo-org.1803224.n2.nabble.com/Produce-KML-KMZ-files-in-GRASS-tp5917938p5921006.html
 Sent from the Grass - Users mailing list archive at Nabble.com.
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: loading kc50 palsar datasets in GRASS

2010-11-25 Thread maning sambale
Follow-up on the, from the website:
*2. HH and HV polarization data are kept in BSQ format. Color
composition of the image data is R=HH, G=HV, B=HH/HV and enhanced.

In gdal docs reference to bsq is in the envi format:
http://www.gdal.org/frmt_various.html#ENVI

On Wed, Nov 24, 2010 at 6:21 PM, maning sambale
emmanuel.samb...@gmail.com wrote:
 Anyone successfully loaded kc50 palsar raw imagery with gdal?


 -- Forwarded message --
 From: maning sambale emmanuel.samb...@gmail.com
 Date: Tue, Nov 23, 2010 at 4:00 PM
 Subject: loading kc50 palsar datasets in GRASS
 To: grass-us...@lists.osgeo.org


 This site provides PALSAR 50m Orthorectified Mosaic Product.
 http://www.eorc.jaxa.jp/ALOS/en/kc_mosaic/kc_mosaic.htm

 I tried running gdalinfo to the directory, img and header files, but
 gdal doesn't recognize the format.  There is pdf (below the webpage)
 but I can't find any reference to the file formats that is supported
 by GDAL and GRASS.  Has anyone used this dataset in GRASS?

 PS.  This is the only dataset without clouds I can find for my study
 area (for obvious reasons :)).
 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --



 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] loading kc50 palsar datasets in GRASS

2010-11-23 Thread maning sambale
This site provides PALSAR 50m Orthorectified Mosaic Product.
http://www.eorc.jaxa.jp/ALOS/en/kc_mosaic/kc_mosaic.htm

I tried running gdalinfo to the directory, img and header files, but
gdal doesn't recognize the format.  There is pdf (below the webpage)
but I can't find any reference to the file formats that is supported
by GDAL and GRASS.  Has anyone used this dataset in GRASS?

PS.  This is the only dataset without clouds I can find for my study
area (for obvious reasons :)).
-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Meta Data Management Suggestions

2010-10-06 Thread maning sambale
For individual layers in GRASS, you can add descriptive information
and processing history using r.support.
To access the information, use r.info.

Many metadata support files in GRASS are simple text files.

A simple shell script can parse that data and reformatted into a human
readable format.



On Thu, Oct 7, 2010 at 10:20 AM, stephen sefick ssef...@gmail.com wrote:
 I would like to keep metadata files for the maps that I import into
 GRASS and also the maps that I create.  Much like the metadata that is
 provided by ... cough... ESRI.  Is there a way to do this in GRASS?
 Should I just store the meta data files as text files in the GRASS
 location and then edit them with the additional GRASS processing
 information?
 thanks,

 --
 Stephen Sefick
 
 | Auburn University                                   |
 | Department of Biological Sciences           |
 | 331 Funchess Hall                                  |
 | Auburn, Alabama                                   |
 | 36849                                                    |
 |___|
 | sas0...@auburn.edu                             |
 | http://www.auburn.edu/~sas0025             |
 |___|

 Let's not spend our time and resources thinking about things that are
 so little or so large that all they really do for us is puff us up and
 make us feel like gods.  We are mammals, and have not exhausted the
 annoying little problems of being mammals.

                                 -K. Mullis

 A big computer, a complex algorithm and a long time does not equal science.

                               -Robert Gentleman

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





-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] MDOW relief shading in GRASS

2010-09-29 Thread maning sambale
Hi,

I am revisiting a previous script I made in producing
Multidirectional, oblique-weighted, shaded-relief.  ArcGis has this
script:
http://blogs.esri.com/Support/blogs/mappingcenter/archive/2008/10/07/updated-hillshade-toolbox.aspx

The commads in ARC 6.0.1 GRID by Rober Mark (USGS) are as follows:

shade225 = hillshade(hawaii, 225, 30, shade, 5)
shade270 = hillshade(hawaii, 270, 30, shade, 5)
shade315 = hillshade(hawaii, 315, 30, shade, 5)
shade360 = hillshade(hawaii, 360, 30, shade, 5)
h00 = resample(hawaii, 1000)
h01 = focalmean(h00)
h02 = focalmean(h01)
h03 = focalmean(h02)
asp = aspect(h03)
asp1 = con(isnull(asp), 293, asp)
w225 = sqr(sin((asp1 - 225) div deg ))
w270 = sqr(sin((asp1 - 270) div deg))
w315 = sqr(sin((asp1 - 315) div deg ))
w360 = sqr(sin(asp1 div deg ))
setcell minof
temp = w225 * shade225 + w270 * shade270 + w315 *
shade315 + w360 * shade360
shade4 = int(temp div 2)

I tried porting in GRASS shell:

g.region rast=$GIS_OPT_input -p

# Compute hillshade at azimuth 225, 270, 315 and 360 at 30 degrees sun
illumination angle
r.shaded.relief map=$GIS_OPT_input shadedmap=shade225 altitude=30
azimuth=225 zmult=1 scale=1
r.shaded.relief map=$GIS_OPT_input shadedmap=shade270 altitude=30
azimuth=270 zmult=1 scale=1
r.shaded.relief map=$GIS_OPT_input shadedmap=shade315 altitude=30
azimuth=315 zmult=1 scale=1
r.shaded.relief map=$GIS_OPT_input shadedmap=shade360 altitude=30
azimuth=360 zmult=1 scale=1

# Resample DEM map
r.neighbors in=$GIS_OPT_input out=res1 method=average size=$GIS_OPT_window

# compute aspect map
r.slope.aspect elevation=res1 format=degrees prec=float aspect=aspect
zfactor=1.0 min_slp_allowed=0.0

# compute weights 225, 270, 315 and 360
r.mapcalc w225 = sin((aspect - 225)*sin(aspect - 225))
r.mapcalc w270 = sin((aspect - 270)*sin(aspect - 270))
r.mapcalc w315 = sin((aspect - 315)*sin(aspect - 315))
r.mapcalc w360 = sin((aspect)*sin(aspect))


#compute weighted
r.mapcalc weightedshade = (((w225 * shade225) + (w270 * shade270) +
(w315 * shade315) + (w360 * shade360))/2)
r.colors map=weightedshade color=grey

But I get not so good results, here's the output using nc_sp_08's
elev_ned_30m (X0=MDOW, X1= default r.shaded.relief module)
http://farm5.static.flickr.com/4148/5037622535_113924312b_b.jpg

Anything wrong with my commands?

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] multi-date image segmentation in grass

2010-09-28 Thread maning sambale
Hi,

I have several Landsat (4-7) images from 1990-2002  which will be used
for change analysis of land cover.  FAO's FRA-RSS [1] developed a
standard methodology for integrating multi-temporal Landsat into a
consistent sample tiles. The article [2] outlined the general method
as:
 - Image subsetting
 - Radiometric normalization
 - Contrast enhancement
 - Multi-date image segmentation

Most of the processing procedures I was able to do in GRASS GIS,
except for the last step.  The process combined the multi-date imagery
where an automated image segmentation was implemented to derive
preliminary polygons to identify changes.  I believe they used
ecognition for the segmentation.  I've used i.smap for preliminary
segmentation on separate dates but not for a combined multi-date
imagery.  Anyone tried this with GRASS? Can i.pca do this or another
tool?

Any idea

[1] http://www.fao.org/forestry/fra/remotesensingsurvey/en/
[2] 
http://www.informaworld.com/smpp/content~content=a923063789~db=all~jumptype=rss

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] [DWE]- Landsat Cloud detection with GRASS

2010-09-28 Thread maning sambale
2010/9/29 Markus Neteler nete...@osgeo.org:
 2010/9/28 António Rocha antonio.ro...@deimos.com.pt:
 Greetings

 I'm trying to find if GRASS includes some-sort of way to detect/eliminatd
 cloud cover and shadows from LANDSAT Images.
 The only thing I found was this:
 http://svn.osgeo.org/grass/grass-addons/imagery/i.landsat.acca/description.html
 But:
 1- has anyone tried/used this? Is this useful and working?

Yes, I've used it in the past using 6.3 svn.  I highly recommend it.
You need to convert your Landsat data to top-of-atmosphere radiance or
top-of-atmosphere reflectance first before running i.landsat.acca:
http://grass.osgeo.org/wiki/GRASS_AddOns#i.landsat.toar

Off-topic:  I hope the landsat tools (i.topo.corr, i.landsat.toar,
i.landsat.acca, i.landsat.dehaze, i.attcorr) gets into the core in the
future.


 It is currently under active development. Useful - hope so :)

 2- can I use it to detect-and-eliminate clouds?

 Should be the idea. The best way is to contact the author who
 is very responsive (see manual page).

 Cheers
 Markus

 Thanks
 Antonio


 __ Information from ESET NOD32 Antivirus, version of virus signature
 database 5486 (20100928) __

 The message was checked by ESET NOD32 Antivirus.

 http://www.eset.com


 ___
 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




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Re: Convert FCELL to CELL

2010-06-30 Thread maning sambale
r.mapcalc new_map=round(original_map)

On 6/30/10, Pedro Venâncio pedrongvenan...@yahoo.com wrote:
 Hello Micha and Eloi,

 This
 procedure only extracts the integer part of the number, it is not rounded.
 For example, 3.8
  becomes 3 instead of 4.

 Is there any way
  to force the rounding?

 Thanks!

 Pedro Venâncio





   Wouldn't:

 r.mapcalc new_map=int(original_map)
 do that?

   Hi,



 How can I convert a raster type FCELL into a raster type CELL?











-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] calculating road density

2010-06-02 Thread maning sambale
Thanks!  I'm trying this route.
On Wed, May 26, 2010 at 11:18 PM, Dylan Beaudette
dylan.beaude...@gmail.com wrote:
 v.rast.stats,

So far its running for 5 days on my whole region!  No segfault error so far. :)


-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] calculating road density

2010-05-26 Thread maning sambale
Hi,

I have two vectors in GRASS:
road network - line
grid - polygon

I want to compute road density and append the attributes to the grid poly.
In QGIS I tried using the vector tools (sum line length) with
shapefiles and it works for small areas bu not for my full region.
Any ideas?

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] r.coin behaving differently

2010-04-14 Thread maning sambale
Hi,

I am using r.coin for two raster maps

map1's range is
Range of data:min = 111  max = 140

map2' range
Range of data:min =   max = 17590


I used to get this result with r.coin
|   cat# |   111 |   112 |   120 |   130 |   140 |
  w cat 0 | w/o cat 0 |
||
|p  |  0.00 |  0.00 |  0.00 |  0.00 |
103.33 |103.33 |103.33 |
|h 1280 |   8376.26 |  51848.35 |395.22 |   5594.52 |
36.80 |  66251.15 |  66251.15 |
|i 1290 |  0.00 |  10521.99 |   3791.04 |636.42 |
162.64 |  15112.09 |  15112.09 |
|l 1330 |  0.00 | 66.13 |577.76 |  0.00 |
141.35 |785.24 |785.24 |
|p 1550 |   2079.69 |  12411.01 |440.71 |   1594.05 |
696.53 |  17221.99 |  17221.99 |
|r 2090 | 51.59 |   3862.45 |  0.00 |  0.00 |
0.00 |   3914.03 |   3914.03 |
|o 2150 | 112366.40 | 197525.87 |  0.00 |  17049.36 |
5081.42 | 332023.05 | 332023.05 |

Now I get a different cat# as row headers
|   cat# |   111 |   112 |   120 |   130 |   140 |
  w cat 0 | w/o cat 0 |
||
|p 1 |  0.00 |  0.00 |  0.00 |  0.00 |103.33 |
   103.33 |103.33 |
|h 3 |   8376.26 |  62370.34 |   4186.26 |   6230.93 |199.44 |
 81363.24 |  81363.24 |
|i 4 |  0.00 | 66.13 |577.76 |  0.00 |141.35 |
   785.24 |785.24 |
|l 7 |   2079.69 |  12411.01 |440.71 |   1594.05 |696.53 |
 17221.99 |  17221.99 |

Any idea why?  I'm running GRASS 6.4.0RC5 from kingchaos mac osx

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: r.coin behaving differently

2010-04-14 Thread maning sambale
Follow-up

On Wed, Apr 14, 2010 at 5:41 PM, maning sambale
emmanuel.samb...@gmail.com wrote:
 map2' range
 Range of data:    min =   max = 17590

The range above is not are ordinal categories, they are simply polygon
codes of administrative regions.

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: r.coin behaving differently

2010-04-14 Thread maning sambale
Never mind, map2 was an Data Type: DCELL, simply convert them in to
CELL and its OK. Doh!

On Wed, Apr 14, 2010 at 5:44 PM, maning sambale
emmanuel.samb...@gmail.com wrote:
 Follow-up

 On Wed, Apr 14, 2010 at 5:41 PM, maning sambale
 emmanuel.samb...@gmail.com wrote:
 map2' range
 Range of data:    min =   max = 17590

 The range above is not are ordinal categories, they are simply polygon
 codes of administrative regions.

 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] claslite modules in GRASS?

2010-04-12 Thread maning sambale
Follow-up inquiry, how do grass use spectral libraries in image
classification?  Google couldn't give any result except for Markus's
paper related to i.spec.unmix

On Fri, Apr 9, 2010 at 12:05 AM, Markus Neteler nete...@osgeo.org wrote:
 i.spec.unmix



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] claslite modules in GRASS?

2010-04-06 Thread maning sambale
claslite modules in GRASS?

I'm wondering if any GRASS users have been able to adapt a similar
workflow of Claslite within GRASS GIS.

http://claslite.ciw.edu/

Claslite is a suite of remote sensing modules developed by Carnegie
Institution for Science for rapid assessment of forest cover using
various satellite platforms.  Unfortunately, the software is not FOSS
but is it possible to adapt the developed tropical forest spectral
library within GRASS?
-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] r.le or r.li?

2010-03-28 Thread maning sambale
 That's right.
 The differences: r.le is the first implementation but unmaintained.
 r.li is a new implementation and somewhat maintained (just fixed
 a bug today). Hopefully more indices are added in future.
Thanks for clarifying. So r.li it is.

 Follow up inquiry, how do I include the the patch stats to a vector
 polygon?  For example, I have a watershed polygon and landcover
 raster.  I want to include as an attribute, the largest interior patch
 within each sub-basin.

 ... sorry, no idea. In any case think raster for this module family...

That's my line of thinking at the moment.



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] r.le or r.li?

2010-03-14 Thread maning sambale
Hi,

I'm doing some patch stats analysis (largest patch, edge and interior
patches).   There are meta modules available in grass r.le and r.li.
The seem to provide similar set of patch statistics calculations.  I
want to know the basic difference of the two.  r.li seems to be very
similar to fragstats or is it?

Follow up inquiry, how do I include the the patch stats to a vector
polygon?  For example, I have a watershed polygon and landcover
raster.  I want to include as an attribute, the largest interior patch
within each sub-basin.


Thanks in advance!
-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: weird results with mapcalc on different region extents

2010-03-11 Thread maning sambale
same issue with grass6.4 in Ubuntu,
runs OK with small region but not in a bigger region


On Thu, Mar 11, 2010 at 3:05 PM, maning sambale
emmanuel.samb...@gmail.com wrote:
 Hi,

 I have this simple mapcalc run that determines the proportion of
 landcover type over the land, given a moving window.

 A = binary map 1=with cover pixels, 0=w/o cover pixels
 B = binary map 1=land pixels 0=non-land pixels
 WINDOW = 5

 r.neighbors input=A output=C method=sum size=$WINDOW --o
 r.neighbors input=B output=D method=sum size=$WINDOW --o

 r.mapcalc  EOF
 E = 1.0 * C
 F = 1.0 * D
 pf = (E/D)
 EOF

 pf should result in a percentage of cover pixels over the land.

 If I run the data in a small region and the results are OK.
 r.info -r pf
 min=0
 max=1

 r.stats -acp pf
  100%
 0-0.003922 83506029.482130 102816  26.75%
 0.039216-0.043137 1436762.431469 1769   0.46%
 0.047059-0.05098 1624.378102 2   0.00%
 0.07-0.070588 2436.567153 3   0.00%
 0.078431-0.082353 916149.249687 1128   0.29%
 0.098039-0.101961 812.189051 1   0.00%
 0.117647-0.121569 987621.886188 1216   0.32%
 0.14902-0.152941 812.189051 1   0.00%
 0.156863-0.160784 724472.633618 892   0.23%
 0.2-0.203922 990058.453341 1219   0.32%
 0.239216-0.243137 842240.046034 1037   0.27%
 0.247059-0.25098 1624.378102 2   0.00%
 0.27-0.270588 2436.567153 3   0.00%
 0.278431-0.282353 616451.489816 759   0.20%
 0.298039-0.301961 1624.378102 2   0.00%
 0.317647-0.321569 604268.654049 744   0.19%
 0.33-0.337255 1624.378102 2   0.00%
 0.356863-0.360784 660309.698578 813   0.21%
 0.4-0.403922 1019297.259182 1255   0.33%
 0.439216-0.443137 621324.624123 765   0.20%
 0.447059-0.45098 1624.378102 2   0.00%
 0.47-0.470588 1624.378102 2   0.00%
 0.478431-0.482353 614827.111714 757   0.20%
 0.498039-0.501961 1624.378102 2   0.00%
 0.517647-0.521569 594522.385435 732   0.19%
 0.53-0.537255 1624.378102 2   0.00%
 0.556863-0.560784 645690.295657 795   0.21%
 0.6-0.603922 1009550.990569 1243   0.32%
 0.639216-0.643137 631883.081788 778   0.20%
 0.647059-0.65098 2436.567153 3   0.00%
 0.67-0.670588 1624.378102 2   0.00%
 0.678431-0.682353 644878.106606 794   0.21%
 0.698039-0.701961 3248.756205 4   0.00%
 0.717647-0.721569 605080.843100 745   0.19%
 0.73-0.737255 1624.378102 2   0.00%
 0.74902-0.752941 812.189051 1   0.00%
 0.756863-0.760784 808940.294937 996   0.26%
 0.8-0.803922 1011175.368671 1245   0.32%
 0.839216-0.843137 745589.548948 918   0.24%
 0.847059-0.85098 3248.756205 4   0.00%
 0.87-0.870588 4873.134307 6   0.00%
 0.878431-0.882353 998180.343853 1229   0.32%
 0.898039-0.901961 1624.378102 2   0.00%
 0.917647-0.921569 933205.219761 1149   0.30%
 0.93-0.937255 8121.890511 10   0.00%
 0.94902-0.952941 8934.079563 11   0.00%
 0.956863-0.960784 1423767.406651 1753   0.46%
 0.996078-1 93443162.522784 115051  29.94%
 * 115074193.521699 141684  36.87%

 But, when I run the data in the full map region, I get a different result

 r.info -r pf
 min=0
 max=25
 r.stats -acp pf
  100%
 0-0.098039 1593170550.199143 1961576  39.19%
 0.098039-0.196078 18625931.509825 22933   0.46%
 0.196078-0.294118 26394519.783980 32498   0.65%
 0.294118-0.392157 13855133.023422 17059   0.34%
 0.392157-0.490196 24355925.265619 29988   0.60%
 0.490196-0.588235 13444165.363544 16553   0.33%
 0.588235-0.686275 24779075.761263 30509   0.61%
 0.686275-0.784314 15370677.792853 18925   0.38%
 0.784314-0.882353 29887744.892934 36799   0.74%
 0.882353-0.980392 25306186.455453 31158   0.62%
 0.980392-1.078431 1858477789.081065 2288233  45.71%
 1.078431-1.176471 77970.148910 96   0.00%
 1.176471-1.27451 51980.099273 64   0.00%
 1.27451-1.372549 25177.860585 31   0.00%
 1.372549-1.470588 21116.915330 26   0.00%
 1.470588-1.568627 28426.616790 35   0.00%
 1.568627-1.67 5685.323358 7   0.00%
 1.67-1.764706 14619.402921 18   0.00%
 1.764706-1.862745 17055.970074 21   0.00%
 1.862745-1.960784 14619.402921 18   0.00%
 1.960784-2.058824 4873.134307 6   0.00%
 2.058824-2.156863 14619.402921 18   0.00%
 2.156863-2.254902 5685.323358 7   0.00%
 2.254902-2.352941 11370.646716 14   0.00%
 2.352941-2.45098 1624.378102 2   0.00%
 2.45098-2.54902 11370.646716 14   0.00%
 2.54902-2.647059 3248.756205 4   0.00%
 2.647059-2.745098 2436.567153 3   0.00%
 2.745098-2.843137 12995.024818 16   0.00%
 2.843137-2.941176 812.189051 1   0.00%
 2.941176-3.039216 7309.701460 9   0.00%
 3.039216-3.137255 10558.457665 13   0.00%
 3.137255-3.235294 812.189051 1   0.00%
 3.235294-3.33 812.189051 1   0.00%
 3.33-3.431373 4060.945256 5   0.00%
 3.431373-3.529412 812.189051 1   0.00%
 3.529412-3.627451 8934.079563 11   0.00%
 3.627451-3.72549 1624.378102 2   0.00%
 3.72549-3.823529 1624.378102 2   0.00%
 3.823529-3.921569 812.189051 1   0.00%
 3.921569-4.019608 3248.756205 4   0.00%
 4.117647-4.215686 10558.457665 13   0.00%
 4.313725-4.411765 812.189051 1   0.00%
 4.411765-4.509804 812.189051 1   0.00%
 4.509804-4.607843 812.189051 1   0.00%
 4.607843-4.705882 812.189051 1   0.00%
 5-5.098039

[GRASS-user] weird results with mapcalc on different region extents

2010-03-10 Thread maning sambale
Hi,

I have this simple mapcalc run that determines the proportion of
landcover type over the land, given a moving window.

A = binary map 1=with cover pixels, 0=w/o cover pixels
B = binary map 1=land pixels 0=non-land pixels
WINDOW = 5

r.neighbors input=A output=C method=sum size=$WINDOW --o
r.neighbors input=B output=D method=sum size=$WINDOW --o

r.mapcalc  EOF
E = 1.0 * C
F = 1.0 * D
pf = (E/D)
EOF

pf should result in a percentage of cover pixels over the land.

If I run the data in a small region and the results are OK.
r.info -r pf
min=0
max=1

r.stats -acp pf
 100%
0-0.003922 83506029.482130 102816  26.75%
0.039216-0.043137 1436762.431469 1769   0.46%
0.047059-0.05098 1624.378102 2   0.00%
0.07-0.070588 2436.567153 3   0.00%
0.078431-0.082353 916149.249687 1128   0.29%
0.098039-0.101961 812.189051 1   0.00%
0.117647-0.121569 987621.886188 1216   0.32%
0.14902-0.152941 812.189051 1   0.00%
0.156863-0.160784 724472.633618 892   0.23%
0.2-0.203922 990058.453341 1219   0.32%
0.239216-0.243137 842240.046034 1037   0.27%
0.247059-0.25098 1624.378102 2   0.00%
0.27-0.270588 2436.567153 3   0.00%
0.278431-0.282353 616451.489816 759   0.20%
0.298039-0.301961 1624.378102 2   0.00%
0.317647-0.321569 604268.654049 744   0.19%
0.33-0.337255 1624.378102 2   0.00%
0.356863-0.360784 660309.698578 813   0.21%
0.4-0.403922 1019297.259182 1255   0.33%
0.439216-0.443137 621324.624123 765   0.20%
0.447059-0.45098 1624.378102 2   0.00%
0.47-0.470588 1624.378102 2   0.00%
0.478431-0.482353 614827.111714 757   0.20%
0.498039-0.501961 1624.378102 2   0.00%
0.517647-0.521569 594522.385435 732   0.19%
0.53-0.537255 1624.378102 2   0.00%
0.556863-0.560784 645690.295657 795   0.21%
0.6-0.603922 1009550.990569 1243   0.32%
0.639216-0.643137 631883.081788 778   0.20%
0.647059-0.65098 2436.567153 3   0.00%
0.67-0.670588 1624.378102 2   0.00%
0.678431-0.682353 644878.106606 794   0.21%
0.698039-0.701961 3248.756205 4   0.00%
0.717647-0.721569 605080.843100 745   0.19%
0.73-0.737255 1624.378102 2   0.00%
0.74902-0.752941 812.189051 1   0.00%
0.756863-0.760784 808940.294937 996   0.26%
0.8-0.803922 1011175.368671 1245   0.32%
0.839216-0.843137 745589.548948 918   0.24%
0.847059-0.85098 3248.756205 4   0.00%
0.87-0.870588 4873.134307 6   0.00%
0.878431-0.882353 998180.343853 1229   0.32%
0.898039-0.901961 1624.378102 2   0.00%
0.917647-0.921569 933205.219761 1149   0.30%
0.93-0.937255 8121.890511 10   0.00%
0.94902-0.952941 8934.079563 11   0.00%
0.956863-0.960784 1423767.406651 1753   0.46%
0.996078-1 93443162.522784 115051  29.94%
* 115074193.521699 141684  36.87%

But, when I run the data in the full map region, I get a different result

r.info -r pf
min=0
max=25
r.stats -acp pf
 100%
0-0.098039 1593170550.199143 1961576  39.19%
0.098039-0.196078 18625931.509825 22933   0.46%
0.196078-0.294118 26394519.783980 32498   0.65%
0.294118-0.392157 13855133.023422 17059   0.34%
0.392157-0.490196 24355925.265619 29988   0.60%
0.490196-0.588235 13444165.363544 16553   0.33%
0.588235-0.686275 24779075.761263 30509   0.61%
0.686275-0.784314 15370677.792853 18925   0.38%
0.784314-0.882353 29887744.892934 36799   0.74%
0.882353-0.980392 25306186.455453 31158   0.62%
0.980392-1.078431 1858477789.081065 2288233  45.71%
1.078431-1.176471 77970.148910 96   0.00%
1.176471-1.27451 51980.099273 64   0.00%
1.27451-1.372549 25177.860585 31   0.00%
1.372549-1.470588 21116.915330 26   0.00%
1.470588-1.568627 28426.616790 35   0.00%
1.568627-1.67 5685.323358 7   0.00%
1.67-1.764706 14619.402921 18   0.00%
1.764706-1.862745 17055.970074 21   0.00%
1.862745-1.960784 14619.402921 18   0.00%
1.960784-2.058824 4873.134307 6   0.00%
2.058824-2.156863 14619.402921 18   0.00%
2.156863-2.254902 5685.323358 7   0.00%
2.254902-2.352941 11370.646716 14   0.00%
2.352941-2.45098 1624.378102 2   0.00%
2.45098-2.54902 11370.646716 14   0.00%
2.54902-2.647059 3248.756205 4   0.00%
2.647059-2.745098 2436.567153 3   0.00%
2.745098-2.843137 12995.024818 16   0.00%
2.843137-2.941176 812.189051 1   0.00%
2.941176-3.039216 7309.701460 9   0.00%
3.039216-3.137255 10558.457665 13   0.00%
3.137255-3.235294 812.189051 1   0.00%
3.235294-3.33 812.189051 1   0.00%
3.33-3.431373 4060.945256 5   0.00%
3.431373-3.529412 812.189051 1   0.00%
3.529412-3.627451 8934.079563 11   0.00%
3.627451-3.72549 1624.378102 2   0.00%
3.72549-3.823529 1624.378102 2   0.00%
3.823529-3.921569 812.189051 1   0.00%
3.921569-4.019608 3248.756205 4   0.00%
4.117647-4.215686 10558.457665 13   0.00%
4.313725-4.411765 812.189051 1   0.00%
4.411765-4.509804 812.189051 1   0.00%
4.509804-4.607843 812.189051 1   0.00%
4.607843-4.705882 812.189051 1   0.00%
5-5.098039 20304.726279 25   0.00%
5.196078-5.294118 812.189051 1   0.00%
5.294118-5.392157 2436.567153 3   0.00%
5.980392-6.078431 5685.323358 7   0.00%
6.176471-6.27451 12995.024818 16   0.00%
6.27451-6.372549 812.189051 1   0.00%
6.67-6.764706 812.189051 1   0.00%
6.960784-7.058824 1624.378102 2   0.00%

Re: [GRASS-user] reproject from lat/long to utm

2009-12-01 Thread maning sambale
Try to include the full path of your location ([/work/projects/Willamette])

On Wed, Dec 2, 2009 at 1:47 PM, Janet Choate jsc@gmail.com wrote:
 HI,
 thanx for any help.  seems this should be a pretty straight forward thing to
 do, but am not having any success.
 i am trying to reproject a raster from lat/long to utm.
 i created a new location with a utm projection, then tried to use r.proj to
 reproject a raster from another location.

 r.proj input=dem location=/work/projects/Willamette mapset=PERMANENT

 however, i get the following error:
 Mapset [PERMANENT] in input location [/work/projects/Willamette] - not found

 thanx for any assistance,
 Janet

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





-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] reproject from lat/long to utm

2009-12-01 Thread maning sambale
On Wed, Dec 2, 2009 at 2:28 PM, Janet Choate jsc@gmail.com wrote:
 I tried using the full path (i.e. the PERMANENT mapset lives within the
 /work/projects/Willamette location), but get the 'not found' error.
 is there any other command i need to run prior to r.proj, or any additional
 options to include, or how to list the path other than the way i am?

just to check from your lat/long mapset, run g.gisenv
It should show you the location, mapset and GRASS database.

This is what you use in r.proj parameters

 the original data i imported was a dem from usgs national elevation dataset
 (ned) seamless data distribution system (sdds).  i had to export it from ned
 as a geotiff to bring it into grass.  i tried exporting the data from ned as
 an arcgrid and importing into grass with r.in.gdal, but was unable to bring
 it into grass.

 so if anyone can tell me the steps to bringing ned seamless data into grass,
 then reprojecting to utm, i would be grateful.  or at this point, just how
 to use r.proj to reproject to utm.
 thanx, janet

 On Tue, Dec 1, 2009 at 10:14 PM, maning sambale emmanuel.samb...@gmail.com
 wrote:

 Try to include the full path of your location
 ([/work/projects/Willamette])

 On Wed, Dec 2, 2009 at 1:47 PM, Janet Choate jsc@gmail.com wrote:
  HI,
  thanx for any help.  seems this should be a pretty straight forward
  thing to
  do, but am not having any success.
  i am trying to reproject a raster from lat/long to utm.
  i created a new location with a utm projection, then tried to use r.proj
  to
  reproject a raster from another location.
 
  r.proj input=dem location=/work/projects/Willamette mapset=PERMANENT
 
  however, i get the following error:
  Mapset [PERMANENT] in input location [/work/projects/Willamette] - not
  found
 
  thanx for any assistance,
  Janet
 
  ___
  grass-user mailing list
  grass-user@lists.osgeo.org
  http://lists.osgeo.org/mailman/listinfo/grass-user
 
 



 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --





-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] reproject from lat/long to utm

2009-12-01 Thread maning sambale
try this:

r.proj input=dem dbase=/work/projects/ location=Willamette mapset=PERMANENT

On Wed, Dec 2, 2009 at 2:45 PM, maning sambale
emmanuel.samb...@gmail.com wrote:
 On Wed, Dec 2, 2009 at 2:28 PM, Janet Choate jsc@gmail.com wrote:
 I tried using the full path (i.e. the PERMANENT mapset lives within the
 /work/projects/Willamette location), but get the 'not found' error.
 is there any other command i need to run prior to r.proj, or any additional
 options to include, or how to list the path other than the way i am?

 just to check from your lat/long mapset, run g.gisenv
 It should show you the location, mapset and GRASS database.

 This is what you use in r.proj parameters

 the original data i imported was a dem from usgs national elevation dataset
 (ned) seamless data distribution system (sdds).  i had to export it from ned
 as a geotiff to bring it into grass.  i tried exporting the data from ned as
 an arcgrid and importing into grass with r.in.gdal, but was unable to bring
 it into grass.

 so if anyone can tell me the steps to bringing ned seamless data into grass,
 then reprojecting to utm, i would be grateful.  or at this point, just how
 to use r.proj to reproject to utm.
 thanx, janet

 On Tue, Dec 1, 2009 at 10:14 PM, maning sambale emmanuel.samb...@gmail.com
 wrote:

 Try to include the full path of your location
 ([/work/projects/Willamette])

 On Wed, Dec 2, 2009 at 1:47 PM, Janet Choate jsc@gmail.com wrote:
  HI,
  thanx for any help.  seems this should be a pretty straight forward
  thing to
  do, but am not having any success.
  i am trying to reproject a raster from lat/long to utm.
  i created a new location with a utm projection, then tried to use r.proj
  to
  reproject a raster from another location.
 
  r.proj input=dem location=/work/projects/Willamette mapset=PERMANENT
 
  however, i get the following error:
  Mapset [PERMANENT] in input location [/work/projects/Willamette] - not
  found
 
  thanx for any assistance,
  Janet
 
  ___
  grass-user mailing list
  grass-user@lists.osgeo.org
  http://lists.osgeo.org/mailman/listinfo/grass-user
 
 



 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --





 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] internal buffer of polygon

2009-10-28 Thread maning sambale
I tried using a negative values as
buffer=-3000 in v.buffer but it doesn't give me internal buffer in a
polygon layer.

Any hints?

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Re: v.out.ogr kml does not add name field

2009-10-15 Thread maning sambale
It works using ogr2ogr.

On Thu, Oct 15, 2009 at 12:45 AM, Hermann Peifer pei...@gmx.eu wrote:
 Works fine for me with Grass 6.4 from SVN.

 Which version are you using?
old grass6.4svn

 Hermann


 maning sambale wrote:

 Hi,

 I have a grass vector of provincial boundaries
 db.describe -c generalized_adminncols: 3

 nrows: 415
 Column 1: cat:INTEGER:11
 Column 2: region:CHARACTER:15
 Column 3: province:CHARACTER:50

 Now, I exported the vector to KML polygon using the province
 attributes as the KML name element.

 v.out.ogr in=generalized_admin format=KML type=area dsn=test_kml.kml
 layer=1 dsco=NameField=province
 Exporting 220 areas (may take some time)...
  100%
 227 features written
 WARNING: 7 features without category written

 The output KML does not include the province names instead the cat
 values are in the KML name element.

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




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] v.out.ogr kml does not add name field

2009-10-14 Thread maning sambale
Hi,

I have a grass vector of provincial boundaries
db.describe -c generalized_adminncols: 3

nrows: 415
Column 1: cat:INTEGER:11
Column 2: region:CHARACTER:15
Column 3: province:CHARACTER:50

Now, I exported the vector to KML polygon using the province
attributes as the KML name element.

v.out.ogr in=generalized_admin format=KML type=area dsn=test_kml.kml
layer=1 dsco=NameField=province
Exporting 220 areas (may take some time)...
 100%
227 features written
WARNING: 7 features without category written

The output KML does not include the province names instead the cat
values are in the KML name element.
-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] v.edit does not copy attribute?

2009-09-25 Thread maning sambale
I tried this:

v.edit map=map2 bgmap=map1 tool=copy bbox=188716,1353547,341637,1499539 -e --o

but the attributes table from map1 are not copied.  Tried adding
attributes  with v.distance an v.category

v.extract doesn't have the bbox feature

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] How to get perimeter?

2009-09-17 Thread maning sambale
v.report map=name option=length units=meters

On 9/17/09, Margherita Di Leo direg...@gmail.com wrote:
 Hi List,

 I'm searching for a way to delineate the perimeter of an area in a
 raster map. The area is like the result of r.clump, so it is
 caracterized by a integer number in the raster map. I have no problem in
 trasforming the raster map in vector with r.to.vect, but the vector is
 an area, not a line. I need to know the lenght of the contour of the
 areas. What is the best way to obtain that?

 Thank you in advance

 Margherita

 --
 Eng. Margherita Di Leo
 Ph.D. Student
 Methods and Technologies for Environmental Monitoring
 Department of Environmental Engineering and Physics (DIFA)

 University of Basilicata
 Campus Macchia Romana
 85100 - Potenza
 Italy

 Office: +39-0971205363
 Fax: +39-0971205160
 E-mail: dileomargherita AT gmail DOT com
 Skype: dileomargherita
 URL:
 http://www.difa.unibas.it/A_Manager_PP.do?azione=visualizzaHomePageid=106





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



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Dsn error in Window$

2009-09-16 Thread maning sambale
afaik,

it should be:
 v.out.ogr -e input=hidr...@permanent dsn=C:/tmp/hidro3D

On 9/17/09, Vincent Blanqué vincent.blan...@oriondata.cl wrote:
 Hi,

 It's my first post here. I am quite happy to join the powerful world of
 GRASS. My first good impression is with the 'v.drape' process.

 With the 'v.out.ogr' command, I get an issue with the dsn parameter.
 I tried various syntax for dsn path  filename but none works fine.

 v.out.ogr -e --verbose input=hidr...@permanent dsn=C:/tmp/hidro3D.shp


 v.out.ogr -e --verbose input=hidr...@permanent dsn=C:/tmp/hidro3D.shp

 ERROR 1: Failed to create directory C:/tmp/hidro3D.shp
 for shapefile datastore.

 v.out.ogr -e input=hidr...@permanent dsn=/c/tmp/hidro3D.shp

 I am using GRASS 6.4.0svn (2009) - Revision: 37101 in Windows Vista
 (very bad choice).

 Do you know what happens ?

 In Linux Ubuntu, everything is great.

 see U,

 Vincent






-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] creating map series purely on ps.map

2009-08-20 Thread maning sambale
Hi,

Basically my workflow for creating paper maps is create the main map
in ps.map.  Export the ps output to pdf or png in a graphics software
(scribus or illustrator).  Edit all map elements legend, anotations,
etc.

What I want is to be bale to create the map purely within grass and
ps.map module.

Several questions:
1. I need to prepare separate mapsheets of landcover data.  I have a
vector layer in grass of all the boundingboxes.  How can I call ps.map
to create mapsheets using the individual boundingboxes as the region
settings?

2. Each mapsheet requires a different locatormap as reference.  How do
I add the locatormap using a bigger region settings around the mapped
region settings?

3. Legend symbology for different landcover categories within a single
polygon vector.  As a workaround, I separate landcover categories as
separate vector layers.

4. Tips and tricks for different annotations.

5. Combine different ps.map output into one map.

Any advice? Thanks in advance!
-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] creating map series purely on ps.map

2009-08-20 Thread maning sambale
Thanks hamish!  I will try you suggestions back at work.  Most of my
queries where answered.
One of my question, I discovered that when I open the .ps file to
Adobe Illustrator.  It seems that the whole vector layer is included
with the a white foreground mask outside the region settings.

Am I correct that the vector is not clipped within the region settings?

On 8/20/09, Hamish hamis...@yahoo.com wrote:
 maning wrote:
 Basically my workflow for creating paper maps is create the
 main map in ps.map.  Export the ps output to pdf or png in a
 graphics software (scribus or illustrator).  Edit all map elements
 legend, anotations, etc.

 What I want is to be bale to create the map purely within
 grass and ps.map module.

 sure, it's possible.

 Several questions:
 1. I need to prepare separate mapsheets of landcover data.  I have a
 vector layer in grass of all the boundingboxes.  How can I call ps.map
 to create mapsheets using the individual boundingboxes as the region
 settings?

 hmmm how to extract the bounding box of each category in a vector map?
 I think there might be a way, but I'm not sure, maybe someone else knows.

 anyway you can write a shell loop to use v.extract to pull out that
 one category and then use g.region vect=extract_cat_$CAT to zoom in
 on that one area. If you want a buffer around it, try like:

 g.region vect=$MAP
 g.region res=5000 -a

 which will grow outwards until it hits a 5km round number.

 while you are there save the region with 'g.region save=' with a name
 derived from that box name so you can jump to it quickly later on.

 something like

 for CAT in `v.category master_bounding_box_map option=print` ; do
v.extract in=master_bounding_box_map out=bbox_$CAT list=$CAT
g.region -d
g.region vect=bbox_$CAT
g.region res=5000 -a
g.region res=10   # back to normal for raster map res
g.region save=bbox_$CAT
 done

 then

 for REGION in `g.mlist region pattern=bbox_*` ; do
g.region region=$REGION
ps.map out=landcover_$REGION in=-  EOF

end
 EOF
ps2pdf landcover_$REGION.ps
 done



 2. Each mapsheet requires a different locatormap as reference.  How
 do I add the locatormap using a bigger region settings around
 the mapped region settings?

 inlays:
 first for each mapsheet save a region with 'g.region save='.
 then (in a loop) create a ps.map .eps file drawing the overview
 map but also using the ps.map region instruction to draw a box
 over the top of it. I like to use a red frame with width of 2,
 but you'll have to experiment with the width as this .eps will be
 much smaller on the page once it is embedded. Because it will be
 so reduced you can probably set the resolution quite coarse, with
 a little math you can figure out how many raster column cells equate
 to 300dpi at your target width (e.g. if you want the inset to be 2
 wide set the region resolution to be so columns is about 600 cells.)

 once you have those just use the ps.map 'eps' instruction to place
 them on the map. Use the 'rectangle' instruction to add a thicker
 border if you like.

 again, if you write the script like this:

 for MAP in `g.mlist rast patter=*.landuse` ; do
case $MAP in
   landuse23)
  MAP_TITLE=Map of landuse category 23
  ;;
landuse24)
  MAP_TITLE= ...
  ;;
esac

ps.map out=$MAP.ps in=-  EOF
   raster $MAP
   eps 10% 10%
 epsfile locator_$MAP.eps
 scale 0.2
 end
   text 65 7365000 Landuse for $MAP_TITLE
 end
   end
 EOF

 then you can use shell variable replacement within the instruction file.
 another way is to have a shell script write a bunch of files instead of
 running ps.map directly; or just make a master instruction file then copy
 it a bunch of times and tweak each individually. that last way is a pain
 if you decide to change anything though, as you have to make the edit
 in every copy of the file.

 you can put the non changing bits in a template file and then do like:

 for MAP in `g.mlist` ; do
   ps.map out=test.ps in=-  EOF
  raster $MAP
  eps ..
epsfile region.$MAP
end
  `cat common_instructions.txt`
  end
 EOF


 3. Legend symbology for different landcover categories within a single
 polygon vector.  As a workaround, I separate landcover categories as
 separate vector layers.

 you can keep them all in the same file, just repeat the vareas command
 but limit by SQL query each time, and set the label for each vareas
 as needed. see example here:
  http://grass.osgeo.org/wiki/Psmap_fill_patterns
  http://grass.osgeo.org/wiki/AreaFillPatterns


 4. Tips and tricks for different annotations.

 you'd have to be more specific about what you mean. many tips here:
   http://grass.osgeo.org/wiki/Ps.map
   http://grass.osgeo.org/wiki/IconSymbols#North_Arrows

 see also ps.output in wiki addons.


 5. Combine different ps.map output into one map.

 see
  

Re: [GRASS-user] change some labels to upper/lower case in vector attributes

2009-08-18 Thread maning sambale
Thanks will try this one.

On Tue, Aug 18, 2009 at 7:36 PM, Moritz
Lennertmlenn...@club.worldonline.be wrote:
 On 17/08/09 08:19, maning sambale wrote:

 One column of my vector is the name of the town all in uppercase.
 What I want is to convert them, i.e.,

 SAN DIEGO  San Diego

 My grass dbase is sqlite

 Something like this should do the trick (untested):

 UPDATE YourTable set town=substr(town,1,1) || lower(substr(town,2))

 Moritz




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] change some labels to upper/lower case in vector attributes

2009-08-17 Thread maning sambale
One column of my vector is the name of the town all in uppercase.
What I want is to convert them, i.e.,

SAN DIEGO  San Diego

My grass dbase is sqlite
-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] updating vector attributes for separate islands

2009-08-17 Thread maning sambale
Thanks! it's working now

On 8/17/09, Markus Metz markus.metz.gisw...@googlemail.com wrote:

 maning sambale wrote:
 Hi,

 I imported a shapefile of provinces in grass.  The shapefile has
 around 84 provinces composed of different islands.  The attribute
 table contains only 84 rows corresponding to 84 provinces.  I used
 v.clean bpol tool to topologically clean the vectors.  Then used
 v.db.droptable and add a new  a table to add island names for each
 polygon.

 The new table created still contain 84 rows and when I edit an island
 centroid the island name is also added to the other polygon covered by
 the same province.

 They all have the same category number, so if you add a row to the
 attribute table, all polygons with that category number will use that
 row. You can set new categories for each polygon separately with
 v.category in a new layer 2, then add a table for layer 2 and edit the
 attribute table.

 Markus M



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] updating vector attributes for separate islands

2009-08-16 Thread maning sambale
Hi,

I imported a shapefile of provinces in grass.  The shapefile has
around 84 provinces composed of different islands.  The attribute
table contains only 84 rows corresponding to 84 provinces.  I used
v.clean bpol tool to topologically clean the vectors.  Then used
v.db.droptable and add a new  a table to add island names for each
polygon.

The new table created still contain 84 rows and when I edit an island
centroid the island name is also added to the other polygon covered by
the same province.

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] bus error in v.in.ogr

2009-08-14 Thread maning sambale
I am importing a polygon shapefile and I get this error.  The
shapefile was created using QGIS.

v.in.ogr layer=forest_p116r048_repair output=forest_p116r048
dsn=/Users/maning/pfua/map_data/forest_2002_digit/forest_p116r048_p117r047_edited
--o
WARNING: Vector map forest_p116r048 already exists and will be
 overwritten
Projection of input dataset and current location appear to match
WARNING: Vector map forest_p116r048 already exists and will be
 overwritten
Layer: forest_p116r048_repair
Importing map 181 features...
WARNING: Area size [7.5e-09], area not imported
-
Building topology for vector map forest_p116r048...
58952 vertices registered
230 primitives registered
Building areas:  100%
225 areas built
225 isles built
Attaching islands:  100%
Attaching centroids:  100%
Topology was built
Number of nodes :   228
Number of primitives:   230
Number of points:   0
Number of lines :   0
Number of boundaries:   230
Number of centroids :   0
Number of areas :   225
Number of isles :   225
Number of incorrect boundaries   :   5
Number of areas without centroid :   225
-
WARNING: Cleaning polygons, result is not guaranteed!
Building topology for vector map forest_p116r048...
Topology was built
Number of nodes :   228
Number of primitives:   230
Number of points:   0
Number of lines :   0
Number of boundaries:   230
Number of centroids :   0
Number of areas :   -
Number of isles :   -
-
Break polygons:
Registering points 46354
All points (vertices): 58943
Registered points (unique coordinates): 46354
Points marked for break:   345
Breaks:   610
-
Remove duplicates:
Duplicates:   230
-
Break boundaries:
Intersections:39
-
Remove duplicates:
Duplicates: 7
-
Clean boundaries at nodes:
Modifications: 0Bus error


-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] bus error in v.in.ogr

2009-08-14 Thread maning sambale
On Fri, Aug 14, 2009 at 9:39 PM, maning
sambaleemmanuel.samb...@gmail.com wrote:
 There was a vector file but not layers so in effect, there was an error.

 On Fri, Aug 14, 2009 at 9:21 PM, Milton Cezar
 Ribeiromiltinho.astrona...@gmail.com wrote:
 Hi Maning,

 O flyed on your email and saw only warnings, not errors.
 It appears to worked.

 bests
 milton

 2009/8/14 maning sambale emmanuel.samb...@gmail.com

 I am importing a polygon shapefile and I get this error.  The
 shapefile was created using QGIS.

 v.in.ogr layer=forest_p116r048_repair output=forest_p116r048

 dsn=/Users/maning/pfua/map_data/forest_2002_digit/forest_p116r048_p117r047_edited
 --o
 WARNING: Vector map forest_p116r048 already exists and will be
         overwritten
 Projection of input dataset and current location appear to match
 WARNING: Vector map forest_p116r048 already exists and will be
         overwritten
 Layer: forest_p116r048_repair
 Importing map 181 features...
 WARNING: Area size [7.5e-09], area not imported
 -
 Building topology for vector map forest_p116r048...
 58952 vertices registered
 230 primitives registered
 Building areas:  100%
 225 areas built
 225 isles built
 Attaching islands:  100%
 Attaching centroids:  100%
 Topology was built
 Number of nodes     :   228
 Number of primitives:   230
 Number of points    :   0
 Number of lines     :   0
 Number of boundaries:   230
 Number of centroids :   0
 Number of areas     :   225
 Number of isles     :   225
 Number of incorrect boundaries   :   5
 Number of areas without centroid :   225
 -
 WARNING: Cleaning polygons, result is not guaranteed!
 Building topology for vector map forest_p116r048...
 Topology was built
 Number of nodes     :   228
 Number of primitives:   230
 Number of points    :   0
 Number of lines     :   0
 Number of boundaries:   230
 Number of centroids :   0
 Number of areas     :   -
 Number of isles     :   -
 -
 Break polygons:
 Registering points 46354
 All points (vertices): 58943
 Registered points (unique coordinates): 46354
 Points marked for break:   345
 Breaks:   610
 -
 Remove duplicates:
 Duplicates:   230
 -
 Break boundaries:
 Intersections:    39
 -
 Remove duplicates:
 Duplicates:     7
 -
 Clean boundaries at nodes:
 Modifications:     0Bus error


 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user





 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] wishlist to add a grass module for Landsat ETM+ SLC-off Gapfill

2009-06-22 Thread maning sambale
Landsat ETM+ SLC-off Gapfill
http://l7gapfill.sourceforge.net/

L7gapfill uses a multi-scale segment model to guide interpolation of
spectral data across gaps in Landsat 7 SLC-off images. Gap pixels are
filled with concurrent spectral data, essential for applications that
require same-day spectral information.

Would be nice to add to the expanding grass LANDSAT modules
-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] calculate speed fron GPX in GRASS

2009-06-10 Thread maning sambale
I want to calculate the speed of my GPS trace in grass.  Converting
the vector points where an attribute speed is added.
Any ideas how to do it?

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] error in ps.map scalebar units argument

2009-05-26 Thread maning sambale
Ahh no units option in my grass version:

/250k_layout.map
+ REGION=surigao_test_reg
+ TITLE=SURIGAO
+ PCGS=2535
+ g.region region=surigao_test_reg res=100
+ ps.map out=surigao_test_reg.ps
where  x y
length  length
height  height
segment no_segemnts
numbers no_labels
fontsize   fontsize
background [Y|n]
enter 'end' when done, 'exit' to quit
+ echo 'finished generating postscript map of SURIGAO'
finished generating postscript map of SURIGAO
[Raster MASK present]


I'll try a new grass version then.


On Tue, May 26, 2009 at 7:35 PM, Hamish hamis...@yahoo.com wrote:

 maning wrote:
 GRASS 6.4.svn
 
 ERROR: units kilometers : illegal request (scalebar)


 that indicates it doesn't understand units (as opposed to km,miles,etc
 qualifier being bad).

 I don't know why it doesn't work for you.


 can you try putting in your script:

  scalebar s
   help
   exit
   end
  end

 and see if units   auto|meters|kilometers|feet|miles|nautmiles is
 there in the output?


 All I can suggest is to 'make clean', 'svn update' to the latest 6.4svn
 (checking for 'C'onflicts),  and rebuild.

 And make sure there are no invisible non-ascii chars instead of spaces
 or tabs on that line.. but I doubt that.  can you email me the script
 or check with a hexdump program?


 worst come to worst you can add the kilometers text in with the
 text x% y% command by hand, then open the PS file in your favourite
 industrial strength text editor and search for the ZB (1) PB,
 (2), (3) m text and edit them to be like (10) (20) (30).
 [a wonderful trick for quickly fixing typos in the hardcopy]
 That should be pretty much near the end of the file.


 ???,
 Hamish









-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: r.neighbors modal filter for two raster layers

2009-04-16 Thread maning sambale
Revisiting this thread,  how do I add in the r.mapcalc command:

where if a pixel in the window contains a NULL value, do not count the
NULL in computing the mode

 For example if map 1 has many NULL values while map2 has valid
values, juts count the mode within map2.  This way I can patch areas
where map1 has NULLs.  And vise-versa if map2 has many NULLs.

The command below is OK when both maps have few NULL values.



On Fri, Jan 30, 2009 at 3:31 PM, maning sambale
emmanuel.samb...@gmail.com wrote:
 Well, it should be ,

 r.mapcalc test_andmode = mode (map1[-1,-1] , map1[-1,0] , map1[-1,1]
 , map1[0,-1] , map1[1,-1] , map1[1,1] , map1[0,0] , map1[0,1] ,
 map1[-1,0] , map2[-1,-1] , map2[-1,0] , map2[-1,1] , map2[0,-1] ,
 map2[1,-1] , map2[1,1] , map2[0,0] , map2[0,1] , map2[-1,0])



 On Fri, Jan 30, 2009 at 2:16 PM, maning sambale
 emmanuel.samb...@gmail.com wrote:
 What operator in mapcalc should I use to evaluate a mode from a given window:
 r.mapcalc test_andmode = mode (map1[-1,-1]  map1[-1,0] 
 map1[-1,1]  map1[0,-1]  \
 map1[1,-1]  map1[1,1]  map1[0,0]  map1[0,1]  map1[-1,0] 
 map2[-1,-1]  map2[-1,0]  \
 map2[-1,1]  map2[0,-1]  map2[1,-1]  map2[1,1]  map2[0,0] 
 map2[0,1]  map2[-1,0])

 Is the  correct to list values?

 On Fri, Jan 30, 2009 at 11:53 AM, maning sambale
 emmanuel.samb...@gmail.com wrote:
 Hi,

 Any ideas in doing this:

 I want to get the mode of raster categorical values for a defined
 window based on 2 raster layers.
 For example, I set a window of 5x5, so I get mode value from 50 pixels
 (25 for raster a and 25 for raster 2).

 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: r.neighbors modal filter for two raster layers

2009-04-16 Thread maning sambale
Again (doh!) answering to my own query.

use r.series method=mode

my grass module of the day!

On Thu, Apr 16, 2009 at 4:51 PM, maning sambale
emmanuel.samb...@gmail.com wrote:
 Revisiting this thread,  how do I add in the r.mapcalc command:

 where if a pixel in the window contains a NULL value, do not count the
 NULL in computing the mode

  For example if map 1 has many NULL values while map2 has valid
 values, juts count the mode within map2.  This way I can patch areas
 where map1 has NULLs.  And vise-versa if map2 has many NULLs.

 The command below is OK when both maps have few NULL values.



 On Fri, Jan 30, 2009 at 3:31 PM, maning sambale
 emmanuel.samb...@gmail.com wrote:
 Well, it should be ,

 r.mapcalc test_andmode = mode (map1[-1,-1] , map1[-1,0] , map1[-1,1]
 , map1[0,-1] , map1[1,-1] , map1[1,1] , map1[0,0] , map1[0,1] ,
 map1[-1,0] , map2[-1,-1] , map2[-1,0] , map2[-1,1] , map2[0,-1] ,
 map2[1,-1] , map2[1,1] , map2[0,0] , map2[0,1] , map2[-1,0])



 On Fri, Jan 30, 2009 at 2:16 PM, maning sambale
 emmanuel.samb...@gmail.com wrote:
 What operator in mapcalc should I use to evaluate a mode from a given 
 window:
 r.mapcalc test_andmode = mode (map1[-1,-1]  map1[-1,0] 
 map1[-1,1]  map1[0,-1]  \
 map1[1,-1]  map1[1,1]  map1[0,0]  map1[0,1]  map1[-1,0] 
 map2[-1,-1]  map2[-1,0]  \
 map2[-1,1]  map2[0,-1]  map2[1,-1]  map2[1,1]  map2[0,0] 
 map2[0,1]  map2[-1,0])

 Is the  correct to list values?

 On Fri, Jan 30, 2009 at 11:53 AM, maning sambale
 emmanuel.samb...@gmail.com wrote:
 Hi,

 Any ideas in doing this:

 I want to get the mode of raster categorical values for a defined
 window based on 2 raster layers.
 For example, I set a window of 5x5, so I get mode value from 50 pixels
 (25 for raster a and 25 for raster 2).

 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] convert r.report output to spreadsheet format

2009-02-24 Thread maning sambale
Hi,

Looking for tips in transforming r.report output into a
spreadsheet like format.

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] convert r.report output to spreadsheet format

2009-02-24 Thread maning sambale
Thanks!
An example of the r.report output:
|-|
|  Category Information   |  | cell|   %  |
| #|description   |  hectares|count| cover|
|-|
| 1|Siquijor  |31,313|   385511|  0.05|
|  |--|--|-|--|
|  |112|Broadleaf - open canopy. . . . . . . . . .|  8265|   101755| 26.39|
|  |140|Mangrove . . . . . . . . . . . . . . . . .|   173| 2125|  0.55|
|  |210|Agriculture-annual . . . . . . . . . . . .|  7370|90737| 23.54|
|  |310|Scrub. . . . . . . . . . . . . . . . . . .|  1063|13081|  3.39|
|  |320|Grassland. . . . . . . . . . . . . . . . .|  7185|88460| 22.95|
|  |400|Built-up . . . . . . . . . . . . . . . . .|   103| 1270|  0.33|
|  |510|Inland water . . . . . . . . . . . . . . .|44|  539|  0.14|
|  |  *|no data. . . . . . . . . . . . . . . . . .|  7111|87544| 22.71|

I want converted to:
Province,land_cover_code,description,hectares,cell_count,percent_cover

I will try the cryptic awk lines.


On Wed, Feb 25, 2009 at 12:12 AM, Giovanni Pasini jynx...@gmail.com wrote:
 maning sambale ha scritto:

 Hi,

 Looking for tips in transforming r.report output into a
 spreadsheet like format.


 Hi,
 I use awk (change NF and the other fields $*):

 r.report -n map=gri_1000_in,dem2_idw_100_slope_rec units=h,p,c | awk 'BEGIN
 {FS=|;OFS=|} !/---/  NF==7 {intest=$2} NF==8 {sub(/[ .]*$/,,$4);
 print intest, $4, $6}'  report_slope.txt

 This from a non awk expert!

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




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] sample bash script for multiple mapsets

2009-02-15 Thread maning sambale
In a single script I would like to:
open a mapset and do vector processing, then,
open another mapset,
import vectors from previous mapset
continue  processing.

Any psuedoscript to do this?

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] r.neighbors modal filter for two raster layers

2009-01-29 Thread maning sambale
Hi,

Any ideas in doing this:

I want to get the mode of raster categorical values for a defined
window based on 2 raster layers.
For example, I set a window of 5x5, so I get mode value from 50 pixels
(25 for raster a and 25 for raster 2).

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: r.neighbors modal filter for two raster layers

2009-01-29 Thread maning sambale
What operator in mapcalc should I use to evaluate a mode from a given window:
r.mapcalc test_andmode = mode (map1[-1,-1]  map1[-1,0] 
map1[-1,1]  map1[0,-1]  \
map1[1,-1]  map1[1,1]  map1[0,0]  map1[0,1]  map1[-1,0] 
map2[-1,-1]  map2[-1,0]  \
map2[-1,1]  map2[0,-1]  map2[1,-1]  map2[1,1]  map2[0,0] 
map2[0,1]  map2[-1,0])

Is the  correct to list values?

On Fri, Jan 30, 2009 at 11:53 AM, maning sambale
emmanuel.samb...@gmail.com wrote:
 Hi,

 Any ideas in doing this:

 I want to get the mode of raster categorical values for a defined
 window based on 2 raster layers.
 For example, I set a window of 5x5, so I get mode value from 50 pixels
 (25 for raster a and 25 for raster 2).

 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: r.neighbors modal filter for two raster layers

2009-01-29 Thread maning sambale
Well, it should be ,

r.mapcalc test_andmode = mode (map1[-1,-1] , map1[-1,0] , map1[-1,1]
, map1[0,-1] , map1[1,-1] , map1[1,1] , map1[0,0] , map1[0,1] ,
map1[-1,0] , map2[-1,-1] , map2[-1,0] , map2[-1,1] , map2[0,-1] ,
map2[1,-1] , map2[1,1] , map2[0,0] , map2[0,1] , map2[-1,0])



On Fri, Jan 30, 2009 at 2:16 PM, maning sambale
emmanuel.samb...@gmail.com wrote:
 What operator in mapcalc should I use to evaluate a mode from a given window:
 r.mapcalc test_andmode = mode (map1[-1,-1]  map1[-1,0] 
 map1[-1,1]  map1[0,-1]  \
 map1[1,-1]  map1[1,1]  map1[0,0]  map1[0,1]  map1[-1,0] 
 map2[-1,-1]  map2[-1,0]  \
 map2[-1,1]  map2[0,-1]  map2[1,-1]  map2[1,1]  map2[0,0] 
 map2[0,1]  map2[-1,0])

 Is the  correct to list values?

 On Fri, Jan 30, 2009 at 11:53 AM, maning sambale
 emmanuel.samb...@gmail.com wrote:
 Hi,

 Any ideas in doing this:

 I want to get the mode of raster categorical values for a defined
 window based on 2 raster layers.
 For example, I set a window of 5x5, so I get mode value from 50 pixels
 (25 for raster a and 25 for raster 2).

 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Announcing the release of QGIS 1.0 'Kore'

2009-01-25 Thread maning sambale
apologies for the cross-post:

Dear QGIS Users, Developers and FOSSGIS devotees

We are extremely pleased to announce the release of QGIS 1.0 and the
shiny new QGIS 1.0 User's Guide. We have also revamped our web site at
http://qgis.org. See the note below from our project chair:

-

This release is the culmination of literally hundreds of hours of work
by a team of developers, translators, documentation writers, and
graphics designers.

QGIS began life in February of 2002, with the first release in June of
the same year. The initial goal was to create a viewer for PostGIS data
that ran on Linux. From those humble beginnings, QGIS has become a true
cross-platform application that runs on all major versions of unix,
Linux, as well as Mac and Windows. It supports editing and map
composition as well as integration with GRASS to provide powerful GIS
capability. QGIS has been translated into many languages by a dedicated
internationalization team.

At 1.0 we provide a stable API from which you can develop custom
solutions in Python or C++. Even though 1.0 is fresh, there are a number
of exciting developments underway in both the core application and
plugins. Although it took nearly 7 years to get to version 1.0, I think
you'll find that this version is the best yet. Thanks for using
QGIS---you, the users, have played a large part in its success.


-- Gary Sherman, Chair, QGIS Project Steering Committee
-

For more detailed release information, visit:

http://blog.qgis.org/node/123

For Gnu/Linux, Mac OS X and Microsoft Windows downloads, visit:

http://download.qgis.org *

* Source code and User's Guide are also available from this site.

--
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] mosaicing classified images

2009-01-22 Thread maning sambale
But the landsat images are of different dates and in some case different years.

I found some mosaicing debates here:
http://www.faqs.org/faqs/sci/Satellite-Imagery-FAQ/part3/section-14.html

2 schools of thought appear:
1. Mosaic image then classify
2. Classify individual images then mosaic

I am following the great grass developer Michael Shapiro's (is he the
grass developer?) advice on option 2.

On Wed, Jan 21, 2009 at 7:16 PM, Markus Neteler nete...@osgeo.org wrote:
 On Wed, Jan 7, 2009 at 4:04 AM, maning sambale
 emmanuel.samb...@gmail.com wrote:
 Hi,
 I have a series of classified rasters from different landsat images
 (various dates).  I want to moasic/stitch them to a one raster.  The
 problem is there are subtle differences in the classified output. Any
 advise in mosaicing them?

 (I am afraid that it is easier to first mosaic (e.g., with OSSIM), then
  do the classification...)

 You could vectorize the raster areas and then try some SQL
 magic to match close polygons. Also v.clean to remove tiny
 areas from the map if needed.

 Markus




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Fwd: mosaicing classified images

2009-01-21 Thread maning sambale
Hi,
I am reposting this hoping for some answers.  Manually reclassifying
conflicting values seem to be a pain.  Any advice?



-- Forwarded message --
From: maning sambale emmanuel.samb...@gmail.com
Date: Wed, Jan 7, 2009 at 11:04 AM
Subject: mosaicing classified images
To: grass mailing list grass-user@lists.osgeo.org


Hi,
I have a series of classified rasters from different landsat images
(various dates).  I want to moasic/stitch them to a one raster.  The
problem is there are subtle differences in the classified output. Any
advise in mosaicing them?

--
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--



-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Why v.extract produces columns of type CHARACTER?

2009-01-14 Thread maning sambale
On Wed, Jan 14, 2009 at 5:02 PM, Glynn Clements
gl...@gclements.plus.com wrote:

 Nikos Alexandris wrote:

  The columns produces by v.extract are of type CHARACTER and v.dissolve
  does not like this. It's an old issue. Can someone explain why it
  becomes CHARACTER since grass' type for strings is varchar?
 
  Regards, Nikos

 Sorry, wrong question - false alarm.

 Actually, the *real* question I have asked in the past but never really
 got a reply is: why sqlitebrowser reports the columns as varchar
 while db.describe reports (in grass-shell) the same columns to be of
 type CHARACTER?

 GRASS' DBMI doesn't distinguish CHARACTER from CHARACTER VARYING (aka
 VARCHAR); the constant DB_SQL_TYPE_CHARACTER is used for both
 (although there is also DB_SQL_TYPE_TEXT for TEXT).

 Also, SQLite doesn't really have column types; it assigns types to
 individual values rather than to columns. Columns may have a type of
 SQLITE_INTEGER, SQLITE_FLOAT, SQLITE_TEXT, SQLITE_BLOB or SQLITE_NULL,
 but it doesn't actually require that the values' types conform to the
 type.
A bit off-topic.  Does this mean I can actually edit vector data
attributes as text even if the column is an INTEGER?
I'm using Sqlite database browser for editing.  A bit dangerous for
me.  Any other way to add checks/controls for editing mistakes?


 --
 Glynn Clements gl...@gclements.plus.com
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Why v.extract produces columns of type CHARACTER?

2009-01-14 Thread maning sambale
 IMHO, the best way to deal with attributes in GRASS + SQLite is to do it via
 the built-in tools in GRASS, and not via the SQLiteBrowser or other tools as
 this can lead to incompatibilities and confusion. The new wxgui has a very
 nice interface for attribute and table management.

Ahh! Thanks for the tip.  I haven't really used GRASS GUI apart from
the quick d.mon, d.vect, d.rast.
I guess now is the time to tryout wxgui.

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: mosaicing classified images

2009-01-12 Thread maning sambale
Hi,

I trying to mosaic my classified images, I have been experimenting
with David Finlayson's r.smoothpatch and r.eucdist.  This looks like
the tools I need.  However, the output map produces FCELL rasters.  My
classified image is an integer where integer values represent a
landcover category.

Is there a way to use this for integer-nominal data?

On Wed, Jan 7, 2009 at 11:04 AM, maning sambale
emmanuel.samb...@gmail.com wrote:
 Hi,
 I have a series of classified rasters from different landsat images
 (various dates).  I want to moasic/stitch them to a one raster.  The
 problem is there are subtle differences in the classified output. Any
 advise in mosaicing them?

 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] mosaicing classified images

2009-01-06 Thread maning sambale
Hi,
I have a series of classified rasters from different landsat images
(various dates).  I want to moasic/stitch them to a one raster.  The
problem is there are subtle differences in the classified output. Any
advise in mosaicing them?

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] osm data in grass

2008-12-11 Thread maning sambale
Hi,

AFAIK the v.in.osm is for the old API (0.4) of OSM.  You can't use it
for the current data. There are utilities to import OSM data to
various formats checkout their wiki:
http://wiki.openstreetmap.org/

But if you want GIS ready data derived from OSM, you can download
shapefiles here:
http://download.geofabrik.de/
and here
http://downloads.cloudmade.com/


On Fri, Dec 12, 2008 at 5:18 AM, Torsten Eckart tore...@arcor.de wrote:
 Hallo,
 could someone give me a hint how to import osm data in to a grass
 location.
 The module v.in.osm what I download from
 http://kripton.kripserver.net/software/v.in.osm/v.in.osm-2007-07-14-00
 gives not a good result:
 Reading from file /home/torsten/Arbeit/Gis/Daten/osm/heim.osm, writing
 to map heim Reading nodes ... [ OK ] (339 nodes)
 Reading segs ... [ OK ] (0 segments)
 Importing way 1/77 (ID: 3213956) can't read segcounter: no such
 variable while executing
 puts -nonewline ($segcounter segments) ...
(while body line 58)
invoked from within
 while {![eof $fchan]} {
  gets $fchan line
  if {[string match *way\ id=* $line]} {
incr waycount
set list [split $line]
foreach elem $li...
(file /usr/lib/grass/scripts/v.in.osm line 185)

 Additional I do not understand what properties the target location has
 to have. Do's v.in.osm transform the data into the right location e.g.:
 projection: 99 (Transverse Mercator)
 zone:   0
 datum:  potsdam
 ellipsoid:  bessel
 north:  5661000
 south:  566
 west:   5411000
 east:   5412000
 nsres:  1
 ewres:  1
 rows:   1000
 cols:   1000
 cells:  100

 Thanks for Your help.
 Torsten
 --
 Chatterreichbarkeit:
 e...@jabber.chaotic.de
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] LANDSAT data archive now free

2008-12-11 Thread maning sambale
We can now access the full 35+ years of archive Landsat 1 through 7 for free.

-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] LANDSAT data archive now free

2008-12-11 Thread maning sambale
Here:
http://glovis.usgs.gov/

or

http://edcsns17.cr.usgs.gov/EarthExplorer/

On Fri, Dec 12, 2008 at 2:11 PM, Rainer M Krug r.m.k...@gmail.com wrote:
 On Fri, Dec 12, 2008 at 7:38 AM, maning sambale
 emmanuel.samb...@gmail.com wrote:
 We can now access the full 35+ years of archive Landsat 1 through 7 for free.

 Sounds great - could you provide a URL?


 --
 cheers,
 maning
 --
 Freedom is still the most radical idea of all -N.Branden
 wiki: http://esambale.wikispaces.com/
 blog: http://epsg4253.wordpress.com/
 --
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user




 --
 Rainer M. Krug, PhD (Conservation Ecology, SUN), MSc (Conservation
 Biology, UCT), Dipl. Phys. (Germany)

 Centre of Excellence for Invasion Biology
 Faculty of Science
 Natural Sciences Building
 Private Bag X1
 University of Stellenbosch
 Matieland 7602
 South Africa




-- 
cheers,
maning
--
Freedom is still the most radical idea of all -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: combining and evaluating vector dbf values

2008-11-25 Thread maning sambale
Nearly there but not quite:

The final column should combine the values from the 2 columns, am I
right that concat doesn't work with dbf?

So this doesn't work
echo UPDATE vectormap SET c_COVER = a_COVER || a_value | db.execute

cheers,
maning
On Tue, Nov 25, 2008 at 2:34 PM, maning sambale
[EMAIL PROTECTED] wrote:
 Hi,

 Answering partially to my problem.

 I found a way to update columns by from a criteria of another column

 like this:
 v.info -c overlay_trainingANDforest2002_sq
 lDisplaying column types/names for database connection of layer 1:
 INTEGER|cat
 INTEGER|a_cat
 INTEGER|a_value
 CHARACTER|a_label
 INTEGER|a_COVER
 INTEGER|b_cat
 CHARACTER|b_ID
 CHARACTER|b_TYPE
 INTEGER|c_COVER

 v.db.update map=overlay_trainingANDforest2002_sql col=c_COVER
 value=112 where=a_Cover=112

 And It's really fast for 100K records!  Now, how do make an sql
 expression by evaluationg multiple columns?

 cheers,
 maning


 On Tue, Nov 25, 2008 at 12:01 PM, maning sambale
 [EMAIL PROTECTED] wrote:
 Hi,

 Any help in doing the following:

 I have a combined vector in grass using v.overlay

 I now have a dbf table with the following
 a_value = integer
 a_value2 = integer
 b_value = string

 Now I want to combine the categories in another column (c_value) like :

 evaluate b_value, a_value, a_value2
 then assign c_value from a lookup table

 the lookup table is a set of integer values based on the various
 combinations of b_value, a_value, a_value2

 Any ideas doing this in GRASS or oodbase?

 cheers,
 maning


 --
 |-|--|
 | __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
 | '-._7' |Freedom is still the most radical idea of all -N.Branden|
 |  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
 |  |  /T  |http://esambale.wikispaces.com/ |
 | _)_/L I http://epsg4253.wordpress.com/ |
 |-|--|




 --
 |-|--|
 | __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
 | '-._7' |Freedom is still the most radical idea of all -N.Branden|
 |  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
 |  |  /T  |http://esambale.wikispaces.com/ |
 | _)_/L I http://epsg4253.wordpress.com/ |
 |-|--|




-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com/ |
| _)_/L I http://epsg4253.wordpress.com/ |
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] combining and evaluating vector dbf values

2008-11-24 Thread maning sambale
Hi,

Any help in doing the following:

I have a combined vector in grass using v.overlay

I now have a dbf table with the following
a_value = integer
a_value2 = integer
b_value = string

Now I want to combine the categories in another column (c_value) like :

evaluate b_value, a_value, a_value2
then assign c_value from a lookup table

the lookup table is a set of integer values based on the various
combinations of b_value, a_value, a_value2

Any ideas doing this in GRASS or oodbase?

cheers,
maning


-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com/ |
| _)_/L I http://epsg4253.wordpress.com/ |
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: combining and evaluating vector dbf values

2008-11-24 Thread maning sambale
Hi,

Answering partially to my problem.

I found a way to update columns by from a criteria of another column

like this:
v.info -c overlay_trainingANDforest2002_sq
lDisplaying column types/names for database connection of layer 1:
INTEGER|cat
INTEGER|a_cat
INTEGER|a_value
CHARACTER|a_label
INTEGER|a_COVER
INTEGER|b_cat
CHARACTER|b_ID
CHARACTER|b_TYPE
INTEGER|c_COVER

v.db.update map=overlay_trainingANDforest2002_sql col=c_COVER
value=112 where=a_Cover=112

And It's really fast for 100K records!  Now, how do make an sql
expression by evaluationg multiple columns?

cheers,
maning


On Tue, Nov 25, 2008 at 12:01 PM, maning sambale
[EMAIL PROTECTED] wrote:
 Hi,

 Any help in doing the following:

 I have a combined vector in grass using v.overlay

 I now have a dbf table with the following
 a_value = integer
 a_value2 = integer
 b_value = string

 Now I want to combine the categories in another column (c_value) like :

 evaluate b_value, a_value, a_value2
 then assign c_value from a lookup table

 the lookup table is a set of integer values based on the various
 combinations of b_value, a_value, a_value2

 Any ideas doing this in GRASS or oodbase?

 cheers,
 maning


 --
 |-|--|
 | __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
 | '-._7' |Freedom is still the most radical idea of all -N.Branden|
 |  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
 |  |  /T  |http://esambale.wikispaces.com/ |
 | _)_/L I http://epsg4253.wordpress.com/ |
 |-|--|




-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com/ |
| _)_/L I http://epsg4253.wordpress.com/ |
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] New user grass help in Australia

2008-11-23 Thread maning sambale
If you can get hold of the GRASS bible by Markus and Helena:

http://www.grassbook.org/

You're on your way to become a great grass user.


On Sun, Nov 23, 2008 at 3:49 PM, Markus Neteler [EMAIL PROTECTED] wrote:
 Hi Matt,

 On Sat, Nov 22, 2008 at 10:55 AM, Matt B [EMAIL PROTECTED] wrote:
 Hi Group,
 I've been plugging away with grass for a few months now and not making
 progress as fast as I would like. Is there a user group (preferably in
 Melbourne though Adelaide is doable) or someone who would be able to sit
 down with me for an hour or two to help me through a few startup problems.

 If all are physically to remote, consider to post here.

 Even having someone on chat (I haven't had much luck on the grass irc
 channel, maybe because of the time difference?) walk me through a few push
 this button steps would be useful.

 The GRASS IRC channel
  http://grass.itc.it/community/support.php#irc
  isn't too crowded, the mailing lists are usually a good place to ask.

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




-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com/ |
| _)_/L I http://epsg4253.wordpress.com/ |
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] GRASS tools for palsar/radar images

2008-11-20 Thread maning sambale
Hi,

I'm looking for some tools in processing radar/palsar data.
pre-processing, image classification, etc.
So far, I found RAT radartools http://srv-43-200.bv.tu-berlin.de/rat/index.php

Any GRASS user doing radar processing?

cheers,
maning

-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com/ |
| _)_/L I http://epsg4253.wordpress.com/ |
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] quickbird imagery: i.fusion.brovey, r.composite...

2008-11-04 Thread maning sambale
Hi,

(Somewhat related but not entirely)

I am looking for cases/howto/examples of using quickbird images in
GRASS.  Please refer some related:

orthorectification of quickbird in GRASS
change detection using quickbird in GRASS
image classification with quickbird (using texture, geometry and
spectra) in GRASS

Link I found so far:
http://geoinformatics.fsv.cvut.cz/wiki/index.php/Change_Detection_with_GRASS_GIS_-_Comparison_of_images_taken_by_different_sensors


cheers,
maning

On Wed, Nov 5, 2008 at 6:12 AM, tommaso [EMAIL PROTECTED] wrote:
 Hi List,

 I'm trying to display a Quickbird imagery with nice real colors.

 I tried with:

 r.composite red=m1.red green=m1.green blue=m1.blue out=composite
 i.fusion.brovey ms1=m1.green ms2=m1.nir ms3=m1.red pan=pan
 outputprefix=fusion -q

 I tried also to set the RGB maps to grey und grey.eq but unfortunately until
 now the results are not very satisfying.

 Any suggestion will be appreciated.

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





-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com|
| _)_/LI
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] d.out.file now does GeoTIFFs

2008-10-29 Thread maning sambale
Hamish,

Just tried it.  This makes a quick and dirty visualization really
fast!  Especially for non GIS colleagues.

I can simply import the d.mon display to geotiff then run
gdal2tile.py.  VoilaI  I have a KML of the output map for viewing in
GoogleEarth.

Thanks!

cheers,
maning

On 10/29/08, Hamish [EMAIL PROTECTED] wrote:
 Nikos Alexandris:
  Well,
  so I did.

 maning sambale wrote:
 How?  I haven't been doing svn update for the past 2
 months because its working OK.
 Would love to try this out also.


 if you are working from a svn source tree it is just:

 cd scripts/d.out.file
 svn up
 make


 you only have to do ./configure again and rebuild everything for deep
 changes. (if you haven't in a few months it is probably a good idea
 thought)


 and whatever you do to make install


 but this one is just a script, you can download the latest version and
 simply replace the file. (remember 'chmod a+x')


 Hamish








-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com|
| _)_/LI
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] d.out.file now does GeoTIFFs

2008-10-29 Thread maning sambale
Hamish,

I just found out that the output from d.out.file geotiff has some distortions.

See this: http://farm4.static.flickr.com/3007/2984908579_b73cdd3b99_o.png

The lighter relief is the original raster from the GRASS dbase and the
darker one with some green vector is the d.out.file output.  I'm using
qgis to visualize.

cheers,
maning

On Thu, Oct 30, 2008 at 5:03 AM, Glynn Clements
[EMAIL PROTECTED] wrote:

 Hamish wrote:

   Unfortunately 'd.info -g' is buggy, it reports exactly the same as
   'g.region -g'. But if that was fixed it would be an easy solution.

 fixed in devbr6 but not trunk; I'm unsure about the changes needed.

 7.0 doesn't have D_get_screen_window(); D_setup() uses R_get_window()
 instead (R_get_window() reports the raster clip region, which is
 initially set from $GRASS_FRAME).

 Also, I would suggest reporting the values returned from D_get_u_*
 rather than using the current region directly (which can always be
 obtained from g.region).

 Although the two should match, it's the values reported by D_get_u_*
 which are actually used, so using D_get_u_* would be more robust in
 the face of any future extensions (e.g. allowing the region to be
 enlarged to fit the window rather than shrinking the window to fit the
 region).

 --
 Glynn Clements [EMAIL PROTECTED]
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user




-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com|
| _)_/LI
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Landsat SLC-Off gap filling

2008-10-20 Thread maning sambale
+1

I had just downloaded new Landsat 7 SLC-Off data.  It still useable
for land use change analysis.  A grass tool would be very helpful.

cheers,
maning




On 10/20/08, Ned Horning [EMAIL PROTECTED] wrote:
 Hi - Is there anyone out there who has worked on an algorithm to fill gaps
 in
 Landsat ETM+ SLC-Off data? Now that the Landsat 7 archive is freely
 available
 I am looking for an open source approach to fill the gaps in the SLC-Off
 data
 and GRASS seems like a good option. I'm most interested in multi-date
 approaches
 but at this stage I'm looking at single-date approaches as well. Any ideas?


 Ned

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



-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com|
| _)_/LI
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Problem grass-6.4.0 spearfish60 aspect raster image looks heavily pixellated

2008-09-30 Thread maning sambale
Elvis,

You region settings says:
EW Res: 149.++
NS Res: 150.++

Set g.region to resolution of the layer
g.region rast=aspect -p

cheers,
maning





On 9/30/08, Elvis Dowson [EMAIL PROTECTED] wrote:
 Hi,
   I just tried to display the spearfish60 aspect raster image,
 but it looks heavily pixellated. What could the problem be?




-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com|
| _)_/LI
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] dbf editor for mac for grass and qgis

2008-09-25 Thread maning sambale
Hi,

My usual workflow is like this:

All my vectors are in GRASS database and I view them with QGIS.
For editing vector attributes, I basically browse the GRASS vectors in
QGIS then edit the dbf attributes with OpenOffice calc or base.
So far so good.  This time I have a vector with more than 100K
records.  Editing with calc is not possible (hitting the limit of 65K
rows) and its painstakingly slow in OO base.

I'm looking for ways to speed things up.  I cannot find a decent dbf
editor for mac.  Any ideas?

cheers,
maning

-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com|
| _)_/LI
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: [Qgis-user] dbf editor for mac for grass and qgis

2008-09-25 Thread maning sambale
Just finished installing macports gnumeric and it does have the 65K row limit

On Thu, Sep 25, 2008 at 4:56 PM, maning sambale
[EMAIL PROTECTED] wrote:
 Thanks for that tip.  I'm installing from macports now.
 But I read somewhere that: GNUmeric could accept more than 256 columns
 and more of 65000 rows (needs to be re-compiled)

 On Thu, Sep 25, 2008 at 4:47 PM, Nikos Alexandris
 [EMAIL PROTECTED] wrote:
 On Thu, 2008-09-25 at 16:38 +0800, maning sambale wrote:
 Hi,

 My usual workflow is like this:

 All my vectors are in GRASS database and I view them with QGIS.
 For editing vector attributes, I basically browse the GRASS vectors in
 QGIS then edit the dbf attributes with OpenOffice calc or base.
 So far so good.  This time I have a vector with more than 100K
 records.  Editing with calc is not possible (hitting the limit of 65K
 rows) and its painstakingly slow in OO base.
 Can you use Gnumeric? It's faster and limitless :-) ...ok, maybe no
 limitless but reads bigger files than OOo

 I'm looking for ways to speed things up.  I cannot find a decent dbf
 editor for mac.  Any ideas?

 cheers,
 maning






 --
 |-|--|
 | __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
 | '-._7' |Freedom is still the most radical idea of all -N.Branden|
 |  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
 |  |  /T  |http://esambale.wikispaces.com|
 | _)_/LI
 |-|--|




-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com|
| _)_/LI
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: [Qgis-user] dbf editor for mac for grass and qgis

2008-09-25 Thread maning sambale
Andreas,

I agree with you.  I actually have postgresql installed on my mac.
But this would just be a temporary vector file (an input to a
satellite classification).  So I am looking (if there are) ways just
to make it usable for now.  I might consider postgis if all quick and
dirty efforts would be fruitless

cheers,
maning

On Thu, Sep 25, 2008 at 5:57 PM, Andreas Neumann [EMAIL PROTECTED] wrote:
 With this many records I would seriously consider moving to a professional
 database, such as PostgreSQL/Postgis or Oracle Spatial. Both should run on a
 Mac. Postgis is of course easier to install and maintain than Oracle. It is
 also cheaper.

 Andreas

 maning sambale wrote:

 Just finished installing macports gnumeric and it does have the 65K row
 limit

 On Thu, Sep 25, 2008 at 4:56 PM, maning sambale
 [EMAIL PROTECTED] wrote:

 Thanks for that tip.  I'm installing from macports now.
 But I read somewhere that: GNUmeric could accept more than 256 columns
 and more of 65000 rows (needs to be re-compiled)

 On Thu, Sep 25, 2008 at 4:47 PM, Nikos Alexandris
 [EMAIL PROTECTED] wrote:

 On Thu, 2008-09-25 at 16:38 +0800, maning sambale wrote:

 Hi,

 My usual workflow is like this:

 All my vectors are in GRASS database and I view them with QGIS.
 For editing vector attributes, I basically browse the GRASS vectors in
 QGIS then edit the dbf attributes with OpenOffice calc or base.
 So far so good.  This time I have a vector with more than 100K
 records.  Editing with calc is not possible (hitting the limit of 65K
 rows) and its painstakingly slow in OO base.

 Can you use Gnumeric? It's faster and limitless :-) ...ok, maybe no
 limitless but reads bigger files than OOo

 I'm looking for ways to speed things up.  I cannot find a decent dbf
 editor for mac.  Any ideas?

 cheers,
 maning




 --
 |-|--|
 | __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
 | '-._7' |Freedom is still the most radical idea of all -N.Branden|
 |  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
 |  |  /T  |http://esambale.wikispaces.com|
 | _)_/LI
 |-|--|





 --

 --
 Andreas Neumann
 Böschacherstrasse 6
 CH-8624 Grüt (Gossau ZH)
 Switzerland
 Phone: ++41-44-2736668
 Email: [EMAIL PROTECTED]

 Web: http://www.carto.net/neumann/
 SVG Examples: http://www.carto.net/papers/svg/samples/
 SVG.Open: http://www.svgopen.org/




-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com|
| _)_/LI
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: [Qgis-user] dbf editor for mac for grass and qgis

2008-09-25 Thread maning sambale
On Thu, Sep 25, 2008 at 6:33 PM, Nikos Alexandris
[EMAIL PROTECTED] wrote:
 On Thu, 2008-09-25 at 17:53 +0800, maning sambale wrote:
 Just finished installing macports gnumeric and it does have the 65K row limit


 Oh! I am so sorry for the mis-information. I remembered a discussion we
 had at some point about those limits. The problem back then was about
 the limit of columns (256).The oss solution for that is 32767 for
 KSpread.

 Sorry,
 Nikos

 No need for apologies,  I am actually pushing the limits to
spreadsheet apps, which shouldn't be.





-- 
|-|--|
| __.-._  |Ohhh. Great warrior. Wars not make one great. -Yoda |
| '-._7' |Freedom is still the most radical idea of all -N.Branden|
|  /'.-c  |Linux registered user #402901, http://counter.li.org/ |
|  |  /T  |http://esambale.wikispaces.com|
| _)_/LI
|-|--|
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


  1   2   >