C:/cmcintos/defOrgs/examples/vesselCrawler/source/Sense_HessianBased.h

#if defined (_MSC_VER) && (_MSC_VER >= 1000)
#pragma once
#endif
#ifndef _Sense_HessianBased
#define _Sense_HessianBased

#include "itkGradientMagnitudeRecursiveGaussianImageFilter.h"
#include "itkGradientRecursiveGaussianImageFilter.h"
#include "itkImageFileWriter.h"
#include "Sensor.h"
namespace mial
{


        template<class DataType, class TInputImage, class TGradientImage, int nDims>
        class Sense_HessianBased: public Sensor
        {
                //****************************************************
                //Member variables
                //****************************************************
        public:

                typedef TInputImage                                                                     InputImageType;
                typedef typename InputImageType::ConstPointer           InputImagePointer;

                typedef itk::Image<DataType,nDims>                                      OutputImageType;
                typedef typename OutputImageType::ConstPointer          OutputImagePointer;


                typedef TGradientImage                                                          GradientImageType;
                typedef typename GradientImageType::Pointer                     GradientImagePointer;

                struct In: public Sensor::sensorIn{
                        //Sensor::sensorIn baseIn;//Import standards
                        //Define all custom here
                        InputImagePointer imageIn;
                        DataType sigma;
                };
                typedef In sensorIn; 
                struct Out: public Sensor::sensorOut{
                        //Sensor::sensorOut baseOut;//Import standards
                        //Define all custom here
                        GradientImagePointer imageOut;

                        //Testing
                        //InputImagePointer imageOut;
                };
                typedef Out sensorOut;
        private:

                //Define all variables specific to this sensor here
                typedef typename itk::GradientMagnitudeRecursiveGaussianImageFilter<InputImageType, OutputImageType> MagnitudeFilterType;
                typedef typename itk::GradientRecursiveGaussianImageFilter<OutputImageType, GradientImageType> GradientFilterType;

                //Testing
                //typedef typename itk::GradientMagnitudeRecursiveGaussianImageFilter<InputImageType, InputImageType> MagnitudeFilterType;

                typename MagnitudeFilterType::Pointer magnitudeFilter;
                typename GradientFilterType::Pointer gradientFilter;

                sensorOut sensorOutput;


                //****************************************************
                //Member functions
                //****************************************************
        public:

                Sense_HessianBased();//default constructor

                //Run the sensor
                virtual void run(void * i);


                //Overide in derived class if custom sensor output desired
                virtual void * getOuput(){ return (void *)&sensorOutput;};

        private:
                //write an image to file
                void writeImageToFile(OutputImagePointer img,char * name);

        };
} // end namespace mial
#include "Sense_HessianBased.cxx"
#endif

---