康秀娟，施 展

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

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

康秀娟，施 展

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

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

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

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

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

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

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

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

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

圖2 圖形用戶界面總體設計

2 圖像預處理

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

2.1 直方圖增強的改進算法

由于直方圖均衡算法經(jīng)常會顯著改變圖像特征[1-4]，本文采用直方圖增強的改進算法,將直方圖規(guī)定化與局部對比度增強相結合，進行圖像增強處理。直方圖規(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為處理窗口內所含像素的平均灰度值。k(f(x,y)-Mx,y)為局部對比度增強后的灰度值，相當于一個高通濾波器，起到強化細節(jié)的作用。k值根據(jù)處理窗口的灰度方差來選取，可表達為

(2)

k=k′(t2-1)

(3)

圖3 k值變化示意圖

2.2 自適應中值濾波的改進算法

圖4 改進的自適應中值濾波算法流程圖

3 圖像幾何校正

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

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

(4)

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

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

圖5 雙線性內插法鄰近點選擇

雙線性內插法利用待求像素周圍的4個鄰近像素在x、y方向上分別做加權線性插值。如圖5所示，設(x′,y′)4個鄰近像素點為A、B、C、D，其坐標分別為(i,j)，(i+1,j)，(i,j+1)，(i+1,j+1)，灰度值分別為g(A)、g(B)、g(C)、g(D)，設α=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)標定與異物定位測量

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

k=P/L

(8)

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

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

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

(9)

圖6 圖像坐標系

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

Z=Z′/sinθ

(10)

可將異物點的三維坐標信息求出。

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

(11)

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

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

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

表1 10 mm標準尺寸測量結果 /mm

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

6 結束語

本文介紹了人體異物圖像定位測量系統(tǒng)的設計，該設計系統(tǒng)具有圖像增強、濾波、系統(tǒng)標定、異物坐標及采點測距等功能。系統(tǒng)誤差分析結果表明本文設計軟件標定誤差較小，尺寸測量比較精確，系統(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] 龍清.直方圖圖像增強技術[J].電腦知識與技術,2011,7(4):883-886.

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

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

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

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

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

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

[9] 郭羽,楊紅,楊照金,等.CCD攝像系統(tǒng)鏡頭的畸變測量[J].應用光學,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] 符祥,郭寶龍.圖像插值技術綜述[J].計算機工程與設計,2009,30(1):141-144.

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

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

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

[15] 劉浩,韓晶.Matlab R2012a完全自學一 本通[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

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

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

TP391.41

A

1007-7820(2018)01-071-04