runtime.cpp

Go to the documentation of this file.

#include <qfiledialog.h>
#include <qimage.h>
#include <typeinfo>
#include <math.h>
#include "runtime.h"
#include "image-dialog.h"
#include "text.h"
#include "line.h"
#include "box.h"
#include "status.h"

/*
Value * Runtime::distributionxz(ImageValue * image)
{
  int as=image->width();
  int bs=256;
  int cs=image->channelcount();
  ImageValue * O = new ImageValue(as,bs,cs);
  int xs = image->width();
  int ys = image->height();
  O->clear(0);
  double zmax = image->maxz();
  double zmin = image->minz();
  double ztot = zmax-zmin;
  for(int c = 0; c < cs ; c ++)
    for(int x = 0; x < xs ; x ++)
      for(int y = 0 ; y < ys ; y++)
        {
          double z = image->get(x,y,c);
          z-=zmin;
          z/=ztot;
          z*=255.0;
          O->inc(x,(int)z,c);
        }
  return O;
}

Value * Runtime::distributionyz(ImageValue * image)
{
  int as=256;
  int bs=image->height();
  int cs=image->channelcount();
  ImageValue * O = new ImageValue(as,bs,cs);
  int xs = image->width();
  int ys = image->height();
  O->clear(0);
  double zmax = image->maxz();
  double zmin = image->minz();
  double ztot = zmax-zmin;
  for(int c = 0; c < cs ; c ++)
    for(int x = 0; x < xs ; x ++)
      for(int y = 0 ; y < ys ; y++)
        {
          double z = image->get(x,y,c);
          z-=zmin;
          z/=ztot;
          z*=255.0;
          O->inc((int)z,y,c);
        }
  return O;
}

Value * Runtime::distribution(ImageValue * image)
{
  if (image->width() > image->height()) 
    return distributionxz(image);
  return distributionyz(image);
}

Value * Runtime::distribution(vector<Value*> v)
{
  if (v.size()==0) return NULL;
  Value * v0 = v[0];
  if (v0==NULL) return NULL;
  if (typeid(*v0)!=typeid(ImageValue)) return NULL;
  ImageValue *I=(ImageValue*)v0;
  if (v.size()==1) return distribution(I);
  return NULL;
}


Value * Runtime::projectxz(ImageValue * image)
{
  int as=image->width();
  int bs=256;
  int cs=image->channelcount();
  ImageValue * O = new ImageValue(as,bs,cs);
  int xs = image->width();
  int ys = image->height();
  O->clear(1);
  double zmax = image->maxz();
  double zmin = image->minz();
  double ztot = zmax-zmin;
  for(int c = 0; c < cs ; c ++)
    for(int x = 0; x < xs ; x ++)
      {
        double z = 0;
        for(int y = 0 ; y < ys ; y++)
          z+=image->get(x,y,c);
        z/=ys;
        z-=zmin;
        z/=ztot;
        z*=255;
        O->set(x,(int)z,c,0);
      }
  return O;
}

Value * Runtime::projectyz(ImageValue * image)
{
  int as=256;
  int bs=image->height();
  int cs=image->channelcount();
  ImageValue * O = new ImageValue(as,bs,cs);
  int xs = image->width();
  int ys = image->height();
  O->clear(1);
  double zmax = image->maxz();
  double zmin = image->minz();
  double ztot = zmax-zmin;
  for(int c = 0; c < cs ; c ++)
    for(int y = 0 ; y < ys ; y++)
      {
        double z = 0;
        for(int x = 0; x < xs ; x ++)
          z += image->get(x,y,c);
        z/=xs;
        z-=zmin;
        z/=ztot;
        z*=255.0;
        O->set((int)z,y,c,0);
      }
  return O;
}

Value * Runtime::project(ImageValue * image)
{
  if (image->width() > image->height()) 
    return projectxz(image);
  return projectyz(image);
}

Value * Runtime::project(vector<Value*> v)
{
  if (v.size()==0) return NULL;
  Value * v0 = v[0];
  if (v0==NULL) return NULL;
  if (typeid(*v0)!=typeid(ImageValue)) return NULL;
  ImageValue *I=(ImageValue*)v0;
  if (v.size()==1) return project(I);
  return NULL;
}

*/

Value * Runtime::apply(QString op, vector<Value*> evaluated)
{
  if (op=="+") return add(evaluated);
  if (op=="-") return sub(evaluated);
  if (op=="/") return div(evaluated);
  if (op=="*") return mul(evaluated);
  if (op=="image") return image(evaluated);
  // selection
  if (op=="select") return select(evaluated);
  if (op=="line") return line(evaluated);
  if (op=="box") return box(evaluated);
  // translation, rotation and affine transformation
  if (op=="center") return center(evaluated);
  if (op=="zoom") return zoom(evaluated);
  if (op=="zoomx") return zoomx(evaluated);
  if (op=="zoomy") return zoomy(evaluated);
  if (op=="straighten") return straighten(evaluated);
  if (op=="crop") return crop(evaluated);
  // projection like operations
  //  if (op=="distribution") return distribution(evaluated);
  //  if (op=="project") return project(evaluated);
#ifdef PSYTRACK
  if (op=="morph") return morph(evaluated);
#endif
  throw new Error("Unknown operator "+op);
}

bool Runtime::requires_user_input(Value * val)
{
  if (!val) return false;
  if (typeid(*val)!=typeid(Expression))
    return false;
  Expression * app = (Expression*)val;
  for(int i = 0 ; i < app->argc() ; i++)
    if (requires_user_input(app->argv(i))) return true;
  QString op = app->op();
  if (op=="image") return true;
  if (op=="select") return true;
  if (op=="line") return true;
  if (op=="box") return true;
  return false;
}

Runtime runtime;

---