CylindersWithLineSource.h

/*
    Copyright (C) 2000- 2002, IRSL
    See STIR/LICENSE.txt for details
*/

#include "stir/Shape/EllipsoidalCylinder.h"
#include "stir/Shape/CombinedShape3D.h"

START_NAMESPACE_STIR

class LineSource_phantom
{
public:
  inline LineSource_phantom();

  inline const shared_ptr<Shape3D>& get_A_ptr() const { return A_ptr; }

  inline shared_ptr<Shape3D> make_union_ptr(const float fraction) const;

  inline void translate(const CartesianCoordinate3D<float>& direction);
  inline void scale(const CartesianCoordinate3D<float>& scale3D);

private:
  shared_ptr<Shape3D> A_ptr;
};

LineSource_phantom::LineSource_phantom()
{
  A_ptr = new EllipsoidalCylinder(150, 0, 0, CartesianCoordinate3D<float>(0, 0, 0), 0, 0, 0);
}

void
LineSource_phantom::translate(const CartesianCoordinate3D<float>& direction)
{
  A_ptr->translate(direction);
}

void
LineSource_phantom::scale(const CartesianCoordinate3D<float>& scale3D)
{
  // check this!!!
  A_ptr->scale_around_origin(scale3D);
}

shared_ptr<Shape3D>
LineSource_phantom::make_union_ptr(const float fraction) const
{
  shared_ptr<Shape3D> full_A = get_A_ptr()->clone();

  return full_A;
}

class CylindersWithLineSource_phantom
{
public:
  inline CylindersWithLineSource_phantom();

  // inline const shared_ptr<Shape3D>& get_A_ptr() const
  //{ return A_ptr; }

  inline const shared_ptr<Shape3D>& get_B_ptr() const { return B_ptr; }

  inline const shared_ptr<Shape3D>& get_C_ptr() const { return C_ptr; }

  // inline const shared_ptr<Shape3D>& get_B1_ptr() const
  //{ return B1_ptr; }

  // inline const shared_ptr<Shape3D>& get_C1_ptr() const
  //{ return C1_ptr; }

  inline shared_ptr<Shape3D> make_union_ptr(const float fraction) const;

  inline void translate(const CartesianCoordinate3D<float>& direction);
  inline void scale(const CartesianCoordinate3D<float>& scale3D);

private:
  // shared_ptr<Shape3D> A_ptr;
  shared_ptr<Shape3D> B_ptr;
  shared_ptr<Shape3D> C_ptr;
  // shared_ptr<Shape3D> B1_ptr;
  // shared_ptr<Shape3D> C1_ptr;
};

CylindersWithLineSource_phantom::CylindersWithLineSource_phantom()
{
  // A_ptr = new EllipsoidalCylinder (150,19,19,
  // CartesianCoordinate3D<float>(0,0,0),
  // 0,0,0);

  // for 966
  // B_ptr = new EllipsoidalCylinder (500,60,60,
  // CartesianCoordinate3D<float>(0,0,0),
  //  0,0,0);

  // two planes only
  /* B_ptr = new EllipsoidalCylinder (1,6,6,
     CartesianCoordinate3D<float>(0,-14,0),
     0,0,0);

   C_ptr = new EllipsoidalCylinder (1,6,6,
     CartesianCoordinate3D<float>(0,14,0),
     0,0,0);*/
  // off centre

  // this is to check the bloody arthifacts in new filters

  /*B_ptr = new EllipsoidalCylinder (1,6,6,
    CartesianCoordinate3D<float>(0,-13,4),
    0,0,0);

  C_ptr = new EllipsoidalCylinder (1,6,6,
    CartesianCoordinate3D<float>(0,13,4),
    0,0,0);*/

  // these are the ones noramly used
  B_ptr = new EllipsoidalCylinder(1, 8, 8, CartesianCoordinate3D<float>(0, -14, 5), 0, 0, 0);

  C_ptr = new EllipsoidalCylinder(1, 8, 8, CartesianCoordinate3D<float>(0, 14, 5), 0, 0, 0);

  // the one that I use for resolution
  /* B_ptr = new EllipsoidalCylinder (100,8,8,
     CartesianCoordinate3D<float>(0,-14,0),
     0,0,0);

   C_ptr = new EllipsoidalCylinder (100,8,8,
     CartesianCoordinate3D<float>(0,14,0),
     0,0,0);*/

  // the one that I use for resolution
  /* B1_ptr = new EllipsoidalCylinder (1,6,6,
    CartesianCoordinate3D<float>(0,0,-14),
    0,0,0);

  C1_ptr = new EllipsoidalCylinder (1,6,6,
    CartesianCoordinate3D<float>(0,0,14),
    0,0,0);*/
  // off centre

  /*B1_ptr = new EllipsoidalCylinder (1,6,6,
    CartesianCoordinate3D<float>(0,-2,-14),
   0,0,0);

  C1_ptr = new EllipsoidalCylinder (1,6,6,
    CartesianCoordinate3D<float>(0,-2,14),
    0,0,0);*/
}
void
CylindersWithLineSource_phantom::translate(const CartesianCoordinate3D<float>& direction)
{
  //  A_ptr->translate(direction);
  B_ptr->translate(direction);
  C_ptr->translate(direction);
  // B1_ptr->translate(direction);
  // C1_ptr->translate(direction);
}

