VtkStationSource.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
// ** VTK INCLUDES **
#include "VtkStationSource.h"
 
#include <vtkCellArray.h>
#include <vtkCellData.h>
#include <vtkDoubleArray.h>
#include <vtkInformation.h>
#include <vtkInformationVector.h>
#include <vtkLine.h>
#include <vtkPointData.h>
#include <vtkPoints.h>
#include <vtkPolyData.h>
#include <vtkProperty.h>
#include <vtkSmartPointer.h>
#include <vtkStreamingDemandDrivenPipeline.h>
 
#include "BaseLib/Logging.h"
#include "GeoLib/StationBorehole.h"
#include "vtkObjectFactory.h"
 
vtkStandardNewMacro(VtkStationSource);
 
VtkStationSource::VtkStationSource()
{
    _removable = false;  // From VtkAlgorithmProperties
    this->SetNumberOfInputPorts(0);
 
    const DataHolderLib::Color c = DataHolderLib::getRandomColor();
    GetProperties()->SetColor(c[0] / 255.0, c[1] / 255.0, c[2] / 255.0);
}
 
void VtkStationSource::PrintSelf(ostream& os, vtkIndent indent)
{
    this->Superclass::PrintSelf(os, indent);
 
    if (_stations->empty())
    {
        return;
    }
 
    os << indent << "== VtkStationSource =="
       << "\n";
 
    int i = 0;
    for (auto station : *_stations)
    {
        os << indent << "Station " << i << " (" << (*station)[0] << ", "
           << (*station)[1] << ", " << (*station)[2] << ")\n";
        i++;
    }
}

int VtkStationSource::RequestData(vtkInformation* request,
                                  vtkInformationVector** inputVector,
                                  vtkInformationVector* outputVector)
{
    (void)request;
    (void)inputVector;
 
    if (!_stations)
    {
        return 0;
    }
    std::size_t nStations = _stations->size();
    if (nStations == 0)
    {
        return 0;
    }
 
    bool useStationValues(false);
    double sValue =
        static_cast<GeoLib::Station*>((*_stations)[0])->getStationValue();
    for (std::size_t i = 1; i < nStations; i++)
    {
        if (static_cast<GeoLib::Station*>((*_stations)[i])->getStationValue() !=
            sValue)
        {
            useStationValues = true;
            break;
        }
    }
 
    bool isBorehole = dynamic_cast<GeoLib::StationBorehole*>((*_stations)[0]);
 
    vtkSmartPointer<vtkInformation> outInfo =
        outputVector->GetInformationObject(0);
    vtkSmartPointer<vtkPolyData> output =
        vtkPolyData::SafeDownCast(outInfo->Get(vtkDataObject::DATA_OBJECT()));
 
    vtkSmartPointer<vtkPoints> newStations = vtkSmartPointer<vtkPoints>::New();
    vtkSmartPointer<vtkCellArray> newVerts =
        vtkSmartPointer<vtkCellArray>::New();
    newVerts->Allocate(nStations);
 
    vtkSmartPointer<vtkCellArray> newLines;
 
    if (isBorehole)
    {
        newLines = vtkSmartPointer<vtkCellArray>::New();
    }
 
    if (outInfo->Get(vtkStreamingDemandDrivenPipeline::UPDATE_PIECE_NUMBER()) >
        0)
    {
        return 1;
    }
 
    vtkSmartPointer<vtkIntArray> station_ids =
        vtkSmartPointer<vtkIntArray>::New();
    station_ids->SetNumberOfComponents(1);
    station_ids->SetName("SiteIDs");
 
    vtkSmartPointer<vtkDoubleArray> station_values =
        vtkSmartPointer<vtkDoubleArray>::New();
    station_values->SetNumberOfComponents(1);
    station_values->SetName("StationValue");
 
    vtkSmartPointer<vtkIntArray> strat_ids =
        vtkSmartPointer<vtkIntArray>::New();
    strat_ids->SetNumberOfComponents(1);
    strat_ids->SetName("Stratigraphies");
 
    std::size_t lastMaxIndex(0);
    std::size_t site_count(0);
 
    // Generate graphic objects
    for (auto station : *_stations)
    {
        vtkIdType sid = newStations->InsertNextPoint(station->data());
        station_ids->InsertNextValue(site_count);
        if (useStationValues)
        {
            station_values->InsertNextValue(
                static_cast<GeoLib::Station*>(station)->getStationValue());
        }
 
        if (!isBorehole)
        {
            newVerts->InsertNextCell(1, &sid);
        }
        else
        {
            std::vector<GeoLib::Point*> profile =
                static_cast<GeoLib::StationBorehole*>(station)->getProfile();
            std::vector<std::string> soilNames =
                static_cast<GeoLib::StationBorehole*>(station)->getSoilNames();
            const std::size_t nLayers = profile.size();
 
            for (std::size_t i = 1; i < nLayers; i++)
            {
                newStations->InsertNextPoint(profile[i]->data());
                station_ids->InsertNextValue(site_count);
                newLines->InsertNextCell(2);
                newLines->InsertCellPoint(
                    lastMaxIndex);  // start of borehole-layer
                newLines->InsertCellPoint(
                    ++lastMaxIndex);  // end of boreholelayer
                strat_ids->InsertNextValue(this->GetIndexByName(soilNames[i]));
                if (useStationValues)
                {
                    station_values->InsertNextValue(
                        static_cast<GeoLib::Station*>(station)
                            ->getStationValue());
                }
            }
            lastMaxIndex++;
        }
        site_count++;
    }
 
    output->SetPoints(newStations);
 
    if (!isBorehole)
    {
        output->SetVerts(newVerts);
        output->GetCellData()->AddArray(station_ids);
        output->GetCellData()->SetActiveAttribute(
            "SiteIDs", vtkDataSetAttributes::SCALARS);
    }
    else
    {
        output->SetLines(newLines);
        // output->GetCellData()->AddArray(station_ids);
        output->GetCellData()->AddArray(strat_ids);
        output->GetCellData()->SetActiveAttribute(
            "Stratigraphies", vtkDataSetAttributes::SCALARS);
    }
    if (useStationValues)
    {
        output->GetPointData()->AddArray(station_values);
    }
 
    output->Squeeze();
 
    return 1;
}
 
int VtkStationSource::RequestInformation(vtkInformation* /*request*/,
                                         vtkInformationVector** /*inputVector*/,
                                         vtkInformationVector* /*outputVector*/)
{
    return 1;
}
 
void VtkStationSource::SetUserProperty(QString name, QVariant value)
{
    Q_UNUSED(name);
    Q_UNUSED(value);
}
 
std::size_t VtkStationSource::GetIndexByName(std::string const& name)
{
    vtkIdType max_key(0);
    for (auto& it : _id_map)
    {
        if (name == it.first)
        {
            return it.second;
        }
        if (it.second > max_key)
        {
            max_key = it.second;
        }
    }
 
    vtkIdType new_index = (_id_map.empty()) ? 0 : (max_key + 1);
    INFO("Key '{:s}' has been assigned index {:d}.", name, new_index);
    _id_map.insert(std::pair<std::string, vtkIdType>(name, new_index));
    return new_index;
}