diff --git a/Modules/DiffusionImaging/Rendering/mitkConnectomicsNetworkMapper3D.cpp b/Modules/DiffusionImaging/Rendering/mitkConnectomicsNetworkMapper3D.cpp index 13eaa56edf..145701b316 100644 --- a/Modules/DiffusionImaging/Rendering/mitkConnectomicsNetworkMapper3D.cpp +++ b/Modules/DiffusionImaging/Rendering/mitkConnectomicsNetworkMapper3D.cpp @@ -1,181 +1,263 @@ /*=================================================================== The Medical Imaging Interaction Toolkit (MITK) Copyright (c) German Cancer Research Center, Division of Medical and Biological Informatics. All rights reserved. This software is distributed WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See LICENSE.txt or http://www.mitk.org for details. ===================================================================*/ #include "mitkConnectomicsNetworkMapper3D.h" +#include +#include "vtkGraphLayout.h" +#include +#include "vtkGraphToPolyData.h" +#include +#include "vtkGlyph3D.h" +#include "vtkGlyphSource2D.h" + mitk::ConnectomicsNetworkMapper3D::ConnectomicsNetworkMapper3D() { //TODO: implement m_NetworkAssembly = vtkPropAssembly::New(); } mitk::ConnectomicsNetworkMapper3D:: ~ConnectomicsNetworkMapper3D() { //TODO: implement m_NetworkAssembly->Delete(); } void mitk::ConnectomicsNetworkMapper3D::GenerateDataForRenderer(mitk::BaseRenderer* renderer) { //TODO: implement if( this->GetInput() == NULL ) { return; } if( this->GetInput()->GetIsModified( ) ) { GenerateData(); } } void mitk::ConnectomicsNetworkMapper3D::GenerateData() { m_NetworkAssembly->Delete(); m_NetworkAssembly = vtkPropAssembly::New(); // Here is the part where a graph is given and converted to points and connections between points... std::vector< mitk::ConnectomicsNetwork::NetworkNode > vectorOfNodes = this->GetInput()->GetVectorOfAllNodes(); std::vector< std::pair< std::pair< mitk::ConnectomicsNetwork::NetworkNode, mitk::ConnectomicsNetwork::NetworkNode > , mitk::ConnectomicsNetwork::NetworkEdge > > vectorOfEdges = this->GetInput()->GetVectorOfAllEdges(); - mitk::Point3D tempWorldPoint, tempCNFGeometryPoint; - - //////////////////////Create Spheres///////////////////////// - for(unsigned int i = 0; i < vectorOfNodes.size(); i++) + // Decide on the style of rendering due to property + if( false ) { - vtkSmartPointer sphereSource = - vtkSmartPointer::New(); + mitk::Point3D tempWorldPoint, tempCNFGeometryPoint; - for(unsigned int dimension = 0; dimension < 3; dimension++) + //////////////////////Create Spheres///////////////////////// + for(unsigned int i = 0; i < vectorOfNodes.size(); i++) { - tempCNFGeometryPoint.SetElement( dimension , vectorOfNodes[i].coordinates[dimension] ); - } + vtkSmartPointer sphereSource = + vtkSmartPointer::New(); - this->GetData()->GetGeometry()->IndexToWorld( tempCNFGeometryPoint, tempWorldPoint ); + for(unsigned int dimension = 0; dimension < 3; dimension++) + { + tempCNFGeometryPoint.SetElement( dimension , vectorOfNodes[i].coordinates[dimension] ); + } - sphereSource->SetCenter( tempWorldPoint[0] , tempWorldPoint[1], tempWorldPoint[2] ); - sphereSource->SetRadius(1.0); + this->GetData()->GetGeometry()->IndexToWorld( tempCNFGeometryPoint, tempWorldPoint ); - vtkSmartPointer mapper = - vtkSmartPointer::New(); - mapper->SetInput(sphereSource->GetOutput()); + sphereSource->SetCenter( tempWorldPoint[0] , tempWorldPoint[1], tempWorldPoint[2] ); + sphereSource->SetRadius(1.0); - vtkSmartPointer actor = - vtkSmartPointer::New(); - actor->SetMapper(mapper); + vtkSmartPointer mapper = + vtkSmartPointer::New(); + mapper->SetInput(sphereSource->GetOutput()); - m_NetworkAssembly->AddPart(actor); - } + vtkSmartPointer actor = + vtkSmartPointer::New(); + actor->SetMapper(mapper); - //////////////////////Create Tubes///////////////////////// - - for(unsigned int i = 0; i < vectorOfEdges.size(); i++) - { + m_NetworkAssembly->AddPart(actor); + } - vtkSmartPointer lineSource = - vtkSmartPointer::New(); + //////////////////////Create Tubes///////////////////////// - for(unsigned int dimension = 0; dimension < 3; dimension++) + for(unsigned int i = 0; i < vectorOfEdges.size(); i++) { - tempCNFGeometryPoint[ dimension ] = vectorOfEdges[i].first.first.coordinates[dimension]; - } - this->GetData()->GetGeometry()->IndexToWorld( tempCNFGeometryPoint, tempWorldPoint ); + vtkSmartPointer lineSource = + vtkSmartPointer::New(); - lineSource->SetPoint1(tempWorldPoint[0], tempWorldPoint[1],tempWorldPoint[2] ); + for(unsigned int dimension = 0; dimension < 3; dimension++) + { + tempCNFGeometryPoint[ dimension ] = vectorOfEdges[i].first.first.coordinates[dimension]; + } - for(unsigned int dimension = 0; dimension < 3; dimension++) - { - tempCNFGeometryPoint[ dimension ] = vectorOfEdges[i].first.second.coordinates[dimension]; - } + this->GetData()->GetGeometry()->IndexToWorld( tempCNFGeometryPoint, tempWorldPoint ); - this->GetData()->GetGeometry()->IndexToWorld( tempCNFGeometryPoint, tempWorldPoint ); + lineSource->SetPoint1(tempWorldPoint[0], tempWorldPoint[1],tempWorldPoint[2] ); - lineSource->SetPoint2(tempWorldPoint[0], tempWorldPoint[1], tempWorldPoint[2] ); + for(unsigned int dimension = 0; dimension < 3; dimension++) + { + tempCNFGeometryPoint[ dimension ] = vectorOfEdges[i].first.second.coordinates[dimension]; + } - vtkSmartPointer tubes = vtkSmartPointer::New(); - tubes->SetInput( lineSource->GetOutput() ); - tubes->SetNumberOfSides( 12 ); - double radiusFactor = 1.0 + ((double) vectorOfEdges[i].second.weight) / 10.0 ; - tubes->SetRadius( std::log10( radiusFactor ) ); + this->GetData()->GetGeometry()->IndexToWorld( tempCNFGeometryPoint, tempWorldPoint ); - vtkSmartPointer mapper2 = - vtkSmartPointer::New(); - mapper2->SetInput( tubes->GetOutput() ); + lineSource->SetPoint2(tempWorldPoint[0], tempWorldPoint[1], tempWorldPoint[2] ); - vtkSmartPointer actor = - vtkSmartPointer::New(); - actor->SetMapper(mapper2); + vtkSmartPointer tubes = vtkSmartPointer::New(); + tubes->SetInput( lineSource->GetOutput() ); + tubes->SetNumberOfSides( 12 ); + double radiusFactor = 1.0 + ((double) vectorOfEdges[i].second.weight) / 10.0 ; + tubes->SetRadius( std::log10( radiusFactor ) ); - double maxWeight = (double) this->GetInput()->GetMaximumWeight(); - double colourFactor = vectorOfEdges[i].second.weight / maxWeight; - actor->GetProperty()->SetColor( colourFactor, colourFactor, colourFactor); + vtkSmartPointer mapper2 = + vtkSmartPointer::New(); + mapper2->SetInput( tubes->GetOutput() ); + vtkSmartPointer actor = + vtkSmartPointer::New(); + actor->SetMapper(mapper2); - m_NetworkAssembly->AddPart(actor); + double maxWeight = (double) this->GetInput()->GetMaximumWeight(); + double colourFactor = vectorOfEdges[i].second.weight / maxWeight; + actor->GetProperty()->SetColor( colourFactor, colourFactor, colourFactor); + + m_NetworkAssembly->AddPart(actor); + + } } + else + { + vtkSmartPointer graph = + vtkSmartPointer::New(); - (static_cast ( GetData() ) )->SetIsModified( false ); + std::vector< vtkIdType > networkToVTKvector; + networkToVTKvector.resize(vectorOfNodes.size()); - //this->UpdateVtkObjects(); + for(unsigned int i = 0; i < vectorOfNodes.size(); i++) + { + networkToVTKvector[vectorOfNodes[i].id] = graph->AddVertex(); + } + for(unsigned int i = 0; i < vectorOfEdges.size(); i++) + { + graph->AddEdge(networkToVTKvector[vectorOfEdges[i].first.first.id], networkToVTKvector[vectorOfEdges[i].first.second.id]); + } + vtkSmartPointer points = + vtkSmartPointer::New(); + for(unsigned int i = 0; i < vectorOfNodes.size(); i++) + { + double x = vectorOfNodes[i].coordinates[0]; + double y = vectorOfNodes[i].coordinates[1]; + double z = vectorOfNodes[i].coordinates[2]; + points->InsertNextPoint( x, y, z); + } + + graph->SetPoints(points); + + vtkGraphLayout* layout = vtkGraphLayout::New(); + layout->SetInput(graph); + layout->SetLayoutStrategy(vtkPassThroughLayoutStrategy::New()); + + vtkGraphToPolyData* graphToPoly = vtkGraphToPolyData::New(); + graphToPoly->SetInputConnection(layout->GetOutputPort()); + + // Create the standard VTK polydata mapper and actor + // for the connections (edges) in the tree. + vtkPolyDataMapper* edgeMapper = vtkPolyDataMapper::New(); + edgeMapper->SetInputConnection(graphToPoly->GetOutputPort()); + vtkActor* edgeActor = vtkActor::New(); + edgeActor->SetMapper(edgeMapper); + edgeActor->GetProperty()->SetColor(0.0, 0.5, 1.0); + + // Glyph the points of the tree polydata to create + // VTK_VERTEX cells at each vertex in the tree. + vtkGlyph3D* vertGlyph = vtkGlyph3D::New(); + vertGlyph->SetInputConnection(0, graphToPoly->GetOutputPort()); + vtkGlyphSource2D* glyphSource = vtkGlyphSource2D::New(); + glyphSource->SetGlyphTypeToVertex(); + vertGlyph->SetInputConnection(1, glyphSource->GetOutputPort()); + + // Create a mapper for the vertices, and tell the mapper + // to use the specified color array. + vtkPolyDataMapper* vertMapper = vtkPolyDataMapper::New(); + vertMapper->SetInputConnection(vertGlyph->GetOutputPort()); + /*if (colorArray) + { + vertMapper->SetScalarModeToUsePointFieldData(); + vertMapper->SelectColorArray(colorArray); + vertMapper->SetScalarRange(colorRange); + }*/ + + // Create an actor for the vertices. Move the actor forward + // in the z direction so it is drawn on top of the edge actor. + vtkActor* vertActor = vtkActor::New(); + vertActor->SetMapper(vertMapper); + vertActor->GetProperty()->SetPointSize(5); + vertActor->SetPosition(0, 0, 0.001); + + m_NetworkAssembly->AddPart(edgeActor); + m_NetworkAssembly->AddPart(vertActor); + } + + (static_cast ( GetData() ) )->SetIsModified( false ); } const mitk::ConnectomicsNetwork* mitk::ConnectomicsNetworkMapper3D::GetInput() { return static_cast ( GetData() ); } void mitk::ConnectomicsNetworkMapper3D::SetDefaultProperties(DataNode* node, BaseRenderer* renderer , bool overwrite) { //TODO: implement // hand it to the superclass for base default properties Superclass::SetDefaultProperties(node, renderer, overwrite); } void mitk::ConnectomicsNetworkMapper3D::ApplyProperties(mitk::BaseRenderer* renderer) { //TODO: implement } void mitk::ConnectomicsNetworkMapper3D::SetVtkMapperImmediateModeRendering(vtkMapper *mapper) { //TODO: implement } void mitk::ConnectomicsNetworkMapper3D::UpdateVtkObjects() { //TODO: implement } vtkProp* mitk::ConnectomicsNetworkMapper3D::GetVtkProp(mitk::BaseRenderer *renderer) { return m_NetworkAssembly; }