OpenCV笔记12 距离变换

一、距离图像 

距离变换于被广泛应用于图像分析、计算机视觉、模式识别等领域,人们利用它来实现目标细化、骨架提取、形状插值及匹配、粘连物体的分离等。距离变换的图像是二值图像针。二值图像可以做事两个部分组成的前景目标和背景,前景目标的像素值为1,背景的像素值为0;距离变换的结果是一幅灰度级图像,即距离图像,图像中每个像素的灰度值为该像素与距其最近的背景像素间的距离。
距离变换算法主要采用两类距离测度:非欧式距离和欧式距离。

1. 欧式距离:其物理意义表示距离小于等于某一个值得像素中心(x,y)且半径为R的圆平面,r

如上图所示,p1h和p2 ,dist=1.4.14

2. 棋盘距离:是指数字栅格上像素按照对角线方向移动形成的距离

如上图所示,p1h和p2 ,dist=2

3. 曼哈顿距离:从起点到终点所需要的最小步数,只允许横向移动和纵向移动

如上图所示,p1h和p2 ,dist=1

 例如:如下1表示前景目标,0表示背景,则曼哈顿(城市街区)距离:

 算子解释

 void cv::distanceTransform(
  InputArray  src,  // cv_8U的单通道的图像
  OutputArray  dst, // 输出图像, cv_8U或者CV_32F
  OutputArray  labels,// 二维的标签数组(离散Voronoi图),与输入图像具有相同的尺寸,数据类型为CV_32S的单通道数据
  int  distanceType, // 计算距离的方法  DIST_USER   -1 ,表示用户自定义距离
                                          CV_DIST_L1,  1, 曼哈顿距离(城市街区距离)  距离的计算是精确的
                                          CV_DIST_L2,  2,  欧式距离 计算有某些相对误差 (5×5 mask 给出更精确的结果) 
                                          CV_DIST_C    3,  棋盘距离  距离的计算是精确的

  int  maskSize,// 掩码大小  :DIST_MASK_3(3×3)和DIST_MASK_5(5×5)
  int  labelType = DIST_LABEL_CCOMP  //labelType == DIST_LABEL_CCOMP时,该函数会自动在输入图像中找到0像素的连通分量,并用不同的标签标记它们。当labelType == DIST_LABEL_CCOMP时,该函数扫描输入图像并用不同的标签标记所有0像素。
  )

 

距离变换算法步骤:

1.将图像进行二值化,子图像值为0,背景为255;

2.利用Mask 1从左向右,从上到下扫描,p点是当前像素点,q点是Mask 1中AL邻域中的点,D为距离计算,包括棋盘距离、城市距离和欧式距离。Fpp为p点的像素值,计算

F(p) = min( F(p),  F(q)+D(p,q) ), 其中,q属于AL.

3.再利用Mask 2从右向左,从下向上扫描

Fpp = minF(pF(p, Fqq+Dp,qp,q ), 其中,q属于BR

4.Fpp 则为距离变换后的图像。

代码详情请看:距离变换 - 芒果的博客 - 芒果的个人博客

  • 3
    点赞
  • 18
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
下面是一个基于opencv的分水岭距离变换的代码实现: ```python import cv2 import numpy as np # 读取像 img = cv2.imread('image.png') gray = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY) # 进行阈值化处理 ret, thresh = cv2.threshold(gray, 0, 255, cv2.THRESH_BINARY_INV+cv2.THRESH_OTSU) # 计算距离变换像 dist_transform = cv2.distanceTransform(thresh, cv2.DIST_L2, 5) # 对距离变换像进行阈值化 ret, sure_fg = cv2.threshold(dist_transform, 0.7*dist_transform.max(), 255, 0) # 对距离变换像求反 sure_fg = np.uint8(sure_fg) unknown = cv2.subtract(thresh, sure_fg) # 对距离变换像进行膨胀操作 kernel = np.ones((3, 3), np.uint8) sure_fg = cv2.dilate(sure_fg, kernel, iterations=3) # 计算标记像 ret, markers = cv2.connectedComponents(sure_fg) # 将标记像加1,未知区域设为0 markers = markers+1 markers[unknown==255] = 0 # 使用分水岭算法进行像分割 markers = cv2.watershed(img, markers) img[markers==-1] = [255,0,0] # 显示结果 cv2.imshow('image', img) cv2.waitKey(0) cv2.destroyAllWindows() ``` 这段代码首先读取一个像,然后将像转换为灰度像,并进行阈值化处理。接着,使用`distanceTransform()`函数计算距离变换像,并对其进行阈值化,得到前景区域。然后,对距离变换像进行膨胀操作,得到局部极大值区域作为分割线。接下来,使用`connectedComponents()`函数计算标记像,然后将标记像加1,未知区域设为0,最后使用`watershed()`函数进行分割。最后,将分割线标记的部分用蓝色进行标记,并显示结果。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值