Re: [GRASS-user] How I can develope code for Polar Diagram(Rose Diagram) in windows

2012-09-05 Thread Marcello Gorini
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. __

Re: [GRASS-user] How to implement this kind of operation?

2012-07-02 Thread Marcello Gorini
DavidRA wrote: I guess I could use the code of r.series, modifying it a bit to obtain a new > function. The problem is that... man, that code is hard to understand. Sorry if I gave you the impression that I am a skilled programmer in C. Actually, I forgot the majority of GRASS modules are in C (

Re: [GRASS-user] Invalid region and coordinates when trying to reproject.

2012-07-02 Thread Marcello Gorini
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

Re: [GRASS-user] Invalid region and coordinates when trying to reproject.

2012-07-02 Thread Marcello Gorini
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

Re: [GRASS-user] Invalid region and coordinates when trying to reproject.

2012-07-01 Thread Marcello Gorini
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 -PROJ_UNITS---

[GRASS-user] Invalid region and coordinates when trying to reproject.

2012-07-01 Thread Marcello Gorini
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 g.reg

Re: [GRASS-user] Correct import of a global grid (0E to 360E) to GRASS (-180W to 180E).

2012-06-30 Thread Marcello Gorini
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 exa

Re: [GRASS-user] Correct import of a global grid (0E to 360E) to GRASS (-180W to 180E).

2012-06-30 Thread Marcello Gorini
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, al

[GRASS-user] Correct import of a global grid (0E to 360E) to GRASS (-180W to 180E).

2012-06-30 Thread Marcello Gorini
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

Re: [GRASS-user] How to implement this kind of operation?

2012-06-29 Thread Marcello Gorini
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. >

Re: [GRASS-user] calculate mode value on moving window (with mapalgebra) with null value cells

2012-06-19 Thread Marcello Gorini
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 > &g

Re: [GRASS-user] calculate mode value on moving window (with mapalgebra) with null value cells

2012-06-19 Thread Marcello Gorini
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

Re: [GRASS-user] Add up several maps, which hold different regions, with mapcalc

2012-06-18 Thread Marcello Gorini
I knew it! Just testing you all :) On Mon, Jun 18, 2012 at 12:54 PM, Paulo van Breugel wrote: > 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

Re: [GRASS-user] Add up several maps, which hold different regions, with mapcalc

2012-06-18 Thread 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 t

Re: [GRASS-user] Euclidean distance (grow.distance) on raster to certain extent?

2012-05-07 Thread 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 o

Re: [GRASS-user] Measuring Distance Along a Line

2012-01-18 Thread Marcello Gorini
On Wed, Jan 18, 2012 at 4:07 PM, Rich Shepard wrote: > 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. module > to do this, I'

Re: [GRASS-user] Subtracting One Vector Map From Another

2012-01-16 Thread Marcello Gorini
On Mon, Jan 16, 2012 at 10:37 PM, Rich Shepard wrote: > 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 what > module will allo

Re: [GRASS-user] Re: Extracting centroids as a raster (and are these centroids really at the center?)

2012-01-09 Thread Marcello Gorini
Markus: Honestly, I tried only with a one area map > and updated the script to the > posted stage after midnight :p > All I can do is to thank you very much for your effort in updating the script, i.e., in helping me. It really did... a lot! Cheers, Marcello. ___

[GRASS-user] Re: Extracting centroids as a raster (and are these centroids really at the center?)

2012-01-09 Thread Marcello Gorini
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

[GRASS-user] Extracting centroids as a raster (and are these centroids really at the center?)

2012-01-03 Thread Marcello Gorini
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 s

Re: [GRASS-user] GRASS citation?

2011-12-15 Thread Marcello Gorini
OK, thanks. And by the way... congrats on the GRASS paper! On Thu, Dec 15, 2011 at 6:22 PM, Markus Neteler wrote: > On Thu, Dec 15, 2011 at 9:10 PM, Marcello Gorini wrote: > > > > Markus: > > > >> > >> > >> Please use > >> "

Re: [GRASS-user] GRASS citation?

2011-12-15 Thread Marcello Gorini
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 thres

Re: [GRASS-user] GRASS citation?

2011-12-14 Thread Marcello Gorini
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) > > * GRA

Re: [GRASS-user] GRASS citation?

