vuML.h

// -*- C++ -*-
#ifndef __MARSCHNERLOBB_H__
#define __MARSCHNERLOBB_H__

// written by Reza Entezari
// adapted to vuVolume by Steven Bergner
// adapted to 775 final proj by Tai Meng

#include <vuVector.h>

#define MAX_INTENSITY 254 //vuVolume doesn't handle 255 properly

class MarschnerLobb {
public:
  MarschnerLobb();
  ~MarschnerLobb();
  
  void initialize(double x, double y, double z);  

  double getValue(const vuVector& p) const;


  double getValueGaussianNoise(const vuVector& p, const double& stdDev) const;

  /*
    Set value to 0 (pepper) with probability p/2.
    Set value to 1 (salt) with probability p/2.
  */
  double getValueSaltPepperNoise(const vuVector&p, const double& percentage) const;

private:
  double xsize, ysize, zsize;
  double fm, alpha;
};

#endif

---