做图像处理这行，要是没跟“阈值分割”打过交道，基本等于没入门。这个操作听起来简单——把灰度图变成黑白的，目标就是背景，背景就是目标，像素值一比较就完事。但真到了实际项目里，光照一变、目标一复杂，固定阈值立刻翻车，背景噪点、目标缺失、边缘断裂全来了。我最近把几种主流阈值分割方法放在同一个测试框架里做了一轮对比实验，用MATLAB实现，跑了几组典型图像，把每种方法的表现、参数敏感度和适用边界都摸了一遍。这篇文章就复盘这次“泾渭分明”微实验的完整思路和实操细节，附核心代码，直接可以拿去改。

这个实验适合刚接触图像分割的初学者，也适合那些跟我一样被光照不均、目标对比度低折磨过的老手。我会把五种方法——固定全局阈值、Otsu大津法、迭代阈值、自适应局部阈值、多阈值分割——的数学原理、MATLAB实现、参数调优和踩坑点全部摊开讲。看完之后你至少能回答三个问题：什么场景该用哪种方法？为什么Otsu在某些图像上会失灵？自适应阈值窗口大小到底怎么定？

1. 实验整体设计与思路拆解

1.1 为什么要做“对比实验”而不是单个方法测试

很多人在学习阈值分割时，会陷入一个误区：拿到图像，直接用graythresh算个阈值，然后二值化，发现效果不好就换一张图，或者调一调参数碰运气。这种做法最大的问题是——你根本不知道当前图像适不适合这个方法，也不知道方法失效是因为参数不对还是原理本身就不匹配。

所以我整个实验的设计核心，是把“方法”、“参数”、“图像特征”三个变量拆开，做交叉对比。测试集刻意选择了四类图像：高对比度合成图、自然场景灰度图、光照不均文档图、低信噪比显微图。这四类基本覆盖了阈值分割在实际项目里遇到的典型情况。每张图都预先人工标注了标准二值掩膜，作为计算IoU、像素精度等指标的“标准答案”。

这个设计看起来简单，但它直接决定实验结论的可靠性。如果你只用一两张图得出结论，很容易以偏概全。比如在合成图上，固定阈值和Otsu几乎没差别，但放到光照不均的文档图上，自适应阈值又和全局方法拉开明显差距。没有多图像交叉验证，你看不到方法之间的真实分水岭。

1.2 五种方法的选型逻辑

选哪几种方法对比，也有一点讲究。我没有一开始就上深度学习方法，因为阈值分割本身就是传统图像处理的基石操作，搞清楚这些基础方法的脾气，才能理解后续更复杂的算法为什么那样设计。

• 固定全局阈值：最朴素的基线方法。所有像素共享一个阈值，大于阈值的记为目标，小于记为背景。实现简单，运行快，但对光照和噪声几乎没有抵抗力。
• Otsu大津法：自动寻找阈值，使得阈值两侧类间方差最大。解决了“阈值怎么定”的主观性问题，是实际项目里最常用的自动全局阈值方法。
• 迭代阈值法：跟Otsu思路不同，它从初始阈值出发，通过“分割→求均值→更新阈值”的方式反复迭代逼近一个稳定值。
• 自适应局部阈值：不再用单一全局阈值，而是在每个像素周围取一个局部窗口，用窗口内的统计信息计算该像素的阈值，专门对付光照不均。
• 多阈值分割：把图像分成多个灰度区间，对应不同类别，本质上是把二值分割扩展到多类分割。

这五条路线从全局到局部、从二值到多类，刚好构成一条完整的阈值分割技术谱系。实验时让它们在同一测试集上跑，各自的优势、短板就全部暴露出来了。

1.3 环境配置与实验流程

我的实验环境是MATLAB R2022a，只用了Image Processing Toolbox，没有额外装第三方库。核心函数涉及imbinarize、graythresh、adaptthresh、multithresh，另外自己手写了迭代阈值算法，方便观察每一步的收敛过程。

实验流程分为四步：先读图并转灰度，归一化到[0,1]区间；然后调用不同方法生成二值掩膜；第三步计算量化指标（IoU、精确率、召回率、F1）；最后输出可视化和指标表格。下面会按这四步逐个展开。

2. 关键原理剖析与参数选择思路

2.1 Otsu的类间方差到底在算什么

