Re: [GRASS-user] using r.resamp.stats- When I reach a tie

2011-02-28 Thread António Rocha

Hi Glynn.
Thanks for your answer. Is there any other alternative, to resample 
data, that uses another approach to solve tie events?

THanks
Antonio

Glynn Clements wrote:

António Rocha wrote:

  
I'm planning on using r.resamp.stats in order to  resample a grid to a 
coarser resolution. Since I'm using this function with classes I will 
use mode method (more frequentclass to fall into a destination cell is 
the chosen). My question is: when I have a tie between two classes, how 
does GRASS decides which class shall be chosen?



In the event of a tie, the lowest-numbered category is chosen.

  




__ Information from ESET NOD32 Antivirus, version of virus signature 
database 5912 (20110227) __

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


Re: [GRASS-user] New User

2011-02-28 Thread Saber Razmjooei
If you are intending to use GIS for visualization and data conversion,
maybe better to use something easy to start off, like QGIS:
http://www.qgis.org/wiki/Download#Release_2

You have GRASS plugin there if needed. But usually you should be able to
do export the files to other OGR/GDAL supported formats from within
QGIS.

Regards
Sab

On Sun, 2011-02-27 at 23:25 -0800, Hamish wrote:
 Thomas wrote:
  I am a new user, just discovered grass and OSGeo. My current
  interest is accessing the files from NaturalEarthData,
  converting them into matlab format, and, using octave to
  manipulate and/or plot.
  
  I use Linux and Unix. Octave 3.4.0, Gnuplot, Grass
  6.4.0.  This workstation has Ubuntu 10.04.
  
  I downloaded the
  10m-admin-0-countries.zip   and
  10m-admin-0-boundary-lines-land.zip files from
  NaturalEarthData.
  
  I want to use grass to display some portion of the data.
  
  Then, use r.out.mat to generate a octave compatible binary
  file.
  
  I started looking at the manual, but, quickly became
  confused.
 
 that would be my fault.
 
  Is there a simple how to for what I want to do?
 
 Could you explain your workflow a bit more? there are many
 choices.   In general it's not so bad, and you can do matlab/oct
 to and from grass without much hassle.
 
 natural earth is polyline data, so vector  thus v.* modules.
 2D array is raster data, so r.* modules.
 see v.to.rast and r.to.vect modules.
 
 
 do you have 2D arrays in matlab/octave from contourf() or similar
 and you want to fill(), area(), or patch() a polygon coastline
 over the top?
 see  http://grass.osgeo.org/wiki/Addons#v.out.ascii.mat
 
 are you using the (free/gratis) m_map toolbox?
   http://www.eos.ubc.ca/~rich/map.html
 or matlab's mapping toolbox? (which is using proj4, same as us)
 
 do you wish to rasterize the natural earth polygons into a 2D
 array?
 see v.to.rast + r.out.mat
 
 do you want to bring your matlab data into grass and plot it
 there?
 see r.in.mat, or v.in.lines, v.in.mapgen.
 
 
 the first step is to bring the data into grass with v.in.ogr.
 
 
 Hamish
 
 
 
   
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user


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


[GRASS-user] Re: How to merge rasters using feathering/blending ?

2011-02-28 Thread kaipi


 You can create a buffer zone between your 2 dataset where they merge and
 set the buffer to null values. Then using r.fillnulls 
 

Thanks. This works fine for elevation data but not for satellite imagery.



 g.region rast=map1,map2
 r.series in=map1,map2 out=both_maps method=average 
 

Thanks. But this is not really a smooth blending/fading/feathering ? Any
alternatives ?

Kaipi


-- 
View this message in context: 
http://osgeo-org.1803224.n2.nabble.com/How-to-merge-rasters-using-feathering-blending-tp6061489p6073011.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


Re: [GRASS-user] Re: How to merge rasters using feathering/blending ?

2011-02-28 Thread Saber Razmjooei
If you are working with Images, maybe try to use GIMP. It has stitching
and smoothing/blending options. I use them quite often for my photos.

Regards
Sab

