- filter->SetOutputSpacing(reader->GetOutput()->GetSpacing());
- const SizeType& input_size = reader->GetOutput()->GetLargestPossibleRegion().GetSize();
- SizeType output_size;
- for (unsigned int i = 0; i < Dimension; i++)
- output_size[i] = input_size[i];
- filter->SetOutputSize(output_size);
+
+ if (!m_ArgsInfo.useSizeAsReference_flag) {
+ filter->SetOutputSpacing(reader->GetOutput()->GetSpacing());
+ if (m_Verbose) {
+ std::cout << "Warning: The image spacing differs between the MHD file and the input dicom series. Performing resampling with default options using spacing as reference (for advanced options, use clitkResampleImage)." << std::endl;
+ std::cout << "MHD -> " << input->GetSpacing() << std::endl;
+ std::cout << "dicom -> " << reader->GetOutput()->GetSpacing() << std::endl;
+ }
+ }
+ else {
+ const SizeType& dicom_size = reader->GetOutput()->GetLargestPossibleRegion().GetSize();
+ SizeType output_size;
+ for (unsigned int i = 0; i < Dimension; i++)
+ output_size[i] = dicom_size[i];
+ filter->SetOutputSize(output_size);
+ if (m_Verbose) {
+ std::cout << "Warning: The image size differs between the MHD file and the input dicom series. Performing resampling with default options using size as reference (for advanced options, use clitkResampleImage)." << std::endl;
+ std::cout << "MHD -> " << input->GetLargestPossibleRegion().GetSize() << std::endl;
+ std::cout << "dicom -> " << reader->GetOutput()->GetLargestPossibleRegion().GetSize() << std::endl;
+ }
+ }
+