<div dir="ltr">Re-opening this and adding back the wider group( I had not noticed that they had dropped off).<div><br></div><div>The question is about saving 3D meshes with volumetric datasets. The context is around creating a python interface and using it for 3D volumetric data.</div><div><br></div><div>There are a number of questions coming out of this work - I will raise individual issues and PRs where appropriate but I wanted to discuss the issues wider:<br><br>1 There seems to be a general lack of drivers capable of writing any sort of data - according to the flags and therefore also according to the C API which checks the flags. UGRID, for instance, was mentioned but the Ugrid driver capabilities do not include the ability to write ANY data. Is that correct? It seems to contradict what has been said elsewhere.<br><br>There are three drivers that can write meshes (2DM, SELAFIN and UGRID). Of these only one can write data - SELAFIN - and only on vertices. This is very limited. Is it correct?</div><div><br></div><div>The flags also say that DataOnFace can be written by the TUFLOW FV driver but that meshes cannot, and that data can be written by the DAT drivers. But none of that helps with creating a new mesh with data of any type (let alone 3D).<br><br>2 There are NO drivers that admit to the ability to write DataOnVolumes datasets. <br><br>Is this true? </div><div><br></div><div>3 The C API MDAL_G_addDataset() method explicitly excludes (i.e. there is a validation check to enforce the exclusion) creation of DataOnVolume datasets. This does make the conversation earlier in the email chain about how to create these datasets incorrect. Is there a reason for this exclusion and is there a way around it.</div><div><br></div><div>Thanks for any pointers / comments</div><div><br></div><div>Paul</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 17 Jun 2021 at 10:22, Paul Harwood <<a href="mailto:runette@gmail.com">runette@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Ok, thanks. Got it now<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 17 Jun 2021 at 10:20, Peter Petrik <<a href="mailto:peter.petrik@lutraconsulting.co.uk" target="_blank">peter.petrik@lutraconsulting.co.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">base net is same as 2D mesh, verticals are just datasets defined in <a href="https://github.com/lutraconsulting/MDAL/blob/42eb61f880f99ca0aa89a72f5afc5193b2f061e2/mdal/api/mdal.h#L648" target="_blank">https://github.com/lutraconsulting/MDAL/blob/42eb61f880f99ca0aa89a72f5afc5193b2f061e2/mdal/api/mdal.h#L648</a> </div><div dir="ltr"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jun 17, 2021 at 11:07 AM Paul Harwood <<a href="mailto:runette@gmail.com" target="_blank">runette@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">To be more precise - I can see how to add volume based datasets but not how to add the vertical structure - i.e. the number of volumes per face of the base net.</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 17 Jun 2021 at 08:18, Paul Harwood <<a href="mailto:runette@gmail.com" target="_blank">runette@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">We probably need to change or clarify the README which says that 3D data is not supported for UGRID. Does that column mean something else ?<div><br></div><div>Is the support read-only? I cannot immediately see how to use the C API to write a layered 3D mesh.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, 17 Jun 2021 at 07:27, Peter Petrik <<a href="mailto:peter.petrik@lutraconsulting.co.uk" target="_blank">peter.petrik@lutraconsulting.co.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><a href="https://ugrid-conventions.github.io/ugrid-conventions/#3d-layered-mesh-topology" target="_blank">https://ugrid-conventions.github.io/ugrid-conventions/#3d-layered-mesh-topology</a> is supported <br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jun 16, 2021 at 8:38 PM Paul Harwood <<a href="mailto:runette@gmail.com" target="_blank">runette@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Sorry - in that last email when I said "is UGRID unstructured 3d MESH supported"? I meant "is UGRID LAYERED 3d mesh supported"!</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, 16 Jun 2021 at 19:35, Paul Harwood <<a href="mailto:runette@gmail.com" target="_blank">runette@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Peter<div><br></div><div>Can I just confirm. I look at the GH readme and it says that only XMDF supports 3d (i.e. I guess TUFLOW). There is an issue about UGRID that mentions a PR. Also the Readme says that XMDF support for 3D is read-only.<br><br>Is UGRID Unstructured 3D mesh supported?</div><div><br></div><div>Is any 3D RW?</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, 16 Jun 2021 at 14:59, Paul Harwood <<a href="mailto:runette@gmail.com" target="_blank">runette@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="auto">That makes sense. Thanks. I think that ugrid is what I want</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, 16 Jun 2021, 14:45 Peter Petrik, <<a href="mailto:peter.petrik@lutraconsulting.co.uk" target="_blank">peter.petrik@lutraconsulting.co.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">I think we only support TUFLOW 3d data and UGRID specification for 3D stuff, which limits your options quite a bit :)<div><br></div><div>P.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Jun 16, 2021 at 3:28 PM Paul Harwood <<a href="mailto:runette@gmail.com" rel="noreferrer" target="_blank">runette@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi all<div><br></div><div>I am thinking about voxel based 3D data in meshes. </div><div><br></div><div>I am aware of this <a href="https://github.com/qgis/QGIS-Enhancement-Proposals/issues/158" rel="noreferrer" target="_blank">https://github.com/qgis/QGIS-Enhancement-Proposals/issues/158</a></div><div><br></div><div>The question I have is a simple one :<br><br>Absent any other constraint, what would be the "best" format to use to hold and transfer voxel based 3D mesh data?</div></div>
_______________________________________________<br>
MDAL-Developer mailing list<br>
<a href="mailto:MDAL-Developer@lists.osgeo.org" rel="noreferrer" target="_blank">MDAL-Developer@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/mdal-developer" rel="noreferrer noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/mdal-developer</a><br>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div>