和聲搜索改進的形態學分析在庫區漂浮物體量預估中應用的研究

花勝強，陳 意，鄭慧娟，袁 帥，高 磊

(南瑞集團(國網電力科學研究院)有限公司，江蘇南京211106)

0 引 言

隨著我國制造產業的迅猛發展，人民群眾生活水平的飛速提升，以及水利水電建設的大規模開展，水電廠、水利樞紐的庫區開始出現和聚集大量的漂浮物，其成分主要包含生活垃圾、工業廢棄物及樹干枝葉等。由于水電廠或水利樞紐一般承擔著防洪、供水、發電、航運及旅游等社會和民生功能，防漂浮物的存在不僅對運管單位的水工設施、交通安全、機電設備等造成危害，更會污染水體，破壞地區生態環境，因此，及時地監控、估算庫區漂浮物，從而合理制定和安排漂浮物清除作業，無疑具有極其重要的經濟和安全效益[1-3]。

當前，對于庫區漂浮物的實時監測主要依賴人工巡檢及機器視覺等方式。人工巡檢實施簡單但耗費人力，也無法實現實時監測，對于當前漂浮物的體量難以量化預估，不利于后續組織安排清除作業船只。基于遠程攝像的機器視覺是一種借助圖形圖像處理方法進行成像中目標物檢測、分割和識別的技術，具備自動化、智能化等特性，是進行庫區漂浮物在線監控和分析的良好途徑。已有文獻提出基于圖像灰度和邊緣曲線的BP神經網絡用于樹枝和漂浮瓶的判別模型，基于Faster R-CNN、R-FCN和SSD等多種深度學習網絡用于落葉和水草的目標檢測，以及基于AlexNet卷積神經網絡用于水面塑料制品的識別等，均取得了一定的識別效果[4-5]。

BP神經網絡及深度學習算法關于復雜非線性映射關系良好的擬合能力，在圖形圖像分析識別中確實具有一定的優勢，但其限制條件也很顯著，如網絡模型的精度是由數據驅動的，需要大樣本的漂浮物圖像集并一一加以標注，數據預處理工作量巨大；網絡模型的結構設計困難且復雜，模型超參數量巨大且不易設定，模型訓練耗時耗資源，導致產出門檻較高等[6]。

針對以上問題，本文提出了一種庫區漂浮物體量預估的方法，首先在完成圖像預處理后，進行樣本圖像中漂浮物目標區域的人工標注；其次，基于形態學方法生成每個樣本圖像的梯度特征空間；然后，基于滑動梯度卷積及和聲搜索算法進行漂浮物識別模型的訓練和優化，并產出庫區漂浮物圖像識別模型；最后，根據預定義的像素面積比、漂浮物平均厚度、漂浮物平均密度等參數，給出庫區漂浮物的面積、體積和質量的預估，以便于后續管理單位的船舶清除調度安排。由于是基于樣本圖像的梯度特征空間進行模型的學習訓練，故該方法只需小樣本量即可生成較高精度的漂浮物識別模型，無論是在樣本獲取、樣本預處理還是模型訓練的時間和空間資源消耗等環節，其工作效率都比較高[6]。

1 基于內容的圖像壓縮

當前圖像進行識別前通常要求壓縮到指定尺寸上，如果原圖與指定尺寸的寬高比差異較大，會導致圖像內容發生形變失真；另外，圖像中的目標區域也可能跟隨著圖像壓縮處理而同步損失，從而丟失了寶貴的目標內容信息。由此，本文中提出一種基于內容的圖像壓縮改進算法，以盡可能保持原圖像中含有重要信息區域的比例，從而避免縮放造成的失真。

1.1 算法原理

當前常用的圖像壓縮算法，是類似于沿著圖像橫向或者縱向做均勻等間隔的降采樣。此法雖然邏輯上簡單明了，但也造成了無差別的內容裁減。由此，本文提出基于圖像能量，而非均勻等間隔選點的策略來確定待壓縮路徑集合。

首先，定義圖像中某點P的能量，即

(1)

其次，定義圖像中的一條路徑為如下所述圖像點的集合，即從圖像的頂部或者左側的任意一點開始，依次向下或者向右相鄰的3個點中擇一移動，如此重復移動，直至達到圖像最底部或者最右側，則所有選擇遍歷過的圖像點即組成一條路徑。

