榮文雅,王仲霞,何婷婷,景 婧,王麗蘋,孫永強(qiáng),朱 云,許文濤,桑秀秀,田 淼,劉力政,冷源銘,楊 偉**,王睿林**
(1.南方醫(yī)科大學(xué)中醫(yī)藥學(xué)院 廣州 510515;2.解放軍總醫(yī)院第五醫(yī)學(xué)中心 北京 100039;3.復(fù)旦大學(xué)工程與應(yīng)用技術(shù)研究院 上海 200433;4.中國中醫(yī)科學(xué)院中醫(yī)臨床基礎(chǔ)醫(yī)學(xué)研究所 北京 100700)
腺病毒是一種DNA病毒,通常會(huì)引起輕微感染,涉及上呼吸道或下呼吸道、胃腸道或結(jié)膜[1],其具有發(fā)病急、感染途徑多、傳播速度快、病死率高和后遺癥發(fā)生率高、易合并肺炎產(chǎn)生重癥病例等特點(diǎn)[2],可引起嚴(yán)重的支氣管堵塞,也可通過血液傳播造成多系統(tǒng)功能損害,重癥肺炎還會(huì)影響電解質(zhì)與酸堿平衡等[3]。腺病毒感染??梢詫?dǎo)致重癥肺炎,尤其是3、7、55型血清型[4]。腺病毒感染所引起的急性呼吸道傳染病可能會(huì)在封閉人群中迅速傳播,擁擠和不良的衛(wèi)生行為可能會(huì)促進(jìn)傳播。近年來,我國關(guān)于腺病毒感染病例報(bào)道也逐漸增多[5-8]。
舌診作為中醫(yī)望診的重要內(nèi)容之一,可以為臨床辨證施治提供有效的依據(jù)[9]。傳統(tǒng)的舌診方法主要是根據(jù)醫(yī)生目測(cè)、語言文字說明、經(jīng)驗(yàn)辨析,并要求在短時(shí)間內(nèi)進(jìn)行判斷,其診斷結(jié)果既受醫(yī)生的知識(shí)水平、思維能力和診斷技能的限制,又受光線、溫度等外部客觀條件的影響,缺乏客觀評(píng)價(jià)標(biāo)準(zhǔn)[10]。而臨床舌診中對(duì)舌質(zhì)、舌苔、舌形與舌態(tài)特征的準(zhǔn)確定量詮釋尤為重要。運(yùn)用現(xiàn)代科學(xué)方法為舌診建立客觀指標(biāo),使其更準(zhǔn)確、客觀地反映人體機(jī)能狀態(tài),對(duì)于提高中醫(yī)舌診的客觀化和精準(zhǔn)度具有重要意義。目前舌象的研究已取得初步成效,如趙靜等[11]通過光譜技術(shù)來進(jìn)行舌象分析,用于疾病快速篩查和中醫(yī)證候辨識(shí)等。周勝等[12]應(yīng)用DS01-A道生四診儀檢測(cè)250例經(jīng)胃鏡診斷為胃潰瘍的舌象,發(fā)現(xiàn)觀察舌苔厚薄指數(shù)、腐膩指數(shù)、剝苔指數(shù)可以間接評(píng)估胃內(nèi)局灶性炎癥反應(yīng)的變化。洪樂等[13]通過舌象分析,綜合利用計(jì)算機(jī)圖像識(shí)別技術(shù)進(jìn)行智能化中醫(yī)體質(zhì)辨識(shí),建立了一個(gè)基于舌象分析的中醫(yī)體質(zhì)辨識(shí)系統(tǒng)。應(yīng)薦等[14]采用自主研發(fā)的PHE201體表紅外光譜儀對(duì)乳腺增生患者舌體各部進(jìn)行紅外輻射光譜測(cè)定,發(fā)現(xiàn)乳腺增生患者舌前、舌中、舌邊等部位紅外輻射光譜形態(tài)相似。
參照以往溫病論治[15],腺病毒感染患者的舌象具有一定特征,然而目前關(guān)于腺病毒感染舌診客觀化研究的報(bào)道較少。本研究以腺病毒感染患者為研究對(duì)象,采集并利用多維圖像分析技術(shù)分析患者舌象,應(yīng)用描述性統(tǒng)計(jì)、統(tǒng)計(jì)檢驗(yàn)和聚類分析討論其特征,使腺病毒中醫(yī)舌診更精準(zhǔn)化、客觀化,為中醫(yī)辨證論治提供支撐。
收集了85例2020年6月-2021年1月就診于解放軍總醫(yī)院第五醫(yī)學(xué)中心肝病醫(yī)學(xué)部中醫(yī)肝病科已確診腺病毒感染患者的舌診圖像:其中腺病毒患者發(fā)熱(A+C)41例、無發(fā)熱(T+P)30例、密切接觸(N)14例(缺1例患者舌診圖像),選取30例與患者年齡、身高及體重相匹配的成年人舌診圖像作為正常對(duì)照。從以上受試者中共采集280張舌象信息,每位受試者最多的采集了舌象7張,最少的采集了舌象1張。由于圖片質(zhì)量問題,共刪除21張舌象,其中,有11張舌象無法識(shí)別(舌象特征取值全為0),有10張舌象識(shí)別不清(多數(shù)舌象特征取值為0),最終納入分析的感染腺病毒患者66例和正常人群30例。然后,對(duì)受試者單日多次拍照的舌象特征數(shù)據(jù)取均值。其中,特征變量“是否檢測(cè)到齒痕”取均值后向上取整。
腺病毒患者均符合2013年頒布的《腺病毒感染診療指南》中腺病毒診斷標(biāo)準(zhǔn)[16]:①咽拭子實(shí)時(shí)定量PCR(Real-time PCR)法檢測(cè)腺病毒特異性核酸陽性;②血清腺病毒特異性IgM抗體陽性;③急性期與恢復(fù)期雙份血清標(biāo)本腺病毒特異性IgG抗體4倍以上升高。
①腺病毒患者:符合2013年頒布的《腺病毒感染診療指南》中腺病毒診斷標(biāo)準(zhǔn)[16];年齡18-80歲,性別不限;愿意配合舌象拍攝并簽署知情同意書者?;颊吲R床資料完整且依從性高;②正常人群:既往無慢性疾病史;3個(gè)月內(nèi)無急性疾病史等;舌象為淡紅舌、薄白苔者;年齡18-70歲,性別不限;愿意配合舌拍攝并簽署知情同意書者。
①因精神因素對(duì)病情描述不清、不能配合采集舌診圖像者;②嚴(yán)重水電解質(zhì)紊亂、酸中毒等急性并發(fā)癥者;③患有其他嚴(yán)重內(nèi)科疾病如腫瘤、免疫系統(tǒng)、血液系統(tǒng)等疾病者;④服用類固醇等影響糖代謝的藥物者;⑤孕婦及哺乳期患者。患兒監(jiān)護(hù)人同意且簽署知情同意書,本研究經(jīng)過醫(yī)院倫理委員會(huì)審批同意實(shí)施(倫理審批編號(hào):2020001D)。
2.1.1 采集環(huán)境與設(shè)備
舌苔圖像采集時(shí)間為早晨的7∶00-9∶00。舌象采集環(huán)境為光線為良好的門診室內(nèi),色溫范圍:4500-7200 K,照度范圍:320-3630 lx。被測(cè)者拍攝前禁食食物,禁飲有色飲料。被檢查者坐位,張嘴伸舌、舌體放松、舌面平展、舌尖自然下垂,色差校正色卡舉在相機(jī)視線內(nèi),拍攝舌圖像。相機(jī)與受試者的拍攝距離為20 cm如需重新拍攝,休息3-5 min后再采集,最后選取最佳者儲(chǔ)存。舌診圖像數(shù)字分析診斷系統(tǒng)用于所有舌像的采集。
自然光色溫照度測(cè)量采用日本柯尼卡美能達(dá)公司的CL-200色溫照度計(jì);舌象采集設(shè)備選用Canon EOS M100數(shù)碼相機(jī),固定于舌象采集專用支架上;校正色卡選用柯達(dá)Q-13灰階卡,拍攝時(shí)放于舌側(cè)。用Canon EOS M100數(shù)碼相機(jī)在自動(dòng)白平衡模式下拍攝舌象以及灰度色卡。選取柯達(dá)Q-13灰階卡的中灰色塊來設(shè)置白平衡,通過Nikon Capture NX軟件進(jìn)行顏色的校正。其原理主要是:軟件根據(jù)數(shù)碼相機(jī)拍攝RAW格式圖像時(shí)記錄的白平衡信息將圖像中的顏色通道通過增益系數(shù)進(jìn)行校正。
2.1.2 舌診圖像目標(biāo)檢測(cè)
對(duì)于舌診圖像目標(biāo)檢測(cè),主要選取了4個(gè)關(guān)鍵目標(biāo)的檢測(cè),包括舌頭區(qū)域、舌苔區(qū)域、齒痕區(qū)域以及裂紋區(qū)域,這里采用了YOLOV4[17]網(wǎng)絡(luò)框架下的舌診圖像目標(biāo)檢測(cè)算法集合,通過對(duì)最新算法的大量實(shí)驗(yàn),得到最優(yōu)的算法組合,提高了卷積神經(jīng)網(wǎng)絡(luò)(CNN)的準(zhǔn)確性。
2.1.3 儀器系統(tǒng)處理
(1)舌苔顏色提?。悍謩e計(jì)算舌苔區(qū)域的RGB分量的顏色均值。R代表紅色值,G代表綠色值,B代表藍(lán)色值。
(2)齒痕判別:齒痕為目標(biāo)檢測(cè)的結(jié)果,判斷齒痕的個(gè)數(shù),置信度取0.6。
(3)舌苔裂紋占比計(jì)算:根據(jù)舌苔上裂紋檢測(cè)的結(jié)果做分割,同時(shí)對(duì)于多個(gè)細(xì)小的裂紋交錯(cuò)的部分做區(qū)域融合,得到裂紋的面積,計(jì)算了裂紋與舌苔的面積占比。
(4)舌苔濕度計(jì)算:舌苔濕度定義為反光區(qū)域,即舌苔部分亮度超過一定閾值的白色點(diǎn),根據(jù)現(xiàn)有提供的圖片樣本,很多舌苔上反光的部分并不是很亮,這里給出的閾值為170,統(tǒng)計(jì)反光區(qū)域的亮度點(diǎn),計(jì)算占比。由于拍攝圖像某些反光并不明顯或圖像偏暗,會(huì)使得某些值偏低。后續(xù)可以考慮統(tǒng)一拍攝環(huán)境及補(bǔ)光標(biāo)準(zhǔn)。
(5)舌苔紋理特征值提取:對(duì)于舌苔部分計(jì)算Tamura紋理特征,Tamura紋理特征的六個(gè)分量對(duì)應(yīng)于心理學(xué)角度上紋理特征的六種屬性分別是粗糙度(Coarseness)、對(duì)比度(Contrast)、方向度(Directionality)、線像度(Linelikeness)、規(guī)整度(Regularity)和粗略度(Roughness)。這里主要選取了粗糙度、對(duì)比度、方向度3分量特征來描述舌苔的紋理特征。
將所有患者基本信息、中醫(yī)舌象特征信息導(dǎo)入數(shù)據(jù)庫,采用R4.0.2和SPSS22.0進(jìn)行統(tǒng)計(jì)分析與驗(yàn)證。各組數(shù)據(jù)均進(jìn)行正態(tài)性檢驗(yàn)和方差齊性檢驗(yàn),計(jì)數(shù)資料采用頻數(shù)表示,使用卡方檢驗(yàn)或Fisher精確檢驗(yàn)進(jìn)行分析。以腺病毒感染者不同程度無發(fā)熱、密切接觸、發(fā)熱和正常為結(jié)局,隨機(jī)森林算法[18]采用Gini指數(shù)確定舌象特征的重要度,模糊K-Means聚類分析采用系統(tǒng)聚類中的R型分析法;計(jì)量資料符合正態(tài)分布的數(shù)據(jù)用均值±標(biāo)準(zhǔn)差(±s)表示,采用t檢驗(yàn)或者方差分析(ANOVA)比較組或多組間的差異,非正態(tài)分布的數(shù)據(jù),用中位數(shù)(M)、四分位間距表示,采用Mann-Whitney U檢驗(yàn)或Kruskall-Wallis H檢驗(yàn)來比較兩組或多組間的差異。以P〈0.05為差異有統(tǒng)計(jì)學(xué)意義。
腺病毒感染不同程度患者與正常人群性別(男69例,女27例)、年齡(23.0(20.0,25.0)歲,最小值18歲,最大值80歲)、身高(172.0±6.9 cm)、體重(65.0±7.3 kg)、基線體溫(36.5(36.5,39.0)℃,最小值36.5℃,最大值40.5℃)等基線情況比較,差異無統(tǒng)計(jì)學(xué)意義(P〉0.05),具有可比性(見表1)。
腺病毒感染不同程度患者與正常人群舌色、舌形、苔色、苔質(zhì)特征的整體比較見表1。處理后的96例舌象數(shù)據(jù)按嚴(yán)重程度分類,針對(duì)每個(gè)變量對(duì)4類人群分析進(jìn)行兩兩比較,結(jié)果如表2所示;其中R代表紅色值,G代表綠色值,B代表藍(lán)色值。四類人群中舌苔顏色R,舌苔顏色G,舌苔顏色B的比較,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05);舌苔顏色R中,發(fā)熱患者與密切接觸患者,密切接觸患者與正常人群,無發(fā)熱患者與正常人群,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05)。舌苔顏色G中,密切接觸患者與正常人群,無發(fā)熱患者與正常人群,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05);舌苔顏色B中,發(fā)熱患者與無發(fā)熱患者,無發(fā)熱患者與正常人群,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05);齒痕有無、舌苔濕度、紋理粗糙度、紋理方向、紋理對(duì)比度的4類人群比較,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05)。其中,齒痕有無檢測(cè)中,發(fā)熱患者與正常人群,密切接觸患者與無發(fā)熱患者,無發(fā)熱患者與正常人群,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05);舌苔濕度中,發(fā)熱患者、密切接觸患者、無發(fā)熱患者與正常人群,差異都具有統(tǒng)計(jì)學(xué)意義(P〈0.05);紋理粗糙度中,發(fā)熱患者、密切接觸患者、無發(fā)熱患者與正常人群,差異都具有統(tǒng)計(jì)學(xué)意義(P〈0.05);紋理方向度中,發(fā)熱患者與正常人群,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05);紋理對(duì)比度中,發(fā)熱患者、密切接觸患者、無發(fā)熱患者與正常人組,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05)。4類人群的典型舌苔圖像見圖1。
表1 舌象特征的組間比較
續(xù)表
圖1 物體探測(cè)儀
以腺病毒感染不同程度:無發(fā)熱、密切接觸、發(fā)熱以及正常人群為結(jié)局,利用隨機(jī)森林算法對(duì)所有舌象特征進(jìn)行建模,用Gini指數(shù)篩選出重要特征:舌苔顏色R均值、舌苔顏色G均值、紋理粗糙度、紋理對(duì)比度,結(jié)合3.2的結(jié)果,納入以上舌象特征,對(duì)感染腺病毒不同程度患者與正常人群進(jìn)行模糊K-Means聚類R型分析,得到如下表3和圖2的聚類結(jié)果。首先每一類都是有代表性的,最具特點(diǎn)的第4類是以正常人為主,結(jié)合圖2,可以看出正常人更多的體現(xiàn)在第2個(gè)主成分上,進(jìn)而可以大致推斷:較正常人群,患病人群在第2主成分上表現(xiàn)稍差;第3類的代表性樣本為發(fā)熱(A+C)和無發(fā)熱(T+P)為主,主要表現(xiàn)在第2主成分上;第2類的代表性樣本以發(fā)熱(A+C)人群更為明顯,其在第1主成分上表現(xiàn)的更好;第1類的代表性樣本分布沒有較大差異,基本上不偏向哪個(gè)主成分。綜上,根據(jù)舌苔顏色R均值、舌苔顏色G均值、紋理粗糙度、紋理對(duì)比度獲取的兩個(gè)主成分,用不同顏色圓實(shí)點(diǎn)畫出的4類人群散點(diǎn)圖具有較好的區(qū)分度。
圖2 4類人群典型舌苔圖像
表3 腺病毒感染的3類不同程度患者與正常人群舌象特征進(jìn)行模糊K-Means聚類R型分析
另外,我們還按舌苔顏色R均值、舌苔顏色G均值、舌苔顏色B均值、齒痕有無、舌苔裂紋占比、舌苔濕度、舌苔紋理(紋理粗糙度,紋理方向度,紋理對(duì)比度)分別畫出4類人群對(duì)應(yīng)的核密度分布圖,從圖3可以看出,舌苔顏色R均值、舌苔顏色G均值、舌苔顏色B均值、紋理對(duì)比度、舌苔濕度、紋理粗糙度對(duì)4類人群具有一定的區(qū)分度。
圖3 腺病毒感染的3類不同程度患者與正常人群舌象特征進(jìn)行模糊K-Means聚類R型分析
舌診作為中醫(yī)臨床診療的重要要素,是一種極具中醫(yī)特色的無創(chuàng)性診斷方法,中醫(yī)溫病學(xué)強(qiáng)調(diào)注重舌象的觀察,有“溫病看苔”的表述[19],腺病毒感染的臨床診療需要參照溫病學(xué)內(nèi)容[20],其舌象分析可以作為溫病舌象觀察的一個(gè)示例。溫病容易化燥傷陰,在溫病后期,邪入營血或深入下焦。則邪熱熏蒸,陰液耗傷嚴(yán)重。溫病變化發(fā)展迅速,而舌象對(duì)于疾病的反映尤為敏感,能夠反映人體氣血津液的盈虧情況。因此通過辨舌可以直觀的判斷津液耗傷的嚴(yán)重程度。溫病對(duì)舌象的觀察包括舌苔、舌質(zhì)、舌態(tài)等內(nèi)容[21]。溫病辨舌可以辨病邪性質(zhì)、病位深淺、病勢(shì)進(jìn)退、津液盈虧等,并且舌象反映的更加直觀便捷,對(duì)于起病急,傳變快的溫病診斷具有重要意義。中醫(yī)舌診的觀察,歷來主要以臨床中醫(yī)師觀察后進(jìn)行語言描述進(jìn)行記錄,舌質(zhì)顏色的紅、絳,舌體的瘦、薄,舌苔的厚、膩,這些表述在不同的中醫(yī)師的觀察往往只能以定性的方式呈現(xiàn),無法加以量化,隨著近代色度學(xué),近代光學(xué)技術(shù),數(shù)字圖像處理技術(shù)和計(jì)算機(jī)處理技術(shù)的發(fā)展,使中醫(yī)舌診從主觀定性的表達(dá)走向客觀化、定量化、標(biāo)準(zhǔn)化成為可能。本研究通過現(xiàn)代圖像分析技術(shù)從舌苔顏色、齒痕有無、舌苔濕度、紋理粗糙度、紋理方向度、紋理對(duì)比度幾個(gè)維度對(duì)腺病毒感染者中醫(yī)舌象觀察內(nèi)容進(jìn)行客觀表述,發(fā)現(xiàn)舌苔顏色R方面,發(fā)熱患者與密切接觸患者、密切接觸患者與正常人群、無發(fā)熱患者與正常人群,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05)。舌苔顏色G方面,密切接觸患者與正常人群、無發(fā)熱患者與正常人群,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05);舌苔顏色B方面,發(fā)熱患者與無發(fā)熱患者、無發(fā)熱患者與正常人群,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05);以上結(jié)果表明:正常人群感染病毒到產(chǎn)生發(fā)熱癥狀的過程中,舌苔顏色發(fā)生了客觀的變化,與中醫(yī)溫病學(xué)認(rèn)為溫病進(jìn)展過程中,舌苔從白到黃,溫?zé)岬牟〕烫攸c(diǎn)是一致的。正常舌象舌質(zhì)淡紅,舌苔薄白;黃苔一般是從白苔轉(zhuǎn)化而來,標(biāo)志邪熱已入氣分,而其厚薄、潤燥又可反映病位的深淺及津傷的程度。《溫?zé)嵴摗氛f:“或黃甚,或如沉香色,或如灰黃色,或老黃色,或中有斷紋,皆當(dāng)下之,如小承氣湯,用檳榔、青皮、枳實(shí)、元明粉、生首烏等?!睖夭∫婞S苔時(shí)應(yīng)當(dāng)以祛邪存陰為主要治療原則。另外,齒痕有無、舌苔濕度、紋理粗糙度、紋理方向、紋理對(duì)比度的四類人群比較,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05)。其中,齒痕有無檢測(cè)方面,發(fā)熱患者與正常人群、密切接觸患者與無發(fā)熱患者、無發(fā)熱患者與正常人群,差異均具有統(tǒng)計(jì)學(xué)意義(P〈0.05)。舌苔濕度方面,發(fā)熱患者與正常人群、密切接觸患者與正常人群、無發(fā)熱患者與正常人群,差異均具有統(tǒng)計(jì)學(xué)意義(P〈0.05);紋理粗糙度方面,發(fā)熱患者與正常人群、密切接觸患者與正常人群、無發(fā)熱患者與正常人群,差異均具有統(tǒng)計(jì)學(xué)意義(P〈0.05);紋理方向度方面,發(fā)熱患者與正常人群,差異具有統(tǒng)計(jì)學(xué)意義(P〈0.05);紋理對(duì)比度方面,發(fā)熱患者與正常人群、密切接觸患者與正常人群、無發(fā)熱患者與正常人群,差異均具有統(tǒng)計(jì)學(xué)意(P〈0.05),將上述的幾個(gè)特征和隨機(jī)森林算法、模糊KMeans聚類分析的圖像分析結(jié)果可與中醫(yī)理論可進(jìn)行互釋:中醫(yī)認(rèn)為齒痕主濕[22],舌苔濕度在中醫(yī)學(xué)中的表述是“苔膩”程度與“濕”,舌苔之膩者,乃為穢濁邪氣入里太深,正氣不能祛邪于外而上泛于口所致,舌苔黃膩而垢,或干澀深黃厚膩是表邪入里之候。《辨舌指南·辨舌之津液》記載:“膩者有形,揩之不去,為穢濁盤踞中宮?!?組人群舌象特征之間的差異增加了濕的客觀精準(zhǔn)辨識(shí);舌苔裂紋的形成在中醫(yī)學(xué)上也常常提示有熱邪傷陰。本研究對(duì)腺病毒不同程度感染者的舌色、齒痕及舌苔紋理、舌苔濕度等方面的動(dòng)態(tài)變化進(jìn)行分析,將主觀的文字表述與客觀的圖像分析結(jié)果進(jìn)行互參互釋,也為腺病毒感染中醫(yī)證候病機(jī)提供客觀依據(jù),使以腺病毒為代表的溫病舌苔診療更加客觀量化,有助于發(fā)揮中醫(yī)藥特色和優(yōu)勢(shì)的診療技術(shù)在腺病毒的防治方面的作用,同時(shí)也有助于對(duì)腺病毒等溫病的臨床療效進(jìn)行客觀評(píng)價(jià),進(jìn)一步的研究可以將不同病毒或細(xì)菌感染屬于中醫(yī)溫病范疇的舌像特征進(jìn)行綜合對(duì)比研究,探討不同溫病之間舌像客觀的差異,為中醫(yī)藥防治新突發(fā)傳染病提供更細(xì)化的技術(shù)支撐。
圖4 腺病毒感染的三類不同程度患者與正常人群的舌象特征核密度分布圖