<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Ken,<div class=""><br class=""></div><div class="">Your comments helped me a lot to see the issue and I appreciate you raising a bug on this.</div><div class=""><br class=""></div><div class="">Thanks again for your prompt replies.</div><div class=""><br class=""></div><div class="">Best,</div><div class=""><br class=""></div><div class="">—manoch</div><div class=""><br class=""></div><div class=""><br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On Jan 18, 2018, at 6:32 AM, Moreland, Kenneth <<a href="mailto:kmorel@sandia.gov" class="">kmorel@sandia.gov</a>> wrote:</div><br class="Apple-interchange-newline"><div class="">
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252" class="">
<div dir="auto" class="">
Manochehr,
<div class=""><br class="">
</div>
<div class="">ParaView always renders in a Cartesian space. So if you want the lat/lon/depth values to look like a sphere, the coordinates have to be converted and you will loose the original values. I can see the point of your use case. It just never came up before.
It seems reasonable to be able to save the original lat/lon values in a field, so I raised a bug for it.</div>
<div class=""><br class="">
</div>
<div class=""><a href="https://gitlab.kitware.com/paraview/paraview/issues/17943" class="">https://gitlab.kitware.com/paraview/paraview/issues/17943</a></div>
<div class=""><br class="">
</div>
<div class="">As far as getting the coordinates in a script, it is possible but not recommended. If you need to access field or coordinate values directly, it is usually best to do so in a programmable filter.</div>
<div class=""><br class="">
</div>
<div class="">-Ken<br class="">
<br class="">
<div class="">Sent from my iPad</div>
<div class=""><br class="">
On Jan 18, 2018, at 7:06 AM, Manochehr Bahavar <<a href="mailto:manoch@iris.washington.edu" class="">manoch@iris.washington.edu</a>> wrote:<br class="">
<br class="">
</div>
<blockquote type="cite" class="">
<div class="">Hi Aashish,
<div class=""><br class="">
</div>
<div class="">Thank you for your reply. That will be great if it gets looked into:) Going back and forth between coordinate systems may be a little confusing, but I will give it a try and see how it works out.</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class="">Could you also point me in the right direction for getting the latitude, longitude and depth arrays? As I mentioned in my reply to Ken (below), using Python, I go as far as getting the extents and bounds but cannot get the actual values.</div>
<div class=""><br class="">
</div>
<div class="">Thank you again for your quick reply and useful suggestions.</div>
<div class=""><br class="">
</div>
<div class="">Best,</div>
<div class=""><br class="">
</div>
<div class="">—manoch</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class=""><br class="">
<blockquote type="cite" class="">
<div class="">On Jan 17, 2018, at 10:29 PM, Aashish Chaudhary <<a href="mailto:aashish.chaudhary@kitware.com" class="">aashish.chaudhary@kitware.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" class="">Hi Manoch,
<div class=""><br class="">
</div>
<div class="">I think what you want is not possible right now but let us get back to you on it as something we might want to look into before we can say it for sure. One workaround would be you load the data (uncheck the sphere coordinate system), apply the
filter, and then write the data back as NetCDF. It is not ideal by any means. You can also look into vtkSphericalTransform (if you want to convert data back to spherical coordinate system) but you will not get the exact shape of earth which is a ellipsoid
(which I am not sure if you would be interested in anyways). </div>
<div class=""><br class="">
</div>
<div class="">Hope this helps. </div>
<br class="">
<div class="gmail_quote">
<div dir="ltr" class="">On Wed, Jan 17, 2018 at 7:12 PM Manochehr Bahavar <<a href="mailto:manoch@iris.washington.edu" target="_blank" class="">manoch@iris.washington.edu</a>> wrote:<br class="">
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word;line-break:after-white-space" class=""><font face="Helvetica" size="2" class="">Ken,</font>
<div class=""><font face="Helvetica" size="2" class=""><br class="">
</font></div>
<div class=""><font face="Helvetica" size="2" class="">Thanks for your reply. The problem with that approach is that all the features I am interested in make sense if they are plotted on a sphere (they approximate the Earth surface). So, there is no builtin
solution for that?</font></div>
<div class=""><font face="Helvetica" size="2" class=""><br class="">
</font></div>
<div class=""><font face="Helvetica" size="2" class="">One more question: How can I access the actual latitude, longitude and depth parameter values using Python. I can see their extent and bounds, but I can not see how I can access the individual values:</font></div>
<div class=""><font face="Helvetica" size="2" class=""><br class="">
</font></div>
<div class=""><font face="Helvetica" size="2" class=""><span style="color:rgb(51,51,51)" class="">temp.GetBounds()</span><br style="color:rgb(51,51,51)" class="">
<span style="color:rgb(51,51,51)" class="">(0.0, 359.0, -89.0, 89.0, -2750.0, 50.0)</span><br style="color:rgb(51,51,51)" class="">
<span style="color:rgb(51,51,51)" class="">temp.GetDimensions()</span><br style="color:rgb(51,51,51)" class="">
<span style="color:rgb(51,51,51)" class="">(360, 179, 57)</span><br style="color:rgb(51,51,51)" class="">
<span style="color:rgb(51,51,51)" class="">temp.GetExtent()</span><br style="color:rgb(51,51,51)" class="">
<span style="color:rgb(51,51,51)" class="">(0, 359, 0, 178, 0, 56)</span><br style="color:rgb(51,51,51)" class="">
<span style="color:rgb(51,51,51)" class="">temp.GetCenter()</span><br style="color:rgb(51,51,51)" class="">
<span style="color:rgb(51,51,51)" class="">(179.5, 0.0, -1350.0)</span><br style="color:rgb(51,51,51)" class="">
<span style="color:rgb(51,51,51)" class="">temp.GetNumberOfPoints()</span><br style="color:rgb(51,51,51)" class="">
<span style="color:rgb(51,51,51)" class="">3673080L</span></font></div>
<div class=""><font color="#333333" face="Helvetica" size="2" class=""><br class="">
</font></div>
<div class=""><font color="#333333" face="Helvetica" size="2" class="">Really appreciate your help</font></div>
<div class=""><font color="#333333" face="Helvetica" size="2" class=""><br class="">
</font></div>
<div class=""><font color="#333333" face="Helvetica" size="2" class="">—manoch</font></div>
</div>
<div style="word-wrap:break-word;line-break:after-white-space" class="">
<div class=""><font color="#333333" face="Lato, Helvetica, Arial, freesans, clean, sans-serif" class=""><span style="font-size:15px" class=""><br class="">
</span></font>
<div class=""><br class="">
<blockquote type="cite" class="">
<div class="">On Jan 17, 2018, at 3:57 PM, Moreland, Kenneth <<a href="mailto:kmorel@sandia.gov" target="_blank" class="">kmorel@sandia.gov</a>> wrote:</div>
<br class="m_5856616774579461781m_-6687049364179838243Apple-interchange-newline">
<div class="">
<div class="">Manochehr,<br class="">
<br class="">
If you uncheck the “Spherical Coordiantes” option, then the longitude, latitude, and depth values will be preserved as x, y, and z coordinates. (Instead of translated to a sphere, the data will appear in a rectangle.)<br class="">
<br class="">
-Ken<br class="">
<br class="">
<br class="">
On 1/17/18, 4:34 PM, "ParaView on behalf of Manochehr Bahavar" <<a href="mailto:paraview-bounces@paraview.org" target="_blank" class="">paraview-bounces@paraview.org</a> on behalf of
<a href="mailto:manoch@iris.washington.edu" target="_blank" class="">manoch@iris.washington.edu</a>> wrote:<br class="">
<br class="">
Hello,<br class="">
<br class="">
I am trying to see if I can visualize my netCDF Earth models (longitude, latitude,depth) using ParaView. I have been able to load the models, with the "Spherical Coordinates" option selected. However it appears that I now lose my coordinates and depth values
and as a result to extract a subset, I have to use the depth index (not the actual depth value). Is there a way to preserve the latitude, longitude and depth values and work with ParaView using the actual values (not their index) directly?<br class="">
<br class="">
Thank you,<br class="">
<br class="">
—manoch<br class="">
_______________________________________________<br class="">
Powered by <a href="http://www.kitware.com/" target="_blank" class="">www.kitware.com</a><br class="">
<br class="">
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank" class="">
http://www.kitware.com/opensource/opensource.html</a><br class="">
<br class="">
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" target="_blank" class="">
http://paraview.org/Wiki/ParaView</a><br class="">
<br class="">
Search the list archives at: <a href="http://markmail.org/search/?q=ParaView" target="_blank" class="">
http://markmail.org/search/?q=ParaView</a><br class="">
<br class="">
Follow this link to subscribe/unsubscribe:<br class="">
<a href="https://paraview.org/mailman/listinfo/paraview" target="_blank" class="">https://paraview.org/mailman/listinfo/paraview</a><br class="">
<br class="">
<br class="">
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
_______________________________________________<br class="">
Powered by <a href="http://www.kitware.com/" rel="noreferrer" target="_blank" class="">
www.kitware.com</a><br class="">
<br class="">
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank" class="">
http://www.kitware.com/opensource/opensource.html</a><br class="">
<br class="">
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" rel="noreferrer" target="_blank" class="">
http://paraview.org/Wiki/ParaView</a><br class="">
<br class="">
Search the list archives at: <a href="http://markmail.org/search/?q=ParaView" rel="noreferrer" target="_blank" class="">
http://markmail.org/search/?q=ParaView</a><br class="">
<br class="">
Follow this link to subscribe/unsubscribe:<br class="">
<a href="https://paraview.org/mailman/listinfo/paraview" rel="noreferrer" target="_blank" class="">https://paraview.org/mailman/listinfo/paraview</a><br class="">
</blockquote>
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</blockquote>
</div>
</div>
</div></blockquote></div><br class=""></div></body></html>