[GRASS-user] Re: r.grow.distance

2010-05-27 Thread fengel
Sorry for that. I hope this one's ok... http://osgeo-org.1803224.n2.nabble.com/file/n5107306/test_case.zip test_case.zip -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/r-grow-distance-tp5100098p5107306.html Sent from the Grass - Users mailing list archive at

Re: Re: [GRASS-user] Multi-scale Segmentation

2010-05-27 Thread Aleksander Yanovskiy
Hi Pedro There is also r.seg GRASS module for doing variational segmentation (http://www.ing.unitn.it/~grass/). Some information about it can be found in http://download.osgeo.org/osgeo/foss4g/2009/SPREP/2Thu/Parkside%20GO4/1500/Thu%20G04%201545%20Zatelli.pdf Regards, Aleksander Hi Nikos

[GRASS-user] Georrectify images without a coordinate system

2010-05-27 Thread Gilbert Ferrara
Greetings I have a couple of images, without any coordinates information, that I need to georrectify using GRASS and a georrectified image. My main concern is regarding r.in.gdal with -o option: If I do this, will it change my image resolution or extent? Thank you Gilbert

[GRASS-user] GRASS GUI assessment

2010-05-27 Thread Jenny Turner
Greetings I have sent an email a couple of weeks ago regarding a quick compare/assessment regarding GRASS versions. Now, I need to do the same but for GRASS GUI's. As far as I can see I can identify 3 GUIS: - Tcl/tk - WxPython - QGIS is there any comparison (Wikis, documents, ...) between those

Re: [GRASS-user] Warnings while importing Kompsat-2 images

2010-05-27 Thread Luis Lisboa
Regarding this issue, at GDAL mailing list they told me: This is not a critical error, only a warning. However, if you are in contact with the data producer you might encourage them to bring their TIFF production into agreement with the specification which requires that the last byte of these

[GRASS-user] Changing Mapset in the same LOCATION without exiting GRASS

2010-05-27 Thread Nikos Dumakis
Greetings I'm trying to find a fast way to change MAPSET, inside the same location, without having to exit GRASS and enter again. Can this be done? ___ grass-user mailing list grass-user@lists.osgeo.org

Re: [GRASS-user] Re: r.grow.distance

2010-05-27 Thread Markus Metz
fine in 7.0, broken in 6.x, this fixes it for me in 6.5: Index: main.c === --- main.c (revision 42353) +++ main.c (working copy) @@ -287,6 +287,12 @@ check(col, 1, -1); } + for (col = 0; col

[GRASS-user] Re: [beginner][r.to.vect] -v? polygonizing only white parts of an image

2010-05-27 Thread baobazz
Hello and thanks, Hamish, for the answer After thinning my image with r.thin, I described my thinned raster r.describe rasterThin * 1 The white parts of my image seemed to be * I launched r.reclass input=rasterThin output=rasterReclass 1=NULL *=1 black But I got everything equal to

[GRASS-user] Extracting grey features from an image

2010-05-27 Thread baobazz
Hello I've got a black and white image with some grey features on it I'd like to create a new image containing only these features so I can vectorize them i think that I have to use r.mask r.mapcalc or r.reclass r.describe raster returns me 0-121 123 125-127 132-255 I created a mask of my

Re: [GRASS-user] Re: [beginner][r.to.vect] -v? polygonizing only white parts of an image

2010-05-27 Thread Hamish
baobazz wrote: Just one question: r.describe rasterThin gave me * 1 Did * mean everything else than 1 or everything? Why did I get a blank/NULL map? * means NULL values. see the help page- http://grass.osgeo.org/grass64/manuals/html64_user/r.describe.html Hamish

Re: [GRASS-user] Changing Mapset in the same LOCATION without exiting GRASS

2010-05-27 Thread Hamish
Nikos Dumakis wrote: I'm trying to find a fast way to change MAPSET, inside the same location, without having to exit GRASS and enter again. Can this be done? http://grass.osgeo.org/grass64/manuals/html64_user/g.mapset.html Hamish ___

Re: [GRASS-user] Extracting grey features from an image

2010-05-27 Thread Micha Silver
On 05/27/2010 12:59 PM, baobazz wrote: Hello I've got a black and white image with some grey features on it I'd like to create a new image containing only these features so I can vectorize them i think that I have to use r.mask r.mapcalc or r.reclass r.describe raster returns me 0-121

Re: [GRASS-user] Multi-scale Segmentation

2010-05-27 Thread Nikos Alexandris
Aleksander Yanovskiy: There is also r.seg GRASS module for doing variational segmentation (http://www.ing.unitn.it/~grass/). Some information about it can be found in http://download.osgeo.org/osgeo/foss4g/2009/SPREP/2Thu/Parkside%20GO4/1500 /Thu%20G04%201545%20Zatelli.pdf Great piece

[GRASS-user] Using a polygon as a Mask

2010-05-27 Thread António Rocha
, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5149 (20100527) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ grass-user mailing list grass-user@lists.osgeo.org http

Re: [GRASS-user] Using a polygon as a Mask

2010-05-27 Thread Christian Schwartze
and only then I can use it in mapcalculator as a mask? Thanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5149 (20100527) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com

Re: [GRASS-user] Extracting grey features from an image

2010-05-27 Thread Christian Schwartze
Use r.mask as already done to identify your grey values only, then apply r.mapcalc new_raster=if(MASK,orig_raster,null()). Christian. -- From: baobazz baob...@hotmail.fr Sent: Thursday, May 27, 2010 11:59 AM To: grass-user@lists.osgeo.org

[GRASS-user] Re: Extracting grey features from an image

2010-05-27 Thread baobazz
Thanks a lot for the answers! I managed to identify the range of values containing the grey category values. So, my command became: r.mapcalc grey_raster=if(tragh.red=132tragh.red=200, 1, null()) It was a little difficult. I had to make it step by step. One way would be to use in addition

[GRASS-user] Re: [beginner][r.to.vect] -v? polygonizing only white parts of an image

2010-05-27 Thread baobazz
thanks for the answer I launched r.reclass input=rasterThin output=rasterReclass 1=NULL *=1 black so as to fill null values I don't understand why my resulting map was a NULL=blank raster.. -- View this message in context:

Re: [GRASS-user] Georrectify images without a coordinate system

2010-05-27 Thread Nikos Alexandris
Gilbert Ferrara: I have a couple of images, without any coordinates information, that I need to georrectify using GRASS and a georrectified image. My main concern is regarding r.in.gdal with -o option: If I do this, will it change my image resolution or extent? No. This flag is to be used if

[GRASS-user] Trying python scripts.

2010-05-27 Thread Marcello Gorini
Dear all, 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 other programming language), so I am trying to run the first examples found in

Re: [GRASS-user] Setting GRASS environment variable in a C program - step two

2010-05-27 Thread rabotin
If (inside a GRASS session), I access with v.db.connect -p map=GU0507 it looks ok : Vector map gu0...@simon9 is connected by: layer 1 table GU0507 in database /home/rabotin/grassdata/Roujan/simon9/dbf/ through driver dbf with key cat Do you know how to launch a grass command (e.g.

Re: [GRASS-user] Trying python scripts.

2010-05-27 Thread Nikos Alexandris
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 other programming language), so I am trying to run the first

Re: [GRASS-user] Using a polygon as a Mask

2010-05-27 Thread António Rocha
to convert it to RASTER and only then I can use it in mapcalculator as a mask? Thanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5149 (20100527) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com

[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. I am a total beginner in python and GRASS (and any

[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 to ANY mailing list, so please forgive me if I don't express myself

[GRASS-user] Re: Get current location projection in a Python Script

2010-05-27 Thread Kim Besson
Michael, regarding the use of regiondict = grass.region() , Does this retrieves Projection EPSG code? Kim 2010/5/25 Kim Besson kimbesson1...@gmail.com Greetings I'm building a GRASS Python Script to automatically import a set of data to my location. At GRASS book (chapter 3.3.3) it uses

Re: [GRASS-user] Using a polygon as a Mask

2010-05-27 Thread Nikos Alexandris
António Rocha wrote: What I am affraid is that when I convert Vector to RASTER I will obtain a mask that is not exacly the polygon in the boundaries of the polygon. For instance, if a pixel is cuttted by the polygons boundaries is it considered part odf the polygon or not? --%--- Antonio, an

Re: [GRASS-user] Multi-scale Segmentation

2010-05-27 Thread Markus Neteler
2010/5/27 Nikos Alexandris nikos.alexand...@felis.uni-freiburg.de: Aleksander Yanovskiy: There is also r.seg GRASS module for doing variational segmentation (http://www.ing.unitn.it/~grass/). Some information about it can be found in

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

2010-05-27 Thread Nikos Alexandris
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

Re: [GRASS-user] Multi-scale Segmentation

2010-05-27 Thread Nikos Alexandris
On Thursday 27 of May 2010 20:27:37 Markus Neteler wrote: 2010/5/27 Nikos Alexandris nikos.alexand...@felis.uni-freiburg.de: Aleksander Yanovskiy: There is also r.seg GRASS module for doing variational segmentation (http://www.ing.unitn.it/~grass/). Some information about it can be found

Re: [GRASS-user] Re: r.grow.distance

2010-05-27 Thread Glynn Clements
Markus Metz wrote: fine in 7.0, broken in 6.x, this fixes it for me in 6.5: Right; this was done for 7.0 in r34429, but it was mixed up with the addition of metric=geodesic, so it got overlooked. -- Glynn Clements gl...@gclements.plus.com ___

Re: [GRASS-user] Setting GRASS environment variable in a C program - step two

2010-05-27 Thread Glynn Clements
rabotin wrote: If (inside a GRASS session), I access with v.db.connect -p map=GU0507 it looks ok : Vector map gu0...@simon9 is connected by: layer 1 table GU0507 in database /home/rabotin/grassdata/Roujan/simon9/dbf/ through driver dbf with key cat Do you know how to launch a grass

Re: [GRASS-user] Re: Get current location projection in a Python Script

2010-05-27 Thread Glynn Clements
Kim Besson wrote: regarding the use of regiondict = grass.region() , Does this retrieves Projection EPSG code? No. You can get projection information as a Python dictionary with e.g.: s = grass.read_command(g.proj, flags='j') kv = grass.parse_key_val(s) In 7.0, use: fgrep

[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 the GRASS prompt, but when I tried to display the map, I found

Re: [GRASS-user] Re: Get current location projection in a Python Script

2010-05-27 Thread Hamish
Glynn wrote: However: I don't know whether the EPSG code is actually stored in the PROJ_INFO file, or whether it just stores the parameters. (it's not) If it's the latter, then there's no way to retrieve an EPSG code. (correct) ___

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

2010-05-27 Thread Nikos Alexandris
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 the outcome somewhere else?