[GRASS-user] Georrectify and nearest neighbor

2010-06-02 Thread Jenny Turner
My name is Jenny Turner and I've been using  georrectify graphical module in
WxpyTHon GRASS6.4. I saw a video of the proceeding but I would like to know
if there is a description or a tutorial of what to do and how. I have
understand what I have to do but I need to explain to other persons
(non-grass users :) ).

Just one more thing: Is it used any resampling to the data or only nearest
neighbor?

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


Re: [GRASS-user] calculating road density

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

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


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


[GRASS-user] Using v.distance between two identical layers...

2010-06-02 Thread devaux

Hi,
do you know if it's possible to use v.distance command when for and  
to layers are

the same layer ?
I have a layer describing ecological areas of interest and I want to  
qualify spatial
relation between all patches. I want to identify and calculate  
shortest relation between

all my patches boundaries.

How can I do as v.distance command doesn't have any option like :  
eliminate identical

objects from the from and to layer ?

I'd like to use this particular command because of the option that  
create lines between

nearest objects.

Thanks

Nicolas


This message was sent using IMP, the Internet Messaging Program.


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


Re: [GRASS-user] modis data

2010-06-02 Thread Benjamin Ducke
On 06/01/2010 09:10 PM, Hamish wrote:
 bharath s wrote:
 iam working on modis data (mydo9q1 250 m) using
 grass (gis) software.  i   want to overlay google earth
 boundary for southindia..  after geocorrection on modis data...  
 but the boundaries are not matching. please help me to solve
 this problem. is the projecton playing any role in it ?
 actually iam  using lat/long  and datum as wgs84.

 I would try to use the natural earth vector boundary instead of
 Google's one.

Yes, indeed. Since Google's shiny 3D toy is getting dangerously close to being 
accepted into mainstream science and research now, one should be aware that 
there are differences between Google's private way of handling lat/lon 
projections and  international geodetic systems:

http://www.sharpgis.net/post/2008/05/SphericalWeb-Mercator-EPSG-code-3785.aspx

This will lead to inaccuracies when trying to overlay Google data and EPSG 4326 
data.
The same problem exists for Microsoft's clone of Google Earth (they do seem to 
take
the copying business serious at MS). The magnitude of error will vary with your
region's scale and location.

NASA's 3D Earth tool is, as may be expected, more competent, so perhaps the
better choice for 3D GIS visualization:

http://worldwindcentral.com/wiki/Map_projection_Plug-in

Best,

Ben


 see 1:10 million, 1:50 million and 1:110million scale maps from
 http://grass.osgeo.org/wiki/Global_datasets#Natural_Earth_imagery

 if 1:10m is too crude for you, there NOAA's coastline extractor
 + v.in.mapgen, or v.in.gshhs from wiki addons; or others listed
 on the above wiki page.

 see also  http://grass.osgeo.org/wiki/MODIS


 Hamish



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





--
Files attached to this email may be in ISO 26300 format (OASIS Open Document 
Format). If you have difficulty opening them, please visit http://iso26300.info 
for more information.

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


[GRASS-user] Using Vector Digitizer in Windows

2010-06-02 Thread Kim Besson
Hi

I'm using Vector Digitizer in MS-Windows (Tcltk since WxGUI is not
available) and I wanted to know How can I put an image as a background? I am
at v.digit form and I cannot see any field regarding that
Thanks
Kim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] get coords from a raster

2010-06-02 Thread Margherita Di Leo
Hi List,

I'm trying to make a script to calculate the slope along a path. I want to
calculate the slope as the average of cell per cell slope. For my purpose I
found r.profile to produce the input file for my script, but I have to give
it the coordinates of the points. As I have a raster and not a vector, I was
wondering if is there a non interactive way like a d.where or something like
that, that i can use for my purpose. I mean, i have a raster in which there
is the path i'm interested in, and the other cells are null.

Thank you in advance

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


Re: [GRASS-user] get coords from a raster