Otsu方法的核心思想是：假设图像灰度直方图有两个峰，分别对应目标和背景。我们要找一个分割阈值，让两类的类间方差最大，也就是让“峰与峰之间的距离”尽可能大。类间方差公式是：

[ \sigma_b^2 = \omega_1(t) \omega_2(t) [\mu_1(t) - \mu_2(t)]^2 ]

其中(\omega_1)和(\omega_2)分别是目标像素和背景像素占全图的比例，(\mu_1)和(\mu_2)是它们的灰度均值。这个式子的直观含义是——当两个类内部越紧凑、彼此中心隔得越远，方差就越大，分裂效果越好。

MATLAB里graythresh返回的就是使类间方差最大的归一化阈值，范围在[0,1]，配合imbinarize(gray, threshold)使用即可。这里有一个很多人踩过的坑：graythresh输入的是uint8原始灰度图或[0,1]范围的double图，结果不同。如果输入uint8，返回的阈值也是0到1之间的归一化小数，而如果你自己写循环遍历阈值，要注意灰度级是0到255，二者需要换算清楚。

2.2 迭代阈值法的收敛逻辑

迭代法是一个经典的“猜测→修正”循环。它先取全图灰度均值或中值作为初始阈值T0，然后按T0把图像分成两组G1和G2，分别计算两组灰度均值m1和m2，新阈值就是(m1+m2)/2。重复这个过程，直到连续两次阈值的差小于某个极小值，比如1e-6，或者达到最大迭代次数。

从实验观察来看，迭代法在大多数情况下的收敛速度很快，通常十几步就能稳定。它比Otsu更“稳健”的一点是，结果不依赖直方图双峰假设的严格成立，哪怕直方图不光滑，迭代过程也会逐渐逼近一个中间态。缺点是它仍然依赖全局统计信息，光照不均的场景下照样无能为力。

2.3 自适应阈值：局部窗口的尺度博弈

自适应阈值计算每个像素的阈值时，用的是一个局部窗口，窗口大小和偏移量C直接决定效果。MATLAB的adaptthresh函数里，最重要的两个参数是NeighborhoodSize和Offset。

NeighborhoodSize的含义是局部统计的邻域范围，推荐值是图像宽度的1/10到1/20。窗口太小，局部统计受噪声影响大，二值化结果会出现大量细碎斑点；窗口太大，局部结构又趋向于全局统计，失去了局部自适应的意义。我在实验里对文档图试过从5到300像素不等的窗口，结论是窗口大约在目标特征尺寸的2到3倍时效果最好。Offset则是一个手工修正项，正数让阈值偏高、分割结果偏暗，负数让阈值偏低、结果偏亮。它本质上是在召回率和精确率之间做权衡——Offset取正数，能滤掉浅色伪影，但也可能丢掉真正的弱目标。

2.4 多阈值分割的使用边界

multithresh函数做的是把灰度直方图分成多段，每段对应一个类别。这个方法的实用性经常被高估。对自然图像，三个或四个灰度区间往往就能覆盖天空、地面、目标物，但一旦图像内的光照变化剧烈，单纯按灰度分层就会把阴影区域和真正的暗目标混为一谈。

我通常建议把它用在那些“目标本身存在不同灰度层次”的场景。比如电路板检测中，焊盘、走线、基板分别处于不同的灰度区间，用三级分割可以同时分离出三类。如果只是想把目标从背景里挖出来，强行用多阈值只会增加后续处理的复杂度。

3. 实操过程与核心代码实现

3.1 读取数据与预处理

实验的第一步是统一灰度归一化。这里要强调一个细节：MATLAB里imread读进来的如果是彩色图，建议用rgb2gray转灰度，而不是直接取某一个颜色通道，因为不同通道的光照响应不同，直接影响后续阈值效果。

img = imread('test_image.png'); if size(img, 3) == 3 gray = rgb2gray(img); else gray = img; end gray = im2double(gray); % 将灰度值归一到[0,1]

这里用im2double而不是double(gray)，是因为im2double会自动把uint8的数据除以255，避免后续算法里单位不统一的问题。归一化之后，所有方法的阈值参数都在同一尺度下，可比较性强。

3.2 各方法的代码与参数设置

固定全局阈值取0.5，作为基线：

bw_fixed = imbinarize(gray, 0.5);

Otsu方法：

