[發明專利]一種基于分水嶺分割的顯微圖像顆粒粒度測量方法

有效

• 申請號: 201410311363.7
• 申請日: 2014-07-02
• 公開（公告）號: CN104075965A
• 公開（公告）日: 2014-10-01
• 發明（設計）人: 李雨洋;王贏元
• 申請（專利權）人: 北京機械設備研究所
• 主分類號: G01N15/02
• 分類號: G01N15/02
• 代理公司: 中國航天科工集團公司專利中心 11024
• 代理人: 岳潔菱

1.一種基于分水嶺分割的顯微圖像顆粒粒度測量方法，其特征在于本方法的具體步驟為：?

第一步構建基于分水嶺分割的顯微圖像顆粒粒度測量系統?

基于分水嶺分割的顯微圖像顆粒粒度測量系統，包括：光學顯微鏡、CCD攝像機、圖像采集卡、圖像預處理模塊、圖像分割模塊和顆粒參數測量模塊；?

圖像預處理模塊的功能為：對獲取的原始顯微圖像進行預處理，去除原始顯微圖像的噪聲，強化顆粒的輪廓、邊緣；?

圖像分割模塊的功能為：計算預處理后顯微圖像的梯度圖像，通過由梯度圖像構建的數據矩陣計算內部標記符集合和外部標記符集合，獲得分割后顯微圖像；?

顆粒參數測量模塊的功能為：計算分割后顯微圖像中顆粒的數目、周長、面積、等效粒徑和圓形度；?

圖像采集卡與CDD攝像機通過RS-422接口連接，圖像采集卡與測量計算機通過PCI總線連接；圖像預處理模塊、圖像分割模塊和顆粒參數測量模塊分別置于測量計算機中；?

第二步原始顯微圖像獲取?

將顆粒樣本放置于測量平臺上，經光學顯微鏡放大后由CCD攝像機成像；圖像采集卡對CCD攝像機輸出的圖像視頻信號進行采集，經A/D轉換后通過PCI總線將原始顯微圖像存儲于測量計算機中；?

第三步圖像預處理模塊對原始顯微圖像進行平滑處理?

圖像預處理模塊接收存儲于測量計算機中的原始顯微圖像，采用滑動窗口為3*3的中值濾波器去除原始顯微圖像中的干擾和噪聲，保留原始顯微圖像中顆粒的輪廓、邊緣信息；?

第四步圖像預處理模塊對平滑后顯微圖像進行二值化處理?

圖像預處理模塊利用平滑后顯微圖像水平方向和垂直方向像素的灰度差計算大小為M×N的顯微圖像全局閾值：?

其中：?

e(x，y)＝max(|e_x|，|e_y|)，?

e_x＝f(x-1，y)-f(x+1，y)，?

e_y＝f(x，y-1)-f(x，y+1)；?

x表示顯微圖像水平方向的坐標，y表示顯微圖像豎直方向的坐標；?

f(x，y)表示顯微圖像像素(x，y)的灰度值；?

f(x-1，y)表示顯微圖像像素(x，y)的左方像素(x-1，y)的灰度值；?

f(x+1，y)表示顯微圖像像素(x，y)的右方像素(x+1，y)的灰度值；?

f(x，y-1)表示顯微圖像像素(x，y)的下方像素(x，y-1)的灰度值；?

f(x，y+1)表示顯微圖像像素(x，y)的上方像素(x，y+1)的灰度值；?

|e_x|表示顯微圖像像素(x，y)左方像素(x-1，y)與右方像素(x+1，y)灰度差的絕對值；?

|e_y|表示顯微圖像像素(x，y)下方像素(x，y-1)與上方像素(x，y+1)的灰度差的絕對值；?

由公式(2)得到二值化顯微圖像：?

第五步圖像分割模塊獲得二值化顯微圖像的梯度圖像?

圖像分割模塊利用3*3的Sobel掩膜計算二值化顯微圖像的梯度圖像；顯微圖像中像素(x,y)的梯度為：?

其中G_x為顯微圖像水平方向的一階導數，G_y為顯微圖像豎直方向的一階導數；?

第六步圖像分割模塊對梯度圖像進行分割?

圖像分割模塊把梯度圖像存儲為數據矩陣，數據矩陣中的每個數值為梯度圖像中對應像素的灰度值，計算梯度圖像的極大值，將極大值點作為分水嶺的內部標記符集合；對內部標記符集合進行距離變換后再進行一次分水嶺變換得到分水嶺脊線，將分水嶺脊線作為分水嶺的外部標記符集合；采用最小強制方式對梯度圖像進行修改，由此得到分割后顯微圖像；?

第七步顆粒參數測量模塊對顆粒數目、大小、形狀進行測量?

顆粒參數測量模塊采用線掃描方式對分割后顯微圖像中各目標區域標號，最大的區域標號即為顯微圖像中顆粒的數目N；?

顆粒參數測量模塊統計目標區域中顆粒邊緣的像素個數，即為此目標區域中顆粒的周長L；?

目標區域中顆粒的面積為顆粒邊緣與邊緣內包含的像素個數之和，此目標區域的顆粒面積為：?

由此計算出顆粒的等效粒徑為：?

圓形度為：?

經過上述各步驟實現了基于分水嶺分割的顯微圖像顆粒粒度測量。?