size = LPR.GetSize();
// origin = LPR.GetIndex(); //in->GetOrigin();
for (unsigned int i=0;i<Dimension;++i)
- {
- origin[i] = LPR.GetIndex()[i]*in->GetSpacing()[i];
- spacing[i] = bbGetInputSpacing()[i];
- double tmp = (LPR.GetSize()[i]*in->GetSpacing()[i]/spacing[i] ) + 0.5;
- size[i] = (long)floor(tmp);
+ {
+ origin[i] = LPR.GetIndex()[i]*in->GetSpacing()[i];
+ spacing[i] = bbGetInputSpacing()[i];
+ double tmp = (LPR.GetSize()[i]*in->GetSpacing()[i]/spacing[i] ) + 0.5;
+ size[i] = (long)floor(tmp);
// size[i] = (long)lrint(LPR.GetSize()[i]*in->GetSpacing()[i]/spacing[i]);
}
// We create an interpolator of the found type
typename InterpolatorType::Pointer interpolator = InterpolatorType::New();
filter->SetInterpolator( interpolator );
-printf("EED ResampleImageFilter NearestNeighbor \n");
}
else if ( bbGetInputInterpolation() == "BSpline") {
-printf("EED ResampleImageFilter BSpline \n");
typedef itk::BSplineInterpolateImageFunction < ImageType, double > InterpolatorType;
// We create an interpolator of the found type
typename InterpolatorType::Pointer interpolator = InterpolatorType::New();
} //end else if
// Interpolation
else { // if ( bbGetInputInterpolation() == "Linear" ) {
-printf("EED ResampleImageFilter Linear \n");
typedef itk::LinearInterpolateImageFunction < ImageType, double > InterpolatorType;
// We create an interpolator of the found type
typename InterpolatorType::Pointer interpolator = InterpolatorType::New();