基于灰度插值直方圖線性化的醫(yī)學(xué)MR圖像處理設(shè)計與實現(xiàn)

摘 要：為了解決醫(yī)學(xué)MR圖像邊緣模糊、亮度不合理和噪聲抑制效果不好等問題，提出了一種灰度插值直方圖線性化算法。該方法先對含混合噪聲的醫(yī)學(xué)MR圖像進行直方圖線性化處理，抑制部分噪聲，使亮度達到合理的范圍，進一步將得到的結(jié)果進行多項式插值運算，消除象素塊間的灰度差異，使圖像達到層次清晰。仿真實驗結(jié)果表明，該算法獲得的圖像亮度、峰值信噪比的值和結(jié)構(gòu)相似性指數(shù)指標(biāo)均最大，其中，峰值信噪比值高出其他方法約為1.6～3.2 dB，結(jié)構(gòu)相似性指數(shù)指標(biāo)高出其他方法約為5%～7%。該方法可有效地降低噪聲，較好地保持了醫(yī)學(xué)MR圖像邊緣和細節(jié)信息，其效果明顯優(yōu)于自適應(yīng)直方圖均衡化和直方圖局部線性化等方法。

關(guān)鍵詞：灰度插值；直方圖；醫(yī)學(xué)MR圖像；線性化

中圖分類號：TP391.4

文獻標(biāo)志碼：A

醫(yī)學(xué)影像圖像線性化是進行圖像分割、特征提取和分類識別的核心［1-2］。通過醫(yī)學(xué)知識與數(shù)學(xué)算法對醫(yī)學(xué)影像圖像進行定量分析和描述解釋使臨床醫(yī)生更好地觀測病變組織和更準(zhǔn)確地診斷疾?。?］，為可靠的治療方案提供證據(jù)，達到精準(zhǔn)治療，對臨床診斷、治療及預(yù)后評估具有重要的意義［4］。許紀(jì)亞等［5］利用基于Laplacian金字塔的二維直方圖均衡化方法對醫(yī)學(xué)CT圖像和MRI進行了對比度增強研究，結(jié)果表明該方法在圖像細節(jié)增強方面有較好的效果；李艷萍等［6］研究了基于限峰分離模糊直方圖均衡化的圖像增強算法，該方法在增強圖像對比度和保留圖像亮度方面比較理想；苑豪杰等［7］報道了自適應(yīng)多子直方圖均衡算法能夠?qū)D像合理增強、抑制過增強的不利，在保持圖像細節(jié)方面有比較好的效果，但易受噪聲影響；王超等［8］探討了自適應(yīng)直方圖均衡化算法對增強有高峰圖像的問題，研究顯示該方法對圖像細節(jié)增強較理想；陳博洋［9］采用局部線性化直方圖對低能彩色遙感圖像處理研究顯示，該方法能提高圖像的可視性能和目標(biāo)識辨率；VAKIL等［10］通過比較顱內(nèi)動脈粥樣硬化斑塊的形態(tài)學(xué)及直方圖紋理特征對癥狀性斑塊特性進行了研究，結(jié)果顯示有較好的區(qū)別度；蘇葉等［11］通過灰度直方圖均衡對由于受標(biāo)尺、偽影、曝光、噪聲等影響的手骨X光片進行預(yù)處理，提高了骨齡預(yù)測正確率；文海瓊等［12］通過多閾值直方圖均衡化方法對低照度圖像進行了增強研究顯示，該方法能增強圖像的對比度和圖像的視覺感受性，但不能較好地保持圖像的邊緣信息；胡小琴等［13］通過灰度直方圖的均值、方差等紋理特性，結(jié)合ROC曲線的簡易聯(lián)合診斷模型對亞實性結(jié)節(jié)的肺腺癌進行研究顯示，該方法能提高肺腺癌病理亞型診斷率；徐茂林等［14］采用乳腺浸潤性導(dǎo)管癌超聲灰度直方圖特性對人表皮生長因子受體-2表達研究顯示，該方法可提高乳腺浸潤性導(dǎo)管癌的識別度，但對含噪聲圖像不理想；管秀紅［15］通過直方圖和多模態(tài)MRI方法對多發(fā)膠質(zhì)瘤影像異質(zhì)性與腫瘤生物學(xué)行為研究顯示，該方法能提高多發(fā)膠質(zhì)瘤診斷率，但對含噪聲圖像的靈敏度低；楊婷婷等［16］利用T2WI灰度直方圖的均值、變異度、偏度、峰度等特征，對原發(fā)性中樞神經(jīng)系統(tǒng)淋巴瘤和腦膠質(zhì)瘤進行了研究，結(jié)果顯示該方法有較高的鑒別度，但對含噪聲圖像不理想；陳倩倩等［17］通過增強T1WI全域灰度直方圖的均值、變異度、峰度、偏度絕對值等參數(shù)和ROC曲線對生殖細胞瘤進行了分析，結(jié)果表明該方法有較高的鑒別功效；潘霞等［18］利用ADC直方圖軟件和Spearman對腫瘤輪廓相應(yīng)的直方圖參數(shù)與原發(fā)腫瘤分期、腫瘤區(qū)域淋巴結(jié)的分期關(guān)系研究表明，該方法有較高的鑒別度；白國杰等［19］利用統(tǒng)計學(xué)方法結(jié)合MaZda軟件將雙參數(shù)磁共振成像灰度直方圖應(yīng)用于移行帶前列腺癌與良性前列腺增生的研究表明，該方法提高了前列腺癌與良性前列腺增生的鑒別度；?？∠嫉龋?0］采用灰度直方圖方法對AT/RT和MB的差異性進行研究顯示，該方法有較高的區(qū)別度，但易受圖像中噪聲的影響。