2010-06-02 Thread Margherita Di Leo
Christian,

thank you for answering. My problem is that I don't want to use the slope
map to get slope, bu i want to calculate it as the
ratio between the difference of the elevation and the distance between two
given points (pixels). To do that, I need to get the coordinates of the
points as a input file for r.profile.

Margherita


2010/6/2 Christian Schwartze christian.schwar...@uni-jena.de

  Why are you not using r.statistics, since it takes one base map (your
 path as a raster) and another map containing data you want to analyze (your
 slope). With the average option/method you will get the desired result, if
 I understand you correctly...

 Christian.

  *From:* Margherita Di Leo direg...@gmail.com
 *Sent:* Wednesday, June 02, 2010 6:17 PM
 *To:* grass-user@lists.osgeo.org
 *Subject:* [GRASS-user] get coords from a raster

 Hi List,

 I'm trying to make a script to calculate the slope along a path. I want to
 calculate the slope as the average of cell per cell slope. For my purpose I
 found r.profile to produce the input file for my script, but I have to give
 it the coordinates of the points. As I have a raster and not a vector, I was
 wondering if is there a non interactive way like a d.where or something like
 that, that i can use for my purpose. I mean, i have a raster in which there
 is the path i'm interested in, and the other cells are null.

 Thank you in advance

 Margherita

 --

 ___
 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] huff model

2010-06-02 Thread Ricardo Rodríguez
hello all, saw some of you to know if there is any implementation of huff
model for grass gis or failing that a script on it in Open Source software
arcgis but I do not want to work with such libraries.

thanks for your answers

Ricardo Rodríguez
Univalle
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] surface cost

2010-06-02 Thread Ricardo Rodríguez
hello all, I have a question that can be evaluated as sustainable urban
accessibility of a zone map go through friction, but the only difference to
the examples in the book that the analysis I want to do is not over in many
respects an attractive location but several points of interest between
themselves, to facilitate the compression of what I want to analyze the
connection (and accessibility) among various warehouses or schools, both
topologically (factor map distance matrix) as the gravitational level ( huff
model), which commands do you recommend for this type of process.

code more than I would like to hear experiences and suggestions on this
topic

thanks for your attention, time and answers
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Using Vector Digitizer in Windows

2010-06-02 Thread Markus Neteler
On Wed, Jun 2, 2010 at 6:10 PM, Kim Besson kimbesson1...@gmail.com wrote:
 Hi
 I'm using Vector Digitizer in MS-Windows (Tcltk since WxGUI is not
 available) and I wanted to know How can I put an image as a background? I am
 at v.digit form and I cannot see any field regarding that

There should be, when starting, the option to select a background display
*command*

http://grass.osgeo.org/grass64/manuals/html64_user/v.digit.html

which could be set to
d.rast map=name_of_raster

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


Re: [GRASS-user] get coords from a raster

2010-06-02 Thread Markus Neteler
On Wed, Jun 2, 2010 at 6:50 PM, Margherita Di Leo direg...@gmail.com wrote:
 Christian,

 thank you for answering. My problem is that I don't want to use the slope
 map to get slope, bu i want to calculate it as the
 ratio between the difference of the elevation and the distance between two
 given points (pixels).

If you treat them as vectors, perhaps this could help?


v.to.db
...
 slope: slope steepness of vector line or boundary
?

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


Re: [GRASS-user] calculating road density

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

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

Perhaps this could be a fast alternative:
http://grass.osgeo.org/wiki/GRASS_AddOns#r.univar.zonal

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


Re: [GRASS-user] huff model

2010-06-02 Thread Markus Neteler
On Wed, Jun 2, 2010 at 7:21 PM, Ricardo Rodríguez
ricardorodo...@gmail.com wrote:
 hello all, saw some of you to know if there is any implementation of huff
 model for grass gis or failing that a script on it in Open Source software
 arcgis but I do not want to work with such libraries.

Do you mean something like this?

