[GRASS-user] Problems w/ v.rectify for 3D Cross-Section

[Phillip Allen]

Hi All,

It has been a few years since I had the time to work with Grass and I must
say the new 7.1 version looks interesting.

I have a vector geological cross-section digitized in the horizontal  with
all Zs=0 meters and I would like to properly Geo-reference it in 3D so it
lines up with my 3D bore-hole data that I am plotting. It is not to scale
so it needs to be scaled, rotated vertically, rotated horizontally,
translated horizontally/vertically.  It is in arbitrary UTM coordinates so
it is in the same location/coordinates.  So the relatively new module
v.rectify seems appropriate.

I have created a ground control coordinate text table but when I run
v.rectify all I get is the RMS table and no new-properly geo referenced.

Is the v.rectify command working properly for 3d polygons, polylines or
points?  My section needs be perfectly vertical so I can imagine that the
vertical polygons will break topology and create problems.  But I can live
with 3D polylines.

Below is the command I used.

Thanks for any help.

Phil

v.rectify -3 -o -r -b --overwrite input=xsec_d_geo1vline_temp3d at PERMANENT
output=xsec_d_geo1vline_gref points=C:\data\company\geology\xsect_d_gcp.txt
index|fwd_dx|fwd_dy|fwd_dz|back_dx|back_dy|back_dz|fwd_RMS|back_RMS
1|2.727261|2.636065|0.347856|8.196966|0.824918|0.428253|3.8089|8.2495
2|0.273510|0.056618|3.913415|0.326901|8.559731|0.518620|3.9234|8.5817
3|0.868377|1.496921|1.443858|3.656405|3.143833|0.838842|2.2538|4.8946
4|0.926665|1.577732|1.121287|3.867425|2.418439|0.865953|2.1460|4.6428
5|1.852593|2.899231|0.228845|7.361002|0.561665|1.352944|3.4482|7.5053
6|3.058322|2.404024|0.771569|8.361038|1.610876|1.298724|3.9659|8.6133
6|1.617788|1.845098|1.304472|5.294956|2.853244|0.883889|3.3467|7.2770
(Tue Oct 20 08:37:10 2015) Command finished (0
sec)


xsect_d_gcp.txt File:
527982.2 3828858.8 0 529227.0102 3834043.5563 1219.2 1
527988.7 3826198.5 0 529227.0102 3834043.5563 0 1
530199 3828873.2 0 529925.888 3833293.767 1219.2 1
530203.1 3828077.2 0 529925.888 3833293.767 853.44 1
529413.5 3828866.3 0 529676.798 3833562.382 1219.2 1
529406.4 3827276.8 0 529676.798 3833562.382 487.68 1
527982.2 3828858.8 200 529085.589 3833902.135 1219.2 0
530199 3828873.2 200 529784.467 3833152.346 1219.2 0

I added the 2 last gcp off the section line to make 3d
triangles/tetrahedrons but this still didn't help so I turned them off in
this example.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20151020/b9b94818/attachment.html>