drawTv_proj3.c

Go to the documentation of this file.

/*********
 * 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_proj3.c,v $
 * Date    :  $Date: 2003/10/01 16:02:47 $
 * Version :  $Revision: 1.2 $
 * CVS Id  :  $Id: drawTv_proj3.c,v 1.2 2003/10/01 16:02:47 tony Exp $
 *
 * Author  : Legacy TINA
 *
 * Notes :
 *
 *********
*/

#include "drawTv_proj3.h"

#if HAVE_CONFIG_H
  #include <config.h>
#endif

#include <math.h>
#include <float.h>
#include <limits.h>
#include <tina/sys/sysDef.h>
#include <tina/math/mathDef.h>
#include <tina/math/mathPro.h>
#include <tinatool/draw/draw_TvDef.h>

static Ipos proj3_orth(Tv * tv, Vec3 v)
{
    Ipos    pos = {Ipos_id};
    double  x, y;

    v = vec3_diff(v, tv->centre3);
    x = tv->cx + tv->scalex * vec3_dot(v, tv->ex3);
    y = tv->cy + tv->scaley * vec3_dot(v, tv->ey3);
    pos.x = ROUND(x);
    pos.y = ROUND(y);
    return (pos);
}

static void ray3_orth(Tv * tv, Ipos pos, Vec3 * p, Vec3 * v)
{
    double  x = (pos.x - tv->cx) / tv->scalex;
    double  y = (pos.y - tv->cy) / tv->scaley;

    *p = vec3_sum3(tv->centre3, vec3_times(x, tv->ex3), vec3_times(y, tv->ey3));
    *v = tv->ez3;
}

static Vec3 backproj3_orth(Tv * tv, Ipos pos)
{
    Vec3    p = {Vec3_id};
    double  x = (pos.x - tv->cx) / tv->scalex;
    double  y = (pos.y - tv->cy) / tv->scaley;

    p = vec3_sum3(tv->centre3, vec3_times(x, tv->ex3), vec3_times(y, tv->ey3));
    return (p);
}

void    tv_set_proj3_orth(Tv * tv)
{
    tv->proj3type = ORTH;
    tv->proj3 = proj3_orth;
    tv->ray3 = ray3_orth;
    tv->backproj3 = backproj3_orth;
}

static Ipos proj3_persp(Tv * tv, Vec3 v)
{
    Ipos    pos = {Ipos_id};
    double  x, y;

    v = vec3_unit(vec3_diff(v, tv->pcentre));
    v = vec3_inter_line_plane(tv->pcentre, v, tv->centre3, tv->ez3);
    v = vec3_diff(v, tv->centre3);
    x = tv->cx + tv->scalex * vec3_dot(v, tv->ex3);
    y = tv->cy + tv->scaley * vec3_dot(v, tv->ey3);
    pos.x = ROUND(x);
    pos.y = ROUND(y);
    return (pos);
}

static void ray3_persp(Tv * tv, Ipos pos, Vec3 * p, Vec3 * v)
{
    double  x = (pos.x - tv->cx) / tv->scalex;
    double  y = (pos.y - tv->cy) / tv->scaley;
    Vec3    q = {Vec3_id};

    q = vec3_sum3(tv->centre3, vec3_times(x, tv->ex3), vec3_times(y, tv->ey3));
    *p = tv->pcentre;
    *v = vec3_unit(vec3_diff(q, tv->pcentre));
}

static Vec3 backproj3_persp(Tv * tv, Ipos pos)
{
    Vec3    p = {Vec3_id};
    double  x = (pos.x - tv->cx) / tv->scalex;
    double  y = (pos.y - tv->cy) / tv->scaley;

    p = vec3_sum3(tv->centre3, vec3_times(x, tv->ex3), vec3_times(y, tv->ey3));
    return (p);
}

void    tv_set_proj3_persp(Tv * tv)
{
    tv->proj3type = PERSP;
    tv->proj3 = proj3_persp;
    tv->ray3 = ray3_persp;
    tv->backproj3 = backproj3_persp;
}

Ipos    tv_proj3(Tv * tv, Vec3 p)
{
    return (tv->proj3(tv, p));
}

void    tv_ray3(Tv * tv, Ipos pos, Vec3 * p, Vec3 * v)
{
    tv->ray3(tv, pos, p, v);
}

Vec3    tv_backproj3(Tv * tv, Ipos pos)
{
    return (tv->backproj3(tv, pos));
}

double  tv_dist3(Tv * tv, Vec3 p)
{
    if (tv == NULL)
        return (FLT_MAX);
    switch (tv->proj3type)
    {
    case ORTH:
        return (vec3_dot(vec3_diff(p, tv->centre3), tv->ez3));
    case PERSP:
        return (vec3_dist(p, tv->pcentre));
    }
    return (FLT_MAX);
}

Ipos    tv_zbuff_proj3(Tv * tv, Vec3 p, int *iz)
{
    double  z;

    z = (tv_dist3(tv, p) - tv->zbuff->zmin) / (tv->zbuff->zmax - tv->zbuff->zmin);
    if (z < 0.0 || z > 1.0)
        *iz = INT_MAX;
    else
        *iz = ROUND(INT_MAX * z);
    return (tv_proj3(tv, p));
}

void    tv_camera3(Tv * tv, Vec3 centre, double radius, double pscale, Vec3 aim, Vec3 down)
{
    tv->scalex = tv->scaley = (float)(0.5 * MIN(tv->width, tv->height) / radius);
    tv->cx = (float)(tv->width / 2.0);
    tv->cy = (float)(tv->height / 2.0);
    tv->centre3 = centre;
    tv->radius3 = (float)radius;
    vec3_basis(aim, down, &tv->ex3, &tv->ey3, &tv->ez3);
    tv->pscale = (float)pscale;
    tv->pcentre = vec3_diff(centre, vec3_times(pscale * radius, tv->ez3));
}

void    tv_aim3(Tv * tv, Vec3 aim)
{
    vec3_basis(aim, vec3_ey(), &tv->ex3, &tv->ey3, &tv->ez3);
}

void    tv_orient3(Tv * tv, Vec3 aim, Vec3 down)
{
    vec3_basis(aim, down, &tv->ex3, &tv->ey3, &tv->ez3);
    tv->pcentre = vec3_diff(tv->centre3,
                       vec3_times(tv->pscale * tv->radius3, tv->ez3));
}

void    tv_set_axis(Tv * tv, Vec3 axis)
{
    tv->axis = vec3_unit(axis);
    tv->axis_set = true;
}

void    tv_unset_axis(Tv * tv)
{
    tv->axis_set = false;
}

double  tv_pscale(Tv * tv, double pscale)
{
    double  oldpscale = tv->pscale;

    tv->pscale = (float)pscale;
    return (oldpscale);
}

void    tv_set_proj3(Tv * tv, int type)
{
    if (tv == NULL)
        return;

    switch (type)
    {
    case ORTH:
        tv_set_proj3_orth(tv);
        break;
    case PERSP:
        tv_set_proj3_persp(tv);
        break;
    }
}

---