<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EstiloDeEmail17
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:8.5in 11.0in;
margin:70.85pt 85.05pt 70.85pt 85.05pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=PT-BR link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US'>Many thanks!<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US'>You saved me!!<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US'><o:p> </o:p></span></p><p class=MsoNormal><b><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri",sans-serif'>De:</span></b><span lang=EN-US style='font-size:11.0pt;font-family:"Calibri",sans-serif'> Vincent Cloarec [mailto:vcloarec@gmail.com] <br><b>Enviada em:</b> Sunday, April 26, 2020 11:35 AM<br><b>Para:</b> Jorge Almerio <jorgealmerio@yahoo.com.br><br><b>Cc:</b> qgis-developer <qgis-developer@lists.osgeo.org>; qgis-developer-bounces@lists.osgeo.org<br><b>Assunto:</b> Re: [QGIS-Developer] Mesh layer<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p><div><div><div><p class=MsoNormal><span lang=EN-US>Hi,<o:p></o:p></span></p></div><div><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p></div><div><p class=MsoNormal><span lang=EN-US>With QGIS master you can use </span><a href="https://qgis.org/pyqgis/master/core/QgsMeshLayer.html#qgis.core.QgsMeshLayer.updateTriangularMesh"><span lang=EN-US>updateTriangularMesh.</span></a><span lang=EN-US><o:p></o:p></span></p></div><div><p class=MsoNormal><span lang=EN-US>Before, there is a workaround that consists to just call </span><a href="https://qgis.org/pyqgis/3.12/core/QgsMeshLayer.html#qgis.core.QgsMeshLayer.createMapRenderer"><span lang=EN-US>createMapRenderer</span></a><span lang=EN-US> without using the renderer.<o:p></o:p></span></p></div><div><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p></div><div><p class=MsoNormal>Vincent<o:p></o:p></p></div></div><p class=MsoNormal><o:p> </o:p></p><div><div><p class=MsoNormal>Le sam. 25 avr. 2020 à 11:32, Jorge Almerio <<a href="mailto:jorgealmerio@yahoo.com.br">jorgealmerio@yahoo.com.br</a>> a écrit :<o:p></o:p></p></div><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0in;margin-bottom:5.0pt'><p class=MsoNormal><span lang=EN-US>Hi everyone,<br><br>Sometimes I got "NaN" return when trying to get scalar value to a point from a Mesh layer. I guess the reason is that it was not rendered previously as said in the documentation.<br><br><</span><a href="https://qgis.org/api/classQgsMeshLayer.html#a07cb62f04b64979500d3873afb580216" target="_blank"><span lang=EN-US>https://qgis.org/api/classQgsMeshLayer.html#a07cb62f04b64979500d3873afb580216</span></a><span lang=EN-US>><br>QgsMeshDatasetValue:<br>"Note:<br> It uses previously cached and indexed triangular mesh and so if the layer has not been rendered previously (e.g. when used in a script) it returns NaN value"<br><br>When I turn ON and OFF the loaded mesh layer in the active project it Works!<br><br>I would like to know if there is a proper way to render the mesh layer, before get values.<br><br>Below is my code:<br><br>proj = QgsProject.instance()<br>layerPto = iface.activeLayer() # a Point layer<br>layer = proj.mapLayersByName("mesh01")[0] # a Mesh layer<br>dataset = QgsMeshDatasetIndex(0,0)<br>for feat in layerPto.getSelectedFeatures():<br> ptoGeo = feat.geometry()<br> pt = ptoGeo.asPoint()<br> value = layer.datasetValue(dataset, pt).scalar()<br> print('elevation={}'.format(value))<br><br><br><br>_______________________________________________<br>QGIS-Developer mailing list<br></span><a href="mailto:QGIS-Developer@lists.osgeo.org" target="_blank"><span lang=EN-US>QGIS-Developer@lists.osgeo.org</span></a><span lang=EN-US><br>List info: </span><a href="https://lists.osgeo.org/mailman/listinfo/qgis-developer" target="_blank"><span lang=EN-US>https://lists.osgeo.org/mailman/listinfo/qgis-developer</span></a><span lang=EN-US><br>Unsubscribe: </span><a href="https://lists.osgeo.org/mailman/listinfo/qgis-developer" target="_blank"><span lang=EN-US>https://lists.osgeo.org/mailman/listinfo/qgis-developer</span></a><span lang=EN-US><o:p></o:p></span></p></blockquote></div></div></div></body></html>