Hamish wrote:
anyway, I use this python one which is quite nice.
http://youarealegend.blogspot.co.nz/2008_09_01_archive.html
http://downloads.sourceforge.net/project/windrose/
Wow, very nice indeed! Thanks for the links!
Cheers,
Marcello.
___
Thanks for the suggestions.
Markus:
Please post
g.proj -p
-PROJ_INFO
-
name : Mercator
proj : merc
datum : wgs84
ellps : wgs84
lon_0 : 0
k : 1
x_0: 0
y_0: 0
no_defs: defined
Makus:
What are the extents of the vector in latlon before projecting into
the mercator location? Note that regions spanning 180 degrees
longitude or more can not be reprojected to mercator locations due to
mathematical constraints.20037508.34278924
That's the result of v.info of the of
Hamish said:
It is impossible to reproject from a simple XY location to a
projected location, or vice versa. Simple XY is just like graph
paper, with no Earth-based geo-* part to it. It is most commonly
used for imagery where x,y are measured in pixels, and manual geo-
referencing must be
Michael wrote:
Which way does it get flipped?
Most NetCDF files will import via r.in.gdal correctly, but they
are displayed as (from L-R) as 0-180E/180W-0. This is the way
that most climate model files are created and then stored.
It's not a climate model, but that's exactly the way
Dear all,
After importing a netCDF file and correcting its bounds using
r.region so that it spans from -180/180 instead of 0/360, a tried to
reproject it in the way I always do, following the GRASS book.
I use v.in.region to get the region as a vector, v.proj in a mercator
location and then
Dear all,
I am pretty sure this must be easy, but I just can't find the correct
answer.
My grid is a netcdf geographic grid that spans from 0E to 360E and from 90N
to -90S.
I use r.in.gdal to import it, but it gets flipped since GRASS region spans
from -180W to 180E.
How to do it right?
Thanks for he fast response
Daniel:
you should be able to define a location with that CRS using either the
correct EPSG code or a WKT string.
I don't have the EPSG code. Only know that it isin lat long and that long
is 0-360.
If need be, you can always create a WKT string as well,
DavidRA wrote:
Hi, here's the problem: I need to implement an operation in GRASS wich takes
an UNDEFINED number of raster layers and, for each position, takes the cell
of each raster, sort those values, makes some calculations and store the
result in the same position of the output raster.
G. Allegri wrote:
I need to assign values to a the cells on the border of a raster. The
inside and the outside are distinguished by having or not having null
values assigned.
I also need to keep the other cell values (internals) untouched.
Hey,
Check out r.grow. I am pretty sure you can
problem is not finding the border. I already have it (with
its own category) but assigning it the value from the surrounding (not
null) value, e.g. max/min/etc.
The problem with kernel filters/moving windows is that they do not filter
out null values...
giovanni
2012/6/19 Marcello Gorini
Büro Seling said:
My problem begins, when I now try to add up all these single maps into one
larger map.
I searched for a solution and tried several things, but nothing worked.
Any tipps or hinds?
Hey,
As always, there is probably a much more elegant solution to this problem,
but this
I knew it! Just testing you all :)
On Mon, Jun 18, 2012 at 12:54 PM, Paulo van Breugel
p.vanbreu...@gmail.comwrote:
First set the region to fit all maps, like suggested by Marcello. Then,
use r.patch to combine the maps.
Cheers,
Paulo
On 06/18/2012 01:36 PM, Marcello Gorini
Tanya wrote:
I am trying to run the grow.distance command to create a raster that
displays the euclidean distance from an 'orchard' outward to 1108m. I am
not trying to measure the distance between 2 non-null features, but simply
wish to have the euclidean distance extend out from the
On Wed, Jan 18, 2012 at 4:07 PM, Rich Shepard rshep...@appl-ecosys.comwrote:
I have a map with a couple of transect lines that intersect points on a
separate map. I want to measure the distance from one end of each transect
line to the point which it overlays. While I thought there was a v.
On Mon, Jan 16, 2012 at 10:37 PM, Rich Shepard rshep...@appl-ecosys.comwrote:
Somehow, my streets map also contains lines for streams. I would like to
subtract the stream lines from the street lines but do not see a module to
do this. Map algebra works well with raster maps (r.mapcalc), but
Marcello:
Dear all,
I need to extract the center of raster areas as another raster.
Since my script is totaly raster-based (and I am not used to using
vectors),
I tried to find some kind of r.thin that would generate the centroids,
instead of trying to force the area into lines.
Markus:
Dear all,
I need to extract the center of raster areas as another raster.
Since my script is totaly raster-based (and I am not used to using vectors),
I tried to find some kind of r.thin that would generate the centroids,
instead of trying to force the area into lines.
I believe there is no
Markus:
Please use
... GRASS GIS (© 1999-2011 GRASS Development Team) version is 6.4.2svn and
the main modules that compose
But the version is 6.4.0svn (2010) actually...
and then when I talk about a specific module I say like
the minimum-watershed-size threshold was set
OK, thanks.
And by the way... congrats on the GRASS paper!
On Thu, Dec 15, 2011 at 6:22 PM, Markus Neteler nete...@osgeo.org wrote:
On Thu, Dec 15, 2011 at 9:10 PM, Marcello Gorini gor...@gmail.com wrote:
Markus:
Please use
... GRASS GIS (© 1999-2011 GRASS Development Team
Humm I think you forgot the thing :) , no?
On Wed, Dec 14, 2011 at 5:22 PM, Stephen Sefick sas0...@auburn.edu wrote:
I saw a paper on this list the other day, but have deleted the original
email. Is this the appropriate thing to cite when referring to GRASS GIS
in a paper?
many thanks,
Stephen Sefick:
I saw a paper on this list the other day, but have deleted the original
email. Is this the appropriate thing to cite when referring to GRASS
GIS in
a paper?
Markus Neteler:
You may cite as this:
(from http://grass.osgeo.org/download/index.php)
* GRASS
Moritz:
As Glynn suggests:
for i = 1 to 10
sum_new = sum + i
g.remove rast=sum
g.rename rast=sum_new,sum
OK, sorry, I got so scared I overlooked Glynn's suggestion. I did another
experiment to see the differences in processing time. Since I iterate 1260
times in my script, I
Marcello:
I have to find the maximum value for each cell in a series
of rasters and generate both a maximum value map and a
categorical map with the raster number which contains the
maximum value.
Hamish:
see also:
Marcello:
But there is no way to avoid r.mapcalc, is there?
Hamish:
I doubt you'd get very much faster or more efficient even by
writing your own C module. (it would probably only take you a
minute to create a r.zero module out of the doc/raster/r.example
code; just change return x; to
Glynn:
Note that using the same map as both input and output results in
undefined behaviour. Even if it happens to work, there's no guarantee
that it will continue to work in future versions.
If you need to implement an iterative algorithm, you should generate a
new map, then replace the
Marcello:
My problem is that I also need to find the *second largest* value and the
corresponding raster number which contains the second largest value.
I am doing that by iterating over all classes through a shell script,
but it
obviously takes much more time than using a simple
Dear all,
I have to find the maximum value for each cell in a series of rasters and
generate both a maximum value map and a categorical map with the raster
number which contains the maximum value.
R.series is, therefore, perfect to accomplish that using
method=maximum,max_raster.
My problem is
Dear all,
Inside a loop in a shell script, I keep incrementing the values of a raster
map using r.mapcalc with some logical expressions.
Therefore, I need to initialize the raster map before the loop itself.
Ideally, I would like to initialize it as a raster with the same dimensions
as my base
:) how difficult, isn't it?
Thanks!
But there is no way to avoid r.mapcalc, is there?
On Thu, Dec 8, 2011 at 8:04 PM, Moritz Lennert mlenn...@club.worldonline.be
wrote:
On 08/12/11 22:26, Marcello Gorini wrote:
Dear all,
Inside a loop in a shell script, I keep incrementing the values
Moritz:
Since you have the max_raster value, can't you just take that map out of
the list you submit to r.series and find the new max_raster and value ?
Unfortunately, I can't. There is not on single map which is the maximum
everywhere. The calculation is done in a cell-by-cell basis, so
Hey Antonio,
The command v.in.region is exactly what you want. I use that all the time.
Cheers,
Marcello.
2011/11/23 António Rocha antonio.ro...@deimos.com.pt
Greetings
I would like to know how can I generate a polygon for the exact
computational region (a square). Is this possible? IF
Moritz:
Hello,
Does anyone see an easy way to count the number of times a specific
category value occurs in a given neighborhood of a central pixel ?
As always, there is surely a better way in GRASS (probably one single
command), but you could binarize your map (1 for your desired category
Hi Helena,
I have come to the same problem and solved it with ugly and very specific
shell scripts. There is most certainly a better way within GRASS and no
doubt you can code it much better in Shell, but anyway, this may get you
going somehow.
Code snippet:
MAP=$1
g.region rast=$MAP
See? I knew there was a much better way :)
On Tue, Nov 15, 2011 at 2:29 PM, Michael Barton michael.bar...@asu.eduwrote:
Use r.report.
This will do exactly what you want and output to a file too.
Michael
C. Michael Barton
Director, Center for Social Dynamics
Marcello Gorini wrote:
Or I guess you could export it with a larger number of pixels by setting
some environment variables.
export GRASS_WIDTH=
export GRASS_HEIGHT=
Glynn:
Those will affect the size of images generated by d.* commands using
the PNG driver. They won't have any
Luisa said:
I'm using r.out.png to generate PNG files from a small patch (like 6x15
pixels) but I'm obtaining a really small PNG file. Sicne I want to
display it a little bit bigger in a website I need to create a bigger
(in size) PNG. What can I do to do this?
Ben said:
You can simply increase
My friend, by your intention of helping Luisa, you ended up helping me a
lot.
I guess I am so attached to the GRASS book that I was conformed to the fact
that I could only export the active display monitor by using the PNG
driver.
D.out.file will help me a lot. Thanks.
Marcello.
On Thu, Nov
As far as I know, v.to.points does not write to an external file, it creates
a point vector map within GRASS.
So you should do output=pointmap and then use another module such as
v.out.ascii to generate an actual text file with the points.
Also, when you do that, do not enclose the output
Raphael wrote:
i have imported the etopo dem (geotifff) into grass.
It appears to be in XY and i would like to convert
these to lat longs.
Markus:
But in any case the ETOPO2 will be better (see wiki for an
article how to import it).
Hamish:
yes for ease of import, but not so much for
## SOLVED ##
The problem was located in-between the computer and my chair :)
Thanks to Hamish for pointing it.
Cheers,
Marcello.
Marcello Gorini wrote:
Sorry for this desperate message, but I am having all kinds of troubles in
the worst moment ever!
My GRASS 6.4.0RC6 on ubuntu 9.10
Glynn Clements wrote:
It looks like you're trying to build a 6.x add-on against 7.0.
Probably because of:
svnurl=http://svn.osgeo.org/grass/grass-addons/
For 7, hits should be
svnurl=https://svn.osgeo.org/grass/grass-addons/grass7
That' s exactly what I am doing, thanks. I
On Mon, Sep 19, 2011 at 6:43 PM, Marcello Gorini gor...@gmail.com wrote:
...
I was advised not to go into grass 7.0 if I wanted a fast result
On Mon, Sep 19, 2011 at 7:24 PM, Markus Neteler nete...@osgeo.org wrote:
... are you sure about this?
Sorry, I must explain better.
I am
Sorry for this desperate message, but I am having all kinds of troubles in
the worst moment ever!
My GRASS 6.4.0RC6 on ubuntu 9.10 (old and in need of update, I know)
suddenly stopped being able to find all my own shell scripts in the scripts
folder and also all add-ons that I had installed. Some
Sorry for this another desperate message, but I am really having all kinds of
troubles in the worst moment ever!
In order to bypass the problem in my previous post, I installed both GRASS
7.0 from svn checkout https://svn.osgeo.org/grass/grass/trunk grass_trunk
and GRASS 6.4.1 from Synaptic in
, such
as Matlab, for instance?
Thanks in advance.
Marcello.
Marcello Gorini wrote:
Dear all,
I am interested in extracting the longest -1,000 meter contour from a
bathymetric DEM and to save it in a text file as points.
By doing:
r.contour in=my_dem out=contour levels=-1000
v.to.points
Dear all,
I am interested in extracting the longest -1,000 meter contour from a
bathymetric DEM and to save it in a text file as points.
By doing:
r.contour in=my_dem out=contour levels=-1000
v.to.points in=contour out=contour_points dmax=my_resolution
v.out.ascii in=contour_points
r.mapcalc new_raster=if(old_raster-0.01 old_raster0.01,0,old_raster)
Cheers,
Marcello.
On Wed, Jul 13, 2011 at 10:10 AM, Rebecca Bennett rabenn...@ymail.comwrote:
Hello Grass users,
I have a (hopefully) quick question - I would like to replace all values in
a raster that fall in the range
Dear all,
This is just a follow-up of my original message as I promissed.
It has been two weeks since I tried to contact the author, but received no
response.
I eventually used a workaround based on r.random.surface itself. I added
random fields to the original continuous parameters instead of
Dear all,
I am using r.fuzzy.system for quite a while now, but still have a persisting
bug/problem that I can't solve. Hopefully it is just me doing something
wrong.
A use a GRASS/Shell script to automatically parse the input .map file (that
contains all the fuzzy sets definitions) and update it
Hamish wrote:
also, for a categorical map you might try r.to.vect + the v.random.cover
module for GRASS 6 from addons:
http://grass.osgeo.org/wiki/Addons#v.random.cover
v.random.cover is a shell script for creating random points constrained
within an irregularly shaped vector area. (v.random
On Wed, Jun 22, 2011 at 2:24 AM, Marcello Gorini lt;gor...@gmail.comgt;
wrote:
Deal all,
I am interested in a module that seems to be discontinued since GRASS 4.
It
is called r.random.model and it is used for Monte Carlo realizations of
categorical maps.
Rainer:
Me to - sounds really
in advance.
Best regards,
Marcello Gorini.
--
View this message in context:
http://osgeo-org.1803224.n2.nabble.com/Monte-Carlo-realizations-of-categorical-maps-tp6502447p6502447.html
Sent from the Grass - Users mailing list archive at Nabble.com.
___
grass
Marcello wrote:
I solved it using:
d.where map_profile
instead of:
r.profile -ig input=map output=map_profile
The only problem is that I can't see the points while I am digitizing.
Hamish wrote:
use the middle mouse button with d.where to show the line.
Yes, thanks, it helped.
Janet Choate said:
i want to display a raster in 3D. i.e. display landscape maps 'draped on
top' of the dem so i can show the terrain in 3D rather than as a flat
image.
do i need to convert a raster to a 3D raster via r.to.rast3 or
r.to.rast3elev, set the region for 3D?
how to display?
Not at
Dear all,
My problem (a very common one) is that I want to digitize a profile and then
sample other maps using the extracted coordinates. In order to do so, I
tried the following;
r.profile -ig input=map output=map_profile
With that, I managed to digitize the profile and to get the following
, instead of sampling the raster.
Still don't know why the method with r.profile doesn't work right though. If
anybody does, I would be glad to hear.
Thanks,
Marcello.
Marcello Gorini wrote:
Dear all,
My problem (a very common one) is that I want to digitize a profile and
then sample other
Luisa Peña wrote:
Greetings
I have a raster map with a size like 25000 x 3 with valid values on
all
its extent. With another raster, I defined a sub-extent of this raster,
much
smaller where only a small area is valid and the rest it NULL() like this:
r.mapcalc output =
Nepomuk Reinhard wrote:
Dear all,
I'm a newbie in grass GIS so I'm sorry, if this question is very simple.
Has anyone any idea how to create a raster map including only the
information of the latitude and longitude of each cell? I want to write a
script to compute the solar elevation
Kim:
I need to calculate the average of a set of raster images. How can I
calculate it without using mapcalc?
THanks
Hello Kim,
You can use the module r.series. Like this, if you have few images:
r.series input=image1,image2,image3... method=average output=averaged_image
Or if you have a
Jenny:
in my base map (called X) I have a few null values and, in Y map I want to
eliminate pixels that match with null values in X. How can I do this?
besides using mapcalc
Hello Jenny,
I don't know if it is a coincidence, but your problem can be solved in the
same way as Kim's problem.
You
Aren Cambre wrote:
I am not getting any responses. :-) Am I doing something really dumb? Is
this the right forum for this kind of question?
Thanks,
Aren
On Sun, Feb 13, 2011 at 8:02 AM, Aren Cambre a...@arencambre.com wrote:
Maybe I could clarify. Here's the command:
*v.kernel
Aren said:
Where is resolution adjusted?
http://grass.fbk.eu/gdp/html_grass64/v.kernel.html doesn't make that clear.
GRASS works with the region concept. You can adjust any region settings with
the command g.region, like:
g.region res=desired_resolution -p
The -p is for you to see the
Hello,
I cannot really answer your question, but I have been working with a
combination of your fuzzy modules for GRASS and r.param.scale in order to
accomplish a fuzzified version of Wood's morphometric feature extraction
method.
In doing so, I noticed that the results of r.param.scale
Timmie says:
An assiciated question:
How do I get a diagram of the frequency distribution of
the floating point values of a raster?
I tried histogram but I cannot get the label on the x and y axis.
Something like a plot of the r.report output.
Hello Timmie,
Not really an answer, but maybe
Pedro Roma:
I want to develop a script that crops a raster in order to only have the
pixels that are located inside a vector polygon.
There is probably a more elegant solution, but this will do:
g.region vect=your_vector
v.to.rast input=your_vector output=your_vector_rasterized use=val
r.mask
Pedro Roma:
I want to develop a script that crops a raster in order to only have the
pixels that are located inside a vector polygon.
Marcello Gorini:
There is probably a more elegant solution, but this will do:
g.region vect=your_vector
v.to.rast input=your_vector output
On Mon, Aug 16, 2010 at 10:20 PM, Marcello Gorini gor...@gmail.com wrote:
Dear all,
By comparing morphometric feature maps obtained from both r.param.scale
in
GRASS and the function sufparam using LandScript in Landserf (Wood,
1996-2009), I noticed great discrepancies.
Actually
Hello Daniel,
I also wanted to use r.area, so I installed it but had an error, as posted
here
http://osgeo-org.1803224.n2.nabble.com/New-modules-in-svn-add-on-r-area-and-r-convergence-td5161926.html#a5199839
Then afterwards Milton Ribeiro sent me another way to get the same results
in a rather
Kwas wrote:
Hello,
I have a raster map (dem) and I want to select a region of the raster
corresponding to a polygon that is in a vector map.nbsp; What is the best
way to do this?nbsp; It looks like r.in.poly might be along the lines of
what I need but in that case I think I would need
António Rocha:
Hi there
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?
Hello
On Mon, May 31, 2010 at 12:11 AM, Marcello Gorini gor...@gmail.com wrote:
OK, right after a wrote that, I thought I was being too bold.
Should I go for:
svn checkout https://svn.osgeo.org/grass/grass/branches/releasebranch_6_4
grass64_release
Markus:
Yes.
Subsequently, it is sufficient
a little
little bit in this list. It is good to at least try to contribute.
Best regards,
Marcello.
2010/5/31 Marcello Gorini gor...@gmail.com
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Nikos Alexandris wrote:
Marcello Gorini:
...
I am a total beginner in python and GRASS (and any other programming
language), so I am trying to run the first examples found in
http://grass.osgeo.org/wiki/GRASS_and_Python.
Marcello,
just fyi: I've added another simple example script
Martin Landa wrote:
Hi,
2010/5/30 Marcello Gorini gor...@gmail.com:
I am a total beginner in python and GRASS (and any other programming
language), so I am trying to run the first examples found in
http://grass.osgeo.org/wiki/GRASS_and_Python.
GRASS 6.4.0RC5 is quite old. Try RC6
OK, right after a wrote that, I thought I was being too bold.
Should I go for:
svn checkout https://svn.osgeo.org/grass/grass/branches/releasebranch_6_4
grass64_release
or this one?
svn checkout
https://svn.osgeo.org/grass/grass/tags/release_20090609_grass_6_4_0RC5
grass640_rc5
Thanks again
On Monday 31 of May 2010 00:11:00 Marcello Gorini wrote:
Should I go for:
svn checkout https://svn.osgeo.org/grass/grass/branches/releasebranch_6_4
grass64_release
or this one?
svn checkout
https://svn.osgeo.org/grass/grass/tags/release_20090609_grass_6_4_0RC5
grass640_rc5
Nikos
hamish-2 wrote:
Marcello wrote:
checking for lex... no
configure: error: *** Unable to locate lex.
...
P.S. I use Ubuntu 9.10 64 bits.
you need to install all the Build-depends packages listed here:
http://svn.debian.org/viewsvn/pkg-grass/packages/grass/trunk/debian/control
(
Marcello Gorini wrote:
hamish-2 wrote:
Marcello wrote:
checking for lex... no
configure: error: *** Unable to locate lex.
...
P.S. I use Ubuntu 9.10 64 bits.
you need to install all the Build-depends packages listed here:
http://svn.debian.org/viewsvn/pkg-grass/packages/grass
Nikos Alexandris wrote:
Marcello Gorini wrote:
Thanks for the try. I tested it again with other rasters, but I got the
same error. Could it possibly be something related to the 64bits
architecture?
I don't think so. I work with Kubuntu Lucid 64 bit :-)
Or could GRASS be sending
Nikos Alexandris wrote:
Marcello Gorini:
...let's go step by step:
ls -l myscript.py
-rwxr-xr-x 1 marcello marcello 1227 2010-05-28 08:13 myscript.py
looks fine
g.gisenv
GISDBASE=/home/marcello/grassdata
LOCATION_NAME=Ocean_floor
MAPSET=level4
MONITOR=x0
GRASS_GUI=tcltk
, Marcello Gorini gor...@gmail.com wrote:
Nikos Alexandris wrote:
Marcello Gorini wrote:
Thanks for the try. I tested it again with other rasters, but I got the
same error. Could it possibly be something related to the 64bits
architecture?
I don't think so. I work with Kubuntu Lucid 64 bit
Nikos:
Just stick on Marcello. Get yourself the GRASS-book, look at the
GRASS-wiki,
invest some time generally and specifically, and you will be addicted ;-)
It
is worthwhile.
I will, don't worry. And I already bought the GRASS book. That is what
really got me going (after Carlos got me
output=outmap raster1=map1 raster2=map2
I am using Ubuntu 9.04 64bits and GRASS 6.4.ORC5.
Can anyone help me?
Best regards,
Marcello Gorini.
--
View this message in context:
http://osgeo-org.1803224.n2.nabble.com/Trying-python-scripts-tp5109105p5109105.html
Sent from the Grass - Users mailing
Hello Nikos,
Marcello Gorini.
Nikos Alexandris wrote:
Hi Marcello,
Marcello Gorini:
This is kind of my first e-mail ever to ANY mailing list, so please
forgive
me if I don't express myself in the best way. I will improve eventually.
I am a total beginner in python and GRASS (and any
Oops, I guess I quoted the text in the wrong way. Sorry for that.
Marcello Gorini wrote:
Hello Nikos,
Marcello Gorini.
Nikos Alexandris wrote:
Hi Marcello,
Marcello Gorini:
This is kind of my first e-mail ever to ANY mailing list, so please
forgive
me if I don't express myself
Marcello Gorini wrote:
[...]
... I am trying to run the first examples found in
http://grass.osgeo.org/wiki/GRASS_and_Python.
--%---
a script was here - look first post !
---%--
It seems to run OK and “100%” is displayed in the GRASS prompt, but
when I tried to display the map, I found
87 matches
Mail list logo