Re: [GRASS-user] Error in GUI's sub-group definition in i.maxlik and others modules
Hello Martin I got this error in GRASS6.4.0.RC6, from two weeks ago in LinuxUbuntu9.10 but also in Windows. Jenny On Fri, Jun 25, 2010 at 9:57 PM, Martin Landa landa.mar...@gmail.comwrote: Hi, 2010/6/24 Jenny Turner jennyturner1...@gmail.com: I have been testing i.maxlik and I'm getting a strange error. When I select the group, in group parameter, automatically, sub-group list includes the available sub-groups but, when i select one, it's selected in the frame but it's not inserted in the command to be runned. In the bottom it's still REQUIRED. I can only run i.maxlik by inserting directly the parameters in the commandline-output. Has anyone having the same problem? GRASS version? OS? Martin -- Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Customization of wxGUI
Hi there I want to change a bit the look of my GRASS's wxGUI (just testing colours, patterns, icons and stuff like that). In which file is all of this defined? Or is it any wiki explaining where/how to change it? Thanks Pedro ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Changing menudata.xml to include an direct processing
Hi Martin, it's not possible in GRASS6.4.0RC6. It still opens the module-panel and parameterizes it. has anyone ever tried to have a module that runs, from the menu, without opening module-panel? Thanks Best regards, Franz On Fri, Jun 18, 2010 at 1:35 PM, Martin Landa landa.mar...@gmail.comwrote: Hi, 2010/6/9 Franz Schiller franzschiller1...@gmail.com: Greetings I'm making some addaptions in menudata.xml in order to have a few functions available more rapidly (for me). I wanted to have a menu-item (I know how to add it) that runs directly without opening any panel. This would mean that all parameterization would be in the xml. Is this possible? yes, it's possible in 6.5+. In 6.4 it will be probably not working. menuitem labelInfo elevation map/label helpShow menudata of elevation raster map/help handlerOnMenuCmd/handler commandr.info map=elevation/command /menuitem Martin -- Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] shell script
Hi I'm trying to run the script below but it's giving me an error saying default region is not set. Which command should I use to set the default region and do I put it? #!/bin/sh #variable to customize: # path to GRASS software main directory GISBASE=/usr/lib/grass64 # path to GRASS database GISDBASE=$HOME/grassdem # nothing to change below MAP=$1 LOCATION=$2 # generate temporal LOCATION: TEMPDIR=FLOODS mkdir -p $GISDBASE/$TEMPDIR/PERMANENT # save existing $HOME/.grassrc6 if test -e $HOME/.grassrc6 ; then mv $HOME/.grassrc6 /tmp/$TEMPDIR.grassrc6 fi echo LOCATION_NAME: $TEMPDIR $HOME/.grassrc6 echo MAPSET:PERMANENT $HOME/.grassrc6 echo DIGITIZER: none $HOME/.grassrc6 echo GISDBASE: $GISDBASE $HOME/.grassrc6 export GISBASE=$GISBASE # Create a WIND file with minimal information and no projection: echo proj: 0 zone: 0 north: 1 south: 0 east: 1 west: 0 cols: 1 rows: 1 e-w resol: 1 n-s resol: 1 top:1 bottom: 0 cols3: 1 rows3: 1 depths: 1 e-w resol3: 1 n-s resol3: 1 t-b resol: 1 $TEMPDIR/$LOCATION_NAME/$MAPSET/WIND # Copy WIND-file to DEFAULT_WIND: cp $TEMPDIR/$LOCATION_NAME/$MAPSET/WIND \ $TEMPDIR/$LOCATION_NAME/$MAPSET/DEFAULT_WIND # Set default database driver: echo DBF_DRIVER: dbf DB_DATABASE : $GISDBASE/$LOCATION_NAME/$MAPSET/dbf/ $TEMPDIR/$LOCATION_NAME/$MAPSET/VAR export PATH=$GISBASE/bin:$GISBASE/scripts:$PATH export LD_LIBRARY_PATH=$GISBASE/lib:$LD_LIBRARY_PATH export GIS_LOCK=$$ export GISRC=$HOME/.grassrc6 db.connect driver=dbf database='$GISDBASE/$LOCATION_NAME/$MAPSET/dbf/' # this should print GRASS version used: g.version g.proj --quiet -c epsg=4326 g.region --quiet -s n=90 s=-90 w=-180 e=180 res=1 # other calculations go here # import rainfall data set. # cd /home/tgumede1/grassdem r.in.gdal input=$HOME/grassdem/TRMMLast1day.tif output=rainfall g.region rast=rainfall -- Kind Regards TS Gumede CSIR, Meraka Institute 072 258 1650 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: shell script
Sandile Gumede wrote: Hi I'm trying to run the script below but it's giving me an error saying default region is not set. Which command should I use to set the default region and do I put it? #!/bin/sh #variable to customize: # path to GRASS software main directory GISBASE=/usr/lib/grass64 # path to GRASS database GISDBASE=$HOME/grassdem # nothing to change below MAP=$1 LOCATION=$2 # generate temporal LOCATION: TEMPDIR=FLOODS mkdir -p $GISDBASE/$TEMPDIR/PERMANENT # save existing $HOME/.grassrc6 if test -e $HOME/.grassrc6 ; then mv $HOME/.grassrc6 /tmp/$TEMPDIR.grassrc6 fi echo LOCATION_NAME: $TEMPDIR $HOME/.grassrc6 echo MAPSET:PERMANENT $HOME/.grassrc6 echo DIGITIZER: none $HOME/.grassrc6 echo GISDBASE: $GISDBASE $HOME/.grassrc6 export GISBASE=$GISBASE # Create a WIND file with minimal information and no projection: echo proj: 0 zone: 0 north: 1 south: 0 east: 1 west: 0 cols: 1 rows: 1 e-w resol: 1 n-s resol: 1 top:1 bottom: 0 cols3: 1 rows3: 1 depths: 1 e-w resol3: 1 n-s resol3: 1 t-b resol: 1 $TEMPDIR/$LOCATION_NAME/$MAPSET/WIND # Copy WIND-file to DEFAULT_WIND: cp $TEMPDIR/$LOCATION_NAME/$MAPSET/WIND \ $TEMPDIR/$LOCATION_NAME/$MAPSET/DEFAULT_WIND Where are the variables $LOCATION_NAME and $MAPSET being set? As far as I can see they seem to be empty - could that be the problem? Paul ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] debian lenny : /usr/lib/grass64/etc/prompt.sh: line 2: g.gisenv: command not found
Hi. I installed the package Version: 6.4.0~rc6-1~gfossit50+1 on debian lenny (Linux version 2.6.26-2-amd64 (Debian 2.6.26-21lenny4) (da...@debian.org) (gcc version 4.1.3 20080704 (prerelease) (Debian 4.1.2-25)) #1 SMP Tue Mar 9 22:29: 32 UTC 2010). When I start grass and select my location (which was created by a previous version of grass 6.4. - RC3) the shell and gui start just fine, but the shell has the error message Cleaning up temporary files ... Starting GRASS ... GISDBASE='/data/z/gis/grass/locations'; LOCATION_NAME='longlatWGS84_EPSG4326_all'; MAPSET='PERMANENT'; __ ___ _____ / / __ \/ | / ___/ ___/ / / _/ ___/ / / __/ /_/ / /| | \__ \\_ \ / / __ / / \__ \ / /_/ / _, _/ ___ |___/ /__/ / / /_/ // / ___/ / \/_/ |_/_/ |_/// \/___/// Welcome to GRASS 6.4.0RC6 (2010) GRASS homepage: http://grass.osgeo.org/ This version running thru: Bash Shell (/bin/bash) Help is available with the command: g.manual -i See the licence terms with: g.version -c If required, restart the GUI with: g.gui oldtcltk When ready to quit enter:exit /usr/lib/grass64/etc/prompt.sh: line 2: g.gisenv: command not found /usr/lib/grass64/etc/prompt.sh: line 2: g.gisenv: command not found /usr/lib/grass64/etc/prompt.sh: line 2: g.gisenv: command not found GRASS 6.4.0RC6 (longlatWGS84_EPSG4326_all):~ No commands are recognised in the shell (as stated in the error message). so e.g. g.list type=vect also fails. BUT strangely, issuing the command from the gui works (oldtcltk and tcltk)! My problem is that I need the shell for scripted gis operations. The info current GRASS environment settings using the gui shows g.gisenv DIGITIZER='none'; GISDBASE='/data/z/gis/grass/locations'; GML_OGR_PATH='/usr/bin/ogr2ogr'; MAPSET='PERMANENT'; LOCATION_NAME='longlatWGS84_EPSG4326_all'; DEBUG='0'; MONITOR='x0'; GRASS_GUI='oldtcltk'; I can't find anything about this problem on Google. Please help. Eduard ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Really dumb script question
It's like this. I'm good with Grass-GIS and not good with complicated scripts. An r.in.wms process downloaded too many tiles, and broke r.in.gdalwarp. So, I would like to break up the downloaded tiles into groups of ten, run r.in.gdalwarp on the groups, then patch the grouped tiles together at the end. I have: Background_ _ 0 to 251.png ie - Background__0.png to Background__251.png What I want to do is for i in group of ten tiles; do r.in.gdalwarp input=grouped filenames output=subset method=nearest s_srs=EPSG:4326 done But how do I break the set of 251 tiles into groups for this? Richard ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Convert FCELL to CELL
Hi, How can I convert a raster type FCELL into a raster type CELL? Thanks, -- Eloi Ribeiro GIS Analyst 39,45º -4,40º http://eloiribeiro.wordpress.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Really dumb script question
Hello Richard, just a suggestion, though I'm not a bash expert : k=10 while [ $k -le 250 ]; do i=0 while [ $i -lt $k ]; do image=`echo Background_$i.png` your stack of operations on $image; ((i++)); done k=$((k+10)); done hmmm, perhaps it's not clean at interval borders, but it may be a good beginning ? Good luck, Vincent. Le lundi 28 juin 2010 à 21:17 +1000, Richard Chirgwin a écrit : It's like this. I'm good with Grass-GIS and not good with complicated scripts. An r.in.wms process downloaded too many tiles, and broke r.in.gdalwarp. So, I would like to break up the downloaded tiles into groups of ten, run r.in.gdalwarp on the groups, then patch the grouped tiles together at the end. I have: Background_ _ 0 to 251.png ie - Background__0.png to Background__251.png What I want to do is for i in group of ten tiles; do r.in.gdalwarp input=grouped filenames output=subset method=nearest s_srs=EPSG:4326 done But how do I break the set of 251 tiles into groups for this? Richard ___ 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] Convert FCELL to CELL
On 28/06/2010 14:40, Eloi Ribeiro wrote: Hi, How can I convert a raster type FCELL into a raster type CELL? Wouldn't: r.mapcalc new_map="int(original_map)" do that? Thanks, -- Eloi Ribeiro GIS Analyst 39,45 -4,40 http://eloiribeiro.wordpress.com This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver http://www.surfaces.co.il/ Arava Development Co. +972-52-3665918 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Really dumb script question
Oops, checking for your reply, I notice an error in my code at line 3 : replace i=0 with i=$((k-10)) Yours, Vincent. Le lundi 28 juin 2010 à 14:06 +0200, Vincent Bain a écrit : Hello Richard, just a suggestion, though I'm not a bash expert : k=10 while [ $k -le 250 ]; do i=0 while [ $i -lt $k ]; do image=`echo Background_$i.png` your stack of operations on $image; ((i++)); done k=$((k+10)); done hmmm, perhaps it's not clean at interval borders, but it may be a good beginning ? Good luck, Vincent. Le lundi 28 juin 2010 à 21:17 +1000, Richard Chirgwin a écrit : It's like this. I'm good with Grass-GIS and not good with complicated scripts. An r.in.wms process downloaded too many tiles, and broke r.in.gdalwarp. So, I would like to break up the downloaded tiles into groups of ten, run r.in.gdalwarp on the groups, then patch the grouped tiles together at the end. I have: Background_ _ 0 to 251.png ie - Background__0.png to Background__251.png What I want to do is for i in group of ten tiles; do r.in.gdalwarp input=grouped filenames output=subset method=nearest s_srs=EPSG:4326 done But how do I break the set of 251 tiles into groups for this? Richard ___ 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
Re: [GRASS-user] Convert FCELL to CELL
2010/6/28 Micha Silver mi...@arava.co.il On 28/06/2010 14:40, Eloi Ribeiro wrote: Hi, How can I convert a raster type FCELL into a raster type CELL? Wouldn't: r.mapcalc new_map=int(original_map) do that? Thanks, -- Eloi Ribeiro GIS Analyst 39,45º -4,40º http://eloiribeiro.wordpress.com This mail was received via Mail-SeCure System. ___ grass-user mailing listgrass-u...@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silverhttp://www.surfaces.co.il/ Arava Development Co. +972-52-3665918 That's it! Thanks Micha! -- Eloi Ribeiro GIS Analyst 39,45º -4,40º http://eloiribeiro.wordpress.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] remote sensing (télédétection)
On Sun, Jun 27, 2010 at 12:10 PM, Etienne DELAY (deust 2010) etienne.de...@etu.unilim.fr wrote: Hello I am a Master (french) student, and I work in Valle d'Aosta on vines issues . I'm a noob to 3 days. I discovered the potential of GRASS with great pleasure. I wanted to know if it was possible to make remote sensing (télédétection) with this fabulous software? R. Zorer, L. Delucchi, M. Neteler, G. Nicolini, 2010: Development of a GRASS-GIS application for the characterization of vineyards in the Province of Trento (Poster). VIII International Terroir Congress, Soave (Verona) Italy. June 14–18, 2010, http://terroir2010.entecra.it/index_en.html Here's the poster link: http://gis.fem-environment.eu/uploads/Zorer2010_et_al_Terroir_Poster.pdf Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Retrieve pixel resolution using Python Script
If you use gdalinfo on an image, it will output the Pixel Size. Hope that helps. Mark 2010/6/28 António Rocha antonio.ro...@deimos.com.pt: Greetings I'm trying to retrieve an image (not yet imported to GRASS) pixel resolution through a Python Script but so far i was unable to figure out a method to do this. Is this possible? I have used g.proj to retrieve its projection and datum but not the Pixel resolution. THanks for your help Best regards, Antonio R. __ Information from ESET NOD32 Antivirus, version of virus signature database 5234 (20100628) __ 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 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Retrieve pixel resolution using Python Script
Olá António Rocha, This works for me, give it a try: ### import os try: from osgeo import gdal from osgeo import osr except: import gdal import osr def gdalinfo(pasta): '''Le metadados com GDAL''' ficheiros = os.listdir(pasta) for ficheiro in ficheiros: if ficheiro[-3:].lower() == 'tif': # abrir imagem filename = pasta+ficheiro dataset = gdal.Open(filename, gdal.GA_ReadOnly) if dataset is None: print 'Nao foi possivel abrir a imagem %s!' % ficheiro RasterXSize = dataset.RasterXSize RasterYSize = dataset.RasterYSize print RasterXSize,RasterYSize print '\nfim' return gdalinfo('/home/eloi/Geo/raster/20790/') ### 2010/6/28 Mark Seibel msei...@gmail.com If you use gdalinfo on an image, it will output the Pixel Size. Hope that helps. Mark 2010/6/28 António Rocha antonio.ro...@deimos.com.pt: Greetings I'm trying to retrieve an image (not yet imported to GRASS) pixel resolution through a Python Script but so far i was unable to figure out a method to do this. Is this possible? I have used g.proj to retrieve its projection and datum but not the Pixel resolution. THanks for your help Best regards, Antonio R. __ Information from ESET NOD32 Antivirus, version of virus signature database 5234 (20100628) __ 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 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Eloi Ribeiro GIS Analyst 39,45º -4,40º http://eloiribeiro.wordpress.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Retrieve pixel resolution using Python Script
Hello Mark I know that I can retrieve that information with gdalinfo. But I need to have it on a variable and get/use it in a Python Script. Mark Seibel wrote: If you use gdalinfo on an image, it will output the Pixel Size. Hope that helps. Mark 2010/6/28 António Rocha antonio.ro...@deimos.com.pt: Greetings I'm trying to retrieve an image (not yet imported to GRASS) pixel resolution through a Python Script but so far i was unable to figure out a method to do this. Is this possible? I have used g.proj to retrieve its projection and datum but not the Pixel resolution. THanks for your help Best regards, Antonio R. __ Information from ESET NOD32 Antivirus, version of virus signature database 5234 (20100628) __ 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 __ Information from ESET NOD32 Antivirus, version of virus signature database 5234 (20100628) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5235 (20100628) __ 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
Re: [GRASS-user] debian lenny : /usr/lib/grass64/etc/prompt.sh: line 2: g.gisenv: command not found
Seems to be broken install. Start GRASS and check output of: env | grep GISBASE It should give something like: GISBASE=/path/to/grass-6.4 If GISBASE is incorrect, all GRASS comands will fail. If it's so, check contents of grass64 file. (/usr/bin/grass64) Maris. 2010/6/28, Eduard Ruffert ruff...@mkm.de: Hi. I installed the package Version: 6.4.0~rc6-1~gfossit50+1 on debian lenny (Linux version 2.6.26-2-amd64 (Debian 2.6.26-21lenny4) (da...@debian.org) (gcc version 4.1.3 20080704 (prerelease) (Debian 4.1.2-25)) #1 SMP Tue Mar 9 22:29: 32 UTC 2010). When I start grass and select my location (which was created by a previous version of grass 6.4. - RC3) the shell and gui start just fine, but the shell has the error message Cleaning up temporary files ... Starting GRASS ... GISDBASE='/data/z/gis/grass/locations'; LOCATION_NAME='longlatWGS84_EPSG4326_all'; MAPSET='PERMANENT'; __ ___ _____ / / __ \/ | / ___/ ___/ / / _/ ___/ / / __/ /_/ / /| | \__ \\_ \ / / __ / / \__ \ / /_/ / _, _/ ___ |___/ /__/ / / /_/ // / ___/ / \/_/ |_/_/ |_/// \/___/// Welcome to GRASS 6.4.0RC6 (2010) GRASS homepage: http://grass.osgeo.org/ This version running thru: Bash Shell (/bin/bash) Help is available with the command: g.manual -i See the licence terms with: g.version -c If required, restart the GUI with: g.gui oldtcltk When ready to quit enter:exit /usr/lib/grass64/etc/prompt.sh: line 2: g.gisenv: command not found /usr/lib/grass64/etc/prompt.sh: line 2: g.gisenv: command not found /usr/lib/grass64/etc/prompt.sh: line 2: g.gisenv: command not found GRASS 6.4.0RC6 (longlatWGS84_EPSG4326_all):~ No commands are recognised in the shell (as stated in the error message). so e.g. g.list type=vect also fails. BUT strangely, issuing the command from the gui works (oldtcltk and tcltk)! My problem is that I need the shell for scripted gis operations. The info current GRASS environment settings using the gui shows g.gisenv DIGITIZER='none'; GISDBASE='/data/z/gis/grass/locations'; GML_OGR_PATH='/usr/bin/ogr2ogr'; MAPSET='PERMANENT'; LOCATION_NAME='longlatWGS84_EPSG4326_all'; DEBUG='0'; MONITOR='x0'; GRASS_GUI='oldtcltk'; I can't find anything about this problem on Google. Please help. Eduard ___ 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] Entering point data with SQLite
Hi I am trying to add point data with SQLite without much success. 1. Is it better to simply create a new table in your sqlite.db or should you create a separate database and use v.in.db? 2. If you are using latitude and longitude coordinates, what is the best Field type to use? I was using the degree minute second format dd:mm:ss.W. eg. 71:43:08.2W. I concede that this may not be the best format. There doesn't seem to be a lot of documentation on this part of database management for GRASS. Also, I have been using SQLite Database Browser 1.3. It doesn't place the rows of field names in the order they are created. In fact there seems to be no order at all and no way to rearrange the order. As always, any help will be greatly appreciated. Kurt ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Really dumb script question
Vincent - thanks twice! I'll let you know how it goes. Richard Vincent Bain wrote: Oops, checking for your reply, I notice an error in my code at line 3 : replace i=0 with i=$((k-10)) Yours, Vincent. Le lundi 28 juin 2010 à 14:06 +0200, Vincent Bain a écrit : Hello Richard, just a suggestion, though I'm not a bash expert : k=10 while [ $k -le 250 ]; do i=0 while [ $i -lt $k ]; do image=`echo Background_$i.png` your stack of operations on $image; ((i++)); done k=$((k+10)); done hmmm, perhaps it's not clean at interval borders, but it may be a good beginning ? Good luck, Vincent. Le lundi 28 juin 2010 à 21:17 +1000, Richard Chirgwin a écrit : It's like this. I'm good with Grass-GIS and not good with complicated scripts. An r.in.wms process downloaded too many tiles, and broke r.in.gdalwarp. So, I would like to break up the downloaded tiles into groups of ten, run r.in.gdalwarp on the groups, then patch the grouped tiles together at the end. I have: Background_ _ 0 to 251.png ie - Background__0.png to Background__251.png What I want to do is for i in group of ten tiles; do r.in.gdalwarp input=grouped filenames output=subset method=nearest s_srs=EPSG:4326 done But how do I break the set of 251 tiles into groups for this? Richard ___ 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