陳吉，劉曉榮，楊佳雯，陳冶秋，王誠，歐夢媛，吳嘉儀，俞尤嘉，李開，陳鵬，陳峰

1.南京醫(yī)科大學(xué)法醫(yī)學(xué)系，江蘇 南京 211166；2.南京市公安局水上分局，江蘇 南京 210036

硅藻檢驗(yàn)是溺死診斷的重要手段[1]。通過對(duì)水中尸體進(jìn)行硅藻檢驗(yàn)，可判斷其死亡原因，并可推斷入水地點(diǎn)[2-4]。硅藻檢驗(yàn)方法有消化法、化學(xué)元素法、分子生物學(xué)檢驗(yàn)法等，其中強(qiáng)酸消化-光學(xué)顯微鏡檢驗(yàn)法因操作簡單、設(shè)備要求低，是基層法醫(yī)學(xué)鑒定實(shí)踐中常用的硅藻檢驗(yàn)方法，但是此方法具有雜質(zhì)消化不全、離心損失大、檢出率低、開放操作易污染等問題[5]。

近年來興起的微波消解-真空抽濾-掃描電子顯微鏡法[6]因具有雜質(zhì)消化完全、操作污染少、硅藻富集率高、陽性檢出率高等優(yōu)勢，已在部分公安檢案中開展應(yīng)用。然而掃描電子顯微鏡產(chǎn)生的圖像數(shù)據(jù)量龐大，不含硅藻的陰性圖像多，人工從上千甚至上萬張圖像中精確識(shí)別硅藻需要耗費(fèi)大量的時(shí)間和精力。硅藻電子顯微鏡圖像的快速、準(zhǔn)確自動(dòng)化識(shí)別已成為一個(gè)亟待解決的問題。

隨著人工智能技術(shù)的不斷發(fā)展，基于深度學(xué)習(xí)進(jìn)行圖像識(shí)別已在多領(lǐng)域應(yīng)用，與人工識(shí)別相比顯著提升了識(shí)別速度并可取得與人工識(shí)別相當(dāng)?shù)臏?zhǔn)確率[7-9]。深度學(xué)習(xí)是模擬人腦信號(hào)傳導(dǎo)機(jī)制構(gòu)建人工神經(jīng)網(wǎng)絡(luò)，通過大量數(shù)據(jù)訓(xùn)練提取數(shù)據(jù)特征，從而完成從特征到目標(biāo)任務(wù)的映射。卷積神經(jīng)網(wǎng)絡(luò)（convolutional neural network，CNN）是深度學(xué)習(xí)中的一種重要模型，特點(diǎn)是利用多層卷積運(yùn)算實(shí)現(xiàn)特征提取、降維和結(jié)果預(yù)測，已廣泛應(yīng)用于計(jì)算機(jī)視覺領(lǐng)域。

目前已有相關(guān)研究涉及人工智能硅藻圖像自動(dòng)識(shí)別[10-12]，但識(shí)別模型較少，識(shí)別速度和性能均有待提升。為更好地解決硅藻電子顯微鏡圖像識(shí)別問題，提高識(shí)別速度和自動(dòng)識(shí)別性能，本研究基于深度學(xué)習(xí)CNN 中的YOLOv3 模型開展硅藻電子顯微鏡圖像自動(dòng)化識(shí)別模型的構(gòu)建，并在實(shí)際硅藻檢驗(yàn)案例中進(jìn)行應(yīng)用，探討該模型的優(yōu)勢和可待改進(jìn)之處。

1 材料與方法

1.1 硅藻圖像來源

本研究使用的所有硅藻電子顯微鏡圖像均由南京市公安局水上分局刑事警察大隊(duì)提供，圖像通過微波消解-真空抽濾-掃描電子顯微鏡法采集于長江流域溺死者的肺、肝、腎組織。圖像采集設(shè)備為Phenom XL G2 臺(tái)式掃描電子顯微鏡、DiatomScope 軟件（荷蘭PhenomWorld 公司），設(shè)置掃描半徑8.0 mm（涵蓋整個(gè)樣本區(qū)域）、放大倍數(shù)1 500×、圖像分辨率1 024×1024，形成的每個(gè)視場大小為0.176mm×0.176mm，完整掃描一個(gè)樣本所有視場生成的圖片總數(shù)為6 392 張。

