sample.h

/***************************************************************************
                          sample.h  -  description
                             -------------------
    begin                : Mon Jul 11 2005
    copyright            : (C) 2001 by Thies H. 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 SAMPLE_H
#define SAMPLE_H

#include <odinpara/jdxblock.h>
#include <odinpara/jdxnumbers.h>
#include <odinpara/jdxarrays.h>
#include <odinpara/odinpara.h>

enum sampleDim { frameDim=0, freqDim, zDim, yDim, xDim, n_sampleDim };




class Sample : public JcampDxBlock {

 public:

  Sample(const STD_string& label="unnamedSample", bool uniformFOV=true, bool uniformT1T2=false);

  Sample(const Sample& ss);

  Sample& operator = (const Sample& ss);

  Sample& set_FOV(float fov);

  Sample& set_FOV(axis direction, float fov);

  float get_FOV(axis direction) const;

  Sample& set_spatial_offset(axis direction, float offs) {offset[direction]=offs; return *this;}

  float get_spatial_offset(axis direction) const {return offset[direction];}

  Sample& set_freqrange(float range) {freqrange=range; return *this;}

  float get_freqrange() const {return freqrange;}

  Sample& set_freqoffset(float offs) {freqoffset=offs; return *this;}

  float get_freqoffset() const {return freqoffset;}


  Sample& set_frame_durations(const dvector& intervals) {frameDurations=intervals; return *this;}

  const dvector& get_frame_durations() const {return frameDurations;}


  Sample& resize(unsigned int framesize, unsigned int freqsize, unsigned int zsize, unsigned int ysize, unsigned int xsize);

  const ndim& get_extent() const {return spinDensity.get_extent();}

  Sample& set_T1(float relaxation_time) {T1=relaxation_time; haveT1map=false; return *this;}

  Sample& set_T2(float relaxation_time) {T2=relaxation_time; haveT2map=false; return *this;}

  Sample& set_T1map(const farray& t1map);

  const farray& get_T1map() const;

  Sample& set_T2map(const farray& t2map);

  const farray& get_T2map() const;

  Sample& set_ppmMap(const farray& ppmmap);

  const farray& get_ppmMap() const;

  Sample& set_spinDensity(const farray& sd);

  const farray& get_spinDensity() const;

  Sample& set_DcoeffMap(const farray& dmap);

  const farray& get_DcoeffMap() const;

  Sample& update();


  // overwriting virtual functions from JcampDxClass
  int load(const STD_string& filename);


 private:
  friend class SeqSimMagsi;

  int append_all_members();

  bool check_and_correct(const char* mapname, const farray& srcmap, farray& dstmap);


  JDXfloat FOVall;
  JDXtriple FOV;
  bool uniFOV;
  JDXtriple offset;

  JDXfloat freqrange;
  JDXfloat freqoffset;

  JDXdoubleArr frameDurations;

  mutable JDXfloatArr spinDensity;

  bool uniT1T2;
  JDXfloat T1;
  JDXfloat T2;
  mutable JDXfloatArr T1map;
  mutable JDXfloatArr T2map;
  mutable bool haveT1map;
  mutable bool haveT2map;

  mutable JDXfloatArr ppmMap;
  mutable bool have_ppmMap;

  mutable JDXfloatArr DcoeffMap;
  mutable bool have_DcoeffMap;

};


#endif

---