[发明专利]一种基于显微图像的工业用精对苯二甲酸粒径分布估计方法

权利要求书

1.一种基于工业用精对苯二甲酸显微图像的新的粒径分布估计方法，其特 征在于包括如下步骤：

1)将体视显微镜成像系统总放大倍数设置为100倍，其中CCD相机的放大 倍数为25倍，物镜镜头的放大倍数为2倍，调节升降调焦旋钮放大倍数为2倍；

2)将测微尺置于载物台上，调节升降调焦旋钮，直到观察到完全清晰的测 微尺图像，利用计算机获取测微尺图像中每一像素所代表的真实长度；

3)将工业用精对苯二甲酸粉末放入样品皿中，置于振动装置中处理1-2分钟 后转移至载物台的黑色背光板上；

4)调节升降调焦旋钮直到完全清晰地观察到工业用精对苯二甲酸颗粒，保 存图像至计算机；

5)将工业用精对苯二甲酸图像转换为灰度图像，进行二值化处理和滤波处 理，提取图像中工业用精对苯二甲酸颗粒的数目、真实面积和最大粒径；

6)将工业用精对苯二甲酸颗粒的真实面积S_actual(i)，i＝1，2，...，n，转换为各工业 用精对苯二甲酸颗粒的当量直径和体积，计算公式如下：

Sactual(i)=π4Dactual2(i)⇒Dactual(i)=4π×Sactual(i)]]>

Vactual(i)=π6Dactual3(i),]]>i＝1，2，...，n

其中D_actual(i)和V_actual(i)分别为第i个工业用精对苯二甲酸颗粒的当量直径和体积；

7)重复步骤3)至步骤6)，得到多幅工业用精对苯二甲酸图像的当量直径 和体积；

8)利用工业用精对苯二甲酸的当量直径获得工业用精对苯二甲酸的粒径分 布和体积分布；

9)利用标准正态分布、偏态分布和对数正态分布对工业用精对苯二甲酸粒 径分布和体积分布结果进行拟合，分别选取误差平方和最小的分布作为工业用 精对苯二甲酸的真实粒径分布和真实体积分布；

所述步骤2)：

(1)将测微尺放在白色背光板上，调节升降调焦旋钮，找到刻度线外的圆 之后，寻找刻度线并将刻度线移至视野中央；调节升降调焦旋钮直到能完全清 晰地观察到各刻度线，保存图像至计算机；

(2)将获得的图像转换为灰度图像，选择参考灰度阈值为50，对灰度图像 进行二值化处理，其中灰度值大于该阈值的为背景区域，否则为目标区域；

(3)利用Matlab的ordfilt2函数对图像进行滤波处理，参数选择为6行1列的 矩阵；

(4)对滤波后的二值图像进行行扫描和列扫描，找出各目标区域，其中刻 度线的标号为K1，K2，...，Kn，各区域代表边界点的坐标值保存在相应的矩阵中，K1 号和Kn号刻度线的中心坐标通过以下公式得到：

X(p)‾=1mΣi=1mX_coordinate(p,i),]]>Y(p)‾=1mΣi=1mY_coordinate(p,i),]]>p＝K1，…Kn

其中是p号刻度线的中心坐标，m是各区域的代表边界点的坐标对 数目，X_coordinate为刻度线的横坐标，Y_coordinate为刻度线的纵坐标；

(5)K1号刻度线和Kn号刻度线坐标中心之间的像素数通过以下公式得到：

Pixels=(X(K1)‾-X(Kn)‾)2+(Y(K1)‾-Y(Kn)‾)2,]]>

每一像素代表的真实长度通过以下公式得到：

1pixel=1Pixelsmm=1000Pixelsμm;]]>

所述步骤5)：

(1)将工业用精对苯二甲酸显微图像转换为灰度图像，选择灰度值阈值为 50，将灰度图像转换为相应的二值图像，其中灰度值大于该阈值的为目标区域， 否则为背景区域；

(2)从二值图像矩阵的第1行、第1列开始每Res行、Res列进行扫描，找出 各工业用精对苯二甲酸颗粒的边界点坐标和像素面积S_pixel(i)，并分别保存；

(3)各工业用精对苯二甲酸颗粒的真实面积S_actual(i)通过以下公式得到： i＝1，2，...，n，其中，S_actual(i)是第i个工业用精对苯二甲酸 颗粒的真实面积，单位为μm²。