康秀娟，施 展

(上海理工大學(xué) 光電信息與計算機工程學(xué)院，上海 200093)

基于圖形處理的人體內(nèi)異物定位測量系統(tǒng)設(shè)計

康秀娟，施 展

(上海理工大學(xué) 光電信息與計算機工程學(xué)院，上海 200093)

臨床醫(yī)學(xué)實踐中，存留的人體異物會給病人身心帶來嚴(yán)重的傷害。因此，快速、準(zhǔn)確地對人體內(nèi)異物進行三維定位，具有重要的研究意義和廣泛的應(yīng)用前景。文中對異物圖像預(yù)處理、幾何校正以及人體異物定位測量等關(guān)鍵技術(shù)進行了深入分析和研究，并設(shè)計開發(fā)了基于Matlab GUI的人體異物圖像測量分析軟件。該軟件系統(tǒng)能實現(xiàn)對醫(yī)學(xué)圖像的預(yù)處理以及對人體內(nèi)異物平面位置和深度信息的定量描述，實驗結(jié)果表明，系統(tǒng)測量誤差不超過0.21 mm。

圖像處理；幾何校正；系統(tǒng)標(biāo)定；異物定位；采點測距

臨床醫(yī)學(xué)中，存留的人體異物對人身心造有極大的危害。雖然目前可以用X線檢查等方法判斷異物位置，但由于圖像不能對異物進行三維定位，經(jīng)常存在定位不夠準(zhǔn)確等問題。因此，快速、準(zhǔn)確判斷傷病員體內(nèi)異物的位置對提高手術(shù)準(zhǔn)確率和安全性具有重要意義。

1 人體異物定位測量系統(tǒng)設(shè)計

為確定異物在人體內(nèi)的平面位置和深度信息，需要至少兩幅不同角度的異物圖像進行測量定位，本文主要通過對兩幅人體異物圖像的處理與分析來確定異物的三維坐標(biāo)、大小尺寸等，定量描述人體內(nèi)異物的平面位置和深度信息。

按照系統(tǒng)的功能要求，本文基于圖形處理的人體內(nèi)異物定位測量系統(tǒng)結(jié)構(gòu)主要分為4部分：圖像預(yù)處理、圖像幾何校正、系統(tǒng)標(biāo)定和人體異物定位測量。系統(tǒng)總體設(shè)計方案如圖1所示。

圖1 人體異物定位測量系統(tǒng)設(shè)計方案

圖形用戶界面總體設(shè)計如圖2所示。界面布局主要分為3部分：圖像顯示區(qū)域、操作區(qū)域及分析結(jié)果輸出區(qū)域。其中顯示區(qū)域用于顯示異物圖像的正側(cè)位投影及處理后圖像，兩側(cè)操作區(qū)域通過功能鍵實現(xiàn)對圖像的各種處理與分析，界面下方為分析結(jié)果輸出欄，用于顯示系統(tǒng)標(biāo)定系數(shù)、異物坐標(biāo)及采點測距的結(jié)果輸出。

圖2 圖形用戶界面總體設(shè)計

2 圖像預(yù)處理

本文采用直方圖增強的改進算法和自適應(yīng)中值濾波的改進算法對圖像進行預(yù)處理，提高圖像清晰度和信噪比。

2.1 直方圖增強的改進算法

由于直方圖均衡算法經(jīng)常會顯著改變圖像特征[1-4]，本文采用直方圖增強的改進算法,將直方圖規(guī)定化與局部對比度增強相結(jié)合，進行圖像增強處理。直方圖規(guī)定化可以修正圖像灰度分布，局部對比度增強算法通過處理窗口對圖像細節(jié)處做局部對比度增強，從而減小對圖像細節(jié)的損壞。得到的圖像在(x,y)處的灰度值可表示為

g′(x,y)=G(x,y)+k(f(x,y)-Mx,y)

(1)

式(1)中，G(x,y)為規(guī)定化處理后的灰度值，f(x,y)為圖像未處理前(x,y)處的灰度值，Mx,y為處理窗口內(nèi)所含像素的平均灰度值。k(f(x,y)-Mx,y)為局部對比度增強后的灰度值，相當(dāng)于一個高通濾波器，起到強化細節(jié)的作用。k值根據(jù)處理窗口的灰度方差來選取，可表達為

(2)

k=k′(t2-1)

(3)

圖3 k值變化示意圖

