Main Page   Groups   Namespace List   Class Hierarchy   Alphabetical List   Compound List   File List   Namespace Members   Compound Members   File Members   Concepts

itkPhasedArray3DSpecialCoordinatesImage.h

Go to the documentation of this file.
00001 /*=========================================================================
00002 
00003   Program:   Insight Segmentation & Registration Toolkit
00004   Module:    $RCSfile: itkPhasedArray3DSpecialCoordinatesImage.h,v $
00005   Language:  C++
00006   Date:      $Date: 2009-02-19 19:41:23 $
00007   Version:   $Revision: 1.22 $
00008 
00009   Copyright (c) Insight Software Consortium. All rights reserved.
00010   See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
00011 
00012      This software is distributed WITHOUT ANY WARRANTY; without even
00013      the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
00014      PURPOSE.  See the above copyright notices for more information.
00015 
00016 =========================================================================*/
00017 #ifndef __itkPhasedArray3DSpecialCoordinatesImage_h
00018 #define __itkPhasedArray3DSpecialCoordinatesImage_h
00019 
00020 #include "itkSpecialCoordinatesImage.h"
00021 #include "itkImageRegion.h"
00022 #include "itkPoint.h"
00023 #include "itkContinuousIndex.h"
00024 #include "vnl/vnl_math.h"
00025 
00026 
00027 namespace itk
00028 {
00029 
00090 template <class TPixel>
00091 class ITK_EXPORT PhasedArray3DSpecialCoordinatesImage :
00092 public SpecialCoordinatesImage<TPixel,3>
00093 {
00094 public:
00096   typedef PhasedArray3DSpecialCoordinatesImage Self;
00097   typedef SpecialCoordinatesImage<TPixel,3>    Superclass;
00098   typedef SmartPointer<Self>                   Pointer;
00099   typedef SmartPointer<const Self>             ConstPointer;
00100   typedef WeakPointer<const Self>              ConstWeakPointer;
00101 
00103   itkNewMacro(Self);
00104 
00106   itkTypeMacro(PhasedArray3DSpecialCoordinatesImage, SpecialCoordinatesImage);
00107 
00110   typedef TPixel PixelType;
00111 
00113   typedef TPixel ValueType;
00114 
00119   typedef TPixel InternalPixelType;
00120 
00121   typedef typename Superclass::IOPixelType   IOPixelType;
00122   
00125   typedef DefaultPixelAccessor< PixelType > AccessorType;
00126 
00130   typedef DefaultPixelAccessorFunctor< Self > AccessorFunctorType;
00131 
00136   itkStaticConstMacro(ImageDimension, unsigned int, 3);
00137 
00139   typedef ImportImageContainer<unsigned long, PixelType> PixelContainer;
00140 
00142   typedef typename Superclass::IndexType      IndexType;
00143   typedef typename Superclass::IndexValueType IndexValueType;
00144 
00146   typedef typename Superclass::OffsetType OffsetType;
00147 
00149   typedef typename Superclass::SizeType   SizeType;
00150 
00154   typedef typename Superclass::RegionType RegionType;
00155 
00162   typedef typename Superclass::SpacingType SpacingType;
00163 
00168   typedef typename Superclass::PointType PointType;
00169 
00171   typedef typename PixelContainer::Pointer      PixelContainerPointer;
00172   typedef typename PixelContainer::ConstPointer PixelContainerConstPointer;
00173 
00178   template<class TCoordRep>
00179   bool TransformPhysicalPointToContinuousIndex(
00180               const Point<TCoordRep, 3>& point,
00181               ContinuousIndex<TCoordRep, 3>& index   ) const
00182     {
00183     RegionType region = this->GetLargestPossibleRegion();
00184     double maxAzimuth =    region.GetSize(0) - 1;
00185     double maxElevation =  region.GetSize(1) - 1;
00186     
00187     // Convert Cartesian coordinates into angular coordinates
00188     TCoordRep azimuth   = vcl_atan(point[0] / point[2]);
00189     TCoordRep elevation = vcl_atan(point[1] / point[2]);
00190     TCoordRep radius    = vcl_sqrt(point[0] * point[0]
00191                               + point[1] * point[1]
00192                               + point[2] * point[2] );
00193     
00194     // Convert the "proper" angular coordinates into index format
00195     index[0] = static_cast<TCoordRep>( (azimuth/m_AzimuthAngularSeparation)
00196                                        + (maxAzimuth/2.0)   );
00197     index[1] = static_cast<TCoordRep>( (elevation/m_ElevationAngularSeparation)
00198                                        + (maxElevation/2.0) );
00199     index[2] = static_cast<TCoordRep>( ( (radius-m_FirstSampleDistance)
00200                                               / m_RadiusSampleSize) );
00201     
00202     // Now, check to see if the index is within allowed bounds
00203     const bool isInside = region.IsInside( index );
00204 
00205     return isInside;
00206     }
00207 
00212   template<class TCoordRep>
00213   bool TransformPhysicalPointToIndex(
00214             const Point<TCoordRep, 3>& point,
00215             IndexType & index                                ) const
00216     {
00217     RegionType region = this->GetLargestPossibleRegion();
00218     double maxAzimuth =    region.GetSize(0) - 1;
00219     double maxElevation =  region.GetSize(1) - 1;
00221 
00222     // Convert Cartesian coordinates into angular coordinates
00223     TCoordRep azimuth   = vcl_atan(point[0] / point[2]);
00224     TCoordRep elevation = vcl_atan(point[1] / point[2]);
00225     TCoordRep radius    = vcl_sqrt(point[0] * point[0]
00226                                 + point[1] * point[1]
00227                                 + point[2] * point[2] );
00228     
00229     // Convert the "proper" angular coordinates into index format
00230     index[0] = static_cast<IndexValueType>(
00231       (azimuth/m_AzimuthAngularSeparation)
00232       + (maxAzimuth/2.0) );
00233     index[1] = static_cast<IndexValueType>(
00234       (elevation/m_ElevationAngularSeparation)
00235       + (maxElevation/2.0) );
00236     index[2] = static_cast<IndexValueType>(
00237       ((radius-m_FirstSampleDistance)
00238        / m_RadiusSampleSize ) );
00239     
00240     // Now, check to see if the index is within allowed bounds
00241     const bool isInside = region.IsInside( index );
00242 
00243     return isInside;
00244     }
00245 
00250   template<class TCoordRep>
00251   void TransformContinuousIndexToPhysicalPoint(
00252             const ContinuousIndex<TCoordRep, 3>& index,
00253             Point<TCoordRep, 3>& point        ) const
00254     {
00255     RegionType region = this->GetLargestPossibleRegion();
00256     double maxAzimuth =    region.GetSize(0) - 1;
00257     double maxElevation =  region.GetSize(1) - 1;
00259 
00260     // Convert the index into proper angular coordinates
00261     TCoordRep azimuth   = ( index[0] - (maxAzimuth/2.0) )
00262                           * m_AzimuthAngularSeparation;
00263     TCoordRep elevation = ( index[1] - (maxElevation/2.0) )
00264                           * m_ElevationAngularSeparation;
00265     TCoordRep radius    = (index[2]*m_RadiusSampleSize)+m_FirstSampleDistance;
00266     
00267     // Convert the angular coordinates into Cartesian coordinates
00268     TCoordRep tanOfAzimuth    = vcl_tan(azimuth);
00269     TCoordRep tanOfElevation  = vcl_tan(elevation);
00270     point[2] = static_cast<TCoordRep>( radius /
00271            vcl_sqrt(1 +
00272                     tanOfAzimuth*tanOfAzimuth +
00273                     tanOfElevation*tanOfElevation));
00274     point[1] = static_cast<TCoordRep>( point[2] * tanOfElevation );
00275     point[0] = static_cast<TCoordRep>( point[2] * tanOfAzimuth );
00276     }
00277 
00283   template<class TCoordRep>
00284   void TransformIndexToPhysicalPoint(
00285                       const IndexType & index,
00286                       Point<TCoordRep, 3>& point ) const
00287     {
00288     RegionType region = this->GetLargestPossibleRegion();
00289     double maxAzimuth =    region.GetSize(0) - 1;
00290     double maxElevation =  region.GetSize(1) - 1;
00292 
00293     // Convert the index into proper angular coordinates
00294     TCoordRep azimuth = 
00295       (static_cast<double>(index[0]) - (maxAzimuth/2.0) )
00296       * m_AzimuthAngularSeparation;
00297     TCoordRep elevation =
00298       (static_cast<double>(index[1]) - (maxElevation/2.0) )
00299       * m_ElevationAngularSeparation;
00300     TCoordRep radius =
00301       (static_cast<double>(index[2]) * m_RadiusSampleSize)
00302       + m_FirstSampleDistance;
00303     
00304     // Convert the angular coordinates into Cartesian coordinates
00305     TCoordRep tanOfAzimuth    = vcl_tan(azimuth);
00306     TCoordRep tanOfElevation  = vcl_tan(elevation);
00307     point[2] = static_cast<TCoordRep>(
00308       radius / vcl_sqrt(
00309         1.0 + tanOfAzimuth*tanOfAzimuth + tanOfElevation*tanOfElevation) );
00310     point[1] = static_cast<TCoordRep>( point[2] * tanOfElevation );
00311     point[0] = static_cast<TCoordRep>( point[2] * tanOfAzimuth );
00312     }
00313   
00314   
00316   itkSetMacro(AzimuthAngularSeparation, double);
00317 
00319   itkSetMacro(ElevationAngularSeparation, double);
00320 
00322   itkSetMacro(RadiusSampleSize, double);
00323 
00325   itkSetMacro(FirstSampleDistance, double);
00326 
00327 #ifdef ITK_USE_ORIENTED_IMAGE_DIRECTION
00328   template<class TCoordRep>
00329   void TransformLocalVectorToPhysicalVector(
00330     const FixedArray<TCoordRep, 3> & inputGradient,
00331           FixedArray<TCoordRep, 3> & outputGradient ) const
00332     {
00333     }
00334 #endif
00335 protected:
00336   PhasedArray3DSpecialCoordinatesImage()
00337     {
00338     m_RadiusSampleSize = 1;
00339     m_AzimuthAngularSeparation   =  1 * (2.0*vnl_math::pi/360.0); // 1 degree
00340     m_ElevationAngularSeparation =  1 * (2.0*vnl_math::pi/360.0); // 1 degree
00341     m_FirstSampleDistance = 0;
00342     }
00343   virtual ~PhasedArray3DSpecialCoordinatesImage() {};
00344   void PrintSelf(std::ostream& os, Indent indent) const;
00345   
00346 private:
00347   PhasedArray3DSpecialCoordinatesImage(const Self&);//purposely not implemented
00348   void operator=(const Self&); //purposely not implemented
00349   
00350   double  m_AzimuthAngularSeparation;   // in radians
00351   double  m_ElevationAngularSeparation; // in radians
00352   double  m_RadiusSampleSize;
00353   double  m_FirstSampleDistance;
00354   
00355 };
00356 } // end namespace itk
00357 
00358 // Define instantiation macro for this template.
00359 #define ITK_TEMPLATE_PhasedArray3DSpecialCoordinatesImage(_, EXPORT, x, y) namespace itk { \
00360   _(1(class EXPORT PhasedArray3DSpecialCoordinatesImage< ITK_TEMPLATE_1 x >)) \
00361   namespace Templates { typedef PhasedArray3DSpecialCoordinatesImage< ITK_TEMPLATE_1 x > \
00362                                          PhasedArray3DSpecialCoordinatesImage##y; } \
00363   }
00364 
00365 #if ITK_TEMPLATE_EXPLICIT
00366 # include "Templates/itkPhasedArray3DSpecialCoordinatesImage+-.h"
00367 #endif
00368 
00369 #if ITK_TEMPLATE_TXX
00370 # include "itkPhasedArray3DSpecialCoordinatesImage.txx"
00371 #endif
00372 
00373 #endif
00374 

Generated at Tue Sep 15 04:22:01 2009 for ITK by doxygen 1.5.8 written by Dimitri van Heesch, © 1997-2000