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

#ifndef _Sense_HessianBased_CXX
#define _Sense_HessianBased_CXX

#include "Sense_HessianBased.h"


namespace mial
{
        //Default constructor
        template<class DataType, class TInputImage, class TGradientImage, int nDims>
                Sense_HessianBased<DataType,TInputImage,TGradientImage,nDims>
                ::Sense_HessianBased()
                : Sensor()
        {
                //Instantiate the filters
                magnitudeFilter = MagnitudeFilterType::New();
                gradientFilter = GradientFilterType::New();
        }

        template<class DataType, class TInputImage, class TGradientImage, int nDims>
                void 
                Sense_HessianBased<DataType,TInputImage,TGradientImage,nDims>
		::writeImageToFile(OutputImagePointer img,char * name)
        {
                typedef typename itk::ImageFileWriter<OutputImageType> WriterType;
                WriterType::Pointer writer = WriterType::New(); 
                writer->SetFileName(name);
                writer->SetInput(img);
                writer->Update();
        }

        template<class DataType, class TInputImage, class TGradientImage, int nDims>
                void 
                Sense_HessianBased<DataType, TInputImage, TGradientImage, nDims>
                ::run(void * i)
        {
                sensorIn in = *(reinterpret_cast<sensorIn *>(i)); //Typecast the input to its desired form

                //Calculate the gradient magnitude at scale sigma
                magnitudeFilter->SetSigma(in.sigma);
                magnitudeFilter->SetNormalizeAcrossScale(1); //Set so magnitude is normalized across scales.
                magnitudeFilter->SetInput(in.imageIn);

                //testing
                //this->writeImageToFile(magnitudeFilter->GetOutput(),"mag.mhd");
                //this->sensorOutput.imageOut = magnitudeFilter->GetOutput();
                //magnitudeFilter->Update();

                //Calculate the gradients
                gradientFilter->SetSigma(in.sigma);
                gradientFilter->SetNormalizeAcrossScale(1); //Set so magnitude is normalized across scales.
                gradientFilter->SetInput(magnitudeFilter->GetOutput());

                this->sensorOutput.imageOut = gradientFilter->GetOutput();

                gradientFilter->Update();
        }
}
#endif

---