On Mon, 2011-02-28 at 04:05 -0800, kaipi wrote:
 
  You can create a buffer zone between your 2 dataset where they merge and
  set the buffer to null values. Then using r.fillnulls 
  
 
 Thanks. This works fine for elevation data but not for satellite imagery.
 
 
 
  g.region rast=map1,map2
  r.series in=map1,map2 out=both_maps method=average 
  
 
 Thanks. But this is not really a smooth blending/fading/feathering ? Any
 alternatives ?
 
 Kaipi
 
 


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


[GRASS-user] Re: How to merge rasters using feathering/blending ?

2011-02-28 Thread kaipi
Gimp is fine if you have a smal number of scenes but I need to process large
amounts of scenes where I need an automated way to feather the seam lines.

-- 
View this message in context: 
http://osgeo-org.1803224.n2.nabble.com/How-to-merge-rasters-using-feathering-blending-tp6061489p6073453.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


Re: [GRASS-user] calculating average value over a base map

2011-02-28 Thread Christian Schwartze
Monica,

try to use r.statistics with the average function/method. Pass you maps through 
the cover= and base= parameters...

Christian.

From: Monica Buescu 
Sent: Monday, February 28, 2011 4:27 PM
To: GRASS user list 
Subject: [GRASS-user] calculating average value over a base map

Greetings 

I have a map with values and a base map created with r.clump. How can I 
calculate the average of values that falls inside each cluster defined by 
r.clump?
Thanks

Best regards,
Monica



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


Re: [GRASS-user] calculating average value over a base map

2011-02-28 Thread Monica Buescu
Greetings Christian

Unfortunely it only works with CELL values. For instance, r.median uses more
than just CELL.
Is there any other alternative? Beucase for mode and median there is
specific functions.

Thamls

Best regards,
Monica

On Mon, Feb 28, 2011 at 3:55 PM, Christian Schwartze 
christian.schwar...@uni-jena.de wrote:

   Monica,

 try to use r.statistics with the average function/method. Pass you maps
 through the cover= and base= parameters...

 Christian.

  *From:* Monica Buescu monicabuescu1...@gmail.com
 *Sent:* Monday, February 28, 2011 4:27 PM
 *To:* GRASS user list grass-user@lists.osgeo.org
 *Subject:* [GRASS-user] calculating average value over a base map

 Greetings

 I have a map with values and a base map created with r.clump. How can I
 calculate the average of values that falls inside each cluster defined by
 r.clump?
 Thanks

 Best regards,
 Monica

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


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


Re: [GRASS-user] calculating average value over a base map

2011-02-28 Thread Christian Schwartze
If you want to deal with raster value types other than CELL, you can use the –c 
flag in r.statistics in order to process the labels. Therfore, you have to copy 
your cat values to the category labels first...

Regards,
Christian.
From: Monica Buescu 
Sent: Monday, February 28, 2011 5:16 PM
To: Christian Schwartze 
Cc: GRASS user list 
Subject: Re: [GRASS-user] calculating average value over a base map

Greetings Christian 

Unfortunely it only works with CELL values. For instance, r.median uses more 
than just CELL.
Is there any other alternative? Beucase for mode and median there is specific 
functions.

Thamls

Best regards,
Monica

On Mon, Feb 28, 2011 at 3:55 PM, Christian Schwartze 
christian.schwar...@uni-jena.de wrote:

  Monica,

  try to use r.statistics with the average function/method. Pass you maps 
through the cover= and base= parameters...

  Christian.

  From: Monica Buescu 
  Sent: Monday, February 28, 2011 4:27 PM
  To: GRASS user list 
  Subject: [GRASS-user] calculating average value over a base map

  Greetings 

  I have a map with values and a base map created with r.clump. How can I 
calculate the average of values that falls inside each cluster defined by 
r.clump?
  Thanks

  Best regards,
  Monica

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

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


Re: [GRASS-user] calculating average value over a base map

2011-02-28 Thread Monica Buescu
Hello Christian.
I believe I didn't quite understand what you have said but I think that is
not the way.
My basemap is a r.clump output in CELL format that ranges from 1 til 1700 or
something
my cover map is NDVI that ranges from -1 (near) til 1 (near) and is a FCELL.

My Objective. calculate for each clump, the average of NDVI pixels

If I do:
r.statistics -c base=new_clump cover=2007_ndvi method=average

I get:
ERROR: This module currently only works for integer (CELL) maps

Is there any alternative to do this?

Thank you Christian and the rest of grass mailing list

Monica