基于上述研究，本文利用灰度插值直方圖線性化算法對醫(yī)學(xué)MR圖像邊緣模糊、亮度不合理和噪聲抑制效果不好等問題進行了探討，并采用Matlab平臺設(shè)計實現(xiàn)了本算法的醫(yī)學(xué)影像圖像處理系統(tǒng)，驗證了本文方法的有效性。對醫(yī)學(xué)MR圖像噪聲抑制、邊緣清晰度的提高、細節(jié)信息保護及亮度合理均衡有實際指導(dǎo)意義。

1 基于灰度插值直方圖線性化的醫(yī)學(xué)MR圖像的數(shù)學(xué)模型分析

醫(yī)學(xué)影像圖像的灰度直方圖線性化是利用灰度變換函數(shù)將原影像圖像的直方圖映射成灰度均勻分布的直方圖。通過函數(shù)灰度點運算和插值使醫(yī)學(xué)影像圖像的直方圖均勻分布，其信息熵達到最大，提高了圖像的清晰度。直方圖能給出圖像的灰度范圍、每個灰度級的頻數(shù)、灰度的分布、整幅圖像的亮度和平均明暗對比度等。每個灰度級出現(xiàn)的概率可表示為［21］

Pr（rk）＝Nk/N，k=0，1，2，…，L－1（1）

式中：Pr（rk）表示第k個灰度級出現(xiàn)的概率；Nk為第k個灰度級出現(xiàn)的頻數(shù)；N為圖像像素總數(shù)；L為圖像中可能的灰度級總數(shù)。由此可得圖像的灰度累積分布函數(shù)sk為

式中：sk為歸一化灰度級。通過直方圖能夠觀察到醫(yī)學(xué)影像圖像整體的灰度數(shù)、灰度級等統(tǒng)計特性。

為了改善圖像的明暗對比度［22］，可用公式（3）表示。

為進一步消除象素塊間的灰度差異，使圖像層次清晰，采用多項式插值法，其公式為

（3）、（4）式中B 為圖像在該處的亮度，α為多項式插值參數(shù)，其變化與圖像亮度關(guān)系如圖1所示，這里α取值為0.6，使亮度達到最佳。

2 灰度插值直方圖線性化在醫(yī)學(xué)MR圖像中的實現(xiàn)

2.1 實驗評價

