<div dir="ltr">Hi Martin and thanks for responding.<div><br></div><div>I should have said that I had studied that doc - and it does not actually answer my question. I also looked at the UGRID CF definitions but they are deliberately vague about this subject as well - suggesting two possible interpretations.<br><br>This problem only really comes to light when you start thinking about creating a 3D visualisation of the data.<br><br>The problem is that the doc has two interpretations, which are not compatible and neither is properly represented by the data - viz:</div><div><br></div><div><img src="cid:ii_kspzn1ft0" alt="v2.png" width="100" height="64">  <img src="cid:ii_kspzn1h01" alt="v1.png" width="100" height="65"> </div><div><br></div><div>I have chosen to interpret the data set as something more like this<br></div><div><br></div><div><img src="cid:ii_kspzn1i52" alt="v3.png" width="100" height="65"></div><div><br></div><div>I think that is correct but it is not clear from the docs and seems a little strange as a data structure<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, 24 Aug 2021 at 10:22, Martin Dobias <<a href="mailto:wonder.sk@gmail.com">wonder.sk@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"><div>Hi Paul</div><div><br></div><div>Probably the best description how the vertical levels in meshes are meant to work is this QGIS enhancement proposal:</div><div><a href="https://github.com/qgis/QGIS-Enhancement-Proposals/issues/158" target="_blank">https://github.com/qgis/QGIS-Enhancement-Proposals/issues/158</a></div><div><br></div><div>In short, within one face, if you have N vertical levels, there will be N-1 volumes in between those levels for that particular face. There are two important considerations that may not be obvious:</div><div>1. for each face, the vertical levels may be at different elevations/depths<br></div><div>2. for each face, there could be a different number of vertical levels</div><div>MDAL_D_volumesCount() returns _total_ number of volumes, not per-face number of volumes. Based on all this, the total number of vertical level entries therefore is "number of volumes" + "number of faces" because for each face we need that one extra vertical level.<br></div><div><br></div><div>Regards</div><div>Martin</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Aug 23, 2021 at 9:23 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">Hello all<div><br></div><div>Another question for you. I am looking at an app that converts from MDAL DataOnVolume dataset to a voxel-based representation.<br><br>That should not be particularly difficult. However, it has made apparent that it is not particularly clear what the "Vertical Levels" data in MDAL actually means.<br><br>My first understanding - that since the data values are volumetric the "vertical levels" represent the centroid of the volume being measured - is slightly undermined by the fact that there are (face count + volume count) values in the array which is too many for that interpretation.<br><br>The size of the array suggests that the interpretation should be one of layered meshes and the volumes are the spaces between the layers. That is fair enough but there is only one value per extruded face cell - when it is clear that the extruded face cell is not necessarily parallel to anything else. Is this vale then the vertical level of the centroid extruded face cell?<br><br>Does anyone have a definitive answer?</div></div>
_______________________________________________<br>
MDAL-Developer mailing list<br>
<a href="mailto:MDAL-Developer@lists.osgeo.org" target="_blank">MDAL-Developer@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/mdal-developer" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/mdal-developer</a><br>
</blockquote></div>
</blockquote></div>