On Mon, Feb 28, 2011 at 5:08 PM, Christian Schwartze 
christian.schwar...@uni-jena.de wrote:

   If you want to deal with raster value types other than CELL, you can use
 the –c flag in r.statistics in order to process the labels. Therfore, you
 have to copy your cat values to the category labels first...

 Regards,
 Christian.
   *From:* Monica Buescu monicabuescu1...@gmail.com
 *Sent:* Monday, February 28, 2011 5:16 PM
 *To:* Christian Schwartze christian.schwar...@uni-jena.de
 *Cc:* GRASS user list grass-user@lists.osgeo.org
 *Subject:* Re: [GRASS-user] calculating average value over a base map

 Greetings Christian

 Unfortunely it only works with CELL values. For instance, r.median uses
 more than just CELL.
 Is there any other alternative? Beucase for mode and median there is
 specific functions.

 Thamls

 Best regards,
 Monica

 On Mon, Feb 28, 2011 at 3:55 PM, Christian Schwartze 
 christian.schwar...@uni-jena.de wrote:

   Monica,

 try to use r.statistics with the average function/method. Pass you maps
 through the cover= and base= parameters...

 Christian.

  *From:* Monica Buescu monicabuescu1...@gmail.com
 *Sent:* Monday, February 28, 2011 4:27 PM
 *To:* GRASS user list grass-user@lists.osgeo.org
 *Subject:* [GRASS-user] calculating average value over a base map

   Greetings

 I have a map with values and a base map created with r.clump. How can I
 calculate the average of values that falls inside each cluster defined by
 r.clump?
 Thanks

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



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


[GRASS-user] Question on Importing E$RI Grid DEM

2011-02-28 Thread Rich Shepard

  I thought I had a saved thread from this list on this subject, but
grepping my mail subdirectory I don't find it. Here's the situation:

  A set of 30m NED (National Elevation Database) DEMs as Lat/Long 1 degree
blocks. Data is in ARC GRID file format. Each directory has three
subdirectores: demgrid/, info/, and metadata/. The demgrid/ subdirectoy
contains: dblbnd.adf, hdr.adf, log, prj.adf, sta.adf, w001001.adf, and
w001001x.adf files.

  I thought I wanted to point the grass65 gui for raster-arc grid import to
the hdr.adf file, but I get errors and no import when I do that. Is there
another file that should be specified for proper import?

TIA,

Rich

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


Re: [GRASS-user] New User

2011-02-28 Thread Thomas D. Dean
On Sun, 2011-02-27 at 23:25 -0800, Hamish wrote:

I know little or nothing about mapping.  I want to understand how the
data files and grass connect together.  Some tutorial using known data
followed by some information on various data formats and how to convert
them into something GRASS can handle

 Could you explain your workflow a bit more? there are many
 choices.   In general it's not so bad, and you can do matlab/oct
 to and from grass without much hassle.

I have data consisting of bearings amd/or distances from known points,
in octave.  I want to plot the bearing lines and the distance arcs on a
wgs84 ellipsoid with national boundaries (outlines).  The end result is
to visually locate the intersections of the various lines with national
outlines for easier identification of RF sources.

 
 natural earth is polyline data, so vector  thus v.* modules.
 2D array is raster data, so r.* modules.
 see v.to.rast and r.to.vect modules.

The data from natural earth is
10m_admin_0_countries.dbf - spreadsheet database file
10m_admin_0_countries.prj - projection
10m_admin_0_countries.shp - shape file
10m_admin_0_countries.shx 

I also have WDB data from
http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO5/BOUNDARY/WDB/

BOUNDS.DAT
COASTS.DAT
ETOPO5.DAT
ETOPO5.DOS
RIVERS.DAT

I have not discovered how to bring the map data into grass.  How do I
arrange the data files so GRASS can find/inport them?  I have a
$home/grassdata directory structure with worldlocation.tar.gz expanded
into it as $home/grassdata/worldlocation.  How do I access/plot this
data in GRASS?

How do I put the ETOPO5 data in the grassdata tree?  I think I need to
convert it to some format that GRASS understands.

How do I put the 10m_admin_0_countries.* files into the grassdata tree?

 
 do you have 2D arrays in matlab/octave from contourf() or similar
 and you want to fill(), area(), or patch() a polygon coastline
 over the top?

