rtkSheppLoganPhantom.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright RTK Consortium
 *
 *  Licensed under the Apache License, Version 2.0 (the "License");
 *  you may not use this file except in compliance with the License.
 *  You may obtain a copy of the License at
 *
 *         https://www.apache.org/licenses/LICENSE-2.0.txt
 *
 *  Unless required by applicable law or agreed to in writing, software
 *  distributed under the License is distributed on an "AS IS" BASIS,
 *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 *  See the License for the specific language governing permissions and
 *  limitations under the License.
 *
 *=========================================================================*/

#ifndef rtkSheppLoganPhantom_h
#define rtkSheppLoganPhantom_h

#include "RTKExport.h"
#include "rtkGeometricPhantom.h"
#include "rtkQuadricShape.h"

namespace rtk
{
class RTK_EXPORT SheppLoganPhantom : public GeometricPhantom
{
public:
  ITK_DISALLOW_COPY_AND_MOVE(SheppLoganPhantom);

  using Self = SheppLoganPhantom;
  using Superclass = itk::DataObject;
  using Pointer = itk::SmartPointer<Self>;
  using ConstPointer = itk::SmartPointer<const Self>;

  using ScalarType = QuadricShape::ScalarType;
  using PointType = QuadricShape::PointType;
  using VectorType = QuadricShape::VectorType;

  itkNewMacro(Self);

#ifdef itkOverrideGetNameOfClassMacro
  itkOverrideGetNameOfClassMacro(SheppLoganPhantom);
#else
  itkTypeMacro(SheppLoganPhantom, GeometricPhantom);
#endif


protected:
  SheppLoganPhantom();
  ~SheppLoganPhantom() override = default;

private:
  void
  SetEllipsoid(ScalarType spax,
               ScalarType spay,
               ScalarType spaz,
               ScalarType centerx,
               ScalarType centery,
               ScalarType centerz,
               ScalarType angle,
               ScalarType density);
};

} // namespace rtk
#endif // rtkSheppLoganPhantom_h