图像细化算法 opencv实现 完整代码(附实验结果)
来源:互联网 发布:北斗手机导航软件 编辑:程序博客网 时间:2024/05/22 12:10
一、关于图像细化的算法可以参看下面两个PDF链接:
http://www.uel.br/pessoal/josealexandre/stuff/thinning/ftp/lam-lee-survey.pdf :总结了几乎所有92年以前的经典细化算法
http://www-prima.inrialpes.fr/perso/Tran/Draft/gateway.cfm.pdf :本文所附代码所参照的算法
二、实验结果
原图:
运行结果(分别为迭代1次、4次、8次、16次、64次的结果图):
三、opencv代码
运行环境:Windows 7,VS 2012,OpenCV 2.3
#include <stdlib.h>#include <stdio.h>#include <math.h>#include <fstream>#include <string>#include <iostream>#include <opencv/cv.h>#include <opencv/highgui.h> using namespace std; void cvThin( IplImage*, IplImage*, int iterations);//使用举例int main(int argc, char* argv[]){ IplImage *pSrc = NULL,*pDst = NULL,*pTmp = NULL;//传入一个灰度图像 pSrc = cvLoadImage("1.bmp",CV_LOAD_IMAGE_GRAYSCALE); if(!pSrc) { return 0; } pTmp = cvCloneImage(pSrc); pDst = cvCreateImage(cvGetSize(pSrc),pSrc->depth,pSrc->nChannels); cvZero(pDst); cvThreshold(pSrc,pTmp,68,1,CV_THRESH_BINARY_INV);//做二值处理,将图像转换成0,1格式 //cvSaveImage("pTmp.bmp",pTmp); cvThin(pTmp,pDst,8);//细化,通过修改iterations参数进一步细化 cvNamedWindow("src",1); cvNamedWindow("dst",1); cvShowImage("src",pSrc); //将二值图像转换成灰度,以便显示 int i = 0,j = 0; CvSize size = cvGetSize(pDst); for(i=0; i<size.height; i++) { for(j=0; j<size.width; j++) { if(CV_IMAGE_ELEM(pDst,uchar,i,j)==1) { CV_IMAGE_ELEM(pDst,uchar,i,j) = 0; } else { CV_IMAGE_ELEM(pDst,uchar,i,j) = 255; } } } cvSaveImage("dst_8.bmp",pDst); cvShowImage("dst",pDst); cvWaitKey(0); cvReleaseImage(&pSrc); cvReleaseImage(&pDst); cvReleaseImage(&pTmp); cvDestroyWindow("src"); cvDestroyWindow("dst"); return 0;}void cvThin( IplImage* src, IplImage* dst, int iterations){ //此时的src是一个二值化的图片 CvSize size = cvGetSize(src); cvCopy(src, dst); int n = 0,i = 0,j = 0; for(n=0; n<iterations; n++)//开始进行迭代 { IplImage* t_image = cvCloneImage(dst); for(i=0; i<size.height; i++) { for(j=0; j<size.width; j++) { if(CV_IMAGE_ELEM(t_image,byte,i,j)==1) { int ap=0; int p2 = (i==0)?0:CV_IMAGE_ELEM(t_image,byte, i-1, j); int p3 = (i==0 || j==size.width-1)?0:CV_IMAGE_ELEM(t_image,byte, i-1, j+1); if (p2==0 && p3==1) { ap++; } int p4 = (j==size.width-1)?0:CV_IMAGE_ELEM(t_image,byte,i,j+1); if(p3==0 && p4==1) { ap++; } int p5 = (i==size.height-1 || j==size.width-1)?0:CV_IMAGE_ELEM(t_image,byte,i+1,j+1); if(p4==0 && p5==1) { ap++; } int p6 = (i==size.height-1)?0:CV_IMAGE_ELEM(t_image,byte,i+1,j); if(p5==0 && p6==1) { ap++; } int p7 = (i==size.height-1 || j==0)?0:CV_IMAGE_ELEM(t_image,byte,i+1,j-1); if(p6==0 && p7==1) { ap++; } int p8 = (j==0)?0:CV_IMAGE_ELEM(t_image,byte,i,j-1); if(p7==0 && p8==1) { ap++; } int p9 = (i==0 || j==0)?0:CV_IMAGE_ELEM(t_image,byte,i-1,j-1); if(p8==0 && p9==1) { ap++; } if(p9==0 && p2==1) { ap++; } if((p2+p3+p4+p5+p6+p7+p8+p9)>1 && (p2+p3+p4+p5+p6+p7+p8+p9)<7) { if(ap==1) { if(!(p2 && p4 && p6)) { if(!(p4 && p6 && p8)) { CV_IMAGE_ELEM(dst,byte,i,j)=0;//设置目标图像中像素值为0的点 } } } } } } } cvReleaseImage(&t_image); t_image = cvCloneImage(dst); for(i=0; i<size.height; i++) { for(int j=0; j<size.width; j++) { if(CV_IMAGE_ELEM(t_image,byte,i,j)==1) { int ap=0; int p2 = (i==0)?0:CV_IMAGE_ELEM(t_image,byte, i-1, j); int p3 = (i==0 || j==size.width-1)?0:CV_IMAGE_ELEM(t_image,byte, i-1, j+1); if (p2==0 && p3==1) { ap++; } int p4 = (j==size.width-1)?0:CV_IMAGE_ELEM(t_image,byte,i,j+1); if(p3==0 && p4==1) { ap++; } int p5 = (i==size.height-1 || j==size.width-1)?0:CV_IMAGE_ELEM(t_image,byte,i+1,j+1); if(p4==0 && p5==1) { ap++; } int p6 = (i==size.height-1)?0:CV_IMAGE_ELEM(t_image,byte,i+1,j); if(p5==0 && p6==1) { ap++; } int p7 = (i==size.height-1 || j==0)?0:CV_IMAGE_ELEM(t_image,byte,i+1,j-1); if(p6==0 && p7==1) { ap++; } int p8 = (j==0)?0:CV_IMAGE_ELEM(t_image,byte,i,j-1); if(p7==0 && p8==1) { ap++; } int p9 = (i==0 || j==0)?0:CV_IMAGE_ELEM(t_image,byte,i-1,j-1); if(p8==0 && p9==1) { ap++; } if(p9==0 && p2==1) { ap++; } if((p2+p3+p4+p5+p6+p7+p8+p9)>1 && (p2+p3+p4+p5+p6+p7+p8+p9)<7) { if(ap==1) { if(p2*p4*p8==0) { if(p2*p6*p8==0) { CV_IMAGE_ELEM(dst, byte,i,j)=0; } } } } } } } cvReleaseImage(&t_image); }}
四、参考链接
[1]http://tech.groups.yahoo.com/group/OpenCV/message/70260
[2]http://blog.csdn.net/byxdaz/article/details/5642669
- 图像细化算法 opencv实现 完整代码(附实验结果)
- 图像细化算法 opencv实现 完整代码(附实验结果)
- 图像拼接实现镜面效果 opencv完整代码 实现(附实验结果)
- 图像融合实现CT与PET图像的融合 OpenCV实现(附实验结果)
- 并行图像细化算法和C代码实现 (转)
- opencv实现图像细化效果
- opencv实现二值图像细化的算法
- opencv实现二值图像细化的算法
- OpenCV实现二值图像细化的算法
- 图像细化matlab代码实现
- 并行图像细化算法和C代码实现
- opencv实现二值图像的细化
- Opencv 图像增强算法 图像检测结果及代码
- 【opencv】图像细化
- opencv-图像细化
- Opencv----图像细化
- Matlab的图像细化算法(详细实现)
- 细化处理(二值图像)Hilditch细化算法
- SimpleDateFormat格式化日期
- jstl fn函数
- 冒泡排序 --跳出循环的break的flag值放错位置
- CCS3.3编译链接库文件与头文件
- c++中static的全部用法
- 图像细化算法 opencv实现 完整代码(附实验结果)
- 双单击
- java List 排序 Collections.sort() 对 List 排序
- VS2010 安装配置 vtk
- iOS核心系统编程最佳实践:绘图
- BitmapFactory decodeResource and decodeStream 's problem
- iOS核心系统编程最佳实践:动画
- JAVA_OPTS参数说明与配置
- Android杂谈--打开文件的Intent及使用