1.連通區(qū)域分析

許多圖形數(shù)據(jù)中，并非只包含一個(gè)對(duì)象（連通區(qū)域）。而在處理這些圖形數(shù)據(jù)時(shí)，有時(shí)需要對(duì)每一個(gè)對(duì)象單獨(dú)處理或者讓其單獨(dú)顯示。比如，利用MarchingCube方法提取三維圖像中的等值面，得到的結(jié)果往往是存在多個(gè)連通的對(duì)象區(qū)域，這是就需要對(duì)圖形數(shù)據(jù)做連通區(qū)域分析，提取每個(gè)連通區(qū)域并計(jì)算其屬性信息，以此來(lái)得到需要的連通區(qū)域。

下面一個(gè)例子來(lái)分析VTK中如何對(duì)圖形數(shù)據(jù)做連通區(qū)域分析：

#include <vtkAutoInit.h> VTK_MODULE_INIT(vtkRenderingOpenGL); VTK_MODULE_INIT(vtkRenderingFreeType); VTK_MODULE_INIT(vtkInteractionStyle);#include <vtkSmartPointer.h> #include <vtkSphereSource.h> #include <vtkConeSource.h> #include <vtkAppendPolyData.h> #include <vtkPolyDataConnectivityFilter.h> #include <vtkPolyDataMapper.h> #include <vtkActor.h> #include <vtkProperty.h> #include <vtkRenderer.h> #include <vtkRenderWindow.h> #include <vtkRenderWindowInteractor.h> int main() {vtkSmartPointer<vtkSphereSource> sphereSource =vtkSmartPointer<vtkSphereSource>::New();sphereSource->SetRadius(10);sphereSource->SetThetaResolution(10);sphereSource->SetPhiResolution(10);sphereSource->Update();vtkSmartPointer<vtkConeSource> coneSource =vtkSmartPointer<vtkConeSource>::New();coneSource->SetRadius(5);coneSource->SetHeight(10);coneSource->SetCenter(25, 0, 0);coneSource->Update();vtkSmartPointer<vtkAppendPolyData> appendFilter =vtkSmartPointer<vtkAppendPolyData>::New();appendFilter->AddInputData(sphereSource->GetOutput());appendFilter->AddInputData(coneSource->GetOutput());appendFilter->Update();vtkSmartPointer<vtkPolyDataConnectivityFilter> connectivityFilter =vtkSmartPointer<vtkPolyDataConnectivityFilter>::New();connectivityFilter->SetInputData(appendFilter->GetOutput());connectivityFilter->SetExtractionModeToCellSeededRegions();connectivityFilter->AddSeed(100);connectivityFilter->Update();vtkSmartPointer<vtkPolyDataMapper> originalMapper =vtkSmartPointer<vtkPolyDataMapper>::New();originalMapper->SetInputConnection(appendFilter->GetOutputPort());originalMapper->Update();vtkSmartPointer<vtkActor> originalActor =vtkSmartPointer<vtkActor>::New();originalActor->SetMapper(originalMapper);vtkSmartPointer<vtkPolyDataMapper> extractedMapper =vtkSmartPointer<vtkPolyDataMapper>::New();extractedMapper->SetInputConnection(connectivityFilter->GetOutputPort());extractedMapper->Update();vtkSmartPointer<vtkActor> extractedActor =vtkSmartPointer<vtkActor>::New();extractedActor->SetMapper(extractedMapper);/double leftViewport[4] = { 0.0, 0.0, 0.5, 1.0 };double rightViewport[4] = { 0.5, 0.0, 1.0, 1.0 };vtkSmartPointer<vtkRenderer> leftRenderer =vtkSmartPointer<vtkRenderer>::New();leftRenderer->SetViewport(leftViewport);leftRenderer->AddActor(originalActor);leftRenderer->SetBackground(1, 0, 0);vtkSmartPointer<vtkRenderer> rightRenderer =vtkSmartPointer<vtkRenderer>::New();rightRenderer->SetViewport(rightViewport);rightRenderer->AddActor(extractedActor);rightRenderer->SetBackground(0, 0, 0);vtkSmartPointer<vtkRenderWindow> renderWindow =vtkSmartPointer<vtkRenderWindow>::New();renderWindow->AddRenderer(leftRenderer);renderWindow->AddRenderer(rightRenderer);renderWindow->SetSize(640, 320);renderWindow->Render();renderWindow->SetWindowName("PolyDataConnectedCompExtract");leftRenderer->ResetCamera();rightRenderer->SetActiveCamera(leftRenderer->GetActiveCamera());vtkSmartPointer<vtkRenderWindowInteractor> interactor =vtkSmartPointer<vtkRenderWindowInteractor>::New();interactor->SetRenderWindow(renderWindow);interactor->Initialize();interactor->Start();return 0; }輸出結(jié)果為：

