Re: [GRASS-user] Add DB entries to vector line (after v.split)
v.to.equidist https://raw.githubusercontent.com/amuriy/GRASS-scripts/master/v.to.equidist addon doesn't work under GRASS 7 (only 6.4.*). But I still have plans to rewrite it in Python :) You could also try the addon v.fixed.segmentpoints, it works similary. 2015-05-08 20:10 GMT+03:00 Ken Nussear knuss...@mac.com: Thanks very much, that is all working well. I did try the addon you mention, but using grass7 it wont install in either linux or mac os I'm having some trouble with the final layer though when trying to sample a raster Using this command v.what.rast map=TSELinesSplit500_w_lay2@PERMANENT type=centroid raster=a250m_ATI@MojaveCovariates_250m column=RC1 WARNING: No record for category 0 in table TSELinesSplit500_w_lay2 Any idea why I have no category 0? Ken Moritz Lennert mlenn...@club.worldonline.be May 7, 2015 at 2:54 AM On 06/05/15 22:12, knussear wrote: Interestingly 4 years later and using V.split in grass7 and I'm still having the same problem tried the v.category solution and the cats were never deleted. Are you sure ? In the original thread, the check for success was wrong [1]: v.db.select checks the attached attribute table (which v.split does not modify), not the categories present in the features. Here's a working example using a simple vector map with three lines: v.split test_line out=line_test_split length=5000 v.category line_test_split op=report Layer/table: 1/line_test_split type countminmax point 0 0 0 line 46 2 4 boundary 0 0 0 centroid 0 0 0 area 0 0 0 face 0 0 0 kernel 0 0 0 all 46 2 4 Now, two options: 1) You want to keep the category values of the original lines, but also add individual category values for the new segements. Here GRASS's layer paradigm comes in handy, allowing you to create a second layer of information: #create second information layer containing new individual category values for the segments v.category line_test_split out=line_test_split_w_lay2 layer=2 op=add v.category line_test_split_w_lay2 op=report Layer/table: 1/line_test_split_w_lay2 type countminmax point 0 0 0 line 46 2 4 boundary 0 0 0 centroid 0 0 0 area 0 0 0 face 0 0 0 kernel 0 0 0 all 46 2 4 Layer/table: 2/line_test_split_w_lay2_2 type countminmax point 0 0 0 line 46 1 46 boundary 0 0 0 centroid 0 0 0 area 0 0 0 face 0 0 0 kernel 0 0 0 all 46 1 46 Now you can add a new attribute table to the second layer: v.db.addtable line_test_split_w_lay2 layer=2 If you want to integrate the old categories (or anything else from the original attribute table) into the attribute table of the new segments: v.db.addcolumn line_test_split_w_lay2 col=old_cat int layer=2 v.to.db line_test_split_w_lay2 layer=2 query_layer=1 op=query query_col=cat col=old_cat v.db.select line_test_split_w_lay2 layer=2 cat|old_cat 1|2 2|2 3|2 4|2 5|2 6|2 7|2 8|2 9|2 10|2 11|2 12|2 13|2 14|2 15|2 16|2 17|3 18|3 19|3 20|3 21|3 etc... 2) You just need the segments with individual category values, but not the old line categories, nor the old attribute table: #delete old cats, don't keep old attribute table to avoid confusion v.category -t line_test_split op=del cat=-1 out=line_test_split_wo_cats #add new cats v.category line_test_split_wo_cats out=line_test_split_w_new_cats op=add #now, if necessary create new attribute table prefilled with cat values v.db.addtable line_test_split_w_new_cats My system will also not install the v.to.equidist add on. Which version of GRASS are you using. The addon was developed for GRASS 6. It is a simple shell script which you can download from [2], but I haven't checked if it will run as-is in GRASS 7 Moritz [1] http://lists.osgeo.org/pipermail/grass-user/2011-October/062318.html [2] https://raw.githubusercontent.com/amuriy/GRASS-scripts/master/v.to.equidist ___ 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] Large vector file processing
Hello, Margherita! I don't know about limitations of raster or memory sizes in GRASS, but the first thought was to cut raster into parts, vectorize them and then patch into new vector map. Of course, with automatic cutting, not manual :) 2015-04-24 10:54 GMT+03:00 Margherita Di Leo direg...@gmail.com: Hi, I'm trying to vectorize a very large file | Type of Map: raster Number of Categories: 255 | | Data Type:CELL | | Rows: 88537 | | Columns: 54410 | | Total Cells: 4817298170 | |Projection: Lambert Azimuthal Equal Area | |N:4239550S:2026125 Res:25 | |E:4284650W:2924400 Res:25 | | Range of data:min = 11 max = 255 using GRASS 7.1.svn (r65027) :~ r.to.vect input=raster_tty output=vector_tty type=area Extracting areas... 100% Writing areas... 100% Building topology for vector map vector_tty@assessment... Registering primitives... Killed what does this mean? :~ v.build vector_tty WARNING: Coor file of vector map vector_tty@assessment is larger than it should be (2403996039 bytes excess) Building topology for vector map vector_tty@assessment... Registering primitives... Killed Question about larger that it should be: What is this limitation? I guess is related to the specs of my laptop? Thanks, -- Best regards, Dr. Margherita DI LEO Scientific / technical project officer European Commission - DG JRC Institute for Environment and Sustainability (IES) Via Fermi, 2749 I-21027 Ispra (VA) - Italy - TP 261 Tel. +39 0332 78 3600 margherita.di-...@jrc.ec.europa.eu Disclaimer: The views expressed are purely those of the writer and may not in any circumstance be regarded as stating an official position of the European Commission. ___ 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] addons web page needs some corrections
Hello! Dear maintainers of GRASS site, could you correct an addons page, please? http://grass.osgeo.org/grass70/manuals/addons/ Descriptions for modules are moved down in some places. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] TIN to raster conversion module - v.tin.to.rast
Thanks for edits, Helmut! I will be more careful with Python scripts parameters. 2015-04-06 0:45 GMT+03:00 Helmut Kudrnovsky hel...@web.de: you use old g.remove interface, it should be: grass.run_command('g.remove', type='raster', name=tmp). [...] I'm not sure about g.region stuff, but I'm ready to upload scripts to grass-addons SVN. change to new g.remove syntax for g7, see: http://lists.osgeo.org/pipermail/grass-commit/2015-April/035894.html - best regards Helmut -- View this message in context: http://osgeo-org.1560.x6.nabble.com/TIN-to-raster-conversion-module-v-tin-to-rast-tp526p5200031.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 mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] TIN to raster conversion module - v.tin.to.rast
Hello all! I have corrected and adapted to recent GRASS versions (6.4.* and 7.0) the script of Antonio Alliegro tin.to.raster. These are 2 separate scripts for 6.4.* and 7.0: https://github.com/amuriy/GRASS-scripts/blob/master/Python/v.tin.to.rast_G64.py https://github.com/amuriy/GRASS-scripts/blob/master/Python/v.tin.to.rast_G70.py Could you test them, please? Maybe some errors or enhancements? Thanks in advance! ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] TIN to raster conversion module - v.tin.to.rast
Thanks for your answer, Anna :) I'm not sure about g.region stuff, but I'm ready to upload scripts to grass-addons SVN. 2015-04-05 21:33 GMT+03:00 Anna Petrášová kratocha...@gmail.com: Hello, On Sun, Apr 5, 2015 at 5:02 AM, Alexander Muriy amu...@gmail.com wrote: Hello all! I have corrected and adapted to recent GRASS versions (6.4.* and 7.0) the script of Antonio Alliegro tin.to.raster. These are 2 separate scripts for 6.4.* and 7.0: https://github.com/amuriy/GRASS-scripts/blob/master/Python/v.tin.to.rast_G64.py https://github.com/amuriy/GRASS-scripts/blob/master/Python/v.tin.to.rast_G70.py Could you test them, please? Maybe some errors or enhancements? Thanks in advance! thank you for your contribution. I tried the G7 version and works great, just fix please the clean up, you use old g.remove interface, it should be: grass.run_command('g.remove', type='raster', name=tmp). Also I would remove stderr = nuldev from all run_commands, not sure why they are there. Have you considered uploading the scripts in the official GRASS addons repository? It would be easier for people to install it then. And it could be maintained by other developers too. http://trac.osgeo.org/grass/wiki/HowToContribute#WriteaccesstotheGRASS-Addons-SVNrepository Cheers, Anna ___ 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] getting from DBF to SQLite
2013/2/14 Markus Neteler nete...@osgeo.org On Mon, Apr 9, 2012 at 1:27 PM, Alexander Muriy amu...@gmail.com wrote: Hi, Richard. Is it is possible to set a single SQLite database for the entire GRASS GIS database or each mapset requires its own sqlite db file? Use db.connect and then v.db.reconnect.all to reconnect from DBF. For my own needs I wrote small shell-script dbf2sqlite that makes one ... For the record: there is now this simplified procedure available: http://grasswiki.osgeo.org/wiki/Convert_all_GRASS_6_vector_maps_to_GRASS_7 (requires current GRASS 7 due to recent bugfixes) Markus Markus, thanks for the tip. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.denoise segfaults on a raster of 1166468401 cells
2012/12/21 grass-user-requ...@lists.osgeo.org Hi Nikos. In my experience, mdenoise utility works well with raster sizes no more than 2000x2000 (regardless of RAM size). That's why the only way to denoise the large rasters ? using of tiles and patching of final results. Of course! Yet, what about differently smoothed tile borders? Any additional hints? In my work, patching smoothed tiles (separate rasters) not cause any serious distortions in result raster. Anyway, it needs some more investigations :) --- Best regards, Alexander. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] i.landsat.trim -- script to trim Landsat images, make footprints and all-coverage mask
2012/12/21 Nikos Alexandris n...@nikosalexandris.net I am testing i.landsat.trim. Nice work! Thanks, Nikos! how do we interpret exactly the rast_buffer parameter? rast_buffer is the service parameter, which is used for several actions with the input raster bands. First of all, it identifies the width of raster buffer zone for making the band mask; also it controls the resolution of region during processing. It was gained experimentally with the best value near 300 meters for Landsat 7 images (sounds not very logical, I understand). Testing on a single L5T band, I get strange (to me) results with various buffer distances (e.g., 100, 150, 200, 5000, 1). Where does the buffer start? Does it start from the outer non-NULL pixels of the image for example? Yes, the buffer start on the border of the band mask, but with the significantly coarser region resolution. 1 Is it expected that some distance (e.g. 100) will result in an (almost) empty image? Can be expected almost anything, that's why the module needs some work to make it more logical and more usable. I already started rewriting it in Python, along the way I'll try to complement it and fix some functions. Best regards, Alexander. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.denoise segfaults on a raster of 1166468401 cells
2012/12/15 grass-user-requ...@lists.osgeo.org The following attempt segfaults on GRASS 6.4.3svn (nearly in the end of the process, if I am not wrong), running on Core i7 with 32GB of RAM. # mdenoise ASTER GDEM2 over Greece r.denoise in=aster_gdem2_ellas out=aster_gdem2_ellas_smoothed iterations=1 threshold=0.8 epsg=2100 The computational area comprises (Rows: 32401, Columns: 36001) 1166468401 cells. Is this too much? Would (re-)running this on GRASS 7.0.svn make any difference? I should better perform the process on a per (ASTER GDEM2) tile(s) basis, shouldn't I? Hi Nikos. In my experience, mdenoise utility works well with raster sizes no more than 2000x2000 (regardless of RAM size). That's why the only way to denoise the large rasters — using of tiles and patching of final results. --- Best regards, Alexander. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] reading AVHRR NDVI composites
I am trying to read the AVHRR 14 day NDVI composites downloaded from Earth Explorer. Each composite contains a .img and a .ddr file and gdal doesn't recognize it. I guess this is more of a gdal related question, but hopefully someone in the list has already done it. AVHRR rasters are in native binary format. Try module r.in.bin (in manual has an example for AVHRR). ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Removing spurious features on vector map
I'm tediously moving one line out of the way, deleting the second line and the point, then moving the first line back again. I don't know that I could specify only certain points to be removed programmatically. Rich, it's an editing hell... Maybe try to use v.select to select points that lie on boundaries and then remove them? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Removing spurious features on vector map
The attached image shows a vector map that has the attribute displayed along the boundaries of the various polygons. It should have that attribute only on the centroid. When I try to delete those spurious features I cannot distinguish them from the boundary lines in the wxpython GUI digitizing window. As I see, these features are zero areas with centroids. You can add polygons areas to attribute table with v.to.db and then remove polygons with the areas smaller than the chosen value (v.edit) or extract large polygons (v.extract). ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.in.ascii error: bad type reading source file
Running 6.5svn on Slackware and trying to import an ASCII text file with lat/lon coordinates. Here are the first two lines of the file: SC-100|41.42188333|-115.99418333 SC|42.30416667|-116.1060 When I run this command: v.in.ascii -n in='/home/rshepard/projects/nevada/queenstake/sites-data/biosites.csv' out=biosites format=standard col='name varchar(10), y double precision, x double precision' x=3 y=2 cat=0 --o Rich, try to import in point format and without -n flag: v.in.ascii in='/home/rshepard/projects/nevada/queenstake/sites-data/biosites.csv' out=biosites col='name varchar(10), y double precision, x double precision' x=3 y=2 cat=0 --o --- Alexander ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] orthorectify satellite images
Hi Josep. Except GRASS you may also try GDAL, OSSIM, ORFEO Toolbox (OTB), wxGIS to orthorectify the imagery. A few links (translated, originals in Russian): http://translate.google.ru/translate?hl=rusl=rutl=enu=http%3A%2F%2Fgis-lab.info%2Fqa%2Fortho-wxgis.html http://translate.google.ru/translate?hl=rusl=rutl=enu=http%3A%2F%2Fgis-lab.info%2Fqa%2Forbview3-ortho-gdal.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] i.landsat.trim -- script to trim Landsat images, make footprints and all-coverage mask
Hello, comrades. I wrote shell script i.landsat.trim, that trims the fringe from the borders of Landsat images, for each band separately or with the MASK where coverage exists for all bands. Optionally saves vector footprints of trimmed rasters and MASK. Works with Landsat 5, Landsat 7 (SLC-on). Link: https://raw.github.com/amuriy/GRASS-scripts/72f039073ff55b006b7aecbaa7870fac193dd9b3/i.landsat.trim Usage: i.landsat.trim [-maf] [input=string] [input_base=string] [input_prefix=string] output_prefix=string [rast_buffer=value] [gener_thresh=value] [--verbose] [--quiet] Flags: -m Trim raster(s) with the MASK where coverage exists for all bands -a Process all bands -f Save vector footprint(s) of trimmed raster bands or coverage MASK --v Verbose module output --q Quiet module output Parameters: input Name of input raster band(s) Example: L5170028_02820070521_B10 input_base Base name of input raster bands Example: L5170028_02820070521 input_prefix Prefix name of input raster bands Example: 'B.' for B.1, B.2, ... output_prefix Prefix for output raster maps Example: 'trim' generates B.1.trim, B.2.trim, ... rast_buffer Distance for raster buffering (in meters) default: 300 gener_thresh Threshold for generalizing of vector footprints or coverage MASK (in meters) default: 3000 Picture with example attached. It was made just for fun, but I hope it would be useful for somebody :) Best regards, Alexander ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Script in Grass - Mapcalc and MODIS
Hi Lucien. Does also someone has some experience with MODIS images in hdf format? Is it possible to work with it with GRASS? A few links: http://grass.osgeo.org/wiki/MODIS http://www.safe.iis.u-tokyo.ac.jp/pdf/MODIS_Image_Processing_using_GRASS_ver1.3.pdf ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to get centroid coordinates
Hi Margherita. I have a vector of polygons and i want to create a vector of points with the centroids of the polygons, with a connected table showing their coordinates. Which is the simplest way to do that on G7? In GRASS 6.4.*: -- extract centroids from polygon map v.extract in=poly out=centr type=centroid -- if needed convert centroids to points v.type in=centr out=centr_pts type=centroid,point -- add attribute table to points v.db.addtable centr_pts col=x double,y double,z double -- insert coords in table v.to.db TMP_centr_pts option=coor col=x,y,z I think in GRASS 7.0 it will be similar. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] creating two new maps from one existing raster
Hi Maris and all. r.split.line is already in my github repository ( https://raw.github.com/amuriy/GRASS-scripts/master/r.split.line) and on GRASS AddOns wiki page ( http://grass.osgeo.org/wiki/GRASS_AddOns#r.split.line). 2012/8/23 Maris Nartiss maris@gmail.com Hello Alexander, if You are not willing to make it into GRASS addon, then, please, submit it to issue tracking system so it doesn't get lost in ML. Best regards for all GRASS'ers, Alexander. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] creating two new maps from one existing raster
Hi, Janet and all. can anyone tell me how to make 2 rasters out of my 1 along a user specified line, or how to successfully digitize this? thank you, I wrote a small shell script r.split.line for splitting the raster map into parts with the line(s). Hope it will be useful. -- ./r.split.line -help Description: Split raster into parts with vector line(s). Lines MUST intersect the whole raster map. Keywords: raster,vector Usage: r.split.line raster=name line=name [--verbose] [--quiet] Flags: --v Verbose module output --q Quiet module output Parameters: raster Name of input raster map line Name of vector map with line(s) --- The script is attached. Spearfish example (raster elevation.10m) with animation in GIF: http://img708.imageshack.us/img708/7835/rsplitline.gif --- r.split.line.gz Description: GNU Zip compressed data ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] create polyines based on an attribute column
Hi Otto. is it possible to create polyines based on an attribute column? I found v.build.polylines using cat values, but that's not, what I need. Could you describe the whole task and your particular data in detail ? I think some scripting is required.. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.select usage
*Rich GRASS don't work with 3D-topology. If your drillholes map is 3D, try to make its 2D copy using v.to.3d module with the flag -r and then try to select drillholes in catchments. v.to.3d -r in=drillholes out=drillholes_2d v.select ain=drillholes_2d bin=catchments out=holesbybasin op=overlap is this not the appropriate module to find the solution to points-in-polygons? Try module v.distance. As written in man page: Point-in-polygon The option *dmax=0* is here important because otherwise for points not falling into any area, the category of the nearest area is recorded. For each point from vector map *pnt*, find the *area* from vector map *ar*in which the individual point falls, and write the related area categories to column *areacat* into the attribute table attached to vector map *pnt*: v.distance from=pnt to=ar dmax=0 upload=cat column=areacat ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.extrude not changing anything
*giannis I have a 2D vector map which is actually a single line (it shows a coastline). I try to give to this fixed height and make it a 3D map (the line to become a wall). I try v.extrude --overwrite input=COASTLINE@final output=coast_3d@finalheight=10 type=line but nothing changing at all (then I export this in .vtk and i check it out in 3D software). Try also module v.to.3d: v.to.3d input=COASTLINE output=coast_3d height=10 type=line ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to split a line using points?
I have points over a line and i need to split the line where the points are located. The points aren't at the same distance to each other. Any idea how to do that? If your points lie exactly on the line, you may do so (in Linux with bash) : v.out.ascii format=point in=pts --q | cut -d'|' -f1,2 | tr '|' ',' | while read COOR; do v.edit map=line tool=break coords=$COOR; done ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to install v.strahler add-on under Windows?
Thanks, Markus and Martin! Now v.strahler is installed under 6.4.* and works fine. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] How to install v.strahler add-on under Windows?
Hi, comrades. Suggest me please, how to compile the module v.strahler under Windows (under Linux I do it easily and with the pleasure)? There are precompiled add-ons here (http://wingrass.fsv.cvut.cz/), but all archives with this module are empty (for example, http://wingrass.fsv.cvut.cz/grass64/addons/grass-6.4.3svn/v.strahler.zip) ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] round floating map
Hi, Salvatore. You can also try the special module r.quanthttp://grass.gis-lab.info/grass64/manuals/html64_user/r.quant.html( http://grass.gis-lab.info/grass64/manuals/html64_user/r.quant.html). I'm trying to round a floating map. I revived an error: Overflow occured in the calculation. Why? I used the command r.mapcalc R_round=round(R_seJ_ha_filled) My goal is to obtain an integer map. Is there another way? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: grass-user Digest, Vol 73, Issue 25
Il 11/05/2012 13:05, Paolo Cavallini ha scritto: Il 11/05/2012 13:03, Alexander Muriy ha scritto: I wrote shell script to find average line(s) of input vector map. BTW: why not adding it to GRASS trunk? Thanks Paolo :) I'll add it soon at add-ons wiki page, and then can add it to AddOns SVN repository. Short manual also needed to be made :) ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: r.finder
I recently published a grass script for focussing surveys. The name is r.finder, hoping it can be of help for somebody... Augusto, great news :) I'll test it in my geomorphological work surely. You would not want to add it to GRASS-Addons wiki page ( http://grass.osgeo.org/wiki/GRASS_AddOns)? One note: in GRASS shell scripts it's usually better to use: -- module g.message instead of standard echo -- special temp files (made with g.tempfile): # create temporary file TMP1=$(g.tempfile pid=$$) if [ $? -ne 0 ] || [ -z $TMP1 ] ; then g.message -e ERROR: Unable to create temporary files. exit 1 fi -- cleanup function and some actions for the case of user exit cleanup() { g.mremove -f vect=MY_TEMP_VECTOR --q rm -f $TMP1* } ## what to do in case of user break: exitprocedure() { echo User break! cleanup exit 1 } ## shell check for user break (signal list: trap -l) trap exitprocedure 2 3 15 --- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.triangle -- module to produce TIN with breaklines using Triangle utility
Hi all. I wrote a shell-script module v.triangle that is front-end for the famous Triangle utility (http://www.cs.cmu.edu/~quake/triangle.html) of J.R. Shewchuk. I'm planning to add it to GRASS-Addons wiki page :) Triangle makes exact Delaunay triangulations, constrained Delaunay triangulations, conforming Delaunay triangulations and high-quality triangular meshes. On GIS terminology, it produces 2D TIN, optionally with breaklines. Link: https://raw.github.com/amuriy/GRASS-scripts/master/v.triangle P.S. It generally works but centroids in final TIN have 0 heights. Now I think about some hack to make the right centroids in all areas of TIN. - Comparison of the one resulting TIN with the TIN produced in ArcView 3D Analyst see in linked PDF ( http://dl.dropbox.com/u/42016908/ArcView_GRASS_TIN_breaklines.pdf) - ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: getting from DBF to SQLite
Hi, Richard. Is it is possible to set a single SQLite database for the entire GRASS GIS database or each mapset requires its own sqlite db file? Use db.connect and then v.db.reconnect.all to reconnect from DBF. For my own needs I wrote small shell-script dbf2sqlite that makes one SQLite DB in mapset (with the mapset name) and copies all tables from DBF to SQLite ( https://raw.github.com/amuriy/GRASS-scripts/f5157c8d2c8a5b3a9149cf207b62ef088a51d1ef/dbf2sqlite ) Best regards, Alexander. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: r.surf.nnbathy and nnbathy installation
Hello Danny. does nobody has an idea what I have done wrong with the implementation of r.surf.nnbathy? :-( I work with r.surf.nnbathy constantly and haven't seen any serious troubles with it. Did you install it according to instruction from the manual ( http://www.sieczka.org/prog/grass/r.surf.nnbathy.html)? -- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: r.surf.nnbathy and nnbathy installation
how did you link the script in your GRASS? It seems that GRASS doesn't start the script. Danny, just put directory with r.surf.nnbathy in your system PATH (or set GRASS_ADDONS_PATH). Of course, /usr/local/bin also should be in PATH. It should be work in any version of GRASS 6.4.* - ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: i.vi.mpi
Hello Ripsime. This is my Makefile MODULE_TOPDIR = /usr/local/grass-6.4.2svn/ CC=mpicc --static PGM = i.vi.mpi LIBES = $(GISLIB) $(GMATHLIB) DEPENDENCIES = $(GISDEP) $(GMATHDEP) include $(MODULE_TOPDIR)/include/Make/Module.make ifneq ($(USE_LARGEFILES),) EXTRA_CFLAGS = -D_FILE_OFFSET_BITS=64 endif default: cmd CC=mpicc -L/usr/local/grass-6.4.2svn/lib/ /usr/local/grass-6.4.2svn/lib/ libgrass_gis.6.4.2svn.so When I do a make I am getting this: make[2]: *** [/usr/local/grass-6.4.2svn/man/man1/i.vi.mpi.1] Error 126 make[2]: Leaving directory `/localuser/ripsime/mpitest' make[1]: *** [mancmd] Error 2 make[1]: Leaving directory `/localuser/ripsime/mpitest' make: *** [cmd] Error 2 You ignored these compilation errors and want the working module ... Why not to use the default Makefile from i.vi.mpi add-on directory? I used command: make MODULE_TOPDIR=/home/user/bin/grass643svn_bin/grass-6.4.3svn/ and i.vi.mpi compiled successfully. So try to recompile module and than run you script again. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.ldm -- add-on module to compute Linear Directional Mean
Dear GRASS'ers! I wrote the shell-script v.ldm that computes Linear Directional Mean of the vector map (as explaned there http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/How_Linear_Directional_Mean_works/005p001r00/) and displays LDM graphics on the X-monitor, optionally save it to vector line and update attribute table with LDM parameters. Links: https://raw.github.com/amuriy/GRASS-scripts/master/v.ldm http://grass.osgeo.org/wiki/GRASS_AddOns#v.ldm I hope it would be useful for somebody :) ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] d.zoom.keys -- shell script for navigation in GRASS monitors from keyboard
Dear GRASS'ers :) Let me introduce (a bit hacky) shell script ( https://raw.github.com/amuriy/GRASS-scripts/master/d.zoom.keys) for navigating in classic X-monitors with the typical keys (arrows and - / +). It works in two modes: from terminal or directly from X-monitor (requires xev and xdotool utilities). Note: tested only with GRASS 6.4.2svn and only in Linux (Ubuntu 10.04). d.zoom.keys -help Description: Allows to change the current geographic region settings interactively, with a keyboard. Use arrow keys to navigate, - to make futher, + (without Shift) to make closer. Keywords: display, zoom Usage: d.zoom.keys [-txr] [mon=name] [magn=number] [--verbose] [--quiet] Flags: -t Use navigation from terminal -x Use navigation in X-monitor (requires xev and xdotool) -r Decrease region resolution during navigation --v Verbose module output --q Quiet module output Parameters: mon Name of graphics monitor to select (default is currently selected monitor). magn Magnification: default is 1 (1/10 of current region). magn 1 gives positive magnification, magn 1 gives negative one. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: grass and palsar processing
Also there is PDF http://gis-lab.info/forum/download/file.php?id=2370, copied from http://www.safe.iis.u-tokyo.ac.jp/, about SAR GRASS. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: GRASS GIS development visualization
To download from Youtube: Firefox with Flash Got plugin or something similar. Also in Linux: youtube-dlhttp://www.arrakis.es/%7Erggi3/youtube-dl/, clive http://clive.sourceforge.net/. Message: 1 Date: Mon, 14 Nov 2011 19:32:10 -0700 From: Michael Barton michael.bar...@asu.edu Subject: [GRASS-user] video To: grass-user@lists.osgeo.org grass-user@lists.osgeo.org Message-ID: 6f076750-0140-4b9d-96cb-b84cab697...@asu.edu Content-Type: text/plain; charset=us-ascii Markus, Is there a way to download the low and hires versions of this video? Michael C. Michael Barton Director, Center for Social Dynamics Complexity Professor of Anthropology, School of Human Evolution Social Change Arizona State University ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: How to generate a bigger PNG file
Luisa, you can use module d.out.file ( http://grass.gis-lab.info/grass64/manuals/html64_user/d.out.file.html) with the *resolution* parameter (double size or quadruple of the image). Date: Thu, 3 Nov 2011 12:16:58 + From: Luisa Pe?a luisapena1...@gmail.com Subject: [GRASS-user] How to generate a bigger PNG file To: GRASS user list grass-user@lists.osgeo.org Message-ID: cah+a188ytff79mfg0ya6f9l9w9ockmxsj6bn1yxxpuguoyr...@mail.gmail.com Content-Type: text/plain; charset=iso-8859-1 Greetings I'm using r.out.png to generate PNG files from a small patch (like 6x15 pixels) but I'm obtaining a really small PNG file. Sicne I want to display it a little bit bigger in a website I need to create a bigger (in size) PNG. What can I do to do this? One idea was to rescale it in a image processing software but I get a low quality image Is there any solution by using grass? Thanks Luisa -- next part -- An HTML attachment was scrubbed... URL: http://lists.osgeo.org/pipermail/grass-user/attachments/2003/6cef341b/attachment-0001.html ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: ps.map arrow head
Vincent, Hamish and all. Maybe it would be easier to do PS with ps.map, convert it to pdf (for example with ps2pdf), then open in vector editor like Inkscape and add arrowheads manually. On this topic I made small shell-script d.arrowheadhttps://raw.github.com/amuriy/GRASS-scripts/master/d.arrowhead that draw arrowheads of vector lines on the graphic monitor using d.graph (in fact, with any custom symbol) and optionally save graph file (to use with d.redraw / d.save etc.) and ps.map file (to use with ps.map read instruction). This command makes red arrows with the custom symbol arrows/arrowhead_1 (created by the script), and saves graph and psmap files d.arrowhead -s test symbol=arrows/arrowhead_1 size=20 width=2 col=red graph=test.graph psmap=test.psmap - Output PNG with d.out.file and PDF with ps.map+ps2pdf: attachment: arrows_test.png arrows_test.pdf Description: Adobe PDF document ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: ps.map arrow head
Hi, Vincent, Hamish and all. Maybe it would be easier to do PS with ps.map, convert it to pdf (for example with ps2pdf), then open in vector editor like Inkscape and add arrowheads manually. On this topic I made small shell-script d.arrowheadhttps://raw.github.com/amuriy/GRASS-scripts/master/d.arrowhead that draw arrowheads of vector lines on the graphic monitor using d.graph (in fact, with any custom symbol) and optionally save graph file (to use with d.redraw / d.save etc.) and ps.map file (to use with ps.map read instruction). This command makes red arrows with the custom symbol arrows/arrowhead_1 (created by the script), and saves graph and psmap files d.arrowhead -s test symbol=arrows/arrowhead_1 size=20 width=2 col=red graph=test.graph psmap=test.psmap - Output PNG with d.out.file and PDF with ps.map+ps2pdf: attachment: arrows_test.png arrows_test.pdf Description: Adobe PDF document ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Add DB entries to vector line (after v.split)
Date: Tue, 25 Oct 2011 12:37:36 -0700 (PDT) From: wittich vwitt...@gmail.com Subject: [GRASS-user] Re: Add DB entries to vector line (after v.split) To: grass-user@lists.osgeo.org Message-ID: 1319571456756-6930154.p...@n2.nabble.com Content-Type: text/plain; charset=us-ascii v.info lines_split ++ | Layer: lines_split | Map format: native | | Type of Map: vector (level: 2) | Number of points: 0 Number of areas: 0 | Number of lines:598 Number of islands:0 | Number of boundaries: 0 Number of faces: 0 | Number of centroids:0 Number of kernels:0 | | Map is 3D: No | Number of dblinks: 1 | | Digitization threshold: 0 | Comments: ++ v.category in=lines_split opt=report Layer/table: 1/lines_split type countminmax point 0 0 0 line 598 1 3 boundary 0 0 0 centroid 0 0 0 area 0 0 0 all 598 1 3 Vector map looks normal. Maybe try to run one more time: v.category in=lines_split out=lines_split_nocats type=line opt=del v.category in=lines_split_nocats out=lines_split_newcats type=line opt=add ?? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Add DB entries to vector line (after v.split)
Message: 5 Date: Mon, 24 Oct 2011 10:34:10 -0700 (PDT) From: wittich vwitt...@gmail.com Subject: [GRASS-user] Re: Add DB entries to vector line (after v.split) To: grass-user@lists.osgeo.org Message-ID: 1319477650371-6925911.p...@n2.nabble.com Content-Type: text/plain; charset=us-ascii thx for your answer. I tried the v.category function but it seams like it doesn't work. And in case that it would work, I would lost the original categories the lines belongs to, don't I? It looks strange. Can you give the output of commands v.info and v.category in=vector opt=report for your segments vector? lines_split_nocats and lines_split_newcats give me the following v.db.select ... cat 1 2 3 It seams like the Option 'del' doesn't work... I tried to delete the table as well, but if I built it again, I've the old categories. But anyway my workaround works. Deletion of the table is unnecessary, because cats in table are depended on categories of the vector objects. So all work is around categories. ps: I'm gonna try your shell-script (v.to.equidist), as soon I'll have installed the GRASS-Addons. Right now it didn't find the grass/gis.h while he compiled it... I hope you're not going to compile a script? :) ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Add DB entries to vector line (after v.split)
Hi, wittich. -- firstly you should delete old categories from your segments' lines (they relate to the initial 3 lines): v.category in=lines_split out=lines_split_nocats opt=del -- then add new cats to segmented lines: v.category in=lines_split_nocats out=lines_split_newcats opt=add *As a promotion*: there is a small shell-script v.to.equidisthttp://grass.osgeo.org/wiki/Addons#v.to.equidist on GRASS-Addons page that generates vector points or line segments along a given vector line(s) with the equal distances (uses v.segment) and populate tables with cats, length etc. Date: Sun, 23 Oct 2011 10:26:06 -0700 (PDT) From: wittich vwitt...@gmail.com Subject: [GRASS-user] Add DB entries to vector line (after v.split) To: grass-user@lists.osgeo.org Message-ID: 1319390766722-6922423.p...@n2.nabble.com Content-Type: text/plain; charset=us-ascii Hi there, I try to find a way to give every line segment in my vector file a db entry. I've create 3 lines and split them into 50m segments. Now I want to do several calculations with it. Therefore I need for each segment a row in the db. *How can I create a entry for each segment?* v.info split_lines Number of lines: 598 v.db.select split_lines cat 1 2 3 I would appreciate if someone can explain me what I do wrong, I tried already everything I could found without success (eg. v.dissolve, v.category, v.centroids, and others) regards Valentin -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Add-DB-entries-to-vector-line-after-v-split-tp6922423p6922423.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] ASTER GDEM Version 2 --- filtering and smoothing
Hello, dear GRASS users. A few words about new ASTER GDEM version 2https://lpdaac.usgs.gov/about/news_archive/monday_october_17_2011 . As was explored, this DEM could be also rather well filtered and smoothed with the Sun's denoising algorithmhttp://personalpages.manchester.ac.uk/staff/neil.mitchell/mdenoise/(using GDAL and free / open source program mdenoise). Also there is GRASS add-on r.denoise http://grass.osgeo.org/wiki/GRASS_AddOns#r.denoise. Experiments showed that the best smoothing of ASTER GDEM 2 is reached with such parameters of mdenoise: threshold = 0.8 iterations = 10-20 Also filtering with r.neighbors by average method and window size =5 is quite usefull to remove some noise from DEM. --- With the best regards from Russia, Alexander Muriy. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Can't edit GRASS-Wiki page normally
Hello all. I would like to make some minor edits in GRASS_AddOns wiki page (r.denoize links update). I'm a registered user. After more than twenty unsuccessful attempts to enter CAPTCHA (all different) I decided to write here. What happened with wiki? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] GDAL/OGR user docs in one PDF
Sorry for some offtopic, but almost all GRASS'ers use GDAL :) NOTE: now all GDAL/OGR user docs are also available as one PDFhttp://gis-lab.info/docs/gdal/gdal_ogr_user_docs.pdfand as a single HTML-page http://gis-lab.info/docs/gdal/gdal_ogr_user_docs.html from GIS-Lab (russian GIS community). ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Cannot create dbf database
Maybe need to check the connection one more time? $ eval $(g.gisenv) $ db.connect driver=dbf database='$GISDBASE/$LOCATION_NAME/$MAPSET/dbf/' $ db.connect -p $ db.tables -p And then: $ v.db.connect -p contours ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Cannot create dbf database
*Craig You used command v.in.ogr -e dsn=contours lo29.shp out=contours2 with the -e flag , this flag is for creating a new location with the imported data. But you don't specify any location.. And may be better to use single quotes with the filenames like 'contours lo29.shp'. --- Alexander ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: grass-user Digest, Vol 65, Issue 49
*Jim Trivially: did you do g.region rast=na_dem_sub_20 before raster computations? And r.mapcalc expression is without any quotes like r.mapcalc na_dem_sub_20 = na_dem_input - 20 or r.mapcalc 'na_dem_sub_20 = na_dem_input - 20' Alexander ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[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
[GRASS-user] [v.to.equidist] -- segments generator for v.segment
Dear GRASS users! I want to introduce you shell-script v.to.equidisthttps://raw.github.com/amuriy/GRASS-scripts/master/v.to.equidist which generates points or line segments along vector lines using module v.segment. It's already posted at GRASS_AddOns wiki page. It's an attempt to have the functionality of perl-scripts v.lmeasure / v.revlmeasure http://grass.osgeo.org/wiki/GRASS_AddOns#v.lmeasure, but in more consistent form. It can do: -- direct and reverse segmenting -- broken line segments instead of continuous -- create attribute table for output vector and populate it with the values of line cats, distances etc. Thanks for your attention. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: v.in.ogr slowness with GRASS 6.4.1 on Linux
If I try to import a shapefile using GRASS 6.4.1 using v.in.ogr, it takes more than an order of magnitude longer than previous versions of GRASS that I built and used. Maybe you're trying to import polygons with a topology cleaning? In that case use -c flag in v.in.ogr (and/or -t flag). What version of GDAL/OGR? (maybe problem in GDAL..) ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: v.in.ogr slowness with GRASS 6.4.1 on Linux
* Did v.in.ogr change from, say, GRASS 6.3 and 6.4.1? In what ways? * Honestly, I don't know about 6.3 and 6.4.1 differences (probably only in GDAL version), but GRASS 7 already ( http://grass.osgeo.org/wiki/GRASS_7_ideas_collection#Vector) is more speed up than 6.4*. Try to search about v.in.ogr in this mail list [http://www.google.com/search?q=v.in.ogras_qdr=allsa=Google+Searchdomains=lists.osgeo.orgsitesearch=lists.osgeo.org%2Fpipermail%2Fgrass-user%2F] Also you may try to install a more ancient version of GRASS from RedHat repositories. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: r.param.scale curvatures formulas
2011/5/31, grass-user-requ...@lists.osgeo.org grass-user-requ...@lists.osgeo.org: Sorry to resurrect a reasonably old topic (following my mail), but I have two questions on the r.param.scale module: 1 - I share Jasiewicz's concerns, the formula do not match the litterature I got [0] Formulas for this module were investigated by Jo Wood in his work http://www.soi.city.ac.uk/~jwo/phd/ So there can be some special cases.. I know other modules, such as v.surf.rst, offer to compute the tangential curvature, but I'd be glad to generate the different curvatures I need using only one command. Also as an idea: it's possible to write custom script (shell, python), which would be compute different types of curvatures by calling different modules (r.param.scale, r.slope.aspect,v.surf.rst) ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: d.frontline Shell script to display frontlines (Alexander Muriy)
Hi all. I corrected my module d.frontline, now it works (I hope) a little better. It's possible to use any user defined symbol (also front/arrow_int and front/arrow_ext added to front/... symbols). Also now it's possible to get output file with ps.map instructions (to use later with ps.map read command). Alexander. d.frontline.gz Description: GNU Zip compressed data ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: d.frontline Shell script to display frontlines (Alexander Mury)
Hello, Hamish. Thank you for your answer. I've looked through it and it is quite well done. Could we add it to the grass-addons repository? It would be great, but first I want to correct it a bit. * consider moving the temp file creation until after all the parsing checks are done, so that they are only generated when we know the module is going to run and empty files are not left behind. * always quote filename $VARIABLES, or else the script will break on MS Windows where path names commonly have spaces in them. Ok, understand. * did you try the geology/strike* symbols before making your own? http://grass.osgeo.org/wiki/IconSymbols#New were they deficient for some reason? (too long bases for the task, or ..?) [no problem, just like feedback if they could be better] Yes, it was my first attempt, but all geology/... symbols have (you're right) too long bases for the task. So I decide to make front symbols == geology symbols without bases. * consider advertising it in the geology + atmospheric science field of study sections in the GRASS wiki: http://grass.osgeo.org/wiki/Applications Ok. * besides left,right, consider a both side of line option. Ok, it's a good idea. the problem of line styles is still not completely solved, hopefully we can do something with that in grass7, but this has been a missing functionality for a long time- so thanks! In GRASS 7 without good-old monitors I can't realize what to do with all of this graph modules (some of them are really useful..) I think it's possible to implement some other line styles in 6.4.* version with the same method as in d.frontline. ps- I just noticed the collection of new symbols by Martin Povolny. lots of them, great. :) http://grass.osgeo.org/wiki/IconSymbols#More_simple_symbols Thanks, --- So, I will correct my module and will send it again. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user