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.
__
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 (
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
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
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---
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
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
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
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
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.
>
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
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
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
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
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
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'
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
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.
___
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
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
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
> >> "
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
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
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
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
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
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
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
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
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
:) 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
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
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
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
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
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.
>
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
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
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
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
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
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
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
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
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
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
## 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
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
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
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
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
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
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
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
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
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
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
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
>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,
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
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
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(
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
>"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).
>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.
>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
> 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
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
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
> 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
>>>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
>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
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,
>>
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
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
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
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
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
>>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
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
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:
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
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
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
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
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
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
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
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
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
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
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
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
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
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.
>
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
96 matches
Mail list logo