使用峰值信噪比和結(jié)構(gòu)相似性指數(shù)指標(biāo)對本文算法抑制噪聲性能高低進行評估［23］。

2.2 實驗流程

實驗采用Matlab程序語言，進行編程實現(xiàn)算法仿真。圖2所示為基于灰度插值直方圖線性化算法的醫(yī)學(xué)影像圖像的處理流程。先讀入原醫(yī)學(xué)影像圖像數(shù)據(jù)，接著統(tǒng)計醫(yī)學(xué)影像圖像灰度值，并計算每級灰度概率，存入緩沖區(qū)，在循環(huán)變量的控制下對直方圖均衡化達到符合的要求后輸出線性化處理后的醫(yī)學(xué)影像圖像。

2.3 實驗結(jié)果分析

該實驗中選用灰暗的原始醫(yī)學(xué)MR圖像如圖3所示，將自適應(yīng)直方圖均衡化、直方圖局部線性化和本文的方法進行了比較研究，驗證了本文方法的有效性。圖4（a）為在原始醫(yī)學(xué)MR圖3中加有零均值、方差為0.30、噪聲密度為0.30的混合噪聲圖。圖4（b）為圖4（a）對應(yīng)的統(tǒng)計直方圖；圖5（a）是圖4（a）經(jīng)過自適應(yīng)直方圖均衡化處理后的圖，圖5（b）為圖5（a）對應(yīng)的統(tǒng)計直方圖；圖6（a）是圖4（a）經(jīng)過直方圖局部線性化處理后的圖，圖6（b）為圖6（a）對應(yīng)的統(tǒng)計直方圖；圖7（a）是圖4（a）經(jīng)過本文算法處理后的圖，圖7（b）為圖7（a）對應(yīng)的統(tǒng)計直方圖。由圖4（a）、（b） 可觀察到：含混合噪聲的原始醫(yī)學(xué)MR圖像整個畫面較暗、模糊不清，灰度值幾乎分布在100以下，其中，灰度值在60以內(nèi)的約占80%，其占有的范圍很小，致使醫(yī)學(xué)MR圖像暗淡、對比度低、層次感不強。通過圖4（a）、圖5（a）、圖6（a）、圖7（a）可以觀察到，自適應(yīng)直方圖均衡化方法雖然一定程度上提高了醫(yī)學(xué)MR圖像的對比度，但對噪聲的抑制效果不理想，圖像存在朦朧感，邊緣不清晰；直方圖局部線性化方法對噪聲抑制較為理想，但圖像存在細節(jié)不清晰、邊緣模糊、對比度不高；本文方法處理后的醫(yī)學(xué) MR圖像變得明亮、層次感得到加強，細節(jié)也明顯變地清晰可辨，圖像中的噪聲得到有效抑制，邊緣相似性得到很好保持，視覺效果明顯比前兩種方法的好。從其對應(yīng)的直方圖4（b）、圖5（b） 、圖6（b）和圖7（b）可看出，本文方法處理之后的圖像灰度值的區(qū)間得以拓寬。

為了更進一步考查三種方法對含混合噪聲的醫(yī)學(xué)MR圖像處理性能的優(yōu)劣，運用亮度、峰值信噪比和結(jié)構(gòu)相似性指數(shù)衡量其性能，表1所列為其處理結(jié)果。由表1可見，本文提出的方法處理后圖像的亮度得到增強，其值比自適應(yīng)直方圖均衡化方法增大了40個單位，比直方圖局部線性化方法增大了7個單位；本文方法處理后圖像的峰值信噪比值最大，其值比自適應(yīng)直方圖均衡化方法高3.2 dB，比直方圖局部線性化方法高1.6 dB；本文方法處理圖像后的結(jié)構(gòu)相似性指數(shù)得以較大的提高，其值比自適應(yīng)直方圖均衡化方法提高7%，比直方圖局部線性化方法提高5%。這表明本文方法處理之后的醫(yī)學(xué)影像圖像既明顯地改善了顯示效果，又最大限度地保持了原始MR圖像的特征?？梢娫摲椒▽︶t(yī)學(xué)影像圖像中主要像素的灰度值得以加寬，而次要像素的灰度值得到壓縮，拓寬了醫(yī)學(xué)影像圖像灰度的動態(tài)范圍，有效抑制了混合噪聲和提高了結(jié)構(gòu)相似性指數(shù)，其效果明顯優(yōu)于自適應(yīng)直方圖均衡化和直方圖局部線性化方法。

