sample.h

/***************************************************************************
                          sample.h  -  description
                             -------------------
    begin                : Mon Jul 11 2005
    copyright            : (C) 2000-2021 by Thies H. Jochimsen
    email                : thies@jochimsen.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/ldrblock.h>
#include <odinpara/ldrnumbers.h>
#include <odinpara/ldrarrays.h>
#include <odinpara/odinpara.h>
 
enum sampleDim { frameDim=0, freqDim, zDim, yDim, xDim, n_sampleDim };
 
 
 
 
class Sample : public LDRblock {
 
 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 LDRbase
  int load(const STD_string& filename, const LDRserBase& serializer=LDRserJDX());
 
 
 private:
  friend class SeqSimMagsi;
 
  int append_all_members();
 
  bool check_and_correct(const char* mapname, const farray& srcmap, farray& dstmap);
 
 
  LDRfloat FOVall;
  LDRtriple FOV;
  bool uniFOV;
  LDRtriple offset;
 
  LDRfloat freqrange;
  LDRfloat freqoffset;
 
  LDRdoubleArr frameDurations;
 
  mutable LDRfloatArr spinDensity;
 
  bool uniT1T2;
  LDRfloat T1;
  LDRfloat T2;
  mutable LDRfloatArr T1map;
  mutable LDRfloatArr T2map;
  mutable bool haveT1map;
  mutable bool haveT2map;
 
  mutable LDRfloatArr ppmMap;
  mutable bool have_ppmMap;
 
  mutable LDRfloatArr DcoeffMap;
  mutable bool have_DcoeffMap;
 
};
 
 
#endif