超音波成像利用探頭發射聲波進入人體內，而與人體內的各種介質發生交互作用，而利用反射回的信號重建人體內影像，陣列超音波系統由於可以進一步調整探頭單元間的時間差而有更好的成像品質與應用而被廣泛地應用在醫學診斷上。然而在目前影像系統中，人體內組織不均勻的特性卻會對成像過程造成影響，舉凡組織聲速、折射、散射、衰減係數…等因素均會隨著不同的組織而有不同的物理效應，進而影響到成像的品質與診斷正確率。本論文將對於偵測平均聲速與亮度衰減補償兩方面提出自動化校正方法，期待藉由提升影像品質而增進診斷正確率。在聲速偵測方面，使用平均聲速成像在超音波成像系統中可降低因組織不均勻而造成影像品質下降的影響，影像品質包括了空間與對比解析度以及病理位置的判定…等，在目前成像系統中，往往採用預設好的聲速(例如:1540m/s)進行成像，而在此預設聲速與真正平均聲速之間的誤差會導致影像品質降低進而影響到診斷效果，因此本論文以平行化架構同時分析不同聲速下的影像品質來進行自動化聲速預估與校正，此分析方法為利用不同位置的子孔徑成像並計算誤差而推估正確聲速。而在亮度衰減補償方面，超音波成像過程之聲波傳遞的過程中，舉凡組織種類、入射角度、有效頻率與孔徑的下降…等均會導致在超音波影像之亮度與動態範圍隨著深度增加而有非規則降低的現象，診斷上可能會造成病理誤判，而在成像系統中，往往以固定預設的線性補償來降低此影響，加上依靠操作人員的診斷與經驗經由調整不同深度的亮度增益與動態範圍才可得到較為真實較為可信的無衰減影像，此過程十分費時與不精確，故此論文亦提出了自動化亮度增益補償方法以改善此人工操作之問題，此補償方法包含了先利用統計分區標定出主要組織以判斷其衰減曲線，接著對其安排了二維查找表與線性組合計算各取樣點的亮度補償，亮度補償後為了處理可能造成信雜比下降的傷害，此方法於最後加上了隨深度變化的選擇性對比加強器以適當調整影像對比，此方法可有效增進影像均勻度使影像品質提升。最後，本論文將此兩種自動化演算法於平行化軟體平台上(CUDA)加以實現，並達到在診斷前以短時間內進行最佳化影像之目標。

In ultrasound array imaging, using correct sound speed is one of the most important tasks to ensure beam-forming accuracy. However, due to variations in body type and tissue inhomogeneities, the actual sound speed ranging from 1400m/s to 1600m/s is hard to be accurately applied during beamforming. In most current systems, the sound velocity is pre-determined as a constant that would influence whether or not the image representation of the anatomy is reliable for diagnosis. In this research, a global mean sound speed estimation method is presented. It utilizes a sub-aperture imaging technique to form two images with different sound speed values, and choose the optimal one with the minimum mean error between these images as the estimated sound speed value. This optimal sound speed effectively provides higher contrast-to-noise ratio on experimental phantom images as well as clinical images. The other part of the research is regarding tissue attenuation that causes brightness changes in the image and thus an adaptive depth gain compensation method is needed. In current systems, incorrect compensation may be presented and the operators needs to adjust the gain and dynamic range manually. Therefore, in this research, an adaptive depth gain compensation algorithm is proposed with the purpose of automatically providing better image uniformity. This compensation algorithm first calculates attenuation curves from main tissue regions of images, and based on these curves, a 2d look-up table for linear composition could be generated to adapt the brightness pixel-by-pixel. Furthermore, in order to decrease influences of signal-to-noise ratio degradation by compensation, an adaptive contrast enhancement is added to last stage of this method. Finally, the computations of these two algorithms can be finished within 1 second on GPU platform.

