vtkPlusRfToBrightnessConvert.h

Go to the documentation of this file.

/*=Plus=header=begin======================================================
  Program: Plus
  Copyright (c) Laboratory for Percutaneous Surgery. All rights reserved.
  See License.txt for details.
=========================================================Plus=header=end*/

#ifndef __vtkPlusRfToBrightnessConvert_h
#define __vtkPlusRfToBrightnessConvert_h

#include "vtkPlusImageProcessingExport.h"
#include "vtkThreadedImageAlgorithm.h"

class vtkPlusImageProcessingExport vtkPlusRfToBrightnessConvert : public vtkThreadedImageAlgorithm
{
public:
  static vtkPlusRfToBrightnessConvert *New();
  vtkTypeMacro(vtkPlusRfToBrightnessConvert,vtkThreadedImageAlgorithm);
  virtual void PrintSelf(ostream& os, vtkIndent indent) VTK_OVERRIDE;
  
  virtual PlusStatus ReadConfiguration(vtkXMLDataElement* rfToBrightnessElement); 

  virtual PlusStatus WriteConfiguration(vtkXMLDataElement* rfToBrightnessElement); 

  vtkSetMacro(ImageType, US_IMAGE_TYPE);
  vtkGetMacro(ImageType, US_IMAGE_TYPE);

  vtkSetMacro(NumberOfHilbertFilterCoeffs, int);
  vtkGetMacro(NumberOfHilbertFilterCoeffs, int);

  vtkSetMacro(BrightnessScale, double);
  vtkGetMacro(BrightnessScale, double);

protected:
  vtkPlusRfToBrightnessConvert();
  ~vtkPlusRfToBrightnessConvert();
  
  virtual int RequestInformation(vtkInformation*,
                                 vtkInformationVector**,
                                 vtkInformationVector* outputVector);

  void ThreadedRequestData( vtkInformation *request,
                            vtkInformationVector **inputVector,
                            vtkInformationVector *outputVector,
                            vtkImageData ***inData, vtkImageData **outData,
                            int outExt[6], int id);

  virtual void ComputeHilbertTransformCoeffs();

  template<typename ScalarType>
  void ThreadedLineByLineHilbertTransform(int inExt[6], int outExt[6], vtkImageData ***inData, vtkImageData **outData, int threadId);

  template<typename ScalarType>
  PlusStatus ComputeHilbertTransform(ScalarType *hilbertTransformOutput, ScalarType *input, int npt);
  
  template<typename ScalarType>
  void ComputeAmplitudeILineQLine(unsigned char *ampl, ScalarType *inputSignal, ScalarType *inputSignalHilbertTransformed, int npt);
  
  template<typename ScalarType>
  void ComputeAmplitudeIqLine(unsigned char *ampl, ScalarType *inputSignal, const int npt);

  double BrightnessScale;

  int NumberOfHilbertFilterCoeffs;

  std::vector<double> HilbertTransformCoeffs;

  US_IMAGE_TYPE ImageType;

private:
  vtkPlusRfToBrightnessConvert(const vtkPlusRfToBrightnessConvert&);  // Not implemented.
  void operator=(const vtkPlusRfToBrightnessConvert&);  // Not implemented.
};

#endif