- typedef itk::JoinSeriesImageFilter<SliceType, ImageType> JointSeriesFilterType;
- typename JointSeriesFilterType::Pointer jointFilter = JointSeriesFilterType::New();
- for(int i=0; i<mInputSlices.size(); i++) {
- jointFilter->SetInput(i, mInputSlices[i]);
+ DD(this->GetIntermediateSpacing());
+ DD(this->GetResampleBeforeRelativePositionFilter());
+ DD("End slice");
+
+ typedef itk::JoinSeriesImageFilter<SliceType, ImageType> JoinSeriesFilterType;
+ typename JoinSeriesFilterType::Pointer joinFilter = JoinSeriesFilterType::New();
+ joinFilter->SetOrigin(input->GetOrigin()[GetDirection()]);
+ joinFilter->SetSpacing(input->GetSpacing()[GetDirection()]);
+ for(unsigned int i=0; i<mInputSlices.size(); i++) {
+ // DD(mInputSlices[i]->GetLargestPossibleRegion().GetIndex());
+// DD(mInputSlices[i]->GetLargestPossibleRegion().GetSize());
+// DD(mInputSlices[i]->GetRequestedRegion().GetIndex());
+// DD(mInputSlices[i]->GetRequestedRegion().GetSize());
+ joinFilter->PushBackInput(mInputSlices[i]);
+ //SetInput(i, mInputSlices[i]);