2011-12-14 Thread Marcello Gorini
Humm I think you forgot "the thing" :) , no? On Wed, Dec 14, 2011 at 5:22 PM, Stephen Sefick 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, > > -- > Stephe

Re: [GRASS-user] How to initialize raster maps.

2011-12-10 Thread Marcello Gorini
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

Re: [GRASS-user] How to find the second largest value in a series of rasters.

2011-12-09 Thread Marcello Gorini
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 s

Re: [GRASS-user] How to initialize raster maps.

2011-12-09 Thread Marcello Gorini
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

Re: [GRASS-user] How to initialize raster maps.

2011-12-09 Thread Marcello Gorini
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 "retur

Re: [GRASS-user] How to find the second largest value in a series of rasters.

2011-12-09 Thread Marcello Gorini
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: > > http://grass.osgeo.org/wiki/Time_series#Common_legends_fo

Re: [GRASS-user] How to find the second largest value in a series of rasters.

2011-12-08 Thread Marcello Gorini
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

Re: [GRASS-user] How to initialize raster maps.

2011-12-08 Thread Marcello Gorini
:) 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 wrote: > On 08/12/11 22:26, Marcello Gorini wrote: > >> Dear all, >> >> Inside a loop in a shell script, I keep incrementing the

[GRASS-user] How to initialize raster maps.

2011-12-08 Thread Marcello Gorini
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 m

[GRASS-user] How to find the second largest value in a series of rasters.

2011-12-08 Thread Marcello Gorini
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 t

Re: [GRASS-user] r.fillnulls for large area

2011-12-01 Thread Marcello Gorini
Nick wrote: > > Maybe I need to split my region into subregions > > and run r.fillnulls on each subregion then patch > > them together. > Hamish: > yes, that should work well. I always wanted to create a script to do that. Since (I think) r.fill.nulls identifies the null locations and create

Re: [GRASS-user] how to count number of occurrences of a specific category value in the neighborhood of a given pixel

2011-11-23 Thread Marcello Gorini
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 catego

Re: [GRASS-user] Generate a box that matches the exact processing region

2011-11-23 Thread Marcello Gorini
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 > Greetings > > I would like to know how can I generate a polygon for the exact > computational region (a square). Is this possible? IF so, how can I do that. >

Re: [GRASS-user] Alternative to r.stats to generate % of surface

2011-11-15 Thread Marcello Gorini
See? I knew there was a much better way :) On Tue, Nov 15, 2011 at 2:29 PM, Michael Barton wrote: > 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 & Complexity > Profess

Re: [GRASS-user] Alternative to r.stats to generate % of surface

2011-11-15 Thread Marcello Gorini
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 max_cats=9

Re: [GRASS-user] Re: GRASS GIS development visualization from 1999 to 2011

2011-11-15 Thread Marcello Gorini
FANTASTIC, AMAZING, INSANE!! Just a few words to describe this video! Congratulations!! On Mon, Nov 14, 2011 at 3:30 PM, Markus Neteler wrote: > On Fri, Nov 11, 2011 at 4:19 PM, Markus Neteler wrote: > > I have created a video of the GRASS GIS development from 1999 to 2011: > > http://grass

Re: [GRASS-user] How to generate a "bigger PNG file"

2011-11-04 Thread Marcello Gorini
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

Re: [GRASS-user] Re: How to generate a "bigger PNG file"

2011-11-03 Thread Marcello Gorini
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 3

Re: [GRASS-user] How to generate a "bigger PNG file"

2011-11-03 Thread Marcello Gorini
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 simpl

Re: [GRASS-user] v.to.points

2011-10-18 Thread Marcello Gorini
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 portion

[GRASS-user] Re: xy to lat long

2011-09-27 Thread Marcello Gorini
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

Re: [GRASS-user] Re: Desperate help needed 2: can' t install add-ons in other versions of GRASS.

2011-09-19 Thread Marcello Gorini
On Mon, Sep 19, 2011 at 6:43 PM, Marcello Gorini 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 wrote: > ... are you sure about this? > > Sorry, I must explain better. I am in a hurry to del

[GRASS-user] Re: Desperate help needed 2: can' t install add-ons in other versions of GRASS.

2011-09-19 Thread Marcello Gorini
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

[GRASS-user] Re: Desperate help needed: GRASS can' t find any add-ons anymore.

