Linux Cross Reference
Tina6/tina-tools/tinatool/draw/drawTv_proj2.c

   /*********
    * Copyright (c) 2003, Division of Imaging Science and Biomedical Engineering,
    * University of Manchester, UK.  All rights reserved.
    * 
    * Redistribution and use in source and binary forms, with or without modification, 
    * are permitted provided that the following conditions are met:
    * 
    *   . Redistributions of source code must retain the above copyright notice, 
    *     this list of conditions and the following disclaimer.
   *    
   *   . Redistributions in binary form must reproduce the above copyright notice,
   *     this list of conditions and the following disclaimer in the documentation 
   *     and/or other materials provided with the distribution.
   * 
   *   . Neither the name of the University of Manchester nor the names of its
   *     contributors may be used to endorse or promote products derived from this 
   *     software without specific prior written permission.
   * 
   * 
   * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" 
   * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 
   * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 
   * ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE 
   * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 
   * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF 
   * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
   * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN 
   * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) 
   * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 
   * POSSIBILITY OF SUCH DAMAGE.
   *
   **********
   *
   * Program :    TINA
   * File    :  $Source: /home/tina/cvs/tina-tools/tinatool/draw/drawTv_proj2.c,v $
   * Date    :  $Date: 2003/10/01 16:02:47 $
   * Version :  $Revision: 1.2 $
   * CVS Id  :  $Id: drawTv_proj2.c,v 1.2 2003/10/01 16:02:47 tony Exp $
   *
   * Author  : Legacy TINA
   *
   * Notes :
   *
   *********
  */
  
  #include "drawTv_proj2.h"
  
  #if HAVE_CONFIG_H
    #include <config.h>
  #endif
  
  #include <stdio.h>
  #include <tina/sys/sysDef.h>
  #include <tina/math/mathDef.h>
  #include <tina/math/mathPro.h>
  #include <tinatool/draw/draw_TvDef.h>
  
  /**2D general projection**/
  
  static Ipos proj2(Tv * tv, Vec2 v)
  {
      Ipos    pos = {Ipos_id};
      double  x, y;
  
      v = vec2_diff(v, tv->centre2);
      x = tv->cx + tv->scalex * vec2_dot(v, tv->ex2);
      y = tv->cy + tv->scaley * vec2_dot(v, tv->ey2);
      pos.x = ROUND(x);
      pos.y = ROUND(y);
      return (pos);
  }
  
  static Vec2 backproj2(Tv * tv, Ipos pos)
  {
      Vec2    p = {Vec2_id};
      double  x = (pos.x - tv->cx) / tv->scalex;
      double  y = (pos.y - tv->cy) / tv->scaley;
  
      p = vec2_sum3(tv->centre2, vec2_times(x, tv->ex2), vec2_times(y, tv->ey2));
      return (p);
  }
  
  void    tv_set_proj2(Tv * tv)
  {
      tv->proj2 = proj2;
      tv->backproj2 = backproj2;
  }
  
  /**generic 2D projection**/
  
  Ipos    tv_proj2(Tv * tv, Vec2 p)
  {
      return (tv->proj2(tv, p));
  }
  
  Vec2    tv_backproj2(Tv * tv, Ipos pos)
  {
      return (tv->backproj2(tv, pos));
 }
 
 Bool    tv_ipos_within(Tv * tv, Ipos pos)
 {
     int     r, c;
 
     if (tv == NULL)
         return (false);
 
     r = ipos_y(pos);
     c = ipos_x(pos);
 
     return (Bool) (r >= 0 && r < tv->height && c >= 0 && c < tv->width);
 }
 
 /**setting and changing 2D cameras**/
 
 void    tv_camera2_rect(Tv * tv, Vec2 v1, Vec2 v2)
 {
     double  x1 = vec2_x(v1);
     double  y1 = vec2_y(v1);
     double  x2 = vec2_x(v2);
     double  y2 = vec2_y(v2);
     double  width, height;
 
     if (tv == NULL || x1 >= x2 || y1 >= y2)
         return;
     width = x2 - x1;
     height = y2 - y1;
     tv->scalex = (float)(tv->width / width);
     tv->scaley = (float)(tv->height / height);
     tv->cx = (float)(tv->width / 2);
     tv->cy = (float)(tv->height / 2);
     tv->centre2 = vec2((x1 + x2) / 2.0, (y1 + y2) / 2.0);
     tv->radius2 = (float)(MAX(width, height));
     tv->ex2 = vec2_ex();
     tv->ey2 = vec2_minus(vec2_ey());
 }
 
 void    tv_camera2_roi(Tv * tv, Imregion * roi)
 {
     double  x1 = roi->lx;
     double  y1 = roi->ly;
     double  x2 = roi->ux;
     double  y2 = roi->uy;
     double  width, height;
 
     if(tv == NULL)
         return;
     if (x1 >= x2 || y1 >= y2)
         return;
     width = x2 - x1;
     height = y2 - y1;
     tv->scalex = tv->scaley = (float)(MIN(tv->width / width, tv->height / height));
     tv->cx = (float)(tv->width / 2);
     tv->cy = (float)(tv->height / 2);
     tv->centre2 = vec2((x1 + x2) / 2.0, (y1 + y2) / 2.0);
     tv->radius2 = (float)(MAX(width, height));
     tv->ex2 = vec2_ex();
     tv->ey2 = vec2_ey();
 }
 
 void    tv_camera2_image(Tv * tv, int width, int height)
 {
     if(tv == NULL)
         return;
     tv->scalex = tv->scaley =
     (float)(MIN(tv->width / (double) width, tv->height / (double) height));
     tv->cx = (float)(tv->width / 2);
     tv->cy = (float)(tv->height / 2);
     tv->centre2 = vec2(width / 2.0, height / 2.0);
     tv->radius2 = (float)(MAX(width, height));
     tv->ex2 = vec2_ex();
     tv->ey2 = vec2_ey();
 }
 
 /* Stretch the camera of a tv to make region of interest of an image
  * (with given aspect ratio) fill the screen. */
 void    tv_camera2_roi_stretch(Tv * tv, Imregion * roi, double aspect_ratio)
 {
     double  width = (double) (roi->ux - roi->lx);
     double  cx = (double) (roi->ux + roi->lx) / 2.0;
     double  cy = (double) (roi->uy + roi->ly) / 2.0;
 
     if(tv == NULL)
         return;
     tv->scalex = (float)(tv->width / width);
     tv->scaley = (float)(tv->scalex * aspect_ratio);
     tv->cx = (float)(tv->width / 2);
     tv->cy = (float)(tv->height / 2);
     tv->centre2 = vec2(cx, cy);
     tv->ex2 = vec2_ex();
     tv->ey2 = vec2_ey();
 }
 
 void    tv_camera2(Tv * tv, Vec2 centre, double radius, Vec2 down)
 {
     if(tv == NULL)
         return;
     tv->scalex = tv->scaley = (float)(0.5 * MIN(tv->width, tv->height) / radius);
     tv->cx = (float)(tv->width / 2);
     tv->cy = (float)(tv->height / 2);
     tv->centre2 = centre;
     tv->radius2 = (float)radius;
     vec2_basis(down, &tv->ex2, &tv->ey2);
 }
 
 void    tv_orient2(Tv * tv, Vec2 down)
 {
     if(tv == NULL)
         return;
     vec2_basis(down, &tv->ex2, &tv->ey2);
 }