接触角的计算----寻找气液固三相的交点

接触角的定义

接触角（contact angle）指在气、液、固三相交点处所作的气-液界面的切线穿过液体与固-液交界线之间的夹角（θ），也是润湿程度的量度，主要测量方法为外形图像分析方法及称重法。————互动百科

具体思路

• 以下图为例：
  

  
  
  

  1. 将原图转为灰度图在转为二值图
     

     
     
     

  2. 去掉连通域中较小的区域。（MATLAB中默认白色为前景区域）
     

     
     
     

  3. 寻找边缘，并找出曲率突变的点即为突变点。

     Note:黄色的线为液滴附着的水平面，一般斜率都不会太大，故在寻找突变点（即交点的时候）可以认为水平面的斜率在一个设置的范围内，在此范围内的点的斜率都置为0，剩下的斜率变化最大的点即为所找的两个突变点，这样利用突变点加上对蓝色弧线（液滴的边缘）的拟合即可求出接触角。

  
  ，
  

  代码实现

%%%%%%%%%%%%接触角交点的自动计算%%%%%%%%%%%%%%
clc
clear
img_orignal =imread('0517.jpg') ;
img_gray =rgb2gray(img_orignal);%变灰度
% figure;imshow(img_gray);title('灰度图像');
level=graythresh(img_gray);
img_bin =im2bw(img_gray,level) ; %灰度图变二值图
figure;imshow(img_bin),title('二值图像');
img_bin = bwmorph(img_bin,'clean',Inf);
% figure;imshow(img_bin);title('clean')
STATS = regionprops(img_bin,'area'); %获得各个连通域的面积
areas =