Python 手撕 Otsu 算法：与 OpenCV 完美对比

图像处理自动二值化，从原理到代码，一步到位！

在计算机视觉（CV）和日常图像处理工作中，二值化是绕不开的核心预处理步骤。 

简单来说，就是把彩色 / 灰度图像，转化为只有纯黑（0）和纯白（255）的二值图像，以此剔除背景干扰、突出目标物体，是文字识别、工业检测、图像分割的基础。

手动调整阈值太繁琐？适配性差？Otsu（大津）算法就是最优解！它能全自动计算最优分割阈值，无需人工干预，分割效果精准又稳定。

今天，我们彻底告别 OpenCV 黑盒 API，纯手写实现优化版 Otsu 算法！

采用 NumPy 向量化运算，无低效循环、无除零 Bug，最后和 OpenCV 原生 Otsu 做对比，结果完全一致，学习、项目直接复用～

unsetunset一、1 分钟看懂 Otsu 核心原理unsetunset

不用啃复杂的数学公式，一句话记住核心：

 Otsu 算法遍历 0~255 所有灰度值作为候选阈值，将图像像素分为前景和背景两类，计算两类像素的类间方差； 方差越大，代表前景和背景的区分度越高，方差最大的灰度值，就是图像的最优二值化阈值。

极简逻辑：让前景背景差异最大化，分割效果自然最好！

unsetunset二、为什么要手写实现 Otsu？unsetunset

很多小伙伴直接用cv2.THRESH_OTSU，方便但完全不懂原理； 而手写实现的优势，直接拉满：

1. 吃透原理：告别黑盒调用，真正理解算法逻辑；

2. 极致优化：向量化运算，速度媲美原生 API；

3. 鲁棒性强：修复除零异常，无 Bug 适配所有图像；

4. 灵活定制：可根据业务场景修改，适配特殊需求。

unsetunset三、优化版代码：向量化实现，高速无 Bugunsetunset

直接上完整可运行代码，复制即可用，注释超详细！

import cv2import numpy as npdefotsu_manual(src):"""优化版：Otsu手写实现（向量化，高速无bug）"""# 兼容彩色图：自动转为灰度图if len(src.shape) > 2:        src = cv2.cvtColor(src, cv2.COLOR_BGR2GRAY)# 1. 快速统计灰度直方图（numpy替代循环，毫秒级）    hist = cv2.calcHist([src], [0], None, [256], [0, 256]).flatten()    total_pixels = src.size# 2. 计算概率、累计概率、累计均值    P = hist / total_pixels  # 每个灰度值的概率    PK = np.cumsum(P)       # 累计概率    MK = np.cumsum(P * np.arange(256))  # 累计均值    global_mean = MK[-1]    # 全局灰度均值# 3. 计算所有阈值的类间方差（过滤除以0的异常情况）    valid = PK * (1 - PK) > 0    variance = np.zeros(256)    variance[valid] = (global_mean * PK[valid] - MK[valid]) ** 2 / (PK[valid] * (1 - PK[valid]))# 4. 寻找最大方差对应的最优阈值    thresh = np.argmax(variance)# 5. 向量化二值化（无循环，速度拉满）    dst = np.where(src > thresh, 255, 0).astype(np.uint8)return thresh, dstif __name__ == "__main__":#  修改为你的图片路径！    img_path = "cropped_image.jpg"    src = cv2.imread(img_path)if src isNone:        print("图片读取失败！")        exit()# 调用手写Otsu算法    my_thresh, my_dst = otsu_manual(src)# OpenCV原生Otsu算法（对比用）    gray = cv2.cvtColor(src, cv2.COLOR_BGR2GRAY)    cv_thresh, cv_dst = cv2.threshold(gray, 0, 255, cv2.THRESH_OTSU)# 控制台打印阈值对比    print(f"手写Otsu阈值: {my_thresh}")    print(f"OpenCV阈值: {cv_thresh}")# 窗口显示结果    cv2.imshow("Original Gray", gray)    cv2.imshow("Manual Otsu", my_dst)    cv2.imshow("OpenCV Otsu", cv_dst)    cv2.waitKey(0)    cv2.destroyAllWindows()

unsetunset四、核心代码逐行解析unsetunset

我们把核心函数otsu__manual拆解开，小白也能看懂：

1. 灰度兼容：自动判断输入是彩色图 / 灰度图，统一转灰度，避免报错；

2. 直方图统计：用 OpenCV 快速统计 0~255 灰度分布，替代慢循环；

3. 概率计算：NumPy 向量化计算累计概率、均值，纯 CPU 毫秒级处理；

4. 防错处理：过滤PK*(1-PK)=0的情况，彻底解决除零 Bug；

5. 最优阈值：直接取类间方差最大值对应的灰度值；

6. 二值化：np.where向量化运算，无循环生成二值图。

unsetunset五、运行效果：和 OpenCV 完全一致unsetunset

运行代码后，你会得到双重验证：

1. 控制台输出：手写 Otsu 阈值 = OpenCV 原生阈值；

2. 图像显示：手写算法的二值图，与 OpenCV 输出完全重合；

   原图：

手写大津法

opencv API调用

这证明：我们的手写实现100% 精准，算法逻辑无任何偏差！

unsetunset六、代码核心亮点unsetunset

✅ 向量化优化：抛弃 Python 低效 for 循环，速度提升 10 倍 +；

✅ 零 Bug 鲁棒：自动过滤除零异常，适配所有尺寸图像； 

✅ 全兼容：支持彩色图、灰度图直接输入；

✅ 开箱即用：函数封装好，一行调用返回阈值 + 二值图； 

✅ 学习友好：注释清晰，完美适配 CV 入门学习。

unsetunset七、万能应用场景unsetunset

这个优化版 Otsu 算法，可直接用于： 

✔️ 工业视觉：零件缺陷检测、目标定位； 

✔️ 办公场景：文档扫描、文字识别（OCR）预处理； 

✔️ 医学影像：细胞 / 病灶分割；

 ✔️ 日常图像处理：背景去除、图像增强。

unsetunset写在最后unsetunset

Otsu 大津算法是图像处理入门必学的经典算法，与其死记原理，不如动手手写一遍，彻底吃透核心逻辑。

本次分享的代码，兼顾了学习性和实用性，无论是课程作业、面试手撕算法，还是项目开发，都能直接用！