mitkCylindricToCartesianFilter.h

Go to the documentation of this file.

/*=========================================================================

Program:   Medical Imaging & Interaction Toolkit
Language:  C++
Date:      $Date$
Version:   $Revision$

Copyright (c) German Cancer Research Center, Division of Medical and
Biological Informatics. All rights reserved.
See MITKCopyright.txt or https://www.mitk.org/copyright.html for details.

This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE.  See the above copyright notices for more information.

=========================================================================*/


#ifndef MITKCYLINDRICTOCARTESIANFILTER_H_HEADER_INCLUDED_C1F48A22
#define MITKCYLINDRICTOCARTESIANFILTER_H_HEADER_INCLUDED_C1F48A22

#include "mitkCommon.h"
#include "MitkExtExports.h"
#include "mitkImageToImageFilter.h"

#ifndef M_PI
  #define M_PI       3.14159265358979323846
#endif

//#endif
namespace mitk {
//##Documentation
//## @brief Filter to convert a cylindric image into a cartesian image
//##
//## @todo position of input image is not yet used to calculate position of the output
//## @todo tag with transducer position not yet transformed.
//## @todo the line limiting the sector is assumed to be available as a pic tag. Is this ok?
//## @ingroup Process
class MitkExt_EXPORT CylindricToCartesianFilter : public ImageToImageFilter
{
public:
  mitkClassMacro(CylindricToCartesianFilter, ImageToImageFilter);

  itkNewMacro(Self);

  //##Documentation
  //## @brief Set background grey level
  itkSetMacro(OutsideValue, float);
  //##Documentation
  //## @brief Get background grey level
  itkGetMacro(OutsideValue, float);

  //##Documentation
  //## @brief Set the size in x-direction of the converted image.
  //## 
  //## A value of 0 (default) means that the filter uses the x-size of the input image.
  itkSetMacro(TargetXSize, unsigned int);
  //##Documentation
  //## @brief Get the size in x-direction of the converted image.
  itkGetMacro(TargetXSize, unsigned int);

  //##Documentation
  //## @brief Build some tables to speed up transformation, e.g., to avoid repeated calulations of the same sinus/cosinus.
  //## @todo Free tables if they were already previously allocated
  void buildTransformShortCuts(int orig_xsize, int orig_ysize, int orig_zsize, int new_xsize, mitkIpPicDescriptor * &rt_pic, mitkIpPicDescriptor * &phit_pic, mitkIpPicDescriptor * &fr_pic, mitkIpPicDescriptor * &fphi_pic, unsigned int * &zt, float * &fz);

  //##Documentation
  //## @brief Build some tables to speed up transformation of images that only contain data within a sector
  //## @todo Free tables if they were already previously allocated
  void buildConeCutOffShortCut(int orig_xsize, int orig_ysize, mitkIpPicDescriptor *rt_pic, mitkIpPicDescriptor *fr_pic, float a, float b, mitkIpPicDescriptor * &coneCutOff_pic);
protected:
  mitkIpPicDescriptor * rt_pic, * phit_pic, *fr_pic, * fphi_pic, * coneCutOff_pic;
  float * fz;
  unsigned int * zt;
  //##Description 
  //## @brief Cut off line for images that only contain data within a sector
  //## 
  //## y=a x+b
  float a;
  //##Description 
  //## @brief Cut off line for images that only contain data within a sector
  //## 
  //## y=a x+b
  float b;

  //##Description 
  //## @brief Time when Header was last initialized
  itk::TimeStamp m_TimeOfHeaderInitialization;
protected:
  virtual void GenerateData();

  virtual void GenerateOutputInformation();

  virtual void GenerateInputRequestedRegion();

  CylindricToCartesianFilter();

  ~CylindricToCartesianFilter();

  float m_OutsideValue;
  unsigned int m_TargetXSize;
};

} // namespace mitk

#endif /* MITKCYLINDRICTOCARTESIANFILTER_H_HEADER_INCLUDED_C1F48A22 */