Detailed Description

A VTK filter that creates a point cloud representing a surface defined by pixel values

Definition at line 16 of file VtkImageDataToSurfacePointsFilter.h.

#include <VtkImageDataToSurfacePointsFilter.h>

Inheritance diagram for VtkImageDataToSurfacePointsFilter:

Collaboration diagram for VtkImageDataToSurfacePointsFilter:

Public Member Functions
	vtkTypeMacro (VtkImageDataToSurfacePointsFilter, vtkPolyDataAlgorithm)
void	PrintSelf (ostream &os, vtkIndent indent) override
	Prints information about itself.
	vtkGetMacro (PointsPerPixel, vtkIdType)
	vtkSetMacro (PointsPerPixel, vtkIdType)

Static Public Member Functions
static VtkImageDataToSurfacePointsFilter *	New ()
	Create a new objects (required because of VTKs reference counting)

Protected Member Functions
	VtkImageDataToSurfacePointsFilter ()
	~VtkImageDataToSurfacePointsFilter () override=default
int	FillInputPortInformation (int port, vtkInformation *info) override
	Sets input port to vtkImageData.
int	RequestData (vtkInformation request, vtkInformationVector inputVector, vtkInformationVector outputVector) override
	Updates the graphical object.

Private Member Functions
void	createPointSurface (vtkSmartPointer< vtkPoints > &points, vtkSmartPointer< vtkCellArray > &cells, std::size_t pnt_idx, MathLib::Point3d const &min_pnt, MathLib::Point3d const &max_pnt, GeoLib::Raster const &raster)
	Creates the point objects based on the parameters set by the user.
double	getRandomNumber (double const &min, double const &max) const
	Returns a random number in [min, max].

Private Attributes
vtkIdType	PointsPerPixel {20}

Constructor & Destructor Documentation

◆ VtkImageDataToSurfacePointsFilter()

VtkImageDataToSurfacePointsFilter::VtkImageDataToSurfacePointsFilter ( )

protecteddefault

Referenced by New(), and vtkTypeMacro().

◆ ~VtkImageDataToSurfacePointsFilter()

VtkImageDataToSurfacePointsFilter::~VtkImageDataToSurfacePointsFilter ( )

overrideprotecteddefault

Member Function Documentation

◆ createPointSurface()

void VtkImageDataToSurfacePointsFilter::createPointSurface	(	vtkSmartPointer< vtkPoints > &	points,
		vtkSmartPointer< vtkCellArray > &	cells,
		std::size_t	pnt_idx,
		MathLib::Point3d const &	min_pnt,
		MathLib::Point3d const &	max_pnt,
		GeoLib::Raster const &	raster )

private

Creates the point objects based on the parameters set by the user.

Definition at line 139 of file VtkImageDataToSurfacePointsFilter.cpp.

{
    auto const n_points(static_cast<std::size_t>(this->GetPointsPerPixel()));
    for (std::size_t i = 0; i < n_points; ++i)
    {
        double p[3];
        p[0] = getRandomNumber(min_pnt[0], max_pnt[0]);
        p[1] = getRandomNumber(min_pnt[1], max_pnt[1]);
        p[2] = raster.interpolateValueAtPoint(
            MathLib::Point3d(std::array<double, 3>{{p[0], p[1], 0}}));
        points->SetPoint(pnt_idx + i, p);
        cells->InsertNextCell(1);
        cells->InsertCellPoint(pnt_idx + i);
    }
}

References getRandomNumber(), and GeoLib::Raster::interpolateValueAtPoint().

Referenced by RequestData().

◆ FillInputPortInformation()

int VtkImageDataToSurfacePointsFilter::FillInputPortInformation	(	int	port,
		vtkInformation *	info )

overrideprotected

Sets input port to vtkImageData.

Definition at line 30 of file VtkImageDataToSurfacePointsFilter.cpp.

{
    info->Set(vtkAlgorithm::INPUT_REQUIRED_DATA_TYPE(), "vtkImageData");
    return 1;
}

◆ getRandomNumber()

double VtkImageDataToSurfacePointsFilter::getRandomNumber	(	double const &	min,
		double const &	max ) const

private

Returns a random number in [min, max].

Definition at line 161 of file VtkImageDataToSurfacePointsFilter.cpp.

{
    return (static_cast<double>(std::rand()) / RAND_MAX) * (max - min) + min;
}

Referenced by createPointSurface().

◆ New()

VtkImageDataToSurfacePointsFilter * VtkImageDataToSurfacePointsFilter::New ( )

