Linux Cross Reference
TINA5/tina-libs/tina/image/imgSeq_imstack.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-libs/tina/image/imgSeq_imstack.c,v $
   * Date    :  $Date: 2003/09/22 16:09:02 $
   * Version :  $Revision: 1.5 $
   * CVS Id  :  $Id: imgSeq_imstack.c,v 1.5 2003/09/22 16:09:02 tony Exp $
   *
   * Author  :  Legacy TINA
   *
   * Notes : image stack handling
   *
   */
  
  /** 
   *  @file
   *  @brief Image stack handling
   *
   */
  
  
  #include "imgSeq_imstack.h"
  
  #if HAVE_CONFIG_H
  #include <config.h>
  #endif
  
  #include <stdio.h>
  #include <tina/sys/sysDef.h>
  #include <tina/sys/sysPro.h>
  #include <tina/image/img_GenDef.h>
  #include <tina/image/img_GenPro.h>
  #include <tina/image/img_SeqDef.h>
  #include <tina/image/imgSeq_voi.h>
  
  /*
  Allocate a slice with space for NVOI volumes of interest.
  */
  static Seqslice *seqslice_alloc(double z)
  {
          Seqslice       *slice = ralloc(sizeof(Seqslice));
          int             i;
          slice->im = NULL;
          for (i = 0; i < NVOI; i++)
                  slice->voi[i] = svoi_alloc();
          slice->z = z;
          return (slice);
  }
  
  /*
  Free a slice.
  */
  static void     seqslice_free(Seqslice * slice)
  {
          int             i;
          if (slice == NULL)
                  return;
          im_free(slice->im);
          for (i = 0; i < NVOI; i++)
                  svoi_free(slice->voi[i]);
          rfree(slice);
  }
  
  /*
  Find power of 2 greater than x
  */
  static int      above2n(int x)
 {
         int             y = 1;
         while (y < x)
                 y *= 2;
         return (y);
 }
 
 /*
 Return an image stack of given image variable type (uchar_v etc.)
 with slice numbers ranging from lz to uz (excluding uz).
 */
 Br_imstack     *seq_imstack_make(Vartype vtype, int lz, int uz, double zscale)
 {
         Br_imstack     *imstack = ralloc(sizeof(Br_imstack));
         int             z;
 
         imstack->vtype = vtype;
         imstack->lx = imstack->ux = imstack->x = 0;
         imstack->ly = imstack->uy = imstack->y = 0;
         imstack->lz = lz;
         imstack->uz = uz;
         imstack->z = (lz + uz) / 2;
         imstack->zmark = 0;
         imstack->slice = tvector_alloc(lz, uz, Seqslice *);
         imstack->zscale = zscale;
         for (z = lz; z < uz; z++)
                 imstack->slice[z] = seqslice_alloc(seq_imstack_zscaled(imstack, z));
         /*
          * imstack->stats = NULL; imstack->tube = NULL;
          */
         return (imstack);
 }
 
 /*
 Free an image stack and associated memory
 */
 void            seq_imstack_free(Br_imstack * imstack)
 {
         int             z;
 
         if (imstack == NULL)
                 return;
         for (z = imstack->lz; z < imstack->uz; z++)
                 seqslice_free(imstack->slice[z]);
         tvector_free(imstack->slice, imstack->lz, Seqslice *);
         /*
          * rfree(imstack->stats); tube_free(imstack->tube);
          */
         rfree(imstack);
 }
 
 /*
 Get true z-coordinate associated with slice number z
 */
 double          seq_imstack_zscaled(Br_imstack * imstack, int z)
 {
         z -= (imstack->lz + imstack->uz) / 2;
         return (imstack->zscale * z);
 }
 
 /*
 Get slice number associated wth true z-coordinate z
 */
 int             seq_imstack_zunscaled(Br_imstack * imstack, double z)
 {
         int             zi;
         z /= imstack->zscale;
         z += (imstack->lz + imstack->uz) / 2;
         zi = ROUND(z);
         return (zi);
 }
 
 /*
 Return transverse image in (y, z) plane at current x-cursor position.
 */
 Imrect         *seq_imstack_xslice(Br_imstack * imstack, int x)
 {
         Imrect         *im;
         Imregion       *roi;
         float          *row;
         int             z, width, height;
 
         if (imstack == NULL)
                 return (NULL);
 
         height = above2n(imstack->uz);
         width = above2n(imstack->uy);
         roi = roi_alloc(imstack->ly, imstack->lz, imstack->uy, imstack->uz);
         im = im_alloc(height, width, roi, imstack->vtype);
         rfree(roi);
 
         row = tvector_alloc(imstack->ly, imstack->uy, float);
         for (z = imstack->lz; z < imstack->uz; z++)
         {
                 im_get_colf(row, imstack->slice[z]->im, x, imstack->ly, imstack->uy);
                 im_put_rowf(row, im, z, imstack->ly, imstack->uy);
         }
         tvector_free(row, imstack->ly, float);
 
         return (im);
 }
 
 /*
 Return transverse image in (x, z) plane at current y-cursor position.
 */
 Imrect         *seq_imstack_yslice(Br_imstack * imstack, int y)
 {
         Imrect         *im;
         Imregion       *roi;
         float          *row;
         int             z, width, height;
 
         if (imstack == NULL)
                 return (NULL);
 
         height = above2n(imstack->uz);
         width = above2n(imstack->ux);
         roi = roi_alloc(imstack->lx, imstack->lz, imstack->ux, imstack->uz);
         im = im_alloc(height, width, roi, imstack->vtype);
         rfree(roi);
 
         row = tvector_alloc(imstack->lx, imstack->ux, float);
         for (z = imstack->lz; z < imstack->uz; z++)
         {
                 im_get_rowf(row, imstack->slice[z]->im, y, imstack->lx, imstack->ux);
                 im_put_rowf(row, im, z, imstack->lx, imstack->ux);
         }
         tvector_free(row, imstack->lx, float);
 
         return (im);
 }