level = graythresh(gray); bw_otsu = imbinarize(gray, level);

迭代法手写实现，方便观察阈值收敛过程：

function [bw, T_final] = iterative_threshold(gray, tol) if nargin < 2 tol = 1e-6; end T0 = mean(gray(:)); while true g1 = gray(gray > T0); g2 = gray(gray <= T0); if isempty(g1) || isempty(g2) break; end m1 = mean(g1); m2 = mean(g2); T_new = (m1 + m2) / 2; if abs(T_new - T0) < tol break; end T0 = T_new; end bw = gray > T0; T_final = T0; end

这段代码里有个很容易被忽略的细节：gray(gray > T0)抽取的是一维线性索引，而不是保留图像的二维结构，但因为我们只用它来计算均值，所以线性索引完全够用。如果后续需要对前景和背景分别做形态学处理，那就需要保留二维掩膜了。

自适应阈值方法：

T_adapt = adaptthresh(gray, 0.4, 'NeighborhoodSize', 2*floor(size(gray,1)/16)+1, 'Offset', 0.05); bw_adapt = imbinarize(gray, T_adapt);

adaptthresh的第一个参数0.4是灵敏度，取值范围[0,1]，它控制的是局部阈值相对于局部均值的偏移程度，数值越小分割越激进，越容易把暗区域判为目标。NeighborhoodSize我取的是图像高度1/8的奇数倍，这样局部统计区域不会太小也不会太大。

多阈值三级分割：

thresh_levels = multithresh(gray, 2); % 2个阈值，划分3类 bw_multi1 = gray > thresh_levels(1) & gray <= thresh_levels(2); bw_multi2 = gray > thresh_levels(2); bw_multi = bw_multi1 | bw_multi2; % 将中间和亮部合并为目标

这里要注意：multithresh返回的阈值向量是升序排列的，第一个阈值把图像分成暗部和亮部，第二个阈值进一步把亮部细分成中间灰和高亮。实验里我把中间灰和高亮都归为目标，把最暗部分作为背景，这样的合并策略是为了和人工标注保持一致的语义。

3.3 评价指标计算方法

量化对比必须用统一指标，我选了四个：IoU（交并比）、Precision（精确率）、Recall（召回率）、F1分数。其中IoU最严格，它同时惩罚误检和漏检，是最常用的分割质量指标。

function [precision, recall, f1, iou] = calc_metric(bw, gt) bw = logical(bw); gt = logical(gt); TP = sum(bw(:) & gt(:)); FP = sum(bw(:) & ~gt(:)); FN = sum(~bw(:) & gt(:)); precision = TP / (TP + FP + eps); recall = TP / (TP + FN + eps); f1 = 2 * precision * recall / (precision + recall + eps); iou = TP / (TP + FP + FN + eps); end

eps在这里用来防止除零错误。我特别提这个是因为真实图像里经常出现某种方法把整个图像都判为背景的情况，TP + FP可能等于0，不加保护直接NaN。

3.4 实验数据观测记录

我在高对比度合成图上跑完五种方法，IoU结果如下表：

这个表格信息量很大。高对比度合成图上五种方法都表现良好，因为图像本身“泾渭分明”；一换到自然图，光照因素开始介入，Otsu和自适应方法有明显优势；光照不均的文档图是自适应阈值的绝对主场，全局方法的IoU全线跌破0.35；低信噪比环境下，自适应阈值没有想象中占绝对优势，只是略好于Otsu。这说明信噪比低时，局部统计同样会被噪声污染，不是加了“局部”两个字就万能了。

3.5 可视化与结果输出

实验中的每张图我都保存了原图灰度图、人工标注、五种方法的分割结果，叠加对比图用imfuse函数生成半透明叠加效果，并在图上标注IoU数值。imfuse的blend模式对肉眼检查分割区域非常管用，比单独看黑白二值图像直观得多。

overlay = imfuse(gray, bw_otsu, 'blend'); figure; imshow(overlay); title(sprintf('Otsu result, IoU = %.3f', iou_otsu));

叠加热力图的方式也推荐一个：先计算分割边界bwperim(bw_otsu)，然后把边界绘制在原始图像上，可以精确看出边缘偏移发生在哪些位置。这一步对分析目标边缘模糊、细碎背景干扰这类问题，比只看IoU数字有用得多。

4. 常见问题与排查技巧实录

