Hi Mathieu :<br><br>You need to override the method AllocateOutputs() and manually allocate each of the outputs.<br><br>The implementation of ImageToImageFilter is such that if the filter is not threaded, (as in it overrides GenerateData()), the implementation is responsible for allocating the outputs. If the filter is threaded (provides ThreadedGenerateData()), the superclass automatically allocates the outputs for you and ends up allocating all outputs of the type specified as the output image template parameter. In your case, the outputs are of different types and hence a memory access crash.
<br><br>The following illustrates the pipeline for threaded filters : <br> * 1) Allocate the output buffer --- (You must override this if the output types are different")<br> * 2) Call BeforeThreadedGenerateData()
<br> * 3) Spawn threads, calling ThreadedGenerateData() in each thread.<br> * 4) Call AfterThreadedGenerateData()<br><br><br>For details, please see itkImageSource.h <br><br>-- <br>Karthik Krishnan<br>R&D Engineer,
<br>Kitware Inc.
<br><br><div><span class="gmail_quote">On 8/22/07, <b class="gmail_sendername">Mathieu Malaterre</b> <<a href="mailto:mathieu.malaterre@gmail.com">mathieu.malaterre@gmail.com</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hello,<br><br><br> I am trying to write an ImagetoImageFilter that would produce two<br>types of output, two output are the same type as the input and the<br>last one (3rd) is floating point type. For some reason the filter is
<br>crashing. I suspect that the allocation did not go correctly for the<br>third output (second type).<br><br>Prototype is:<br><br>template <class TImage><br>class ITK_EXPORT MyImageFilter:<br> public ImageToImageFilter<TImage, TImage>
<br>{<br>...<br> /** Get the error image output of this process object. */<br> typedef Image<float,3> ErrorImageType;<br> ErrorImageType::Pointer GetErrorOutput()<br> { return static_cast<ErrorImageType *>(this->ProcessObject::GetOutput(2)); }
<br><br> /** Set the error image output of this process object. */<br> void SetErrorOutput(ErrorImageType *output)<br> { this->SetNthOutput(2, output); }<br>...<br>}<br><br>And in the ThreadedGenerateData :<br><br>
ThreadedGenerateData(const ImageRegionType& outputRegionForThread,<br> int threadId)<br>{<br> ErrorImageType::Pointer outputErrorPtr = this->GetErrorOutput();<br> OutputIterator outIt(outputPtr, outputRegionForThread);
<br> ...<br> //outIt.Set( residual ); // everything works if I comment out the 'Set'<br> ...<br>}<br><br>Finally here is the cstor:<br><br>::MyImageFilter()<br>{<br> this->ProcessObject::SetNumberOfRequiredOutputs(3);
<br> typename TImage::Pointer output = TImage::New();<br> this->ProcessObject::SetNthOutput(1, output.GetPointer());<br> ErrorImageType::Pointer e1 = ErrorImageType::New();<br> this->ProcessObject::SetNthOutput(2,
e1.GetPointer());<br>}<br><br>Comments welcome,<br><br><br>--<br>Mathieu<br>_______________________________________________<br>Insight-users mailing list<br><a href="mailto:Insight-users@itk.org">Insight-users@itk.org</a>
<br><a href="http://www.itk.org/mailman/listinfo/insight-users">http://www.itk.org/mailman/listinfo/insight-users</a><br></blockquote></div>