VtkCompositeImageToPointCloudFilter.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "VtkCompositeImageToPointCloudFilter.h"
 
#include <vtkPointData.h>
#include <vtkSmartPointer.h>
 
#include <QMap>
#include <QString>
#include <QVariant>
 
#include "VtkImageDataToPointCloudFilter.h"
 
VtkCompositeImageToPointCloudFilter::VtkCompositeImageToPointCloudFilter(
    vtkAlgorithm* inputAlgorithm)
    : VtkCompositeFilter(inputAlgorithm)
{
    init();
}
 
void VtkCompositeImageToPointCloudFilter::init()
{
    _inputDataObjectType = VTK_IMAGE_DATA;
    _outputDataObjectType = VTK_POLY_DATA;
 
    VtkImageDataToPointCloudFilter* point_cloud_filter =
        VtkImageDataToPointCloudFilter::New();
    point_cloud_filter->SetInputConnection(_inputAlgorithm->GetOutputPort());
    _inputAlgorithm->Update();
 
    QList<QVariant> n_points_range_list;
    n_points_range_list.push_back(
        point_cloud_filter->GetMinNumberOfPointsPerCell());
    n_points_range_list.push_back(
        point_cloud_filter->GetMaxNumberOfPointsPerCell());
    (*_algorithmUserVectorProperties)["Number of points range"] =
        n_points_range_list;
    QList<QVariant> vertical_extent_list;
    vertical_extent_list.push_back(point_cloud_filter->GetMinHeight());
    vertical_extent_list.push_back(point_cloud_filter->GetMaxHeight());
    (*_algorithmUserVectorProperties)["Vertical extent"] = vertical_extent_list;
    (*_algorithmUserProperties)["Logarithmic interpolation"] =
        !point_cloud_filter->GetIsLinear();
    (*_algorithmUserProperties)["Gamma value"] = point_cloud_filter->GetGamma();
 
    point_cloud_filter->Update();
    _outputAlgorithm = point_cloud_filter;
}
 
void VtkCompositeImageToPointCloudFilter::SetUserProperty(QString name,
                                                          QVariant value)
{
    VtkAlgorithmProperties::SetUserProperty(name, value);
 
    if ((name == "Gamma value") && (value.toDouble() > 0))
    {
        static_cast<VtkImageDataToPointCloudFilter*>(_outputAlgorithm)
            ->SetGamma(value.toDouble());
    }
    if (name == "Logarithmic interpolation")
    {
        if (value.toBool())
        {
            double const gamma =
                VtkAlgorithmProperties::GetUserProperty("Gamma value")
                    .toDouble();
            if (gamma > 0)
            {
                static_cast<VtkImageDataToPointCloudFilter*>(_outputAlgorithm)
                    ->useLogarithmicInterpolation(gamma);
            }
        }
        else
        {
            static_cast<VtkImageDataToPointCloudFilter*>(_outputAlgorithm)
                ->useLinearInterpolation();
        }
    }
}
 
void VtkCompositeImageToPointCloudFilter::SetUserVectorProperty(
    QString name, QList<QVariant> values)
{
    VtkAlgorithmProperties::SetUserVectorProperty(name, values);
 
    if (name == "Number of points range")
    {
        if (values[0].toInt() >= 0 && values[1].toInt() >= 0 &&
            values[0].toInt() <= values[1].toInt())
        {
            static_cast<VtkImageDataToPointCloudFilter*>(_outputAlgorithm)
                ->SetMinNumberOfPointsPerCell(values[0].toInt());
            static_cast<VtkImageDataToPointCloudFilter*>(_outputAlgorithm)
                ->SetMaxNumberOfPointsPerCell(values[1].toInt());
        }
    }
    else if (name == "Vertical extent")
    {
        if (values[0].toDouble() <= values[1].toDouble())
        {
            static_cast<VtkImageDataToPointCloudFilter*>(_outputAlgorithm)
                ->SetMinHeight(values[0].toDouble());
            static_cast<VtkImageDataToPointCloudFilter*>(_outputAlgorithm)
                ->SetMaxHeight(values[1].toDouble());
        }
    }
}