I have these in octave 3.4.0.

 see  http://grass.osgeo.org/wiki/Addons#v.out.ascii.mat
 
 are you using the (free/gratis) m_map toolbox?
   http://www.eos.ubc.ca/~rich/map.html
 or matlab's mapping toolbox? (which is using proj4, same as us)
 
 do you wish to rasterize the natural earth polygons into a 2D
 array?
 see v.to.rast + r.out.mat
 
 do you want to bring your matlab data into grass and plot it
 there?
 see r.in.mat, or v.in.lines, v.in.mapgen.

This would be Ok, assuming I can export the resulting map in pdf format.
My first thought was to use GRASS to generate a MAT file from map data,
to load into octave.

 
 
 the first step is to bring the data into grass with v.in.ogr.

I can start grass and get the startup screen.  I select worlddata and
PERMANENT and click Enter GRASS.  What is next?

tomdean

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


Re: [GRASS-user] Question on Importing E$RI Grid DEM -- UPDATE

2011-02-28 Thread Rich Shepard

On Mon, 28 Feb 2011, Rich Shepard wrote:


 I thought I wanted to point the grass65 gui for raster-arc grid import
to the hdr.adf file, but I get errors and no import when I do that. Is
there another file that should be specified for proper import?


  How strange. It works from the command line:

GRASS 6.5.svn (Nevada-ll):~/grassdata  r.in.gdal -e --overwrite
input=/home/rshepard/GIS/data/Nevada/DEM/area2/demgrid/hdr.adf
output=dem_area2 title='DEM NED Area 2'
Projection of input dataset and current location appear to match
 100%
r.in.gdal complete. Raster map dem_area2 created.

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


Re: [GRASS-user] Question on Importing E$RI Grid DEM -- UPDATE

2011-02-28 Thread Martin Landa
2011/2/28 Rich Shepard rshep...@appl-ecosys.com:

Just a note: please keep the subject otherwise thread is broken and
it's hard to track it...

Martin

-- 
Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


RE: [GRASS-user] Re: How to merge rasters using feathering/blending ?

2011-02-28 Thread Moskovitz, Bob
I haven't tried it yet, but what about this:
http://grass.osgeo.org/wiki/GRASS_AddOns#r.smoothpatch

There are other useful addons, so you might want to browse around that
page.

Bob

Bob Moskovitz
Seismic Hazards Mapping Program
California Geological Survey
801 K Street   MS-1231
Sacramento, CA  95814
ph. (916) 323-0043

CONFIDENTIALITY NOTICE: This communication is intended only for the use
of the individual or entity to which it is addressed. This message
contains information from the State of California, California Geological
Survey, which may be privileged, confidential and exempt from disclosure
under applicable law, including the Electronic Communications Privacy
Act. If the reader of this communication is not the intended recipient,
you are hereby notified that any dissemination, distribution, or copying
of this communication is strictly prohibited.



 -Original Message-
 From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-
 boun...@lists.osgeo.org] On Behalf Of kaipi
 Sent: Monday, February 28, 2011 6:27 AM
 To: grass-user@lists.osgeo.org
 Subject: [GRASS-user] Re: How to merge rasters using
 feathering/blending ?
 
 Gimp is fine if you have a smal number of scenes but I need to process
 large
 amounts of scenes where I need an automated way to feather the seam
 lines.
 
 --
 View this message in context: http://osgeo-
 org.1803224.n2.nabble.com/How-to-merge-rasters-using-feathering-
 blending-tp6061489p6073453.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


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


Re: [GRASS-user] New User

2011-02-28 Thread Rich Shepard

On Mon, 28 Feb 2011, Thomas D. Dean wrote:


I know little or nothing about mapping.  I want to understand how the
data files and grass connect together.  Some tutorial using known data
followed by some information on various data formats and how to convert
them into something GRASS can handle


Thomas,

  First, there is a tutorial of sorts using the Spearfish, ND data sets.
Second, GRASS is a tool for spatial analyses and modeling, and as such
requires users to know how to properly use the tool. I do not intend to
sound condescending, but I highly recommend you read a book or information
found on the Web about GIS, data formats, spatial analyses, and related
subjects.


In general it's not so bad, and you can do matlab/oct to and from grass
without much hassle.


  MatLab (and Octave) are totally different tools. I've used Octave to
