Hi David.<div><br></div><div>Thank you for replying. I fixed the problem replacing the line "cast -> SetOutputScalarTypeToDouble();" to "cast -> SetOutputScalarTypeToFloat();".</div><div><br></div>
<div>Not sure why float type is required, but it works.</div><div><br></div><div>Best,</div><div><div><div><font color="#555555" face="Helvetica, Arial, sans-serif"><span style="line-height:16px"><b>Matheus Viana</b></span></font></div>
<div><font color="#555555" face="Helvetica, Arial, sans-serif"><span style="line-height:16px"><b>Postdoctoral Research Employee</b></span></font></div><div><b><span style="color:rgb(85,85,85);line-height:16px;font-family:Helvetica,Arial,sans-serif;font-size:13px">Developmental and Cell Biology</span></b></div>
<div><font color="#555555" face="Helvetica, Arial, sans-serif"><span style="line-height:16px"><b>University of California Irvine</b></span></font></div><div><font color="#555555" face="Helvetica, Arial, sans-serif"><span style="line-height:16px"><b><br>
</b></span></font><b><i><br></i></b></div></div><br>
<br><br><div class="gmail_quote">2012/8/15 David Doria <span dir="ltr"><<a href="mailto:daviddoria@gmail.com" target="_blank">daviddoria@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
On Wed, Aug 15, 2012 at 6:37 PM, matheus_viana <<a href="mailto:vianamp@gmail.com">vianamp@gmail.com</a>> wrote:<br>
> Hi Folks,<br>
><br>
> I'm trying to sum two ImageDatas, but I'm getting a blank image as result.<br>
> This is my code:<br>
><br>
> vtkImageCast *cast = vtkImageCast::New();<br>
> cast -> SetInput(ImageData);<br>
> cast -> SetOutputScalarTypeToDouble();<br>
> cast -> Update();<br>
><br>
> vtkImageGaussianSmooth *Gauss = vtkImageGaussianSmooth::New();<br>
> Gauss -> ReleaseDataFlagOff();<br>
> Gauss -> SetInput(cast -> GetOutput());<br>
> Gauss -> SetDimensionality(3);<br>
> Gauss -> SetRadiusFactors(1,1,1);<br>
> Gauss -> SetStandardDeviations(1,1,1);<br>
> Gauss -> Update();<br>
><br>
> vtkImageVariance3D *Variance = vtkImageVariance3D::New();<br>
> Variance -> SetInput(cast->GetOutput());<br>
> Variance -> SetKernelSize(3,3,3);<br>
> Variance -> Update();<br>
><br>
> vtkImageMathematics *mat = vtkImageMathematics::New();<br>
> mat -> SetOperationToAdd();<br>
> mat -> SetInput1(Gauss->GetOutput());<br>
> mat -> SetInput2(Variance->GetOutput());<br>
> mat -> Update();<br>
><br>
> You guys have any idea?<br>
<br>
Please make a self-contained, compilable example of the problem. You<br>
can start with this as a shell:<br>
<br>
<a href="http://www.vtk.org/Wiki/VTK/Examples/Cxx/Images/ImageMathematics" target="_blank">http://www.vtk.org/Wiki/VTK/Examples/Cxx/Images/ImageMathematics</a><br>
<span class="HOEnZb"><font color="#888888"><br>
David<br>
</font></span></blockquote></div><br></div>