-
区域生长法
2018-06-16 21:07:56区域生长的基本思想是将具有相似性质的像素集合起来构成区域。具体先对每个需要分割的区域找一个种子像素作为生长的起点,然后将种子像素周围邻域中与种子像素具有相同或相似性质的像素(根据某种事先确定的生长或...区域生长的基本思想是将具有相似性质的像素集合起来构成区域。具体先对每个需要分割的区域找一个种子像素作为生长的起点,然后将种子像素周围邻域中与种子像素具有相同或相似性质的像素(根据某种事先确定的生长或相似准则来判定)合并到种子像素所在的区域中。将这些新像素当做新的种子像素继续进行上面的过程,直到再没有满足条件的像素可被包括进来,这样,一个区域就长成了。
区域生长是指从某个像素出发,按照一定的准则,逐步加入邻近像素,当满足一定的条件时,区域生长终止。区域生长的好坏决定于1.初始点(种子点)的选取。2.生长准则。3.终止条件。区域生长是从某个或者某些像素点出发,最后得到整个区域,进而实现目标的提取。
简单来说下三个法则,对出需要分割的图像:1、选取图像中的一点为种子点(种子点的选取需要具体情况具体分析)。2、在种子点处进行8邻域或4邻域扩展,判定准则是:如果考虑的像素与种子像素灰度值差的绝对值小于某个门限T,则将该像素包括进种子像素所在的区域。3、当不再有像素满足加入这个区域的准则时,区域生长停止。
区域生长实现的步骤如下:
1. 对图像顺序扫描!找到第1个还没有归属的像素, 设该像素为(x0, y0);
2. 以(x0, y0)为中心, 考虑(x0, y0)的8邻域像素(x, y),如果(x,, y)满足生长准则, 将(x, y)与(x0, y0)合并(在同一区域内), 同时将(x, y)压入堆栈;
3. 从堆栈中取出一个像素, 把它当作(x0, y0)返回到步骤2;
4. 当堆栈为空时!返回到步骤1;
5. 重复步骤1 - 4直到图像中的每个点都有归属时。生长结束。
-
4.1图像分割之区域生长法
2018-08-04 16:37:39其中最基础的是区域生长法。 区域生长法 区域生长是根据事先定义的准则将像素或者子区域聚合成更大区域的过程。其基本思想是从一组生长点开始(生长点可以是单个像素,也可以是某个小区域),将与该...我们将讨论以区域为基础的图像分割处理技术。传统的区域分割方法有区域生长和区域分裂与合并, 其中最基础的是区域生长法。
区域生长法
区域生长是根据事先定义的准则将像素或者子区域聚合成更大区域的过程。其基本思想是从一组生长点开始(生长点可以是单个像素,也可以是某个小区域),将与该生长点性质相似的相邻像素或者区域与生长点合并,形成新的生长点,重复此过程直到不能生长为止。生长点和相似区域的相似性判断依据可以是灰度值、纹理、颜色等图像信息。所以区域生长算法关键有三个:
1、选择合适的生长点
2、确定相似性准则即生长准则
3、确定生长停止条件
下面给出一个区域生长的实例:图(a)为原始图像,数字表示像素的灰度。以灰度值为8的像素为初始的生长点,记为f(i,j)。在8邻域内,生长准则是待测点灰度值与生长点灰度值相差为1或0.那么图(b)是第一次区域生长后,f(i-1,j)、f(i,j-1)、f(i,j+1)和生长点灰度值相差都是1,因而被合并。图©是第二次生长后,f(i+1,j)被合并。图(d)为第三次生长后,f(i+1,j-1)、f(i+2,j)被合并,至此,已经不存在满足生长准则的像素点,生长停止。
示例演示
我们实现一个功能,鼠标选择一个初始生长点,生长准则是待测点灰度值与初始生长点灰度值差的绝对值不大于8。示例代码中regionGrow函数参考了《数字图像处理与机器视觉》第九章第五节,regionGrowFast函数是优化过的实现,代码如下:
cv::Mat MainWindow::regionGrowFast(const cv::Mat &src, const cv::Point2i seed, int throld) { //convert src to gray for getting gray value of every pixel cv::Mat gray; cv::cvtColor(src,gray, cv::COLOR_RGB2GRAY); // set every pixel to black cv::Mat result = cv::Mat::zeros(src.size(), CV_8UC1); if((seed.x < 0) || (seed.y < 0)) return result; result.at<uchar>(seed.y, seed.x) = 255; //gray value of seed int seed_gray = gray.at<uchar>(seed.y, seed.x); //grow direction sequenc int grow_direction[8][2] = {{-1,-1}, {0,-1}, {1,-1}, {1,0}, {1,1}, {0,1}, {-1,1}, {-1,0}}; //seeds collection std::vector<cv::Point2i> seeds; seeds.push_back(seed); //start growing while(! seeds.empty()){ //get a seed cv::Point2i current_seed = seeds.back(); seeds.pop_back(); for(int i = 0; i < 8; ++i){ cv::Point2i neighbor_seed(current_seed.x + grow_direction[i][0], current_seed.y + grow_direction[i][1]); //check wether in image if(neighbor_seed.x < 0 || neighbor_seed.y < 0 || neighbor_seed.x > (gray.cols-1) || (neighbor_seed.y > gray.rows -1)) continue; int value = gray.at<uchar>(neighbor_seed.y, neighbor_seed.x); if((result.at<uchar>(neighbor_seed.y, neighbor_seed.x) == 0) && (abs(value - seed_gray) <= throld)){ result.at<uchar>(neighbor_seed.y, neighbor_seed.x) = 255; seeds.push_back(neighbor_seed); } } } return result; }
运行结果:
-
区域生长法MATLAB版本程序代码
2021-01-14 11:00:31区域生长法MATLAB版本程序代码,需要把图像名字改成程序里一致的 区域生长法MATLAB版本程序代码,需要把图像名字改成程序里一致的 -
区域生长法分割
2020-06-24 10:48:10区域分割区域生长法 区域生长法 区域生长法的基本思想就是将具有相似性的像素集合起来形成一个区域。 具体做法就是,首先选择一个种子点,通过比较种子点邻域的相似性,将邻域中满足相似性准则的像素归入种子点所在...区域分割
区域生长法
区域生长法的基本思想就是将具有相似性的像素集合起来形成一个区域。
具体做法就是,首先选择一个种子点,通过比较种子点邻域的相似性,将邻域中满足相似性准则的像素归入种子点所在的区域,然后将这新的像素点定为种子点,重复上述过程,直到没有满足相似性准则的新的邻域像素点产生为止。通过区域生长,一个区域就形成了。
有三个关键性的问题:
1.种子点的选择问题,如何选择代表所需区域的种子点?
借助具体问题的特点来进行选择。
2.生长过程中停止生长的条件是什么?
一般是再没有满足生长准则的像素产生为止。
3.如何判断相邻两个像素是否是相似的呢?
①基于区域灰度差的方法。
1)对所有的像素进行扫描,找出尚没有归属的像素点。
2)对这像素点的所有邻域像素进行灰度差,找出灰度差满足某一阈值的邻域像素,将他们合并。
3)以新合并的像素为新的中心,进行步骤2)。检查新像素的邻域,直到区域不能继续扩张为止。
4)返回步骤1),继续扫描,直到所有像素都有归属为止。
改进–>该方法对于生长起点的选择具有很大的依赖性。将灰度差阈值设为0,这样所有灰度值相同的像素点就形成了一个区域。然后比较相邻区域之间的平均灰度差,将平均灰度差小于阈值的区域进行合并。问题:对于变化缓慢的相邻区域,会产生错分割。
改进–>在进行区域生长时,不使用新像素点的灰度值与其邻域进行比较,而是使用新像素点所在的区域的平均灰度值与其邻域的像素点的灰度值进行比较,选择满足阈值的点进行合并。
②以灰度分布的相似性为准则。
1)将像素分成互不重叠的小区域。
2)比较相邻区域的累计灰度直方图,即比较灰度分布的相似性。参考https://blog.csdn.net/combfish/article/details/45056239
其实累积直方图就是对直方图分布进行归一化处理。将小于阈值的两个相邻区域进行合并。
3)设定终止准则,通过反复步骤2),依次将各个区域进行合并。区域生长法的优点是:计算简单,比较适合分割均匀的小结构,往往与其他分割方法联合使用。
缺点是:对初始种子点的依赖性,而且对噪声比较敏感,出现分割空洞或者过分割。 -
区域生长法图像分割
2021-01-18 17:13:30传统的区域分割方法有区域生长和区域分裂与合并,其中最基础的是区域生长法。 区域生长法 区域生长是根据事先定义的准则将像素或者子区域聚合成更大区域的过程。其基本思想是从一组生长点开始(生长点可以是单个像素...我们将讨论以区域为基础的图像分割处理技术。传统的区域分割方法有区域生长和区域分裂与合并,其中最基础的是区域生长法。
区域生长法
区域生长是根据事先定义的准则将像素或者子区域聚合成更大区域的过程。其基本思想是从一组生长点开始(生长点可以是单个像素,也可以是某个小区域),将与该生长点性质相似的相邻像素或者区域与生长点合并,形成新的生长点,重复此过程直到不能生长为止。生长点和相似区域的相似性判断依据可以是灰度值、纹理、颜色等图像信息。所以区域生长算法关键有三个:
1、选择合适的生长点
2、确定相似性准则即生长准则
3、确定生长停止条件下面给出一个区域生长的实例:图(a)为原始图像,数字表示像素的灰度。以灰度值为8的像素为初始的生长点,记为f(i,j)。在8邻域内,生长准则是待测点灰度值与生长点灰度值相差为1或0.那么图(b)是第一次区域生长后,f(i-1,j)、f(i,j-1)、f(i,j+1)和生长点灰度值相差都是1,因而被合并。图©是第二次生长后,f(i+1,j)被合并。图(d)为第三次生长后,f(i+1,j-1)、f(i+2,j)被合并,至此,已经不存在满足生长准则的像素点,生长停止。
示例代码:#include <opencv2/opencv.hpp> //头文件 #include <opencv2/highgui.hpp> #include <iostream> #include <map> using namespace cv; using namespace std; /*************************************************************************************** Function: 区域生长算法 Input: src 待处理原图像 pt 初始生长点 th 生长的阈值条件 Output: 肺实质的所在的区域 实质区是白色,其他区域是黑色 Description: 生长结果区域标记为白色(255),背景色为黑色(0) Return: NULL Others: NULL ***************************************************************************************/ void RegionGrow(cv::Mat& src, cv::Mat& matDst, cv::Point2i pt, int th = 40) { cv::Point2i ptGrowing; //待生长点位置 int nGrowLable = 0; //标记是否生长过 int nSrcValue = 0; //生长起点灰度值 int nCurValue = 0; //当前生长点灰度值 matDst = cv::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 } }; std::vector<cv::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 < 8; ++i) { ptGrowing.x = pt.x + DIR[i][0]; ptGrowing.y = pt.y + DIR[i][1]; //检查是否是边缘点 if (ptGrowing.x < 0 || ptGrowing.y < 0 || ptGrowing.x >(src.cols - 1) || (ptGrowing.y > src.rows - 1)) continue; nGrowLable = matDst.at<uchar>(ptGrowing.y, ptGrowing.x); //当前待生长点的灰度值 if (nGrowLable == 0) //如果标记点还没有被生长 { nCurValue = src.at<uchar>(ptGrowing.y, ptGrowing.x); if (abs(nSrcValue - nCurValue) < th) //在阈值范围内则生长 { matDst.at<uchar>(ptGrowing.y, ptGrowing.x) = 255; //标记为白色 vcGrowPt.push_back(ptGrowing); //将下一个生长点压入栈中 } } } } } void on_MouseHandle(int event, int x, int y, int flags, void* param) { cv::Mat& src = *(cv::Mat*) param; cv::Mat src_gray, dst; if (src.channels() > 1) cv::cvtColor(src, src_gray, CV_RGB2GRAY); else src_gray = src.clone(); cv::Point2i pt; switch (event) { //左键按下 case cv::EVENT_LBUTTONDOWN: { //x:列 y:行 pt = cv::Point2i(x, y); std::cout << "(x,y)=" << "(" << x << "," << y << ")" << std::endl; } break; //左键放开 char str[16]; case cv::EVENT_LBUTTONUP: { //cv::circle(src, cv::Point2i(x, y), 1, cv::Scalar(0, 0, 255), -1, CV_AA); //sprintf_s(str, "(%d,%d)", x, y); //cv::putText(src, str, cv::Point2i(x, y), 3, 1, cv::Scalar(150, 200,0), 2, 8); pt = cv::Point2i(x, y); RegionGrow(src_gray, dst, pt); //区域生长 cv::bitwise_and(src_gray, dst, dst); //与运算 //imshow("src", src); imshow("dst", dst); } break; } } int main() { cv::Mat src = cv::imread("toux.jpg", 0); if (src.empty()) { return -1; } /* int d: 表示在过滤过程中每个像素邻域的直径范围。如果这个值是非正数,则函数会从第五个参数sigmaSpace计算该值。 double sigmaColor: 颜色空间过滤器的sigma值,这个参数的值月大,表明该像素邻域内有越宽广的颜色会被混合到一起,产生较大的半相等颜色区域。 (这个参数可以理解为值域核w_r的\sigma_r) double sigmaSpace: 坐标空间中滤波器的sigma值,如果该值较大,则意味着越远的像素将相互影响,从而使更大的区域中足够相似的颜色获取相同的颜色。 当d>0时,d指定了邻域大小且与sigmaSpace无关,否则d正比于sigmaSpace. (这个参数可以理解为空间域核w_d的\sigma_d) int borderType=BORDER_DEFAULT: 用于推断图像外部像素的某种边界模式,有默认值BORDER_DEFAULT. */ int d = 7; double sigmaColor = 140; double sigmaSpace = 140; //cv::imshow("原图", src); cv::Mat gauss_src; cv::bilateralFilter(src, gauss_src, d, sigmaColor, sigmaSpace, cv::BORDER_DEFAULT); cv::namedWindow("双边模糊", CV_WND_PROP_AUTOSIZE);//定义一个img窗口 cv::namedWindow("dst", CV_WND_PROP_AUTOSIZE);//定义一个dst窗口 imshow("双边模糊", gauss_src); cv::setMouseCallback("双边模糊", on_MouseHandle, (void*)&gauss_src);//调用回调函数 cv::waitKey(0); }
结果:
-
区域生长法实现matlab代码
2017-12-03 19:40:41利用区域生长法来实现图像分割,本人已经用肝部CT图像来测试过,非常好用。 -
区域生长法MATLAB
2020-07-08 20:01:09区域生长算法的设计主要由以下三点:生长种子点的确定,区域生长的条件,区域生长停止的条件。 代码采用交互的方式,用鼠标选取合适的种子点,搜索种子点8邻域内的像素。 -
区域生长法 讨论
2020-01-14 06:44:29下面是《数字图像处理与机器视觉VisualC++与Matlab实现》(张铮等编著)配套代码中的区域生长法,大家能不能发现其中的缺陷,此外,其中只能选择一个种子点,你能修改为可选择多个种子点吗?并且加一个算法实现种子... -
区域生长法VS实现。C++
2016-04-11 20:25:40区域生长法VS实现。C++ -
基于c++的区域生长法算法
2018-04-08 18:26:42用visual studio2015编译实现的区域生长法去实现图像分割。注意opencv的环境配置。 -
基于改进的区域生长法的气管与支气管分割
2021-02-20 15:36:43基于改进的区域生长法的气管与支气管分割 -
多点区域生长法 matlab实现
2018-04-28 15:38:40关于区域生长法网上已经有很多算法可以实现了,但是关于多点区域生长法却没什么资料,多点可以很大程度上增加生长速度,经过琢磨才发现在单点区域生长的算法上很容易就能实现多点生长,因此写出来供大家参考。... -
区域生长法的图像分割matlab程序
2015-05-27 13:09:53本程序主要功能为实现基于区域生长法的图像分割,选取种子点,并将灰度差值小于阈值的像素点进行合并生长。 -
matlab图像处理--区域生长法
2018-05-10 16:48:50区域生长法:区域生长是按照事先定义的生长准则将一个像素或者子区域逐步聚合成一个完整独立的连通区域过程。对于图像感兴趣目标区域R,z为区域R上事先发现的种子点,按照规定的生长准则逐步将与种子点z一定邻域内... -
基于区域生长法的图像分割matlab程序
2013-10-19 19:28:31本程序主要功能为实现基于区域生长法的图像分割,选取种子点,并将灰度差值小于阈值的像素点进行合并生长。 -
区域生长图像分割算法_基于综合颜色特征和区域生长法对病斑图像进行分割
2020-12-17 00:47:52点击蓝色字免费订阅,每天收到这样的好资讯本文提出了一种利用综合颜色特征(ACCF)和区域生长法对叶片病斑进行分割的方法,植物表型资讯介绍如下:在田间条件下获取叶片病斑图像主要面临两大挑战,一是背景杂乱,二是... -
区域生长法的C++实现(opencv)
2021-03-26 15:43:10区域生长法的C++实现(opencv) 区域生长法是一种图像分割的常用算法,能够将目标区域精确得分割出来。我的程序是将图像最中央的白色圆形区域提取出来,因此第一个种子点选择在了图像中央。 待处理图片: 提取结果... -
基于彩色图像区域生长法的变电站红外图像分割
2020-06-26 17:45:57本文介绍了变电站红外图像的特点,并通过对现有分割方法的分析,提出了适合实际图像特点的彩色图像区域生长法。实验结果表明,该方法能较好的实现红外图像的分割,有利于后期的图像识别和特征数据的提取,适合于在变电站... -
图像分割——区域生长法
2016-04-08 22:38:20原文地址:[转载]图像分割——区域生长法作者:小腹黑zju 原文出处:http://blog.csdn.net/cay22/article/details/5666091 区域生长方法是根据同一物体区域内象素的相似性质来聚集象素点的方法,从... -
基于C++的区域生长法实现
2018-04-08 18:16:52图像分割是一种重要的数字图像处理技术,生长区域法是选择一个种子点(图像中的某一像素点),对二维图像来言,先遍历... 区域生长法的算法步骤: 1.选择一个种子点,将种子点push到stack。 2.从stack中取出后... -
论文研究-基于形态学的遮挡目标提取的区域生长法研究.pdf
2019-07-22 19:43:22采用改进的区域生长法进行红外热图中遮挡目标的提取。首先利用阈值化方法提取出一个包含全部目标信息的矩形子区域,后续处理仅针对该矩形区域进行处理。然后证明了数学形态学的膨胀梯度定理(dilation gradient ... -
二维码 区域上色 区域生长法
2014-04-01 17:30:15参照csdn上一位大牛对二维码区域上色的思路,自己写的使用区域增长法二维码上色。希望有人在这基础上完成二维码圆角化。 -
二值图像区域生长法+区域面积计算
2016-08-22 16:34:13这个帖子,区域生长法是借鉴现有算法,区域面积计算是根据自己的需求自行加上的,暂且归为原创,如有雷同,请联系我。 这里,前半部分,区域生长法,可以参见此博客:... -
转载]图像分割——区域生长法
2016-05-14 16:51:17转载]图像分割——区域生长法 (2012-10-12 09:02:55) 转载▼ 标签: 图像分割 区域生长法 杂谈 分类: 学术 原文出处:http://blog.csdn.net/cay22/article/details/5666091 ... -
Matlab图像分割(三)区域生长法
2020-04-21 16:02:37Matlab图像分割——区域生长法分割代码 [fileName,pathName] = uigetfile('*.*','Please select an image');%文件筐,选择文件 if(fileName) fileName = strcat(pathName,fileName); fileName = lower(fileName);%... -
区域生长算法 python代码_OpenCV:区域生长法实现
2020-12-15 18:20:04示例程序:OpenCV没有自带区域生长函数;为了说明区域是如何生长的,#include #include #include #include using namespace std;using namespace cv;void RegionGrowing(Mat srcImg,Mat& dstImg,Point pt,int th...