基于DBSCAN的二氧化硅融化過程中質心定位方法

王浩然 楊煉鑫 張劍書

摘要：二氧化硅的融化行為可以在很大程度上代表鐵尾礦的融化行為。通過CCD攝像機可以以非接觸式的方法采集到二氧化硅融化過程的時序圖像，利用圖像分割、邊緣檢測等一系列圖像分析算法從圖像中提取原始圖像數據的輪廓特征，并在此基礎上通過DBSCAN密度聚類算法提取二氧化硅的輪廓，以發現圖像數據中二氧化硅的質心的位置數據?；诙趸璧膶崟r質心位置數據可以進一步發現二氧化硅的融化速率，為提高鐵尾礦的利用率提供理論依據。

關鍵詞：圖像分割; 邊緣檢測; DBSCAN算法; 二氧化硅融化過程;二氧化硅質心位置

中圖分類號：TP18? ? ? ? 文獻標識碼：A

文章編號：1009-3044（2020）23-0001-03

Abstract：The melting behavior of silicon dioxide can largely represent the melting behavior of iron tailings. The time sequence image of silicon dioxide melting process can be collected by CCD camera in a non-contact way. A series of image analysis algorithms such as image segmentation and edge detection can be used to extract the contour features of the original image data. In order to find the centroid position of silicon dioxide， it is necessary to extract the contour of silicon dioxide by DBSCAN density clustering algorithm. Based on the real-time centroid position data of silicon dioxide， the melting rate of silicon dioxide can be further found， which provides a theoretical basis for improving the utilization rate of iron tailings.

Key words：image segmentation; edge detection; DBSCAN algorithm; the melting process of silicon dioxide; the centroid position of silicon dioxide

1 引 言

高爐渣是高爐煉鐵過程中產生的固體廢料，鐵尾礦是選礦過程中產生的殘渣，這些廢棄物如不能很好地處理會造成大量的資源浪費和環境污染[1]。目前，利用高爐渣和鐵尾礦來制造礦渣棉可以將這些材料變廢為寶[2]，其實現關鍵在于如何監測高溫熔池中鐵尾礦的融化行為?？紤]到高溫熔池中的溫度高達1500℃，這對于監測設備來說是非常大的考驗，而采用CCD攝像機可以在非接觸的方式下采集融化過程的圖像數據，基于此分析其融化過程是一個新的研究思路。

二氧化硅是鐵尾礦的主要成分[3]，同時也是其最難融化的成分，因此二氧化硅的融化行為可以間接地表示鐵尾礦的融化行為?；诙趸璧娜诨袨闃嫿ㄨF尾礦融化模型可以為渣棉制備過程中尾礦添加和熱補償提供了指導，從而間接改善了高爐渣直接成纖維技術。本文基于DBSCAN聚類算法提出了一種鐵尾礦融化過程中二氧化硅的質心定位算法，可以從CCD圖像數據中提取出二氧化硅的質心，為鐵尾礦融化模型的構建奠定基礎。

2 基于DBSCAN的二氧化硅融化過程中質心定位

我們研究的原始數據來源于通過CCD設備以1秒/幀的采樣頻率采集到的高溫熔池中二氧化硅的熔融過程，共有114張圖片。我們將首先通過裁剪和灰度化操作對原始圖像數據進行預處理，然后選擇合適的算子提取圖像的邊緣特征，并通過二值化和密度聚類的方法對邊緣特征進行降噪處理，在此基礎上做進一步地聚類提取二氧化硅的輪廓信息，最后通過質心計算公式獲取二氧化硅的質心坐標。

2.1數據預處理

CCD攝像機采集到的原始圖像數據如圖1（a）所示，其中左側黑色的部分為坩堝容器外的背景，右側高亮的圓形盆狀物是容器坩堝。由于在數據采集的過程中CCD攝像機和坩堝的位置都是固定的，因此可以對每張圖片進行裁剪的操作，從而除去左側無用的部分從中提取出僅包含坩堝和二氧化硅的圖像，如圖1（b）所示，裁剪后的圖片中中間偏左的不規則塊狀物即為二氧化硅。

