臺灣博碩士論文加值系統

光學醫療照影技術在生醫界中扮演重要的角色，得以在不破壞組織的情況下進行組織特定物質的探測以及影像重建，像是大腦皮質的血氧分佈、淺層組織之癌細胞分佈或其他載色體物質分佈等等，並做為疾病發現之重要診斷來源。本論文提出一個可攜式之擴散光學斷層掃描(diffuse optical tomography)及二維血氧濃度百分比 (2D-oxygen concentration distribution) 處理器，另外也發展了不同的光學前向模型 (optical forward-model)做為重建演算法的驗證以及開發，並透過此模型設計出下一代新穎性的擴散光學斷層掃描重建演算法。
本論文所提出的整個系統是基於指令集架構 (instruction set architecture) 做特殊應用之處理器，每個指令皆包含一個控制訊號和觸發信號處理 (signal processing) 。另外也制訂一套外部命令 (external command) ，用於與外部設備溝通 (external communication) ，可以藉由外部命令進行最高優先權的系統修改。整個系統控制部分做各模組的調配處理 (deployment processing) ，並配合前端電路控制 (front-end control circuit) 以及資料輸出輸入 (input/output interface) ，與外部控制訊號與前端做完整的配合，並包含多項操作模式 (operation mode) 。此整合型腫瘤檢測系統 (integrated tumor detection system) 設計使用台灣積體電路製造股份有限公司90奈米 CMOS技術下線晶片 (TSMC) 90-nm CMOS technology。
另外本論文提出使用蒙地卡羅法 (Monte Carlo method) 的光子行為描述，基於較偏數學型式上做光子行為深入的模擬及探討，並藉此提供重建演算法開發上較準確且公平的驗證方式，更由此模型探討到的結果研發出一套非線性斷層重建演算法 (nonlinear iterative tomography reconstruction algorithms) ，利用光線在不同擴散距離的行為透過修改過的比爾-朗伯定律(modified Beer–Lambert law)簡化出數學前向公式 (forward equation) ，並利用數學疊代重建技術 (iterative method) 找出公式的解而求得組織的異物分佈。最後也使用蒙地卡羅法前向模型做精確的驗證與比較並證明此演算法的實現可能性。

The optical medical imaging technologies play an important role in medical profession, it can detected specific substances and rebuild image in tissue under the situation without body damaging, such as oxygen distribution in cerebral cortex, cancer location of superficial tissue or other distribution of chromophores, which can be an important diagnosis source of disease findings. In this research, we design a portable processor integrating the diffuse optical tomography (DOT) and 2D-oxygen concentration distribution. Beside, also develop different optical forward-model as a verification method of reconstruction algorithms, and design a new generation novelty reconstruction algorithms through this forward-model.
The system based on the instruction set architecture, processing control signals and trigger signals. To facilitate external communication with peripheral equipment, we make a set of the external commands to implement the system modification with highest priority. To make deployment processing, we integrate front-end control circuit and input/output interface operation mode to coordinate with the external control signal. Moreover, various operation modes are also included, as well. The chip in our integrated tumor detection system is tapped out by Taiwan Semiconductor Manufacturing Company (TSMC) 90-nm CMOS technology.
In other hand, this article proposed an optical forward-model of DOT by Monte Carlo method which is adopted to help simulate the photon propagating behavior with closer inspection. On the basis of mathematical analyses, we hope that the accurate and impartial verification as a pioneer way for new-type reconstruction algorithm. The new generation reconstruction method is nonlinear iterative tomography reconstruction algorithms. Use of light behavior under different diffusion distance to simplify the mathematical formula of Forward equation and using an iterative method to resolve the inverse problem and the result obtain the distribution of distribution of foreign matter. Finally the Monte Carlo optical forward-model also uses to verify this algorithm and proof the implemented possibility. Also, the comparison with conventional works is given and shows the strengths and weaknesses of this innovative reconstruction method.

