简介本资源是一套面向材料科学、表面物理及实验数据分析领域的MATLAB工具包聚焦接触角测量与固体表面能反演这一核心问题适用于高校研究生、科研人员及工程技术人员开展润湿性表征与界面性质研究。压缩包共含10个.m文件总大小仅6KB以主程序Contactangle.m为核心辅以多个fpush/fpop系列函数构成完整的图像边缘提取、液滴轮廓拟合与接触角计算流程支持基于Young-Laplace方程及Wenzel模型的表面能估算。已有576人学习下载体现了其在实验数据处理中的实用价值。用户可直接调用脚本完成从原始图像读取、预处理、接触点识别到角度输出与表面能推算的全流程分析代码结构清晰、模块分工明确且内置基础错误处理机制便于二次开发与适配不同实验条件下的液滴图像分析需求。 做材料表面分析的人几乎都遇到过这类场景从某个渠道下载了一个压缩包文件名写着contact angle.rar里面躺着一个contact_angle.m附带一两张液滴照片。网络热词里也总能看到matlab、contact angle这对组合说明大家都在找能直接跑通的接触角计算程序。但根据我的经验这类程序十有八九第一次运行不是报错就是结果离谱角度算出来 90 度上下浮动连自己都怀疑是不是量角器都比它准。这篇文章不是给你一份现成代码让你复制粘贴就完事而是把这套程序从能运行变成算得准的全过程拆开讲清楚文件里到底有什么、图像预处理哪些步骤省不得、接触角的几何定义是怎么落到atan和polyfit这些函数上的、以及跑完程序之后你凭什么相信那个数字。内容面向正在做接触角测量实验、需要批量处理液滴图片的工程师和科研人员你不需要精通图像处理但看完应该能独立把一套 MATLAB 程序改造得服服帖帖。1. 解压后的第一件事先搞清这个程序到底在算什么1.1 接触角测量到底测的是什么接触角contact angle指的是在固、液、气三相交界处从固液界面经过液滴内部到气液界面切线所夹的角度。这个角度直接反映固体表面的润湿性小于 90 度说明表面亲液大于 90 度说明疏液。测量方式看着简单——把一滴液体滴在固体表面从侧面拍张照片量一下就行了。但真到程序实现时那张照片里的信息量远比想象中多液滴轮廓是一条曲线固体表面是一条直线接触角是两者交点的局部几何关系而不是随便量一个整体角度。很多网上下载的 MATLAB 程序核心逻辑其实就三步读图、找轮廓、算切线。只不过这三步每一步都有不少隐藏细节。contact_angle.m这个脚本通常不是一个人写的它可能合并了多种方法甚至留有原作者实验条件下的硬编码参数。你直接跑很可能算的是某个特定背光条件下、某个特定固体表面的角度换一张图就失效。1.2 文件构成与运行环境检查解压contact angle.rar之后典型的文件结构是这样的contact_angle.m主脚本compute_contact_angle.m子函数负责从边缘点坐标计算角度droplet_sample.png或.bmp示例液滴图片可能还有一个README.txt或说明文档运行前先检查环境这步能省掉一半的报错。contact_angle.m普遍依赖 Image Processing Toolbox用到imread、rgb2gray、edge、bwareaopen等函数和 Curve Fitting Toolbox用到polyfit、polyval。在命令行敲ver查看工具箱列表缺了就去装否则一运行直接弹Undefined function edge这种错误很多人卡在这第一关。另一个容易忽略的是路径。把整个文件夹放在纯英文路径下运行别放桌面带中文的文件夹里。MATLAB 对中文路径的支持时好时坏imread偶尔抽风读不出文件排查起来非常浪费时间。我一般还会在脚本开头加上clc; clear; close all;避免上次运行的变量残留在工作区里干扰结果。提示如果解压后没有示例图片先别急着跑程序。找一张自己拍的液滴照片背景尽量简单、液滴轮廓清晰否则你根本分不清是代码问题还是图片问题。2. 图像预处理这步决定成败灰度化、降噪、背景分离2.1 为什么要先转灰度再做边缘检测相机拍出来的液滴照片是彩色图每个像素有 R、G、B 三个通道。直接对彩色图做边缘检测不是不行但三个通道分别响应边缘会变得极不均匀——红光下液滴边缘清楚蓝光下可能糊成一片。标准做法是先转成灰度图I imread(droplet_sample.png); if size(I, 3) 3 I rgb2gray(I); end这一步把三维矩阵压缩成二维之后的阈值分割、边缘检测全部基于亮度信息计算量也小得多。实测下来一张 1200 万像素的照片灰度化后处理速度能提升三倍左右。2.2 滤波窗口选多大才不会弄丢边缘转灰度之后紧接着是降噪。液滴照片的噪声来源通常有两类一是相机传感器的高斯噪声二是背景里的微小颗粒或灰尘产生的椒盐噪声。很多人一上来就用高斯滤波做平滑结果液滴边缘被抹得圆乎乎接触角算出来系统性地偏大或偏小。我的习惯是先做中值滤波medfilt2窗口取 3×3最多 5×5。中值滤波对椒盐噪声特别有效而且能较好地保留边缘的陡峭程度。高斯滤波imgaussfilt不是不能用但 sigma 超过 1.5 就容易把液滴顶部的小弧度细节吃掉接触角测量对轮廓局部的切线方向极其敏感边缘偏移哪怕一两个像素角度就能差出好几度。I_filt medfilt2(I, [3 3]);这里有个容易被忽略的点滤波窗口越大边缘越钝。做了 5×5 中值滤波之后液滴边缘线会被轻微腐蚀算出来的接触角往往比真实值偏大。所以滤波的原则是能去掉噪声就行窗口尽量小。如果原图噪声不大直接跳过滤波也行。2.3 自动阈值和 ROI 到底该信谁接下来是二值化把液滴从背景里分离出来。最简单的办法是imbinarize配合graythresh自动阈值BW imbinarize(I_filt, graythresh(I_filt));graythresh用的是 Otsu 方法它会自动找一个阈值让前景和背景的类间方差最大。这个方法对背景均匀、液滴和背景对比明显的照片效果很好。问题在于实际拍摄环境很难这么理想。如果光线从侧面打过来液滴表面会有镜面反射反光区域亮度接近背景Otsu 会把那一块误判成背景导致液滴轮廓缺个口子。这种情况下我不会完全信任自动阈值而是先用imcrop把图像裁剪到液滴附近的区域。裁剪有两个好处一是去掉画面边缘无关的高亮区域让 Otsu 的阈值更准二是大幅减少后续边缘检测的计算量。裁剪时不要裁得太紧液滴下方要留出足够的固体表面区域方便后面确定基线。裁剪之后再次imbinarize或者edge成功率会明显提升。如果背景依然不均匀比如有渐变阴影可以试试imbinarize(I_filt, adaptive)它按局部区域计算阈值对光照不均的图片比全局阈值稳得多。代价是速度慢而且参数ForegroundPolarity需要根据液滴是亮还是暗来调整。3. 从二值图到液滴轮廓边缘检测与轮廓筛选3.1 Canny 阈值怎么调才靠谱二值图拿到手下一步是提取液滴边缘。MATLAB 的edge函数支持多种算子算接触角我强烈推荐 Canny。Sobel 算子边缘粗、对噪声敏感算切线方向容易被局部毛刺带偏Canny 做了非极大值抑制和双阈值边缘细而且连续性好。edges edge(BW, Canny, [0.1 0.2]);双阈值的理解很简单低阈值决定边缘检测的灵敏度高阈值决定保留哪些强边缘。[0.1 0.2]这个组合是大多数液滴图片的常用起点。如果检测出来的边缘断裂严重说明高阈值太高降到[0.05 0.15]如果边缘旁边跟着一堆杂线说明低阈值太低往上提到[0.15 0.3]。有个经验技巧在调阈值之前先把BW显示出来肉眼确认液滴在二值图里是不是完整的一块。如果二值图本身就有缺口边缘检测调参怎么调都白搭。这时候应该回去调阈值分割或者用形态学操作把缺口补上。3.2 形态学处理与轮廓去杂二值图里的液滴经常不是完美的一块反光点会在液滴内部形成空洞灰尘会在液滴旁边形成孤立的小白点。直接用bwboundaries提取轮廓你会得到一串乱七八糟的边界只有其中一条才是液滴轮廓。形态学处理能解决大部分问题BW_filled imfill(BW, holes); % 填充液滴内部空洞 BW_clean bwareaopen(BW_filled, 500); % 去掉小于500像素的孤立区域 BW_closed imclose(BW_clean, strel(disk, 3)); % 连接断开的边缘bwareaopen里的 500 是经验值根据你的图片分辨率调整。如果液滴占了画面的大部分这个值可以设成 2000 甚至更大如果液滴本身就很小设 500 可能把液滴主体给删了。判断标准很简单处理后画面里应该只剩下一个最大的连通区域。轮廓提取用bwboundariesB bwboundaries(BW_closed, noholes);返回的B是一个元胞数组每个元素是一圈边界的坐标点N×2 的矩阵第一列是行号 y第二列是列号 x。液滴轮廓通常是其中最长的那个边界。找出来后用plot叠在原图上看一眼确认轮廓是不是完整贴合液滴的外形。3.3 基线固体表面线的确定方法接触角的定义里基线就是固体表面线。如果液滴是滴在一块平整的固体表面上那么照片里固体表面就是一条接近水平的直线。这条线的位置直接决定了接触点的坐标偏差哪怕半度都会让最终角度失真。确定基线有三种方法我按推荐程度排个序方法一霍夫变换检测直线。用houghlines在边缘图里找线段筛选角度接近 0 度水平的那条。这个方法自动但容易误检液滴底部的弧线在某些角度下会被误判成水平线。方法二手动指定两个点。用ginput在图像上点击固体表面的左右两个位置连成直线。这个最稳但不适合批量处理。方法三拟合液滴轮廓的最低点区间。这个方法我用得最多——液滴与固体表面接触的那一段轮廓与基线几乎重合。取轮廓点里 y 坐标最大的那部分点大概是底部 5% 到 10% 的点用polyfit(x, y, 1)做线性拟合拟合出的直线就是基线。idx find(B_pts(:,1) prctile(B_pts(:,1), 90)); p_base polyfit(B_pts(idx,2), B_pts(idx,1), 1);方法三的好处是完全自动且物理意义明确液滴底部被固体表面压平的那一段天然就是基线的一部分。唯一的风险是液滴太小、底部弧线太弯这时拟合出来的线会有点倾斜需要人工检查。基线确定后轮廓与基线的交点就是左右两个接触点。这两个点的 x 坐标一个是最小值、一个是最大值取轮廓点里离基线最近且 x 坐标分别最小和最大的点即可。4. 接触角计算的几何本质与三种实现方案4.1 图像坐标系下的角度陷阱接触角的核心计算其实就是一个初中几何问题曲线上某一点的切线斜率决定了该点切线的角度。但 MATLAB 图像坐标系有一个大坑——原点在左上角x 轴向右y 轴向下。这跟我们熟悉的笛卡尔坐标系正好相反导致atan算出来的角度方向全部翻转。举个例子液滴左侧轮廓从接触点向上走x 坐标逐渐增大接触点在左往右到液滴顶部y 坐标逐渐减小向上在图像坐标系里斜率 dx/dy 是负的。用atan算出来的角度是负值但是真实的物理角度是正的。这就是为什么很多新手程序算出来接触角一直在 0 到 10 度之间徘徊以为固体表面超亲液其实是坐标系符号没处理。正确处理思路先明确你要的角度是接触点处液滴轮廓切线与水平基线之间的夹角然后画个图理清几何关系。用atan2而不是atan它能正确处理四个象限。算出切线向量后再根据左右侧分别调整角度定义。4.2 方案一高宽比法θ/2 法——最快但有限制高宽比法也叫 θ/2 法是接触角计算里最朴素的方法。它假设液滴是一个球冠理想情况下液滴足够小重力影响可忽略此时接触角只跟液滴高度 h 和接触宽度 w 有关theta 2 * atan(2 * h / w)h 是液滴顶点到基线的垂直距离w 是左右接触点之间的水平距离。用 MATLAB 实现就是先找到轮廓的最高点再算高度和宽度。这个方法的优点是代码量十几行计算量几乎为零而且对噪声不敏感。缺点是假设条件太苛刻液滴必须足够小一般直径小于 2 mm才能忽略重力对形状的拉拽而且液滴必须左右对称。实际实验里液滴一旦稍微大一点轮廓就被重力压扁不再是球冠高宽比法算出来的角度偏小。我一般把 θ/2 法当快速预估工具用批量处理几十张图先跑一遍看看数据分布是否合理再对关键样本用更精确的方法复核。4.3 方案二多项式拟合与切线法——工程界的主力方案这是目前工程和科研里最常用的接触角算法也是contact_angle.m这类程序里最可能采用的方法。核心思路是把液滴左右两侧的轮廓分别提取出来用多项式拟合轮廓曲线然后在接触点处求导得到切线斜率再换算成角度。以左侧轮廓为例代码逻辑是这样的% 提取左轮廓点x 最小的那段轮廓 left_mask (B_pts(:,2) min(B_pts(:,2))) | ... (B_pts(:,2) prctile(B_pts(:,2), 5)); % 这里只是示意实际要按连通区间取左边缘点 x_left B_pts(left_mask, 2); y_left B_pts(left_mask, 1); % 二次多项式拟合 y a*x^2 b*x c p polyfit(x_left, y_left, 2); % 在接触点处求导dy/dx 2*a*x b x_contact min(x_left); slope polyval(polyder(p), x_contact); % 切线角度弧度转角度 theta_rad atan(slope); % 图像坐标系翻转修正左接触角 theta_left 180 - abs(rad2deg(theta_rad)) 90; % 按实际几何调整这里有几个细节决定成败。第一拟合阶数二次还是三次液滴轮廓在小范围内曲率变化不大二次多项式够用如果液滴偏大、轮廓明显不对称用三次多项式更贴合。但阶数不是越高越好四阶以上对噪声敏感拟合曲线会故意穿过每个点的波动导致接触点处的切线抖得很离谱。第二拟合点的范围。不是把整条左轮廓都拿去拟合而是取接触点附近的一段——经典做法是取接触点往上约 20% 到 30% 的轮廓高度。范围选太长远处的轮廓形状会影响切线方向选太短拟合不稳定。第三左右接触角要分开算。液滴不可能绝对对称左右接触角有差异是正常的。左右各用各的接触点和轮廓段拟合分别计算。这也是判断程序算得准不准的一个线索如果左右角相差超过 10 度大概率是基线定歪了或者液滴滴歪了。4.4 方案三Young-Laplace 数值拟合——精度天花板如果对精度要求极高比如发论文、做标准测试那就得用 Young-Laplace 方程拟合了。这个方程描述的是液滴外形在表面张力和重力共同作用下的平衡形状是一个二阶常微分方程。思路是把测得的液滴轮廓与方程数值解做非线性最小二乘拟合直接拟合出接触角。MATLAB 实现通常用lsqcurvefit需要 Optimization Toolbox或fitnlm需要 Statistics Toolbox。核心难点在于初值选择和方程参数化拟合变量通常包括液滴底部半径、接触角、表面张力与密度之比毛细管长度等。由于问题高度非线性初值给得不好拟合直接发散。Young-Laplace 拟合对液滴大小没有限制对大液滴特别准因为它的物理模型完整考虑了重力对轮廓的压扁作用。但代价是代码量大、计算慢、参数敏感。我的建议是工程应用和常规实验用方案二多项式拟合已经足够了系数调好之后一致性很好只有在需要追溯真值的场景才上 Young-Laplace。5. 跑完程序之后结果验证与工程化改造5.1 怎么判断程序算得准不准程序跑出数字不等于程序算对了。我拿到一个新接触角程序测试靠谱程度只做三件事。第一把计算出的切线、基线和接触点画在原图上。这一步特别直观看那条切线是否真的贴合液滴轮廓的走向基线是否跟固体表面重合接触点是否在固液气三相的交界处。代码里加几行figure; imshow(I); hold on; plot(B_pts(:,2), B_pts(:,1), g-, LineWidth, 2); plot([x_base1 x_base2], [y_base1 y_base2], r--, LineWidth, 2); plot(x_contact_left, y_contact_left, bo, MarkerSize, 8); % 画切线沿斜率方向画一小段线段只要看到的图跟物理直觉一致那程序基本没问题。如果切线斜得离谱或者切在轮廓外面赶紧回头检查拟合点范围。第二拿 ImageJ 手工量角对比。ImageJ 的接触角插件如 DropSnake 或 LB-ADSA是很多实验室的参考工具拿同一张图在 ImageJ 里量一遍如果 MATLAB 算出来的结果跟 ImageJ 差在 2 度以内说明程序可信。差 3 度以上就说明算法里某个环节的系统性偏差没消除。第三重复性测试。同一张图跑十次结果应该完全一致。批量处理时关注统计量而不是单张值异常值比如角度小于 5 度或大于 175 度标记出来人工复核而不是粗暴删掉。5.2 批量处理从单张到几百张图片实验里往往是一次拍几十上百张图不可能一张张手动跑。把单图脚本改造成批量版本结构很简单files dir(fullfile(droplet_images, *.png)); results table(); for k 1:length(files) img imread(fullfile(files(k).folder, files(k).name)); [theta_left, theta_right] compute_contact_angle(img); results [results; table(files(k).name, theta_left, theta_right, ... VariableNames, {Filename, ThetaLeft, ThetaRight})]; end writetable(results, contact_angle_results.xlsx);批量处理里最容易翻车的是基线拟合。单张图手动确认没问题但换一张光线稍暗的图液滴底部拟合点范围选得不对基线就歪了。我的经验是给批量脚本加一个结果筛选条件左右角差超过 15 度或者角度超出合理范围的样本自动保存原图和处理后的轮廓图到另一个文件夹等全部跑完再人工检查这些异常样本。另外批量处理前建议先把所有图片统一尺寸。不同相机、不同分辨率拍出来的图bwareaopen的阈值和形态学结构元素大小都应该按比例调整。如果图片太大比如 4000×3000可以先imresize到宽 1000 像素左右接触角本质是角度量宽高比不变的情况下缩放不会改变角度值但处理速度能快十几倍。5.3 实操中踩过的坑与几个小提醒最后说几个我在调试这类程序时踩过的坑都是短视频教程不会提的细节。一个是镜头的水平问题。相机必须正对液滴镜头光轴与固体表面平行。如果相机有一点点俯仰拍出来的液滴会变形左右接触角一个偏大一个偏小。拍摄前用水平尺校准比事后用算法修正靠谱得多。另一个是液滴体积的选取问题。用微量注射器滴液滴不同实验之间尽量保持体积一致。液滴太大重力影响强轮廓偏离球冠液滴太小蒸发速度快拍摄时液滴可能在持续收缩。控制变量这件事程序帮不了你只能靠实验设计。还有一个出乎意料的小坑液滴内部的反光。如果光源位置不好液滴表面会出现高光斑块图像处理时这个反光区域的边缘可能被识别成第二层轮廓。形态学里imfill能填掉大部分反光造成的空洞但如果是很大面积的高光建议调整光源位置或者加扩散片从源头上解决。我自己在实际使用中养成的习惯是所有接触角算法都不做黑盒使用。不管程序是从网上下载的还是自己写的第一次跑通后一定花半小时把每一步的中间结果画出来看一遍——灰度图、二值图、边缘图、轮廓图、最终切线图。这五张图全部确认没问题批处理才有意义。接触角测量这件事算法的上限取决于图像质量程序只是把图像里的信息尽可能准确地提取出来。把样本拍好、把中间过程盯住比到处找更高级的拟合算法要管用得多。本文还有配套的精品资源点击获取
简介本资源是一套面向材料科学、表面物理及实验数据分析领域的MATLAB工具包聚焦接触角测量与固体表面能反演这一核心问题适用于高校研究生、科研人员及工程技术人员开展润湿性表征与界面性质研究。压缩包共含10个.m文件总大小仅6KB以主程序Contactangle.m为核心辅以多个fpush/fpop系列函数构成完整的图像边缘提取、液滴轮廓拟合与接触角计算流程支持基于Young-Laplace方程及Wenzel模型的表面能估算。已有576人学习下载体现了其在实验数据处理中的实用价值。用户可直接调用脚本完成从原始图像读取、预处理、接触点识别到角度输出与表面能推算的全流程分析代码结构清晰、模块分工明确且内置基础错误处理机制便于二次开发与适配不同实验条件下的液滴图像分析需求。 做材料表面分析的人几乎都遇到过这类场景从某个渠道下载了一个压缩包文件名写着contact angle.rar里面躺着一个contact_angle.m附带一两张液滴照片。网络热词里也总能看到matlab、contact angle这对组合说明大家都在找能直接跑通的接触角计算程序。但根据我的经验这类程序十有八九第一次运行不是报错就是结果离谱角度算出来 90 度上下浮动连自己都怀疑是不是量角器都比它准。这篇文章不是给你一份现成代码让你复制粘贴就完事而是把这套程序从能运行变成算得准的全过程拆开讲清楚文件里到底有什么、图像预处理哪些步骤省不得、接触角的几何定义是怎么落到atan和polyfit这些函数上的、以及跑完程序之后你凭什么相信那个数字。内容面向正在做接触角测量实验、需要批量处理液滴图片的工程师和科研人员你不需要精通图像处理但看完应该能独立把一套 MATLAB 程序改造得服服帖帖。1. 解压后的第一件事先搞清这个程序到底在算什么1.1 接触角测量到底测的是什么接触角contact angle指的是在固、液、气三相交界处从固液界面经过液滴内部到气液界面切线所夹的角度。这个角度直接反映固体表面的润湿性小于 90 度说明表面亲液大于 90 度说明疏液。测量方式看着简单——把一滴液体滴在固体表面从侧面拍张照片量一下就行了。但真到程序实现时那张照片里的信息量远比想象中多液滴轮廓是一条曲线固体表面是一条直线接触角是两者交点的局部几何关系而不是随便量一个整体角度。很多网上下载的 MATLAB 程序核心逻辑其实就三步读图、找轮廓、算切线。只不过这三步每一步都有不少隐藏细节。contact_angle.m这个脚本通常不是一个人写的它可能合并了多种方法甚至留有原作者实验条件下的硬编码参数。你直接跑很可能算的是某个特定背光条件下、某个特定固体表面的角度换一张图就失效。1.2 文件构成与运行环境检查解压contact angle.rar之后典型的文件结构是这样的contact_angle.m主脚本compute_contact_angle.m子函数负责从边缘点坐标计算角度droplet_sample.png或.bmp示例液滴图片可能还有一个README.txt或说明文档运行前先检查环境这步能省掉一半的报错。contact_angle.m普遍依赖 Image Processing Toolbox用到imread、rgb2gray、edge、bwareaopen等函数和 Curve Fitting Toolbox用到polyfit、polyval。在命令行敲ver查看工具箱列表缺了就去装否则一运行直接弹Undefined function edge这种错误很多人卡在这第一关。另一个容易忽略的是路径。把整个文件夹放在纯英文路径下运行别放桌面带中文的文件夹里。MATLAB 对中文路径的支持时好时坏imread偶尔抽风读不出文件排查起来非常浪费时间。我一般还会在脚本开头加上clc; clear; close all;避免上次运行的变量残留在工作区里干扰结果。提示如果解压后没有示例图片先别急着跑程序。找一张自己拍的液滴照片背景尽量简单、液滴轮廓清晰否则你根本分不清是代码问题还是图片问题。2. 图像预处理这步决定成败灰度化、降噪、背景分离2.1 为什么要先转灰度再做边缘检测相机拍出来的液滴照片是彩色图每个像素有 R、G、B 三个通道。直接对彩色图做边缘检测不是不行但三个通道分别响应边缘会变得极不均匀——红光下液滴边缘清楚蓝光下可能糊成一片。标准做法是先转成灰度图I imread(droplet_sample.png); if size(I, 3) 3 I rgb2gray(I); end这一步把三维矩阵压缩成二维之后的阈值分割、边缘检测全部基于亮度信息计算量也小得多。实测下来一张 1200 万像素的照片灰度化后处理速度能提升三倍左右。2.2 滤波窗口选多大才不会弄丢边缘转灰度之后紧接着是降噪。液滴照片的噪声来源通常有两类一是相机传感器的高斯噪声二是背景里的微小颗粒或灰尘产生的椒盐噪声。很多人一上来就用高斯滤波做平滑结果液滴边缘被抹得圆乎乎接触角算出来系统性地偏大或偏小。我的习惯是先做中值滤波medfilt2窗口取 3×3最多 5×5。中值滤波对椒盐噪声特别有效而且能较好地保留边缘的陡峭程度。高斯滤波imgaussfilt不是不能用但 sigma 超过 1.5 就容易把液滴顶部的小弧度细节吃掉接触角测量对轮廓局部的切线方向极其敏感边缘偏移哪怕一两个像素角度就能差出好几度。I_filt medfilt2(I, [3 3]);这里有个容易被忽略的点滤波窗口越大边缘越钝。做了 5×5 中值滤波之后液滴边缘线会被轻微腐蚀算出来的接触角往往比真实值偏大。所以滤波的原则是能去掉噪声就行窗口尽量小。如果原图噪声不大直接跳过滤波也行。2.3 自动阈值和 ROI 到底该信谁接下来是二值化把液滴从背景里分离出来。最简单的办法是imbinarize配合graythresh自动阈值BW imbinarize(I_filt, graythresh(I_filt));graythresh用的是 Otsu 方法它会自动找一个阈值让前景和背景的类间方差最大。这个方法对背景均匀、液滴和背景对比明显的照片效果很好。问题在于实际拍摄环境很难这么理想。如果光线从侧面打过来液滴表面会有镜面反射反光区域亮度接近背景Otsu 会把那一块误判成背景导致液滴轮廓缺个口子。这种情况下我不会完全信任自动阈值而是先用imcrop把图像裁剪到液滴附近的区域。裁剪有两个好处一是去掉画面边缘无关的高亮区域让 Otsu 的阈值更准二是大幅减少后续边缘检测的计算量。裁剪时不要裁得太紧液滴下方要留出足够的固体表面区域方便后面确定基线。裁剪之后再次imbinarize或者edge成功率会明显提升。如果背景依然不均匀比如有渐变阴影可以试试imbinarize(I_filt, adaptive)它按局部区域计算阈值对光照不均的图片比全局阈值稳得多。代价是速度慢而且参数ForegroundPolarity需要根据液滴是亮还是暗来调整。3. 从二值图到液滴轮廓边缘检测与轮廓筛选3.1 Canny 阈值怎么调才靠谱二值图拿到手下一步是提取液滴边缘。MATLAB 的edge函数支持多种算子算接触角我强烈推荐 Canny。Sobel 算子边缘粗、对噪声敏感算切线方向容易被局部毛刺带偏Canny 做了非极大值抑制和双阈值边缘细而且连续性好。edges edge(BW, Canny, [0.1 0.2]);双阈值的理解很简单低阈值决定边缘检测的灵敏度高阈值决定保留哪些强边缘。[0.1 0.2]这个组合是大多数液滴图片的常用起点。如果检测出来的边缘断裂严重说明高阈值太高降到[0.05 0.15]如果边缘旁边跟着一堆杂线说明低阈值太低往上提到[0.15 0.3]。有个经验技巧在调阈值之前先把BW显示出来肉眼确认液滴在二值图里是不是完整的一块。如果二值图本身就有缺口边缘检测调参怎么调都白搭。这时候应该回去调阈值分割或者用形态学操作把缺口补上。3.2 形态学处理与轮廓去杂二值图里的液滴经常不是完美的一块反光点会在液滴内部形成空洞灰尘会在液滴旁边形成孤立的小白点。直接用bwboundaries提取轮廓你会得到一串乱七八糟的边界只有其中一条才是液滴轮廓。形态学处理能解决大部分问题BW_filled imfill(BW, holes); % 填充液滴内部空洞 BW_clean bwareaopen(BW_filled, 500); % 去掉小于500像素的孤立区域 BW_closed imclose(BW_clean, strel(disk, 3)); % 连接断开的边缘bwareaopen里的 500 是经验值根据你的图片分辨率调整。如果液滴占了画面的大部分这个值可以设成 2000 甚至更大如果液滴本身就很小设 500 可能把液滴主体给删了。判断标准很简单处理后画面里应该只剩下一个最大的连通区域。轮廓提取用bwboundariesB bwboundaries(BW_closed, noholes);返回的B是一个元胞数组每个元素是一圈边界的坐标点N×2 的矩阵第一列是行号 y第二列是列号 x。液滴轮廓通常是其中最长的那个边界。找出来后用plot叠在原图上看一眼确认轮廓是不是完整贴合液滴的外形。3.3 基线固体表面线的确定方法接触角的定义里基线就是固体表面线。如果液滴是滴在一块平整的固体表面上那么照片里固体表面就是一条接近水平的直线。这条线的位置直接决定了接触点的坐标偏差哪怕半度都会让最终角度失真。确定基线有三种方法我按推荐程度排个序方法一霍夫变换检测直线。用houghlines在边缘图里找线段筛选角度接近 0 度水平的那条。这个方法自动但容易误检液滴底部的弧线在某些角度下会被误判成水平线。方法二手动指定两个点。用ginput在图像上点击固体表面的左右两个位置连成直线。这个最稳但不适合批量处理。方法三拟合液滴轮廓的最低点区间。这个方法我用得最多——液滴与固体表面接触的那一段轮廓与基线几乎重合。取轮廓点里 y 坐标最大的那部分点大概是底部 5% 到 10% 的点用polyfit(x, y, 1)做线性拟合拟合出的直线就是基线。idx find(B_pts(:,1) prctile(B_pts(:,1), 90)); p_base polyfit(B_pts(idx,2), B_pts(idx,1), 1);方法三的好处是完全自动且物理意义明确液滴底部被固体表面压平的那一段天然就是基线的一部分。唯一的风险是液滴太小、底部弧线太弯这时拟合出来的线会有点倾斜需要人工检查。基线确定后轮廓与基线的交点就是左右两个接触点。这两个点的 x 坐标一个是最小值、一个是最大值取轮廓点里离基线最近且 x 坐标分别最小和最大的点即可。4. 接触角计算的几何本质与三种实现方案4.1 图像坐标系下的角度陷阱接触角的核心计算其实就是一个初中几何问题曲线上某一点的切线斜率决定了该点切线的角度。但 MATLAB 图像坐标系有一个大坑——原点在左上角x 轴向右y 轴向下。这跟我们熟悉的笛卡尔坐标系正好相反导致atan算出来的角度方向全部翻转。举个例子液滴左侧轮廓从接触点向上走x 坐标逐渐增大接触点在左往右到液滴顶部y 坐标逐渐减小向上在图像坐标系里斜率 dx/dy 是负的。用atan算出来的角度是负值但是真实的物理角度是正的。这就是为什么很多新手程序算出来接触角一直在 0 到 10 度之间徘徊以为固体表面超亲液其实是坐标系符号没处理。正确处理思路先明确你要的角度是接触点处液滴轮廓切线与水平基线之间的夹角然后画个图理清几何关系。用atan2而不是atan它能正确处理四个象限。算出切线向量后再根据左右侧分别调整角度定义。4.2 方案一高宽比法θ/2 法——最快但有限制高宽比法也叫 θ/2 法是接触角计算里最朴素的方法。它假设液滴是一个球冠理想情况下液滴足够小重力影响可忽略此时接触角只跟液滴高度 h 和接触宽度 w 有关theta 2 * atan(2 * h / w)h 是液滴顶点到基线的垂直距离w 是左右接触点之间的水平距离。用 MATLAB 实现就是先找到轮廓的最高点再算高度和宽度。这个方法的优点是代码量十几行计算量几乎为零而且对噪声不敏感。缺点是假设条件太苛刻液滴必须足够小一般直径小于 2 mm才能忽略重力对形状的拉拽而且液滴必须左右对称。实际实验里液滴一旦稍微大一点轮廓就被重力压扁不再是球冠高宽比法算出来的角度偏小。我一般把 θ/2 法当快速预估工具用批量处理几十张图先跑一遍看看数据分布是否合理再对关键样本用更精确的方法复核。4.3 方案二多项式拟合与切线法——工程界的主力方案这是目前工程和科研里最常用的接触角算法也是contact_angle.m这类程序里最可能采用的方法。核心思路是把液滴左右两侧的轮廓分别提取出来用多项式拟合轮廓曲线然后在接触点处求导得到切线斜率再换算成角度。以左侧轮廓为例代码逻辑是这样的% 提取左轮廓点x 最小的那段轮廓 left_mask (B_pts(:,2) min(B_pts(:,2))) | ... (B_pts(:,2) prctile(B_pts(:,2), 5)); % 这里只是示意实际要按连通区间取左边缘点 x_left B_pts(left_mask, 2); y_left B_pts(left_mask, 1); % 二次多项式拟合 y a*x^2 b*x c p polyfit(x_left, y_left, 2); % 在接触点处求导dy/dx 2*a*x b x_contact min(x_left); slope polyval(polyder(p), x_contact); % 切线角度弧度转角度 theta_rad atan(slope); % 图像坐标系翻转修正左接触角 theta_left 180 - abs(rad2deg(theta_rad)) 90; % 按实际几何调整这里有几个细节决定成败。第一拟合阶数二次还是三次液滴轮廓在小范围内曲率变化不大二次多项式够用如果液滴偏大、轮廓明显不对称用三次多项式更贴合。但阶数不是越高越好四阶以上对噪声敏感拟合曲线会故意穿过每个点的波动导致接触点处的切线抖得很离谱。第二拟合点的范围。不是把整条左轮廓都拿去拟合而是取接触点附近的一段——经典做法是取接触点往上约 20% 到 30% 的轮廓高度。范围选太长远处的轮廓形状会影响切线方向选太短拟合不稳定。第三左右接触角要分开算。液滴不可能绝对对称左右接触角有差异是正常的。左右各用各的接触点和轮廓段拟合分别计算。这也是判断程序算得准不准的一个线索如果左右角相差超过 10 度大概率是基线定歪了或者液滴滴歪了。4.4 方案三Young-Laplace 数值拟合——精度天花板如果对精度要求极高比如发论文、做标准测试那就得用 Young-Laplace 方程拟合了。这个方程描述的是液滴外形在表面张力和重力共同作用下的平衡形状是一个二阶常微分方程。思路是把测得的液滴轮廓与方程数值解做非线性最小二乘拟合直接拟合出接触角。MATLAB 实现通常用lsqcurvefit需要 Optimization Toolbox或fitnlm需要 Statistics Toolbox。核心难点在于初值选择和方程参数化拟合变量通常包括液滴底部半径、接触角、表面张力与密度之比毛细管长度等。由于问题高度非线性初值给得不好拟合直接发散。Young-Laplace 拟合对液滴大小没有限制对大液滴特别准因为它的物理模型完整考虑了重力对轮廓的压扁作用。但代价是代码量大、计算慢、参数敏感。我的建议是工程应用和常规实验用方案二多项式拟合已经足够了系数调好之后一致性很好只有在需要追溯真值的场景才上 Young-Laplace。5. 跑完程序之后结果验证与工程化改造5.1 怎么判断程序算得准不准程序跑出数字不等于程序算对了。我拿到一个新接触角程序测试靠谱程度只做三件事。第一把计算出的切线、基线和接触点画在原图上。这一步特别直观看那条切线是否真的贴合液滴轮廓的走向基线是否跟固体表面重合接触点是否在固液气三相的交界处。代码里加几行figure; imshow(I); hold on; plot(B_pts(:,2), B_pts(:,1), g-, LineWidth, 2); plot([x_base1 x_base2], [y_base1 y_base2], r--, LineWidth, 2); plot(x_contact_left, y_contact_left, bo, MarkerSize, 8); % 画切线沿斜率方向画一小段线段只要看到的图跟物理直觉一致那程序基本没问题。如果切线斜得离谱或者切在轮廓外面赶紧回头检查拟合点范围。第二拿 ImageJ 手工量角对比。ImageJ 的接触角插件如 DropSnake 或 LB-ADSA是很多实验室的参考工具拿同一张图在 ImageJ 里量一遍如果 MATLAB 算出来的结果跟 ImageJ 差在 2 度以内说明程序可信。差 3 度以上就说明算法里某个环节的系统性偏差没消除。第三重复性测试。同一张图跑十次结果应该完全一致。批量处理时关注统计量而不是单张值异常值比如角度小于 5 度或大于 175 度标记出来人工复核而不是粗暴删掉。5.2 批量处理从单张到几百张图片实验里往往是一次拍几十上百张图不可能一张张手动跑。把单图脚本改造成批量版本结构很简单files dir(fullfile(droplet_images, *.png)); results table(); for k 1:length(files) img imread(fullfile(files(k).folder, files(k).name)); [theta_left, theta_right] compute_contact_angle(img); results [results; table(files(k).name, theta_left, theta_right, ... VariableNames, {Filename, ThetaLeft, ThetaRight})]; end writetable(results, contact_angle_results.xlsx);批量处理里最容易翻车的是基线拟合。单张图手动确认没问题但换一张光线稍暗的图液滴底部拟合点范围选得不对基线就歪了。我的经验是给批量脚本加一个结果筛选条件左右角差超过 15 度或者角度超出合理范围的样本自动保存原图和处理后的轮廓图到另一个文件夹等全部跑完再人工检查这些异常样本。另外批量处理前建议先把所有图片统一尺寸。不同相机、不同分辨率拍出来的图bwareaopen的阈值和形态学结构元素大小都应该按比例调整。如果图片太大比如 4000×3000可以先imresize到宽 1000 像素左右接触角本质是角度量宽高比不变的情况下缩放不会改变角度值但处理速度能快十几倍。5.3 实操中踩过的坑与几个小提醒最后说几个我在调试这类程序时踩过的坑都是短视频教程不会提的细节。一个是镜头的水平问题。相机必须正对液滴镜头光轴与固体表面平行。如果相机有一点点俯仰拍出来的液滴会变形左右接触角一个偏大一个偏小。拍摄前用水平尺校准比事后用算法修正靠谱得多。另一个是液滴体积的选取问题。用微量注射器滴液滴不同实验之间尽量保持体积一致。液滴太大重力影响强轮廓偏离球冠液滴太小蒸发速度快拍摄时液滴可能在持续收缩。控制变量这件事程序帮不了你只能靠实验设计。还有一个出乎意料的小坑液滴内部的反光。如果光源位置不好液滴表面会出现高光斑块图像处理时这个反光区域的边缘可能被识别成第二层轮廓。形态学里imfill能填掉大部分反光造成的空洞但如果是很大面积的高光建议调整光源位置或者加扩散片从源头上解决。我自己在实际使用中养成的习惯是所有接触角算法都不做黑盒使用。不管程序是从网上下载的还是自己写的第一次跑通后一定花半小时把每一步的中间结果画出来看一遍——灰度图、二值图、边缘图、轮廓图、最终切线图。这五张图全部确认没问题批处理才有意义。接触角测量这件事算法的上限取决于图像质量程序只是把图像里的信息尽可能准确地提取出来。把样本拍好、把中间过程盯住比到处找更高级的拟合算法要管用得多。本文还有配套的精品资源点击获取