本研究在上述所有樣本采集圖像中選取25 000張經(jīng)2 名法醫(yī)病理學(xué)專家肉眼確認(rèn)含完整硅藻數(shù)量大于1 的硅藻電子顯微鏡圖像，作為建立數(shù)據(jù)集所使用的樣本圖像集。同時(shí)選取樣本圖像集以外8例溺死者的肺、肝、腎組織的106 272 張硅藻電子顯微鏡圖像，以檢驗(yàn)?zāi)Ｐ蛯?shí)際應(yīng)用效果。

本研究僅對(duì)南京市公安局水上分局刑事警察大隊(duì)提供的相關(guān)硅藻電子顯微鏡圖像進(jìn)行處理和使用，符合醫(yī)學(xué)倫理學(xué)相關(guān)條款規(guī)定。

1.2 數(shù)據(jù)集構(gòu)建

由硅藻檢驗(yàn)經(jīng)驗(yàn)豐富的2名法醫(yī)病理學(xué)專家使用LabelImg v1.4.3 軟件（Tzutalin 開發(fā)；http://tzutalin.github.io/labelImg/）共同對(duì)樣本圖像集中的硅藻進(jìn)行標(biāo)注（標(biāo)注時(shí)所有完整硅藻和碎片硅藻均視為硅藻）。為了讓模型達(dá)到更好的訓(xùn)練效果，將標(biāo)注所得圖像數(shù)據(jù)進(jìn)行歸一化處理，另外使用旋轉(zhuǎn)、平移、縮放、裁剪、變形等方式進(jìn)行數(shù)據(jù)增強(qiáng)，從而擴(kuò)大訓(xùn)練圖像樣本量。

1.3 模型訓(xùn)練與測試

選取CNN 中的YOLOv3 模型作為識(shí)別模型的基礎(chǔ)。YOLOv3 模型是一種使用“端對(duì)端”方式訓(xùn)練的單階段目標(biāo)檢測模型[13-15]（圖1）。輸入的圖像經(jīng)過全卷積網(wǎng)絡(luò)（Darknet-53）進(jìn)行特征提取和特征融合后，輸出3 種依據(jù)不同張量劃分的特征圖[15]，特征圖中涉及目標(biāo)對(duì)象的每一個(gè)網(wǎng)格根據(jù)模型先前有監(jiān)督學(xué)習(xí)的訓(xùn)練結(jié)果預(yù)測存在對(duì)象的3 個(gè)邊界框的相應(yīng)信息，包括位置和大?。ㄖ行淖鴺?biāo)、框的寬度和高度）、置信度以及對(duì)象類別數(shù)。對(duì)于同一個(gè)目標(biāo)對(duì)象對(duì)應(yīng)的各邊界框按照各自的置信度進(jìn)行排序，以置信度最高的邊界框作為該目標(biāo)對(duì)象所在的位置。以識(shí)別硅藻電子顯微鏡圖像為例，模型將待識(shí)別的圖像轉(zhuǎn)化為416×416 個(gè)網(wǎng)格，經(jīng)過特征提取后輸出包括13×13、26×26、52×52 這3 種張量的特征圖。特征圖中涉及硅藻的每一個(gè)網(wǎng)格將預(yù)測3 個(gè)可能存在硅藻的邊界框，所有邊界框中置信度最高者為該硅藻在所屬圖像中的位置，該置信度即為模型預(yù)測邊界框內(nèi)對(duì)象為硅藻的可能性。

圖1 YOLOv3 模型原理示意圖Fig.1 The schematic diagram of the YOLOv3 model

模型識(shí)別時(shí)可設(shè)置不同的閾值，該閾值是模型最終判定是否為目標(biāo)對(duì)象的置信度標(biāo)準(zhǔn)，若識(shí)別對(duì)象置信度高于該閾值，則判定該對(duì)象為目標(biāo)對(duì)象硅藻，并顯示相應(yīng)的邊界框和置信度，閾值越高，識(shí)別結(jié)果是真正硅藻的可能性越大。

