normalise.h

Go to the documentation of this file.

/* Copyright (c) 2008-2022 the MRtrix3 contributors.
 *
 * This Source Code Form is subject to the terms of the Mozilla Public
 * License, v. 2.0. If a copy of the MPL was not distributed with this
 * file, You can obtain one at http://mozilla.org/MPL/2.0/.
 *
 * Covered Software is provided under this License on an "as is"
 * basis, without warranty of any kind, either expressed, implied, or
 * statutory, including, without limitation, warranties that the
 * Covered Software is free of defects, merchantable, fit for a
 * particular purpose or non-infringing.
 * See the Mozilla Public License v. 2.0 for more details.
 *
 * For more details, see http://www.mrtrix.org/.
 */
 
#ifndef __image_filter_normalise3D_h__
#define __image_filter_normalise3D_h__
 
#include "image.h"
#include "algo/threaded_copy.h"
#include "adapter/normalise3D.h"
#include "filter/base.h"
 
namespace MR
{
  namespace Filter
  {
    class Normalise : public Base { MEMALIGN(Normalise)
 
      public:
        template <class HeaderType>
        Normalise (const HeaderType& in) :
            Base (in),
            extent (1,3) {
          datatype() = DataType::Float32;
        }
 
        template <class HeaderType>
          Normalise (const HeaderType& in, const std::string& message) :
            Base (in, message),
            extent (1,3) {
              datatype() = DataType::Float32;
            }
 
        template <class HeaderType>
        Normalise (const HeaderType& in, const vector<uint32_t>& extent) :
            Base (in),
            extent (extent) {
          datatype() = DataType::Float32;
        }
 
        template <class HeaderType>
          Normalise (const HeaderType& in, const std::string& message, const vector<uint32_t>& extent) :
            Base (in, message),
            extent (extent) {
              datatype() = DataType::Float32;
            }
 
        void set_extent (const vector<uint32_t>& ext) {
          for (size_t i = 0; i < ext.size(); ++i) {
            if (!(ext[i] & uint32_t(1)))
              throw Exception ("expected odd number for extent");
          }
          extent = ext;
        }
 
        template <class InputImageType, class OutputImageType>
        void operator() (InputImageType& in, OutputImageType& out) {
          Adapter::Normalise3D<InputImageType> normalise (in, extent);
          if (message.size())
            threaded_copy_with_progress_message (message, normalise, out);
          else
            threaded_copy (normalise, out);
        }
 
    protected:
        vector<uint32_t> extent;
    };
  }
}
 
 
#endif