口試委員會審定書 #
誌謝 i
摘要 ii
ABSTRACT iii
目錄 iv
圖目錄 vi
表目錄 ix
Chapter 1 緒論 1
1.1 超音波陣列成像系統 1
1.2 聲速與成像 2
1.3 亮度衰減與深度增益補償 4
1.4 研究動機與目標 6
1.5 論文架構 7
Chapter 2 自動化平均聲速偵測 9
2.1 組織聲速誤差與影像品質 9
2.1.1 空間解析度 9
2.1.2 對比解析度 11
2.2 聲速偵測原理及方法 12
2.2.1 先前技術偵測方法與原理 14
2.2.2 雙子孔徑成像法偵測聲速與比較 16
2.3 多聲速掃描系統架構 20
2.4 聲速偵測系統架構 22
2.4.1 雙子孔徑與多聲速成像 23
2.4.2 影像誤差趨勢與最佳聲速偵測 25
Chapter 3 自動化深度增益補償演算法 28
3.1 深度增益補償方法及原理 28
3.1.1 衰減曲線判定 29
3.1.2 補償器設計 30
3.1.3 對比度放大器設計 32
3.2 自動化深度補償系統架構 33
3.2.1 影像分群衰減曲線判定 33
3.2.2 影像深度增益補償器 34
3.2.3 選擇性偵測對比加強器 35
Chapter 4 結果與討論 37
4.1 偵測聲速效率與誤差 37
4.1.1 標準仿體影像 37
4.1.2 加入不同相位偏移影像 40
4.1.3 人體內影像 43
4.2 聲速偵測其他方法與比較 48
4.3 自動化深度增益補償 51
4.3.1 標準仿體影像 51
4.3.2 人體內影像 53
4.4 亮度增益補償其他方法與比較 55
4.4.1 組織分區判斷衰減之相關方法 55
4.4.2 自動化補償方法 57
4.5 自動化偵測聲速於CUDA平行化軟體系統平台之呈現 58
4.6 自動化深度增益補償於CUDA軟體系統之呈現 58
Chapter 5 結論與未來工作 60
5.1 自動化偵測聲速 60
5.2 自動化深度增益補償 61
REFERENCE 63

