癌細胞形態學分析

clc;

clear;
i=imread('cancer02.bmp');       % 注意必須保證二值影象中，細胞區域為白色區域或者畫素點值為「1」
[x,y]=size(i);
bw = bwperim(i,8);              % 檢測細胞的邊緣跟蹤，用於計算周長 
figure(3);
imshow(bw);
%檢測垂直方向連讀的周長畫素點%
p1=0;
p2=0;
ny=0;                           % 記錄垂直方向連續周長畫素點的個數
for i=1:x
for j=1:y
if (bw(i,j)>0)
p2=j;
if ((p2-p1)==1)    % 判斷是否為垂直方向連續的周長畫素點
ny=ny+1;
endp1=p2;
endend
end    
%檢測水平方向連讀的周長畫素點%
p1=0;
p2=0;
nx=0;                           % 記錄水平方向連續周長畫素點的個數
for j=1:y
for i=1:x
if (bw(i,j)>0)
p2=i;
if ((p2-p1)==1)    % 判斷是否為水平方向連續的周長畫素點
nx=nx+1;
endp1=p2;
endend
end 
sn=sum(sum(bw));            % 計算周長畫素點的總數
nd=sn-nx-ny;                % 計算奇數碼的鏈碼數目
h=max(sum(i));              % 計算細胞的高度 
w=max(sum(i'));             % 圖象i經矩陣轉置後，計算寬度
l=sqrt(2)*nd+nx+ny;         % 計算周長
%====四個形態特徵值計算===%
a=bwarea(i);                % 計算細胞的面積
c=4*pi*a/(l*l);             % 計算圓度
r=a/(h*w);                  % 計算矩形度
e=min(h,w)/max(h,w);        % 計算伸長度
%==設定相關閾值，識別癌細胞==%
ath1=10000; ath2=50000;
cth=0.5;      rth=0.5;   eth=0.8;
if ((a>=ath1)&&(a=cth)&&(r<=rth)&&(e>eth))
cancer_right=1                     % 結論為癌細胞
end    
else if (a>=ath2)
cancer_right=2                     % 結論為可疑小細胞癌細胞      
else
cancer_right=0                         % 結論為正常細胞 
endend
figure(1);
imshow(i);
figure(2);
imshow(bw);

形態學重建

在形態學梯度影象的基礎上，利用形態學的開閉重建運算對梯度影象進行重建，在保留重要區域倫敦的同時去除細節和雜訊。分水嶺變換存在過分割現象，原因在於檢測的區域性極值過多，造成極值過多的原因在於影象中的非規則灰度擾動和雜訊。對於好的分水嶺影象分割方法，不僅能消除過分割現象，而且應保證分割後的區域倫敦邊緣具...

形態學點滴

令a為z中的集合，如果a a1，a2 是a中的元素，則 a a 如果a不是a中的元素，則 a a matlab函式 ismember a,a 似乎還有個函式in與此類似 另外 判斷集合中的元素是否按公升序排列 issorted a 想想issorted a end 1 1 是什麼意思？刪除集合中的重...

形態學濾波

include pch.h include include include include using namespace std using namespace cv 形態學 mat g srcimage,g dstimage int g nelementshape morph rect int ...