complement of v.prune?

Susan Huse sue at ced.berkeley.edu
Thu Apr 7 19:42:42 EDT 1994


> From grass-lists-owner at max.cecer.army.mil Thu Apr  7 16:27:52 1994
> Date: Fri, 8 Apr 94 08:58:55 EST
> From: simon at artemis.earth.monash.edu.au (Simon Cox)
> Sender: grass-lists-owner at max.cecer.army.mil
> Reply-To: grassu-list at max.cecer.army.mil
> To: grassu-list at amber.cecer.army.mil
> Subject: complement of v.prune?
> Content-Length: 1195
> 
> No responses first time around, so I'll ask again!
> -----
> 
> 
> Is there a way to do the opposite of v.prune - ie ADD vertices
> along lines to keep them to a set spacing?
> 
> I suspect this
> could be done by exporting the file using v.to.sites using the
> -i option, and then re-importing the sites file into a vector
> (I do not look forward to editting the ascii file, however!) -
> except that the -i flag in v.to.sites doesn't work for me!
> (I've filed a bug report).
> 
> Thanks    Simon Cox
> 
> 
> 
> ----
> ______________________________________________________________________________
> from May 1st (or thereabouts)        __  L     Dr Simon Cox                   
> CSIRO Div. Exploration & Mining   ,~'  L_|\    VIEPS Dept of Earth Sciences,  
> 39 Fairway, PO Box 437,        ,-'         \   Monash University, Clayton     
> Nedlands, WA  6009,  Australia (   soon     \  Vic 3168,  Australia           
>       Phone +61 9 389 8421     X<~~~~~~~~~  /         Phone +61 3 905 5762    
>       Fax   +61 9 389 1906      L,~'   "\_x/          Fax   +61 3 905 5062    
>     s.cox at dem.csiro.au                    u  simon at artemis.earth.monash.edu.au
> ______________________________________________________________________________
> 

The only way I know would be to mess with the ascii, sorry. Use awk.
It may be a pain, but if v.to.sites -i doesn't work I don't know of a
better way.  

You could do it as a sites file or as a dig_ascii file.


If sqrt ( (x2 - x1)^2 + (y2-y1)^2 ) > threshold then newE=0.5(x2+x1);
newN=0.5(y2+y1) printf("%s %s\n%s %s\n", newE, newN, x2, y2) >
"filename".   Or something to that effect.  You'd also need a 
test loop to be sure that the distance between point2 and newE/N is
less than your threshold.

-sue




More information about the grass-user mailing list