http://www.directionsmag.com/article.php?article_id=896

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


Re: [GRASS-user] get coords from a raster

2010-06-02 Thread Margherita Di Leo
Markus,

Thank you. I'm pretty sure this is the way. So I run:

r.to.vect input=mainchannel output=points feature=point --o
v.to.db map=points option=coor col=x,y (like the man page suggests)

but I get a number of errors like:

Column 'x' not found
Error in db_execute_immediate()

where am I wrong now?



2010/6/2 Markus Neteler nete...@osgeo.org


 If you treat them as vectors, perhaps this could help?


 v.to.db
 ...
 slope: slope steepness of vector line or boundary


I calculate slope as

S = (z2-z1)/(sqrt((x2-x1)^2+(y2-y1)^2))*100

where x1,y1 are the coords of point 1, x2,y2 are the coords of point 2, and
z1, z2 are the elevations in the raster map associate at the points 1 and 2
respectively.

Margherita


 ?

 Markus

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


Re: [GRASS-user] get coords from a raster

2010-06-02 Thread Markus Neteler
On Wed, Jun 2, 2010 at 10:07 PM, Margherita Di Leo direg...@gmail.com wrote:
 Markus,

 Thank you. I'm pretty sure this is the way. So I run:

 r.to.vect input=mainchannel output=points feature=point --o
 v.to.db map=points option=coor col=x,y (like the man page suggests)

 but I get a number of errors like:

 Column 'x' not found
 Error in db_execute_immediate()

 where am I wrong now?

Two times v.db.addcol is missing for the x and y columns (before v.to.db).

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


Re: [GRASS-user] Using v.distance between two identical layers...

2010-06-02 Thread Hamish
devaux wrote:
 do you know if it's possible to use v.distance command when
 for and to layers are
 the same layer ?
 I have a layer describing ecological areas of interest and
 I want to qualify spatial
 relation between all patches. I want to identify and
 calculate shortest relation between
 all my patches boundaries.
 
 How can I do as v.distance command doesn't have any option
 like : eliminate identical
 objects from the from and to layer ?
 
 I'd like to use this particular command because of the
 option that create lines between
 nearest objects.


use g.copy or v.extract as a first step to separate/duplicate the
features?


Hamish


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


Re: [GRASS-user] Select random cells from map

2010-06-02 Thread Hamish
  António wrote:
  I need to randomly select non-null cells from a given map.
  I tried r.random.cells but I cannot define a maximum number
  random cells and some selected random cells are not valid
  values in my map. is there any other function to randomly
  select cells from a raster?

Hamish:
  r.to.vect to get a vector area coverage of your raster (do
  not use the -s flag), then v.random.cover from wiki addons
  then v.to.rast, and finally r.mapcalc to transfer values from
  source map to random cell map? (probably we/I should make a
  new wrapper script to do that automatically)
 
  or you might skip the last steps and just use v.random.cover's
  raster= option to sample the raster map directly.
 
  probably the best idea is to fix r.random's cover= option to
  loop until it has enough points.
 
 
  also if you will be doing a lot of geostats it is probably
  worth your while to learn about the GRASS-R interface:
   http://grass.osgeo.org/wiki/R_statistics
 
  ... then the world is open to you.

António  wrote:
 Thanks for your answer. I tried r.random with my input data and the 
 number of cells and it worked but I do have one question: I didn't 
 understand the objective of cover_map. I used the map from I want to 
 extract random cells as input but what is the objective of cover_map?

r.random creates random points all over the current rectangular region
block.

Adding a cover map makes those random points only occur within the
irregular area covered by the cover map. (but r.random cover= is a bit
buggy right now; see new ticket in the bug trac'er)


Hamish




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


[GRASS-user] In Which OS should I develop: Linux or Windows

2010-06-02 Thread katrin eggert
Greetings

I'm GRASS user in Linux and in MS-Windows and Now I need to start
developing a few things. In which one should I start? In which GRASS is
more stable and robust?

Thanks

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