Hi Liam,<br><br>this is how I usually add arrays to image data<br><br> vtkDataArray *Velocity;<br> Velocity = Velocity->CreateDataArray(VTK_DOUBLE);<br> Velocity->SetNumberOfComponents(3);<br> Velocity->SetNumberOfTuples(grid->GetNumberOfPoints());<br>
Velocity->SetName("vertexVelocity");<br> grid->GetPointData()->AddArray(Velocity);<br><br> vtkDataArray *Density;<br>
Density =Density ->CreateDataArray(VTK_DOUBLE);<br>
Density->SetNumberOfComponents(1);<br>
Density ->SetNumberOfTuples(grid->GetNumberOfPoints());<br>
Density->SetName("Density");<br>
grid->GetPointData()->AddArray(Density);<br><br>I checked my arrays, they look ok.<br>I was able to calculate avg just using for-loop,<br>but I'd still like to understand/use vtk functions for<br>these tasks.<br>
<br>Thank you,<br>ula <br>
<br><br><br><div class="gmail_quote">On Fri, Apr 29, 2011 at 3:07 AM, Liam Kurmos <span dir="ltr"><<a href="mailto:quantum.leaf@gmail.com">quantum.leaf@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
hi Ula,<br>
<br>
can i see the code for how you add the 2 fields to the imagedata?<br>
<br>
thanks,<br>
<br>
Liam<br>
<div><div></div><div class="h5"><br>
<br>
On Fri, Apr 29, 2011 at 7:33 AM, Ula Popov <<a href="mailto:ula.popov@gmail.com">ula.popov@gmail.com</a>> wrote:<br>
> Hello all,<br>
><br>
> I have an object of class vtkImageData.<br>
> It has two arrays: density (scalar) and velocity (vector).<br>
><br>
> I want to know mean and median values of the density field.<br>
><br>
> I tried to use vtkImageAccumulate to achieve this goal, but<br>
> my program crashes on my input (actually, I don't know how<br>
> to specify on which of the arrays I want to apply. How<br>
> should I do it?)<br>
> When I tried to specify origin and spacing, I got similar result (seg fault)<br>
><br>
> vtkImageAccumulate *Accumulate = vtkImageAccumulate::New();<br>
><br>
> Accumulate->SetInput(grid);<br>
><br>
> Accumulate->Update();<br>
> <------- crashes here<br>
> double *mean = Accumulate->GetMean();<br>
><br>
><br>
> Thank you!<br>
> ula<br>
><br>
</div></div>> _______________________________________________<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 VTK FAQ at:<br>
> <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
><br>
><br>
</blockquote></div><br>