vtkPlusUsScanConvertCurvilinear.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 __vtkPlusUsScanConvertCurvilinear_h
#define __vtkPlusUsScanConvertCurvilinear_h

#include "vtkPlusImageProcessingExport.h"
#include "vtkPlusUsScanConvert.h"

class vtkPlusImageProcessingExport vtkPlusUsScanConvertCurvilinear : public vtkPlusUsScanConvert
{
public:
  static vtkPlusUsScanConvertCurvilinear* New();
  vtkTypeMacro(vtkPlusUsScanConvertCurvilinear,vtkPlusUsScanConvert);
  virtual void PrintSelf(ostream& os, vtkIndent indent) VTK_OVERRIDE;

  virtual const char* GetTransducerGeometry()
  {
    return "CURVILINEAR";
  }

  virtual PlusStatus ReadConfiguration(vtkXMLDataElement* scanConversionElement);

  virtual PlusStatus WriteConfiguration(vtkXMLDataElement* scanConversionElement);

  virtual vtkImageData* GetOutput();

  struct InterpolatedPoint
  {
    double weightCoefficients[4];
    int inputPixelIndex;
    // Position of the output pixel (in the image matrix) */
    int outputPixelIndex;
  };

  const std::vector<InterpolatedPoint>& GetInterpolatedPointArray()
  {
    return this->InterpolatedPointArray;
  };

  vtkSetMacro(RadiusStartMm, double);
  vtkGetMacro(RadiusStartMm, double);
  vtkSetMacro(RadiusStopMm, double);
  vtkGetMacro(RadiusStopMm, double);
  vtkSetMacro(ThetaStartDeg, double);
  vtkSetMacro(ThetaStopDeg, double);
  vtkSetMacro(OutputImageStartDepthMm, double);

  PlusStatus GetScanLineEndPoints(int scanLineIndex, double scanlineStartPoint_OutputImage[4],double scanlineEndPoint_OutputImage[4]);

  virtual double GetDistanceBetweenScanlineSamplePointsMm();

protected:
  vtkPlusUsScanConvertCurvilinear();
  virtual ~vtkPlusUsScanConvertCurvilinear();

  virtual int SplitExtent(int splitExt[6], int startExt[6], int num, int total);

  virtual int RequestInformation(vtkInformation*, vtkInformationVector**, vtkInformationVector*);

  virtual int RequestUpdateExtent (vtkInformation*, vtkInformationVector**, vtkInformationVector*);
  
  virtual void AllocateOutputData(vtkImageData *output, vtkInformation* outInfo, int *uExtent);

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

  double OutputImageStartDepthMm;

  double RadiusStartMm;
  double RadiusStopMm;

  double ThetaStartDeg;
  double ThetaStopDeg;

  double OutputIntensityScaling;

  std::vector<InterpolatedPoint> InterpolatedPointArray;

  int InterpInputImageExtent[6];
  double InterpRadiusStartMm;
  double InterpRadiusStopMm;
  double InterpThetaStartDeg;
  double InterpThetaStopDeg;
  int InterpOutputImageExtent[6];
  double InterpOutputImageSpacing[3];
  double InterpTransducerCenterPixel[2];
  double InterpIntensityScaling;

  void ComputeInterpolatedPointArray(
    int* inputImageExtent, double radiusStartMm, double radiusStopMm, double thetaStartDeg, double thetaStopDeg,
    int* outputImageExtent, double* outputImageSpacing, double* transducerCenterPixel, double intensityScaling
  );

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

#endif