中文摘要 i
Abstract ii
誌謝 iii
Contents ix
List of Figures xii
List of Tables xv
Chapter 1 Introduction 1
1.1 Preface 1
1.2 Diffusion Optical Tomography 2
1.3 Previous Studies 3
1.4 Motivation 5
1.5 Thesis Organization 6
Chapter 2 Integrated Processor of DOT &; NIR 7
2.1 Hardware Architecture 8
2.1.1 System Control Center 9
2.1.2 Digital Signal Processor 11
2.1.3 Front-end Sensor Controller 13
2.1.4 Input/Output Controller 14
2.2 Data Transmission Protocol 16
2.2.1 External Input Command 17
2.2.2 Input Command Decode 18
2.2.3 Internal Output Command 21
2.2.4 Output Command Encode 22
2.3 Signal Process Flow 23
2.3.1 System State 23
2.3.2 Program Register 25
2.3.3 Execution Flow 26
2.4 System Implementation 28
2.4.1 Chip Tape Out 28
2.4.2 Chip Specification 29
2.5 Test Results 31
2.5.1 Functional Verification 31
2.5.2 Power Consumption Measurement 33
2.5.3 System Integration 37
Chapter 3 Diffuse Optical Verification Model 39
3.1 Diffuse Optical Forward Model 39
3.1.1 Model Synopsis 40
3.1.2 Related Literature 40
3.1.3 Photon Behavior 43
3.2 Methods 45
3.2.1 Monte Carlo Methods 45
3.2.2 Smirnov Transform 47
3.2.3 Mathematical Expression 48
3.2.4 Propagation Rules 49
3.3 Simulation Result 51
3.3.1 Simulation Image 51
3.3.2 Analytical Solution Comparison 52
3.4 Accelerated Monte Carlo Method 54
3.4.1 Source-Detector Operation Module 54
3.4.2 Annular Receiving 56
3.4.3 Paths Database Establishment 57
3.5 Application Results 58
3.5.1 Situation 58
3.5.2 Comparison 59
3.5.3 Discussion 61
Chapter 4 Nonlinear Reconstruction Algorithm 62
4.1 Structure 62
4.1.1 Geometric Position 62
4.1.2 Detection Scale 63
4.2 Algorithm Concept 63
4.2.1 Light Propagation 63
4.2.2 Variable Definitions 65
4.3 Formula Derivation 65
4.3.1 Top Layer 66
4.3.2 Middle Layer 68
4.3.3 Bottom Layer 70
4.3.4 Reconstruction Process 71
4.3.5 Parameter Table 72
4.4 Simulation Results 74
4.4.1 Reconstructed Images 74
4.4.2 Performance 76
4.4.3 Comparison 79
Chapter 5 Conclusions 80
Reference 81

