void mitk::ConnectomicsNetworkCreator::CreateNetworkFromFibersAndSegmentation()
{
//empty graph
m_ConNetwork->clear();
m_LabelToVertexMap.clear();
m_LabelToNodePropertyMap.clear();
vtkSmartPointer<vtkPolyData> fiberPolyData = m_FiberBundle->GetFiberPolyData();
vtkSmartPointer<vtkCellArray> vLines = fiberPolyData->GetLines();
vLines->InitTraversal();
int numFibers = m_FiberBundle->GetNumFibers();
for( int fiberID( 0 ); fiberID < numFibers; fiberID++ )
{
vtkIdType numPointsInCell(0);
vtkIdType* pointsInCell(NULL);
vLines->GetNextCell ( numPointsInCell, pointsInCell );
TractType::Pointer singleTract = TractType::New();
for( int pointInCellID( 0 ); pointInCellID < numPointsInCell ; pointInCellID++)
{
// push back point
PointType point = GetItkPoint( fiberPolyData->GetPoint( pointsInCell[ pointInCellID ] ) );
singleTract->InsertElement( singleTract->Size(), point );
}
//MappingStrategy strategy = EndElementPosition;
//MappingStrategy strategy = JustEndPointVerticesNoLabel;
MappingStrategy strategy = EndElementPositionAvoidingWhiteMatter;
if ( singleTract && ( singleTract->Size() > 0 ) )
{
AddConnectionToNetwork(
ReturnAssociatedVertexPairForLabelPair(
ReturnLabelForFiberTract( singleTract, strategy )
)
);
}
}
// Prune unconnected nodes
m_ConNetwork->PruneUnconnectedSingleNodes();
// provide network with geometry
m_ConNetwork->SetGeometry( m_Segmentation->GetGeometry() );
m_ConNetwork->UpdateBounds();
m_ConNetwork->SetIsModified( true );
MBI_INFO << mitk::ConnectomicsConstantsManager::CONNECTOMICS_WARNING_INFO_NETWORK_CREATED;
}
void mitk::ConnectomicsNetworkCreator::CreateNetworkFromFibersAndSegmentation()
{
//empty graph
m_ConNetwork = mitk::ConnectomicsNetwork::New();
m_LabelToVertexMap.clear();
m_LabelToNodePropertyMap.clear();
idCounter = 0;
vtkSmartPointer<vtkPolyData> fiberPolyData = m_FiberBundle->GetFiberPolyData();
int numFibers = m_FiberBundle->GetNumFibers();
for( int fiberID( 0 ); fiberID < numFibers; fiberID++ )
{
vtkCell* cell = fiberPolyData->GetCell(fiberID);
int numPoints = cell->GetNumberOfPoints();
vtkPoints* points = cell->GetPoints();
TractType::Pointer singleTract = TractType::New();
for( int pointInCellID( 0 ); pointInCellID < numPoints ; pointInCellID++)
{
// push back point
PointType point = GetItkPoint( points->GetPoint( pointInCellID ) );
singleTract->InsertElement( singleTract->Size(), point );
}
if ( singleTract && ( singleTract->Size() > 0 ) )
{
AddConnectionToNetwork(
ReturnAssociatedVertexPairForLabelPair(
ReturnLabelForFiberTract( singleTract, m_MappingStrategy )
)
);
m_AbortConnection = false;
}
}
// Prune unconnected nodes
//m_ConNetwork->PruneUnconnectedSingleNodes();
// provide network with geometry
m_ConNetwork->SetGeometry( dynamic_cast<mitk::BaseGeometry*>(m_Segmentation->GetGeometry()->Clone().GetPointer()) );
m_ConNetwork->UpdateBounds();
m_ConNetwork->SetIsModified( true );
MBI_INFO << mitk::ConnectomicsConstantsManager::CONNECTOMICS_WARNING_INFO_NETWORK_CREATED;
}