[1]Y. F. Li and P. C. Li, “Software beamforming: comparison between a phased array and synthetic transmit aperture,” Ultrasonic Imaging, vol. 33, pp. 108-119, 2011.
[2]M. E. Anderson, M.S. McKeag, and G. E. Trahey, ”The impact of sound speed errors on medical ultrasound imaging”, J. Acoust. Soc. Am., vol. 107, pp. 3540-3548, 2000.
[3]M. H. Cho, L.H. Kang, J.S. Kim, and S.Y. Lee, "An efficient sound speed estimation method to enhance image resolution in ultrasound imaging," Ultrasonics, vol. 49, pp. 774-778 , 2009.
[4]J. F. Krucker, J. B. Fowlkes, and P. L. Carson, “Sound speed estimation using automatic ultrasound image registration”, IEEE Transaction on Sonic and Ultrasonics, vol. 51 , no. 9 , pp. 1095-1106, 2004.
[5]C. H. Yoon, Y. H. Lee, J. H. Chang, T. K. Song, and Y. M. Yoo, "In vitro estimation of mean sound speed based on minimum average phase variance in medical ultrasound imaging," Ultrasonics, vol. 51, pp. 795-802, 2011.
[6]X. L. Qu, T. Azuma, J. T. Liang, and Y. Nakajima, “Average sound speed estimation using speckle analysis of medical ultrasound data, ” Int. J. CARS, vol. 7, no. 6, pp. 891-899, 2012.
[7]M. W. Tang, and D.C. Liu, “Rationalized gain compensation for ultrasound imaging,” APCMBE, IFMBE Proceedings, vol. 19, pp. 282-285, 2008.
[8]M. L. Li, and P. C. Li, “A new adaptive imaging technique using generalized coherence factor,” 2002 IEEE Ultrasonics Symposium Proceedings, pp.1627-1630, 2002.
[9]P. C. Li and M. L. Li, “Adaptive imaging using the generalized coherence factor,” IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control, vol. 50, pp.128-141, 2003.
[10]R. F. Wagner, S. W. Smith, J. M. Sandrik, and H. Lopez, "Statistics of speckle in ultrasound B-scans," IEEE Transaction on Sonic and Ultrasonics, vol. 30, no. 3 , pp. 156-163 , 1983.
[11]D. Lee, Y. S. Kim and J. B. Ra, “Automatic time gain compensation and dynamic range control in ultrasound imaging systems,” SPIE Medical Imaging, vol. 6147, pp. 614708-1 – 614708-9, 2006.
[12]S. W. Huang and P. C. Li, “Ultrasound computed tomography reconstruction of the attenuation coefficient using a linear array,” IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control, vol. 52, no. 11, pp. 2011-2022, 2005.
[13]S. W. Huang and P. C. Li, “Computed tomography sound velocity reconstruction using incomplete data,” IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control, vol. 51, no. 3, pp. 329-342, 2004.
[14]S. W. Huang and P. C. Li, “Experimental investigation of computed tomography sound velocity reconstruction using incomplete data,” IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control, vol. 51, no. 9, pp. 1072-1081, 2004.
[15]M. Krueger, V. Burow, K. M. Hiltawsky, and H. Ermert, “Limited angle ultrasonic transmission tomography of the compressed female breast,” IEEE Ultrasonics Symposium Proceedings, pp. 1345-1348, 1998.
[16]H. Kim, J. A. Zagzebski, and T.Varghese, “Estimation of ultrasound attenuation from broadband echo-signals using bandpass filtering,” IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control, vol. 55, no. 5, pp. 1153-1159, 2008.
[17]K. Horsch, M. L. Giger, L. A.Venta, and C. J. Vyborny, “Computerized diagnosis of breast lesions on ultrasound,” Medical Physics, vol. 29, no. 2, pp. 157-164, 2002.
[18]G. Kossoff, E. K. Fry, and J. Jellins, ”Average velocity of ultrasound in the human female breast,” J. Acoust. Soc. Am., vol. 93, no. 3, pp. 1609-1612, 1993.
[19]L. Landini, R. Sarnelli, and F. Squartini, “Frequency-dependent attenuation in breast tissue characterization,” Ultrasound in Medicine and Biology Society, pp. 395-396, 1989.
[20]K. Hensel, G. Li, and G. Schmitz, “Evaluation of the local speed of sound estimation for the correction of ultrasound compound imaging by speckle analysis,” Medical Physics and Engineering World Congress, 2009.
[21]Q. Zhu and B. D. Steinberg, “Large-transducer measurement of wavefront distortion in the female breast,” Ultrasonics Imaging, vol. 14, no. 3, pp. 276-299, 1992.
[22]J. S.Schreiman, J. J. Gisvold, J. F. Greenleaf, and R. C. Bahn, “Ultrasound transmission computed tomography of the breast,” Radiology, vol. 150, no. 2, pp. 523-530, 1984.
[23]R. F. Chang, W. J. Wu, and D. R. Chen, “Automatic ultrasound segmentation and morphology based diagnosis of solid breast tumors,” Breast Cancer Research and Treatment, vol. 89, pp. 179-185, 2005.
[24]J. A. Noble, and D. Boukerroi, “Ultrasound image segmentation: a survey,” IEEE Transactions on Medical Imaging, vol. 25, no. 8, pp. 987-1010, 2006.
[25]M. C. Axelsen, K. F. Poeboe, and J. A. Jensen, “Evaluation of automatic time gain compensated in-vivo ultrasound sequences”, 2010 IEEE Ultrasonics Symposium Proceedings, pp. 1640-1643, 2010.