在VTK中连通区域分析常用到的类是：vtkPolyDataConnectivityFilter。从名字可以看出来其输入类型为vtkPolyData；了解VTK的人都知道这是图形的基本数据结构保存了图形的几何结构和拓扑结构。如果你要进行连通分析的对象是图形那自然就没问题。但如果需要对图像进行连通区域分析，那该类就不能直接用了，要先进行一步数据类型转换。根据参考文1，我们可以以一种非严谨方式将VTK数据处理算子（Filter）分为一下几类：

    1、图像处理   其输入为vtkImageData(图像基本数据类型)，输出也为vtkImageData

    2、几何提取   其输入为vtkImageData，输出为vtkPolyData

    3、图形处理   其输入为vtkPolyData，输出也为vtkPolyData，例如前面的vtkPolyDataConnectivityFilter

    4、渲染       如Mapper和Actor 来完成图形图像的显示

    为了实现vtkImageData==>vtkPolyData  我们需要使用一种几何提取处理器，参考文献一提到了一种vtkContourFilter(轮廓提取)来完成。vtkContourFilter使用较为简单，只需要设置一个参数等值线，其可以通过函数SetValue或GenerateVale()来完成设置，其具体的不同见参考文献2.

     完成数据类型转换我们就可以使用vtkPolyDataConnectivityFilter来提取连通区域。该类在提取连通区域时候有如下6中模式：

    1）提取数据集中的最大（最多点）连接区域：SetExtractionModeToLargestRegion();

    2）提取指定区域号：SetExtractionModeToSpecifiedRegions(); 

    3）提取共享指定点ids的所有区域：SetExtractionModeToPointSeededRegions(); 

    4）提取共享指定单元ID的所有区域：SetExtractionModeToCellSeededRegions(); 

    5）提取最靠近指定点的区域：SetExtractionModeToClosestPointRegion(); 

    6）提取所有区域（用于着色区域）：SetExtractionModeToAllRegions()

   区域标记号被以点的标量属性数据来保存。为了标记所有区域还必须配合函数ColorRegionOn()

  下面演示代码： 模式为SetExtractionModeToAllRegions() 获取连通区域后 遍历所有点来查看其所属于的连通区域号

//[1] --read image
    QString str="F:/lunkuo5.jpg";
    QByteArray ba=str.toLocal8Bit();
    const char *fileName=ba.data();
    vtkSmartPointer<vtkJPEGReader> reader
            =vtkSmartPointer<vtkJPEGReader>::New();
    reader->SetFileName(fileName);
    reader->Update();
    //[1]

    //[2] --灰度化
    vtkSmartPointer<vtkImageLuminance> luninanceFilter=
            vtkSmartPointer<vtkImageLuminance>::New();
    luninanceFilter->SetInputConnection(reader->GetOutputPort());
    luninanceFilter->Update();
    //[2]

    //[3] --图像二值化
    vtkSmartPointer<vtkImageThreshold> thresholdFilter=
            vtkSmartPointer<vtkImageThreshold>::New();
    thresholdFilter->SetInputConnection(luninanceFilter->GetOutputPort());
    thresholdFilter->ThresholdByUpper(100);//大于该灰度范围的区间为有效范围
    thresholdFilter->SetInValue(255);//有效范围内的输出值
    thresholdFilter->SetOutValue(0);//无效范围的输出值
    thresholdFilter->UpdateDataObject();
    thresholdFilter->UpdateInformation();
    thresholdFilter->Update();
    //[3]

    //[4] --exactor contour
    vtkSmartPointer<vtkContourFilter> contourFilter=
            vtkSmartPointer<vtkContourFilter>::New();
    contourFilter->SetInputConnection(thresholdFilter->GetOutputPort());
    contourFilter->ComputeNormalsOn();
    contourFilter->SetValue(0,100);
    contourFilter->Update();
    //[4]

    //[4] --连通区域分析
    vtkSmartPointer<vtkPolyDataConnectivityFilter> connectFilter1=
            vtkSmartPointer<vtkPolyDataConnectivityFilter>::New();
    connectFilter1->SetInputConnection(contourFilter->GetOutputPort());
    connectFilter1->SetExtractionModeToAllRegions();
    connectFilter1->ColorRegionsOn();
    connectFilter1->MarkVisitedPointIdsOn();
    connectFilter1->Update();
    int regionNum=connectFilter1->GetNumberOfExtractedRegions();
    qDebug()<<"regionNum:"<<regionNum;
    //[4]
  
    for(int i=0;i<connectFilter1->GetOutput()->GetNumberOfPoints();++i)
    {
        int regionId=connectFilter1->GetOutput()->GetPointData()->GetScalars()->GetTuple1(i);
        double *xyz=connectFilter1->GetOutput()->GetPoint(i);

    }

参考文献：

 1  第06章-图像处理及可视化（1）

 2 vtkContourFilter部分成员函数的使用说明

 3 vtkPolyDataConnectivityFilter

 4 ExtrenalContour

 5 vtk-example