networkobject.cpp

Go to the documentation of this file.

//=============================================================================================================
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "networkobject.h"
 
#include <connectivity/network/networknode.h>
#include <connectivity/network/networkedge.h>
#include <disp/plots/helpers/colormap.h>
 
#include <QQuaternion>
#include <QDebug>
#include <cmath>
 
using namespace CONNECTIVITYLIB;
using namespace DISPLIB;
using namespace Eigen;
 
namespace DISP3DLIB
{
 
//=============================================================================================================
// DEFINE MEMBER METHODS
//=============================================================================================================
 
NetworkObject::NetworkObject() = default;
 
//=============================================================================================================
 
NetworkObject::~NetworkObject() = default;
 
//=============================================================================================================
 
InstancedStreams NetworkObject::nodeStreams() const
{
    return {{m_nodeVertexData.constData(), static_cast<quint32>(m_nodeVertexData.size()), m_nodeGeometryGeneration},
            {m_nodeIndexData.constData(), static_cast<quint32>(m_nodeIndexData.size()), m_nodeGeometryGeneration},
            {m_nodeInstanceData.constData(), static_cast<quint32>(m_nodeInstanceCount * sizeof(InstanceData)), m_nodeInstanceGeneration},
            m_nodeIndexCount,
            m_nodeInstanceCount};
}
 
InstancedStreams NetworkObject::edgeStreams() const
{
    return {{m_edgeVertexData.constData(), static_cast<quint32>(m_edgeVertexData.size()), m_edgeGeometryGeneration},
            {m_edgeIndexData.constData(), static_cast<quint32>(m_edgeIndexData.size()), m_edgeGeometryGeneration},
            {m_edgeInstanceData.constData(), static_cast<quint32>(m_edgeInstanceCount * sizeof(InstanceData)), m_edgeInstanceGeneration},
            m_edgeIndexCount,
            m_edgeInstanceCount};
}
 
//=============================================================================================================
 
void NetworkObject::load(const Network& network, const QString& sColormap)
{
    m_network = network;
    m_colormap = sColormap;
 
    createNodeGeometry();
    createEdgeGeometry();
    buildNodeInstances();
    buildEdgeInstances();
}
 
//=============================================================================================================
 
void NetworkObject::setThreshold(double dThreshold)
{
    m_network.setThreshold(dThreshold);
    buildNodeInstances();
    buildEdgeInstances();
}
 
//=============================================================================================================
 
void NetworkObject::setColormap(const QString& sColormap)
{
    m_colormap = sColormap;
    buildNodeInstances();
    buildEdgeInstances();
}
 
//=============================================================================================================
 
void NetworkObject::createNodeGeometry()
{
    if (!m_nodeVertexData.isEmpty())
        return;
 
    // Create an icosphere (subdivision level 1) for nodes
    const int subdivisions = 1;
    const float radius = 1.0f; // Unit sphere, scaled per-instance
 
    // Start with icosahedron
    const float t = (1.0f + std::sqrt(5.0f)) / 2.0f;
 
    std::vector<QVector3D> vertices = {
        QVector3D(-1, t, 0).normalized() * radius,
        QVector3D(1, t, 0).normalized() * radius,
        QVector3D(-1, -t, 0).normalized() * radius,
        QVector3D(1, -t, 0).normalized() * radius,
        QVector3D(0, -1, t).normalized() * radius,
        QVector3D(0, 1, t).normalized() * radius,
        QVector3D(0, -1, -t).normalized() * radius,
        QVector3D(0, 1, -t).normalized() * radius,
        QVector3D(t, 0, -1).normalized() * radius,
        QVector3D(t, 0, 1).normalized() * radius,
        QVector3D(-t, 0, -1).normalized() * radius,
        QVector3D(-t, 0, 1).normalized() * radius,
    };
 
    std::vector<uint32_t> indices = {
        0,
        11,
        5,
        0,
        5,
        1,
        0,
        1,
        7,
        0,
        7,
        10,
        0,
        10,
        11,
        1,
        5,
        9,
        5,
        11,
        4,
        11,
        10,
        2,
        10,
        7,
        6,
        7,
        1,
        8,
        3,
        9,
        4,
        3,
        4,
        2,
        3,
        2,
        6,
        3,
        6,
        8,
        3,
        8,
        9,
        4,
        9,
        5,
        2,
        4,
        11,
        6,
        2,
        10,
        8,
        6,
        7,
        9,
        8,
        1,
    };
 
    // Subdivide
    for (int s = 0; s < subdivisions; ++s) {
        std::vector<uint32_t> newIndices;
        std::map<uint64_t, uint32_t> midpointCache;
 
        auto getMidpoint = [&](uint32_t i0, uint32_t i1) -> uint32_t {
            uint64_t key = static_cast<uint64_t>(std::min(i0, i1)) << 32 | std::max(i0, i1);
            auto it = midpointCache.find(key);
            if (it != midpointCache.end())
                return it->second;
 
            QVector3D mid = ((vertices[i0] + vertices[i1]) / 2.0f).normalized() * radius;
            uint32_t idx = static_cast<uint32_t>(vertices.size());
            vertices.push_back(mid);
            midpointCache[key] = idx;
            return idx;
        };
 
        for (size_t i = 0; i < indices.size(); i += 3) {
            uint32_t a = indices[i], b = indices[i + 1], c = indices[i + 2];
            uint32_t ab = getMidpoint(a, b);
            uint32_t bc = getMidpoint(b, c);
            uint32_t ca = getMidpoint(c, a);
 
            newIndices.insert(newIndices.end(), {a, ab, ca});
            newIndices.insert(newIndices.end(), {b, bc, ab});
            newIndices.insert(newIndices.end(), {c, ca, bc});
            newIndices.insert(newIndices.end(), {ab, bc, ca});
        }
 
        indices = std::move(newIndices);
    }
 
    // Build vertex data with normals (normal = normalized position for sphere)
    std::vector<VertexData> vd;
    vd.reserve(vertices.size());
    for (const auto& v : vertices) {
        QVector3D n = v.normalized();
        vd.push_back({v.x(), v.y(), v.z(), n.x(), n.y(), n.z()});
    }
 
    m_nodeIndexCount = static_cast<int>(indices.size());
 
    m_nodeVertexData.resize(vd.size() * sizeof(VertexData));
    memcpy(m_nodeVertexData.data(), vd.data(), m_nodeVertexData.size());
 
    m_nodeIndexData.resize(indices.size() * sizeof(uint32_t));
    memcpy(m_nodeIndexData.data(), indices.data(), m_nodeIndexData.size());
 
    ++m_nodeGeometryGeneration;
}
 
//=============================================================================================================
 
void NetworkObject::createEdgeGeometry()
{
    if (!m_edgeVertexData.isEmpty())
        return;
 
    // Create a unit cylinder along Y axis (height=1, radius=1, scaled per-instance)
    const int segments = 8;
    const float radius = 1.0f;
    const float halfHeight = 0.5f;
 
    std::vector<VertexData> vertices;
    std::vector<uint32_t> indices;
 
    // Top center (0)
    vertices.push_back({0, halfHeight, 0, 0, 1, 0});
    // Bottom center (1)
    vertices.push_back({0, -halfHeight, 0, 0, -1, 0});
 
    // Side vertices: top ring (2..2+segments-1), bottom ring (2+segments..2+2*segments-1)
    for (int i = 0; i < segments; ++i) {
        float angle = 2.0f * static_cast<float>(M_PI) * i / segments;
        float x = radius * std::cos(angle);
        float z = radius * std::sin(angle);
 
        QVector3D normal(x, 0, z);
        normal.normalize();
 
        // Top side vertex
        vertices.push_back({x, halfHeight, z, normal.x(), normal.y(), normal.z()});
        // Bottom side vertex
        vertices.push_back({x, -halfHeight, z, normal.x(), normal.y(), normal.z()});
    }
 
    // Top cap vertices (for proper normals)
    int topCapStart = static_cast<int>(vertices.size());
    for (int i = 0; i < segments; ++i) {
        float angle = 2.0f * static_cast<float>(M_PI) * i / segments;
        float x = radius * std::cos(angle);
        float z = radius * std::sin(angle);
        vertices.push_back({x, halfHeight, z, 0, 1, 0});
    }
 
    // Bottom cap vertices
    int botCapStart = static_cast<int>(vertices.size());
    for (int i = 0; i < segments; ++i) {
        float angle = 2.0f * static_cast<float>(M_PI) * i / segments;
        float x = radius * std::cos(angle);
        float z = radius * std::sin(angle);
        vertices.push_back({x, -halfHeight, z, 0, -1, 0});
    }
 
    // Side faces
    int sideStart = 2;
    for (int i = 0; i < segments; ++i) {
        int next = (i + 1) % segments;
        int topCur = sideStart + i * 2;
        int botCur = sideStart + i * 2 + 1;
        int topNext = sideStart + next * 2;
        int botNext = sideStart + next * 2 + 1;
 
        indices.insert(indices.end(), {static_cast<uint32_t>(topCur), static_cast<uint32_t>(topNext), static_cast<uint32_t>(botCur)});
        indices.insert(indices.end(), {static_cast<uint32_t>(botCur), static_cast<uint32_t>(topNext), static_cast<uint32_t>(botNext)});
    }
 
    // Top cap
    for (int i = 0; i < segments; ++i) {
        int next = (i + 1) % segments;
        indices.push_back(0); // center
        indices.push_back(topCapStart + i);
        indices.push_back(topCapStart + next);
    }
 
    // Bottom cap
    for (int i = 0; i < segments; ++i) {
        int next = (i + 1) % segments;
        indices.push_back(1); // center
        indices.push_back(botCapStart + next);
        indices.push_back(botCapStart + i);
    }
 
    m_edgeIndexCount = static_cast<int>(indices.size());
 
    m_edgeVertexData.resize(vertices.size() * sizeof(VertexData));
    memcpy(m_edgeVertexData.data(), vertices.data(), m_edgeVertexData.size());
 
    m_edgeIndexData.resize(indices.size() * sizeof(uint32_t));
    memcpy(m_edgeIndexData.data(), indices.data(), m_edgeIndexData.size());
 
    ++m_edgeGeometryGeneration;
}
 
//=============================================================================================================
 
void NetworkObject::buildNodeInstances()
{
    if (m_network.isEmpty()) {
        m_nodeInstanceCount = 0;
        ++m_nodeInstanceGeneration;
        return;
    }
 
    const auto& nodes = m_network.getNodes();
    qint16 iMaxDegree = m_network.getMinMaxThresholdedDegrees().second;
    if (iMaxDegree == 0)
        iMaxDegree = 1;
 
    VisualizationInfo vizInfo = m_network.getVisualizationInfo();
 
    std::vector<InstanceData> instances;
    instances.reserve(nodes.size());
 
    for (int i = 0; i < nodes.size(); ++i) {
        qint16 degree = nodes[i]->getThresholdedDegree();
        if (degree == 0)
            continue;
 
        const RowVectorXf& vert = nodes[i]->getVert();
        QVector3D pos(vert(0), vert(1), vert(2));
 
        // Scale: nodes with higher degree are larger
        // Range: 0.0006 to 0.005 (same as disp3D)
        float scaleFactor = (static_cast<float>(degree) / static_cast<float>(iMaxDegree)) * (0.005f - 0.0006f) + 0.0006f;
 
        QMatrix4x4 m;
        m.translate(pos);
        m.scale(scaleFactor);
 
        InstanceData inst;
        const float* mPtr = m.constData();
        for (int j = 0; j < 16; ++j)
            inst.model[j] = mPtr[j];
 
        // Color: colormap-based or fixed
        if (vizInfo.sMethod == "Map") {
            float normalized = static_cast<float>(degree) / static_cast<float>(iMaxDegree);
            QRgb rgb = ColorMap::valueToColor(normalized, vizInfo.sColormap.isEmpty() ? m_colormap : vizInfo.sColormap);
            QColor color(rgb);
            float alpha = std::pow(normalized, 4.0f); // Same as disp3D
            inst.color[0] = color.redF();
            inst.color[1] = color.greenF();
            inst.color[2] = color.blueF();
            inst.color[3] = alpha;
        } else {
            inst.color[0] = vizInfo.colNodes[0] / 255.0f;
            inst.color[1] = vizInfo.colNodes[1] / 255.0f;
            inst.color[2] = vizInfo.colNodes[2] / 255.0f;
            inst.color[3] = vizInfo.colNodes[3] / 255.0f;
        }
        inst.isSelected = 0.0f;
 
        instances.push_back(inst);
    }
 
    m_nodeInstanceCount = static_cast<int>(instances.size());
    m_nodeInstanceData.resize(m_nodeInstanceCount * sizeof(InstanceData));
    if (m_nodeInstanceCount > 0) {
        memcpy(m_nodeInstanceData.data(), instances.data(), m_nodeInstanceData.size());
    }
    ++m_nodeInstanceGeneration;
 
    qDebug() << "NetworkObject: Built" << m_nodeInstanceCount << "node instances";
}
 
//=============================================================================================================
 
void NetworkObject::buildEdgeInstances()
{
    if (m_network.isEmpty()) {
        m_edgeInstanceCount = 0;
        ++m_edgeInstanceGeneration;
        return;
    }
 
    const auto& edges = m_network.getThresholdedEdges();
    const auto& nodes = m_network.getNodes();
 
    double dMaxWeight = m_network.getMinMaxThresholdedWeights().second;
    double dMinWeight = m_network.getMinMaxThresholdedWeights().first;
    double dWeightRange = dMaxWeight - dMinWeight;
    if (dWeightRange == 0.0)
        dWeightRange = 1.0;
 
    VisualizationInfo vizInfo = m_network.getVisualizationInfo();
 
    std::vector<InstanceData> instances;
    instances.reserve(edges.size());
 
    for (int i = 0; i < edges.size(); ++i) {
        auto& edge = edges[i];
        if (!edge->isActive())
            continue;
 
        int iStart = edge->getStartNodeID();
        int iEnd = edge->getEndNodeID();
 
        if (iStart < 0 || iStart >= nodes.size() || iEnd < 0 || iEnd >= nodes.size())
            continue;
 
        const RowVectorXf& vStart = nodes[iStart]->getVert();
        const RowVectorXf& vEnd = nodes[iEnd]->getVert();
 
        QVector3D startPos(vStart(0), vStart(1), vStart(2));
        QVector3D endPos(vEnd(0), vEnd(1), vEnd(2));
 
        if (startPos == endPos)
            continue;
 
        double dWeight = std::fabs(edge->getWeight());
        if (dWeight == 0.0)
            continue;
 
        QVector3D diff = endPos - startPos;
        QVector3D midPoint = startPos + diff / 2.0f;
        float length = diff.length();
 
        // Build transform: translate to midpoint, rotate Y-axis to diff direction, scale
        float normalizedWeight = static_cast<float>(std::fabs((dWeight - dMinWeight) / dWeightRange));
 
        // Cylinder radius: proportional to weight, range 0.0001 to 0.001
        float edgeRadius = 0.0001f + normalizedWeight * 0.0009f;
 
        QMatrix4x4 m;
        m.translate(midPoint);
        m.rotate(QQuaternion::rotationTo(QVector3D(0, 1, 0), diff.normalized()));
        m.scale(edgeRadius, length, edgeRadius);
 
        InstanceData inst;
        const float* mPtr = m.constData();
        for (int j = 0; j < 16; ++j)
            inst.model[j] = mPtr[j];
 
        // Color
        if (vizInfo.sMethod == "Map") {
            float normalized = (dMaxWeight != 0.0) ? static_cast<float>(std::fabs(dWeight / dMaxWeight)) : 0.0f;
            QRgb rgb = ColorMap::valueToColor(normalized, vizInfo.sColormap.isEmpty() ? m_colormap : vizInfo.sColormap);
            QColor color(rgb);
            float alpha = std::pow(normalized, 1.5f); // Same as disp3D
            inst.color[0] = color.redF();
            inst.color[1] = color.greenF();
            inst.color[2] = color.blueF();
            inst.color[3] = alpha;
        } else {
            inst.color[0] = vizInfo.colEdges[0] / 255.0f;
            inst.color[1] = vizInfo.colEdges[1] / 255.0f;
            inst.color[2] = vizInfo.colEdges[2] / 255.0f;
            inst.color[3] = vizInfo.colEdges[3] / 255.0f;
        }
        inst.isSelected = 0.0f;
 
        instances.push_back(inst);
    }
 
    m_edgeInstanceCount = static_cast<int>(instances.size());
    m_edgeInstanceData.resize(m_edgeInstanceCount * sizeof(InstanceData));
    if (m_edgeInstanceCount > 0) {
        memcpy(m_edgeInstanceData.data(), instances.data(), m_edgeInstanceData.size());
    }
    ++m_edgeInstanceGeneration;
 
    qDebug() << "NetworkObject: Built" << m_edgeInstanceCount << "edge instances";
}
 
//=============================================================================================================
 
} // namespace DISP3DLIB