再次，定義路徑的能量，即為該路徑上所有像素點的能量之和。

最后，可得出本算法的目標為，每次找到能量最小的一條路徑，然后刪除此路徑上所有像素點，如此循環迭代，直至達到壓縮尺寸為止。

1.2 算法流程

該算法流程如下：

(1)根據原圖尺寸和壓縮后的目標尺寸，確定橫向與縱向各需要循環的次數。

(2)計算圖像中每個像素點的能量。

(3)基于求解矩陣最短路徑的動態規劃算法，找出圖像當前橫向或縱向方向上的能量最小的路徑。

(4)移除最小路徑上所有像素點，得到新圖像。

(5) 重復步驟(2)至步驟(4)，直到新圖像橫向與縱向都達到目標尺寸為止。

2 形態學分析

2.1 分析流程

形態學分析通常基于數學工具從圖像中提取出表達和描繪區域形狀特征的圖像分量，其基本運算包括腐蝕和膨脹、開運算和閉運算、形態學梯度、空間變換等等。由此，本文提出基于形態學分析的庫區漂浮物識別方法流程如下：

(1)能量矩陣生成。基于公式(1)定義的圖像點能量定義，生成與原圖像素點一一對應的二維能量矩陣。顯然，能量矩陣的尺寸與原圖一致，矩陣中每個元素值代表了原圖中相應位置像素點的能量。

(2)能量矩陣分塊與選擇。首先，基于給定尺寸的滑動窗口，將能量矩陣分割成多個窗口大小的矩陣塊，并計算每個矩陣塊的能量，其中，矩陣塊能量即為塊內各個元素值之和；然后，基于給定的能量閾值，選擇所有能量不低于閾值的矩陣塊并加以標注。

(3)目標矩陣塊形態學處理。將上步中標注的目標矩陣塊，基于指定的尺寸分別進行腐蝕、膨脹處理，以消除目標矩陣塊所構成的整體識別區域內部中的氣泡、區域邊緣的毛刺等，使得整體的識別區域得到更好的內部連通性與邊緣平滑性。

(4)目標塊聯通化識別與過濾。首先，將各目標矩陣進行聯通化分析和識別，以將所有相連的矩陣塊合并成一個整體區域。其次，基于指定的最小區域尺寸閾值，去除尺寸過小的識別區域，以減少圖像中噪音部分的干擾，提高識別的準確率。

2.2 基礎參數

根據上述流程可知，形態學分析共含有6個基礎參數，包括分割能量矩陣的滑動窗口像素尺寸(窗口高WH、寬WL)、選擇矩陣塊的最小能量閾值E、對矩陣塊進行腐蝕操作的像素尺寸C、對矩陣塊進行膨脹操作的像素尺寸S、單堆漂浮物的最小區域像素尺寸A。顯而易見，這些參數的取值設定會直接決定圖像中漂浮物區域識別的精度高低，而各個參數取值的排列組合的復雜度，導致尋找出最優參數取值組合近似為NP問題，故考慮采用基于群體仿生智能的啟發式搜索算法，來尋找基礎參數的最優取值組合。

3 和聲搜索算法

3.1 算法原理

和聲搜索算法是一種借鑒了樂隊進行音樂創作過程的啟發式全局搜索算法，在多個組合優化領域中得到了成功應用，并展示出了較傳統群體智能仿生算法更好的尋優能力和運算性能。樂隊創作時，每位樂師演奏一種樂器，所有樂器的演奏加起來對應一組和聲，樂師們會不斷調整各自演奏的音調，以求達到整體最好的和聲效果。和聲搜索算法正是基于此原理，將各樂器視為優化問題中的決策變量，各樂器發出的音調視為變量取值，各樂器音調組合成的和聲視為一組解向量，樂師們對和聲的評價視為優化問題中目標函數對解向量好壞的評估。和聲算法的最優化過程，就是對每一組和聲進行評價，如果達不到要求，各樂器就繼續各自調整音調，直到得到一組滿意的和聲為止[7]。

3.2 算法流程