4.1 Otsu在单峰直方图上完全失效

这是实验中最常见的问题。当图像中目标面积占比极小时，灰度直方图几乎是单峰——大量背景像素堆在一起，目标像素占比太低，对类间方差的贡献可以忽略。此时graythresh算出来的阈值会偏向背景峰，把目标和一部分背景一起切成前景或背景。

我在一组细胞显微图像上就遇到了这种情况。目标细胞在整个视野中面积占比不到5%，Otsu分割后大部分细胞都被漏掉，IoU只有0.35左右。解决办法是先用形态学顶帽变换imtophat抑制背景，再执行Otsu分割，或者改用局部自适应方法。顶帽变换能去除大面积缓慢变化的背景光照，效果立竿见影。具体参数上，结构元素大小取目标直径的1.5到2倍左右比较合适。

4.2 自适应阈值窗口大小的选择逻辑

窗口选多大的问题，实验里我从5像素试到500像素，经历了三个明显阶段：窗口太小时，图像内部纹理被误判为不同区域，分割结果变成“胡椒盐”噪声图；窗口适中时，目标和背景被干净分开；窗口过大时，局部统计逐渐退化为全局统计，光照不均问题重新出现。

判断窗口是否合适的经验法则是：用目标在图像中的最大尺寸作为参考，窗口大小取该尺寸的2到4倍。如果目标本身没有固定尺度，就先用形态学开运算估计背景照度分布，再用背景照度平滑尺度作为窗口参考。还有个小技巧是，NeighborhoodSize一定要设置为奇数，否则MATLAB会报错。

4.3 噪声敏感与预处理取舍

阈值分割本身没有去噪能力，它只会忠实地把噪声像素也分到某一类。实验里低信噪比图的对比最明显——固定阈值方法直接在噪声点周围产生大量孤立前景像素，IoU掉到0.3附近。

要预处理还是不要预处理，取决于噪声的类型和后续应用。如果是高斯噪声，先跑一个imgaussfilt(gray, 1.5)能显著提升分割稳定性；如果是脉冲噪声（椒盐噪声），用medfilt2中值滤波更有效，窗口3×3就够了。但要注意，滤波会损失边缘细节，对边缘有严格要求的测量场景，宁可接受噪声导致的零星误检，也不要预先平滑。这里没有银弹，必须根据下游任务取舍。

4.4 代码层面的几个隐蔽坑

imbinarize和graythresh配合时，imbinarize默认是global方法，用graythresh返回的阈值是没有问题的。如果你直接写gray > level进行比较，千万确认gray和level的数据类型一致——gray是double，level是double，没有问题；但如果gray是uint8，level是小数值，两者直接比较会触发数据类型转换，结果完全不是你想要的。

此外，用adaptthresh返回的阈值矩阵给imbinarize作为第二参数时，imbinarize会自动识别为自适应分割，不需要额外指定方法。但如果你把adaptthresh的返回结果另存后，修改了矩阵尺寸再传入imbinarize，会出现维度不匹配的报错。调试时改来改去就容易踩这种坑。

5. 实验之外的一些经验补充

这个实验做完后，我最大的感触是：阈值分割不存在“最好的方法”，只存在“最合适的方法”。在实际项目里，决定用哪种方法之前，先花几分钟判断图像的光照分布是否均匀、目标占比是否过小、噪声是否严重、纹理是否复杂。这个判断做对了，方法选择就是水到渠成的事。

还有一个经验：把adaptthresh的Offset从0开始按0.02步进扫描，观察IoU曲线的走向，经常能帮你找到当前图像的最佳修正量。这个思路比拍脑袋设参数靠谱得多，而且完全可以写成一个自动调参的循环。

代码仓库里我保留了实验脚本和测试图像，后续还想补充两部分内容：一是把multithresh的类别数自动选择加进来——用类间方差或BIC准则对类别数做评估；二是把三种方法封装成统一的函数接口，方便直接部署到实时图像处理流程中。做实验就要做成能复用的工具，而不是跑完就扔的脚本。

我个人的体会是，这种“小实验”最大的价值不是得出结论本身，而是通过亲手调参，建立对算法行为的直觉。有了这种直觉，再看那些复杂的分割网络，你就能理解它为什么要这样设计，以及它在什么情况下会重蹈传统方法的覆辙。这个过程，比背一百篇论文都管用。