[1] Herman, G. T., Fundamentals of computerized tomography: Image reconstruction from projection, 2nd edition, Springer, 2009
[2] Hendee, et al., Magnetic Resonance Imaging Part I—Physical Principles, West J Med., 1984, 141 (4): 491–500
[3] Gary Strangman, David A Boas, Jeffrey P Sutton, Non-invasive neuroimaging using near-infrared light, Biological Psychiatry, Volume 52, Issue 7, 1 October 2002, Pages 679-693
[4] Pogue, et al., Comparison of imaging geometries for diffuse optical tomography of tissue, Opt. Express 4, 270-286 ,1999
[5] H Sato, A Study of Optical Topography for Imaging Human Brain Functions: Evaluation of Activation Signals and Its Methodological Improvement, March 2006
[6] JongKwan Choi, MinGyu Choi, Hyeon-Min Bae, An efficient data extraction method for high-temporal-and-spatial-resolution near infrared spectroscopy (NIRS) systems, Circuits and Systems (ISCAS), 2012 IEEE International Symposium on , vol., no., pp.560,563, 20-23 ,May 2012
[7] D. Joseph, T. Huppert, M. Franceschini, D. Boas, Diffuse optical tomography system to image brain activation with improved spatial resolution and validation with functional magnetic resonance imaging, Appl. Opt. 45, 8142-8151 ,2006
[8] Shih Kang, et al., Portable SoC design for CW diffusion optical tomography, Consumer Electronics (ICCE), 2012 IEEE International Conference on , vol., no., pp.341-342, 13-16 Jan. 2012
[9] Ishimaru, A., Wave Propagation and Scattering in Random Media, Academic Press,San Diego, 1978
[10] Patterson, M.S., Chance, B., Wilson, B.C., Time resolved reflectance and transmittance for the non-invasive measurement of tissue optical properties, Appl.Optics, 28, 2331, 1989
[11] Haskell, R.C. et al., Boundary conditions for the diffusion equation in radiative transfer, J. Opt. Soc. Am. A, 11, 2727, 1994
[12] A. C. Kak, M. Slaney, Principles of Computerized Tomographic Imaging, IEEE Press, New York 1988
[13] Weiwei Ma, et al., An FPGA-Based Singular Value Decomposition Processor, Electrical and Computer Engineering, 2006. CCECE '06. Canadian Conference on , vol., no., pp.1047-1050, May 2006
[14] Shih Kang, et al., A VLSI design of singular value decomposition processor for portable continuous-wave diffusion optical tomography systems, Life Science Systems and Applications Workshop (LiSSA), 2011 IEEE/NIH , vol., no., pp.100-103, 7-8 April 2011
[15] Yuan-Huang Hsu, et al., A VLSI-inspired image reconstruction algorithm for continuous-wave diffuse optical tomography systems, Life Science Systems and Applications Workshop, 2009. LiSSA 2009. IEEE/NIH pp.88-91, 9-10 April 2009
[16] Ching-Ju Cheng, et al., A SoC design for portable 2-dimension oximeter image system, Circuits and Systems (ISCAS), 2012 IEEE International Symposium, Page(s): 556 - 559, 2012
[17] S. Chandrasekhar, Radiative Transfer, Dover, New York,1960
[18] Nicholas Metropolis, S. Ulam, The Monte Carlo Method, Journal of the American Statistical Association , Vol. 44, No. 247, Sep., 1949, pp. 335-341
[19] Reddy, J. N., Introduction to the Finite Element Method, 2nd edn. McGraw-Hill, New York, 1992
[20] Li, Zhi, Joint use of early and quasi-continuous wave photons in time-gated tomographic imaging, Electrical and Computer Engineering Master's Thesis, 2011 Paper 64
[21] L. Wang, S. L. Jacques, L. Zheng, MCML—Monte Carlo modeling of light transport in multi-layered tissues, Comput. Methods Programs Biomed. 47, 131–146, 1995
[22] Tarvainen T., Computational methods for light transport in optical tomography, Doctoral Dissertation, Department of Physics, University of Kuopio, Finland, 2006
[23] J. Brock, Miriam Leeser, Heterogeneous Tasks and Conduits Framework for Rapid Application Portability and Deployment, Innovative Parallel Computing. Foundations &; Applications of GPU, Manycore, and Heterogeneous Systems, May 2012
[24] A. Liemert, A. Kienle, Analytical solution of the radiative transfer equation for the infinite-space fluence, Phys. Rev. A83, 015804, 2011
[25] A. Joshi, et al. Plane-wave fluorescence tomography with adaptive finite elements, Opt Lett, 31 ,2006, pp. 193–195
[26] E.E. Graves, et al., Singular-value analysis and optimization of experimental parameters in fluorescence molecular tomography, Opt Soc Am A, 21 ,2004, pp. 231–241
[27] Ishimaru A., Wave propagation and scattering in random media, vol. 1, Single Scattering and Transport Theory. New York: Academic Press; 1978
[28] Tarvainen T., Computational methods for light transport in optical tomography, Doctoral Dissertation, Department of Physics, University of Kuopio, Finland, 2006
[29] van Gemert, et al., Skin optics, Biomedical Engineering, IEEE Transactions on , vol.36, no.12, pp.1146-1154, Dec. 1989
[30] David A. Boas, et al., The Accuracy of Near Infrared Spectroscopy and Imaging during Focal Changes in Cerebral Hemodynamics, NeuroImage, Volume 13, Issue 1, January 2001, Pages 76-90
[31] A. Sassaroli, S. Fantini, Comment on the modified Beer-Lambert law for scattering media, Phys. Med. Biol. 49, 2004
[32] L. C. Henyey, J. L. Greenstein, Diffuse radiation in the galaxy, Astrophys. J. 93, 70–83 ,1941
[33] W. M. Irvine, Multiple scattering by large particles, Astrophys. J. 4, 1563–1575 ,1965
[34] Jacod, J., A. N. Shiryaev, Limit Theorems for Stochastic Processes ,Second Ed,. Berlin : Springer Verlag. ,2003
[35] Born, M., and E. Wolf, Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light, Sixth corrected Ed., Pergamon Press ,1986
[36] Jaeckel, P., Monte Carlo methods in finance, Stochastic Dynamics, 2001
[37] Shih-Yang Wu , Shih Kang ,Wai-chi Fang, The Forward Problem Algorithm Based on Monte Carlo Method for Diffusion Optical Tomography, World Congress in Computer Science, Computer Engineering, and Applied Computing, 2011. Worldcomp'11, July 18-21, 2011
[38] T. J. Farrell, M. S. Patterson, B. C. Wilson, A diffusion theory model of spatially resolved, steady-state diffuse reflectance for the noninvasive determination of tissue optical properties in vivo, Med. Phys. 19, 879–888 ,1992
[39] A. Liemert, A. Kienle, Analytical solution of the radiative transfer equation for the infinite-space fluence, Phys. Rev. A83, 015804 ,2011
[40] Qianqian Fang, et al., Nonlinear image reconstruction algorithm for diffuse optical tomography using iterative block solver and automatic mesh generation from tomosynthesis images. Proc. SPIE 6081, Multimodal Biomedical Imaging, 60810O ,February 13, 2006
[41] Boas D A, et al., The accuracy of near infrared spectroscopy and imaging during focal changes in cerebral hemodynamics, Neuroimage 13 76–90, 2001