2011-09-19 Thread Marcello Gorini
## 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

[GRASS-user] Desperate help needed 2: can' t install add-ons in other versions of GRASS.

2011-09-18 Thread Marcello Gorini
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 ano

[GRASS-user] Desperate help needed: GRASS can' t find any add-ons anymore.

2011-09-18 Thread Marcello Gorini
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

[GRASS-user] Re: How to extract one single contour.

2011-08-01 Thread Marcello Gorini
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_d

[GRASS-user] How to extract one single contour.

2011-07-31 Thread Marcello Gorini
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 out=contour_p

Re: [GRASS-user] Mapcalc to replace a range of values

2011-07-13 Thread Marcello Gorini
r.mapcalc "new_raster=if(old_raster>-0.01 && old_raster<0.01,0,old_raster)" Cheers, Marcello. On Wed, Jul 13, 2011 at 10:10 AM, Rebecca Bennett wrote: > Hello Grass users, > > I have a (hopefully) quick question - I would like to replace all values in > a raster that fall in the range -0.01 to

[GRASS-user] R.fuzzy.system bug/limitation or just me ?

2011-07-05 Thread Marcello Gorini
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

[GRASS-user] Re: Monte Carlo realizations of categorical maps.

2011-07-05 Thread Marcello Gorini
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

[GRASS-user] Re: Monte Carlo realizations of categorical maps.

2011-06-23 Thread Marcello Gorini
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.r

[GRASS-user] Re: Monte Carlo realizations of categorical maps.

2011-06-22 Thread Marcello Gorini
On Wed, Jun 22, 2011 at 2:24 AM, Marcello Gorini <gor...@gmail.com> 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 >> cate

[GRASS-user] Monte Carlo realizations of categorical maps.

2011-06-21 Thread Marcello Gorini
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

[GRASS-user] Re: create 3D raster and display

2011-04-06 Thread Marcello Gorini
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? N

[GRASS-user] Re: Duplication of lines in output of r.profile.

2011-04-06 Thread Marcello Gorini
>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,

[GRASS-user] Re: Duplication of lines in output of r.profile.

2011-04-05 Thread Marcello Gorini
only gets the coordinates, 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) i

[GRASS-user] Duplication of lines in output of r.profile.

2011-04-05 Thread Marcello Gorini
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 f

[GRASS-user] Re: Region definition

2011-03-26 Thread Marcello Gorini
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 = if(

[GRASS-user] Re: Creating raster of geographic coordinates

2011-03-03 Thread Marcello Gorini
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 elevati

[GRASS-user] Re: third raster -> as a function of two input raster (r.series?)

2011-02-26 Thread Marcello Gorini
>"Gabriele N." says: > >I have two raster. > >I would like in output a third raster that contains only the areas that have high values of the first >raster >and with high values of the second raster (and with the areas that have low values of the first >raster >and low values of the second raster).

[GRASS-user] Re: Eliminating values based on a null

2011-02-16 Thread Marcello Gorini
>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.

[GRASS-user] Re: Calculate averagfe of raster maps without mapcalc

2011-02-16 Thread Marcello Gorini
>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

[GRASS-user] Re: v.kernel -- how to get useful results

2011-02-15 Thread Marcello Gorini
> 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 resul

[GRASS-user] Re: v.kernel -- how to get useful results

2011-02-15 Thread Marcello Gorini
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 wrote: > >> Maybe I could clarify. Here's the command: >> *v.kernel --verbo

[GRASS-user] Re: r.param.scale curvatures: are really OK

2010-11-02 Thread Marcello Gorini
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 param=f

[GRASS-user] Re: reclassify raster based on distribution

2010-10-27 Thread Marcello Gorini
> 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

[GRASS-user] Re: A Script that only selects raster inside a polygon

2010-09-24 Thread Marcello Gorini
>>>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 ve

[GRASS-user] Re: A Script that only selects raster inside a polygon

2010-09-24 Thread Marcello Gorini
>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.m

[GRASS-user] Re: Feature extraction in GRASS and Landserf.

2010-08-25 Thread Marcello Gorini
On Mon, Aug 16, 2010 at 10:20 PM, Marcello Gorini 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, >>

[GRASS-user] Feature extraction in GRASS and Landserf.

2010-08-16 Thread Marcello Gorini
erf Version 2.3. I thank you all very much in advance. Best regards, Marcello Gorini. -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Feature-extraction-in-GRASS-and-Landserf-tp5429470p5429470.html Sent from the Grass - Users mailing l

[GRASS-user] Re: Installing r.area

2010-08-01 Thread Marcello Gorini
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 f

[GRASS-user] Re: r.terracost and gem

2010-06-25 Thread Marcello Gorini
Benjamin Ducke wrote: > > See here: > > http://grass.osgeo.org/wiki/GRASS_AddOns#Building_and_installing_Addons > > GEM is installed autaomatically as part of GRASS 6.4. > But it will not work for installing r.terracost. > > Ben > > - Oryginalna wiadomość - >> hello all, I have a pro

[GRASS-user] Re: New modules in svn add-on: r.area and r.convergence

2010-06-19 Thread Marcello Gorini
Jarek Jasiewicz wrote: > > Hi > > I added new two modules to svn add-on. > > r.area > > Simple module available in almost every raster processing GIS > application. Calculates area size of clumped areas and ascribe these > values to these areas. Also allow remove small areas with size lower

[GRASS-user] Re: how to get raster region corresponding to vector polygon region?

2010-06-17 Thread Marcello Gorini
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.  What is the best > way to do this?  It looks like r.in.poly might be along the lines of > what I need but in that case I think I would need to

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

2010-06-01 Thread Marcello Gorini
>>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 rast

Re: [GRASS-user] Re: Trying python scripts.

2010-05-31 Thread Marcello Gorini
lse) could try them for you. It would be good for my learning as well. >Thanks! >Daniel No, I thank you, since you are the first person I might be helping a little little bit in this list. It is good to at least try to contribute. Best regards, Marcello. 2010/5/31 Marcello Gorini ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Re: Trying python scripts.