2.2 自適應(yīng)中值濾波的改進算法

圖4 改進的自適應(yīng)中值濾波算法流程圖

3 圖像幾何校正

由于成像系統(tǒng)具有非線性以及拍攝角度的問題，會使生成的圖像存在一定的幾何畸變[7]。本文通過空間坐標(biāo)變換和對非整數(shù)處坐標(biāo)進行雙線性插值兩個步驟校正畸變。

(1)假設(shè)圖像為雙線性畸變，校正公式為

(4)

式中共有8個參數(shù)，可通過圖像中的4個點對將系數(shù)解出；

(2)空間變化后的像素坐標(biāo)會產(chǎn)生非整數(shù)，因此，坐標(biāo)變換后還需確定各像素的灰度值。像素灰度內(nèi)插的方法主要有最近鄰點法、雙線性內(nèi)插法和三次內(nèi)插法3種[8-10]。本文采用雙線性內(nèi)插法對非整數(shù)處的坐標(biāo)進行插值。雙線性內(nèi)插法又稱一階插值法[11]，其需要經(jīng)3 次插值才能獲得最終結(jié)果，是對最近鄰插值法的一種改進[12]。

圖5 雙線性內(nèi)插法鄰近點選擇

雙線性內(nèi)插法利用待求像素周圍的4個鄰近像素在x、y方向上分別做加權(quán)線性插值。如圖5所示，設(shè)(x′,y′)4個鄰近像素點為A、B、C、D，其坐標(biāo)分別為(i,j)，(i+1,j)，(i,j+1)，(i+1,j+1)，灰度值分別為g(A)、g(B)、g(C)、g(D)，設(shè)α=x′-i，β=y′-j。首先計算出E、F兩處的灰度值g(E)和g(F)

g(E)=β[g(C)-g(A)]+g(A)

(5)

g(F)=β[g(D)-g(B)]+g(B)

(6)

則(x′,y′)點的灰度值g(x′,y′)為

g(x′,y′)=α[g(F)-g(E)]+g(E)

(7)

g(x′,y′)即為校正后圖像中(x′,y′)處的灰度值。

4 系統(tǒng)標(biāo)定與異物定位測量

系統(tǒng)標(biāo)定[13]可以實現(xiàn)圖像像素點數(shù)與實際尺寸的轉(zhuǎn)換。這里將像素與實際尺寸之間的比例關(guān)系定義為標(biāo)定系數(shù)k

k=P/L

(8)

式中，L表示被測物體的實際尺寸，單位mm，P表示被測物的像素個數(shù)，k表示系統(tǒng)標(biāo)定系數(shù)，單位pixel/mm。

本文通過選取一個標(biāo)準(zhǔn)件的方法，將待測量標(biāo)準(zhǔn)件的精確尺寸與經(jīng)系統(tǒng)處理后的數(shù)字圖像的像素個數(shù)作比對，根據(jù)公式(8)多次實驗取平均值，從而計算出系統(tǒng)的標(biāo)定系數(shù)[14-13]。

獲取標(biāo)定系數(shù)后，還需建立坐標(biāo)系，實現(xiàn)像素坐標(biāo)與物理坐標(biāo)之間的換算。本文建立圖像坐標(biāo)系如圖6所示，將原點規(guī)定位于圖像的左下角，設(shè)O1在u，v坐標(biāo)系中的像素坐標(biāo)為(u0,v0)，系統(tǒng)標(biāo)定系數(shù)為k，則O1點的物理坐標(biāo)(x,y)可表示為

(9)

圖6 圖像坐標(biāo)系

假設(shè)兩幅人體異物圖像的相位夾角為θ(一般為90°)，異物點側(cè)位投影的橫坐標(biāo)為Z′，根據(jù)幾何變換關(guān)系

Z=Z′/sinθ

(10)

可將異物點的三維坐標(biāo)信息求出。

根據(jù)兩點間距離公式

(11)

通過采點測距的方式可實現(xiàn)圖像中任意兩點間距離的計算[15]。

5 系統(tǒng)誤差分析

本文選取分辨率為96 dpi，尺寸大小為981×698像素的標(biāo)準(zhǔn)尺寸圖像，在圖像進行預(yù)處理、幾何校正及系統(tǒng)標(biāo)定后，利用采點測距功能對10 mm標(biāo)準(zhǔn)尺寸進行重復(fù)測量，10次測量結(jié)果和誤差分析如表1所示。

