hi, I have a small remark concerning the itkResampleImageFilter: In its ThreadedGenerateData method an itk::ImageRegionIterator is used to walk the output region. Since for every pixel the index is needed as well, it would be more efficient to use the itk::ImageRegionIteratorWithIndex, wouldn't it? regards, Martin