以7∶2∶1 的比例將1.2 節(jié)中構(gòu)建的數(shù)據(jù)集樣本劃分為訓(xùn)練集、驗(yàn)證集和測試集，基于python 3.6.8 進(jìn)行YOLOv3 框架搭建和模型訓(xùn)練。設(shè)置初始學(xué)習(xí)率為5×10-4，最小學(xué)習(xí)率為5×10-7；每次訓(xùn)練批次大小為2；每迭代200 次評(píng)估1 次訓(xùn)練總損失值——總損失值是對(duì)象類別、置信度、邊界框中心坐標(biāo)以及邊界框?qū)捄透? 項(xiàng)數(shù)值損失值的總和，體現(xiàn)模型誤差大小。最后根據(jù)模型收斂情況和損失值的最低值確定模型最終訓(xùn)練的迭代次數(shù)。

根據(jù)訓(xùn)練模型在驗(yàn)證集、測試集上的查準(zhǔn)率（precision rate，PR）和召回率（recall rate，RR）和各閾值下的F1 分?jǐn)?shù)實(shí)時(shí)評(píng)價(jià)模型訓(xùn)練效果（閾值取值范圍為0～1，以0.05 遞增）。PR是指在所有被預(yù)測為正例（即模型檢測的目標(biāo)對(duì)象）的樣本中，實(shí)際為正例的樣本比例；RR是指所有實(shí)際為正例的樣本中，被預(yù)測為正例的樣本比例；F1 分?jǐn)?shù)是PR和RR的加權(quán)平均，體現(xiàn)了模型在各閾值下的識(shí)別效能。最后結(jié)合模型在不同閾值下的RR和PR計(jì)算平均精度均值（mean average precision，mAP）評(píng)價(jià)模型性能。公式如下：

式中，TP為真陽性，即模型預(yù)測為硅藻而實(shí)際也為硅藻的數(shù)量；FP為假陽性，即模型預(yù)測為硅藻而實(shí)際非硅藻的數(shù)量；FN為假陰性，即模型預(yù)測為非硅藻而實(shí)際為硅藻的數(shù)量；AP為平均精度，即模型在每個(gè)閾值下的PR的加權(quán)平均值，權(quán)重為RR的增加值，平均精度相當(dāng)于某一類別PR-RR曲線（P-R曲線）下的面積，k為總類別數(shù)。

由公式可得，當(dāng)RR越高，模型識(shí)別硅藻的敏感性越高；當(dāng)PR越高，模型識(shí)別硅藻的假陽性率越低，模型特異性越高；當(dāng)mAP越高時(shí)，模型檢測性能越好。

1.4 模型應(yīng)用

從1.1 節(jié)中選取的8例溺死者的肺、肝、腎組織的106 272 張硅藻電子顯微鏡圖像用于測試該模型的實(shí)際應(yīng)用效果。其中，肝、腎組織樣本由于含硅藻數(shù)量較少（整個(gè)樣本中可少至個(gè)位數(shù)），需要對(duì)完整掃描整個(gè)樣本所得的所有電子顯微鏡圖像進(jìn)行人工識(shí)別，因此本研究對(duì)每例肝、腎組織樣本掃描所得的6 392 張圖像進(jìn)行識(shí)別研究；肺組織樣本由于含硅藻數(shù)量較多（整個(gè)樣本中可達(dá)幾千至幾萬個(gè)），在實(shí)際檢案中通常在6 392 個(gè)視場中隨機(jī)選取200～300 個(gè)視場掃描圖像進(jìn)行人工識(shí)別，在本研究中為了更好地探究模型識(shí)別肺組織硅藻的RR和PR，對(duì)每例肺組織樣本掃描所得的500 張電子顯微鏡圖像進(jìn)行識(shí)別研究。

該識(shí)別模型的運(yùn)行流程如圖2 所示，使用者只需選擇硅藻電子顯微鏡圖像所在的文件夾，并設(shè)置篩選閾值，模型會(huì)將識(shí)別出硅藻置信度大于閾值的目標(biāo)對(duì)象框出并顯示相應(yīng)的置信度，有目標(biāo)對(duì)象和沒有目標(biāo)對(duì)象的圖像將被分別存放至相應(yīng)的文件夾。

圖2 硅藻電子顯微鏡圖像識(shí)別模型運(yùn)行流程圖Fig.2 The flow chart of the diatom identification model of scanning electron microscope images