表1 10 mm標(biāo)準(zhǔn)尺寸測量結(jié)果 /mm

根據(jù)誤差理論，極限誤差取標(biāo)準(zhǔn)差的3倍時，極限誤差δ=0.202 2 mm ，滿足人體異物定位測量時的實際要求。

6 結(jié)束語

本文介紹了人體異物圖像定位測量系統(tǒng)的設(shè)計，該設(shè)計系統(tǒng)具有圖像增強、濾波、系統(tǒng)標(biāo)定、異物坐標(biāo)及采點測距等功能。系統(tǒng)誤差分析結(jié)果表明本文設(shè)計軟件標(biāo)定誤差較小，尺寸測量比較精確，系統(tǒng)整體測量誤差不超過0.21 mm，滿足實際要求。

[1] Wang L.Enhancement of medical ultrasonic image based on gray-level histogram equalization[J].Journal of Sichuan University,2002,34(1):105-108.

[2] 龍清.直方圖圖像增強技術(shù)[J].電腦知識與技術(shù),2011,7(4):883-886.

[3] 趙秋宇,王曉紅,張德喜.可增強圖像細節(jié)的改進的直方圖均衡化算法[J].信陽師范學(xué)院報,2008,21(4):601-602.

[4] 萬輝，楊有.一種改進的局部對比度圖像增強算法[J].重慶文理學(xué)院報:自然科學(xué)版,2007,26(3):29-32.

[5] 茍中魁,張少軍,李忠富,等.一種基于極值的自適應(yīng)中值濾波算法[J].紅外與激光,2005,34(1):98-101.

[6] 郭海霞，解凱.一種改進的自適應(yīng)中值濾波算法[J].中國圖像圖形學(xué)報,2007,12(7):1185-1188.

[7] 尹青山.數(shù)字圖像幾何校正方法[J].丹東師專學(xué)報,2003,25(3):148-149.

[8] 崔洪州,孔淵,周起勃,等.基于畸變率的圖像幾何校正[J].應(yīng)用光學(xué),2006,27(3):183-185.

[9] 郭羽,楊紅,楊照金,等.CCD攝像系統(tǒng)鏡頭的畸變測量[J].應(yīng)用光學(xué),2008,29(2):279-282.

[10] Lin J M,Xing M L,Sha D G.Distortion measurement of CCD imaging system with short focal length and large-field objective[J].Optics and Lasers in Engineering,2005(43):1137-1144.

[11] 符祥,郭寶龍.圖像插值技術(shù)綜述[J].計算機工程與設(shè)計,2009,30(1):141-144.

[12] 唐勛俊.基于圖像處理的擠出胎面尺寸控制應(yīng)用[J].電子科技,2015,28(1):140-144.

[13] 王永皎,郭力爭.基于圖像處理技術(shù)的尺寸測量中邊緣定位算法[J].信陽師范學(xué)院報,2014,27(1):135-138.

[14] 吳德剛.基于圖像處理的零件二維幾何尺寸測量算法研究[D].鄭州:鄭州大學(xué)，2011.

[15] 劉浩,韓晶.Matlab R2012a完全自學(xué)一 本通[M].北京:電子工業(yè)出版社,2014.

A System of Foreign Object Localization Based on Graphics Processing

KANG Xiujuan,SHI Zhan

(School of Optical-Electrical and Computer Engineering,University of Shanghai for Science and Technology,Shanghai 200093,China)

In the clinical medical practice, small foreign object in the human body has brings to the patient a great impact of the physical and mental. Therefore quickly and accurately to three- dimensional location the foreign object in human body, has significant research value and application prospect. This paper focused on image preprocessing, geometric correction and human body localization carried on the thorough analysis and research，and designed the foreign object measurement and analysis software based on GUI MATLAB. The software system can realize the pretreatment of medical image and the quantitative description of the plane position and depth information of the foreign object. The experimental results show that the system measurement error less than 0.21mm.

image processing; geometric correction; system calibration; foreign object localization; ranging

2017- 02- 23

上海市科學(xué)技術(shù)委員會科研基金(15110502300)

康秀娟(1993-)，女，碩士研究生。研究方向：數(shù)字圖像處理。施展(1963-)，女，副教授。研究方向：測試?yán)碚撆c測試技術(shù)。

TP391.41

A

1007-7820(2018)01-071-04