网站域名为个人的公司能备案,宜兴建设公司网站,深圳网站设计公司龙岗,郑州市建设安全管理协会网站Hi#xff0c;大家好#xff0c;我是半亩花海。本实验基于 OpenCV 实现了二值图像数字矩阵的距离变换算法。首先生成一个 480x480 的黑色背景图像#xff08;定义黑色为0#xff0c;白色为1#xff09;#xff0c;在其中随机选择了三个白色像素点作为距离变换的原点… Hi大家好我是半亩花海。本实验基于 OpenCV 实现了二值图像数字矩阵的距离变换算法。首先生成一个 480x480 的黑色背景图像定义黑色为0白色为1在其中随机选择了三个白色像素点作为距离变换的原点利用 OpenCV 中 distanceTransform 等相关函数计算并输出这些原点到其他像素点的欧氏距离、D4 距离和 D8 距离及其相应的距离矩阵并将距离变换结果可视化。 文章目录 一、导入必要库二、初始化输入图像和变换结果图像三、根据二值图计算并输出距离矩阵四、将距离矩阵转换为可视化图片 一、导入必要库
导入必要的库cv2、numpy、matplotlib.pyplot为后续的图像处理任务做准备。
#!/usr/bin/env python
# coding: utf-8
import cv2
import numpy as np
import matplotlib.pyplot as plt# 用来设置字体样式(黑体)以正常显示中文标签
plt.rcParams[font.sans-serif] [SimHei]
plt.rcParams[axes.unicode_minus] False二、初始化输入图像和变换结果图像
自定义 480x480 二值图随机生成图案初始值都是 0黑色并显示原始图像。 图像数字化 通过传感器获得的图像是平面坐标 ( x , y ) (x,y) (x,y) 的连续函数 f ( x , y ) f(x,y) f(x,y)它的值图像对应位置的亮度。为了能够让计算机来处理需要对图像进行采样并且对亮度值进行量化。 采样 对连续函数 f ( x , y ) f(x,y) f(x,y) 进行采样就是分别对 x x x 轴和 y y y 轴按照固定间隔取值得到平面坐标上的 M × N M×N M×N 个点将其函数值作为元素生成 M M M 行 N N N 列的矩阵。 量化亮度值。 将 f ( x , y ) f(x,y) f(x,y) 的值转化为等价的整数值的过程称为量化量化的级别越高图像越细致。通常将亮度值表示为 0255 之间的整数。
# 创建一个大小为 480x480 的全黑图像确保数据类型为 8 位无符号整数
mat np.zeros((480, 480), dtypenp.uint8)# 给输入图像指定三个像素点设置为白色像素作为距离变换原点区域块
mat[100, 200] 1
mat[200, 100] 1
mat[300, 300] 1# 显示原始图像和距离变换结果
plt.figure(figsize(15, 5))
plt.scatter([100, 200, 300], [200, 100, 300], colorwhite, markero) # 添加三个白点表示出三个像素点的坐标
plt.imshow(mat, cmapgray)
plt.title(原始图像, fontsize16)# 将输入图像中 1白色和 0黑色调换使得原点距离为 0
mat 1 - mat三、根据二值图计算并输出距离矩阵
距离 距离是描述图像两点像素之间的远近关系的度量常见的度量距离有欧式距离Euchildean distance、城市街区距离City block distance、棋盘距离Chessboard distance。以下以两坐标点 a ( i , j ) a (i, j) a(i,j) 和 b ( k , l ) b (k, l) b(k,l) 的距离为例来说明各种距离的定义方式。
1欧式距离 D e D_e De 欧式距离的定义源于经典的几何学与我们数学中所学的简单几何的两点之间的距离一致为两个像素点坐标值的平方根。欧式距离的优点在于其定义非常地直观是显而易见的但缺点在于平方根的计算是非常耗时的。 D e ( a , b ) ( ( i − k ) 2 ) ( j − l ) 2 D_e(a, b)\sqrt{\left((i-k)^2\right)(j-l)^2} De(a,b)((i−k)2)(j−l)2
2城市街区距离 D 4 D_4 D4 距离描述的是只允许像素坐标系平面中横向和纵向的移动距离4表示在这种定义下像素点是 4 邻接的即每个点只与它的上、下、左、右相邻的 4 个点之间的距离为 1。 D 4 ( a , b ) ∣ i − k ∣ ∣ j − l ∣ D_4(a, b)|i-k||j-l| D4(a,b)∣i−k∣∣j−l∣
3棋盘距离 D 8 D_8 D8 如果允许在图像坐标系中像素点的对角线方向的移动就可以得到棋盘距离8 表示在这种定义下像素点是 8 邻接的即每个点只与它的上、下、左、右、四个对角线方向相邻的 8 个点之间的距离为 1。 D 8 ( a , b ) max { ∣ i − k ∣ , ∣ j − l ∣ } D_8(a, b)\max \{|i-k|,|j-l|\} D8(a,b)max{∣i−k∣,∣j−l∣}
距离变换
距离变换也叫作距离函数或者斜切算法。它是距离概念的一个应用图像处理的一些算法以距离变换为基础。距离变换描述的是图像中像素点与某个区域块的距离区域块中的像素点值为 0临近区域块的像素点有较小的值离它越远值越大。
以二值图像为例其中区域块内部的像素值为 1其他像素值为 0。距离变换给出每个像素点到最近的区域块边界的距离区域块内部的距离变换结果为0。输入图像如图 1 所示 D 4 D_4 D4 距离的距离变换结果如图 2 所示。 下面来讨论距离变换算法其核心是利用两个小的局部掩膜遍历图像。第一遍利用掩模1左上角开始从左往右从上往下。第二遍利用第二个掩模右下角开始从右往左从下往上。掩模形状如下图所示 按照某种距离如 D 4 D_4 D4 距离或 D 8 D_8 D8 距离对大小为 M × N M×N M×N 的图像中的区域块作距离变换算法过程如下
1 建立一个大小为 M × N M×N M×N 的数组 F F F作如下的初始化将区域块中的元素设置为 0其余元素设置为无穷
2 利用掩模1mask1左上角开始从左往右从上往下遍历数组将掩模中P点对应的元素的值作如下更新 F ( P ) min q ∈ mask1 { F ( P ) , D ( P , q ) F ( q ) } F(P)\min _{q \in \operatorname{mask1}}\{F(P), D(P, q)F(q)\} F(P)q∈mask1min{F(P),D(P,q)F(q)}
3 利用掩模2mask2右下角开始从右往左从下往上遍历数组将掩模中P点对应的元素的值作如下更新 F ( P ) min q ∈ mask2 { F ( P ) , D ( P , q ) F ( q ) } F(P)\min _{q \in \operatorname{mask2}}\{F(P), D(P, q)F(q)\} F(P)q∈mask2min{F(P),D(P,q)F(q)}
最终得到的更新后的数组即为距离变换的结果。
这个算法过程在图像的边界处需要做出调整因为在边界处掩模不能全部覆盖图像这时可以将掩模中没有对应元素的位置的值当作 0 来处理即maskSize0。
在 OpenCV 中distanceTransform 函数是用于计算二进制图像中每个非零像素到最近零像素的距离的函数。这个函数通常用于图像处理中的形态学操作和特征提取。下面是 distanceTransform 函数的一般形式
dist_transform cv2.distanceTransform(src, distanceType, maskSize) src 是输入的二进制图像该图像应该是一个 8 位单通道图像。 distanceType 是指定距离度量的类型。 maskSize 是指定计算距离时使用的卷积核的大小。
# 分别利用欧式距离、D4 距离和 D8 距离作距离变换
transMatE cv2.distanceTransform(mat, distanceTypecv2.DIST_L2, maskSize0) # 计算欧氏距离变换
transMatD4 cv2.distanceTransform(mat, distanceTypecv2.DIST_L1, maskSize0) # 计算 D4 距离变换
transMatD8 cv2.distanceTransform(mat, distanceTypecv2.DIST_C, maskSize0) # 计算 D8 距离变换# 输出欧式、D4 和 D8 的距离矩阵
print(欧氏距离的变换矩阵:\n, transMatE)
print(城区距离-D4的变换矩阵:\n, transMatD4)
print(棋盘距离-D8的变换矩阵:\n, transMatD8)四、将距离矩阵转换为可视化图片
因为经过距离矩阵变换之后变换结果的数据类型为 float3232 位浮点数而在 OpenCV 中采用 imshow 函数显示图像时需要使用 uint8 数据类型8 位无符号整数使得像素值的范围是从 0255可以表示灰度图像中的所有可能像素值。而 Matplotlib 可以直接处理 float32 类型的图像数据因此并不需要将图像数据类型转换为 uint8 类型。
plt.figure(figsize(15, 5))plt.subplot(1, 3, 1)
plt.imshow(transMatE, cmapgray)
plt.title(欧氏距离, fontsize16)
plt.colorbar(shrink0.8)plt.subplot(1, 3, 2)
plt.imshow(transMatD4, cmapgray)
plt.title(城区距离-D4, fontsize16)
plt.colorbar(shrink0.8)plt.subplot(1, 3, 3)
plt.imshow(transMatD8, cmapgray)
plt.title(棋盘距离-D8, fontsize16)
plt.colorbar(shrink0.8)plt.tight_layout()
plt.show()
