HOWTO: 如何利用Avizo或Amira計算孔隙率（Porosity）

在作三維可視化數據處理過程當中，咱們常常要提取的一個基本信息就是孔隙率。在今天的文章中咱們要分享兩個信息，一個是如何利用Avizo或Amira進行孔隙率計算；另外是關於Avizo 8.0中孔隙率計算異常的問題應對方法。

Avizo或Amira中如何計算孔隙率

孔隙率（Porosity），簡單理解就是材料中孔隙體積與材料在天然狀態下整體積的百分比。在Avizo和Amira中是利用Volume Fraction計算模塊來提取孔隙率信息，基本步驟以下：

1. 分割獲取孔隙（Foam示例中我只是簡單利用交互閾值分割方式獲取了孔隙，以下圖）：

　　通常這裏我會將得到的孔隙重命名爲Pores。

2. 獲取材料總體分割（相似步驟1，選擇所有閾值範圍，就能夠得到整個材料的分割，不過注意，這裏材料正好是正方形，若是是不規則材料，要考慮如何獲取材料總體形態）

      

　　通常這裏我會將分割得到的材料重命名爲Mark。

3. 利用Volume Fraction進行孔隙率計算。

　　對Pores數據應用Volume Fraction，默認第一個參數就是Pores，第二個Input Image Mask咱們選擇Mark，第三個參數Interpretation選擇3D（若是選擇XY，咱們能夠得到面孔率），應用以後獲取的結果點擊Show，在Tables中展現出得到的孔隙率，以下圖：

        

Avizo 8.0計算孔隙率異常的問題

　　使用Avizo 8.0的用戶要注意，按照前面的方法計算孔隙率極可能是錯誤的，Avizo廠商在8.1的Release Notes中描述了這個問題「 23292 Volume Fraction (Avizo Fire)

Volume Fraction could give incorrect values in some case when used with a mask.」

　　也就是說這個問題在8.1以後版本已經修正，因此這也提醒Avizo或Amira用戶，軟件及時更新很是重要。

　　那這裏還提供一個臨時方案，參考前邊完成孔隙和材料總體的分割，而後利用Material Statistics對材料進行統計計算，獲取的計算結果表中提取Volume列對應材料的體積數，而後相除，也能計算得到孔隙率。