另外， CCD攝像機采集到的原始數據既有彩色圖像也有灰度圖像，為了保證數據的一致性，同時考慮到顏色特征對于質心的定位影響不大，因此首先通過加權平均法對彩色圖像進行灰度化處理，處理前后的二氧化硅融化圖像如圖2所示。

2.2邊緣檢測

為了從灰度圖像中提取二氧化硅的質心，需要首先從圖像中提取出二氧化硅的輪廓[4]，對于圖2中的灰度圖像分別使用了Canny[5]、Scharr[6]和Kirsch[7]三種常用的邊緣檢測算子進行處理，處理的結果如圖3所示。

從圖3可以看出用Canny算子處理后，輪廓特征丟失比較嚴重，二氧化硅的輪廓不是很清晰;而在用Kirsch算子處理的結果中，對比原圖像，二氧化硅輪廓偏離正確位置。因此，相比較而言Scharr算子可以較好地提取出圖像中二氧化硅的輪廓信息。

2.3特征去噪

從圖3中可以發現經過邊緣檢測處理后得到的邊緣特征圖像中仍然存在大量的噪聲，我們無法直接基于此獲取到二氧化硅的邊緣像素坐標，因此還需要對邊緣特征圖像進行去噪處理。

2.3.1二值化處理

一般，屬于圖像輪廓特征的像素點的顏色會比較深，而噪聲則大多是顏色比較淺，因此可以通過設定閾值θ對灰度圖像進行二值化處理，以消除大部分的噪聲。為了選取合適的閾值θ進行二值化處理，設定閾值初始值為70，間隔為10，進行二值化處理，處理后的效果圖如圖4所示。

從圖4可以發現，閾值較低時，雖然圖像中的輪廓特征明顯，但噪聲抑制的較差，不能達到特征去噪的效果，而當閾值設置較高時，噪聲得到明顯抑制的同時輪廓特征也同時被抑制。當閾值設定為90時，輪廓特征完整，且噪聲得到的一定的抑制，處理效果比較好。此時，將二值化處理后圖像中像素為0的點的像素坐標記錄下來，構成數據集Dn以待后續分析。

2.3.2密度聚類

通過二值化處理后的圖像中還存在部分噪聲，這些噪聲對于二氧化硅輪廓的提取仍會造成一定的干擾。同時，考慮到屬于圖像輪廓特征的像素點，在距離上一般呈現連續的特征，而噪聲則大多是在距離上離散的像素點，因此我們將采用基于密度的聚類方法（DBSCAN算法）對邊緣特征圖像做進一步地去噪處理。

DBSCAN是一種密度聚類算法，通過使用一組關于“鄰域”的參數（ε，MinPts）來描述樣本分布的緊密程度。與劃分和層次聚類方法不同，它聚類的時候不需要預先指定簇的個數，最終的簇的個數也不定。DBSCAN算法需要先設定在ε鄰域內的半徑Eps和點的數量MinPts，之后可將數據點分為三類，分別是核心點、邊界點和噪聲點，核心點是指在半徑Eps內含有的點數量超過MinPts數目的那些點，邊界點是指在半徑Eps內點的數量小于MinPts，但落在核心點鄰域內的點，噪聲點是指不在任何鄰域內的點。然后刪除噪音點。接下來為距離在Eps之內的所有核心點之間賦予一條邊，每組連通的核心點組成一個簇，最后將每個邊界點分類到與之關聯的核心點的簇中。

我們設定DBSCAN的掃描半徑（eps）為10，最小包含點數（minPts）為6，以2.3.1節中得到的輪廓特征點的集合Dn為例，將15960點共分為150類，并將每個類中數據點數量小于300的類刪除，得到第二次特征去噪的結果即為新的輪廓特征點的集合Dn，如圖5所示。

2.4輪廓提取

