Re: [GRASS-user] r.basin matplotlib error
Dear Margherita, Thanks! On Thu, Apr 30, 2015 at 4:35 PM, Margherita Di Leo 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] r.basin matplotlib error
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 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 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
[GRASS-user] SRTM 1-arc (~30m) available outside USA and Africa
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
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:
[GRASS-user] GRASS71 development packages for Ubuntu 13.10 and 14.04
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
Re: [GRASS-user] 10day accumulation which t.rast module to use?
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
Re: [GRASS-user] 10day accumulation which t.rast module to use?
> > 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?
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 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] 10day accumulation which t.rast module to use?
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] temporal mapcalc 3hour precipitation to hourly
Thanks Moritz. I renamed all my maps and its working now. On Tue, Dec 2, 2014 at 8:59 PM, Moritz Lennert wrote: > On 02/12/14 13:30, maning sambale wrote: >> >> 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! > > > AFAIK, r.mapcalc has difficulties dealing with map names that contain a dash > in the name as it interprets it as a minus sign. > > Try renaming your maps. > > 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] temporal mapcalc 3hour precipitation to hourly
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:
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 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
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
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 wrote: > On Mon, Nov 3, 2014 at 10:18 AM, maning sambale > 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
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
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?
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
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] Global Land Cover Characterization data type
On Fri, Feb 15, 2013 at 11:40 AM, Yann Chemin wrote: > (maybe in r.in.gdal, maybe another more specific. Or, r.in.bin -- 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
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
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
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
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 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
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 K&C 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 K&C 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 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 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
Re: [GRASS-user] grass and palsar processing
Sweet! Thanks markus. On Feb 10, 2012 7:16 PM, "Markus Neteler" wrote: > On Fri, Feb 10, 2012 at 12:15 PM, Martin Landa > wrote: > > 2012/2/10 Markus Neteler : > >> 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
[GRASS-user] grass and palsar processing
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] measuring area on a 3d polygon
Thanks, will try this.
On Thu, Nov 17, 2011 at 3:41 AM, Benjamin Ducke 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"
> 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"
> 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
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
[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
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 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
r.mapcalc "integer_map = round(fp_maps)" On Thu, May 19, 2011 at 7:04 AM, Michael Morin 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
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 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
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 wrote: > Anyone successfully loaded kc50 palsar raw imagery with gdal? > > > -- Forwarded message ------ > From: maning sambale > 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
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
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 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
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
Re: [GRASS-user] [DWE]- Landsat Cloud detection with GRASS
2010/9/29 Markus Neteler : > 2010/9/28 António Rocha : >> 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
[GRASS-user] multi-date image segmentation in grass
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] Re: Convert FCELL to CELL
r.mapcalc new_map="round(original_map)" On 6/30/10, Pedro Venâncio 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
Thanks! I'm trying this route. On Wed, May 26, 2010 at 11:18 PM, Dylan Beaudette 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
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] Re: r.coin behaving differently
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 wrote: > Follow-up > > On Wed, Apr 14, 2010 at 5:41 PM, maning sambale > 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
[GRASS-user] Re: r.coin behaving differently
Follow-up On Wed, Apr 14, 2010 at 5:41 PM, maning sambale 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] r.coin behaving differently
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
Re: [GRASS-user] claslite modules in GRASS?
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 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?
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?
> 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?
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
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 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 406
[GRASS-user] weird results with mapcalc on different region extents
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
Re: [GRASS-user] Help starting grass for first time
If you can buy/borrow the GRASS Book, you have everything you need. :) http://www.grassbook.org/ On Wed, Dec 9, 2009 at 10:06 AM, Ross Benisch wrote: > Hi, I have downloaded and unzip all the data for the 6.4 version of grass > for windows vista and I can not get the program to run. I try and set the > projection and all i get is an error. If someone would describe the process > of how to start grass for the first time for my setup, I would greatly > appreciate it because it is very frustrating to not even be able to open the > software I want to start using it. I have tried all of the help and start > up guides and they are not solving the problem. Please help. > > Ross > > > Windows Live Hotmail gives you a free,exclusive gift. Click here to > download. > ___ > 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
try this: r.proj input=dem dbase=/work/projects/ location=Willamette mapset=PERMANENT On Wed, Dec 2, 2009 at 2:45 PM, maning sambale wrote: > On Wed, Dec 2, 2009 at 2:28 PM, Janet Choate 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 >> wrote: >>> >>> Try to include the full path of your location >>> ([/work/projects/Willamette]) >>> >>> On Wed, Dec 2, 2009 at 1:47 PM, Janet Choate 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
Re: [GRASS-user] reproject from lat/long to utm
On Wed, Dec 2, 2009 at 2:28 PM, Janet Choate 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 > wrote: >> >> Try to include the full path of your location >> ([/work/projects/Willamette]) >> >> On Wed, Dec 2, 2009 at 1:47 PM, Janet Choate 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
Try to include the full path of your location ([/work/projects/Willamette]) On Wed, Dec 2, 2009 at 1:47 PM, Janet Choate 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
[GRASS-user] internal buffer of polygon
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
It works using ogr2ogr. On Thu, Oct 15, 2009 at 12:45 AM, Hermann Peifer 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 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 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
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 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 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?
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?
v.report map=name option=length units=meters On 9/17/09, Margherita Di Leo 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=visualizzaHomePage&id=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$
afaik, it should be: v.out.ogr -e input=hidr...@permanent dsn=C:/tmp/hidro3D On 9/17/09, Vincent Blanqué 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
Re: [GRASS-user] creating map series purely on ps.map
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 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/IconS
[GRASS-user] creating map series purely on ps.map
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] change some labels to upper/lower case in vector attributes
Thanks will try this one. On Tue, Aug 18, 2009 at 7:36 PM, Moritz Lennert 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
Re: [GRASS-user] updating vector attributes for separate islands
Thanks! it's working now On 8/17/09, Markus Metz 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] change some labels to upper/lower case in vector attributes
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
[GRASS-user] updating vector attributes for separate islands
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
Re: [GRASS-user] bus error in v.in.ogr
On Fri, Aug 14, 2009 at 9:39 PM, maning sambale 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 > Ribeiro 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 >>> >>> 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 already exists and will be >>> overwritten >>> Projection of input dataset and current location appear to match >>> WARNING: Vector map 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 ... >>> 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 ... >>> 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] bus error in v.in.ogr
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 already exists and will be overwritten Projection of input dataset and current location appear to match WARNING: Vector map 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 ... 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 ... 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
[GRASS-user] wishlist to add a grass module for Landsat ETM+ SLC-off Gapfill
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
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
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 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
Re: [GRASS-user] error in ps.map scalebar units argument
GRASS 6.4.svn (pfua_combine):~/pfua/pfua_layout/250K_layout > ./250k_layout.map ERROR: units kilometers : illegal request (scalebar) projection: 1 (UTM) zone: 51 datum: wgs84 ellipsoid: wgs84 north: 1108130.80088945 south: 995437.83819297 west: 718326.7686683 east: 886033.01829335 nsres: 99.99375572 ewres: 100.00372667 rows: 1127 cols: 1677 cells: 1889979 On Tue, May 26, 2009 at 5:04 PM, Hamish wrote: > >> > maning wrote: >> >> scalebar s >> >> where 15.0 20.70 >> >> length 3o >> >> height 0.12 >> >> segment 6 >> >> numbers 2 >> >> units kilometers >> >> fontsize 15 >> >> end >> > >> > length "3o" -> "30"? > > maning wrote: >> Oops. corrected this one still the same error >> >> Disabling the units instruction gives me scalebar but with >> the default unit (in this case meters, I want the scale in KM) > > > what is the exact error message? (cut&paste) > what version of grass do you use? > what paper size? > what are the region extents? > > > (it works for me) > > 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
Re: [GRASS-user] error in ps.map scalebar units argument
Oops. corrected this one still the same error Disabling the units instruction gives me scalebar but with the default unit (in this case meters, I want the scale in KM) On Tue, May 26, 2009 at 3:53 PM, Hamish wrote: > > maning wrote: >> scalebar s >> where 15.0 20.70 >> length 3o >> height 0.12 >> segment 6 >> numbers 2 >> units kilometers >> fontsize 15 >> end > > > length "3o" -> "30"? > > > 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] error in ps.map scalebar units argument
scalebar s where 15.0 20.70 length 3o height 0.12 segment 6 numbers 2 units kilometers fontsize 15 end -- 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] use rgbcolumn for vector line in vareas ps.map
Hamish, Thanks for pointing out the obvious color none! Follow-up Q, how do I exclude polygon holes in vareas? I'm starting to really like ps.map for batch processing map sheets :) On Fri, May 22, 2009 at 2:33 PM, Hamish wrote: > > maning wrote: >> I have an vector polygon where an attribute column defines the rgb >> color for ps.map. The problem is, vareas options does not allow to >> use the rgbcolumn for line color. Setting the width to very small >> doesn't help either. Is there a way to disable line render in vareas. > > to turn off boundaries set "color none" > > rgbcolumn also works for pattern color by the way. > > you can use a second vlines command with "type boundary" to draw > boundary lines, but typically those will not have category numbers > (and thus no GRASSRGB entry) because it is ambiguous as to which area > the boundary belongs to, i.e. take the color from the left or right > side of the line? > > > for vpoints see also http://trac.osgeo.org/grass/ticket/342 > > >> Or should I split the vector polys into separate layer of polygon >> categories (too many). > > where cat IS >= x AND cat IS <= y > ? (or similar SQL query for repeated vareas calls) > > > > 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] use rgbcolumn for vector line in vareas ps.map
Hi, I have an vector polygon where an attribute column defines the rgb color for ps.map. The problem is, vareas options does not allow to use the rgbcolumn for line color. Setting the width to very small doesn't help either. Is there a way to disable line render in vareas. Or should I split the vector polys into separate layer of polygon categories (too many). example below: vareas forestcover_2002 masked y width .001 rgbcolumn rgb cats 1-999 label Forest end Any other elegant solution (within ps.map)? Using GRASS6.4 svn compiled ages ago. -- 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] installing grass trunk to jaunty
mitch,
This works for me without a hitch! thanks!
On Thu, Apr 30, 2009 at 3:53 AM, mitch wrote:
> hi maning, i haven't been playing so much with it (busy)
> but i got grass7 svn work on jaunty
>
> gdal and proj from synaptic
>
>
> DEPENDENCIES
>
>
> sudo apt-get install build-essential \
> flex \
> lesstif2-dev \
> libfftw3-dev \
> libglu1-mesa-dev \
> libglu1-xorg-dev \
> libpng12-dev \
> libpq-dev \
> unixodbc-dev \
> doxygen \
> libmysqlclient15-dev \
> graphviz \
> libcairo2-dev \
> bison \
> libncurses5-dev \
> zlib1g-dev \
> libjpeg62-dev \
> libgdal1-dev \
> libtiff4-dev \
> libgcc1 \
> tcl8.4-dev \
> tk8.4-dev \
> fftw3-dev \
> libfreetype6-dev \
> libavcodec-dev \
> libxmu-dev \
> libreadline5 \
> libreadline5-dev \
> python2.5-dev \
> python-wxversion \
> swig \
> libwxbase2.8-0 \
> libwxbase2.8-dbg \
> libwxbase2.8-dev \
> libwxgtk2.8-0 \
> libwxgtk2.8-dbg \
> libwxgtk2.8-dev \
> python-wxgtk2.8 \
> wx2.8-doc \
> wx2.8-examples \
> wx2.8-headers \
> wx2.8-i18n \
> python-numpy \
> sqlite3 \
> libsqlite3-dev \
> python-opengl
>
>
> svn checkout https://svn.osgeo.org/grass/grass/trunk grass_trunk
> (wherever you want...i did it in /usr/local/src)
>
>
>
> CFLAGS="-g -Wall" ./configure \
> --with-gdal=/usr/local/bin \
> -with-tcltk-includes=/usr/include/tcl8.4 \
> --with-sqlite \
> --with-cxx \
> --with-fftw \
> --with-freetype=yes \
> --with-freetype-includes=/usr/include/freetype2 \
> --with-readline \
> --with-opengl-includes=/usr/include/GL \
> --with-python=/usr/bin/python2.5-config \
> --with-wxwidgets=/usr/bin/wx-config \
> --enable-largefile=yes \
> --with-proj-includes=/usr/include \
> --with-proj-share=/usr/share/proj \
> --with-proj-libs=/usr/lib
>
>
>
>
> GRASS is now configured for: i686-pc-linux-gnu
>
> Source directory: /usr/local/src/grass_trunk
> Build directory: /usr/local/src/grass_trunk
> Installation directory: ${prefix}/grass-7.0.svn
> Startup script in directory: ${exec_prefix}/bin
> C compiler: gcc -g -Wall
> C++ compiler: c++ -g -O2
> Building shared libraries: yes
> 64bit support: no
> OpenGL platform: X11
>
> MacOSX application: no
>
> NVIZ: yes
>
> BLAS support: no
> C++ support: yes
> Cairo support: yes
> DWG support: no
> FFMPEG support: no
> FFTW support: yes
> FreeType support: yes
> GDAL support: yes
> GEOS support: no
> JPEG support: yes
> LAPACK support: no
> Large File support (LFS): yes
> MySQL support: no
> NLS support: no
> ODBC support: no
> OGR support: yes
> OpenGL support: yes
> PNG support: yes
> PostgreSQL support: no
> Python support: yes
> Readline support: yes
> SQLite support: yes
> Tcl/Tk support: yes
> wxWidgets support: yes
> TIFF support: yes
> X11 support: yes
> Regex support: yes
> POSIX thread support: no
>
>
>
> sudo make
>
>
> Errors in:
> No errors detected.
> (hopefully :-))
>
> sudo make install
>
>
> this is what i did (indeed you can change parameters in ./configure)
> btw community, AM I RIGHT?...did i miss something or include to
> much? :-) (suggestion appreciated)
>
>
> another question:
> what do i have to do exactly to update?
>
> svn up
> ./configure (including the same parameters of the previous
> installation?)
> make
> make install (again?)
>
> ..or there is another way to do it?
>
>
> cheers,
> mitch
>
>
>
> On Thu, 2009-04-30 at 01:28 +0800, maning sambale wrote:
>> Compiling from source following:
>>
>> http://grass.osgeo.org/wiki/Compile_and_Install#Ubuntu
>>
>> I am somehow stuck here:
>> ./configure --enable-gpl --enable-libvorbis --enable-libtheora \
>> --enable-liba52 --enable-libdc1394 --enable-libgsm \
>> --enable-libmp3lame --enable-libfaad --enable-libfaac \
>> --enable-libxvid --enable-pthreads --enable-libx264 \
>> --enable-shared
>>
>> Unknown option "--enable-liba52".
>> See ./configure --help for available options.
>>
>>
>> ps.
>> in
[GRASS-user] installing grass trunk to jaunty
Compiling from source following: http://grass.osgeo.org/wiki/Compile_and_Install#Ubuntu I am somehow stuck here: ./configure --enable-gpl --enable-libvorbis --enable-libtheora \ --enable-liba52 --enable-libdc1394 --enable-libgsm \ --enable-libmp3lame --enable-libfaad --enable-libfaac \ --enable-libxvid --enable-pthreads --enable-libx264 \ --enable-shared Unknown option "--enable-liba52". See ./configure --help for available options. ps. in the GRASS dependencies: I think "bpng12-dev" should be "libpng12-dev" sudo apt-get install grass build-essential flex bison libncurses5-dev zlib1g-dev \ libjpeg62-dev libgdal1-dev libtiff4-dev libgcc1 bpng12-dev tcl8.4-dev tk8.4-dev fftw3-dev \ libfreetype6-dev libavcodec-dev libxmu-dev gdal-bin libreadline5 libreadline5-dev \ make python-dev python-wxversion swig -- 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 going to python scripting (was Re: [GRASS-user] new module: r.surf.volcano)
Thanks Michael > GRASS will continue to work with bash on platforms that support > it--primarily Linux and Unix. This gives me a sigh of relief at the moment > However, we are using Python for the GUI now and Python has a number of > advantages over bash for scripting. It is a full-featured programming > language and hence more powerful than bash. It is also very cross-platform > and will run on more OS platforms than bash. So we'd like to encourage users > to begin writing scripts in Python. I learned shell scripting primarily because of GRASS. Probably now is the time to start learning Python (for GRASS, QGIS and R). A beginners tutorial would be an excellent material to lessen the barrier. Just like Chapter 9.2 of the GRASS book. -- 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 going to python scripting (was Re: [GRASS-user] new module: r.surf.volcano)
another possible stupid questions: Am I right that GRASS 7 will use Python as the scripting language, therefore, unix shell scripting won't be supported anymore? Should I start learning Python now? Should I port my scripts to Python? Will a new GRASS book include Python script samples? On Tue, Apr 21, 2009 at 4:36 AM, Hamish wrote: > > Dylan Beaudette wrote: >> Possible stupid question: do we have a standardized Python >> script template, something like the old shell-style template? > > Possible stupid answer: because unix scripting and python both > use '#' as the comment character, 'g.module --script' will create > a script template for both exactly the same. just change the > shebang line at the top from /bin/sh to /bin/python or whatever > that has to be and also replace the "are we already in grass?" > shell-style check at the end. > > and of course reuse from the other python scripts in (mostly) > grass7's scripts/ dir as a good guide for setting up temp files > and other common but subtle tasks. > > > Hamish > > > > > > ___ > 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: r.neighbors modal filter for two raster layers
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 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 > 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 >> 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 >>> 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] Re: r.neighbors modal filter for two raster layers
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 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 > 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 >> 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
Re: [GRASS-user] convert r.report output to spreadsheet format
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 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] convert r.report output to spreadsheet format
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
[GRASS-user] sample bash script for multiple mapsets
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] Re: r.neighbors modal filter for two raster layers
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 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 > 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] Re: r.neighbors modal filter for two raster layers
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 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] r.neighbors modal filter for two raster layers
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] Announcing the release of QGIS 1.0 'Kore'
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
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 wrote: > On Wed, Jan 7, 2009 at 4:04 AM, maning sambale > 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
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 Date: Wed, Jan 7, 2009 at 11:04 AM Subject: mosaicing classified images To: grass mailing list 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"?
> 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
Re: [GRASS-user] Why v.extract produces columns of type "CHARACTER"?
On Wed, Jan 14, 2009 at 5:02 PM, Glynn Clements 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 > ___ > 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: mosaicing classified images
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 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
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] LANDSAT data archive now free
Here: http://glovis.usgs.gov/ or http://edcsns17.cr.usgs.gov/EarthExplorer/ On Fri, Dec 12, 2008 at 2:11 PM, Rainer M Krug wrote: > On Fri, Dec 12, 2008 at 7:38 AM, maning sambale > 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] LANDSAT data archive now free
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] osm data in grass
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 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 *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
Re: [GRASS-user] Re: combining and evaluating vector dbf values
Geez!!! 2 days of banging my head I just did with sqlite databrowser in 15 minutes (including create new mapset, download install sqlite, and run concat to the database) Thanks again! And sqlite will be my default db driver for now. cheers, maning On Wed, Nov 26, 2008 at 6:12 PM, Markus Neteler <[EMAIL PROTECTED]> wrote: > On Wed, Nov 26, 2008 at 10:45 AM, maning sambale > <[EMAIL PROTECTED]> wrote: >> Thanks for bringing some light into this dark FOSS GIS moments. > > :) > Note that SQLite is the default driver in GRASS 7, so those dark moments > shouldn't happen any more in future (mhh, maybe there will be others...). > > Markus > >> On Wed, Nov 26, 2008 at 5:42 PM, Markus Neteler <[EMAIL PROTECTED]> wrote: >>> On Wed, Nov 26, 2008 at 10:40 AM, maning sambale >>> <[EMAIL PROTECTED]> wrote: >>>> So, it's really a limitation inherent to the dbf driver? >>> >>> Yes. DBF != SQL... we hacked in a couple of SQL statements but not entire >>> SQL92. >>> >>>> I guess sqlite it is. Is it possible to g.copy again the vector from the >>>> sqlite Mapset to the dbf Mapset? >>> >>> Sure! >>> >>> Best >>> Markus >>> >>>> cheers, >>>> maning >>>> >>>> >>>> On Wed, Nov 26, 2008 at 5:36 PM, Markus Neteler <[EMAIL PROTECTED]> wrote: >>>>> On Tue, Nov 25, 2008 at 9:07 AM, maning sambale >>>>> <[EMAIL PROTECTED]> wrote: >>>>>> 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 >>>>> >>>>> Did you try the SQLite driver? >>>>> >>>>> Just create a new mapset, define the connection with "db.connect" (see >>>>> manual >>>>> for copy-paste command) and "g.copy" the map from the DBF mapset >>>>> into the sqlite mapset. >>>>> Then the joys of SQL are available. >>>>> >>>>> Markus >>>>> >>>> >>>> >>>> >>>> -- >>>> |-|--| >>>> | __.-._ |"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/ | >>>> |-|--| >>>> >>> >>> >>> >>> -- >>> Open Source Geospatial Foundation >>> http://www.osgeo.org/ >>> http://www.grassbook.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/ | >> | _)_/L I http://epsg4253.wordpress.com/ | >> |-|--| >> > > > > -- > Open Source Geospatial Foundation > http://www.osgeo.org/ > http://www.grassbook.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/ | | _)_/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
After struggling with this without any solution for the past 2 dyas. I will try r.cross On Tue, Nov 25, 2008 at 4:07 PM, maning sambale <[EMAIL PROTECTED]> wrote: > 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/ | > |-|--| > -- |-|--| | __.-._ |"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
