[GRASS-user] using g.mlist to set g.region
Dear all, I have a set of raster maps and I need to run a r.series with them. But each raster have different extent. How can I use g.mlist to get the extent that include all rasters to redefine my working regions? I am running grass under MSYS/Win7 best -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Calculate average of 3 highest values for a set of raster maps
Dear all, Thanks for all replies. I will check it out thoughout the weekend. Just clarifying something, I need to do the estimations focusing on each pixels, and not for the full map (I was not that clear before). Also as my raster maps are very large (~50,000x60,000 pixels) maybe R (at least on w7) maybe will not handle the seven rasters. Best wishes miltinho 2013/3/15, Paulo van Breugel p.vanbreu...@gmail.com: Sorry, some typo's (but the given code should be ok): set all values 4/7 quantile to null(), not 0, and you want to calculate the average, not the sum On Fri 15 Mar 2013 10:05:51 AM CET, Paulo van Breugel wrote: Hi Miltinho I assume you want the mean of the three highest values for each raster cell? If so, you can do this using another step with r.series 1) If your 7 maps are a1 .. a7: mapcalc a1 = 1; r.mapcalc a2 = 2; r.mapcalc a3 = 3; r.mapcalc a4 = 4; r.mapcalc a5 = 5; r.mapcalc a6 = 6; r.mapcalc a7 = 7 2) calculate the 4/7 quantile. Actually, it should be a little bit smaller to make sure you include your 3rd largest layer. In this example below, 0.57 as approximation is fine r.series --overwrite input=a1,a2,a3,a4,a5,a6,a7 quantile=0.57 output=res1 method=quantile 2) Next, for all 7 layers, set all values 4/7 quantile to 0 r.mapcalc b1 = if(a1res1,a1,null()) r.mapcalc b2 = if(a2res1,a2,null()) r.mapcalc b3 = if(a3res1,a3,null()) r.mapcalc b4 = if(a4res1,a4,null()) r.mapcalc b5 = if(a5res1,a5,null()) r.mapcalc b6 = if(a6res1,a6,null()) r.mapcalc b7 = if(a7res1,a7,null()) 3) Now, calculate the sum of the 7 new layers r.series input=b1,b2,b3,b4,b5,b6,b7 method=average output=res2 --overwrite Your output res2 should have a value 6 for all raster cells. Cheers, Paulo On Fri 15 Mar 2013 03:57:35 AM CET, Milton Cezar Ribeiro wrote: I have a set of 7 raster maps. I am able to calculate several stats with these mapas using r.series. Now I need to calculate the mean of the tree highest values of my 7 maps. Any idea of how can I do that? cheers -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Calculate average of 3 highest values for a set of raster maps
I have a set of 7 raster maps. I am able to calculate several stats with these mapas using r.series. Now I need to calculate the mean of the tree highest values of my 7 maps. Any idea of how can I do that? cheers -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.proj - how to use different dbase?
Dear all, I have two grass spatial database, each on different directory and with different projection/datum. How can I import a raster map while reprojectingit? Just suppose I am connected on C:\Mydata\GrassSDB_target newLocation PERMANENT and the source data are on C:\users\GrassSDB_source newLocation PERMANENT mymap_UTM_sad69 I tryed something like r.proj input=mymap_UTM_sad69 location=newLocation mapset=PERMANENT dbase=C:\users\GrassSDB_source But I get ERRO:Mapset PERMANENT in input location newLocation - not found Any idea of how to solve this? cheers milton -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.proj - how to use different dbase?
Hi Markus Thanks. I suggest, if it is a know but, put this on manual page of GRASS while not fixed :-) Cheers miltinho 2013/2/16, Markus Metz markus.metz.gisw...@gmail.com: On Sat, Feb 16, 2013 at 4:50 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Dear all, I have two grass spatial database, each on different directory and with different projection/datum. How can I import a raster map while reprojectingit? Just suppose I am connected on C:\Mydata\GrassSDB_target newLocation PERMANENT and the source data are on C:\users\GrassSDB_source newLocation PERMANENT mymap_UTM_sad69 I tryed something like r.proj input=mymap_UTM_sad69 location=newLocation mapset=PERMANENT dbase=C:\users\GrassSDB_source But I get ERRO:Mapset PERMANENT in input location newLocation - not found Any idea of how to solve this? r.proj tests first for the same location name, which in your case is identical. Try to rename one and only one of your newLocation locations, then try r.proj again. This is a bug in r.proj. HTH, Markus M -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] to move or not move ? grass64 - ? - 7.0
Dear all I have used grass gis 6.4 and I am happy with this version. But I was wondering which are the advantages or disadvantages to move or not grass 7.0 Any thoughts are very welcome. cheers -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] to move or not move ? grass64 - ? - 7.0
Thanks Rashad, I understand that grass 7 have a lot of new features (or improved/fixed ones), including the powerful contribution of python. But I would like to hear from grass7 users about how stable are this version. Anyone survive using only grass7 our we need to switch between them? By the way, a grass GIS DB/Mapset can be used for both versions ? Regards miltinho 2013/1/10, Rashad M mohammedrasha...@gmail.com: This can help you to make a decision[1]. [1] http://trac.osgeo.org/grass/wiki/Grass7/NewFeatures On Thu, Jan 10, 2013 at 11:13 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Dear all I have used grass gis 6.4 and I am happy with this version. But I was wondering which are the advantages or disadvantages to move or not grass 7.0 Any thoughts are very welcome. cheers -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Regards, Rashad -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Multi-path corridors in grass?
Dear all, Anyone using GRASS to map multi-path corridors between paired or multiple habitat patches in landscape ecology studies? best -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.drain problem when flow reach the border
Hi Markus Thanks for all your help. Please, which svn command I need to run to get the daily trunk subversion? best miltinho 2012/12/31, Markus Metz markus.metz.gisw...@gmail.com: Markus Metz wrote: I am going to backport r.cost/r.walk/r.drain from 6.5 to 6.4. Done in r54470-1. Happy new year, Markus M best wishes miltinho 2012/12/30, Markus Metz markus.metz.gisw...@gmail.com: On Sun, Dec 30, 2012 at 8:59 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Hi Rashad I tested with the trunk svn, but apparently it still not worked to me. Are you able to test with my datased? https://www.yousendit.com/download/UW16RE9zTkxlM1JBSXNUQw What are the commands you used? I found r.drain input=custo_aux_cost output=custo_aux_cost_drain vector_points=map_target_pt You need to also provide a direction map with the indir option of r.drain. Markus M Thanks! miltinho 2012/12/27, Mohammed Rashad mohammedrasha...@gmail.com: Hi Milton, I think this fix will solve the problem[0]. For elevation data there might be cases where no path exists ie, only single pixel is the selected path. I had tested it by myself and it worked. Please apply the patch and let me know. If any problem comes in compilation or doing r.drain let me know. It will be great if you can share the test data you are using in case not working after applying the patch [0] http://trac.osgeo.org/grass/attachment/ticket/1840/r.drain.diff [1] http://trac.osgeo.org/grass/attachment/ticket/1840/drain.c On Thu, Dec 27, 2012 at 2:31 AM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Hi Mohamed and Markus Mohamed, thanks for your reply, I will wait for a solution within GRASS :-) Markus, thanks for your reply,. Sorry but I disagree a tiny bit, because is expected that between a source and target I will find a solution with the Least Cost path. Do you think that is there any situation where a least cost between source and target will not be found? Maybe if other software provide the soluction, although is not the bestest option, I prefere receive one complete output than a broken one. Happy new year! miltinho 2012/12/23, Markus Metz markus.metz.gisw...@gmail.com: On Sun, Dec 23, 2012 at 7:07 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Dear all, I am running a least cost path analysis combining r.cost and r.drain. But on some cases, when the least cost path between source and target reach the border of my map, r.drain stop on that position, and target point isn´t reached. This is correct because at the borders, the direction can not be determined, the path could lead out of the current region, you can not tell. Running on other softwares the algorithms the source and target is connected using edges without problem. If in doubt, I would call it a bug in that other software. Any hint are welcome. Try to expand the current region Markus M -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Regards, Rashad -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http
Re: [GRASS-user] r.drain problem when flow reach the border
Hi all, thanks again for the help with this. I tried several versions of 6.4X (6.42 weekly snapshot, 6.43.. ) and 6.5 but in none of them, when I try to see the parameters for create output direction on r.cost man page I can't see how can I generate direction for input on r.drain. Including the suggestion of M.Neterler svn checkout https://svn.osgeo.org/grass/grass/branches/develbranch_6 grass6_devel Sorry to do this so trouble, but maybe I am not finding where can I get an updated version of r.cost, r.drain and r.walk. Happy 2013! miltinho 2012/12/31, Markus Neteler nete...@osgeo.org: On Mon, Dec 31, 2012 at 3:05 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Hi Markus Thanks for all your help. Please, which svn command I need to run to get the daily trunk subversion? See here: http://trac.osgeo.org/grass/wiki/Release/6.4.3RC2-News -- SVN Checkout latest 6.4 release branch (including unreleased bugfixes): Best MarkusN -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.drain problem when flow reach the border
Hi there Everything worked fine. Congrats for all your great job on GRASS-GIS solution, and again many thanks for the effert on deal with this in the transition 2012/2013 (yes! 2013 in Australia and other places around world!). best wishes miltinho 2012/12/31, Markus Neteler nete...@osgeo.org: On Mon, Dec 31, 2012 at 4:41 PM, Markus Metz markus.metz.gisw...@gmail.com wrote: see http://grass.osgeo.org/grass65/manuals/r.cost.html http://grass.osgeo.org/grass65/manuals/r.walk.html For r.drain, the option to use directions as input is indir, see http://grass.osgeo.org/grass65/manuals/r.drain.html Just FYI: I forced the creation of the 6.4 manual an hour ago, so that's also already updated (usually it happens on Saturdays). MarkusN -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.drain problem when flow reach the border
Hi Rashad I tested with the trunk svn, but apparently it still not worked to me. Are you able to test with my datased? https://www.yousendit.com/download/UW16RE9zTkxlM1JBSXNUQw Thanks! miltinho 2012/12/27, Mohammed Rashad mohammedrasha...@gmail.com: Hi Milton, I think this fix will solve the problem[0]. For elevation data there might be cases where no path exists ie, only single pixel is the selected path. I had tested it by myself and it worked. Please apply the patch and let me know. If any problem comes in compilation or doing r.drain let me know. It will be great if you can share the test data you are using in case not working after applying the patch [0] http://trac.osgeo.org/grass/attachment/ticket/1840/r.drain.diff [1] http://trac.osgeo.org/grass/attachment/ticket/1840/drain.c On Thu, Dec 27, 2012 at 2:31 AM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Hi Mohamed and Markus Mohamed, thanks for your reply, I will wait for a solution within GRASS :-) Markus, thanks for your reply,. Sorry but I disagree a tiny bit, because is expected that between a source and target I will find a solution with the Least Cost path. Do you think that is there any situation where a least cost between source and target will not be found? Maybe if other software provide the soluction, although is not the bestest option, I prefere receive one complete output than a broken one. Happy new year! miltinho 2012/12/23, Markus Metz markus.metz.gisw...@gmail.com: On Sun, Dec 23, 2012 at 7:07 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Dear all, I am running a least cost path analysis combining r.cost and r.drain. But on some cases, when the least cost path between source and target reach the border of my map, r.drain stop on that position, and target point isn´t reached. This is correct because at the borders, the direction can not be determined, the path could lead out of the current region, you can not tell. Running on other softwares the algorithms the source and target is connected using edges without problem. If in doubt, I would call it a bug in that other software. Any hint are welcome. Try to expand the current region Markus M -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Regards, Rashad -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.drain problem when flow reach the border
Hi Markus I am using this commands. r.cost -k input=friction_landuse_roads output=cost_map start_points=map_source_pt stop_points=map_target_pt --o r.drain input=cost_map output=drain_map vector_points=map_target_pt --o I have not a direction map, and as for other source-target-friction maps everything works fine without direction map I thougth that this should work for all. If you have the chance, please run the command and take a look that the drain line finish on right border. best wishes miltinho 2012/12/30, Markus Metz markus.metz.gisw...@gmail.com: On Sun, Dec 30, 2012 at 8:59 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Hi Rashad I tested with the trunk svn, but apparently it still not worked to me. Are you able to test with my datased? https://www.yousendit.com/download/UW16RE9zTkxlM1JBSXNUQw What are the commands you used? I found r.drain input=custo_aux_cost output=custo_aux_cost_drain vector_points=map_target_pt You need to also provide a direction map with the indir option of r.drain. Markus M Thanks! miltinho 2012/12/27, Mohammed Rashad mohammedrasha...@gmail.com: Hi Milton, I think this fix will solve the problem[0]. For elevation data there might be cases where no path exists ie, only single pixel is the selected path. I had tested it by myself and it worked. Please apply the patch and let me know. If any problem comes in compilation or doing r.drain let me know. It will be great if you can share the test data you are using in case not working after applying the patch [0] http://trac.osgeo.org/grass/attachment/ticket/1840/r.drain.diff [1] http://trac.osgeo.org/grass/attachment/ticket/1840/drain.c On Thu, Dec 27, 2012 at 2:31 AM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Hi Mohamed and Markus Mohamed, thanks for your reply, I will wait for a solution within GRASS :-) Markus, thanks for your reply,. Sorry but I disagree a tiny bit, because is expected that between a source and target I will find a solution with the Least Cost path. Do you think that is there any situation where a least cost between source and target will not be found? Maybe if other software provide the soluction, although is not the bestest option, I prefere receive one complete output than a broken one. Happy new year! miltinho 2012/12/23, Markus Metz markus.metz.gisw...@gmail.com: On Sun, Dec 23, 2012 at 7:07 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Dear all, I am running a least cost path analysis combining r.cost and r.drain. But on some cases, when the least cost path between source and target reach the border of my map, r.drain stop on that position, and target point isn´t reached. This is correct because at the borders, the direction can not be determined, the path could lead out of the current region, you can not tell. Running on other softwares the algorithms the source and target is connected using edges without problem. If in doubt, I would call it a bug in that other software. Any hint are welcome. Try to expand the current region Markus M -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Regards, Rashad -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http
[GRASS-user] recompile grass 6.42 within window / Msys
Dear all, I need to recompile current grass 6.42 within windows. Where can I find the tips to deal with this? -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.drain problem when flow reach the border
Hi Markus You are right. Thanks for your comments. In fact I don´t thougth on null cell between Source and Target, and agree that extent influences one archive or not a real least cost path solution. Cheers miltinho 2012/12/27, Markus Metz markus.metz.gisw...@gmail.com: On Wed, Dec 26, 2012 at 10:01 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Hi Mohamed and Markus Mohamed, thanks for your reply, I will wait for a solution within GRASS :-) Markus, thanks for your reply,. Sorry but I disagree a tiny bit, because is expected that between a source and target I will find a solution with the Least Cost path. Do you think that is there any situation where a least cost between source and target will not be found? Maybe if other software provide the soluction, although is not the bestest option, I prefere receive one complete output than a broken one. Ah, you are right, if the input is output from r.cost, the cost directions should be used. Depending on the input data and the region extents, the least cost path found may not be the real least cost path, but there should be one solution. There are situations where there is no connection between source and target if source and target are separated by NULL cells and no NULL cost is specified to r.cost. Markus M Happy new year! miltinho 2012/12/23, Markus Metz markus.metz.gisw...@gmail.com: On Sun, Dec 23, 2012 at 7:07 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Dear all, I am running a least cost path analysis combining r.cost and r.drain. But on some cases, when the least cost path between source and target reach the border of my map, r.drain stop on that position, and target point isn´t reached. This is correct because at the borders, the direction can not be determined, the path could lead out of the current region, you can not tell. Running on other softwares the algorithms the source and target is connected using edges without problem. If in doubt, I would call it a bug in that other software. Any hint are welcome. Try to expand the current region Markus M -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.drain problem when flow reach the border
Hi Mohamed and Markus Mohamed, thanks for your reply, I will wait for a solution within GRASS :-) Markus, thanks for your reply,. Sorry but I disagree a tiny bit, because is expected that between a source and target I will find a solution with the Least Cost path. Do you think that is there any situation where a least cost between source and target will not be found? Maybe if other software provide the soluction, although is not the bestest option, I prefere receive one complete output than a broken one. Happy new year! miltinho 2012/12/23, Markus Metz markus.metz.gisw...@gmail.com: On Sun, Dec 23, 2012 at 7:07 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Dear all, I am running a least cost path analysis combining r.cost and r.drain. But on some cases, when the least cost path between source and target reach the border of my map, r.drain stop on that position, and target point isn´t reached. This is correct because at the borders, the direction can not be determined, the path could lead out of the current region, you can not tell. Running on other softwares the algorithms the source and target is connected using edges without problem. If in doubt, I would call it a bug in that other software. Any hint are welcome. Try to expand the current region Markus M -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.drain problem when flow reach the border
Dear all, I am running a least cost path analysis combining r.cost and r.drain. But on some cases, when the least cost path between source and target reach the border of my map, r.drain stop on that position, and target point isn´t reached. Running on other softwares the algorithms the source and target is connected using edges without problem. Any hint are welcome. Merry Xmas -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] WinGrass under MSYS
Great, Thanks milton 2012/10/20, Martin Landa landa.mar...@gmail.com: Hi, 2012/10/20 Milton Cezar Ribeiro miltinho.astrona...@gmail.com: But I noticed that on the recent versions of Standalone WinGRASS, althoug MSYS is installed, we loose the capability of copy and past sets of grass command (like was in GRASS 6.4.2RC2 / MSYS). Any idea of how can I run Wingrass within Msys environment again? you didn't loose any capability, go to Start menu GRASS GIS version | GRASS GIS version GUI with MSYS You can add it to your desktop on your own if you prefer. Martin -- Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] WinGrass under MSYS
On the previous version of Standalone WinGRASS, when installing I always got a icon on my desktop that allowed me to start grass inside Msys, so I was able to run scripts by copy and pasting several command lines directly on Msys console. But I noticed that on the recent versions of Standalone WinGRASS, althoug MSYS is installed, we loose the capability of copy and past sets of grass command (like was in GRASS 6.4.2RC2 / MSYS). Any idea of how can I run Wingrass within Msys environment again? Best -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 Depto Ecologia http://www.rc.unesp.br/ib/ecologia/ PG ECO BIODIV http://www.rc.unesp.br/ib/ecologia/posbiodiversidade/index.php CV http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4792988H6mostrarNroCitacoesISI=truemostrarNroCitacoesScopus=true Google citations http://scholar.google.com/citations?user=OWX_2eAJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] unsupervised image classification
Excelent! thanks miltinho 2012/3/28, Markus Neteler nete...@osgeo.org: On Tue, Mar 27, 2012 at 11:03 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Dear all, I need to classify a Landsat image on GRASS but I would like to automatically define the classes, not to define a priori. Is there a way of I set the input bands on a processing, and the algorithm find the signatures within the images and define the number of classes? Yes, please check http://grass.osgeo.org/wiki/Image_classification#Unsupervised_classification Markus -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil http://www.rc.unesp.br/ib/ecologia/ Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] unsupervised image classification
Dear all, I need to classify a Landsat image on GRASS but I would like to automatically define the classes, not to define a priori. Is there a way of I set the input bands on a processing, and the algorithm find the signatures within the images and define the number of classes? In this moment I am running I.SMAP, which maybe will produce an interesting result, but as I am using multi-temporal bands (for 3 different decades) it is not an easy task to defini the multi-transitional classes. Thanks in advance -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil http://www.rc.unesp.br/ib/ecologia/ Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.out.gdal : exporting category in HFA format
Dear all, I have a raster map with its own category information. Now I am using r.out.gdal to export to HFA (.img) format. But I would like to add the labels (category) in the output file, intead of have only VALUEs. Any idea? best -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil http://www.rc.unesp.br/ib/ecologia/ Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] get values for each point within lines
Dear Grass gurus I have a vector map with lines, and I need to get raster values each 5 meters. So I used v.split to spring my lines in segments of 5 meters, and after that I used v.to.points to generate points on each node of my line segments. Now I am trying to overlap my points against a elevation map, but when I run v.what.rast vect=lines_split_05m_points rast=dem col=elevation the numbers of records on my v.db.select is the same of my lines instead of my points. The question is how to get the elevation values and keep the other attributes from my source line. best wishes -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil http://www.rc.unesp.br/ib/ecologia/ Fone: +55 19 3526-9647 (office) 19 3526-9680 (lab) Cel: 19 9853-3220 / 19 9853-5430 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] zoom errors in grass WinGRASS-6.4.2RC2-1-Setup
Dear all, I am running grass 6.4.2 which I just downloaded and installed from OSGEO server. When I try to zoon in or out, I get the following error under Windows/7 64bits Traceback (most recent call last): File c:/osgeo4w/usr/src/grass-6.4.2RC2/dist.i686-pc- mingw32/etc/wxpython/gui_modules/mapdisp_window.py, line 109, in OnMotion KeyError : 'proj' In fact, I had the same error with previously installed WinGRASS-6.4.2RC1-1-Setup version. Any hint? cheers -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil http://www.rc.unesp.br/ib/ecologia/ Fone: +55 19 3526-9647 Cel: 19 9853-3220 / 19 9853-5430 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.to.rast problem
Hola José, Send also db.describe of your map and the complete command of v.to.rast that you are using. Maybe you need to use something like: v.to.rast input=yourmap output=yourmap_rast use=attr column=PUT_A_NUMERIC_COLUMN best wishes miltinho 2011/12/18, José Carlos Guerrero Antúnez jcgantu...@gmail.com: Hi, I am working with GRASS 6.4.1(2011) on a OSX 10.6.8 (MaBookPro). When I do the process for importing with v.to.rast, According to different Source of raster values: - attr Column parameter missing (or use value parameter)) - cat (Loading data, Reading areas, Reading features, Writing raster map...) the process becomes eternal - val (loading data, converted areas: 1 of 1, Converted points/lines: 0 of 0, t.to.rast complete, Reading features, Writing raster map,...) The process becomes eternal - z (Loading data) Vector map Uruguay@PERMANENT is not 3D - dir (Loading data, Reading features, Writing raster map, Converted areas: 0 of 1, Converted pointes/lines: 0 of 0, v.to.rast complete) It is not created. Vector file data are showed below: ++ | Layer: Uruguay@PERMANENT | | Mapset: PERMANENT | | Location:Sudamérica | | Database:/Users/joseca/grass | | Title: | | Map scale: 1:1 | | Map format: native | | Name of creator: joseca | | Organization: | | Source date: Sun Dec 18 19:31:15 2011 | || | Type of Map: vector (level: 2) | | | | Number of points: 0 Number of areas: 1 | | Number of lines:0 Number of islands:1 | | Number of boundaries: 1 Number of faces: 0 | | Number of centroids:1 Number of kernels:0 | | | | Map is 3D: No | | Number of dblinks: 1 | | | | Projection: x,y | | N: 30:06:02.013245SS: 34:58:40.005798S | | E: 53:05:42.004395WW: 58:26:19.037476W | | | | Digitization threshold: 0 | | Comments: | | | ++ Thank you, José Guerrero -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil http://www.rc.unesp.br/ib/ecologia/ Fone: +55 19 3526-9647 Cel: 19 9853-3220 / 19 9853-5430 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] delineate Watersheds
Hi Ruth, Try g.region rast=YourMap r.terraflow... good look milton 2011/11/11, Ruth Sos Del Diego rsosdeldi...@googlemail.com: Hello list, I trying to delineate watersheds in a STRM-data(DGM.tif) | Data Type:CELL | | Rows: 1809 | | Columns: 1413 | | Total Cells: 2556117 I´m trying with r.terraflow and the error is: Cell file DGM.tif resolution differs from current region and too I try with r.watershed and finally the run the basin is empty. (0.0) I´m using the nightly build- SEXTANTE (gvSig 1.11) Any idea what is wrong? Any suggestions are welcome! Cheers! Ruth -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil http://www.rc.unesp.br/ib/ecologia/ Fone: +55 19 3526-9647 Cel: 19 9853-3220 / 19 9853-5430 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] removing empty files from grass mapset
Dear all, I just remove 5,000 maps from a mapset using g.remove -f map=mymap, and observed that, although the command worked fine, when I go under the directories (like .../Mapset/cell) I see there a lot of empty files (0 bytes), and the removed raster map filenames are still there. How can I remove that? best wishes milton -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil http://www.rc.unesp.br/ib/ecologia/ Fone: +55 19 3526-9647 Cel: 19 9853-3220 / 19 9853-5430 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] g.remove from another mapset
Dear all, I need to remove raster maps from another mapset instead PERMANENT As I am putting this in a python-grass script, I tried something like g.remove -f rast=mymap@OTHERMAPSET but I get the error message that the raster map is not in the current mapset (PERMANENT). Any suggestion? best wishes -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil http://www.rc.unesp.br/ib/ecologia/ Fone: +55 19 3526-9647 Cel: 19 9853-3220 / 19 9853-5430 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] g.remove from another mapset
Hi Martin, Worked fine! best milton 2011/10/28, Martin Landa landa.mar...@gmail.com: Hi, 2011/10/28 Milton Cezar Ribeiro miltinho.astrona...@gmail.com: I need to remove raster maps from another mapset instead PERMANENT As I am putting this in a python-grass script, I tried something like g.remove -f rast=mymap@OTHERMAPSET but I get the error message that the raster map is not in the current mapset (PERMANENT). you are allowed to modify data only from your current mapset. In other words you cannot delete a map from another mapsets. So first you need to switch to OTHERMAPSET g.mapset OTHERMAPSET and then to delete the map g.remove -f rast=mymap Martin -- Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil http://www.rc.unesp.br/ib/ecologia/ Fone: +55 19 3526-9647 Cel: 19 9853-3220 / 19 9853-5430 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Linear regression between 2 rasters (old and new)
Dear Antonio, Depending on your raster sizes you can load both maps into R and compute the stats over there. best miltinho 2011/9/30, António Rocha antonio.ro...@deimos.com.pt: I I have been using r.series to calculate a linear regression slope for an old and new raster map using r.series. But now I need to also obtain the error. It seems not to be possible to use r.series. Is this true? if yes, what alternative do I have? THanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6506 (20110930) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Miltinho - m...@rc.unesp.br Laboratório de Ecologia Espacial e Conservação - LEEC Depto de Ecologia - UNESP - Rio Claro Av. 24A, 1515- Bela Vista 13506-900 Rio Claro, SP, Brasil http://www.rc.unesp.br/ib/ecologia/ Fone: +55 19 3526-9647 Cel: 19 9853-3220 / 19 9853-5430 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Fwd: Area Extraction tools
Hi Yada, If you have a raster map with your area of interest you can use r.mask ... r.mapcalc DEM_clip=DEM or g.region rast=mask_file g.region res=DEFINE DESIRED RESOLUTION r.mapcalc DEM_clip=DEM I hope this helps. milton 2011/7/27 Yada Elisabeth yadaelisab...@gmail.com Dear GRASS users, I download a DEM file( SRTM geotiff) that I opened in GRASS 6.4.1 One problem I am facing is that I want to extract a part (my area of interest) of this DEM I don't know how to proceed Yada Thank You in advance! ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Miminum Spanning-Tree (MST)
Dear all, Is there a way of I obtaion a polygon- or raster-based Minimum Spanning Tree (MST) on GRASS? It is like to join every polygons to the neighbours using edges (linear links) without allow the edges overlap each other. cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] voronoi tesselation for vector/raster patches
Dear all, Can I generate voronoi tesselation for polygons (instead of points) of raster-based patches? cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: raster stats within vector
Hi Paulo, Thanks a lot. cheers milton 2011/7/17 Paulo van Breugel p.vanbreu...@gmail.com ** v.rast.stats (menu: Vector-Update area attribute from raster) will calculate n, min, max, range, mean, stddev, variance, coeff_var, sum and optionally the 1st quartile, median, 3rd quartile, and given percentile. Columns will be created on the fly, you only have to give the 'base name' of the columns. I.e., no need to create a column yourself. On 07/16/2011 10:14 PM, Milton Cezar Ribeiro wrote: BTW, my vector maps have polygons. milton 2011/7/16 Milton Cezar Ribeiro miltinho.astrona...@gmail.com Dear all, I have a raster map with continuous values, and need to extract some stats (like mean, std etc) assigning these zonal values as a vector column. I just added the empty columns on the my vector map using : v.db.addcol map=frag_ts columns=FLOW_STD DOUBLE PRECISION How can I do this on grass? cheers milton ___ grass-user mailing listgrass-user@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.rast.stats stucked
Dear all, I have a vector map with only 7 features (polygons). I use the command below to extract raster info. v.rast.stats vector=my_polygons raster=my_raster colprefix=LG As my polygons have holes, I receive several messages as the below: Area without centroid (may be OK for island) But I don´t know why grass get stucked and don´t finish the process. By the way, I am using MSYS grass 6.4.1 version under windows. cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: raster stats within vector
BTW, my vector maps have polygons. milton 2011/7/16 Milton Cezar Ribeiro miltinho.astrona...@gmail.com Dear all, I have a raster map with continuous values, and need to extract some stats (like mean, std etc) assigning these zonal values as a vector column. I just added the empty columns on the my vector map using : v.db.addcol map=frag_ts columns=FLOW_STD DOUBLE PRECISION How can I do this on grass? cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] raster stats within vector
Dear all, I have a raster map with continuous values, and need to extract some stats (like mean, std etc) assigning these zonal values as a vector column. I just added the empty columns on the my vector map using : v.db.addcol map=frag_ts columns=FLOW_STD DOUBLE PRECISION How can I do this on grass? cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] calculating perimeter per patch
Hi Markus, Thanks for the reply. In fact a vector-based should solve this. But, unfortulatelly, as I am working with a very large raster map (48,000 x 60,000 pixels; 250,000 clumps), and GRASS fail when I try r.to.vect under Windows/Vista. Is there any other options of convert large raster maps to vector, instead of r.to.vect? cheers milton 2011/2/12 Markus Metz markus.metz.gisw...@googlemail.com On Sat, Feb 12, 2011 at 5:11 AM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Dear all, I have a raster map where each patch (i.e. forest fragment) have different ID. I need to estimate the perimeter per patch, but R.le.patch can't do the job. Have you tried a vector approach? Then the perimeter can be simply uploaded with v.to.db option=perimeter. (Watch out for the fractal problem, both r.le.patch and v.to.db) Markus M GRASS 6.4 r.le.patch map=fragments sam=w per=p1 siz=s1 out=STATS.txt PARAMETER CHOICES: MAP: fragments SAMPLE: whole map TRACING: 8 neighbor SIZE MEASURES: mean patch size PERIMETER MEASURES: sum of perims R.LE.PATCH IS WORKING; ERROR: G_calloc: unable to allocate 24806 * 8 bytes at alloc_cell.c:83 Any hint? bests milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] calculating perimeter per patch
Hi Markus, I will try that. The GRASS_VECTOR_LOWMEM=1 works only for G7? best milton 2011/2/12 Markus Metz markus.metz.gisw...@googlemail.com On Sat, Feb 12, 2011 at 11:58 AM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Hi Markus, Thanks for the reply. In fact a vector-based should solve this. But, unfortulatelly, as I am working with a very large raster map (48,000 x 60,000 pixels; 250,000 clumps), and GRASS fail when I try r.to.vect under Windows/Vista. Is there any other options of convert large raster maps to vector, instead of r.to.vect? No, but you could try a recent version of GRASS 7 and after starting GRASS, set at the commandline export GRASS_VECTOR_LOWMEM=1 then try r.to.vect It will take some time, but it should finish eventually. Markus M Markus Metz: Milton Cezar Ribeiro: Dear all, I have a raster map where each patch (i.e. forest fragment) have different ID. I need to estimate the perimeter per patch, but R.le.patch can't do the job. Have you tried a vector approach? Then the perimeter can be simply uploaded with v.to.db option=perimeter. (Watch out for the fractal problem, both r.le.patch and v.to.db) Markus M GRASS 6.4 r.le.patch map=fragments sam=w per=p1 siz=s1 out=STATS.txt PARAMETER CHOICES: MAP: fragments SAMPLE: whole map TRACING: 8 neighbor SIZE MEASURES: mean patch size PERIMETER MEASURES: sum of perims R.LE.PATCH IS WORKING; ERROR: G_calloc: unable to allocate 24806 * 8 bytes at alloc_cell.c:83 Any hint? bests milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.neighbors - annulus (torus) neighbor
Hi Jarek and r.neighbor maintenators, I suggest put a flag on r.neighbor for annulus (torus), so in future versions the module will be always update. Just a suggestion. Jarek, congrats. cheers milton 2010/11/15 Jasiewicz Jarosław jar...@amu.edu.pl Hi all! while I need annulus (or torus) neighbor to some of my work I added easily it to neighbors (grass 65 delev) module. Because I do not won't to create a fork of that module I put it on that url: http://www.staff.amu.edu.pl/~jarekj/r.neighbors2.tar.gzhttp://www.staff.amu.edu.pl/%7Ejarekj/r.neighbors2.tar.gz Annulus neighbor is a circular neighbor but with exuding nearest neighborhood defined by inner radius. Maybe someones will find this useful. regards Jarek ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.stream : vetor of horton order
Dear all, I extracted the drainage and direction from a DEM using r.stream.extract. r.stream.extract --o elevation=aster_dem_wgs84_UTM threshold=50 stream_vect=aster_dem_wgs84_UTM_vect_th00050 stream_rast=aster_dem_wgs84_UTM_rast_th00050 direction=aster_dem_wgs84_UTM_rast_th00050_direction How can I obtain a vector of Horton order?! cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] grass 64b in windows7
Dear all, How can I have a 64bit windows grass under a 6Gb Windows7 dell machine? I am using OSGeo4W/MSYS version, but many of my processing I get memory errors. cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problems starting Grass6.3 on Windows
Hi Vishal, Is there a why of you use grass 6.3? How about you try grass 6.4R6? Have you tryed grass with a sample (Northcaroline or Spearfish) mapset? cheers milton 2010/5/11 Vishal Mehta vishalm1...@gmail.com I'm having some recent trouble getting Grass to start on Windows XP. when i try to start Grass i get the following message: error setting region (Problem with g.region): child killed:unknown signal and then a message box with lots of errors starting with cant read parts(w): no such variable ... ... I've tried and failed with: - creating a new grass database directory without spaces in the path (the default is in My Documents but it used to work before); - unistalling and re-installing Grass - found some reference to gis_mlock in the archives but did not understand a word of it... Other info: - I have Grass 6.3.0 - Used to work until about a week ago - Only difference since then is that i've installed QGIS , also postgis shaepefile loader (not sure if that would cause the trouble, perhaps a different gdal version or something??). It would be nice to fix this so that I can use Grass commands from QGIS.. Any tips? Thanks, Vishal -- Vishal K. Mehta, PhD Scientist Stockholm Environment Institute - US 133 D St Suite F Davis CA 95616 www.sei-us.org ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Spgrass6 and R using GRASS Bash Scripts
Hi Gilbert, There ist more than one way. I will tell you the way that works for me: 0. Be sure your Rgui.exe are in the windows path like C:\Program Files (x86)\R\R-2.9.0\bin 1. Start grass/MSYS (sorry, I am on windows). 2. Connect to you mapseet 3. inside MSYS shell type: Rgui.exe 4. Inside R you can run your scripts like: #this is on R require(spgrass6) gmeta6() myMap=readRAST6(...) summary(myMap) system(g.list.exe rast) -- But I know that the way I do are not updated :-) cheers milton 2010/4/16 Gilbert Ferrara gilbertferrara1...@gmail.com Good Morning I need to do something similar to what Neteler/Mitasova do in their book with spgrass6 and R. Since I will have to perform this a couple of times I want to build a simple Script with input/output GUI. But I have never seen a GRASS Bash script that calls R and spgrass6. So my question is? Does anyone have/knows one? If not, I just need to send to the command line exacly the same commands? e.g. R , summary(data) ) Thanks Gilbert ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: r.mapcalc and nested if
Hi there, may be the last command need double == r.mapclac 'final=if(first_map==1,0,first_map' cheers milton 2010/4/14 Jarek Jasiewicz jar...@amu.edu.pl thedok78 pisze: Sorry, this is the expression: r.mapcalc 'map2=if(map ==1 ,if(mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1])==1,0,mode(map[-1,-1],map[-1,0],map[-1,1],map[0,1],map[1,1],map[1,0],map[1,-1],map[0,-1]),map))' Thanks Luca seem that much easier will be: r.neighbors input=your_map method=mode size=3 output=mode_output r.mapcalc 'first_map=if(your_map==1,output_mode,your_map)' r.mapclac 'final=if(first_map=1,0,first_map' ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Running Python scripts inside GRASS
Hi There, I am happy running python inside a OSGeo4W grass instalallation. I start MSYS, start grass, and then ryn python 2.5 which are run from C:\OSGeo4W\apps\Python25. But I needed to copy the grass/python library to C:\OSGeo4W\apps\Python25\Lib\site-packages. Good luck, milton 2010/3/9 Hamish hamis...@yahoo.com Sonja wrote: Is there any other place where a wrong PYTHONPATH could be stored??? if you have the grass source code try from the the msys prompt: grep -ri python25 mswindows/* wingrass is a bit of a tangle of startup scripts, especially the osgeo4w build. the above searches through the custom/lesser known wingrass ones. you can also try a recursive grep in the GRASS-64-SVN\ directory to try and find and. If doing that I'd add the -I flag to skip searching in binaries. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Region Definition in GRASS
Pedro, If you run a r.mapcalc, the extent and resolution will follow the current region. So, check g.region -p before run r.mapcalc. You can use g.region rast=MapLandsat to set the current region to use MapLandsat dimensions as reference. After you can reset your extention and resolution using g.region -d . By the way, read g.region --help. best milton 2010/2/5 Pedro Roma pedroroma1...@gmail.com Ok but for instance. If my region is defined with a Spatial Resolution of 1kilometer and, if I do some r.mapcalc with two LANDSAT images, the output image will have which one of the resolutions? the Region or the 2 landsat images? Thank you Glynn for the overview. And thank you Hamish for the weblinks On Thu, Feb 4, 2010 at 7:14 PM, Glynn Clements gl...@gclements.plus.comwrote: Pedro Roma wrote: I'm a kind of a new GRASS user and I'm a bit lost on regarding something. REGION (or spatial extent). For instance the Location North-Carolina. When the location was created a certain extent, projkection a Spatial Resolution was defined. But, inside this location there is a lot of data with different extent and regional cover. And, inside GRASS I can re-define region properties (g.region) and use a certain data to define region. And, it's possible to have a region with a small extent than a vectorial that exists in the Location. My question is, what is the region and what is its purpose? Because I don't see any impact of region over my image/data processing or image visualization. Or is it? The region is mostly used for raster processing and graphical display. Most commands which read or write raster maps will use the bounds and resolution from the current region. Any input maps will be cropped, padded and/or resampled according to the region settings. Display commands translate geographic or cartographic coordinates to display coordinates such that the current region just fits inside the display window. Vector commands are normally unaffected by the current region. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Raster Map subset
Hi Stephen, I suggest you use spgrass6 package on R, set your g.region using system(g.region n= s= w= e=) on R and use readRAST6 to load your map as SGDF. Of course, if you are under windows, you know that depending the size of your grid you will get troubles :-) best milton 2010/2/5 stephen sefick ssef...@gmail.com Is there a way to get a subset of a larger raster map. I have zoomed into an area and would like to make a raster map only containing that area. Can I do this? I am going to use this as input into R for some topmodel runs. thanks, -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Interpreting Terrestrial Geomorphometry
Hi There, I put my 2cents on the Markus' suggestion. This book changed my life :-) Thanks Tom! best regards milton 2010/1/27 Rich Shepard rshep...@appl-ecosys.com On Wed, 27 Jan 2010, Markus Neteler wrote: Let me recommend this book: http://geomorphometry.org/book Thank you, Markus. I'll ask my local library to borrow a copy for me. (of course it has a GRASS chapter :) Well, duh! GRASS grows on soils doesn't it? Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Problems Running r.flow
Hi Rich, try g.region nsres=5 ewres=5, may be you can bypass this. best regards milton 2010/1/27 Rich Shepard rshep...@appl-ecosys.com In -6.4svn from last Friday's weekly snapshot. After setting g.region rast=aber5m I tried running r.flow elevin=aber5m aspin=aber5aspect flout=aberFlowline lgout=aberFlowlength and ran into resolution problems: Reading input files: elevation ERROR: Elevation file's resolution differs from current region resolution g.region -p reports: nsres: 4.99982001 ewres: 4.99979979 g.region rast=aber5m -p reports: nsres: 4.99982001 ewres: 4.99979979 I also get the same error message about the aspect map (derived from the elevation map) which has the same resolution as above. Any ideas what's going on? More importantly, how do I get his module running here? Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Problems Running r.flow
*Hum* Interesting -a flag. But g.region res=5 and g.region nsres=5 ewres=5 produce the same result? I never tryed. Rich, another way is you use n= s= w= e= options, that may be will result something like -a, I suppose. cheers milton 2010/1/27 Rich Shepard rshep...@appl-ecosys.com On Wed, 27 Jan 2010, Milton Cezar Ribeiro wrote: try g.region nsres=5 ewres=5, may be you can bypass this. milton, That's how I defined the resolution, but it was returned as slightly less for some reason. See: GRASS 6.4.0svn (Oregon):/usr4/grassbase g.region res=5 -p projection: 99 (Lambert Conformal Conic) zone: 0 datum: nad83 ellipsoid: grs80 north: 1334419.25160578 south: 1279151.24118496 west: 769192.9282895 east: 819255.9236 nsres: 4.99982001 ewres: 4.99979979 rows: 11054 cols: 10013 cells: 110683702 Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Problems Running r.flow
Rich, If your are using a resolution of ~5m, I suggest you adjust your region ~1m. I think that the final result will not suffer so much changes. If so, may be you can try finer resolution (4m?). g.region n=1334420 s=1279150 w=769190 e=819255 res=5 good luck milton 2010/1/27 Rich Shepard rshep...@appl-ecosys.com On Wed, 27 Jan 2010, Milton Cezar Ribeiro wrote: But g.region res=5 and g.region nsres=5 ewres=5 produce the same result? I never tryed. Milton, Yes. res is equal to setting both nsres and weres to the same values. Rich, another way is you use n= s= w= e= options, that may be will result something like -a, I suppose. That changes the bounds. I don't want to do that other than the small increment necessary to get an integer resolution value. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Problem with r.watershed in Spearfish
Margherita, Check with g.region -p and see how many 100x100 pixels you have. May be it is very small for your threshold. best milton 2010/1/25 Margherita Di Leo direg...@gmail.com Hi list! I often used r.watershed without problems, but today i had a very strange one. Please can anyone verify that? I used Spearfish mapset, map: elevation.10m: g.region res=100 r.watershed elevation=elevation@permanent threshold=100 accumulation=acc_10 drainage=dir_10 stream=str_10 convergence=5 memory=300 The stream map has many streams with 0 value. The problem is only with res=100. Does r.watershed fail? Thank you -- Eng. Margherita Di Leo Ph.D. Candidate Methods and Technologies for Environmental Monitoring Department of Environmental Engineering and Physics (DIFA) University of Basilicata Campus Macchia Romana 85100 - Potenza Italy Office: +39-0971205363 Fax: +39-0971205160 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] forcing raster accept only integer 0
Hi Glynn, Thanks for the reply. By the way, where I can have access to the list of datatypes for a raster map? I ran on my raster map and got datatype=CELL. I was expecting to get byte,float,integer... bests milton 2010/1/22 Glynn Clements gl...@gclements.plus.com Milton Cezar Ribeiro wrote: I am trying to work around with a problem of map size under Vista (2G). I was wondering if I use the command below I will get a raster map with values integer and =0. r.mapcalc Dem.new=if(Dem=0,int(Dem),null()) That will produce an integer map with no negative values. As my original map have values positive and negative, as well as it is fload, my question is if are there a way of I force my output map accept only integers=0. By the way, is there a way of I know if a map is float or integer? r.info, or r.info -t if you only want the type. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] change the cell value of a geotiff
Hey someone, don't push.. r.mapcalc NewMap=if(OldMap==999,null(),OldMap) bests milton 2010/1/20 maven apache apachemav...@gmail.com Hi: I have a geotiff now, and some cells have a not-expected value, so I want to Interpolation them, and I notice that there is a command named r.fillnull, however the cell to be changed is not null, so I want to : 1) change the special cell value to null 2) use the r.fillnull to Interpolation. The problem is I do not know how to do the step 1 with grass, anyone can do me a favor? Thanks. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] change the cell value of a geotiff
Hi Maven, do you think that if you have min = -1267650600228229401496703205376.00 is the same value you used on r.null setnull=-1.26765060022823e+030 ? I think following the Hamish's suggestion you will get there :-) bests milton 2010/1/20 maven apache apachemav...@gmail.com 2010/1/21 Hamish hamis...@yahoo.com maven apache wrote: I have a geotiff now, and some cells have a not-expected value, so I want to Interpolation them, and I notice that there is a command named r.fillnull, however the cell to be changed is not null, so I want to : 1) change the special cell value to null 2) use the r.fillnull to Interpolation. The problem is I do not know how to do the step 1 with grass, anyone can do me a favor? try r.null. Hamish Hi; It seems that the r.null does not work. 1) I use the r.info to get the info of the tiff which I have imported to grass and I got that: Hi: (Thu Jan 21 09:56:46 2010) r.info map=...@permanent ++ | Layer:a...@permanent Date: Tue Jan 12 16:02:39 2010| | Mapset: PERMANENT Login of Creator: Administrator | | Location: hdf5loc | | DataBase: D:\GrassDataBase | | Title: ( aaa ) | | Timestamp: none| || | | | Type of Map: raster Number of Categories: 499 | | Data Type: FCELL | | Rows: 720| | Columns: 1440 | | Total Cells: 1036800| |Projection: Latitude-Longitude | |N:90NS:90S Res: 0:15 | |E: 180EW: 180W Res: 0:15 | | Range of data:min = -1267650600228229401496703205376.00 max = 4 | | | | Data Description:| |generated by r.in.gdal | | | | Comments:| |r.in.gdal input=D:\data\hdf5\aaa.tif output=aaa | | | ++ (Thu Jan 21 09:56:47 2010) Command finished (0 sec) (Thu Jan 21 09:56:50 2010) 2)then I use r.null to set the null value to null: r.null setnull=-1.26765060022823e+030 3) I use the r.info again to see if it is changed : r.info -r map=...@permanent min=-1.26765060022823e+030 max=498.5 (Thu Jan 21 09:56:50 2010) Command finished (0 sec) It seems that the valuse are not changed to null. I wonder why? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] how to know if LFS flag is enable
Hi Glynn, sorry my delay on answer. So under windows we are really limitted to 2Gb independent of compile using LFS flag? :-( Are there plans to by pass this!? I hope windows users could also be favoured of GRASS for work with large maps on near future. bests milton 2010/1/19 Glynn Clements gl...@gclements.plus.com Milton Cezar Ribeiro wrote: I am using r.series to join several large raster maps, but the south part of my map is completely empty, and whem I try r.what I get a error message can't read. I was wondering if it is a problem with LFS (large file system) flag. I compiled grass 6.4.0svn two weeks ago using OSGeo4W/Msys under Vista, in a dell/64bit computer with 6G ram. Any help are welcome. LFS currently doesn't work on Windows. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: compiling grass 64bit under vista?
Dear Helmut and Markus I checked out the Helmut's suggestion, and in fact are there a MinGW-w64 that is a mix between 32-64b. I think it is a very promising solution for windows-based users, mainly because many of new computers that we buy for GIS end-users area, in general, with large memory (4Gb+) and 64bit. But as Helmut advice, I can get in trouble if we mix 32 64b. Addictionally, I have no idea of how many users compile grass from source under MSYS/Windows using the easy-to-follow and valuable intructions found at http://trac.osgeo.org/grass/wiki/CompileOnWindows I think rare people will be able to bypass the difficulties of try compile grass under windows using other configurations/libs different from the instructions above. So, if someone feel confortable to try, and could prepare a *CompileOnWindowsW64 *it will be great. By the way, the new python-based grass7 will still have the same limitations under windows (not use LFS and only run as 32b)? bests niltinho 2010/1/19 Helmut Kudrnovsky hel...@web.de hi, Dear all, I have used grass under Vista, and compiled it using Msys/OSgeo4W version. If I am not wrong I will get only a 32bit grass, because Msys/gcc and some of libs downloadable from the links http://trac.osgeo.org/grass/wiki/CompileOnWindows deal with 32bit instead 64. Is there a way of I prepare and environment under Vista to compile it with 64b? May be someone could help to preparar a *CompileOnWindows64* :-) Any hint are welcome. there seems to be a working mingw64-environment: http://sourceforge.net/projects/mingw-w64/ The mingw-w64 project is a complete runtime environment for gcc to support binaries native to Windows 64-bit and 32-bit operating systems. maybe mixing 32bit-dependecies from osgeo4w and 64bit GRASS isn't a very good idea, so maybe it would be necessary to build a new 64bit building environment with mingw64 and 64bit dependecies like gdal, proj, tcltk, etc. but it would be a nice idea ... :o) best regards Helmut ___ Preisknaller: WEB.DE http://web.de/ DSL Flatrate für nur 16,99 Euro/mtl.! http://produkte.web.de/go/02/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] how to know if LFS flag is enable
Dear all, I am using r.series to join several large raster maps, but the south part of my map is completely empty, and whem I try r.what I get a error message can't read. I was wondering if it is a problem with LFS (large file system) flag. I compiled grass 6.4.0svn two weeks ago using OSGeo4W/Msys under Vista, in a dell/64bit computer with 6G ram. Any help are welcome. To compile this I used the instructions under http://trac.osgeo.org/grass/wiki/CompileOnWindows I checked out ./mswindows/osgeo4w/package.sh and found a flag --enable-largefile \ but I am not sure if it define that my system support or not LFS. Any hint are welcome milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] compiling grass 64bit under vista?
Dear all, I have used grass under Vista, and compiled it using Msys/OSgeo4W version. If I am not wrong I will get only a 32bit grass, because Msys/gcc and some of libs downloadable from the links http://trac.osgeo.org/grass/wiki/CompileOnWindows deal with 32bit instead 64. Is there a way of I prepare and environment under Vista to compile it with 64b? May be someone could help to preparar a *CompileOnWindows64* :-) Any hint are welcome. bests miltinho ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] exporting to KML
Dear all, I have some vector files generated using r.stream.extract, and now I need to convert it to KML. But my mapset is on Albers projection. Is there a way of I export automatically my vector map using v.out.ogr withour projection transformation? v.out.ogr.exe --o input=myvect dsn=myvect_ln.KML type=line format=KML cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] slope unit in r.param.scale
Dear friends, I generated a map with slopes using: r.param.scale input=srtm output=srtm_SLOPE_ws007 param=slope size=07 The output range from r.info is min = 0.00 max = 3.930307 May I consider that all values 1 is equivalent to 100% (or 45 degrees) of steepnees? Case not, how can I identify the pixels with 100% of steepness? By the way, my maps projection is in meters. cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Clipping Raster Map
Hi Rich, Have you tried r.mask ? bests milton 2010/1/1 Rich Shepard rshep...@appl-ecosys.com While the raster DEM displays quickly in the GUI it is unacceptably slow trying to display in the CLI monitor. Does the large number of colors explain the time difference? GRASS 6.4.0svn (Oregon):/usr4/grassbase/Oregon d.mon x0 using default visual which is TrueColor ncolors: 16777216 Graphics driver [x0] started I want to extract only those sells within the project bounds. I don't see an r.overlay or r.select equivalent to the vector commands. Running 'g.region -c map=dem vector=abernethy' did nothing. I don't know if there's a conflict between the '-c' option (use current region) and the definition of the vector boundary map. When I removed the '-c' option and ran the command the result appears to be the entire original map squeezed into the smaller region's bounds. Not what's needed. Do I need to build a raster version of the vector project area, and use that to mask everything external to zero? I cannot display the DEM using a script because it would take a _very_ long time for that one map to display. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] meaning of g.region zoom?
Hi Paulo, Have you tried g.region -d ? ciao milton 2009/12/30 Paolo Cavallini cavall...@faunalia.it Hi all. If I zoom to a small map full of nulls, with only a non-nulla area on one side, and zoom back to the full extent, I do not see the difference: === g.region zoom=los_175_...@paolo g.region -p projection: 99 (Transverse Mercator) zone: 0 datum: rome40 ellipsoid: international north: 4834913.78536118 south: 4821192.33020234 west: 1668983.89916058 east: 1685595.88381969 nsres: 100.156607 ewres: 100.07219674 rows: 137 cols: 166 cells: 22742 g.region zoom=...@permanent g.region -p projection: 99 (Transverse Mercator) zone: 0 datum: rome40 ellipsoid: international north: 4834913.78536118 south: 4821192.33020234 west: 1668983.89916058 east: 1685595.88381969 nsres: 100.156607 ewres: 100.07219674 rows: 137 cols: 166 cells: 22742 === Where am I wrong? All the best. -- Paolo Cavallini: http://www.faunalia.it/pc ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: mswindows/osgeo4w/mklibs.sh: dumpbin: command not found
Hi Helmut, Thanks for the reply. Apparently evething is working fine, thanks. Have happy holidays, milton 2009/12/30 Helmut Kudrnovsky hel...@web.de Hi Milton, Dear all, I am still trying to compile the grass under Vista/MSYS, and following the instructions available at: When I run ./mswindows/osgeo4w/package.sh I get the following error message: [...] Mon Dec 28 21:56:25 EST 2009: STARTING cleanup Mon Dec 28 21:56:26 EST 2009: STARTING building vc libraries libgrass_I.6.4.0svn.dll = grass_I mswindows/osgeo4w/mklibs.sh: dumpbin: command not found please comment out from line 94 ( http://trac.osgeo.org/grass/browser/grass/branches/releasebranch_6_4/mswindows/osgeo4w/package.sh#L94 ) to line 136 ( http://trac.osgeo.org/grass/browser/grass/branches/releasebranch_6_4/mswindows/osgeo4w/package.sh#L136 ) in the package.sh. these lines are there for building under MS VisualC and not for building under MSYS. the wiki (http://trac.osgeo.org/grass/wiki/CompileOnWindows) was updated a while ago according to a commit to the source (http://lists.osgeo.org/pipermail/grass-commit/2009-December/011019.html) for simplifying the package.sh-file, but commenting out the lines 94-136 was not in the commit. so this should be done in the source. best regards Helmut __ GRATIS für alle WEB.DE-Nutzer: Die maxdome Movie-FLAT! Jetzt freischalten unter http://movieflat.web.de ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] extract stream heads
Hey Markus, You are the same guy that is working on the Sandro's (UFZ) project and that wrote the r.stream grass program!? That small world. Good to know. cheers miltinho 2009/12/29 Markus Metz markus.metz.gisw...@googlemail.com Milton Cezar Ribeiro wrote: Dear all, I am (finally :-) extracting the drainage from a DEM using r.stream.extract. How can I extract also the stream heads? In the vector output of r.stream.extract, stream heads are all points with stream_type = 'start' in layer 1 or cat = 0 in layer 2. See also Options, stream_vect in the manual ;-) Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] drainage extraction - so homogeneous output across different relief conditions
Dear all, I extracted drainage from a large number of tiles of SRTM DEMs, for several regions of Brazil. The regions have very different relief conditions, ranging from almost flat to very inclined. I have some cartographic maps at scale 1:50,000, and the (reference) densities ranges from 7 to 16%. But I was surprise that using the command below the estimated drainage is (on my point of view) very homogeneous, ranging from 3 to 4% only: r.watershed.exe -m --o memory=2048 elevation=srtm_map threshold=250 stream=srtm_map_drainage Why my output map is so homogeneous accross regions with very different terrain conditions? Any hint are welcome. cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Compiling addons
Dear all, I have compiled grass using the instructions available at http://trac.osgeo.org/grass/wiki/CompileOnWindows But now I need to compile some Addons, and I must confess that I have no idea of how to do that. I know that I need to download the addons sources, but I can't figure out what steps I need to do to get it compiled and installed on the proper direction. I am using OSGeo4W/MSYS under Vista. I need to compile the full grass again or I can just compile the add-ons? Case the answer is to compile only the add-ons, what steps/command I need to run under Msys to get it compiled/installed. Thanks in advance, milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Pre-built binaries (bc and readline): where find
Dear colleagues, I am just updating my grass system, to compile it under Vista. I am following the the instructions available at http://trac.osgeo.org/grass/wiki/CompileOnWindows. But when trying to download the Pre-built binaries bc and readline, I can't find/download it from sourceforge. Any idea? bests milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Pre-built binaries (bc and readline): where find
Dear all, It was solved. I don't know why my internet explorer (under vista) not download the files, but Firefox do. cheers milton 2009/12/28 Milton Cezar Ribeiro miltinho.astrona...@gmail.com Dear colleagues, I am just updating my grass system, to compile it under Vista. I am following the the instructions available at http://trac.osgeo.org/grass/wiki/CompileOnWindows. But when trying to download the Pre-built binaries bc and readline, I can't find/download it from sourceforge. Any idea? bests milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] extract stream heads
Dear all, I am (finally :-) extracting the drainage from a DEM using r.stream.extract. How can I extract also the stream heads? bests milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] fast updating of a large map from several small ones.
Dear all, I have a large raster project (~60,000 x 48,000 pixels) and I extracted the drainage from DEM for ~1,000 subbasin. Now I would like to join all the pieces of drainage, but whem I tryed : g.region rast=`g.mlist pat=ANA_COD5_drenbin_* sep=,` r.series input=`g.mlist pat=ANA_COD5_drenbin_* sep=,` output=ANA_COD5_drenbin_join method=maximum --o I get error because the (I suppose) the shell not support my large number of raster names on rast= and input= at g.region and r.series, respectivelly. As my base map is very large, if I try something like r.mapcalc ANA_COD5_drenbin_joinAUX=if(AnyMAP0, AnyMAP, ANA_COD5_drenbin_join) g.rename rast=ANA_COD5_drenbin_joinAUX,ANA_COD5_drenbin_join it is very time expensive, because of the size of my base/large map. Is there a way of I update only part of a map, using g.region or r.region for each subbasing, update my base map without loss other information already stored on base map? A little bit late, but still wishing Merry Xtmas. miltinho ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Landscape connectivity
Hi Markus, So, if I understood well the practical theory you need is not that one described by Urban Keitt (2001), because they use nodes and edges, and you want a patch-based graph theory, isn't it? May be Fall et all (2007) is a better technical approach? Andrew Fall, Marie-Josee Fortin, Micheline Manseau, and DanO’Brien Spatial Graphs: Principles and Applications for Habitat Connectivity Ecosystems (2007) 10: 448–461 - DOI: 10.1007/s10021-007-9038-7 bests miltinho brazil=toronto 2009/12/2 Markus Metz markus.metz.gisw...@googlemail.com Hi all, for a raster map with habitat patches separated by gaps in between, I want to know for each habitat patch 1) how far away is the nearest other habitat patch, 2) what's the id of the nearest habitat patch (each contiguous patch must have a unique ID), 3) how do I get to the nearest habitat patch, i.e. what is the shortest path connecting one patch and its nearest other patch. Then I need to get clusters of patches, each cluster consisting of patches connected by paths not longer than max_distance, i.e. I can visit all patches within a cluster by traversing gaps not broader than max_distance (hop from one patch to the next, but I can hop only so far...) That would be interesting for meta-populations, dispersal etc, habitats could be islands, forest fragments, trees of a particular species, any habitat type with a patchy spatial distribution and one or more species dependent on that habitat. Any idea how to do that in grass or some other (fragment analysis-specific) application? Using the habitat map twice as input to r.distance would be a start with r.distance maps=habitat_map,habitat_map, but r.distance is too slow for my taste (running for 2 hours by now), and I still have to get the nearest patches, shortest paths and clusters connected with max_distance. Running r.cost or r.walk for each patch as starting points and all other patches as stop points is not an option because I have ~200,000 patches in a 300 million cell raster. Converting the raster habitat map to a vector with areas and using v.net.* modules doesn't work because AFAIK these modules work with points as graph nodes, and some of the areas are quite large, using centroids would give wrong distances. The theory of what I want to do is described in Urban and Keitt (2001) [1] Thanks, Markus M [1] link to pdf: http://www.keittlab.org/~tkeitt/papers/urban-keitt-2001.pdf PS: The grass module I would like to have for this kind of analysis would use a habitat map as input, have as optional output raster paths with distances as cell values, vector paths with costs and ids of the two patches connected by each path, one cluster map for each maxdistance=max1,max2,max3,..., and an ascii file in the same format like r.distance for all paths connecting habitat patches. X-mas is getting close ;-) ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] % of cover class upstream
Hi guys, sorry for my delay. I am trying your suggestions. cheers milton 2009/11/14 Rich Shepard rshep...@appl-ecosys.com On Sat, 14 Nov 2009, Milton Cezar Ribeiro wrote: I have a DEM, a drainage and a Cover mapa, and I need to estimate - for ANY pixel of my drainage - what is the amount of my cover classes that contribute (or are upstream) of that pixel. milton, Not all my data are imported and defined to the critical area so I cannot offer direct experience. However, Single Flow Direction (SFD, D8) moves flow into a single downslope cell and is used by r.watershed and r.terraflow. D-infinity, or vector grid approach, is used by r.flow. Multiple Flow Directions (MFD, what you probably want to use) partitions the flow into two or more downslope directions and is used by r.terraflow, r.topmodel, and r.mapcalc flow implementation in the r.mapcalc tutorial. Also check out r.hydro.CASC2D. HTH, Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] % of cover class upstream
Dear all, I have a DEM, a drainage and a Cover mapa, and I need to estimate - for ANY pixel of my drainage - what is the amount of my cover classes that contribute (or are upstream) of that pixel. Any help are VERY welcome. Cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Data Import and Coordinate Information
Hi Rich, Are you using v.in.ogr with -o flag? bests milton 2009/11/14 Rich Shepard rshep...@appl-ecosys.com For my current project, almost all data come from the Oregon GIS office (now called the Geographic Enterprise Office). All the data available for download (.shp and .e00 primarily) are supposed to be in the state's own coordinate system, Oregon Conformal Conic. This is actually Lambert Conformal Conic with specific bounds, parallels, etc. Since they're all supposed to have the same bounds and other parameters I should be able to import them to the same location and map set. But I cannot do this. If I tell the import wizard to use the existing coordinates it tells me that they don't match what are in the file. So I have to import each theme into a different location. My first question is what might be causing this to happen? I'll reproject them all to a common location/map set when I get them all in, but wonder why there's an issue. More importantly, when I look at the location/PERMANENT/DEFAULT_WIND and WIND files I find the coordinates are far off what they should be, and what I defined them to be in the GUI. For example, instead of North and South bounds being 46.2N and 41.9N the files show them as 63N and 39N. Second question: what might be happening here? In the meantime, I'll create new locations for all themes then read up on how to reproject them so they're all recognized as using the same coordinates and parameters and will overlay properly. As an aside, I don't see a menu item (under File, I suppose it should be) to create a new location and map set for a new import. As a result, I exit grass and restart it. That's too Microsoftish to be correct. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Programación en GIS
Hi Flavio, how about Python? http://grass.osgeo.org/wiki/GRASS_and_Python saludos milton 2009/11/4 Flavio ESPECHE NIEVA flavio.espe...@gmail.com De mi mayor consideración: Soy nuevo en la lista y estoy comenzando a trabajar con GRASS. Necesito desarrollar aplicaciones (en C, Java o PHP) que utilicen GRASS. ¿pueden recomendarme un buen manual o tutorial sobre la API y librerías de este sistema, para comenzar? No importa si está en inglés o en castellano (español). Desde ya, muchas gracias!! -- Flavio ESPECHE NIEVA ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] internal buffer of polygon
Hi Maning, Sometimes I have the same issue, but I don't solve this on vector. Just supose you have a map where 1 is your raster-based-polygon and NULL is outside of polygon. You can use r.mapcalc to invert it, in the maner that 1 come to be your non-polygon, and NULL your polygon. So you can run r.buffer with positive value and at the end may be you need to do some kind of crossing between original and output maps. Bests. miltinho 2009/10/29 maning sambale emmanuel.samb...@gmail.com I tried using a negative values as buffer=-3000 in v.buffer but it doesn't give me internal buffer in a polygon layer. Any hints? -- cheers, maning -- Freedom is still the most radical idea of all -N.Branden wiki: http://esambale.wikispaces.com/ blog: http://epsg4253.wordpress.com/ -- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] stream extraction from ASTER dem
Hi Hamish, I will try without the memory saving flag... Computer with more memory? The mine is 6Gb (!), and I tryed with 1/4 of the image (5000x3000) but still get error. Any other hint are welcome. cheers milton 2009/9/28 Hamish hamis...@yahoo.com try a smaller region (or try on a computer with more memory) (yes, I saw the seg part, hey try without that memory saving flag too!) H --- On Mon, 28/9/09, Markus Metz markus.metz.gisw...@googlemail.com wrote: From: Markus Metz markus.metz.gisw...@googlemail.com Subject: Re: [GRASS-user] stream extraction from ASTER dem To: Milton Cezar Ribeiro miltinho.astrona...@gmail.com Cc: grass-user grass-user grass-user@lists.osgeo.org Received: Monday, 28 September, 2009, 7:14 PM Milton Cezar Ribeiro wrote: Dear all, I have a DEM with dimensions 10,000 x 6,000. Now I need to extract all the stream or drainage, as well as generated some other stream/flow related features like flow direction, stream lenght, etc. I work on a computer with 6Gb under Vista, and the grass version is the OSGeo4W one. But when I try to run the command below I get one error relate to r.watershed.seg.exe: r.watershed elevation=aster_dem stream=aster_dem_stream threshold=150 -m memory=2048 Any idea of how to solve that? Maybe it is the Vista UAC problem, there were recently other posts related to this. Just a guess, Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] stream extraction from ASTER dem
Hi Markus, May be you are right (Vista issue!). Two month ago I was trying to run the same steps on grass, but I was expecting that new OSGeo4W release included some fix for that. :-( Anyway, thanks for your time. Any colleague that could suggest a way of I extract some geomorphometry features (sensu Tom Hengl!) from large DEM (10,000x6,000), please let me know. bests milton 2009/9/28 Markus Metz markus.metz.gisw...@googlemail.com Milton Cezar Ribeiro wrote: Dear all, I have a DEM with dimensions 10,000 x 6,000. Now I need to extract all the stream or drainage, as well as generated some other stream/flow related features like flow direction, stream lenght, etc. I work on a computer with 6Gb under Vista, and the grass version is the OSGeo4W one. But when I try to run the command below I get one error relate to r.watershed.seg.exe: r.watershed elevation=aster_dem stream=aster_dem_stream threshold=150 -m memory=2048 Any idea of how to solve that? Maybe it is the Vista UAC problem, there were recently other posts related to this. Just a guess, Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] stream extraction from ASTER dem
Dear all, I have a DEM with dimensions 10,000 x 6,000. Now I need to extract all the stream or drainage, as well as generated some other stream/flow related features like flow direction, stream lenght, etc. I work on a computer with 6Gb under Vista, and the grass version is the OSGeo4W one. But when I try to run the command below I get one error relate to r.watershed.seg.exe: r.watershed elevation=aster_dem stream=aster_dem_stream threshold=150 -m memory=2048 Any idea of how to solve that? cheers milton brazil=toronto ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] set color for negative and positive values
Dear all, I have a set of 10,000 distance maps. On this map I have negative (-300 to 0) and posivie (0 to +300) values. I would like to attrib a color table to those maps from yellow-cyan-dark green when the negative values goes from 0 to -300, and from yellow-orange-red when values goes from 0 to +300 Any help are welcome. Cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] MAP Calculator mean of eight neighbors to interpolate 0 value in raster map
May be something like: (please, check several values!) r.mapcalc B=if(A==0, sum(A[-1,-1], A[0, -1], A[1,-1], A[-1,0], A[1,0], A[-1,1], A[0, 1], A[1,1])/8, A) cheers milton 2009/8/28 stephen sefick ssef...@gmail.com if(0, (sum([-1,1], [0, 1], [1,1], [-1,0], [0,0], [1,0],[-1,-1], [0,-1],[1,-1])/8)) I want to use mapcalc to do the above. I am missing a syntax step. if a value is zero take the average of the eight nearest neighbors and replace the value with this. Otherwise, keep the value unchanged -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] MAP Calculator mean of eight neighbors to interpolate 0 value in raster map
Hi stephen, just to remember not include yourmap[0,0] bests milton 2009/8/28 stephen sefick ssef...@gmail.com I think so- let me try it On Fri, Aug 28, 2009 at 12:04 PM, achima...@jupiter.uni-freiburg.de wrote: r.mapcalc newmap = if ( yourmap==0 , ( sum( yourmap[-1,1] , yourmap[0,1] , ...)/8 ) , yourmap) is that what you mean? stephen sefick schrieb: if(0, (sum([-1,1], [0, 1], [1,1], [-1,0], [0,0], [1,0],[-1,-1], [0,-1],[1,-1])/8)) I want to use mapcalc to do the above. I am missing a syntax step. if a value is zero take the average of the eight nearest neighbors and replace the value with this. Otherwise, keep the value unchanged -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] MAP Calculator mean of eight neighbors to interpolate 0 value in raster map
Hi Stephen, may be you can share some picture? bests milton 2009/8/28 stephen sefick ssef...@gmail.com This does not quite work for my purposes. The problem is the the 0 values are in a 2 or 3 grid wide seam in between patched raster maps, so the value is lower than it should be using r.mapcalc benning_dem=if(patch_benning==0, (patch_benning[-1,-1]+patch_benning[1,-1]+patch_benning[-1,0]+ patch_benning[1,0]+patch_benning[-1,1]+patch_benning[1,1])/6, patch_benning) any thoughts, Stephen Sefick On Fri, Aug 28, 2009 at 12:08 PM, Milton Cezar Ribeiromiltinho.astrona...@gmail.com wrote: Hi stephen, just to remember not include yourmap[0,0] bests milton 2009/8/28 stephen sefick ssef...@gmail.com I think so- let me try it On Fri, Aug 28, 2009 at 12:04 PM, achima...@jupiter.uni-freiburg.de wrote: r.mapcalc newmap = if ( yourmap==0 , ( sum( yourmap[-1,1] , yourmap[0,1] , ...)/8 ) , yourmap) is that what you mean? stephen sefick schrieb: if(0, (sum([-1,1], [0, 1], [1,1], [-1,0], [0,0], [1,0],[-1,-1], [0,-1],[1,-1])/8)) I want to use mapcalc to do the above. I am missing a syntax step. if a value is zero take the average of the eight nearest neighbors and replace the value with this. Otherwise, keep the value unchanged -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Removing 0 values from raster map
Hi stephen, May be this? g.region rast=`g.mlist pat=srtm* sep=,` r.series input=`g.mlist pat=srtm* sep=,` output=srtm_leste method=maximum of you can play with r.mapcalc to identify where is zero, and calculate the mean using [-1,-1][-1,0], [1,1] (not use [0,0]) to average (sum/8). good luck, miltinho 2009/8/21 Markus Neteler nete...@osgeo.org On Sat, Aug 22, 2009 at 12:33 AM, stephen sefickssef...@gmail.com wrote: I have four maps that I would like to patch together. The problem is that on the boundary there are zero values after patching. I would like to get rid of these values- interpolate the values based on the surrounding values. I can not firgure out how to do this. Any help is greatly appreciated You could try this: http://grass.osgeo.org/grass64/manuals/html64_user/r.fillnulls.html Or r.neighbors with mode and a large moving window? Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Reading list of vector maps into file or array
It is amazing!!! Thanks hamish! milton 2009/8/20 Hamish hamis...@yahoo.com Mark: If I do 'g.list vect walker.list', the list is not in one column and seems to be separated by tabs, with 4 columns and a bunch of rows of vector files. g.mlist I'm building an animated GIF of the walkers over time, on a shaded DEM. The method I am using is tedious and am seeking to automate the process. I load all the layers (shaded DEM, walkers at various intervals) and turn on one walker vector, and output a .png from the display monitor. Wash, rinse, repeat until all walkers have been captured as an image, then animate them with the 'convert' command. see loop shell scripts in http://grass.osgeo.org/wiki/Movies Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Reading list of vector maps into file or array
Hi Mark, Have you tried something like grass.mlist_grouped(type = 'rast') (may be you need change to vect) on python? good luck milton 2009/8/20 M S msei...@gmail.com Greetings All: How would one read the vector maps in a current mapset into a file, or a list in an array (using python)? If I do 'g.list vect walker.list', the list is not in one column and seems to be separated by tabs, with 4 columns and a bunch of rows of vector files. Would a simple solution be to begin processing that redirected output? As background... After running the r.sim.water module, I get a bunch of walker vector maps as outputs for given intervals. Works well thus far. I'm building an animated GIF of the walkers over time, on a shaded DEM. The method I am using is tedious and am seeking to automate the process. I load all the layers (shaded DEM, walkers at various intervals) and turn on one walker vector, and output a .png from the display monitor. Wash, rinse, repeat until all walkers have been captured as an image, then animate them with the 'convert' command. If anyone has some suggestions on how to get vector maps into a list, please let me know. Also, if there is a better way to solve this problem, any feedback would be greatly appreciated. Mark ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] magic?
Hi Achim, could you, please, share with us your region -p ? bests milton 2009/8/20 achim a...@jupiter.uni-freiburg.de Hi, I tried a lot. And I found a solution: I made a buffer around the area with 0 distance. Then v.to.rast works. If anyone knows why, it would be very interesting for me. Moreover, when this happens in a big dataset, the error will probably not recognized. Thanks for attention, Achim achim schrieb: Hi, using a new location: no success export and import as shape-file: no success reproduce it with spearfish dataset: no success Has anybody an idea? I attach the problem-polygon as a shape... hope somebody has the same problem with it ...and an idea. Thanks, Achim achim schrieb: Hello, I used v.to.rast v.to.rast --overwrite input=inmap output=outmap use=cat to convert 29 areas to raster. It says Converted areas: 29 of 29 but one area is missing!!! Its the one with cat number 16. Doing the same just for the selected 16, conversion results in am empty raster map. the same with use=attr or usr=val. Another vector map also does not convert some areas. The missing areas are not special What happens? Magic? I tried it with grass64 and newly compiled grass65. (OpenSuse 11.1 64bit) Achim PS: Layer: 1 Typ AnzahlMinMax Punkt 0 0 0 Linie 0 0 0 Grenze 0 0 0 Zentroid 29 1 29 Fläche 0 0 0 Alle 29 1 29 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] sqlite3 and grass64 OSGeoW4 installation
Dear Grass-gurus, I am running grass64 available on OSGeoW4 Express installation. Now I need to play with grass python and sqlite3. I gone on OSGeoW4 installer again, and when I try to choose the sqlite3 lib to install the installer (available at http://trac.osgeo.org/grass/wiki/CompileOnWindows) day that sqlite3 lib are installed. How can I run a simple test to check if sqlite3 are really installed on my system? When I start Msys/bash and type sqlite3 and press TAB the msys shell suggest me sqlite3.dll. Any help are welcome. cheers milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] bus error in v.in.ogr
Hi Maning, O flyed on your email and saw only warnings, not errors. It appears to worked. bests milton 2009/8/14 maning sambale emmanuel.samb...@gmail.com I am importing a polygon shapefile and I get this error. The shapefile was created using QGIS. v.in.ogr layer=forest_p116r048_repair output=forest_p116r048 dsn=/Users/maning/pfua/map_data/forest_2002_digit/forest_p116r048_p117r047_edited --o WARNING: Vector map forest_p116r048 already exists and will be overwritten Projection of input dataset and current location appear to match WARNING: Vector map forest_p116r048 already exists and will be overwritten Layer: forest_p116r048_repair Importing map 181 features... WARNING: Area size [7.5e-09], area not imported - Building topology for vector map forest_p116r048... 58952 vertices registered 230 primitives registered Building areas: 100% 225 areas built 225 isles built Attaching islands: 100% Attaching centroids: 100% Topology was built Number of nodes : 228 Number of primitives: 230 Number of points: 0 Number of lines : 0 Number of boundaries: 230 Number of centroids : 0 Number of areas : 225 Number of isles : 225 Number of incorrect boundaries : 5 Number of areas without centroid : 225 - WARNING: Cleaning polygons, result is not guaranteed! Building topology for vector map forest_p116r048... Topology was built Number of nodes : 228 Number of primitives: 230 Number of points: 0 Number of lines : 0 Number of boundaries: 230 Number of centroids : 0 Number of areas : - Number of isles : - - Break polygons: Registering points 46354 All points (vertices): 58943 Registered points (unique coordinates): 46354 Points marked for break: 345 Breaks: 610 - Remove duplicates: Duplicates: 230 - Break boundaries: Intersections:39 - Remove duplicates: Duplicates: 7 - Clean boundaries at nodes: Modifications: 0Bus error -- cheers, maning -- Freedom is still the most radical idea of all -N.Branden wiki: http://esambale.wikispaces.com/ blog: http://epsg4253.wordpress.com/ -- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] passing raster map to Python object
Hi Hamish, Thanks for your encourager's email. I started to play with numpy for other purposes, but I will give a look in deep for this. If I get success I will be happy report it here. bests milton 2009/8/12 Hamish hamis...@yahoo.com Milton Cezar Ribeiro wrote: I am building a simulation python code that read a raster maps on grass mapset. As my maps is not so large (512x512 pixels) I dump it on a text file (using r.out.ascii), after I read the text file on python, parse it using a function that I wrote, and after create a matrix-like object on Python. I would like to know if I can read my raster map on python wihout save it on external text file. Case yes, a short example are welcome. you might look at NumPy http://numpy.scipy.org/ + NumPtr python extension. NumPtr can be found in the swig/python/ dir in the grass source code. (also examples/m.distance + rasteraccess.py there) NumPy is also used by the wxGUI profile tool. SciPy is another thing to check out. please post back to the list with what you learn :) if you succeed maybe add a page to the grass wiki with steps? Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] python script: grass.list_grouped for other mapset
Dear all, I am building some codes on python that read grass informations. If I start msys - grass - select my location mapset, I am able to both get the list of rasters or read a raster map. To get the list of raster maps I use something like: landscape_list=grass.list_grouped('rast')['PERMANENT'] This works fine. But my application will need to read information for more than one mapset at same time, and the are not stored on PERMANENT. As I am working with landscape simulations, for each landscape (I have about 30,000 landscapes with 512x512 pixels) I generated about 5 new maps, and for convenience (including backup) echo set of maps are stored on different mapsets. I really need to stay conected on one mapset and get the list of rasters as well as the raster maps from other mapsets, and this need to be done on the fly of my script. Any hint are welcome. cheers milton brazil=toronto ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] passing raster map to Python object
Dear all, I am building a simulation python code that read a raster maps on grass mapset. As my maps is not so large (512x512 pixels) I dump it on a text file (using r.out.ascii), after I read the text file on python, parse it using a function that I wrote, and after create a matrix-like object on Python. I would like to know if I can read my raster map on python wihout save it on external text file. Case yes, a short example are welcome. cheers milton brazil=toronto ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] osgeo4w grass instalation : Tkinter error on python script
Dear all, [Windows 64 / XP] Sorry if this is not the best place to ask for this. I installed the Osgeo4W grass, using the installer available at: http://trac.osgeo.org/grass/wiki/CompileOnWindows After startup the installer, I choose Express instalation, and installed grass from it. Grass is running find. Now I am building a script on python that import grass.script as grass import Tkinter Unfortunatelly I get the following error: $python landscape.py Traceback (most recent call last): File landscape.py, line 6, in module import Tkinter File C:\OSGeo4W\apps\Python25\lib\lib-tk\Tkinter.py, line 38, in module import _tkinter # If this fails your Python may not be configured for Tk ImportError: No module named _tkinter In fact, I have loaded a raster map on python variable (a list of lists), and now I would like to plot it on the screen. Other suggestions different from Tkinter are welcome. bests milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] can I access mapset outside of grass, by using python?
Dear friends, I can run some simple scripts on python if I start Msys - grass64 - python myscrypt.py But now I would like to know if I can build a python code outside of a msys/grass session, and get access on a mapset withour stay running grass. thanks a lot, milton ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] convert a raster map on a matrix on python.
Dear all, how can I read all lines of a raster map and assing it on a matrix (list of lists) on python. I am able to get all informations about a raster map, using for example grass.raster_info('MYMAP') but I need to load the grass raster map as a grid on python. Any help are welcome. Cheers milton brazil=toronto ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] grass 6.4 for win - old tcltk GUI
Hi Edmond, Any error message?! XP / VISTA / 32 / 64 ? MSYS or Standalone version? I can't help, but I think this information could be useful for helpers. Cheeers milton 2009/8/6 Edmondo Elisei edmondo.eli...@gmail.com Dear Grass lists, the old GUI doesn't start. It happens with 6.4 RC5 and 6.4 SVN both (for WIN) thanks for your help regards Edmondo -- -BEGIN PGP PUBLIC KEY BLOCK- Version: GnuPG v1.4.0 (MingW32) - GPGshell v3.47 mQGiBEL7K/4RBACPWlwG4/4nkjfVpHRBr/1tPjTWzs5Lgddtc32XzbvwW6Kb0RzH eWWdTxXLJ+otNT6BaJGQNNGj5UaBgoMlRcjAewHlvslWJfc81xc0kE8jr5QZ+5Zz 3I2OSYddhIvHw5oY86kVZJ1uIB1baRkQbTYo1hiPaj92CU3WKbOA4ckyfwCgz91L EEbEpxlNJYnarj5At5MVihcD/iTwavn90R/qf/AOzizYAcCqNFzcVqOS1K14PHYD O3Y8N9XXPbpCXmQiMPYFfwJDfl6wjOLcdrnGqMH7P61uVCShm+vhzD1pDPXBom2W 8VfcoYtHctbgG9+H2Y9kPTkeIsNmfMAFi7OJeWfaF8P/IXeWS2hof2O61EXZGk1J S2g0A/4oTupCd4IW7wuk7kG/30kZBqGy8LQRK5UkR/gXYbhvXyzrrfvPVoVJ56tp T2gH+sgmI7W42cAMJLMcHvP0uaQ4iFqLlzVraan5Ap7r2Iu4K9mxKy1Xdzv1zcWY q+OrweOwv9AnqykpDzloP6CRWlte7SinuzS3nFQhnAK4i+nSwLQpZWRtb25kbyBl bGlzZWkgPGVkbW9uZG8uZWxpc2VpQGdtYWlsLmNvbT6IXgQTEQIAHgUCQvsr/gIb AwYLCQgHAwIDFQIDAxYCAQIeAQIXgAAKCRC0eTcIfCWG2jwXAJ4mkNv1yh6xo/kA VbntEZZ1bVVtQgCgptslr3dCWykYivo04QlknGyEnLu5AQ0EQvssARAEALCTTEbp +9kMtkDh/LoOlow+lNDZ3w8DdgikKrst1DC/TaoL01LoUyCVa8FkJ4eS1G8i77Hz hsyfGEtTDB4uv65+HRamLh9RlSdFw3gZ5LOU7Syn7vC7zZ4peW9gpiuNi0nh8d1z SaMK+b1esnUB1DdbD174UpqdPgaZQRP8/i43AAMFA/47MFVyh2V6/KF3nOzFabgb 402admgZ3/z/rsxlwOGYERk8PYHgw7BXx09xu8I+By7NM+NZNajjskcuPRxbWZeJ u3Q4aOK9/0LSKQDJ71NJm5BC9nRFxK/vJgLtOm8Gwdx/SZ0CgjKK5U7g7ArMS08w MDr3hydn62YWZotCc5nQ04hJBBgRAgAJBQJC+ywBAhsMAAoJELR5Nwh8JYba1kUA n2gLdxJMSvXu2Evs1g6d5It9Bx3kAJ90oIbRVsh+6cJtZUzKhbR6dZDAfw== =F04X -END PGP PUBLIC KEY BLOCK- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] g.list not working on GUI after edit .grassrc6
Hi Nikos Micha, Nikos, I tryed: (it is little different from my post because there I give a example) $ /c/OSGeo4W/apps/grass/bin/grass64 -gui C:/Data/Tibagi/vegPR_probio_proatlan_maio2009/tibagi_Probio_grass_sdb/ne wLocation/PERMANENT Cleaning up temporary files ... Starting GRASS ... c:/Users/famiglia/C:/Data/Tibagi/vegPR_probio_proatlan_maio2009/tibagi_Probio_grass_sdb/newLocation/PERMANENT: Not a valid GRASS location see that besides I grass try start my local from c:/Data/... grass add c:/Users/famiglia first. Famiglia is the name of my computer. May be we have a issue to solve here. Anyway Following the suggestion by Micha (remove .grassrc6) everything is ok now!!! Thanks a lot Nikos Micha! milton 2009/8/5 Micha Silver mi...@arava.co.il Nikos Alexandris wrote: On Tue, 2009-08-04 at 19:33 -0400, Milton Cezar Ribeiro wrote: Hi Nikos. Thanks for the reply. Yes, location is included on .grassrc6. $ cat .grassrc6 GISDBASE: C:/Mymapset LOCATION_NAME: newLocation MAPSET: PERMANENT GRASS_GUI: tcltk Any other idea?! cheers milton Milton, from your previous post in the other thread it seems that you do not define the location/mapset [1]. Did you try to start like: $ /c/OSGeo4W/apps/grass/bin/grass64 -gui C:/Mymapset/newLocation/PERMANENT or $ /c/OSGeo4W/apps/grass/bin/grass64 -text C:/Mymapset/newLocation/PERMANENT Also, I think you can just blow away the .grass6rc file and start fresh. This way, on starting, GRASS should ask where the GISDBASE is, and offer to create locations and mapsets... ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] can I run python code in standalone grass 6.4?
Dear all, Is this possible? I installed standalone grass 6.4. from http://grass.itc.it/grass64/binary/mswindows/native/WinGRASS-6.4.0RC5-1-Setup.exe After start this mswindows version I can start MSYS/grass64, and also run python test.py. Now I would like to run some grass-python code but when I try import grass.script as grass or import core from grass or from grass.script import core as grass I get error messages like: ImportError: No module named grass.script any help are welcome. miltinho ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user