ExampleBSplineCurve.cpp

/*
This file is part of the BIAS library (Basic ImageAlgorithmS).

Copyright (C) 2003-2009, 2005    (see file CONTACT for details)
  Multimediale Systeme der Informationsverarbeitung
  Institut fuer Informatik
  Christian-Albrechts-Universitaet Kiel


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

BIAS is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU Lesser General Public License for more details.

You should have received a copy of the GNU Lesser General Public License
along with BIAS; if not, write to the Free Software
Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
*/

/** @example ExampleBSplineCurve.cpp
    @relates BSplineCurve
    @ingroup g_examples    
    @brief   Example draws:
   - regular open quadratic bspline into ExampleBSplineCurve0.mip
   - regular closed quadratic bspline into ExampleBSplineCurve1.mip
   - regular open cubic bspline into ExampleBSplineCurve2.mip
   - regular closed cubic bspline into ExampleBSplineCurve3.mip
   @author MIP
*/


#include <Base/ImageUtils/ImageDraw.hh>
#include <Base/Image/ImageIO.hh>
#include <Base/Math/Matrix.hh>
#include <Base/Math/Vector.hh>
#include <Base/Math/Vector2.hh>
#include <MathAlgo/BSplineCurve.hh>

using namespace BIAS;
using namespace std;


BIAS::Image<unsigned char> img;
unsigned int imgCounter=0;
unsigned char value[1];


void drawQuadratic(BIAS::BSplineCurve bspline){
  img.FillImageWithConstValue(value);
  double t=0;
  double d=0.01;//stepsize
  Vector2<double> tmp;//remember last result
  Vector2<double> res;

  bspline.QuadraticUniformBSpline(res,t);//first result
  tmp=res;
  t+=d;
  while(t<1.+d){
    bspline.QuadraticUniformBSpline(res,t);
    ImageDraw<unsigned char>::Line(img,(unsigned int)tmp[0],
                                   (unsigned int)tmp[1],
                                   (unsigned int)res[0],
                                   (unsigned int)res[1],
                                   255);
    tmp=res;
    t+=d;
  }
  char title[255];
  sprintf(title, "ExampleBSplineCurve%01i.mip",imgCounter++);

  BIAS::ImageIO::Save(title, img);
}

void drawCubic(BIAS::BSplineCurve bspline){
  img.FillImageWithConstValue(value);
  double t=0;
  double d=0.01;//stepsize
  Vector2<double> tmp;//remember last result
  Vector2<double> res;

  bspline.CubicUniformBSpline(res,t);//first result
  tmp=res;
  t+=d;
  while(t<1.+d){
    bspline.CubicUniformBSpline(res,t);
    ImageDraw<unsigned char>::Line(img,(unsigned int)tmp[0],
                                   (unsigned int)tmp[1],
                                   (unsigned int)res[0],
                                   (unsigned int)res[1],
                                   255);
    tmp=res;
    t+=d;
  }
  char title[255];
  sprintf(title, "ExampleBSplineCurve%01i.mip",imgCounter++);

  BIAS::ImageIO::Save(title, img);
}


