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

#ifndef _Sense_ProjectiveSpherical_CXX
#define _Sense_ProjectiveSpherical_CXX

#include "Sense_ProjectiveSpherical.h"
#define PI 3.141592654

namespace mial
{
        //Default constructor
        template<class DataType, class TInputImage, class TGradientImage, int nDims>
                Sense_ProjectiveSpherical<DataType,TInputImage,TGradientImage,nDims>
                ::Sense_ProjectiveSpherical()
                : Sensor()
        {
                //Instantiate the filters
        }

        template<class DataType, class TInputImage, class TGradientImage, int nDims>
                void 
                Sense_ProjectiveSpherical<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_ProjectiveSpherical<DataType, TInputImage, TGradientImage, nDims>
                ::run(void * i)
        {
                //TODO change to base class of sensorIn
                sensorIn in = *(reinterpret_cast<sensorIn *>(i)); //Typecast the input to its desired form

                //Build the hemi-sphere (a rotated, scaled, and translated unit sphere)
                int N = ceil(sqrt(1000));// //The number of points on its surface
                DataType theta =0;
                DataType phi =0;
                MatrixType points(1000,3);
                VectorType pnt(3);
                MatrixType M(3,3);
                VectorType trans(3);
                //TODO fill in r*M,trans
                int cnt =0;
                for(int i=0; i < N; i++)
                {
                        theta = -PI+2*i/(N-1)*PI;
                        for(int j=0; j < N; j++)
                        {
                                phi = (j/(N-1))*(PI/2);
                                pnt(1) = cos(phi)*cos(theta);
                                pnt(2) = cos(phi)*sin(theta);
                                pnt(3) = sin(phi);
                                points.set_row(cnt) = M*pnt + trans;
                                cnt++;
                        }
                }
                //Obtain the spherical slice

                //Project onto the plane

                //Calculate centroids
                
                //Calculate the gradients
                
                //Set the sensorOutput
                this->sensorOutput.imageOut = gradientFilter->GetOutput();

                gradientFilter->Update();
        }
}
#endif

---