接觸角的計算 尋找氣液固三相的交點

2021-08-20 06:20:28 字數 4072 閱讀 2819

接觸角(contact angle)指在氣、液、固三相交點處所作的氣-液介面的切線穿過液體與固-液交界線之間的夾角(θ),也是潤濕程度的量度,主要測量方法為外形影象分析方法及稱重法。————互動百科

%%%%%%%%%%%%接觸角交點的自動計算%%%%%%%%%%%%%%

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 = [stats.area];

img_change = bwareaopen(img_bin,max(areas));%去掉小於最大連通域面積的小連通域

figure

imshow(img_change)

title('img_change')

% imshow(label2rgb(l, @jet, [.5 .5 .5]))%顯示影象

img_edge = bwperim(img_change,8) ; %查詢二值影象的邊緣

figure;

handle = imshow(img_edge);

title('img_edge');

[m,n] = size(img_change);

img_coordinate = zeros(m,n);

p = 0;

forj = 100:n-100

fori = m:-1:1

if img_change(i,j) == 1

img_coordinate(i,j) = 1;

p = p+1;

coor(p,1) = i;

coor(p,2) = j;

break;

endendend

% figure;imshow(img_coordinate),title('img_coordinate');

% x,y座標分離

coor_x = coor(:,2);

coor_y = coor(:,1);

fori = 2:size(coor,1)

if coor_x(i)~= coor_x(i-1) %防止斜率出現無窮大的情況,保證程式可以正常執行

k(i-1) = (coor_y(i)-coor_y(i-1))/(coor_x(i)-coor_x(i-1)); %相鄰兩點之間的斜率

endend

forj = 1:size(k,2)

if k(j)<8&&k(j)>-8

%小於或者大於指定的斜率k認為是固體平面的斜率,並將其置0

k(j) = 0;

endend[kv,index] = sort(k,'descend'); %kv:相連點的斜率,index:kv值對應的索引

catastrophe_a = coor_x(index(1));

catastrophe_b = coor_x(index(end));

if catastrophe_a>catastrophe_b

left_p = index(end); %左突變點

right_p = index(1); %右突變點

else

left_p = index(1); %左突變點

right_p = index(end);%右突變點

end%%%%%%%%%%%%%%%直線座標與圓弧座標點的分離%%%%%%%%%%%%%

line_num = 0;%計數器

circle_num = 0;%計數器

fori = 1:size(coor,1)

ifi2||i>right_p+2

line_num = line_num+1;

line(line_num,1) = coor_x(i); %line第一列是x,第二列是y

line(line_num,2) = coor_y(i);

else

circle_num = circle_num+1;

circle(circle_num,1) = coor_x(i); %circle第一列是x,第二列是y

circle(circle_num,2) = coor_y(i);

endendhold on

[b,l] = bwboundaries(img_coordinate,'noholes');%尋找邊緣,不包括孔

imshow(label2rgb(l, @jet, [.5

.5.5]))%顯示影象

for k = 1:length(b) %length(b)=4 即:b中有4個封閉區域

boundary = b;

handle = plot(boundary(:,2),boundary(:,1), 'g', 'linewidth',0.00001); %boundary(:,2)代表索引矩陣第二列

end%整個迴圈表示的是描邊

title('邊緣提取後的影象');

str = ('交點');

msize = 8;

handle = plot(coor_x(left_p),coor_y(left_p),'o','linewidth',2,'color','r');

text(coor_x(left_p),coor_y(left_p)+50,str,'fontsize',msize);

handle = plot(coor_x(right_p+1),coor_y(right_p+1),'o','linewidth',2,'color','r');

text(coor_x(right_p+1),coor_y(right_p+1)+50,str,'fontsize',msize);

msize = 8;

handle = plot(coor_x(left_p),coor_y(left_p),'o','linewidth',2,'color','r');

text(coor_x(left_p),coor_y(left_p)+50,str,'fontsize',msize);

handle = plot(coor_x( right_p + 1 ),coor_y( right_p + 1 ),'o','linewidth',2,'color','r');

text(coor_x(right_p+1),coor_y(right_p+1)+50,str,'fontsize',msize);`

變成

msize = 8;

handle = plot(coor_x(left_p),coor_y(left_p),'o','linewidth',2,'color','r');

text(coor_x(left_p),coor_y(left_p)+50,str,'fontsize',msize);

handle = plot(coor_x( right_p ),coor_y( right_p ),'o','linewidth',2,'color','r');

text(coor_x(right_p+1),coor_y(right_p+1)+50,str,'fontsize',msize);

得到的交點圖

如何形成歐姆接觸 肖特基勢壘的計算金屬源漏接觸

金屬 半導體接觸 金屬與半導體的接觸分為兩類 第一類是歐姆接觸,指的是其接觸面的電阻值遠小於半導體本身的電阻,在器件工作時,大部分的電壓降在有效區域 active region 而不在接觸面,且不會使半導體內部的平衡載流子濃度發生顯著的改變。形成歐姆接觸的條件有兩種 金屬與半導體間有低的勢壘高度或者...