static

Create a new objects (required because of VTKs reference counting)

References VtkImageDataToSurfacePointsFilter().

Referenced by VtkCompositeImageToSurfacePointsFilter::init().

◆ PrintSelf()

void VtkImageDataToSurfacePointsFilter::PrintSelf	(	ostream &	os,
		vtkIndent	indent )

override

Prints information about itself.

Definition at line 25 of file VtkImageDataToSurfacePointsFilter.cpp.

{
    this->Superclass::PrintSelf(os, indent);
}

◆ RequestData()

int VtkImageDataToSurfacePointsFilter::RequestData	(	vtkInformation *	request,
		vtkInformationVector **	inputVector,
		vtkInformationVector *	outputVector )

overrideprotected

Updates the graphical object.

Definition at line 37 of file VtkImageDataToSurfacePointsFilter.cpp.

{
    vtkDebugMacro(<< "Executing VtkImageDataToSurfacePointsFilter");
 
    vtkInformation* input_info = inputVector[0]->GetInformationObject(0);
    vtkInformation* output_info = outputVector->GetInformationObject(0);
    vtkImageData* input = vtkImageData::SafeDownCast(
        input_info->Get(vtkDataObject::DATA_OBJECT()));
    vtkPolyData* output = vtkPolyData::SafeDownCast(
        output_info->Get(vtkDataObject::DATA_OBJECT()));
 
    void* pixvals = input->GetScalarPointer();
    int n_comp = input->GetNumberOfScalarComponents();
    if (n_comp < 1)
    {
        vtkDebugMacro("Error reading number of components.");
    }
    if (n_comp > 2)
    {
        vtkDebugMacro(
            "RGB colours detected. Using only first channel for computation.");
    }
 
    auto const n_points = static_cast<std::size_t>(input->GetNumberOfPoints());
    if (n_points == 0)
    {
        vtkDebugMacro("No data found!");
        return -1;
    }
 
    double spacing[3];
    input->GetSpacing(spacing);
    int dimensions[3];
    input->GetDimensions(dimensions);
    double origin[3];
    input->GetOrigin(origin);
    MathLib::Point3d const ll(
        std::array<double, 3>{{origin[0], origin[1], origin[2]}});
 
    std::vector<double> pixels;
    pixels.reserve(n_points);
    for (std::size_t i = 0; i < n_points; ++i)
    {
        if ((n_comp == 2 || n_comp == 4) &&
            ((static_cast<float*>(pixvals))[(i + 1) * n_comp - 1] <
             0.00000001f))
        {
            pixels.push_back(-9999);
        }
        else
        {
            pixels.push_back((static_cast<float*>(pixvals))[i * n_comp]);
        }
    }
    GeoLib::Raster const* const raster(new GeoLib::Raster(
        {static_cast<std::size_t>(dimensions[0]),
         static_cast<std::size_t>(dimensions[1]), 1, ll, spacing[0], -9999},
        pixels.begin(),
        pixels.end()));
 
    vtkSmartPointer<vtkPoints> new_points = vtkSmartPointer<vtkPoints>::New();
    new_points->SetNumberOfPoints(PointsPerPixel * n_points);
    vtkSmartPointer<vtkCellArray> cells = vtkSmartPointer<vtkCellArray>::New();
    cells->Allocate(PointsPerPixel * n_points);
 
    double const half_cellsize(spacing[0] / 2.0);
    std::size_t pnt_idx(0);
    for (std::size_t i = 0; i < static_cast<std::size_t>(n_points); ++i)
    {
        // Skip transparent pixels
        if (n_comp == 2 || n_comp == 4)
        {
            if ((static_cast<float*>(pixvals))[(i + 1) * n_comp - 1] <
                0.00000001f)
            {
                continue;
            }
        }
 
        double p[3];
        input->GetPoint(i, p);
        MathLib::Point3d min_pnt{std::array<double, 3>{
            {p[0] - half_cellsize, p[1] - half_cellsize, 0}}};
        MathLib::Point3d max_pnt{std::array<double, 3>{
            {p[0] + half_cellsize, p[1] + half_cellsize, 0}}};
        createPointSurface(
            new_points, cells, pnt_idx, min_pnt, max_pnt, *raster);
        pnt_idx += PointsPerPixel;
    }
 
    output->SetPoints(new_points);
    output->SetVerts(cells);
    output->Squeeze();
 
    vtkDebugMacro(<< "Created " << new_points->GetNumberOfPoints()
                  << " points.");
    return 1;
}