calculate eigenvectors of symetrical matrices, but would not use it for any
spatial analyses.


I have data consisting of bearings amd/or distances from known points, in
octave.  I want to plot the bearing lines and the distance arcs on a wgs84
ellipsoid with national boundaries (outlines).  The end result is to
visually locate the intersections of the various lines with national
outlines for easier identification of RF sources.


  This GRASS can do.


The data from natural earth is 10m_admin_0_countries.dbf - spreadsheet
database file 10m_admin_0_countries.prj - projection
10m_admin_0_countries.shp - shape file 10m_admin_0_countries.shx


  What you want is File - Import vector data - top menu choice for ESRI
.shp files. Select the .shp as the file to open. Or, use v.in.ogr from the
command line.


I have not discovered how to bring the map data into grass.  How do I
arrange the data files so GRASS can find/inport them?  I have a
$home/grassdata directory structure with worldlocation.tar.gz expanded
into it as $home/grassdata/worldlocation.  How do I access/plot this
data in GRASS?


  Put your shape files in there, too. Read the quick introduction to GRASS
on the Web site as it explains the relationships among the root database,
locations, and mapsets. It's all there and very well written.

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


Re: [GRASS-user] Question on Importing E$RI Grid DEM -- UPDATE

2011-02-28 Thread Martin Landa
Hi,

2011/2/28 Rich Shepard rshep...@appl-ecosys.com:
 On Mon, 28 Feb 2011, Martin Landa wrote:

 Just a note: please keep the subject otherwise thread is broken and
 it's hard to track it...

 Martin,

  Does the postfix affect the thread? I thought it did not which is why many
 of us add FIXED or SOLVED to close threads.

ops, it seems that postfix doesn't affect the thread [1]. Anyway
thread is broken at least in gmail web client. So sorry for the noise.

Martin

[1] http://lists.osgeo.org/pipermail/grass-user/2011-February/thread.html

-- 
Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] 6.4svn graphical modeller instead of vector or raster add

2011-02-28 Thread Stephen Sefick
I can't add raster or vector layer to grass 6.4svn because when I click
on the GUI for adding layers the graphic modeller comes up.  What should
I do?

Stephen Sefick

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


Re: [GRASS-user] 6.4svn graphical modeller instead of vector or raster add

2011-02-28 Thread Martin Landa
Hi,

2011/2/28 Stephen Sefick ssef...@gmail.com:
 I can't add raster or vector layer to grass 6.4svn because when I click
 on the GUI for adding layers the graphic modeller comes up.  What should
 I do?

fresh SVN installation? OS?

Martin

-- 
Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] 6.4svn graphical modeller instead of vector or raster add

2011-02-28 Thread Stephen Sefick
Yes, fresh svn installation, and Ubuntu 10.04.

On Mon, 2011-02-28 at 20:22 +0100, Martin Landa wrote:
 Hi,
 
 2011/2/28 Stephen Sefick ssef...@gmail.com:
  I can't add raster or vector layer to grass 6.4svn because when I click
  on the GUI for adding layers the graphic modeller comes up.  What should
  I do?
 
 fresh SVN installation? OS?
 
 Martin
 


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


Re: [GRASS-user] Re: How to merge rasters using feathering/blending ?

2011-02-28 Thread Michael Perdue
In GRASS, one thing you can try is to normalize your data so the variations 
from scene to scene are minimized. ie, run a linear regression between 
overlapping maps and then use the coefficients to normalize one of the data 
sets. However, it is tedious and doesn't work well if the two images arn't well 
correlated to begin with.

A second option is to switch to OSSIM (another open source project) for the 
mosaicing task. It has a number of tools for feathering, blending and histogram 
matching. It's all fairly automated and straight forward to use.

A third option is enblend. It uses multi-resolution splines to blend images 
together. From my own experience, it does the best job at blending Ortho's 
together, but it's not geographically aware, which means that you need to do 
some hacking to develop a workflow that you can use to maintain/re-establish a 
spatial reference system.

--- On Mon, 2/28/11, kaipi mapcoll...@gmx.net wrote:

From: kaipi mapcoll...@gmx.net
Subject: [GRASS-user] Re: How to merge rasters using feathering/blending ?
To: grass-user@lists.osgeo.org
Received: Monday, February 28, 2011, 9:27 AM

