<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'>
So, I´m really fed up with this:<br>Why does this filter link without problems :<br><br>vtkImageHorizontalAverage.h:<br>#ifndef __vtkImageHorizontalAverage_h<br>#define __vtkImageHorizontalAverage_h<br><br>// additional macro for paraview running under windows<br>//#define vtkImaging_EXPORTS<br><br>#include "vtkSimpleImageToImageFilter.h"<br><br>class VTK_IMAGING_EXPORT vtkImageHorizontalAverage: public vtkSimpleImageToImageFilter<br>{<br>public:<br> static vtkImageHorizontalAverage *New();<br> vtkTypeRevisionMacro(vtkImageHorizontalAverage,vtkSimpleImageToImageFilter);<br>// void PrintSelf(ostream& os, vtkIndent indent);<br><br> // Set/Get the isNormalize value.<br> vtkSetMacro(Normalize,int);<br> vtkGetMacro(Normalize,int);<br> vtkBooleanMacro(Normalize, int);<br> <br>protected:<br> vtkImageHorizontalAverage();<br> ~vtkImageHorizontalAverage() {};<br><br> int Normalize;<br> <br> virtual void SimpleExecute(vtkImageData* input, vtkImageData* output);<br><br>private:<br> vtkImageHorizontalAverage(const vtkImageHorizontalAverage&); // Not implemented.<br> void operator=(const vtkImageHorizontalAverage&); // Not implemented.<br>};<br><br>#endif<br><br><br>vtkImageHorizontalAverage.cpp:<br>#include "vtkImageHorizontalAverage.h"<br><br>#include "vtkImageData.h"<br>#include "vtkObjectFactory.h"<br><br>vtkCxxRevisionMacro(vtkImageHorizontalAverage, "$Revision: 1.10 $");<br>vtkStandardNewMacro(vtkImageHorizontalAverage);<br><br><br>// constructor<br>vtkImageHorizontalAverage::vtkImageHorizontalAverage()<br>{<br> this->Normalize = 0;<br>}<br><br>// The switch statement in Execute will call this method with<br>// the appropriate input type (IT). Note that this example assumes<br>// that the output data type is the same as the input data type.<br>// This is not always the case.<br>template <class IT><br>void vtkImageHorizontalAverageExecute(vtkImageHorizontalAverage *self,<br> vtkImageData* input,<br> vtkImageData* output,<br> IT* inPtr, IT* outPtr)<br>{<br> int i,j,k,dims[3];<br> input->GetDimensions(dims);<br> if (input->GetScalarType() != output->GetScalarType())<br> {<br> vtkGenericWarningMacro(<< "Execute: input ScalarType, " << input->GetScalarType()<br> << ", must match out ScalarType " << output->GetScalarType());<br> return;<br> }<br><br> double *horizontal = new double[dims[0]];<br> double *horizontal_sigma = new double[dims[0]];<br> int index;<br> // compute the horizontal average<br> for (i=0; i<dims[0]; i++) { <br> horizontal[i] = 0.0;<br> for (j=0; j<dims[1]; j++) {<br> for (k=0; k<dims[2]; k++) {<br> index = k*(dims[1]*dims[0]) + j*dims[0] + i;<br> horizontal[i] += inPtr[index];<br> }<br> }<br> horizontal[i] = horizontal[i]/(dims[1]*dims[2]);<br> }<br><br> //subtract horizontal average<br> for (i=0; i<dims[0]; i++) { <br> for (j=0; j<dims[1]; j++) {<br> for (k=0; k<dims[2]; k++) {<br> index = k*(dims[1]*dims[0]) + j*dims[0] + i;<br> outPtr[index] = inPtr[index] - horizontal[i];<br> }<br> }<br> }<br> <br> // normalize with sigma <br> if (self->GetNormalize()) {<br> // compute sigma<br> for (i=0; i<dims[0]; i++) { <br> horizontal_sigma[i] = 0.0;<br> for (j=0; j<dims[1]; j++) {<br> for (k=0; k<dims[2]; k++) {<br> index = k*(dims[1]*dims[0]) + j*dims[0] + i;<br> horizontal_sigma[i] += (inPtr[index])*(inPtr[index]);<br> }<br> }<br> horizontal_sigma[i] /= (dims[1]*dims[2]-1);<br> horizontal_sigma[i] = sqrt(horizontal_sigma[i]); <br> }<br> <br> // now normalize<br> for (i=0; i<dims[0]; i++) { <br> for (j=0; j<dims[1]; j++) {<br> for (k=0; k<dims[2]; k++) {<br> index = k*(dims[1]*dims[0]) + j*dims[0] + i;<br> outPtr[index] /= horizontal_sigma[i];<br> }<br> }<br> }<br> }<br> <br> delete [] horizontal;<br> delete [] horizontal_sigma;<br>}<br><br>void vtkImageHorizontalAverage::SimpleExecute(vtkImageData* input,<br> vtkImageData* output)<br>{<br> // copy information<br> output->DeepCopy(input);<br><br> void* inPtr = input->GetScalarPointer();<br> void* outPtr = output->GetScalarPointer();<br><br> switch(output->GetScalarType())<br> {<br> // This is simple a #define for a big case list. It handles<br> // all data types vtk can handle.<br> vtkTemplateMacro5(vtkImageHorizontalAverageExecute, this, input, output, <br> (VTK_TT *)(inPtr), (VTK_TT *)(outPtr));<br> default:<br> vtkGenericWarningMacro("Execute: Unknown input ScalarType");<br> return;<br> }<br>}<br><br>and the other one does not???????????????????????????? (The other one - vtkHorizontalAverage - is posted in the post before)<br>Both should do the same thing, compute the horizontal average of the points of the grid, the only difference should be the data type of the input file! one is structured_points, the other one structured_grid!!!!!!<br><br>I´d be very grateful for any comments!<br>thx,<br>NH<br><br><br /><hr />Express yourself instantly with MSN Messenger! <a href='http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/' target='_new'>MSN Messenger</a></body>
</html>