<div dir="ltr">One solution is to export, pipe through invproj, and re-import:<div><br></div><div><br></div><div><div>v.out.ascii Z layer=-1 -c format=standard output=Zproj --o</div><div><br></div><div>rm Zrot</div><div>regex="[-]\d*\.\d*\s*\d*\.\d*"</div><div>while read line; do</div><div> if [[ $line =~ $regex ]]; then<br></div><div> echo $line | invproj -f "%f" +proj=ob_tran +o_proj=latlon +o_lon_p=-200 +o_lat_p=18 +lon_0=180 -m 57.295779506 >> Zrot</div><div> else</div><div> echo $line >> Zrot</div><div> fi</div><div>done < Zproj</div></div><div><br></div><div> -k.</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 13, 2018 at 3:06 PM, Ken Mankoff <span dir="ltr"><<a href="mailto:mankoff@gmail.com" target="_blank">mankoff@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Tue, Mar 13, 2018 at 1:05 PM, Ken Mankoff <span dir="ltr"><<a href="mailto:mankoff@gmail.com" target="_blank">mankoff@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 13, 2018 at 12:21 PM, Helmut Kudrnovsky <span dir="ltr"><<a href="mailto:hellik@web.de" target="_blank">hellik@web.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Ken Mankoff wrote<span><br><span>> Is there a way to take an arbitrary vector from<br>
> GRASS and apply an arbitrary transformation to each element? v.edit has a<br>
> move feature, but this is not what I want.<br>
><br>
<br>
</span></span><span>v.transform - Performs an affine transformation (shift, scale and rotate) on<br>
vector map.<br>
<a href="https://grass.osgeo.org/grass74/manuals/v.transform.html" rel="noreferrer" target="_blank">https://grass.osgeo.org/grass7<wbr>4/manuals/v.transform.html</a><br>
<br>
not sure it is what you're looking for.<br></span></blockquote><div><br></div><div>Yes that would do it just fine. Unfortunately I'm not sure how to get the transformation matrix that invproj uses. I'll try to figure this out...</div><span class="m_1664017099963134802HOEnZb"><font color="#888888"><div><br></div></font></span></div></div></div></blockquote><div><br></div></span><div>Actually this does not work because I'm in EPSG:4326 (lon,lat) location and there are scaling issues. I think I need to access proj directly for all vertices...</div><span class="HOEnZb"><font color="#888888"><div><br></div><div> -k. </div></font></span></div></div></div>
</blockquote></div><br></div>