和聲搜索算法包含如下幾個基本參數：①和聲庫大小，表征用于存儲候選最優和聲的記憶庫的大小。②和聲庫取值概率，表征每次迭代時候直接從和聲庫中選擇的可能性，否則就按照各決策變量的取值范圍構成的解空間，隨機生成一組和聲。③微調概率，表征當直接從記憶庫中選擇和聲時候，對選中的和聲中各個變量取值再進行微調的概率。④微調帶寬，表征當要對選中的和聲進行微調時的調整幅度大小。⑤最大迭代次數，表征上述整個尋優調整過程最大迭代的次數。和聲搜索算法流程如圖1所示。

圖1 和聲搜索算法流程

3.3 算法參數分析與設置

和聲搜索算法基本參數的取值，直接關系到后續算法尋優的精度和性能，而其中和聲庫取值概率、微調概率、微調帶寬更是尤為重要，現將此3個參數的取值策略分析如下：

(1)和聲庫取值概率，當取值較大時，有利于局部解空間內的最優值搜索收斂，而取值較小時，有利于增加候選最優和聲的多樣性，防止陷入局部最優的早熟，故可以采取動態更新策略，算法啟動時取值較小以增大搜索的多樣性，隨著迭代次數增加而逐漸擴大，以便更聚焦最優局部空間內尋優，其概率取值為

(2)

式中，m為最大迭代次數；i為當前尋優迭代的次數；Pi為第i次時的和聲庫取值概率；Ps為初始的和聲庫取值概率取值；Pe為最終的和聲庫取值概率取值。

(2)微調概率，與和聲庫取值概率恰好相反，取值較大時利于搜索的多樣性，較小時利于局部尋優收斂，故仍可以采取動態更新策略，算法啟動時取較大值以增大搜索的多樣性，隨著迭代次數增加而逐漸減小，其概率取值為

(3)

式中，Qi為第i次時的微調概率；Qs為初始微調概率值；Qe為最終微調概率值。其他含義同上。

(3)微調帶寬，原始和聲搜索算法中一般取固定值，即各個決策變量的微調帶寬取值初始化設定后，運行過程中將不再改變，本文仍然考慮動態取值策略，其帶寬取值為

(4)

式中，Ri為第i次時的微調帶寬；R為初始設定帶寬；Ran(-1，1)為位于區間[-1，1]內的隨機浮點數。

4 方法流程

4.1 圖像預處理

和聲搜索改進的形態學分析在庫區漂浮物體量預估中的應用流程詳細步驟說明如下[8-9]：

(1)圖像壓縮。基于本文提出的壓縮算法，將所有樣本圖像壓縮至指定尺寸，以便于后續統一的圖像分析處理，并極大地降低分析計算量。

(2)灰度轉換與目標區域標注。將壓縮后的樣本圖像統一變換為灰度圖，變換方法如下

Gray=(38R+75G+15B)?7

(5)

式中，Gray為圖像中某像素點轉換后的灰度值；R為該像素點在原圖像RGB空間中的紅色分量值；G為綠色分量值；B為藍色分量值。同時，將圖像中的所有漂浮物區域進行標注，并記錄下區域坐標集合，以便后續識別學習時進行精度判定。

4.2 模型訓練

(1)基于和聲搜索算法進行形態學分析模型的訓練。對于形態學分析中涉及到的基礎算法參數，則基于和聲搜索算法進行各參數最優化取值組合的自動化尋優；作為優化目標的損失函數，定義為形態學分析出的漂浮物區域，與標注的漂浮物區域相差異的總像素點數。顯然，差異點總數越小，當前的形態學分析基礎參數取值組合的分析精度越高。

(2)得到和聲搜索改進的庫區漂浮物形態學分析識別模型并進行體量預估基于和聲搜索算法的訓練完成后，得到最終的形態學分析模型，模型輸出內容即為輸入圖片中判定為漂浮物區域的總像素點數N，以及相應的像素點坐標集合。

4.3 漂浮物體量預估

分析識別時，將待處理的庫區圖片按照4.1的步驟(1)和(2)預處理完成后，代入步驟(3)產出的形態學分析模型，得到分析識別出來的漂浮物區域信息和總像素點數N，然后，按照設定好的3個預估參數：①像素-面積轉換比A，m3/像素。該參數表征了圖像中1個像素點代表了水庫庫區中的實際面積大小。②漂浮物平均深度H，m。該參數代表了庫區水上漂浮物的平均厚度，包括水面上下的厚度。③漂浮物平均密度P，kg/m3。根據以上參數，得到漂浮物總體積V為

