<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=iso-8859-1">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:Helvetica;
panose-1:2 11 6 4 2 2 2 2 2 4;}
@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;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 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;}
p
{mso-style-priority:99;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
span.apple-tab-span
{mso-style-name:apple-tab-span;}
span.apple-style-span
{mso-style-name:apple-style-span;}
span.EmailStyle20
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
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=EN-US link=blue vlink=purple style='word-wrap: break-word;
-webkit-nbsp-mode: space;-webkit-line-break: after-white-space'>
<div class=WordSection1>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Thanks all! I believe we have found numerous solutions.<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'>Alan<o:p></o:p></span></p>
<p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";
color:#1F497D'><o:p> </o:p></span></p>
<div>
<div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in'>
<p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>From:</span></b><span
style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'> Eric E. Monson
[mailto:emonson@cs.duke.edu] <br>
<b>Sent:</b> Friday, August 20, 2010 9:26 AM<br>
<b>To:</b> Scott, W Alan<br>
<b>Cc:</b> David E DeMarle; Aurélien Marsan; paraview@paraview.org list<br>
<b>Subject:</b> Re: [Paraview] nan<o:p></o:p></span></p>
</div>
</div>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>If you can use numpy, and you're using a recent-enough
version of ParaView, then you can also use the
super-spiffy-numpy-hidden-behind-the-scenes programmable filter API to replace
all of the attribute arrays with zero'd versions:<o:p></o:p></p>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<div>
<div>
<p class=MsoNormal>pdi = self.GetInputDataObject(0,0)<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal>pdo = self.GetOutputDataObject(0)<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal>pdo.CopyStructure(pdi)<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal style='-qt-paragraph-type:empty;-qt-block-indent:0'><o:p> </o:p></p>
</div>
<div>
<p class=MsoNormal>for att_name in inputs[0].PointData.keys():<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><span class=apple-tab-span> </span>naninc =
inputs[0].PointData[att_name]<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><span class=apple-tab-span> </span>zerod =
numpy.nan_to_num(naninc)<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><span class=apple-tab-span> </span>output.PointData.append(zerod,
att_name)<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<div>
<p class=MsoNormal>Talk to you later,<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal>-Eric<o:p></o:p></p>
</div>
<div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
<div>
<div>
<p class=MsoNormal><span style='font-size:9.0pt;font-family:"Helvetica","sans-serif";
color:black'>------------------------------------------------------<o:p></o:p></span></p>
</div>
<div>
<p class=MsoNormal><span style='font-size:9.0pt;font-family:"Helvetica","sans-serif";
color:black'>Eric E Monson<o:p></o:p></span></p>
</div>
<div>
<p class=MsoNormal><span style='font-size:9.0pt;font-family:"Helvetica","sans-serif";
color:black'>Duke Visualization Technology Group<o:p></o:p></span></p>
</div>
<div>
<p class=MsoNormal><span style='font-size:9.0pt;font-family:"Helvetica","sans-serif";
color:black'><o:p> </o:p></span></p>
</div>
</div>
<p class=MsoNormal><o:p> </o:p></p>
<div>
<div>
<p class=MsoNormal>On Aug 20, 2010, at 10:00 AM, Aurélien Marsan wrote:<o:p></o:p></p>
</div>
<p class=MsoNormal><br>
<br>
<o:p></o:p></p>
<p class=MsoNormal style='margin-bottom:12.0pt'>Hi, <br>
<br>
If you're using python, and if you can use the numpy library, you can use the
function numpy.nan_to_num too. <br>
<a
href="http://docs.scipy.org/doc/numpy/reference/generated/numpy.nan_to_num.html#numpy.nan_to_num">http://docs.scipy.org/doc/numpy/reference/generated/numpy.nan_to_num.html#numpy.nan_to_num</a><br>
<br>
Regards, <br>
<br>
Aurélien<o:p></o:p></p>
<div>
<p class=MsoNormal>2010/8/20 David E DeMarle <<a
href="mailto:dave.demarle@kitware.com">dave.demarle@kitware.com</a>><o:p></o:p></p>
<p class=MsoNormal>That said, the same trick may work in the standard
calculator filter<br>
with the expression:<br>
<br>
if(val=val, val, 0.0)<o:p></o:p></p>
<div>
<p class=MsoNormal style='margin-bottom:12.0pt'><br>
David E DeMarle<br>
Kitware, Inc.<br>
R&D Engineer<br>
28 Corporate Drive<br>
Clifton Park, NY 12065-8662<br>
Phone: 518-371-3971 x109<br>
<br>
<br>
<o:p></o:p></p>
</div>
<div>
<div>
<p class=MsoNormal>On Fri, Aug 20, 2010 at 9:26 AM, David E DeMarle<br>
<<a href="mailto:dave.demarle@kitware.com">dave.demarle@kitware.com</a>>
wrote:<br>
> You might write a python filter that iterates over all floating point
arrays and<br>
> replaces nan's with 0.<br>
><br>
> According to <a
href="http://stackoverflow.com/questions/944700/how-to-check-for-nan-in-python"
target="_blank">http://stackoverflow.com/questions/944700/how-to-check-for-nan-in-python</a><br>
> The most py version robust way to check for nan is:<br>
><br>
> def isNaN(num):<br>
> return num != num<br>
><br>
> David E DeMarle<br>
> Kitware, Inc.<br>
> R&D Engineer<br>
> 28 Corporate Drive<br>
> Clifton Park, NY 12065-8662<br>
> Phone: 518-371-3971 x109<br>
><br>
><br>
><br>
> On Thu, Aug 19, 2010 at 11:24 PM, Moreland, Kenneth <<a
href="mailto:kmorel@sandia.gov">kmorel@sandia.gov</a>> wrote:<br>
>> This might not be the best solution, but you can use the threshold
filter to<br>
>> remove NANs. A NAN will always fall outside the threshold range.<br>
>><br>
>> -Ken<br>
>><br>
>><br>
>> On 8/19/10 6:28 PM, "Scott, W Alan" <<a
href="mailto:wascott@sandia.gov">wascott@sandia.gov</a>> wrote:<br>
>><br>
>> Is there a way to convert nan’s to zeros in ParaView? I
have a user that is<br>
>> trying to use the integrate data filter, and it is having troubles
with<br>
>> NANs.<br>
>><br>
>> Thanks,<br>
>><br>
>> Alan<br>
>><br>
>><br>
>><br>
>><br>
>><br>
>> _______________________________________________<br>
>> Powered by <a href="http://www.kitware.com/" target="_blank">www.kitware.com</a><br>
>><br>
>> Visit other Kitware open-source projects at<br>
>> <a href="http://www.kitware.com/opensource/opensource.html"
target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
>><br>
>> Please keep messages on-topic and check the ParaView Wiki at:<br>
>> <a href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
>><br>
>> Follow this link to subscribe/unsubscribe:<br>
>> <a href="http://www.paraview.org/mailman/listinfo/paraview"
target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a><br>
>><br>
>><br>
><br>
_______________________________________________<br>
Powered by <a href="http://www.kitware.com/" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a
href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ParaView Wiki at: <a
href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.paraview.org/mailman/listinfo/paraview" target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a><o:p></o:p></p>
</div>
</div>
</div>
<p class=MsoNormal><br>
_______________________________________________<br>
Powered by <a href="http://www.kitware.com">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a
href="http://www.kitware.com/opensource/opensource.html">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ParaView Wiki at: <a
href="http://paraview.org/Wiki/ParaView">http://paraview.org/Wiki/ParaView</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.paraview.org/mailman/listinfo/paraview">http://www.paraview.org/mailman/listinfo/paraview</a><o:p></o:p></p>
</div>
<p class=MsoNormal><o:p> </o:p></p>
</div>
</div>
</body>
</html>