Re: [GRASS-user] Add DB entries to vector line (after v.split)

2015-05-08 Thread Alexander Muriy
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

2015-04-24 Thread Alexander Muriy
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

2015-04-07 Thread Alexander Muriy
​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

2015-04-05 Thread Alexander Muriy
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

2015-04-05 Thread Alexander Muriy
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

2015-04-05 Thread Alexander Muriy
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-02-13 Thread Alexander Muriy
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-20 Thread Alexander Muriy
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-20 Thread Alexander Muriy
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-19 Thread Alexander Muriy
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

2012-11-28 Thread Alexander Muriy
 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

2012-10-10 Thread Alexander Muriy

  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

2012-10-09 Thread Alexander Muriy
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

2012-10-06 Thread Alexander Muriy
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

2012-09-14 Thread Alexander Muriy
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

2012-09-03 Thread Alexander Muriy
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

2012-08-29 Thread Alexander Muriy
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

2012-08-24 Thread Alexander Muriy
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

2012-08-23 Thread Alexander Muriy
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

2012-08-22 Thread Alexander Muriy
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

2012-07-13 Thread Alexander Muriy
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

2012-06-28 Thread Alexander Muriy
*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

2012-06-28 Thread Alexander Muriy
*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?

2012-06-14 Thread Alexander Muriy



 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?

2012-06-12 Thread Alexander Muriy
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?

2012-06-10 Thread Alexander Muriy
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

2012-06-07 Thread Alexander Muriy
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

2012-05-11 Thread Alexander Muriy
 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

2012-05-11 Thread Alexander Muriy
   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

2012-04-20 Thread Alexander Muriy
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

2012-04-09 Thread Alexander Muriy
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

2012-03-26 Thread Alexander Muriy
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

2012-03-26 Thread Alexander Muriy

 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

2012-03-15 Thread Alexander Muriy
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

2012-03-08 Thread Alexander Muriy
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

2012-02-13 Thread Alexander Muriy
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

2012-02-10 Thread Alexander Muriy
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

2011-11-15 Thread Alexander Muriy
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

2011-11-03 Thread Alexander Muriy
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

2011-10-30 Thread Alexander Muriy
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

2011-10-30 Thread Alexander Muriy
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)

2011-10-26 Thread Alexander Muriy
 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)

2011-10-25 Thread Alexander Muriy
 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)

2011-10-24 Thread Alexander Muriy
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

2011-10-24 Thread Alexander Muriy
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

2011-10-23 Thread Alexander Muriy
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

2011-10-01 Thread Alexander Muriy
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

2011-09-26 Thread Alexander Muriy
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

2011-09-26 Thread Alexander Muriy
*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

2011-09-26 Thread Alexander Muriy
*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)

2011-09-21 Thread 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

2011-09-16 Thread Alexander Muriy
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

2011-06-05 Thread Alexander Muriy
 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

2011-06-05 Thread Alexander Muriy
* 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-05-31 Thread Alexander Muriy
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)

2011-05-27 Thread 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)

2011-05-25 Thread Alexander Muriy
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