在轮廓分析中，获取轮廓的中心点（质心）是目标定位、尺寸测量的核心步骤。本文通过 OpenCV 的图像矩（cv.moments()）实现轮廓中心点的精准计算与可视化，新手可直接复用完整流程。
核心代码实现
python
import cv2 as cv
import numpy as np
# 1. 读取图像并校验
src = cv.imread('./image/21.bmp')
if src is None:
print('图像读取失败，请检查路径！')
exit()
# 2. 预处理：高斯滤波去噪 + 二值化
dst = cv.GaussianBlur(src, (3, 3), 15) # 高斯滤波去除噪点，避免伪轮廓干扰
grays = cv.cvtColor(dst, cv.COLOR_BGR2GRAY) # 转灰度图
# 二值化（分离前景轮廓与背景，阈值200可根据图像调整）
ret, binary = cv.threshold(grays, 200, 255, cv.THRESH_BINARY)
# 3. 查找轮廓（CCOMP模式：检索内外层轮廓，适配有孔洞的目标）
contours, hierarchy = cv.findContours(binary, cv.RETR_CCOMP, cv.CHAIN_APPROX_SIMPLE)
# 4. 遍历轮廓，计算并绘制每个轮廓的中心点
for i, contour in enumerate(contours):
# 绘制轮廓（蓝色，线宽5，便于直观查看轮廓范围）
cv.drawContours(src, contours, i, (255, 0, 0), 5)
# 计算轮廓的图像矩（核心：通过矩获取中心点）
M = cv.moments(contour) # 修正原代码：使用当前轮廓contour，而非固定contours[1]
# 避免除零错误（过滤无效小轮廓，m00为零表示轮廓无有效像素）
if M["m00"] != 0:
# 计算中心点坐标（x = m10/m00，y = m01/m00）
center_x = int(M["m10"] / M["m00"])
center_y = int(M["m01"] / M["m00"])
# 绘制中心点（绿色圆形，半径10，线宽2，醒目易见）
cv.circle(src, (center_x, center_y), 10, (0, 255, 0), 2)
# 打印中心点坐标，便于数据分析
print(f"轮廓{i} 中心点坐标：x={center_x}, y={center_y}")
else:
print(f"轮廓{i} 为无效轮廓，跳过中心点计算")
# 5. 显示结果图像
cv.namedWindow('src', cv.WINDOW_NORMAL)
cv.resizeWindow('src', 600, 600)
cv.imshow('src', src)
cv.waitKeyEx(0)
cv.destroyAllWindows()关键知识点解析
1. 核心流程拆解
| 步骤 | 核心 API | 作用说明 |
|---|---|---|
| 图像矩计算 | cv.moments(contour) |
计算轮廓的几何矩，包含轮廓的位置、面积等信息 |
| 中心点求解 | center_x = m10/m00、center_y = m01/m00 |
m00是轮廓面积矩，m10/m01是一阶矩，通过比值求解中心点 |
| 中心点可视化 | cv.circle() |
用圆形标注中心点，直观验证计算结果 |
| 无效轮廓过滤 | M["m00"] != 0 |
避免小噪点轮廓导致的除零报错，提升代码健壮性 |
2. 核心概念说明
- 图像矩（Moments） ：描述图像轮廓几何特征的数值集合，OpenCV 中
cv.moments()返回的矩包含多个关键字段，核心用于中心点计算的仅有 3 个：M["m00"]：零阶矩，对应轮廓的面积（非零表示轮廓有效）；M["m10"]：一阶矩，用于计算中心点 x 坐标；M["m01"]：一阶矩，用于计算中心点 y 坐标。
- 轮廓索引 ：原代码固定使用
contours[1]，仅能计算第 2 个轮廓的中心点，改为contour遍历每个轮廓，适配多轮廓场景。
3. 避坑与优化技巧
- 除零错误 ：必须添加
M["m00"] != 0的判断，否则小噪点轮廓会导致ZeroDivisionError； - 阈值调整：二值化阈值 200 需根据图像亮度调整，暗轮廓可降低阈值（如 150），亮背景可提高阈值（如 220）；
- 多轮廓适配 ：遍历
enumerate(contours)而非固定索引，支持图像中多个目标轮廓的中心点批量计算； - 精度优化 ：若需更高精度的中心点，可保留浮点数（去除
int()转换），适用于精密测量场景。
总结
- 轮廓中心点计算的核心是图像矩（
cv.moments()） ，通过m10/m00和m01/m00求解 x、y 坐标； - 必须添加
M["m00"] != 0的判断，避免除零错误，提升代码健壮性； - 遍历所有轮廓而非固定索引，才能实现多目标轮廓的中心点批量提取与可视化。