由于隨著閾值設(shè)定的改變，模型的RR和PR會(huì)相應(yīng)地改變，通常設(shè)定的閾值越高，RR越低、PR越高。為了確定能夠使模型取得最佳性能的閾值大小，選取0.4、0.6、0.8 3 種閾值進(jìn)行識(shí)別效果評(píng)價(jià)。每例樣本的硅藻電子顯微鏡圖像分別接受2 名法醫(yī)病理學(xué)專家的肉眼確認(rèn)和該識(shí)別模型的自動(dòng)識(shí)別，記錄識(shí)別速度和識(shí)別結(jié)果，并將人工識(shí)別和模型識(shí)別的結(jié)果相對(duì)比，計(jì)算模型識(shí)別在實(shí)際應(yīng)用中的RR和PR。

模型識(shí)別將完整硅藻和碎片硅藻統(tǒng)一識(shí)別為硅藻，實(shí)際應(yīng)用中人工統(tǒng)計(jì)時(shí)需要將兩者區(qū)分開。對(duì)于碎片硅藻的識(shí)別，采用如下規(guī)則：面積大于50%的不完整硅藻視為完整硅藻，面積小于50%的不完整硅藻視為碎片硅藻，完整硅藻和碎片硅藻分2 類計(jì)數(shù)。由于模型識(shí)別出的假陽性硅藻無法判斷是對(duì)于完整硅藻還是碎片硅藻的假陽性，故分別計(jì)算完整硅藻、碎片硅藻和總體PR時(shí)將所有假陽性硅藻數(shù)包括在內(nèi)。

1.5 統(tǒng)計(jì)分析

人工識(shí)別和模型識(shí)別所得的硅藻陽性圖像數(shù)、不同組織中不同閾值下模型識(shí)別的RR、PR和F1 分?jǐn)?shù)用均數(shù)±標(biāo)準(zhǔn)差（）表示，采用SPSS 22.0軟件（美國IBM公司）進(jìn)行t檢驗(yàn)，檢驗(yàn)水準(zhǔn)α=0.05。

2 結(jié)果

2.1 模型訓(xùn)練情況和性能評(píng)價(jià)

圖3 為訓(xùn)練過程中總損失值的收斂曲線。圖中隨著迭代次數(shù)增加，模型迅速收斂并將損失值保持在0.2 以下，當(dāng)?shù)螖?shù)為128 600 次（即紅線標(biāo)注處）時(shí)總損失值達(dá)到最低，以此作為模型的最大迭代次數(shù)以取得理想的訓(xùn)練效果。

圖3 模型訓(xùn)練的收斂曲線Fig.3 Convergence curve of model training

訓(xùn)練完成后的模型在驗(yàn)證集和測試集上各閾值的RR、PR及F1 分?jǐn)?shù)見表1。模型在驗(yàn)證集和測試集的mAP分別為94.8%和94.3%；各閾值下的平均RR分別為81.2%和81.5%。

表1 各閾值下識(shí)別模型的RR、PR 和F1 分?jǐn)?shù)Tab.1 The precision rate，recall rate and F1 score of the identification model under different thresholds

2.2 模型應(yīng)用效果

2.2.1 識(shí)別速度和硅藻陽性圖像檢出量

對(duì)于每例肝、腎組織樣本的硅藻電子顯微鏡圖像，人工對(duì)6 392 張圖像逐一識(shí)別，需要花大于3 h，而使用該模型約15 min 內(nèi)完成，且不同的設(shè)定閾值下所用的識(shí)別時(shí)間相同，每個(gè)樣本模型識(shí)別和人工復(fù)核的時(shí)間總和不超過20 min，圖像識(shí)別速度可達(dá)0.143 s/張，所用時(shí)間是人工識(shí)別的1/9。人工識(shí)別每例肝、腎組織樣本6 392 張電子顯微鏡圖像檢出的硅藻陽性圖像（即含有硅藻的圖像）均值分別為（4.3±1.7）、（9.0±8.6）張，模型在0.4、0.6 和0.8 閾值下檢出肝組織硅藻陽性圖像均值分別為（269.4±197.8）、（218.5±165.9）、（110.4±85.6）張，檢出腎組織硅藻陽性圖像均值分別為（289.8±103.8）、（202.6±60.3）、（117.0±29.4）張，t檢驗(yàn)顯示與人工識(shí)別檢出圖像數(shù)之間差異均有統(tǒng)計(jì)學(xué)意義（P＜0.05）。