void
CylindersWithLineSource_phantom::scale(const CartesianCoordinate3D<float>& scale3D)
{
  // check this!!!
  // A_ptr->scale_around_origin(scale3D);
  B_ptr->scale_around_origin(scale3D);
  C_ptr->scale_around_origin(scale3D);
  // B1_ptr->scale_around_origin(scale3D);
  // C1_ptr->scale_around_origin(scale3D);
}

shared_ptr<Shape3D>
CylindersWithLineSource_phantom::make_union_ptr(const float fraction) const
{
  // shared_ptr<Shape3D> full_A = get_A_ptr()->clone();
  shared_ptr<Shape3D> full_B = get_B_ptr()->clone();
  shared_ptr<Shape3D> full_C = get_C_ptr()->clone();

  // shared_ptr<Shape3D> full_B1 = get_B1_ptr()->clone();
  // shared_ptr<Shape3D> full_C1 = get_C1_ptr()->clone();

  shared_ptr<Shape3D> AB_union = new CombinedShape3D<logical_and<bool>>(full_C, full_B);

  // shared_ptr<Shape3D>  AB1_union =
  // new CombinedShape3D <logical_and<bool> >( full_C1,full_B1);

  // shared_ptr<Shape3D>  AB1_AB_union =
  // new CombinedShape3D <logical_and<bool> >( AB_union, AB1_union);

  //  shared_ptr<Shape3D>  ABC_union =
  // new CombinedShape3D <logical_and<bool> >( AB_union,full_C);

  // return AB1_AB_union;
  return AB_union;
  // return full_B;
}

// OLD
#if 0

class CylindersWithLineSource_phantom
{
public:


  inline CylindersWithLineSource_phantom();
 
  inline const shared_ptr<Shape3D>& get_A_ptr() const
  { return A_ptr; }
  
  //inline const shared_ptr<Shape3D>& get_B_ptr() const
 // { return B_ptr; }

  /* inline const shared_ptr<Shape3D>& get_C_ptr() const
  { return C_ptr; }

   inline const shared_ptr<Shape3D>& get_D_ptr() const
  { return D_ptr; }*/

 inline shared_ptr<Shape3D> make_union_ptr(const float fraction) const;
   
  inline void translate(const CartesianCoordinate3D<float>& direction);
  inline void scale(const CartesianCoordinate3D<float>& scale3D);

private: 
  shared_ptr<Shape3D> A_ptr;
 // shared_ptr<Shape3D> B_ptr;
 // shared_ptr<Shape3D> C_ptr;
 //shared_ptr<Shape3D> D_ptr;
  


};


CylindersWithLineSource_phantom::CylindersWithLineSource_phantom()
{

  A_ptr = new EllipsoidalCylinder (100,1,1,
    CartesianCoordinate3D<float>(0,0,0),
    0,0,0);


 //  A_ptr = new EllipsoidalCylinder (170,10,10,
   // CartesianCoordinate3D<float>(0,15,-10),
  //  0,0,0);

 // B_ptr = new EllipsoidalCylinder (100,10,10,
   // CartesianCoordinate3D<float>(0,-15,-10),
   // 0,0,0);

  /* A_ptr = new EllipsoidalCylinder (100,90,90,
    CartesianCoordinate3D<float>(0,15,-25),
    0,0,0);

  B_ptr = new EllipsoidalCylinder (150,1,1,
    CartesianCoordinate3D<float>(0,0,0),
    0,0,0);*/


 /* C_ptr = new EllipsoidalCylinder (170,1,1,
    CartesianCoordinate3D<float>(0,0,-25),
    0,0,0);*/

 /* D_ptr = new EllipsoidalCylinder (170,1,1,
    CartesianCoordinate3D<float>(0,0,20),
    0,0,0);*/


}


void 
CylindersWithLineSource_phantom::translate(const CartesianCoordinate3D<float>& direction)
{
  
  A_ptr->translate(direction);
 // B_ptr->translate(direction);
 // C_ptr->translate(direction);
//  D_ptr->translate(direction);

}

void 
CylindersWithLineSource_phantom::scale(const CartesianCoordinate3D<float>& scale3D)
{
  //check this!!!  
  A_ptr->scale_around_origin(scale3D);
 // B_ptr->scale_around_origin(scale3D);
  //C_ptr->scale_around_origin(scale3D);
 // D_ptr->scale_around_origin(scale3D);

}

shared_ptr<Shape3D>
CylindersWithLineSource_phantom::make_union_ptr(const float fraction) const
{
   shared_ptr<Shape3D> full_A = get_A_ptr()->clone();
  // shared_ptr<Shape3D> full_B = get_B_ptr()->clone();
   //shared_ptr<Shape3D> full_C = get_C_ptr()->clone();
  // shared_ptr<Shape3D> full_D = get_D_ptr()->clone();
   
   //shared_ptr<Shape3D>  AB_union =
   
   return full_A ;
     
  // new CombinedShape3D <logical_and<bool> >( full_A,full_B);

  
  
   /* shared_ptr<Shape3D>  ABC_union = 
     new CombinedShape3D <logical_and<bool> >( AB_union,full_C);*/

   // return 
  //  new CombinedShape3D< logical_and<bool> >( AB_union,full_D);

}

#endif
END_NAMESPACE_STIR