這個(gè)例子夠早了一個(gè)含有多個(gè)連通區(qū)域的模型數(shù)據(jù)。vtkAppendPolyData可以實(shí)現(xiàn)vtkPolyData的合并，使用該類可以方便地構(gòu)造含有多個(gè)連通區(qū)域的數(shù)據(jù)，該類型接收兩個(gè)或者多個(gè)vtkPolyData數(shù)據(jù)輸入，合并結(jié)果包含輸入數(shù)據(jù)的所有幾何和拓?fù)鋽?shù)據(jù)。若輸入為兩個(gè)或者多個(gè)數(shù)據(jù)都含有點(diǎn)屬性數(shù)據(jù)，則將其存儲(chǔ)值輸出結(jié)果中；對(duì)于單元屬性數(shù)據(jù)亦是如此。

2.VTKPolyDataConnectivityFilter類解析

vtk中的vtkPolyDataConnectivityFilter類可以用于實(shí)現(xiàn)連通區(qū)域分析，該類接受vtkPolyData數(shù)據(jù)作為輸入。集體使用如下： vtkSmartPointer<vtkPolyDataConnectivityFilter> connectivityFilter =vtkSmartPointer<vtkPolyDataConnectivityFilter>::New();connectivityFilter->SetInputData(appendFilter->GetOutput());connectivityFilter->SetExtractionModeToCellSeededRegions();connectivityFilter->AddSeed(100);connectivityFilter->Update();SetExtractionModeToLargestRegion()：用于提取具有最多點(diǎn)數(shù)的連通區(qū)域； SetExtractionModeToAllRegions()：該模式主要用于連通區(qū)域標(biāo)記，配合函數(shù)ColorRegionsOn()使用，在連通區(qū)域像是的同時(shí)，生成一個(gè)名為RegionId的點(diǎn)屬性數(shù)據(jù)。
SetExtractionModeToSpecifiedRegions()：該模式用于提取一個(gè)或多個(gè)連通區(qū)域，在該模式下，需要通過AddSpecifiedRegion()來(lái)添加西藥提取的區(qū)域號(hào)，區(qū)域號(hào)從零開始。 SetExtractionModeToClosestPointRegion()：該模式需要使用SetClosestPoint()函數(shù)設(shè)置一個(gè)空間點(diǎn)坐標(biāo)，執(zhí)行結(jié)果為離該點(diǎn)最近的連通區(qū)域。 SetExtractionModeToPointSeededRegions()：該模式下需要使用AddSeed()函數(shù)添加種子點(diǎn)，提取種子點(diǎn)所在的區(qū)域。 SetExtractionModeToCellSeededRegions()：該模式下需要使用AddSeed()函數(shù)調(diào)價(jià)種子單元，提取種子單元所在的區(qū)域。

3.參看資料

1.《C++ primer》
2.《The VTK User’s Guide – 11thEdition》
3. ?張曉東, 羅火靈. VTK圖形圖像開發(fā)進(jìn)階[M]. 機(jī)械工業(yè)出版社, 2015.

總結(jié)

以上是生活随笔為你收集整理的VTK修炼之道51：图形基本操作进阶_连通区域分析的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。