[GRASS-user] Strategies to Import Geological Cross Sections

Phillip Allen phillipjallen830 at gmail.com
Tue Jan 12 14:37:49 PST 2016


Hi all,
I have digitized a few geological cross-sections each consisting of a 1
polygon & 1 polyline layer (representing lithology & faults respectively).
The layers were digitized in plan view with arbitrary coordinates in meters
and imported into PostgreSQL/PostGIS. I then created views using
ST_Translate(ST_RotateZ(ST_RotateX(ST_Translate(geom, ))) to correctly plot
the cross-section polygons & polylines in 3D space. One day I need to learn
how to multiply matrices so I can do it all in a single step!

I just tried to import them into GRASS 7.0.3RC1 using v.in.ogr and got the
following error message:

ERROR 1: Reading EWKB with 4-dimensional coordinates (XYZM) is not supported

Can GRASS not import PostGIS 3D data? Is it the vertical Polygon with zero
area topology problem?

What other options do I have to get 3D-cross-section data into GRASS?

I tried to use the ST_Extrude (to avoid topology problems) command in
PostGIS and Grass's v.in.ogr didn't like that either. Is there another way
to get 3D objects into GRASS?

I want to use v.hull on the vertical cross-sectional data and symbolize
some geochemical drill data in Grass's 3D viewer.

Regards,

Phil
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20160112/cd7f2611/attachment.html>


More information about the grass-user mailing list