医学影像处理：区域生长算法详解

/*************************************************************************************** 
Function:  区域生长算法 
Input:     src 待处理原图像 pt 初始生长点 th 生长的阈值条件 
Output:    肺实质的所在的区域 实质区是白色，其他区域是黑色 
Description: 生长结果区域标记为白色(255),背景色为黑色(0) 
Return:    Mat 
Others:    NULL 
***************************************************************************************/  
Mat RegionGrow(Mat src, Point2i pt, int th)  
{  
Point2i ptGrowing;                      //待生长点位置  
int nGrowLable = 0;                             //标记是否生长过  
int nSrcValue = 0;                              //生长起点灰度值  
int nCurValue = 0;                              //当前生长点灰度值  
Mat matDst = Mat::zeros(src.size(), CV_8UC1);   //创建一个空白区域，填充为黑色  
//生长方向顺序数据  
int DIR[8][2] = {{-1,-1}, {0,-1}, {1,-1}, {1,0}, {1,1}, {0,1}, {-1,1}, {-1,0}};    
Vector<Point2i> vcGrowPt;                     //生长点栈  
vcGrowPt.push_back(pt);                         //将生长点压入栈中  
matDst.at<uchar>(pt.y, pt.x) = 255;               //标记生长点  
nSrcValue = src.at<uchar>(pt.y, pt.x);            //记录生长点的灰度值  

while (!vcGrowPt.empty())                       //生长栈不为空则生长  
{  
pt = vcGrowPt.back();                       //取出一个生长点  
vcGrowPt.pop_back();                          

//分别对八个方向上的点进行生长  
for (int i = 0; i<9; ++i)  
{