考慮到2.3節中得到的去噪后的輪廓特征集合Dn中不僅包含二氧化硅的輪廓還包含坩堝容器的輪廓信息，因此我們還需要對其進行二次聚類，依舊采用DBSCAN算法，調整其聚類參數，將掃描半徑（eps）設置為25，最小包含點數（minPts）設置為10，得到的二次DBSCAN聚類的前9類聚類點的集合結果如圖6所示。

從圖6可以看出，通過二次DBSCAN聚類可以將特征圖像中屬于二氧化硅的輪廓的像素點和屬于坩堝容器的像素點區分開，聚類結果中的第4類中包含的像素點集合DSiO?即為圖像中二氧化硅的輪廓。

2.5質心搜索

針對2.4節中提取出的二氧化硅輪廓特征點的集合DSiO?，通過以下公式計算二氧化硅的形心：

其中（xi，yi）為二氧化硅輪廓特征點的集合DSiO?中的像素點，N為二氧化硅輪廓特征點的集合DSiO?中的像素點的個數。

對于2.4節中得到的特征點的集合，計算得到的二氧化硅形心如圖7（a）所示，其中黑色的像素點即為計算得到的二氧化硅的形心。采用上述方法，求出用CCD攝像機采集到的二氧化硅的熔融過程的圖像序列中每一張圖像中二氧化硅的形心，利用曲線擬合這些形心點坐標，即可獲得形心運動軌跡，結果如圖7（b）所示，其中紅色的折線為二氧化硅熔融過程中的形心運動軌跡。

3 結論

針對用CCD攝像機獲取到的高溫熔融過程中的二氧化硅的圖像數據，本文提出了一種基于DBSCAN聚類的輪廓特征提取方法，可以有效地從圖像中提取二氧化硅的輪廓，并在此基礎上檢測到二氧化硅在熔融過程中質心的位置變化過程?；诖?，可以進一步地提取二氧化硅的形狀、周長、面積等數據來構建二氧化硅的邊緣輪廓特征指數，并以此評估二氧化硅的融化速率，為鐵尾礦的融化過程提供參考依據，提高材料的利用率。

參考文獻：

[1] 劉愛平，曹苗佳，張雄濤.鐵尾礦規?；醚芯楷F狀及其利用途徑[J].中國金屬通報，2019（11）：279-280.

[2] 邢宏偉，劉振超，梁偉，等.鐵尾礦在高爐熔渣調質過程中熱力學行為的研究[J].鋼鐵釩鈦，2017，38（5）：90-95，100.

[3] 于洪浩.鞍山高硅鐵尾礦的增值化利用研究[D].沈陽：東北大學，2010.

[4] 張晗，錢育蓉，王躍飛，等.基于圖像邊緣增強與弱化的邊緣檢測[J].計算機工程與設計，2019，40（11）：3106-3110.

[5] Canny J.A computational approach to edge detection[J].IEEE Transactions on Pattern Analysis and Machine Intelligence，1986，PAMI-8（6）：679-698.

[6] Scharr H，Weickert J.An anisotropic diffusion algorithm with optimized rotation invariance[M]//Informatik aktuell. Berlin， Heidelberg： Springer Berlin Heidelberg， 2000： 460-467.

[7] Kirsch R A.Computer determination of the constituent structure of biological images[J].Computers and Biomedical Research， 1971，4（3）：315-328.

[8] Ester M，Kriegel H P，Sander J，et al.A density-based algorithm for discovering clusters in large spatial databases with noise[EB/OL].1996.

[9] 王光，林國宇.改進的自適應參數DBSCAN聚類算法[J/OL].計算機工程與應用：1-8[2020-02-07].

[10] Gan J H，Tao Y F.DBSCAN revisited：mis-claim，un-fixability，and approximation[C]//Proceedings of the 2015 ACM SIGMOD International Conference on Management of Data - SIGMOD '15. May 27-June 4， 2015. Melbourne， Victoria， Australia. New York： ACM Press， 2015： 519-530.

【通聯編輯：唐一東】