Re: [GRASS-user] Removing small areas using v.clean
On Wed, Sep 21, 2011 at 12:51 AM, Pierre Roudier pierre.roud...@gmail.com wrote: Markus, Sorry for this late reply as I was away from keyboard the last few days, I encountered this problem while testing the latest grass_svn (grass 7). I can give it a shot using a stable version if you want, You can try, but I would be surprised if the result is any different. Can you make the vector available for testing, or can you reproduce this with a vector of the sample datasets? Markus M 2011/9/17 Markus Metz markus.metz.gisw...@googlemail.com: On Fri, Sep 16, 2011 at 6:28 AM, Pierre Roudier pierre.roud...@gmail.com wrote: Dear list, I got a problem as I'm trying to remove small areas of a vector layer using v.clean with the rmarea tool: some of these small areas that are on the boundary of the layer are getting erased rather than being merged with the closest polygon. Is that an expected behaviour? No, they should always be merged with a neighbouring area along the shared boundary. Areas are supposed to be erased only if they are isolated, i.e. do not share a boundary with another area. What GRASS version are you using? There was a bug which was fixed 2 years ago. Markus M -- Scientist Landcare Research, New Zealand ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] simultaneous raster maps export
Hi to all, i try to export, together in one file, 2 or 3 raster maps displayed simultaneously in Grass, so they will remain referenced together and furtherly used seperately or together. Does anybody has any idea if this can happen and how (also to which format)? In ArcGis it can be done in vrml2 (.wrl) format easily. Thanks. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] albedo map
Hi, I want to generate an annual albedo map, but I never work with satellite images…is it possible with GRASS? In which way can I start? Please help me…___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Open Source GIS Consultancy
Dear sir/madam; I am writing to you for my company which is a large group on GIS sector in Turkey. We are looking for an open source gis expert. We need a consultant to guide us the best way for open source GIS. (we want to work with grass, mapserver and oracle. ) We already have systems on mapinfo and oracle. And willing to change it to open source so we can also contribute to the open source GIS universe. Kindest Regards Serdar SerdarBuyuktemiz image001.jpg___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
RE: [GRASS-user] Open Source GIS Consultancy
Don't forget QGIS and Postgres/PostGIS. Jim Burton Promap/GIS Business Services National Grid Property Ltd National Grid Workspace: A1 Property National Grid House, Warwick Technology Park Gallows Hill, Warwick, CV34 6DA Tel: (Internal) 7474 5432 (External) +44 (0)1926 65 5432 Mob: Address: jim.bur...@uk.ngrid.com From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of Serdar BUYUKTEMIZ Sent: Wednesday, September 21, 2011 2:21 PM To: grass-user@lists.osgeo.org Subject: [GRASS-user] Open Source GIS Consultancy Dear sir/madam; I am writing to you for my company which is a large group on GIS sector in Turkey. We are looking for an open source gis expert. We need a consultant to guide us the best way for open source GIS. (we want to work with grass, mapserver and oracle. ) We already have systems on mapinfo and oracle. And willing to change it to open source so we can also contribute to the open source GIS universe... Kindest Regards Serdar This e-mail, and any attachments are strictly confidential and intended for the addressee(s) only. The content may also contain legal, professional or other privileged information. If you are not the intended recipient, please notify the sender immediately and then delete the e-mail and any attachments. You should not disclose, copy or take any action in reliance on this transmission. You may report the matter by calling us on 0800 085 4806. Please ensure you have adequate virus protection before you open or detach any documents from this transmission. The Group Companies do not accept any liability for viruses. An e-mail reply to this address may be subject to monitoring for operational reasons or lawful business practices. For the registered information on the UK operating companies within the National Grid group please use the attached link: http://www.nationalgrid.com/corporate/legal/registeredoffices.htmimage001.jpg___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] GRASS and GRASS-DEV Package
Hi, I have GRASS 6.4.1 and Ubuntu 10.10 (Maverick). I want to install an addons called 'r.seg'. I have tried to install this via SVN but I doesn't work. After various google research, I've understand that I don't have grass-dev package installed!I have searched around but didn't found anything, just this: http://packages.ubuntu.com/maverick/grass-dev http://packages.ubuntu.com/maverick/grass-dev I've tried to install but is not compatible!Can't install! Do you know where I can find the grass-dev package for my GRASS and Ubuntu version? One more stupid thing; If I install an addons on GRASS, after that, can I use the same on the GRASS Plugin in QGIS? Thank you very much! Solimyr -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/GRASS-and-GRASS-DEV-Package-tp6816343p6816343.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
[GRASS-user] Error in MAN page in G.manual
Hi I'm trying to launch a manual page using g.manual. It works ok but If I swelect -m flag I get : ERROR: no MAN page entry for r.in.gdal. Why do I get this MAN page error for r.in.gdal if r.in.gdal.html is available? Thanks Luisa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error in MAN page in G.manual
2011/9/21 Luisa Peña luisapena1...@gmail.com: I'm trying to launch a manual page using g.manual. It works ok but If I swelect -m flag I get : ERROR: no MAN page entry for r.in.gdal. Why do I get this MAN page error for r.in.gdal if r.in.gdal.html is available? version ? OS? Martin -- Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Where GRASS_DOC_BASE variable is stored?
Greetings I would like to know where GRASS_DOC_BASE is stored and/or how is builted? THanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6482 (20110921) __ 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] problems with r.to.rast3elev
Thanks Soeren, As a side note, it did run on my Mac, using Kygschaos package (6.4.1, I guess)... here is the output if g.region -p3: GRASS 7.0.svn (santa_catarina@dunas2):~ g.region -p3 projection: 1 (UTM) zone: -22 datum: wgs84 ellipsoid: wgs84 north: 6901400 south: 6899300 west: 731900 east: 732900 top:60. bottom: 0. nsres: 0.5 nsres3: 0.5 ewres: 0.5 ewres3: 0.5 tbres: 0.5 rows: 4200 rows3: 4200 cols: 2000 cols3: 2000 depths: 120 cells: 840 cells3: 100800 and here of gdb: GRASS 7.0.svn (santa_catarina@dunas2):~ r.to.rast3elev --overwrite input=mask_volume_garopaba elevation=garopaba_SOLO_bicubic output=garopaba_vol05m lower=1 Creating 3D raster map Segmentation fault GRASS 7.0.svn (santa_catarina@dunas2):~ gdb `which r.to.rast3elev` GNU gdb (Ubuntu/Linaro 7.2-1ubuntu11) 7.2 Copyright (C) 2010 Free Software Foundation, Inc. License GPLv3+: GNU GPL version 3 or later http://gnu.org/licenses/gpl.html This is free software: you are free to change and redistribute it. There is NO WARRANTY, to the extent permitted by law. Type show copying and show warranty for details. This GDB was configured as i686-linux-gnu. For bug reporting instructions, please see: http://www.gnu.org/software/gdb/bugs/... Reading symbols from /usr/local/grass-7.0.svn/bin/r.to.rast3elev...(no debugging symbols found)...done. (gdb) run input=mask_volume_garopaba elevation=garopaba_SOLO_bicubic output=garopaba_vol05m lower=1 Starting program: /usr/local/grass-7.0.svn/bin/r.to.rast3elev input=mask_volume_garopaba elevation=garopaba_SOLO_bicubic output=garopaba_vol05m lower=1 Creating 3D raster map 0% Program received signal SIGSEGV, Segmentation fault. 0xb7fb5ae4 in Rast_set_d_null_value () from /usr/local/grass-7.0.svn/lib/ libgrass_raster.7.0.svn.so (gdb) bt full #0 0xb7fb5ae4 in Rast_set_d_null_value () from /usr/local/grass-7.0.svn/lib/libgrass_raster.7.0.svn.so No symbol table info available. #1 0xb7fd3611 in Rast3d_set_null_value () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so No symbol table info available. #2 0xb7fd6a92 in Rast3d_set_null_tile_type () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so No symbol table info available. #3 0xb7fd6e3e in Rast3d_read_tile () from /usr/local/grass-7.0.svn/lib/ libgrass_g3d.7.0.svn.so No symbol table info available. #4 0xb7fc9288 in ?? () from /usr/local/grass-7.0.svn/lib/ libgrass_g3d.7.0.svn.so No symbol table info available. #5 0xb7fc89fb in Rast3d_cache_elt_ptr () from /usr/local/grass-7.0.svn/lib/ libgrass_g3d.7.0.svn.so No symbol table info available. #6 0xb7fd642d in Rast3d_get_tile_ptr () from /usr/local/grass-7.0.svn/lib/ libgrass_g3d.7.0.svn.so No symbol table info available. #7 0xb7fd4761 in Rast3d_put_double () from /usr/local/grass-7.0.svn/lib/ libgrass_g3d.7.0.svn.so No symbol table info available. #8 0x08049880 in elev_raster_to_g3d () No symbol table info available. #9 0x08049e6e in main () No symbol table info available. (gdb) l No symbol table is loaded. Use the file command. (gdb) frame 2 #2 0xb7fd6a92 in Rast3d_set_null_tile_type () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so (gdb) l No symbol table is loaded. Use the file command. (gdb) best, Carlos 2011/9/21 Sören Gebbert soerengebb...@googlemail.com Hi Carlos, can you please provide more information about the grass7 segfault? Can you please print your region settings (g.region -p3)? Can you please use gdb or valgrind to check at which line in the code the segfault appears? I will try to fix this issue in grass7, but i cant reproduce it on my system. Best regards Soeren 2011/9/20 Carlos Grohmann carlos.grohm...@gmail.com: Hello all, I'm experiencing some issues with r.to.rast3elev , that I didn't had before. I'm trying to reprocess some volume calculation on dune fields in southern Brazil (see this: http://geomorphometry.org/Grohmann2011) so I can evolve what I presented in the Geomorphometry Meeting into a full paper, but at this point I can't calculate the volumes any more! (BTW, Helena, I'm thinking about comparing the volumes from r.volume and r.to.rast3elev, since I think there will be differences..) OS: Linux, Ubuntu 11.04 In GRASS 6.4.1, installed from Ubuntu repositories, I got this (running inside a python session): grass.run_command('r.to.rast3elev', input=clump, elevation=dem, output=dem3d, lower=1, overwrite=True) Creating 3D raster map ERROR: G3d_cache_hash_remove_name: name not in hashtable and in GRASS 7.0svn updated today, I got a segmentation fault. the region is set to match the 'clump' map (which is the same as the mask, to limit the calculations to the active dunes). any help is appreciated Carlos -- Prof. Carlos Henrique Grohmann - Geologist D.Sc. Institute of Geosciences - Univ. of São Paulo, Brazil
[GRASS-user] Re: GRASS and GRASS-DEV Package (Alexander Muriy)
*Solimyr* wrote: http://packages.ubuntu.com/maverick/grass-dev http://packages.ubuntu.com/maverick/grass-dev I've tried to install but is not compatible!Can't install! Do you know where I can find the grass-dev package for my GRASS and Ubuntu version? From my modest experience: to install add-on normally you, in general, must install GRASS from source and then install add-on as noted on http://grass.osgeo.org/wiki/Compile_and_Install#Addons page. Other things like grass-dev packages and so on -- IMHO just wishes how it will be.. If I'm wrong, correct me. One more stupid thing; If I install an addons on GRASS, after that, can I use the same on the GRASS Plugin in QGIS? In simple way, as I know, NO. But it's possible to add modules to GRASS Plugin in QGIS with the special method. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problems with r.to.rast3elev
Hi Carlos, please try svn revision 48398 of grass7. I think the reason was the wrong computation of the tile size, because the chosen variable types are not large enough. This should be fixed now. I hope it works now. Best regards Soeren 2011/9/21 Carlos Grohmann carlos.grohm...@gmail.com: Thanks Soeren, As a side note, it did run on my Mac, using Kygschaos package (6.4.1, I guess)... here is the output if g.region -p3: GRASS 7.0.svn (santa_catarina@dunas2):~ g.region -p3 projection: 1 (UTM) zone: -22 datum: wgs84 ellipsoid: wgs84 north: 6901400 south: 6899300 west: 731900 east: 732900 top: 60. bottom: 0. nsres: 0.5 nsres3: 0.5 ewres: 0.5 ewres3: 0.5 tbres: 0.5 rows: 4200 rows3: 4200 cols: 2000 cols3: 2000 depths: 120 cells: 840 cells3: 100800 and here of gdb: GRASS 7.0.svn (santa_catarina@dunas2):~ r.to.rast3elev --overwrite input=mask_volume_garopaba elevation=garopaba_SOLO_bicubic output=garopaba_vol05m lower=1 Creating 3D raster map Segmentation fault GRASS 7.0.svn (santa_catarina@dunas2):~ gdb `which r.to.rast3elev` GNU gdb (Ubuntu/Linaro 7.2-1ubuntu11) 7.2 Copyright (C) 2010 Free Software Foundation, Inc. License GPLv3+: GNU GPL version 3 or later http://gnu.org/licenses/gpl.html This is free software: you are free to change and redistribute it. There is NO WARRANTY, to the extent permitted by law. Type show copying and show warranty for details. This GDB was configured as i686-linux-gnu. For bug reporting instructions, please see: http://www.gnu.org/software/gdb/bugs/... Reading symbols from /usr/local/grass-7.0.svn/bin/r.to.rast3elev...(no debugging symbols found)...done. (gdb) run input=mask_volume_garopaba elevation=garopaba_SOLO_bicubic output=garopaba_vol05m lower=1 Starting program: /usr/local/grass-7.0.svn/bin/r.to.rast3elev input=mask_volume_garopaba elevation=garopaba_SOLO_bicubic output=garopaba_vol05m lower=1 Creating 3D raster map 0% Program received signal SIGSEGV, Segmentation fault. 0xb7fb5ae4 in Rast_set_d_null_value () from /usr/local/grass-7.0.svn/lib/libgrass_raster.7.0.svn.so (gdb) bt full #0 0xb7fb5ae4 in Rast_set_d_null_value () from /usr/local/grass-7.0.svn/lib/libgrass_raster.7.0.svn.so No symbol table info available. #1 0xb7fd3611 in Rast3d_set_null_value () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so No symbol table info available. #2 0xb7fd6a92 in Rast3d_set_null_tile_type () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so No symbol table info available. #3 0xb7fd6e3e in Rast3d_read_tile () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so No symbol table info available. #4 0xb7fc9288 in ?? () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so No symbol table info available. #5 0xb7fc89fb in Rast3d_cache_elt_ptr () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so No symbol table info available. #6 0xb7fd642d in Rast3d_get_tile_ptr () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so No symbol table info available. #7 0xb7fd4761 in Rast3d_put_double () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so No symbol table info available. #8 0x08049880 in elev_raster_to_g3d () No symbol table info available. #9 0x08049e6e in main () No symbol table info available. (gdb) l No symbol table is loaded. Use the file command. (gdb) frame 2 #2 0xb7fd6a92 in Rast3d_set_null_tile_type () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so (gdb) l No symbol table is loaded. Use the file command. (gdb) best, Carlos 2011/9/21 Sören Gebbert soerengebb...@googlemail.com Hi Carlos, can you please provide more information about the grass7 segfault? Can you please print your region settings (g.region -p3)? Can you please use gdb or valgrind to check at which line in the code the segfault appears? I will try to fix this issue in grass7, but i cant reproduce it on my system. Best regards Soeren 2011/9/20 Carlos Grohmann carlos.grohm...@gmail.com: Hello all, I'm experiencing some issues with r.to.rast3elev , that I didn't had before. I'm trying to reprocess some volume calculation on dune fields in southern Brazil (see this: http://geomorphometry.org/Grohmann2011) so I can evolve what I presented in the Geomorphometry Meeting into a full paper, but at this point I can't calculate the volumes any more! (BTW, Helena, I'm thinking about comparing the volumes from r.volume and r.to.rast3elev, since I think there will be differences..) OS: Linux, Ubuntu 11.04 In GRASS 6.4.1, installed from Ubuntu repositories, I got this (running inside a python session): grass.run_command('r.to.rast3elev', input=clump, elevation=dem, output=dem3d, lower=1, overwrite=True) Creating 3D raster map ERROR:
Re: [GRASS-user] Removing small areas using v.clean
Thanks Markus, Here's an archive containing the input vector layer, along with v.clean output: https://df.bestgrid.org/quickshare/fb0bb3af778c8a03/vectors.tar.gz Here's the result of v.info -h on the output: COMMAND: r.to.vect input=c_means output=foo feature=area GISDBASE: /data/GRASSDATA LOCATION: NZTM2000 MAPSET: BH USER: roudierp DATE: Thu Sep 15 17:29:53 2011 - COMMAND: v.clean input=foo layer=1 type=point,line,boundary,centroid,area,face,kernel output=bar tool=rmar ea thresh=900 GISDBASE: /data/GRASSDATA LOCATION: NZTM2000 MAPSET: BH USER: roudierp DATE: Thu Sep 15 17:31:45 2011 [Raster MASK present] Cheers, Pierre 2011/9/21 Markus Metz markus.metz.gisw...@googlemail.com: On Wed, Sep 21, 2011 at 12:51 AM, Pierre Roudier pierre.roud...@gmail.com wrote: Markus, Sorry for this late reply as I was away from keyboard the last few days, I encountered this problem while testing the latest grass_svn (grass 7). I can give it a shot using a stable version if you want, You can try, but I would be surprised if the result is any different. Can you make the vector available for testing, or can you reproduce this with a vector of the sample datasets? Markus M 2011/9/17 Markus Metz markus.metz.gisw...@googlemail.com: On Fri, Sep 16, 2011 at 6:28 AM, Pierre Roudier pierre.roud...@gmail.com wrote: Dear list, I got a problem as I'm trying to remove small areas of a vector layer using v.clean with the rmarea tool: some of these small areas that are on the boundary of the layer are getting erased rather than being merged with the closest polygon. Is that an expected behaviour? No, they should always be merged with a neighbouring area along the shared boundary. Areas are supposed to be erased only if they are isolated, i.e. do not share a boundary with another area. What GRASS version are you using? There was a bug which was fixed 2 years ago. Markus M -- Scientist Landcare Research, New Zealand -- Scientist Landcare Research, New Zealand ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] GRASS and GRASS-DEV Package
Solimyr wrote: I have GRASS 6.4.1 and Ubuntu 10.10 (Maverick). I want to install an addons called 'r.seg'. I have tried to install this via SVN but I doesn't work. After various google research, I've understand that I don't have grass-dev package installed!I have searched around but didn't found anything, just this: http://packages.ubuntu.com/maverick/grass-dev http://packages.ubuntu.com/maverick/grass-dev I've tried to install but is not compatible!Can't install! Do you know where I can find the grass-dev package for my GRASS and Ubuntu version? just open a terminal and run: sudo apt-get install grass-dev but it is broken right now; unless you are a Makefile guru you need the full grass source code if you want to compile GRASS addons in Ubuntu currently. If it were a python or shell script you could just copy the script into the right place and make it executable, but as r.seg is a C program it will need to be compiled. In the last few days I have been working on fixing this in the DebianGIS package building rules and the g.extension(.sh) script and am very very close to having it function properly there, docs and all. (and thus for future ubuntus too) @devs: current problem is that Grass.make's INST_NOW=y works, but wants to run make a second time. the first time it all builds ok in the user's .tmp dir, but the (unneeded?!) 2nd pass it breaks looking for /usr/.../dist.x86_64-unknown-linux-gnu/ which it doesn't have permission to create/write to. One more stupid thing; If I install an addons on GRASS, after that, can I use the same on the GRASS Plugin in QGIS? yes, but not automatically. You need to write a .qgm file and stash it in the right place. The only grass addon which has that support AFAIK is v.surf.icw, although I haven't tested with a modern version of QGIS. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Where GRASS_DOC_BASE variable is stored?
António wrote: I would like to know where GRASS_DOC_BASE is stored and/or how is builted? It is only in grass 6.x scripts/g.manual, for most cases it will just be the same as $GISBASE. do you get an error that the manuals are not found? if so, what does the error look like? If using debian or ubuntu, is the grass-doc package installed? Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Removing small areas using v.clean
BTW I confirm I got similar problems with the 6.4 series, Pierre 2011/9/22 Pierre Roudier pierre.roud...@gmail.com: Thanks Markus, Here's an archive containing the input vector layer, along with v.clean output: https://df.bestgrid.org/quickshare/fb0bb3af778c8a03/vectors.tar.gz Here's the result of v.info -h on the output: COMMAND: r.to.vect input=c_means output=foo feature=area GISDBASE: /data/GRASSDATA LOCATION: NZTM2000 MAPSET: BH USER: roudierp DATE: Thu Sep 15 17:29:53 2011 - COMMAND: v.clean input=foo layer=1 type=point,line,boundary,centroid,area,face,kernel output=bar tool=rmar ea thresh=900 GISDBASE: /data/GRASSDATA LOCATION: NZTM2000 MAPSET: BH USER: roudierp DATE: Thu Sep 15 17:31:45 2011 [Raster MASK present] Cheers, Pierre 2011/9/21 Markus Metz markus.metz.gisw...@googlemail.com: On Wed, Sep 21, 2011 at 12:51 AM, Pierre Roudier pierre.roud...@gmail.com wrote: Markus, Sorry for this late reply as I was away from keyboard the last few days, I encountered this problem while testing the latest grass_svn (grass 7). I can give it a shot using a stable version if you want, You can try, but I would be surprised if the result is any different. Can you make the vector available for testing, or can you reproduce this with a vector of the sample datasets? Markus M 2011/9/17 Markus Metz markus.metz.gisw...@googlemail.com: On Fri, Sep 16, 2011 at 6:28 AM, Pierre Roudier pierre.roud...@gmail.com wrote: Dear list, I got a problem as I'm trying to remove small areas of a vector layer using v.clean with the rmarea tool: some of these small areas that are on the boundary of the layer are getting erased rather than being merged with the closest polygon. Is that an expected behaviour? No, they should always be merged with a neighbouring area along the shared boundary. Areas are supposed to be erased only if they are isolated, i.e. do not share a boundary with another area. What GRASS version are you using? There was a bug which was fixed 2 years ago. Markus M -- Scientist Landcare Research, New Zealand -- Scientist Landcare Research, New Zealand -- Scientist Landcare Research, New Zealand ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user