【作 者】莊偉龍，魏錫安，賈春榕，洪衍界，石錦添

福建省級機關醫(yī)院，福建，福州，350003

1 故障現(xiàn)象

我院VCT 在做頭顱掃描后,進行拆薄(32*0.625或者64*0.625) 時,在前兩層會出現(xiàn)如圖1所示的環(huán)狀偽影。用水模直接進行薄層掃描（32*0.625或者64*0.625）時，也能看到前兩層出現(xiàn)環(huán)狀偽影，如圖2所示。

圖1 正常頭顱掃描時的環(huán)狀偽影

圖2 水模測試時的環(huán)狀偽影

2 分析與檢修

出現(xiàn)環(huán)狀偽影的可能原因很多，我們只能根據(jù)經(jīng)驗及設備手冊進行排查分析。查看error log，機器并沒有出現(xiàn)打火現(xiàn)象及其它報錯提示。檢查Mylar及球管的射線出口以及探測器表面，也未發(fā)現(xiàn)有污垢，可見并不是因為球管氣泡或者射線被臟污影響而導致的偽影?？赡苁切饰募霈F(xiàn)錯誤，進行Detail calibration，發(fā)現(xiàn)偽影現(xiàn)象沒有改變。重新安裝系統(tǒng)軟件現(xiàn)象不變，由此可排除軟件故障的可能性。

這時我們推測，應該是在接受X-ray后處理出現(xiàn)的問題，即DAS或Recon部分的問題。DAS主要的功能是接收到穿透人體后到達Detector Module的X-ray，由FRDM將其轉化為電信號，通過DIFB再傳送到主機的Recon部分。

我們先開始做一個Diagnostic Data Collection，然后做Scan Analysis，分析DAS搜集到的數(shù)據(jù)是否存在干擾，發(fā)現(xiàn)在Detector Row 1B,2B,3B上搜集到的信號存在干擾（信號波形上有明顯毛刺），如圖3所示，位于1B的291channel，2B 的291channel，3B的291、243channel位置。

經(jīng)過數(shù)據(jù)分析，信號干擾對應的DIFB是7、9兩塊，F(xiàn)RDM（Detector module）是17、19兩塊（DIFB板共有10塊，F(xiàn)RDM板共有57塊）。將有可能存在問題的兩塊DIFB和FRDM板分別與其它位置的同型號板進行交換，交換7#、9#兩塊DIFB板后，發(fā)現(xiàn)故障沒有改變；而交換了FRDM17#、19#兩塊板后，重新做Diagnostic Data Collection 和Scan Analysis，發(fā)現(xiàn)信號波形的干擾還在，但是對應的channel位置發(fā)生了變化，變成1B的243channel和2B 的243channel，3B的243、291channel。

由此可以判定，F(xiàn)RDM 17#、19#兩塊板在處理數(shù)據(jù)時出現(xiàn)問題，導致信號失真，出現(xiàn)了偽影。在更換了兩塊FRDM板后，重新進行Detail calibration，故障現(xiàn)象消失。血壓信號卻相反，其P值大于1，即3～5 Hz的能量反而比1～3 Hz高。高血壓信號的5 Hz以上的高頻能量比正常信號高很多，這是公認的事實。而從表1中可以大膽推論，高頻部分的能量增加不是突然發(fā)生的，而是5 Hz后半段能量不斷增加“抬高”，溢出5 Hz的分界線而造成的。因此，研究5Hz以內(nèi)的能量分布，也能判斷是否發(fā)生了高血壓病變。

圖6 第4層細節(jié)信號的功率譜密度Fig.6:the 4th detail signal’s power spectral density

3 結論

通過脈象頻譜研究，發(fā)現(xiàn)正常人的動脈信號的大部分能量集中在10 Hz以內(nèi)，而病理性脈搏信號則有相當一部分能量集中在高頻段10 Hz以上[1]。然而在5 Hz以內(nèi)，病理信號的能量分布與健康信號相比也存在差異。通過小波細節(jié)分解把5 Hz以內(nèi)的頻段進行分段，得出高血壓信號在5 Hz后半段3～5 Hz頻段內(nèi)能量顯著提高，與1～3 Hz段的能量比值增加到大于1。因此本方法可以通過對5 Hz內(nèi)的能量分布來分析是否有病變發(fā)生。

此外，由于高血壓可以分成很多種類，每個種類的主導病變的頻率肯定有所不同，不能同一而論，所以針對某類高血壓的能量分布情況還需要進一步研究。

[1] 羅志昌, 張松, 楊益民. 脈搏波的工程分析與臨床應用[M]. 北京:科學出版社, 2006.

[2] Chun T Lee, Ling Y Wei. Spectrum Analysis of Human Pluse[J].IEEE Trans.on BME, 1983, 30(6): 348-351.

[3] 廖云朋, 劉曼. 二次樣條小波的P波檢測方法[J]. 北京生物醫(yī)學工程, 2010, 29(3): 266-269.

[4] 許瑞慶. 脈搏檢測系統(tǒng)設計及其信號處理算法研究[D]. 南京: 南京信息工程大學, 2008

[5] 姜潤翔, 史建偉, 龔沈光. 基于尺度—小波能量譜的船舶水壓場信號檢測[J]. 電子學報, 2009, 37(1): 122-125.

[6] 程軍圣, 于德介, 鄧乾旺, 等.時間—小波能量譜在滾動軸承故障診斷中的應用[J]. 振動與沖擊，2004, 23(2): 34-36.

[7] 葛哲學, 陳仲生.Matlab時頻分析技術及其應用[M]. 北京: 人民郵電出版社, 2006.

[8] 柴新禹, 吳朝霞. 時頻分析方法及其在醫(yī)學信號處理中的應用[J].生物醫(yī)學工程學雜志, 2001, 18(1): 138-144.