grid.h

/***************************************************************************
                          grid.h  -  description
                             -------------------
    begin                : Wed Feb 21 2007
    copyright            : (C) 2001 by Thies Jochimsen
    email                : jochimse@cns.mpg.de
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/

#ifndef RECOGRID_H
#define RECOGRID_H

#include "step.h"



template<int NDim>
class RecoGrid : public RecoStep {

  // implementing virtual functions of RecoStep
  STD_string label() const {return "grid"+postfix();}
  STD_string description() const {return "Grid signal data with arbitrary k-space trajectory in "+itos(NDim)+" dimensions";}
  bool process(RecoData& rd, RecoController& controller);
  RecoCoord input_coord() const {return RecoCoord::coord_with_mode(RecoIndex::collected, readout);}
  void modify_coord(RecoCoord& coord) const {if(NDim==2) coord.set_mode(RecoIndex::collected, line);}
  bool query(RecoQueryContext& context);
  RecoStep* allocate() const {return new RecoGrid<NDim>;}
  void init();

  JDXbool   deapo_shift;
  JDXfilter kernel;
  JDXfloat  kernelwidth;
  JDXfloat  gridScale;

  static STD_string postfix() {return itos(NDim)+"d";}

  static TinyVector<int,NDim> gridshape(const RecoCoord& coord, RecoController& controller);

  const Gridding<STD_complex,NDim>& get_gridder(const RecoCoord& coord, RecoController& controller);

  // cache gridders
  typedef STD_map<const KspaceTraj*, Gridding<STD_complex,NDim> > GridderMap;
  GridderMap gridmap;
  Mutex gridmutex;

};




template<int NDim>
class RecoDeapodize : public RecoStep {

  // implementing virtual functions of RecoStep
  STD_string label() const {return "deapodize"+postfix();}
  STD_string description() const {return "De-apodize gridded data in "+itos(NDim)+" dimensions";}
  bool process(RecoData& rd, RecoController& controller);
  RecoCoord input_coord() const {return RecoCoord::coord_with_mode(RecoIndex::collected, line3d, line, readout);}
  void modify_coord(RecoCoord& coord) const {}
  RecoStep* allocate() const {return new RecoDeapodize;}
  void init() {}

  static STD_string postfix() {return itos(NDim)+"d";}

};

class RecoGridCut : public RecoStep {

  // implementing virtual functions of RecoStep
  STD_string label() const {return "gridcut"; }
  STD_string description() const {return "Cut out central part of image to remove effects from gridding on a denser grid";}
  bool process(RecoData& rd, RecoController& controller);
  RecoCoord input_coord() const {return RecoCoord::coord_with_mode(RecoIndex::collected, line3d, line, readout);}
  void modify_coord(RecoCoord& coord) const {}
  RecoStep* allocate() const {return new RecoGridCut;}
  void init() {}
};


#endif

---