3 結(jié)論

論文較詳細地討論了醫(yī)學(xué)影像圖像灰度插值直方圖線性化算法的模型。采用Matlab平臺設(shè)計實現(xiàn)了基于灰度插值直方圖線性化的醫(yī)學(xué)MR圖像處理系統(tǒng)。將該系統(tǒng)進行了仿真實驗，并對醫(yī)學(xué)MR圖像的亮度、峰值信噪比和結(jié)構(gòu)相似性指數(shù)等進行了分析，研究表明本文方法拓寬了醫(yī)學(xué)MR圖像灰度的動態(tài)范圍，有效抑制了混合噪聲和提高了結(jié)構(gòu)相似性指數(shù)，增加了醫(yī)學(xué)圖像的亮度、對比度、清晰度。為醫(yī)學(xué)圖像處理的實際應(yīng)用和理論研究提供借鑒。

參考文獻：

[1]黃渝萍， 李偉生. 醫(yī)學(xué)圖像融合方法綜述［J］. 中國圖象圖形學(xué)報， 2023， 28（1）： 118-143.

［2］ 陳軍. 基于小波分析的醫(yī)學(xué)影像圖像除噪設(shè)計及仿真研究［J］. 貴州大學(xué)學(xué)報（自然科學(xué)版）， 2020， 37（5）： 78-81， 88.

［3］ 楊宏宇， 楊帆. 基于圖像去噪和圖像生成的對抗樣本檢測方法［J］. 湖南大學(xué)學(xué)報（自然科學(xué)版）， 2023， 50（8）： 72-81.

［4］ THANH D N H， PRASATH V B S， SON N V， et al. An adaptive image inpainting method based on the modified mumford-shah model and multiscale parameter estimation［J］. Computer Optics， 2019， 43（2）： 251-257.

［5］ 許紀(jì)亞， 孫佳寧， 喬雙. 基于多尺度二維直方圖均衡化的醫(yī)學(xué)圖像增強方法［J］. 東北師大學(xué)報（自然科學(xué)版）， 2020， 52（1）： 88-91.

［6］ 李艷萍， 寧躍飛， 楊偉. 基于限峰分離模糊直方圖均衡化的圖像增強算法［J］. 西南師范大學(xué)學(xué)報（自然科學(xué)版）， 2018， 43（3）： 42-49.

［7］ 苑豪杰， 劉昌祿， 許建平， 等. 基于自適應(yīng)多子直方圖均衡的圖像增強算法［J］. 指揮控制與仿真， 2017， 39（5）： 45-49.

［8］ 王超， 孫玉秋， 徐石瑤， 等. 自適應(yīng)直方圖均衡化圖像增強算法研究［J］. 長江大學(xué)學(xué)報（自然科學(xué)版）， 2018， 15（1）： 55-59.

［9］ 陳博洋. 彩色遙感圖像的亮度直方圖局部線性化增強［J］. 光學(xué)精密工程， 2017， 25（2）： 502-508.

［10］VAKIL P， ELMOKA A H， SYED F H， et al. Quantifying intracranial plaque permeability with dynamic contrast-enhanced MRI： a pilot study［J］. American Journal of Neuroradiology， 2017， 38（2）： 243-249.

［11］蘇葉， 李婧， 徐寅林. 手骨X光片骨齡預(yù)測中圖像預(yù)處理的研究［J］. 南京師范大學(xué)學(xué)報（工程技術(shù)版）， 2021， 21（2）： 54-59.

［12］文海瓊， 李建成. 基于直方圖均衡化的自適應(yīng)閾值圖像增強算法［J］. 中國集成電路， 2022， 31（3）： 38-42， 71.