Gimp is fine if you have a smal number of scenes but I need to process large
amounts of scenes where I need an automated way to feather the seam lines.

-- 
View this message in context: 
http://osgeo-org.1803224.n2.nabble.com/How-to-merge-rasters-using-feathering-blending-tp6061489p6073453.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


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


Re: [GRASS-user] How to merge rasters using feathering/blending ?

2011-02-28 Thread Jamie Adams
I wrote a script to do something like this at one point, though I can't find
it now.  It was intended for terrain, but could be used per band on imagery.
 It used r.buffer to create a distance gradient and then r.recode to convert
this gradient to a percentage.  I then used r.mapcalc to merge the 2 rasters
based on the percentage.

The script automated the buffer distances and the creation of the recode.txt
file, but went something like this:

## make a null raster
r.mapcalc temp_null = if(isnull(source),1,null)

## create a distance gradient from the null raster.  The output will have
int values 1-7 where 1 equals the input null raster
r.buffer in=temp_null out=temp_buf distances=10,20,30,40,50,60 units=meters

## fill nulls in the distance gradient with 0 for recoding
r.mapcalc temp_buf_fill=if(isnull(temp_buf),0,temp_buf)

## recode to a percentage
r.recode in=temp_buf_fill out=temp_recode rules=recode.txt

where recode.txt would have the range from temp_buf_fill expressed as
percentages:
1:1:100:100
2:2:86:86
3:3:71:71
4:4:57:57
5:5:43:43
6:6:29:29
7:7:14:14
0:0:0:0