對(duì)于肺組織硅藻電子顯微鏡圖像，模型對(duì)每張圖像的識(shí)別速度與肝、腎組織相同。人工識(shí)別每例肺組織樣本的500 張電子顯微鏡圖像檢出的硅藻陽性圖像均值為（368.5±176.5）張。模型在0.4、0.6 和0.8 閾值下檢出硅藻陽性圖像均值分別為（355.5±179.2）、（330.4±176.9）、（310.0±188.2）張，t檢驗(yàn)顯示與人工識(shí)別檢出圖像數(shù)之間差異均無統(tǒng)計(jì)學(xué)意義（P＞0.05）。

2.2.2 不同閾值下模型的檢測效能

從表2 可見，在0.4 閾值下，8例案例各組織的大部分硅藻（包括碎片硅藻）都能被正確檢出，并且肝、腎組織的硅藻RR達(dá)到了100%。隨著閾值的增加，硅藻RR在各組織均有所下降。從圖4 可見，當(dāng)閾值從0.4 上升至0.6 時(shí)，相同視野識(shí)別出的硅藻個(gè)數(shù)下降。

圖4 設(shè)定不同閾值時(shí)肺組織的硅藻識(shí)別結(jié)果Fig.4 Diatom identification results under different thresholds

表2 不同組織中不同閾值下模型識(shí)別的RR、PR 和F1 分?jǐn)?shù)Tab.2 The recall rate，precision rate，and F1 score of the identification model under different thresholds in different tissues （n=8，）

表2 不同組織中不同閾值下模型識(shí)別的RR、PR 和F1 分?jǐn)?shù)Tab.2 The recall rate，precision rate，and F1 score of the identification model under different thresholds in different tissues （n=8，）

在各閾值下，肺組織的PR均較高；肝、腎組織的PR均較低，假陽性率高。隨著閾值的增加，PR在各組織呈上升趨勢。

從各組織各閾值的RR均值和PR均值綜合計(jì)算所得的F1 分?jǐn)?shù)來看，肺組織高于肝、腎組織，模型在肺組織的F1 分?jǐn)?shù)隨著閾值升高而降低，在肝、腎組織的F1 分?jǐn)?shù)隨著閾值升高而升高。

3 討論

硅藻檢驗(yàn)是目前溺死診斷的金標(biāo)準(zhǔn)。近年來快速發(fā)展的微波消解-真空抽濾-掃描電子顯微鏡法有著一系列的優(yōu)勢：微波消解使雜質(zhì)消化得更為徹底，真空抽濾技術(shù)提高了硅藻的富集率，電子顯微鏡的高放大倍數(shù)使硅藻圖像分辨率大大提高，更利于分辨硅藻和進(jìn)行種屬分類，然而電子顯微鏡圖像數(shù)量龐大、陰性圖像多、人工識(shí)別耗時(shí)耗力。

為解決硅藻電子顯微鏡圖像識(shí)別問題，本研究訓(xùn)練了一個(gè)基于深度學(xué)習(xí)YOLOv3 模型的識(shí)別模型。YOLO 是REDMON 等[13]于2016 年提出的一個(gè)單階段目標(biāo)檢測模型，該模型將目標(biāo)檢測變?yōu)橐粋€(gè)回歸問題，比起傳統(tǒng)的二階段目標(biāo)檢測算法的先定位再分類，YOLO 僅經(jīng)過一個(gè)CNN 就能實(shí)現(xiàn)多個(gè)目標(biāo)的定位和識(shí)別，大大提升了識(shí)別速度，其檢測速度是R-CNN的1 000 倍、Fast R-CNN 的100 倍[9]。早期的YOLO 和YOLOv2 在識(shí)別精度上與二階段檢測器相比稍顯不足[14]，而2018 年推出的YOLOv3 大幅度提高了檢測的精度，并且在準(zhǔn)確率相當(dāng)?shù)那闆r下其檢測速度是其他模型的3～4 倍[15]。此外，該模型將圖像作為整體輸入特征提取網(wǎng)絡(luò)進(jìn)行訓(xùn)練，能夠?qū)W習(xí)到背景與目標(biāo)對(duì)象之間的更多信息，更適于復(fù)雜背景下小物體的目標(biāo)檢測。