V=NAH

(6)

由此，漂浮物總質量M為

M=VP

(7)

5 實例及分析

以湖北省某大型水庫上游庫區遠程攝像頭所拍圖像為例，除去因天氣、攝像頭遮擋等導致質量不佳的成像外，共計得到63張有效樣本圖像，統一按照流程(1)、(2)進行預處理后，將其中的48張作為訓練樣本集，15張作為測試樣本集。

5.1 基于內容的圖像壓縮示例

基于圖像內容的壓縮效果如圖2所示。從圖2可以看出，本文算法對于目標區域的損失極小，基本上保留了漂浮物區域的原始細節信息。

圖2 基于內容的圖像壓縮算法對比

5.2 和聲搜索算法改進的形態學分析

本實例中和聲搜索算法的基礎參數設定如下：和聲庫大小為200，和聲庫取值概率起始值為0.10，最終值為0.90，微調概率起始值為0.80，最終值為0.25，微調帶寬取值為1.0，最大迭代次數為6 000。

此外，和聲庫初始化時形態學分析的基礎參數取值組合時，以滑動窗口高WH=6、寬WL=6、矩陣塊最小能量閾值E=180、矩陣塊腐蝕尺寸C=5、矩陣塊膨脹尺寸S=5、單堆漂浮物的最小區域像素尺寸A=150作為種子，隨機生成各種參數的取值組合。

訓練完成后，以產出的形態學分析模型對測試集圖像進行漂浮物區域識別，并輸出識別區域，效果如圖3、圖4所示。其中，能量矩陣塊圖為經過能量閾值選擇過的矩陣塊集，形態學處理圖為經過腐蝕、膨脹、矩陣塊聯通化處理的圖像。

圖3 和聲搜索改進的形態學分析模型識別效果示例1

圖4 和聲搜索改進的形態學分析模型識別效果示例2

圖3和圖4的形態學分析模型在訓練集中的整體精度達到95.61%左右，在測試集中的平均識別精度也達到了94.37%左右。以圖5為例，左上為測試集中的原始庫區圖像，右上為達到能量閾值的矩陣塊集合，左下為經過腐蝕、膨脹、聯通化和最小面積區域過濾的識別圖像，右下為最終識別出的區域與原圖疊加對比的效果圖。從圖5可以看出，本方法基本以較好的精度達到了庫區漂浮物識別的目標。最后，水庫單位根據多年來運行管理經驗，基于漂浮物的常見構成部分和比例設定了漂浮物平均深度、平均密度，最終得出了庫區當前漂浮物的體積與質量的估計值。

圖5 和聲搜索改進的形態學分析模型整體效果對比

本文提出的庫區漂浮物形態學分析識別方法已在湖北省某大型水庫上得到了實際應用，投運以來的運行實踐表明，其計算精度良好，魯棒性強，為水庫的庫區漂浮物監控管理工作提供了實時、可靠的決策輔助支持。

6 結 語

本文提出了基于和聲搜索算法改進的庫區漂浮物形態學分析識別方法，其核心在基于數學形態學的基礎上提出了漂浮物識別的方法論、演算函數以及損失函數定義，然后利用和聲搜索算法在形態學分析參數的取值組合空間內尋求最優解，從而得到最終的分析識別模型。相比較目前流行的深度學習圖像識別，該方法邏輯清晰簡明，易于實現、部署與迭代，計算資源需求極小，且更適用于小樣本量下的庫區漂浮物分析識別處理。

庫區漂浮物形態學分析識別方法在體量預估上還有持續改進的空間，包括雨雪、霧霾、強光、暗光、陰影等條件下，圖像識別前需考慮增加銳化、增強、濾波等預處理操作；此外，體量預估中人工設定像素-面積轉換比、漂浮物平均深度、漂浮物平均密度等參數的精度難以把握，如何繼續基于圖形圖像處理技術給出更精確的估值等，都是下一步工作的重點與方向。