Re: [GRASS-user] Long post - Remote Sensing - Classifications
On Wed, 2009-02-25 at 20:46 -0800, Hamish wrote: ok so from the command line 'd.mon x0' should work to open a monitor, and you could try the steps first from the command line (keep it simple). If all works there you can move to trying launching it from a GUI. i.class will always want to run in an xmon + console though. this gets i.class working for me using the sample Spearfish imagery dataset: g.region rast=spot.comp d.mon x0 When I do that, I get the following warning using the wxPython GUI. http://n2.nabble.com/file/n2389455/dx0.png dx0.png d.rast spot.comp # need to do following in this order to get the r,g,b bands in correct # order in the group/$GROUP/REF file (??) I'm not really sure. # create group i.group gr=spotspear in=spot.ms.1,spot.ms.2,spot.ms.3 # create subgroup i.group gr=spotspear subgr=spot.test in=spot.ms.1,spot.ms.2,spot.ms.3 Yes. # add another map for fun i.group gr=spotspear in=spot.comp # run it i.class map=spot.comp group=spotspear1 subgroup=spot.test outsig=test.txt it doesn't seem to mind that I was running it from the user1 mapset and the images were stored in PERMANENT. and then work your way though the menus. the result is a signature file in $MAPSET/group/$GROUPNAME/subgroup/$SUBGRNAME/sig/sig.txt I'm not sure what good results are supposed to look like. That's the point. Isn't supposed to be a monitor display? ..divided into three parts. In the upper right corner you see the image. In the lower right corner, zoomed map portions will be displayed when using the zoom function. In the left section of the monitor, histograms for the selected training areas for all channels will be displayed. The training areas can be digitized (Neteler, M., H. Mitasova 2008, pp. 320-321) ...right? It's not happening to me. because I'm runnig this magic on Lenny ;-) (I actually got the 6.4 version from Sid repository). nice, I hadn't realized it has been uploaded there yet. is there any way in GRASS to link different map displays I am not sure what you mean by that. zoom in one moves both? Yes... and have a brushing (update while dragging) d.rast.edit? No. They are general simple windows displaying info regarding coordinates and Digital Number values like: Diplay: # R:# G:# B:# Proj: UTM ZONE Map: #, # meters Lat Long: DD, DD Moving the cursor on 2 linked displays you have the same geo-info (you are supposed to link 2 linkable displays) but different RGB values (generally 2 -- or more -- different bands or color composition are loaded on different displays). The meaning of this is allowing a visual inspection (in terms of either different color composition or shades of gray) of the same area, supported by a numerical info (DN -- therefore a precise representation of the reflectance behavior of different surfaces). I don't know if I made myself clear..I can show you a couple of screen-shots as soon as I get to the department (it's going to be Envi and Erdas though...I hope you guys don't mind commercial images on this forum :-)) or even a static d.what.rast? I am not sure what you mean by that. maybe r.what does what you want? In order to use r.what you are supposed to use a separate file containing the points you are querying, right? No, this is not exactly what I'm talking about...or better, the outcome is but it's not the interaction that I'm looking for. I'm talking about a pointer based like d.what.rast one. Hamish mitch -- View this message in context: http://n2.nabble.com/Long-post---Remote-Sensing---Classifications-tp2347807p2389455.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
On 26/02/09 13:58, mitch_TX wrote: On Wed, 2009-02-25 at 20:46 -0800, Hamish wrote: ok so from the command line 'd.mon x0' should work to open a monitor, and you could try the steps first from the command line (keep it simple). If all works there you can move to trying launching it from a GUI. i.class will always want to run in an xmon + console though. this gets i.class working for me using the sample Spearfish imagery dataset: g.region rast=spot.comp d.mon x0 When I do that, I get the following warning using the wxPython GUI. http://n2.nabble.com/file/n2389455/dx0.png dx0.png Hamish meant running these modules from the command line. That's the point. Isn't supposed to be a monitor display? ..divided into three parts. In the upper right corner you see the image. In the lower right corner, zoomed map portions will be displayed when using the zoom function. In the left section of the monitor, histograms for the selected training areas for all channels will be displayed. The training areas can be digitized (Neteler, M., H. Mitasova 2008, pp. 320-321) Again, this is when run from the command line using the x-monitors (i.e. d.mon). No. They are general simple windows displaying info regarding coordinates and Digital Number values like: Diplay: # R:# G:# B:# Proj: UTM ZONE Map: #, # meters Lat Long: DD, DD Moving the cursor on 2 linked displays you have the same geo-info (you are supposed to link 2 linkable displays) but different RGB values (generally 2 -- or more -- different bands or color composition are loaded on different displays). The meaning of this is allowing a visual inspection (in terms of either different color composition or shades of gray) of the same area, supported by a numerical info (DN -- therefore a precise representation of the reflectance behavior of different surfaces). I don't know if I made myself clear..I can show you a couple of screen-shots as soon as I get to the department (it's going to be Envi and Erdas though...I hope you guys don't mind commercial images on this forum :-)) In the wxGUI, if you display the two raster in the same Map Display, select both of them in the layer manager and use the query tool in the Display, you get the numerical information on both maps. But you have to open another Map Display if you want to _see_ both maps. or even a static d.what.rast? I am not sure what you mean by that. maybe r.what does what you want? In order to use r.what you are supposed to use a separate file containing the points you are querying, right? No, this is not exactly what I'm talking about...or better, the outcome is but it's not the interaction that I'm looking for. I'm talking about a pointer based like d.what.rast one. You have the query tool in the Map Display of the GUI. Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
On Wed, 2009-02-25 at 21:28 -0800, Hamish wrote: Hamish: i.class map=spot.comp group=spotspear subgroup=spot.test outsig=test.txt hmmm, if I try the same from gis.m in grass65 I get: = [...] required: NO Enter the name of an existing input file Hit RETURN to cancel request [$GISBASE]/etc/grass-run.sh: line 29: 20596 Segmentation fault $@ ERROR: i.class exited abnormally. Press enter to continue. = and the xterm closes. Yes, this is what I was talking about in one of my previous post while using the wxPython GUI I get this http://n2.nabble.com/file/n2390051/x0.png x0.png but we already said that was under development. Ok, it's a bug: the menu item wasn't opening an xmon for you, and wasn't disabling GRASS_RENDER_IMMEDIATE for i.class. For the Tcl/Tk GIS.m GUI now fixed in svn, rev 36103. cd gui/tcltk/ svn up make I'm sorry Hamish, could you be a bit more for dummies for me ? :-) I'm still not so familiar with this. Do I have to go in the tcltk directory and run svn up - make? Or do I have to download something somewhere. Again, I got the 6.4 version from Sid but now I switched back on Lenny...would that work the same? Besides, are you suggesting to switch back to the tcltk GUI to make it work?..I kind of liked the wxPython one. Hamish Cheers, mitch -- View this message in context: http://n2.nabble.com/Long-post---Remote-Sensing---Classifications-tp2347807p2390051.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
On Thu, 2009-02-26 at 14:25 +0100, Moritz Lennert wrote: On 26/02/09 13:58, mitch_TX wrote: On Wed, 2009-02-25 at 20:46 -0800, Hamish wrote: ok so from the command line 'd.mon x0' should work to open a monitor, and you could try the steps first from the command line (keep it simple). If all works there you can move to trying launching it from a GUI. i.class will always want to run in an xmon + console though. this gets i.class working for me using the sample Spearfish imagery dataset: g.region rast=spot.comp d.mon x0 When I do that, I get the following warning using the wxPython GUI. http://n2.nabble.com/file/n2389455/dx0.png dx0.png Hamish meant running these modules from the command line. That's the point. Isn't supposed to be a monitor display? ..divided into three parts. In the upper right corner you see the image. In the lower right corner, zoomed map portions will be displayed when using the zoom function. In the left section of the monitor, histograms for the selected training areas for all channels will be displayed. The training areas can be digitized (Neteler, M., H. Mitasova 2008, pp. 320-321) Again, this is when run from the command line using the x-monitors (i.e. d.mon). The funny thing is that I thought I did all of the above (from CL)...I'll try again. No. They are general simple windows displaying info regarding coordinates and Digital Number values like: Diplay: # R:# G:# B:# Proj: UTM ZONE Map: #, # meters Lat Long: DD, DD Moving the cursor on 2 linked displays you have the same geo-info (you are supposed to link 2 linkable displays) but different RGB values (generally 2 -- or more -- different bands or color composition are loaded on different displays). The meaning of this is allowing a visual inspection (in terms of either different color composition or shades of gray) of the same area, supported by a numerical info (DN -- therefore a precise representation of the reflectance behavior of different surfaces). I don't know if I made myself clear..I can show you a couple of screen-shots as soon as I get to the department (it's going to be Envi and Erdas though...I hope you guys don't mind commercial images on this forum :-)) In the wxGUI, if you display the two raster in the same Map Display, select both of them in the layer manager and use the query tool in the Display, you get the numerical information on both maps. But you have to open another Map Display if you want to _see_ both maps. or even a static d.what.rast? I am not sure what you mean by that. maybe r.what does what you want? In order to use r.what you are supposed to use a separate file containing the points you are querying, right? No, this is not exactly what I'm talking about...or better, the outcome is but it's not the interaction that I'm looking for. I'm talking about a pointer based like d.what.rast one. You have the query tool in the Map Display of the GUI. Moritz I'll try this too..I gotta go now though. cheers, mitch ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
Hamish: hmmm, if I try the same from gis.m in grass65 I get: = [...] required: NO Enter the name of an existing input file Hit RETURN to cancel request [$GISBASE]/etc/grass-run.sh: line 29: 20596 Segmentation fault $@ ERROR: i.class exited abnormally. Press enter to continue. = and the xterm closes. (fixed in svn) Yes, this is what I was talking about in one of my previous post while using the wxPython GUI I get this http://n2.nabble.com/file/n2390051/x0.png x0.png anything in the terminal window? or does that launch an xterm for you? but we already said that was under development. that's what we're doing now. :) Ok, it's a bug: the menu item wasn't opening an xmon for you, and wasn't disabling GRASS_RENDER_IMMEDIATE for i.class. For the Tcl/Tk GIS.m GUI now fixed in svn, rev 36103. cd gui/tcltk/ svn up make I'm sorry Hamish, could you be a bit more for dummies for me ? :-) I'm still not so familiar with this. Do I have to go in the tcltk directory and run svn up - make? Or do I have to download something somewhere. Again, I got the 6.4 version from Sid but now I switched back on Lenny... argh sorry forgot that. the instructions are cutpaste from the command line from the svn checkout source code dir. would that work the same? you could try to paste the files in, but it's a bit sloppy. Besides, are you suggesting to switch back to the tcltk GUI to make it work?..I kind of liked the wxPython one. well it should now work from the tcltk gui (latest unreleased svn), no idea from the wx GUI. (I still haven't updated from etch here so no official wx2.8 or gui for me except from winGrass or Mac packages) Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
On Mon, 2009-02-23 at 01:05 -0800, Hamish wrote: mitch wrote: GUI: I got stuck at the i.cluster window. How do I set the sampling interval values? I'm not allow to use any kind of character but numbers...how do I differentiate row and column value? http://n2.nabble.com/file/n2369388/sampling.png If I leave 1 value I get this error: ERROR: option sample must be provided in multiples of 2 You provided 1 items: 50 that's a bug, it should let you give two values. --help says: [sample=row_interval,col_interval] got it. I've notice that i.group window has been changed: shouldn't “Name of imagery sub-group” and “Name of the raster map(s) to include in group” be under the Required tab rather than Optional? 'group=name' is the only required option. you can list maps in the group without adding any, so the subgroup and input raster maps are not needed sometimes. also subgroup really is optional. btw (for better or worse) those required/optional groupings happen automatically based on module requirements. I'm not really a fan of that but others think it's a useful feature. shrug I see what you're saying (don't get me wrong, I don't want to criticize in a “bad” way, it's because maybe you receive a lot of feedbacks about other parts and this module is not so usedI haven't seen so many Remote Sensing questions lately...that's why). polite and researched feedback is always welcome! once you are comfortable that an issue is really a bug, please file a ticket, one for each issue. if just left in the mailing list it'll more easily be forgotten in the noise. -Supervised- I couldn't figure out how to make it work. ...but nothing more than that. Shouldn't I be prompted to a display where I can draw my region of interest? (I don't know, maybe I'm expecting something that it's not supposed to be). for sure old xmon based interactive monitors are gone for wxPy on Windows. not sure about what to expect from these modules, but supervised says to me some sort of map interaction is required... GUI: Imagery → Classify image → Interactive input for supervised classification (requires Xterm) it doesn't just need an terminal window for the text based menu data entry, in needs an Xmonitor as well, which is unavailable in the native MS Windows build. a blank Monitor display and a new shell appear and nothing else. we need a better sorry, won't work on MS Windows (so far) error message there. Is the wxPyhton GUI supposed to be run only on Windows? I guess not..right? http://svn.osgeo.org/grass/grass/trunk/gui/wxpython/README because I'm runnig this magic on Lenny ;-) (I actually got the 6.4 version from Sid repository). Are all of the GUI issues due to this http://grass.osgeo.org/wiki/WxPython-based_GUI_for_GRASS#Imagery_tools (Status: development not started yet) perhaps better stated as redevelopment not started in earnest yet, Some behind-the-scenes code has been prepared, but not much. Volunteers welcome! I would if I could but I don't have the required knowledge (I planned to learn a bit of Pyhton though, so maybe one of these days). Hamish Since I'm here I'll use the same post to ask you another thing: is there any way in GRASS to link different map displays and have a brushing (update while dragging) or even a static d.what.rast? Thanks, mitch ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
Hamish: i.class map=spot.comp group=spotspear subgroup=spot.test outsig=test.txt hmmm, if I try the same from gis.m in grass65 I get: = [...] required: NO Enter the name of an existing input file Hit RETURN to cancel request [$GISBASE]/etc/grass-run.sh: line 29: 20596 Segmentation fault $@ ERROR: i.class exited abnormally. Press enter to continue. = and the xterm closes. Ok, it's a bug: the menu item wasn't opening an xmon for you, and wasn't disabling GRASS_RENDER_IMMEDIATE for i.class. For the Tcl/Tk GIS.m GUI now fixed in svn, rev 36103. cd gui/tcltk/ svn up make Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
mitch wrote: GUI: I got stuck at the i.cluster window. How do I set the sampling interval values? I'm not allow to use any kind of character but numbers...how do I differentiate row and column value? http://n2.nabble.com/file/n2369388/sampling.png If I leave 1 value I get this error: ERROR: option sample must be provided in multiples of 2 You provided 1 items: 50 that's a bug, it should let you give two values. --help says: [sample=row_interval,col_interval] I've notice that i.group window has been changed: shouldn't “Name of imagery sub-group” and “Name of the raster map(s) to include in group” be under the Required tab rather than Optional? 'group=name' is the only required option. you can list maps in the group without adding any, so the subgroup and input raster maps are not needed sometimes. also subgroup really is optional. btw (for better or worse) those required/optional groupings happen automatically based on module requirements. I'm not really a fan of that but others think it's a useful feature. shrug (don't get me wrong, I don't want to criticize in a “bad” way, it's because maybe you receive a lot of feedbacks about other parts and this module is not so usedI haven't seen so many Remote Sensing questions lately...that's why). polite and researched feedback is always welcome! once you are comfortable that an issue is really a bug, please file a ticket, one for each issue. if just left in the mailing list it'll more easily be forgotten in the noise. -Supervised- I couldn't figure out how to make it work. ...but nothing more than that. Shouldn't I be prompted to a display where I can draw my region of interest? (I don't know, maybe I'm expecting something that it's not supposed to be). for sure old xmon based interactive monitors are gone for wxPy on Windows. not sure about what to expect from these modules, but supervised says to me some sort of map interaction is required... GUI: Imagery → Classify image → Interactive input for supervised classification (requires Xterm) it doesn't just need an terminal window for the text based menu data entry, in needs an Xmonitor as well, which is unavailable in the native MS Windows build. a blank Monitor display and a new shell appear and nothing else. we need a better sorry, won't work on MS Windows (so far) error message there. Are all of the GUI issues due to this http://grass.osgeo.org/wiki/WxPython-based_GUI_for_GRASS#Imagery_tools (Status: development not started yet) perhaps better stated as redevelopment not started in earnest yet, Some behind-the-scenes code has been prepared, but not much. Volunteers welcome! Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
Hi, let me start saying that I do like the wxpython GUI...it's very nice and usable. It really is!...so, thank you guys for all the work you've been doing for the community! As far as the classification: -Unsupervised- Command: it works either @user or @permanent. GUI: I got stuck at the i.cluster window. How do I set the sampling interval values? I'm not allow to use any kind of character but numbers...how do I differentiate row and column value? http://n2.nabble.com/file/n2369388/sampling.png sampling.png If I leave 1 value I get this error: ERROR: option sample must be provided in multiples of 2 You provided 1 items: 50 I've notice that i.group window has been changed: shouldn't “Name of imagery sub-group” and “Name of the raster map(s) to include in group” be under the Required tab rather than Optional? (don't get me wrong, I don't want to criticize in a “bad” way, it's because maybe you receive a lot of feedbacks about other parts and this module is not so usedI haven't seen so many Remote Sensing questions lately...that's why). -Supervised- I couldn't figure out how to make it work. Command: i.group group=l_sup sub=l_sup in=lsat7_2002_10,lsat7_2002_20,lsat7_2002_30,lsat7_2002_40,lsat7_2002_50,lsat7_2002_70 Adding raster map lsat7_2002...@permanent to group Adding raster map lsat7_2002...@permanent to group Adding raster map lsat7_2002...@permanent to group Adding raster map lsat7_2002...@permanent to group Adding raster map lsat7_2002...@permanent to group Adding raster map lsat7_2002...@permanent to group Adding raster map lsat7_2002...@permanent to subgroup Adding raster map lsat7_2002...@permanent to subgroup Adding raster map lsat7_2002...@permanent to subgroup Adding raster map lsat7_2002...@permanent to subgroup Adding raster map lsat7_2002...@permanent to subgroup Adding raster map lsat7_2002...@permanent to subgroup i.group complete. g.region -p rast=lsat7_2002_40 projection: 99 (Lambert Conformal Conic) zone: 0 datum: nad83 ellipsoid: a=6378137 es=0.006694380022900787 north: 228513 south: 214975.5 west: 629992.5 east: 645012 nsres: 28.5 ewres: 28.5 rows: 475 cols: 527 cells: 250325 Command finished (0 sec) i.class group=l_sup sub=l_sup map=lsat7_2002_40 outsig=outsig Command finished (0 sec) ...but nothing more than that. Shouldn't I be prompted to a display where I can draw my region of interest? (I don't know, maybe I'm expecting something that it's not supposed to be). GUI: Imagery → Classify image → Interactive input for supervised classification (requires Xterm) a blank Monitor display and a new shell appear and nothing else. Are all of the GUI issues due to this http://grass.osgeo.org/wiki/WxPython-based_GUI_for_GRASS#Imagery_tools (Status: development not started yet) Cheers, mitch -- View this message in context: http://n2.nabble.com/Long-post---Remote-Sensing---Classifications-tp2347807p2369388.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
Hi, it's in Raster - Manage map colors - Create RGB there was a wish a while ago to build a menu tree for the commands, not sure if it went anywhere. Maybe something could be built into tools/module_synopsis.sh (add a line of text under each module descr with italicized a-b-c menu location?) todo. hmmm... MODULE=r.composite xml2 gui/wxpython/xml/menudata.xml | grep -B5000 $MODULE | \ sed -e 's+^/menudata/menubar/menu/++' | tac | \ grep '/label=\|^label=' | grep -B1000 '^label=' -m 1 | ... ? may have to call in awk.. ..thanks, Moritz told me where to find it...my bad! if you know it's there but can't find it, what hope does a new user have of learning that it exists? It indicates a non-intuitive placement or where the label could be improved. This sort of feedback is rather valuable IMO. mitch: I'm generally familiar with getting to know new GIS sw and new GUIs..that was my bad (probably not too much attention). Commercial sw give you the possibility to save your visualization once loaded (menu file-save as...and stuff like that)...d.rast doesn't. Maybe I was expecting to find composite as a key word in a first level menu rather that Create RGB in a second level...but again, I wasn't paying too much attention. Once you start understanding that in the GRASS GUI you have to look for a word related to the module, you've got it (i.*it's got to be in that Imagery menu somewhere, doesn't it?) I have to be honest though, the amount of mickeys using GRASS is pretty relevant but, giving its CL nature (virtually everything could be done just typing), this comment doesn't really make sense (unless the direction is really toward a usable GUI as I read in a paper somewhere). Hamish: I would say mitigated not corrected. The imagery modules generally still don't like to look at maps in @othermapsets, and generally don't like @currentmapset tacked on map names either. But hopefully it breaks less now and exits with an understandable error message. The imagery libs need some TLC. Are you saying that it's not a matter of version? no. it should be better in 6.4. Just that even with the new fixes the imagery modules (library) are still not as solid as the raster modules in this regard. Would I have the same problem with the 6.4? I don't know. Hopefully it is fixed, but it might not be. (feedback from users is really the only way to know) I'm sure it doesn't work when I copy the images in the user mapset... would that work working in PERMANENT? try with the imagery maps in the same mapset. perhaps with the maps in PERMANENT you don't need @mapset (it's in the search path), so it would be ok. shrug. I did try to include the images in the group using the graphic way but deleting (in the box) @mapset...no luck. you can look in the $MAPSET/groups/ files too, they are just plain text. for all of the above -- mitch: I installed 6.4 0rc3 and the I realized that for a couple of times unsupervised was working entering @user and using permanent as a pool. Then, same thing (unable to find the signature file...when all the files are there in the group!). So I entered @permanent and worked there...fine (I still don't understand the inconsistency). I entered @user, copied the images and it worked. BTW I used only the GUI (I was testing it) so, a bunch of double clicks generating a bunch of @here and @there in the file naming...that wasn't the problem (at least in my case). But still, how about the Supervised? Wouldn't that work and tell me that THERE IS a sub-group since I used it (the one created via command) for the Usupervised? (I tried to create another group -- via command -- and run the i.class but when I have to define the subgroup...nothing!..it doesn't find it. I don't work with any of the classification stuff, so am not really familiar with them, but if you can recreate the problem(s) from the command line, using the North Carolina or Spearfish sample datasets, please post the commands to trigger it to the bug tracker and we can try and reproduce it at our end. mitch: supervised...oh man, my black ship :-) from GUI it seemed to be working the first time: OPTION: Name of raster map to be displayed key: map format: name required: YES Enter the name of an existing raster file Enter 'list' for a list of existing raster files Hit RETURN to cancel request lsat7_2002_40 lsat7_2002_40 OPTION: Name of input imagery group key: group format: name required: YES Enter the name of an existing group file Enter 'list' for a list of existing group files Hit RETURN to cancel request landsuper landsuper OPTION: Name of input imagery subgroup key: subgroup format: name required: YES enter option landsuper You have chosen: subgroup=landsuper Is this correct? (y/n) [y] OPTION: File to contain result signatures key: outsig format: name required: YES Enter a new output file name Hit RETURN to cancel request land_sup_class land_sup_class
Re: [GRASS-user] Long post - Remote Sensing - Classifications
Hamish wrote: [finding r.composite] mitch_TX wrote: it's in Raster - Manage map colors - Create RGB there was a wish a while ago to build a menu tree for the commands, not sure if it went anywhere. Maybe something could be built into tools/module_synopsis.sh (add a line of text under each module descr with italicized a-b-c menu location?) todo. hmmm... MODULE=r.composite xml2 gui/wxpython/xml/menudata.xml | grep -B5000 $MODULE | \ sed -e 's+^/menudata/menubar/menu/++' | tac | \ grep '/label=\|^label=' | grep -B1000 '^label=' -m 1 | ... ? may have to call in awk.. I don't know what you're trying to do, but I'm 95% certain that it would be easier to do it in Python, probably by cannibalising gui/wxpython/gui_modules/menudata.py. -- 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] Long post - Remote Sensing - Classifications
Oh man...you guys have to excuse me...I've just started the wxpython GUI and it looks pretty cool =)..so, forget all I posted before and let me try it. I'll maybe solve my problems...thanks. Sorry again, mitch Glynn Clements wrote: Hamish wrote: [finding r.composite] mitch_TX wrote: it's in Raster - Manage map colors - Create RGB there was a wish a while ago to build a menu tree for the commands, not sure if it went anywhere. Maybe something could be built into tools/module_synopsis.sh (add a line of text under each module descr with italicized a-b-c menu location?) todo. hmmm... MODULE=r.composite xml2 gui/wxpython/xml/menudata.xml | grep -B5000 $MODULE | \ sed -e 's+^/menudata/menubar/menu/++' | tac | \ grep '/label=\|^label=' | grep -B1000 '^label=' -m 1 | ... ? may have to call in awk.. I don't know what you're trying to do, but I'm 95% certain that it would be easier to do it in Python, probably by cannibalising gui/wxpython/gui_modules/menudata.py. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- View this message in context: http://n2.nabble.com/Long-post---Remote-Sensing---Classifications-tp2347807p2366095.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
Hi Hamish, GRASS cd $GISBASE/etc/gm/ GRASS less gmmenu.tcl /composite [Enter] up up up it's in Raster - Manage map colors - Create RGB there was a wish a while ago to build a menu tree for the commands, not sure if it went anywhere. Maybe something could be built into tools/module_synopsis.sh (add a line of text under each module descr with italicized a-b-c menu location?) [note to self: update the synopsis PDF for 6.4.0] ..thanks, Moritz told me where to find it...my bad! perhaps the computational region is different than the map's? config - region - change region - set region to match this raster map ??? I did that but I'll try it again (I'm not completely sure that I checked the region while I was trying to do it via GUI -- I'm sure I did it via command...in fact good output). I would say mitigated not corrected. The imagery modules generally still don't like to look at maps in @othermapsets, and generally don't like @currentmapset tacked on map names either. But hopefully it breaks less now and exits with an understandable error message. The imagery libs need some TLC. Are you saying that it's not a matter of version? Would I have the same problem with the 6.4? I'm sure it doesn't work when I copy the images in the user mapset...would that work working in PERMANENT? I did try to include the images in the group using the graphic way but deleting (in the box) @mapset...no luck. But still, how about the Supervised? Wouldn't that work and tell me that THERE IS a sub-group since I used it (the one created via command) for the Usupervised? (I tried to create another group -- via command -- and run the i.class but when I have to define the subgroup...nothing!..it doesn't find it. Cheers, mitch. Da: Hamish (via Nabble) ml-user+64140-560766...@n2.nabble.com A: mitch_TX reflex...@yahoo.it Inviato: Mercoledì 18 febbraio 2009, 20:25:52 Oggetto: Re: [GRASS-user] Long post - Remote Sensing - Classifications mitch wrote: How do I actually save a color composite image as a raster? I couldn't find r.composite :working:..too much work. GRASS cd $GISBASE/etc/gm/ GRASS less gmmenu.tcl /composite [Enter] up up up it's in Raster - Manage map colors - Create RGB there was a wish a while ago to build a menu tree for the commands, not sure if it went anywhere. Maybe something could be built into tools/module_synopsis.sh (add a line of text under each module descr with italicized a-b-c menu location?) [note to self: update the synopsis PDF for 6.4.0] The other question is something that has been getting me crazy for a couple of days: classifications. (GUI: Imagery → Classify Image → Clustering input for unsupervised classification) output Reading image ... Iteration 1: % Convergence: 80.00 (0s elapsed, 0s left) Iteration 2: % Convergence: 100.00 (0s elapsed, 0s left) (Pretty short output, innit?) perhaps the computational region is different than the map's? config - region - change region - set region to match this raster map ??? Moritz: There were issues in 6.3 with the GUI appending @MAPSET to names whereas the i.* commands couldn't deal with such names. This was corrected since (but probably after 6.3 came out). I would say mitigated not corrected. The imagery modules generally still don't like to look at maps in @othermapsets, and generally don't like @currentmapset tacked on map names either. But hopefully it breaks less now and exits with an understandable error message. The imagery libs need some TLC. Hamish ___ grass-user mailing list grass-u...@... http://lists.osgeo.org/mailman/listinfo/grass-user This email is a reply to your post @ http://n2.nabble.com/Long-post---Remote-Sensing---Classifications-tp2347807p2350763.html You can reply by email or by visting the link above. Passa a Yahoo! Mail. La webmail che ti offre GRATIS spazio illimitato, antispam e messenger integrato. http://it.mail.yahoo.com/ -- View this message in context: http://n2.nabble.com/Long-post---Remote-Sensing---Classifications-tp2347807p2354014.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
[finding r.composite] mitch_TX wrote: it's in Raster - Manage map colors - Create RGB there was a wish a while ago to build a menu tree for the commands, not sure if it went anywhere. Maybe something could be built into tools/module_synopsis.sh (add a line of text under each module descr with italicized a-b-c menu location?) todo. hmmm... MODULE=r.composite xml2 gui/wxpython/xml/menudata.xml | grep -B5000 $MODULE | \ sed -e 's+^/menudata/menubar/menu/++' | tac | \ grep '/label=\|^label=' | grep -B1000 '^label=' -m 1 | ... ? may have to call in awk.. ..thanks, Moritz told me where to find it...my bad! if you know it's there but can't find it, what hope does a new user have of learning that it exists? It indicates a non-intuitive placement or where the label could be improved. This sort of feedback is rather valuable IMO. Hamish: I would say mitigated not corrected. The imagery modules generally still don't like to look at maps in @othermapsets, and generally don't like @currentmapset tacked on map names either. But hopefully it breaks less now and exits with an understandable error message. The imagery libs need some TLC. Are you saying that it's not a matter of version? no. it should be better in 6.4. Just that even with the new fixes the imagery modules (library) are still not as solid as the raster modules in this regard. Would I have the same problem with the 6.4? I don't know. Hopefully it is fixed, but it might not be. (feedback from users is really the only way to know) I'm sure it doesn't work when I copy the images in the user mapset... would that work working in PERMANENT? try with the imagery maps in the same mapset. perhaps with the maps in PERMANENT you don't need @mapset (it's in the search path), so it would be ok. shrug. I did try to include the images in the group using the graphic way but deleting (in the box) @mapset...no luck. you can look in the $MAPSET/groups/ files too, they are just plain text. But still, how about the Supervised? Wouldn't that work and tell me that THERE IS a sub-group since I used it (the one created via command) for the Usupervised? (I tried to create another group -- via command -- and run the i.class but when I have to define the subgroup...nothing!..it doesn't find it. I don't work with any of the classification stuff, so am not really familiar with them, but if you can recreate the problem(s) from the command line, using the North Carolina or Spearfish sample datasets, please post the commands to trigger it to the bug tracker and we can try and reproduce it at our end. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Long post - Remote Sensing - Classifications
Hi community, This semester I'm teaching an Intro Remote Sensing course using commercial softwares. Since I'm a free spirit, my goal at the end of the semester is to show the students that everything could be done using GRASS rather than Envi,Erdas or Idrisi. Given the fact that they are familiar with GUIs, I'd like to use “graphical” commands as much as possible (I don't want to scare them) I have a few questions: How do I actually save a color composite image as a raster? I know how to create it and visualize it (GUI: Add RGB or HIS layer → I choose the bands → Display active layer on the Map Display) (Command: d.rgb red=lsat5_1987...@landsat green=lsat5_1987...@landsat blue=lsat5_1987...@landsat -o) ...but I don't know how to save it, say, to use it as a base image to draw ground truth for a supervised classification. The other question is something that has been getting me crazy for a couple of days: classifications. (GUI: Imagery → Develop images and groups → Create/edit group. In the i.group window I type the name of the group, sub-group and I include rasters in the group) (Command: i.group group=lsat_try subgroup=lsat_try input=lsat5_1987...@landsat,lsat5_1987...@landsat,lsat5_1987...@landsat,lsat5_1987...@landsat,lsat5_1987...@landsat,lsat5_1987...@landsat) and the output is: group [lsat_try] - does not yet exist. Creating... Adding files to group [lsat_try] Adding raster map [lsat5_1987_10] Adding raster map [lsat5_1987_20] Adding raster map [lsat5_1987_30] Adding raster map [lsat5_1987_40] Adding raster map [lsat5_1987_50] Adding raster map [lsat5_1987_70] Writing group REF Done. Adding files to subgroup [lsat_try] Adding raster map [lsat5_1987_10] Adding raster map [lsat5_1987_20] Adding raster map [lsat5_1987_30] Adding raster map [lsat5_1987_40] Adding raster map [lsat5_1987_50] Adding raster map [lsat5_1987_70] Writing subgroup REF Done. So far so good..right? Now.. (GUI: Imagery → Classify Image → Clustering input for unsupervised classification) In the i.cluster window I type lsat_try for the group, lsat_try for the sub-group, lsat_sign for the signature file, say 15 classes, 50,50 for the sampling, lsat_rep for the report.txt (the rest is by default) (Command: i.cluster group=lsat_try subgroup=lsat_try sigfile=lsat_sign classes=15 sample=50,50 iterations=30 convergence=98.0 separation=0.0 min_size=17 reportfile=lsat_rep.txt) output Reading image ... Iteration 1: % Convergence: 80.00 (0s elapsed, 0s left) Iteration 2: % Convergence: 100.00 (0s elapsed, 0s left) (Pretty short output, innit?) Now.. (GUI: Imagery → Classify Image →Maximum likelyhood classification) I fill the boxes with the required info (Command: i.maxlik group=lsat_try subgroup=lsat_try sigfile=lsat_sign class=class_try reject=class_try_rej) but my output is the following: Can't read signature file [lsat_sign]. The thing that is literally driving me crazy is that sometimes it works (few times) others it doesn't (I generally I get error messages saying that the sub-group doesn't exist in the mapset). I delete everything and I start over using only the command line: i.group group=lsat_try sub=lsat_try in=lsat5_1987_10,lsat5_1987_20,lsat5_1987_30,lsat5_1987_40,lsat5_1987_50,lsat5_1987_70 output: group [lsat_try] - does not yet exist. Creating... Adding files to group [lsat_try] Adding raster map [lsat5_1987_10] Adding raster map [lsat5_1987_20] Adding raster map [lsat5_1987_30] Adding raster map [lsat5_1987_40] Adding raster map [lsat5_1987_50] Adding raster map [lsat5_1987_70] Writing group REF Done. Adding files to subgroup [lsat_try] Adding raster map [lsat5_1987_10] Adding raster map [lsat5_1987_20] Adding raster map [lsat5_1987_30] Adding raster map [lsat5_1987_40] Adding raster map [lsat5_1987_50] Adding raster map [lsat5_1987_70] Writing subgroup REF Done. i.cluster group=lsat_try sub=lsat_try sig=lsat_sign classes=15 report=lsat_rep.txt output Reading image ... Iteration 1: % Convergence: 65.36 (0s elapsed, 0s left) Iteration 2: % Convergence: 71.92 (0s elapsed, 0s left) Iteration 3: % Convergence: 86.67 (0s elapsed, 0s left) Iteration 4: % Convergence: 90.74 (0s elapsed, 0s left) Iteration 5: % Convergence: 93.32 (0s elapsed, 0s left) Iteration 6: % Convergence: 95.15 (0s elapsed, 0s left) Iteration 7: % Convergence: 96.00 (0s elapsed, 0s left) Iteration 8: % Convergence: 96.29 (0s elapsed, 0s left) Iteration 9: % Convergence: 96.70 (0s elapsed, 0s left) Iteration 10: % Convergence: 96.31 (0s elapsed, 0s left) Iteration 11: % Convergence: 96.70 (0s elapsed, 0s left) Iteration 12: % Convergence: 96.68 (0s elapsed, 0s left) Iteration 13: % Convergence: 96.75 (0s elapsed, 0s left) Iteration 14: % Convergence: 96.84 (0s elapsed, 0s left) Iteration 15: % Convergence: 96.98 (0s elapsed, 0s left) Iteration 16: % Convergence: 97.27 (0s elapsed, 0s left) Iteration 17: % Convergence: 97.54 (0s elapsed, 0s left) Iteration 18: %
Re: [GRASS-user] Long post - Remote Sensing - Classifications
I forgot: Debian Lenny fresh installation GRASS 6.3 fresh installation Cheers mitch mitch_TX wrote: Hi community, This semester I'm teaching an Intro Remote Sensing course using commercial softwares. Since I'm a free spirit, my goal at the end of the semester is to show the students that everything could be done using GRASS rather than Envi,Erdas or Idrisi. Given the fact that they are familiar with GUIs, I'd like to use “graphical” commands as much as possible (I don't want to scare them) I have a few questions: How do I actually save a color composite image as a raster? I know how to create it and visualize it (GUI: Add RGB or HIS layer → I choose the bands → Display active layer on the Map Display) (Command: d.rgb red=lsat5_1987...@landsat green=lsat5_1987...@landsat blue=lsat5_1987...@landsat -o) ...but I don't know how to save it, say, to use it as a base image to draw ground truth for a supervised classification. The other question is something that has been getting me crazy for a couple of days: classifications. (GUI: Imagery → Develop images and groups → Create/edit group. In the i.group window I type the name of the group, sub-group and I include rasters in the group) (Command: i.group group=lsat_try subgroup=lsat_try input=lsat5_1987...@landsat,lsat5_1987...@landsat,lsat5_1987...@landsat,lsat5_1987...@landsat,lsat5_1987...@landsat,lsat5_1987...@landsat) and the output is: group [lsat_try] - does not yet exist. Creating... Adding files to group [lsat_try] Adding raster map [lsat5_1987_10] Adding raster map [lsat5_1987_20] Adding raster map [lsat5_1987_30] Adding raster map [lsat5_1987_40] Adding raster map [lsat5_1987_50] Adding raster map [lsat5_1987_70] Writing group REF Done. Adding files to subgroup [lsat_try] Adding raster map [lsat5_1987_10] Adding raster map [lsat5_1987_20] Adding raster map [lsat5_1987_30] Adding raster map [lsat5_1987_40] Adding raster map [lsat5_1987_50] Adding raster map [lsat5_1987_70] Writing subgroup REF Done. So far so good..right? Now.. (GUI: Imagery → Classify Image → Clustering input for unsupervised classification) In the i.cluster window I type lsat_try for the group, lsat_try for the sub-group, lsat_sign for the signature file, say 15 classes, 50,50 for the sampling, lsat_rep for the report.txt (the rest is by default) (Command: i.cluster group=lsat_try subgroup=lsat_try sigfile=lsat_sign classes=15 sample=50,50 iterations=30 convergence=98.0 separation=0.0 min_size=17 reportfile=lsat_rep.txt) output Reading image ... Iteration 1: % Convergence: 80.00 (0s elapsed, 0s left) Iteration 2: % Convergence: 100.00 (0s elapsed, 0s left) (Pretty short output, innit?) Now.. (GUI: Imagery → Classify Image →Maximum likelyhood classification) I fill the boxes with the required info (Command: i.maxlik group=lsat_try subgroup=lsat_try sigfile=lsat_sign class=class_try reject=class_try_rej) but my output is the following: Can't read signature file [lsat_sign]. The thing that is literally driving me crazy is that sometimes it works (few times) others it doesn't (I generally I get error messages saying that the sub-group doesn't exist in the mapset). I delete everything and I start over using only the command line: i.group group=lsat_try sub=lsat_try in=lsat5_1987_10,lsat5_1987_20,lsat5_1987_30,lsat5_1987_40,lsat5_1987_50,lsat5_1987_70 output: group [lsat_try] - does not yet exist. Creating... Adding files to group [lsat_try] Adding raster map [lsat5_1987_10] Adding raster map [lsat5_1987_20] Adding raster map [lsat5_1987_30] Adding raster map [lsat5_1987_40] Adding raster map [lsat5_1987_50] Adding raster map [lsat5_1987_70] Writing group REF Done. Adding files to subgroup [lsat_try] Adding raster map [lsat5_1987_10] Adding raster map [lsat5_1987_20] Adding raster map [lsat5_1987_30] Adding raster map [lsat5_1987_40] Adding raster map [lsat5_1987_50] Adding raster map [lsat5_1987_70] Writing subgroup REF Done. i.cluster group=lsat_try sub=lsat_try sig=lsat_sign classes=15 report=lsat_rep.txt output Reading image ... Iteration 1: % Convergence: 65.36 (0s elapsed, 0s left) Iteration 2: % Convergence: 71.92 (0s elapsed, 0s left) Iteration 3: % Convergence: 86.67 (0s elapsed, 0s left) Iteration 4: % Convergence: 90.74 (0s elapsed, 0s left) Iteration 5: % Convergence: 93.32 (0s elapsed, 0s left) Iteration 6: % Convergence: 95.15 (0s elapsed, 0s left) Iteration 7: % Convergence: 96.00 (0s elapsed, 0s left) Iteration 8: % Convergence: 96.29 (0s elapsed, 0s left) Iteration 9: % Convergence: 96.70 (0s elapsed, 0s left) Iteration 10: % Convergence: 96.31 (0s elapsed, 0s left) Iteration 11: % Convergence: 96.70 (0s elapsed, 0s left) Iteration 12: % Convergence: 96.68 (0s elapsed, 0s left) Iteration
Re: [GRASS-user] Long post - Remote Sensing - Classifications
On 18/02/09 17:17, mitch_TX wrote: Hi community, This semester I'm teaching an Intro Remote Sensing course using commercial softwares. Since I'm a free spirit, my goal at the end of the semester is to show the students that everything could be done using GRASS rather than Envi,Erdas or Idrisi. Given the fact that they are familiar with GUIs, I'd like to use “graphical” commands as much as possible (I don't want to scare them) I have a few questions: How do I actually save a color composite image as a raster? I know how to create it and visualize it (GUI: Add RGB or HIS layer → I choose the bands → Display active layer on the Map Display) (Command: d.rgb red=lsat5_1987...@landsat green=lsat5_1987...@landsat blue=lsat5_1987...@landsat -o) ...but I don't know how to save it, say, to use it as a base image to draw ground truth for a supervised classification. r.composite (Raster-Manage map colors-Create RGB) The other question is something that has been getting me crazy for a couple of days: classifications. [...] Is this a bug or am I “the bug” :-) ?? I know that if the command line works, you all may not “see what the problem is”? But as I said I'd like to use the GUI. There were issues in 6.3 with the GUI appending @MAPSET to names whereas the i.* commands couldn't deal with such names. This was corrected since (but probably after 6.3 came out). You might want to try with 6.4 (IIUC, there should be a Lenny package available in the gfoss.it repository soon). Let's try with the Supervised using the group and sub-group that worked before: (GUI: Imagery → Classify image → Interactive input for supervised classifications) i.class window, say I want to use the NIR band to draw my truths: lsat5_1987_40 (good) group lsat_try (the group that worked before) (good) sub (if I type list to see the available sub-groups I get: no subgruop available in the current mapset) the same happens if I create a group and sub in the user mapset using the copied images. What am I doing wrong? :,( Can't help you on this one, but maybe same issue of file naming ? Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
Thanks Moritz, I couldn't find r.composite :working:..too much work. As far as the 6.4, I'm looking forward to giving it a shot. Thanks again, mitch Moritz Lennert wrote: On 18/02/09 17:17, mitch_TX wrote: Hi community, This semester I'm teaching an Intro Remote Sensing course using commercial softwares. Since I'm a free spirit, my goal at the end of the semester is to show the students that everything could be done using GRASS rather than Envi,Erdas or Idrisi. Given the fact that they are familiar with GUIs, I'd like to use “graphical” commands as much as possible (I don't want to scare them) I have a few questions: How do I actually save a color composite image as a raster? I know how to create it and visualize it (GUI: Add RGB or HIS layer → I choose the bands → Display active layer on the Map Display) (Command: d.rgb red=lsat5_1987...@landsat green=lsat5_1987...@landsat blue=lsat5_1987...@landsat -o) ...but I don't know how to save it, say, to use it as a base image to draw ground truth for a supervised classification. r.composite (Raster-Manage map colors-Create RGB) The other question is something that has been getting me crazy for a couple of days: classifications. [...] Is this a bug or am I “the bug” :-) ?? I know that if the command line works, you all may not “see what the problem is”? But as I said I'd like to use the GUI. There were issues in 6.3 with the GUI appending @MAPSET to names whereas the i.* commands couldn't deal with such names. This was corrected since (but probably after 6.3 came out). You might want to try with 6.4 (IIUC, there should be a Lenny package available in the gfoss.it repository soon). Let's try with the Supervised using the group and sub-group that worked before: (GUI: Imagery → Classify image → Interactive input for supervised classifications) i.class window, say I want to use the NIR band to draw my truths: lsat5_1987_40 (good) group lsat_try (the group that worked before) (good) sub (if I type list to see the available sub-groups I get: no subgruop available in the current mapset) the same happens if I create a group and sub in the user mapset using the copied images. What am I doing wrong? :,( Can't help you on this one, but maybe same issue of file naming ? Moritz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- View this message in context: http://n2.nabble.com/Long-post---Remote-Sensing---Classifications-tp2347807p2348306.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Long post - Remote Sensing - Classifications
mitch wrote: How do I actually save a color composite image as a raster? I couldn't find r.composite :working:..too much work. GRASS cd $GISBASE/etc/gm/ GRASS less gmmenu.tcl /composite [Enter] up up up it's in Raster - Manage map colors - Create RGB there was a wish a while ago to build a menu tree for the commands, not sure if it went anywhere. Maybe something could be built into tools/module_synopsis.sh (add a line of text under each module descr with italicized a-b-c menu location?) [note to self: update the synopsis PDF for 6.4.0] The other question is something that has been getting me crazy for a couple of days: classifications. (GUI: Imagery → Classify Image → Clustering input for unsupervised classification) output Reading image ... Iteration 1: % Convergence: 80.00 (0s elapsed, 0s left) Iteration 2: % Convergence: 100.00 (0s elapsed, 0s left) (Pretty short output, innit?) perhaps the computational region is different than the map's? config - region - change region - set region to match this raster map ??? Moritz: There were issues in 6.3 with the GUI appending @MAPSET to names whereas the i.* commands couldn't deal with such names. This was corrected since (but probably after 6.3 came out). I would say mitigated not corrected. The imagery modules generally still don't like to look at maps in @othermapsets, and generally don't like @currentmapset tacked on map names either. But hopefully it breaks less now and exits with an understandable error message. The imagery libs need some TLC. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user