本研究顯示，該模型在驗(yàn)證集和測試集上取得了94%以上的mAP和81%以上的平均RR，具有良好的性能。實(shí)際應(yīng)用中，在識(shí)別速度方面，該模型識(shí)別硅藻電子顯微鏡圖像的速度可達(dá)0.143 s/張，與人工識(shí)別相比顯著縮短了識(shí)別時(shí)間。在識(shí)別所用圖像數(shù)量和硅藻陽性圖像檢出量方面，肝、腎組織樣本由于硅藻數(shù)量少，需對(duì)所有掃描電子顯微鏡圖像進(jìn)行識(shí)別，識(shí)別模型檢出的硅藻陽性圖像數(shù)為300 張以下，顯著減少了需要人工復(fù)核的可能存在硅藻的圖像數(shù)；對(duì)于肺組織，人工識(shí)別和模型識(shí)別硅藻陽性圖像檢出量差異無統(tǒng)計(jì)學(xué)意義，但由于各樣本硅藻含量不同，不同案例肺組織檢出硅藻陽性圖像數(shù)有較大差異，因此識(shí)別肺組織硅藻所用圖像數(shù)量需要根據(jù)樣本實(shí)際硅藻含量決定，若樣本硅藻含量較少，則需要識(shí)別更多的圖片。

在識(shí)別RR方面，模型在0.4 的識(shí)別閾值下各組織的硅藻RR均較高，隨著識(shí)別閾值的增加，各組織的硅藻RR均有所下降。在識(shí)別PR方面，模型在0.4 的識(shí)別閾值下肺組織硅藻PR均值達(dá)87.8%，肝、腎組織PR均值不到5%，這可能是由于肝、腎組織實(shí)際含硅藻數(shù)量少，模型預(yù)測出的假陽性硅藻數(shù)量多，導(dǎo)致計(jì)算出的PR數(shù)值偏低。隨著閾值的增加，各組織硅藻PR均有所升高。

從F1 分?jǐn)?shù)看，模型在肺組織的識(shí)別效能顯著高于肝、腎組織，這可能是由于肝、腎組織硅藻含量較少，而肺組織硅藻含量較多，使得模型訓(xùn)練時(shí)所用的硅藻樣本圖像大多來自肺組織，導(dǎo)致模型識(shí)別肝、腎組織時(shí)預(yù)測的假陽性硅藻多，使模型在肝、腎組織的查準(zhǔn)率低、綜合識(shí)別效能低。F1 分?jǐn)?shù)顯示，模型在肺組織的識(shí)別效能隨著閾值增高而降低，結(jié)合RR和PR的變化來看，這可能是由于隨著閾值的升高，肺組織的PR較高，無太多上升空間，而肺組織RR下降明顯，使總體識(shí)別效能下降；模型在肝、腎組織的識(shí)別效能則隨著閾值增高而升高，這可能是由于肝、腎組織每張圖像上的硅藻數(shù)量少，模型預(yù)測每個(gè)目標(biāo)對(duì)象的置信度都較高，所以適當(dāng)提高閾值對(duì)真正硅藻的RR影響較小，并且能降低假陽性率，提高PR，取得更好的識(shí)別效能。

綜合以上應(yīng)用效果，筆者認(rèn)為，目前該模型應(yīng)用達(dá)到最佳效果需要注意以下幾點(diǎn)：（1）對(duì)于肺組織硅藻電子顯微鏡圖像識(shí)別，應(yīng)選擇較低的閾值（如0.4），在此閾值下硅藻RR和PR都能保持在較高水平。（2）對(duì)于肝、腎組織硅藻電子顯微鏡圖像識(shí)別，應(yīng)選擇中高閾值（如0.6、0.8），以在保持高RR的基礎(chǔ)上提高PR，取得更好的識(shí)別效能。（3）若肝、腎組織樣本硅藻含量極低，為避免漏檢，識(shí)別閾值設(shè)置為0.4較為適宜。

