I don&#39;t have any problem running Median or Gradient on image of type<br>double. If you are allowed to share your data, I may give a try here.<br><br>Maybe are you running out of memory? You can try to extract a subset<br>
of your data and see if the segfault occurs.<br><br>Jerome<br><br><div class="gmail_quote">2011/5/20 Dominik Szczerba <span dir="ltr">&lt;<a href="mailto:dominik@itis.ethz.ch">dominik@itis.ethz.ch</a>&gt;</span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hi Jerome,<br>
<br>
Thanks for the references, but there does not to be any solution in there?<br>
And yes, &quot;strange extents&quot; are also taking place in some cases, i.e.,<br>
extents are actually correct, but the content of the image data is<br>
smaller if the new spacing is finer than the original.<br>
Due to the lack of detailed description of this whole mechanism I can<br>
not say if this is my lack of knowledge or a bug.<br>
<br>
Regarding segfaults, the only non-standard thing I can say about my<br>
input image data is that it contains a scalar type &quot;double&quot;. The said<br>
filters will work with image generated by Sources/Wavelet, but the<br>
type of that one is float. I am afraid something is hardcoded<br>
somewhere.<br>
<br>
Regards,<br>
<font color="#888888">Dominik<br>
</font><div><div></div><div class="h5"><br>
On Fri, May 20, 2011 at 8:25 AM, Jérôme &lt;<a href="mailto:jerome.velut@gmail.com">jerome.velut@gmail.com</a>&gt; wrote:<br>
&gt; Hi,<br>
&gt; regarding the vtkImageSlice exposition in ParaView, you should take<br>
&gt; a look at these threads:<br>
&gt; - issue report:<br>
&gt;    <a href="http://www.mail-archive.com/paraview@paraview.org/msg02225.html" target="_blank">http://www.mail-archive.com/paraview@paraview.org/msg02225.html</a><br>
&gt;    <a href="http://www.mail-archive.com/paraview@paraview.org/msg08260.html" target="_blank">http://www.mail-archive.com/paraview@paraview.org/msg08260.html</a><br>
&gt; - issue deep explanation:<br>
&gt;    <a href="http://www.cmake.org/pipermail/paraview/2010-June/017738.html" target="_blank">http://www.cmake.org/pipermail/paraview/2010-June/017738.html</a><br>
&gt;<br>
&gt; regarding the segfault, I feel surprise ! I use Median and Gradient very<br>
&gt; often, and I never noticed a bug. Can you give the step to reproduce ?<br>
&gt;<br>
&gt; Best regards,<br>
&gt; Jerome<br>
&gt;<br>
&gt;<br>
&gt; 2011/5/19 Dominik Szczerba &lt;<a href="mailto:dominik@itis.ethz.ch">dominik@itis.ethz.ch</a>&gt;<br>
&gt;&gt;<br>
&gt;&gt; I think my theory that it has to do with threading may be correct...<br>
&gt;&gt;<br>
&gt;&gt; I just tried instead of my custom plugin to run the built-in Gradient<br>
&gt;&gt; or Median on my vtkImageData and Paraview (3.11) silently segfaults.<br>
&gt;&gt; Both these built-in filters derive from vtkThreadedImageAlgorithm.<br>
&gt;&gt; Built-in vtkExtractVOI. however, works correctly, just as my other<br>
&gt;&gt; custom plugins deriving from vtkImageAlgorithm.<br>
&gt;&gt;<br>
&gt;&gt; Any pointers in this bug hunting are appreciated.<br>
&gt;&gt;<br>
&gt;&gt; Dominik<br>
&gt;&gt;<br>
&gt;&gt; On Thu, May 19, 2011 at 9:30 PM, Dominik Szczerba &lt;<a href="mailto:dominik@itis.ethz.ch">dominik@itis.ethz.ch</a>&gt;<br>
&gt;&gt; wrote:<br>
&gt;&gt; &gt; So far I was able to add a few simple custom plugins, but am slowly<br>
&gt;&gt; &gt; losing the grip.<br>
&gt;&gt; &gt; I am trying to resample an image like this:<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; &lt;SourceProxy name=&quot;ImageResample&quot; class=&quot;vtkImageResample&quot;<br>
&gt;&gt; &gt; label=&quot;ImageResample&quot;&gt;<br>
&gt;&gt; &gt; &lt;InputProperty name=&quot;Input&quot; command=&quot;SetInputConnection&quot;&gt;<br>
&gt;&gt; &gt;        &lt;ProxyGroupDomain name=&quot;groups&quot;&gt;<br>
&gt;&gt; &gt;                &lt;Group name=&quot;sources&quot;/&gt;<br>
&gt;&gt; &gt;                &lt;Group name=&quot;filters&quot;/&gt;<br>
&gt;&gt; &gt;        &lt;/ProxyGroupDomain&gt;<br>
&gt;&gt; &gt;        &lt;DataTypeDomain name=&quot;input_type&quot;&gt;<br>
&gt;&gt; &gt;                &lt;DataType value=&quot;vtkImageData&quot;/&gt;<br>
&gt;&gt; &gt;        &lt;/DataTypeDomain&gt;<br>
&gt;&gt; &gt; &lt;/InputProperty&gt;<br>
&gt;&gt; &gt; &lt;DoubleVectorProperty name=&quot;Magnification Factor&quot;<br>
&gt;&gt; &gt; command=&quot;SetAxisMagnificationFactor&quot; number_of_elements=&quot;3&quot;<br>
&gt;&gt; &gt; number_of_elements_per_command=&quot;1&quot; repeat_command=&quot;1&quot; use_index=&quot;1&quot;<br>
&gt;&gt; &gt; default_values=&quot;1 1 1&quot;/&gt;<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; I get these 2 problems:<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; 1) There appear two (!) outputs (named Output-0 and Output-1). I get<br>
&gt;&gt; &gt; this message:<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; ERROR: In<br>
&gt;&gt; &gt; /home/dsz/pack/ParaView/VTK/Filtering/vtkDemandDrivenPipeline.cxx,<br>
&gt;&gt; &gt; line 827<br>
&gt;&gt; &gt; vtkPVPostFilterExecutive (0x3aa4f60): Input for connection index 0 on<br>
&gt;&gt; &gt; input port index 0 for algorithm vtkPVPostFilter(0x2cfb850) is of type<br>
&gt;&gt; &gt; vtkImageStencilData, but a vtkDataSet is required.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; 2) The first output has the expected Spacing, but it does not contain<br>
&gt;&gt; &gt; the same scalar array as in the input. Instead, it contains a new<br>
&gt;&gt; &gt; array ImageScalars ranging [0,0]. I get this message:<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; Warning: In<br>
&gt;&gt; &gt; /home/dsz/pack/ParaView/ParaViewCore/VTKExtensions/vtkTexturePainter.cxx,<br>
&gt;&gt; &gt; line 179<br>
&gt;&gt; &gt; vtkTexturePainter (0x339fa70): Failed to locate selected scalars. Will<br>
&gt;&gt; &gt; use image scalars by default.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; I have been so far successful with other image algorithms, however<br>
&gt;&gt; &gt; they were all deriving from vtkImageAlgorithm. This one derives from<br>
&gt;&gt; &gt; vtkThreadedImageAlgorithm - could this have anything to do here?<br>
&gt;&gt; &gt; Could anyone point me into the right direction? After inspection of<br>
&gt;&gt; &gt; the xml files in the ParaView source tree I am still clueless.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; Thanks and regards,<br>
&gt;&gt; &gt; Dominik<br>
&gt;&gt; &gt;<br>
&gt;&gt; _______________________________________________<br>
&gt;&gt; Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
&gt;&gt;<br>
&gt;&gt; Visit other Kitware open-source projects at<br>
&gt;&gt; <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
&gt;&gt;<br>
&gt;&gt; Please keep messages on-topic and check the ParaView Wiki at:<br>
&gt;&gt; <a href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
&gt;&gt;<br>
&gt;&gt; Follow this link to subscribe/unsubscribe:<br>
&gt;&gt; <a href="http://www.paraview.org/mailman/listinfo/paraview" target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a><br>
&gt;<br>
&gt;<br>
</div></div></blockquote></div><br>