［13］胡小琴， 黃原義， 萬兵， 等. CT灰度直方圖在鑒別小于20 mm亞實性結(jié)節(jié)的肺腺癌病理亞型的價值［J］. 臨床放射學(xué)雜志， 2021， 40（4）： 712-716.

［14］徐茂林， 李芳， 曾書娥， 等. 乳腺浸潤性導(dǎo)管癌超聲灰度直方圖特征與HER-2表達的相關(guān)性［J］. 中國醫(yī)學(xué)影像學(xué)雜志， 2022， 30（12）： 1224-1229.

［15］管秀紅. 多模態(tài)MRI技術(shù)聯(lián)合直方圖分析多發(fā)膠質(zhì)瘤影像異質(zhì)性與腫瘤生物學(xué)行為的關(guān)系［J］. 中國CT和MRI雜志， 2023， 21（1）： 22-24.

［16］楊婷婷， 占鳴， 謝春梅， 等. T2WI灰度直方圖定量分析在原發(fā)性中樞神經(jīng)系統(tǒng)淋巴瘤和腦膠質(zhì)瘤鑒別診斷中的價值［J］. 醫(yī)學(xué)信息， 2021， 34（13）： 4-6.

［17］陳倩倩， 張焱， 程敬亮， 等. 增強T1WI全域灰度直方圖在睪丸生殖細胞瘤鑒別診斷中的價值［J］. 鄭州大學(xué)學(xué)報（醫(yī)學(xué)版）， 2020， 55（4）： 551-555.

［18］潘霞， 嚴(yán)陳晨， 劉松， 等. ADC全容積直方圖在直腸癌分期中的應(yīng)用［J］. 醫(yī)學(xué)影像學(xué)雜志， 2018， 28（1）： 124-129.

［19］白國杰， 萬業(yè)達. Bp-MRI灰度直方圖在鑒別移行帶前列腺癌與良性前列腺增生中的應(yīng)用價值［J］. 天津醫(yī)科大學(xué)學(xué)報， 2021， 27（1）： 47-51.

［20］牛俊霞， 冉云彩， 陳銳， 等. 基于磁共振表觀擴散系數(shù)及對比增強T1WI直方圖分析鑒別診斷非典型性畸胎樣/橫紋肌樣腫瘤與髓母細胞瘤［J］. 磁共振成像， 2022， 13（10）： 98-102， 120.

［21］向巍， 鐘魁松， 張振博. 一種面向城市道路的視覺圖像去霧方法［J］. 貴州大學(xué)學(xué)報（自然科學(xué)版）， 2023， 40（4）： 66-72.

［22］ 張震， 張照崎， 苗志濱， 等. 基于Harris-Hist的特征匹配及目標(biāo)定位算法［J］. 吉林大學(xué)學(xué)報（理學(xué)版）， 2021， 59（2）： 333-341.

［23］牟希農(nóng). 基于閾值鄰域均值的醫(yī)學(xué)超聲圖像去噪算法［J］. 貴州大學(xué)學(xué)報（自然科學(xué)版）， 2023， 40（1）： 75-78.

Design and Implementation of Medical MR Image Based on the

Gray Interpolation Histogram Linearization

Abstract：

In order to solve the problems of blurred edges， unreasonable brightness and poor noise suppression of medical MR image， a gray interpolation histogram linearization algorithm is proposed. This method firstly linearizes the histogram of the medical MR image with mixed noise， so as to removing part of the noise and make the brightness reach a reasonable range， and then performs polynomial interpolation operation on the obtained results to eliminate the gray difference between pixel blocks and make the image clear. The simulation results show that the brightness， peak signail-to-noise ratio（PSNR） and struetural similarity index method（SSIM） of the image obtained by this algorithm are the largest， and the PSNR is about 1.6～3.2 dB higher than other denoising methods， and the SSIM is about 5%～7% higher. This algorithm can effectively reduce the noise and keep the edge and detail information of medical MR image， and its effect is obviously better than that of adaptive histogram equalization and histogram local linearization.

Key words：

gray interpolation; histogram; medical MR image; linearization