## merge the output
r.mapcalc
merged=if(isnull(source),patch,(patch*(temp_recode/100.0))+(source*(1.0-(temp_recode/100.0

Sorry if that doesn't work, but hopefully you get the idea.

- Jamie

On Thu, Feb 24, 2011 at 10:53 AM, kaipi mapcoll...@gmx.net wrote:


 I need to patch rasters containing NULL areas with alternative rasters. I
 can
 do this with r.patch but due to different quality/resolution of the rasters
 this produces prominent seamlines. Is there a way to feather/blend between
 merged rasters using Grass GIS?

 Thanks,

 Kaipi
 --
 View this message in context:
 http://osgeo-org.1803224.n2.nabble.com/How-to-merge-rasters-using-feathering-blending-tp6061489p6061489.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

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


[GRASS-user] calculate raster statistics for vector polygons

2011-02-28 Thread Kim Besson
Greetings I have a vector file with polygons and I want to obtain a raster
map with statistics  of an input raster for each polygon. What can I use to
do this?

Thanks

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


Re: [GRASS-user] calculate raster statistics for vector polygons

2011-02-28 Thread Saber Razmjooei
http://grass.fbk.eu/grass64/manuals/html64_user/v.rast.stats

If you have lots of polygons, you probably need to look into 
v.rast.stats2 , an addon

Cheers
Sab

On Mon, 2011-02-28 at 21:32 +, Kim Besson wrote:
 Greetings I have a vector file with polygons and I want to obtain a
 raster map with statistics  of an input raster for each polygon. What
 can I use to do this?
 
 
 Thanks
 
 
 Kim
 ___
 grass-user mailing list
 grass-user@lists.osgeo.org
 http://lists.osgeo.org/mailman/listinfo/grass-user


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


Re: [GRASS-user] 6.4svn graphical modeller instead of vector or raster add

2011-02-28 Thread Martin Landa
[back to ML]

2011/2/28 Stephen Sefick ssef...@gmail.com:
 Yes whenever I click Add raster map layer I get a modeller window.

it's a bug in silk icon set. Martin

-- 
Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] 6.4svn graphical modeller instead of vector or raster add

2011-02-28 Thread stephen sefick
How do I get rid of the silk icon set?  Or how do I not use it in the
build for both trunk and 6.4?
Thanks,

Stephen

On Mon, Feb 28, 2011 at 3:48 PM, Martin Landa landa.mar...@gmail.com wrote:
 [back to ML]

 2011/2/28 Stephen Sefick ssef...@gmail.com:
 Yes whenever I click Add raster map layer I get a modeller window.

 it's a bug in silk icon set. Martin

 --
 Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa




-- 
Stephen Sefick

| Auburn University                                         |
| 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


Re: [GRASS-user] 6.4svn graphical modeller instead of vector or raster add

2011-02-28 Thread Martin Landa
2011/2/28 stephen sefick ssef...@gmail.com:
 How do I get rid of the silk icon set?  Or how do I not use it in the
 build for both trunk and 6.4?

from Layer Manager menu - Settings | Preferences | Advanced - set Icon
Theme to 'grass2', save settings and restart GUI.

Martin

-- 
Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] 6.4svn graphical modeller instead of vector or raster add

2011-02-28 Thread Martin Landa
Hi,

2011/2/28 Martin Landa landa.mar...@gmail.com:
 2011/2/28 stephen sefick ssef...@gmail.com:
 How do I get rid of the silk icon set?  Or how do I not use it in the
 build for both trunk and 6.4?

 from Layer Manager menu - Settings | Preferences | Advanced - set Icon
 Theme to 'grass2', save settings and restart GUI.

just not that trunk contains only one icon theme 'grass' which is
based on 'grass2' from GRASS 6.

Martin

-- 
Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] v.rast.stats2 bug?

2011-02-28 Thread Mark Seibel
I just checked out v.rast.stats2 (Feb.28) and when I run it I get the
following error.

Do I need to file a ticket, or am I missing something obvious?  While
it says 'done', the table is unpopulated.


awk: {printf \nUPDATE acoe_jd SET dist_n
awk: ^ unterminated string
sed: couldn't write 176 items to stdout: Broken pipe
Restoring previous MASK...
Statistics calculated from raster map streams_distance_ft@PERMANENT
and uploaded to attribute table of vector map acoe_jd.
Done.


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


[GRASS-user] Tutorial Not Same as GRASS 6.4

2011-02-28 Thread Thomas D. Dean
I have GRASS 6.4.0RC5+39438 (2009).  I followed the tutorial, with lots
of mental corrections and managed to get the example displayed.  The
tutorial is good, just needs some tweeking to match the GRASS display.

A quick introduction to the 6.4 wxGUI and raster maps,
http://grass.osgeo.org/wiki/Quick_wx_tutorial,
Does not correspond to the GRASS display.

For example,

 Once inside add a raster map layer such as elevation from the
  PERMANENT mapset. To do this click on the checkerboard toolbar
 button with a + on it and then select the map name from the
 pull-down list which appears, and click [Ok].

There is no icon checkerboard with a + on it.  Rather, there is a
green, red, orange icon with an arrow in it.  However, it does have the
description Add raster layer.  And, there is not a pull down list.
Rather, a layer is added, labeled raster1 and clicking the i button
brings up a pop-up that allows selecting the map.

I am not good at describing things.  I don't know what the author
intended the icons to be.

How do I help feeding back the corrections?  Or, is this tutorial for a
newer version?

tomdean

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


[GRASS-user] Points, Distances, and, Bearings.

2011-02-28 Thread Thomas D. Dean
I can get the Natural Earth data to display.  Neat.

How do I enter a point, a bearing from that point, and a distance from
that point, plot the point, the bearing line, and, the distance arc?

tomdean

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


Re: [GRASS-user] Points, Distances, and, Bearings.

2011-02-28 Thread Hamish
Thomas wrote:
 I can get the Natural Earth data to
 display.  Neat.
 
 How do I enter a point, a bearing from that point, and a
 distance from that point, plot the point, the bearing line,
 and, the distance arc?

check out the d.rhumbline, d.geodesic, r.transect, r.profile,
and m.cogo modules.

also install the proj-bin package (assuming linux) and you
get the geod program which will calculate the start and end
points for you using a great-circle calculation.

are you using a lat/lon location? how large an area will your
plot cover? continental?


Hamish




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


Re: [GRASS-user] Points, Distances, and, Bearings.

2011-02-28 Thread Thomas D. Dean
On Mon, 2011-02-28 at 19:36 -0800, Hamish wrote:

 are you using a lat/lon location? how large an area will your
 plot cover? continental?
 
Thanks for the reply.  I will be away for a couple days, then, will
check out the d.rhumbline, d.geodesic, r.transect, r.profile,
and m.cogo modules.

The points are located by lat/lon, wgs84.  I also want a lat/lon output
that is the centroid of the intersection of several lines.

The plot will cover up to 1/4 of the globe.  Then, zoom down to a
country.

tomdean

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