//int main(int args, char * argv[])
int main()
{
  //prepage image for dawing
  value[0]=0;
  img.Init(800,800);
  img.FillImageWithConstValue(value);

  //bspline constructor
  BSplineCurve bsplineQuadraticOpen,bsplineQuadraticClosed,
    bsplineCubicOpen,bsplineCubicClosed;

  //set controlpoints for simple open bspline (cubic and quadratic)
  std::vector<BIAS::Vector2<double> > cPnt;
  BIAS::Vector2<double> item;
  item[0]=0.;
  item[1]=0.;
  cPnt.push_back(item);
  item[0]=1345.996845;
  item[1]=1263.704834;
  cPnt.push_back(item);
  item[0]=1911.155259;
  item[1]=1785.161187;
  cPnt.push_back(item);
  item[0]=2385.32832;
  item[1]=2306.61754;
  cPnt.push_back(item);

  item[0]=2850.048898;
  item[1]=2828.073893;
  cPnt.push_back(item);

  item[0]=3375.122;
  item[1]=3349.530246;
  cPnt.push_back(item);

  item[0]=3934.636532;
  item[1]=3870.986599;
  cPnt.push_back(item);

  item[0]=4446.42884;
  item[1]=4392.442952;
  cPnt.push_back(item);

  item[0]=4990.218638;
  item[1]=4913.899306;
  cPnt.push_back(item);

  item[0]=5604.887952;
  item[1]=5435.355659;
  cPnt.push_back(item);

  item[0]=6120.227702;
  item[1]=5956.812012;
  cPnt.push_back(item);

  item[0]=8000;
  item[1]=8000;
  cPnt.push_back(item);


  for(unsigned int i=0;i<cPnt.size();i++)
    cPnt[i].DivideIP(10);

//   item[0]=200.;
//   item[1]=150.;
//   cPnt.push_back(item);
//   item[0]=100.;
//   item[1]=100.;
//   cPnt.push_back(item);
//   item[0]=100.;
//   item[1]=200.;
//   cPnt.push_back(item);
//   item[0]=200.;
//   item[1]=300.;
//   cPnt.push_back(item);
//   item[0]=200.;
//   item[1]=300.;
//   cPnt.push_back(item);
//   item[0]=300.;
//   item[1]=200.;
//   cPnt.push_back(item);
//   item[0]=300.;
//   item[1]=100.;
//   cPnt.push_back(item);
//   item[0]=200.;
//   item[1]=150.;
//   cPnt.push_back(item);
//   item[0]=200.;
//   item[1]=200.;
//   cPnt.push_back(item);


  bsplineQuadraticOpen.SetControlPoints(cPnt);
  bsplineCubicOpen.SetControlPoints(cPnt);


  //set controlpoints for simple closed bspline
  //cPnt[n],cPnt[0],...,cPnt[n],cPnt[0] //quadratic case
  //cPnt[n],cPnt[0],...,cPnt[n],cPnt[0],cPnt[1] //cubic case
  cPnt.clear();
  item[0]=200.;
  item[1]=200.;
  cPnt.push_back(item);
  item[0]=200.;
  item[1]=150.;
  cPnt.push_back(item);
  item[0]=100.;
  item[1]=100.;
  cPnt.push_back(item);
  item[0]=100.;
  item[1]=200.;
  cPnt.push_back(item);
  item[0]=200.;
  item[1]=300.;
  cPnt.push_back(item);
  item[0]=200.;
  item[1]=300.;
  cPnt.push_back(item);
  item[0]=300.;
  item[1]=200.;
  cPnt.push_back(item);
  item[0]=300.;
  item[1]=100.;
  cPnt.push_back(item);
  item[0]=200.;
  item[1]=150.;
  cPnt.push_back(item);
  item[0]=200.;
  item[1]=200.;
  cPnt.push_back(item);
  item[0]=200.;
  item[1]=150.;
  cPnt.push_back(item);
  bsplineQuadraticClosed.SetControlPoints(cPnt);
  item[0]=100.;
  item[1]=100.;
  cPnt.push_back(item);
  bsplineCubicClosed.SetControlPoints(cPnt);

  //init bsplines (compute coefficient matrices for segments)
  bsplineQuadraticOpen.InitQuadraticUniformBSpline();
  bsplineQuadraticClosed.InitQuadraticUniformBSpline();
  bsplineCubicOpen.InitCubicUniformBSpline();
  bsplineCubicClosed.InitCubicUniformBSpline();


  //draw open quadratic bspline into ExampleBSplineCurve0.mip
  drawQuadratic(bsplineQuadraticOpen);

  //draw closed quadratic bspline into ExampleBSplineCurve1.mip
  drawQuadratic(bsplineQuadraticClosed);

  //draw open cubic bspline into ExampleBSplineCurve2.mip
  drawCubic(bsplineCubicOpen);

  //draw closed cubic bspline into ExampleBSplineCurve3.mip
  drawCubic(bsplineCubicClosed);
  return 0;
}