本研究構(gòu)建的模型能夠快速識(shí)別硅藻電子顯微鏡圖像并在0.4 閾值時(shí)取得各器官較高的硅藻RR和肺組織較高的硅藻PR，對(duì)于碎片硅藻也能很好地檢出。人工識(shí)別和模型識(shí)別肺組織硅藻陽性圖像檢出量差異無統(tǒng)計(jì)學(xué)意義，結(jié)合RR及PR可知，模型對(duì)于肺組織硅藻的泛化能力較好。雖然模型在肝、腎組織的硅藻PR不高，但由于肝、腎組織硅藻含量較低，電子顯微鏡圖像大多為無硅藻的圖片，模型能夠在盡量不遺漏肝、腎組織所有真正硅藻的情況下大量減少需要人工復(fù)核的圖片數(shù)量，使人工識(shí)別量從6 392 張降至300 張以下，很大程度上節(jié)省了識(shí)別時(shí)間。

然而，該模型也存在一些不足，值得改進(jìn)：（1）該模型對(duì)于肝、腎組織的PR不高，雖然已大大減少了待識(shí)別的圖像數(shù)量但仍需要硅藻檢驗(yàn)經(jīng)驗(yàn)豐富的專家進(jìn)行復(fù)核。這與模型訓(xùn)練使用的樣本圖像集中肺組織圖像較多有關(guān)，后續(xù)將加入大量肝、腎圖像樣本，使模型提高該類組織的PR，降低假陽性率，使沒有硅藻檢驗(yàn)經(jīng)驗(yàn)的人員也能使用該模型進(jìn)行識(shí)別。（2）目前該模型僅支持硅藻檢出，并不支持硅藻種屬分類識(shí)別。后期也將在訓(xùn)練圖像集中增加硅藻種屬類別標(biāo)簽進(jìn)行訓(xùn)練，嘗試?yán)迷撃Ｐ瓦M(jìn)行硅藻種屬自動(dòng)識(shí)別研究。

近年來，利用數(shù)字化方法識(shí)別硅藻越來越受人們關(guān)注，有多項(xiàng)研究已涉及基于深度學(xué)習(xí)的硅藻圖像識(shí)別。ZHOU 等[10]基于GoogLeNet 模型建立了強(qiáng)酸消化-光學(xué)顯微鏡檢驗(yàn)法的硅藻自動(dòng)識(shí)別模型，模型受試者操作特征（receiver operator characteristic，ROC）曲線的曲線下面積（area under the curve，AUC）達(dá)0.982，并在實(shí)際案例應(yīng)用中取得了92.45%的總體PR[11]。YU 等[12]利用RetinaNet 模型實(shí)現(xiàn)了硅藻電子顯微鏡圖像的自動(dòng)識(shí)別，模型在0.5 閾值下的PR為0.82、RR為0.88。在實(shí)際應(yīng)用中，YOLOv3 模型識(shí)別能夠達(dá)到較高的RR，并大幅降低了需要人工復(fù)核的硅藻陽性圖像數(shù)量，對(duì)于肺、肝、腎組織的完整硅藻和碎片硅藻都能有效檢出，在微波消解-真空抽濾-掃描電子顯微鏡法檢案中具有應(yīng)用可行性。

目前利用深度學(xué)習(xí)模型進(jìn)行硅藻種屬識(shí)別的研究也在不斷進(jìn)行[16-19]，相信未來基于深度學(xué)習(xí)模型進(jìn)行硅藻識(shí)別與分類的研究將會(huì)更加成熟與完善，該技術(shù)將逐漸應(yīng)用于法醫(yī)學(xué)實(shí)踐中提升硅藻檢驗(yàn)的效率。

綜上所述，本研究基于深度學(xué)習(xí)中的YOLOv3 模型構(gòu)建了一個(gè)硅藻電子顯微鏡圖像識(shí)別模型，模型在訓(xùn)練階段取得了良好的性能，并在實(shí)際案例應(yīng)用中展現(xiàn)了可行性。該模型能以極快的速度進(jìn)行識(shí)別，并在一定的閾值下取得良好的硅藻RR和肺組織硅藻PR，顯著降低了人工識(shí)別硅藻的工作量。未來計(jì)劃進(jìn)一步提高模型在肝、腎組織中的硅藻PR，完善硅藻種屬識(shí)別功能，擴(kuò)展該模型的應(yīng)用前景。