2010-05-31 Thread Marcello Gorini
On Mon, May 31, 2010 at 12:11 AM, Marcello Gorini 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:

[GRASS-user] Re: Trying python scripts.

2010-05-30 Thread Marcello Gorini
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 al

[GRASS-user] Re: Trying python scripts.

2010-05-30 Thread Marcello Gorini
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

[GRASS-user] Re: Trying python scripts.

2010-05-30 Thread Marcello Gorini
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://sv

[GRASS-user] Re: Trying python scripts.

2010-05-30 Thread Marcello Gorini
stephen sefick wrote: > > This is a script that I wrote this morning to grab the latest 6.4 svn > sources and the addons, and then compile the addons that I want and > install grass. the sources are checked out to > ~/sources/grass/grass64_release and ../grass-addons. I don't know if > you are

[GRASS-user] Re: Trying python scripts.

2010-05-30 Thread Marcello Gorini
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

[GRASS-user] Re: Trying python scripts.

2010-05-30 Thread Marcello Gorini
Martin Landa wrote: > > Hi, > > 2010/5/30 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_P

[GRASS-user] Re: Trying python scripts.

2010-05-30 Thread Marcello Gorini
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, > &g

[GRASS-user] Re: Trying python scripts.

2010-05-28 Thread Marcello Gorini
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 intere

[GRASS-user] Re: Trying python scripts.

2010-05-28 Thread Marcello Gorini
that? > > best > > Carlos > > > > On Fri, May 28, 2010 at 08:23, Marcello Gorini wrote: >> >> >> Nikos Alexandris wrote: >>> >>> Marcello Gorini wrote: >>>> Thanks for the try. I tested it again with other rasters, but I

[GRASS-user] Re: Trying python scripts.

2010-05-28 Thread Marcello Gorini
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

[GRASS-user] Re: Trying python scripts.

2010-05-28 Thread Marcello Gorini
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

[GRASS-user] Re: Trying python scripts.

2010-05-27 Thread Marcello Gorini
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 t

[GRASS-user] Re: Trying python scripts.

2010-05-27 Thread Marcello Gorini
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 t

[GRASS-user] Re: Trying python scripts.

2010-05-27 Thread Marcello Gorini
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. >

[GRASS-user] Trying python scripts.

2010-05-27 Thread Marcello Gorini
the script I simply create it (say “myscript.py”), give permissions and run it in a GRASS session with the necessary options: > ./myscript.py 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 t