<div dir="ltr"><div><div><div><div><div><div><div>Markus,<br><br></div>Thanks for shedding some light on this problem.<br><br></div>Yesterday I tried to work around it by extending my sql view with:<br><br>ST_Boundary(ST_Force3d(ST_Force2d(geom)))::geometry(MultiLineStringZ,32611)<br><br>so it would eliminate the M values (also I had a few z errors in the cross-section) and convert the polygons to polyLines.  I then exported them to shapefiles.  After the closed polyline-shapefiles were imported into Grass I then converted them to faces using:<br>v.type input=xsec_c_lithL output=xsec_c_lithL_face from_type=line to_type=face<br><br></div>This seemed to work so I would have an enclosed area to color.  But now I cannot seem to color the cross-section faces?  It is just black.  I have a text column with rgb values (255:255:255 format) but when I use it for "use color for thematic mapping"  it keeps reporting <br>DBMI-SQLite driver error: error in sqlite3_prepare(): unrecognized token error "'255:255:115"<br><br></div>What is odd is that I was able to color the polylines....I must admit that I have never really had much luck thematically mapping vectors in Grass.<br><br></div>So now the question is how to color the thematically mapped cross-section Faces?<br></div><div>I am using the latest version 7.0.3RC1.<br><br></div>regards,<br></div>Phil<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 14, 2016 at 6:55 AM, Markus Neteler <span dir="ltr"><<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<span class=""><br>
On Tue, Jan 12, 2016 at 11:37 PM, Phillip Allen<br>
<<a href="mailto:phillipjallen830@gmail.com">phillipjallen830@gmail.com</a>> wrote:<br>
> Hi all,<br>
><br>
> I have digitized a few geological cross-sections each consisting of a 1<br>
> polygon & 1 polyline layer (representing lithology & faults respectively).<br>
> The layers were digitized in plan view with arbitrary coordinates in meters<br>
> and imported into PostgreSQL/PostGIS. I then created views using<br>
> ST_Translate(ST_RotateZ(ST_RotateX(ST_Translate(geom, ))) to correctly plot<br>
> the cross-section polygons & polylines in 3D space. One day I need to learn<br>
> how to multiply matrices so I can do it all in a single step!<br>
><br>
> I just tried to import them into GRASS 7.0.3RC1 using v.in.ogr and got the<br>
> following error message:<br>
><br>
> ERROR 1: Reading EWKB with 4-dimensional coordinates (XYZM) is not supported<br>
><br>
> Can GRASS not import PostGIS 3D data?<br>
<br>
</span>This error (given the fact that it is "ERROR 1" must come from OGR and<br>
not from GRASS GIS as we don't use numbers).<br>
Hence, ogr2ogr would likely report the same.<br>
<span class=""><br>
> Is it the vertical Polygon with zero area topology problem?<br>
><br>
> What other options do I have to get 3D-cross-section data into GRASS?<br>
><br>
> I tried to use the ST_Extrude (to avoid topology problems) command in<br>
> PostGIS and Grass's v.in.ogr didn't like that either. Is there another way<br>
> to get 3D objects into GRASS?<br>
<br>
</span>Which OGR version is installed? As soon as it works there I am<br>
confident that it also works with v.in.ogr.<br>
<br>
I see some discussion about the problem also here (random search result):<br>
<br>
<a href="https://lists.osgeo.org/pipermail/gdal-dev/2015-February/041056.html" rel="noreferrer" target="_blank">https://lists.osgeo.org/pipermail/gdal-dev/2015-February/041056.html</a><br>
<span class=""><br>
> I want to use v.hull on the vertical cross-sectional data and symbolize some<br>
> geochemical drill data in Grass's 3D viewer.<br>
<br>
</span>IMHO it needs to be solved in OGR first (maybe it is in the latest<br>
version, I didn't verify).<br>
<br>
Best,<br>
Markus<br>
</blockquote></div><br></div>