隨著數位醫學影像的發展，藉由電腦科技將醫學影像做特殊的分析與適當的處理後，可將特定組織或病灶保留或凸顯出來，使醫生在術前能做出正確的診斷與規劃手術的路徑，術中也可做為手術定位的圖像界面。然而光從單一種類的影像上無法獲得足夠的資訊，所以將多種類的醫學影像融合在一起是有其重要性。
本研究是針對頭部兩種不同特性的醫學影像: 電腦斷層影像(Computed Tomography, CT)與核磁共振影像（Magnetic resonance images, MRI）做影像處理。首先，將電腦斷層利用三角網格法(Marching Cubes)重建而核磁共振影像則是利用實體描繪法(Volume Rendering)重建，接著利用反覆疊代法(Iterative Closest Point)進行影像對位(Registration)接著利用共同資訊法(Mutual Information)進行兩種立體影像之融合(Fusion)，利用模擬軟體將最後融合後之立體圖進行軸向(Axial)切片、矢向(Sagittal)切片與冠向(Coronal)切片之距離的量測，提供醫生在術前能更精確地知道病灶的空間位置，也可輔助醫生在手術前之規劃，提升手術之成功率。

With the development of digital medical images, after special analysis and appropriate processing of medical images based on computer technology, the specific tissues or foci can be reserved or highlighted. As a result, the doctors can make correct diagnosis and plan surgical approach before operation, and an image interface for surgical positioning is provided in operation. However, a single type of image cannot provide adequate information, thus it is important to fuse diversified medical images.
This study processed two kinds of medical images with different characteristics of head, computed tomography and magnetic resonance images. First, the computed tomography was reconstructed by using Marching Cubes, and the magnetic resonance image was reconstructed by using magnetic resonance image. Afterwards, the Iterative Closest Point was used for image registration, and then the Mutual Information was used to fuse the two kinds of 3D images. The distances of axial slice, sagittal slice and coronal slice of the fused 3D image were measured by using simulation software, so that the doctors could know the spatial location of foci accurately before operation. The proposed method can assist the doctors in pre-operational planning, thus increasing the success rate of operation.

[1] Available: http://www.doh.gov.tw/cht/index.aspx
[2]Steve W., “The physics of three dimensional radiation therapy: conformal radiotherapy, Radiosurgery and Treatment Planning (Medical Science) ,” Taylor & Francis, (1993).
[3]Qiu W., Titley-Peloquin D, Soleimani M., “Blockwise conjugate gradient methods for image reconstruction in volumetric CT,” Computer methods and programs in biomedicine , Vol. 108, Issue 2, pp.669-678,(2012).
[4]Matteo, N., Domenico, S., Luigi, M. C., Joaquim, E., Isam, A.,Guadalupe, S., Joan, B. G., Enrique, F., and Alberto P. G., “The use of a three-dimensional novel computer-based model for analysis of the endonasal endoscopic approach to the midline skull base,’’ World Neurosurgery, Vol. 75, No. 1, pp.106-113 ,(2011).
[5]Annette M., Pham, A., Rafi M., Metzger, J., and Strong, E.B. “Computer modeling and intraoperative navigation in maxillofacial surgery,” Otolaryngology - Head and Neck Surgery, Vol. 137, No. 4, pp. 624-631, (2007).
[6]Georg E., Hitomi S., Gavin K., and Joachim M., “The accuracy of image guided surgery based on cone beam computer tomography image data,”Oral Surg Oral Med Oral Pathol Oral Radiol Endod, pp.41-48,(2009).

[7]Haralambidis A., Ari-Demirkaya A., Acar A., Kucukkele N., Ate M., and Ozkaya S., “Morphologic changes of the nasal cavity induced by rapid maxillary expansion: A study on 3-dimensional computed tomography models,” American Journal of Orthodontics and Dentofacial Orthopedics, Vol. 136, No. 6, pp.815- 821, (2009).
[8]Biedert R., Sigg A., Gal I., and Gerber H., “3D representation of the surface topography of normal and dysplastic trochlea using MRI,’’ The Knee, Vol. 18, No 5, pp. 340-346 , (2011).
[9]Bae K., Shim H., Tao C., Chang S., Wang J., Boudreau R., and Kwoh C., “Intra-and inter-observer reproducibility of volume measurement of knee cartilage segmented from the OAI MR image set using a novel semi-automated segmentation method,” Osteoarthritis and Cartilage, (2009).
[10] Wang R., Yang D., and Li S., “Three-dimensional virtual model and animation of penile lengthening surgery,” Journal of Plastic, Reconstructive & Aesthetic Surgery, Available online 2 May, (2012).
[11]Clatz O., Delingette H., Talos I. F., Golby A. J., Kikinis R.,. Jolesz F. A, et al, “Robust non-Rigid registration to capture brain shift form intra-operative MRI”, in proceedings of MICCAI’05, pp. 295-302, (2007).
[12]Brem M.H., Lang P.K., Neumann G., Schlechtweg P.M., Schneider E., Jackson R., Yu J., Eaton C.B., Hennig F.F., Yoshioka H., Pappas G., and Duryea J., “Magnetic resonance image segmentation using semi-automated software for quantification of knee articular cartilage—initial evaluation of a technique for paired scans,” Skeletal Radiology, Vol. 38, pp. 505-511, (2009).
[13]Emre H., Erkan U., Mumcuoglu., Hamcan S., “Automatic segmentation of human facial tissue by MRI–CT fusion: A feasibility study,” Computer Methods and Programs in Biomedicine ,Vol. 108, Issue 3, pp. 1106-1120, (2012).
[14]Sjoberg C., and Ahnesjo A., “Multi-atlas based segmentation using probabilistic label fusion with adaptive weighting of image similarity measures,” Computer Methods and Programs in Biomedicine , Vol. 110, Issue 3, pp.308-319, (2013).
[15]Stein D., Fritzsche K.H., Nolden M., Meinzer H.P., “The extensible open-source rigid and affine image registration module of the Medical Imaging Interaction Toolkit (MITK),” Computer Methods and Programs in Biomedicine, Vol. 100, Issue 1, pp.79-86, (2010).
[16]Dura E., Domingo J., Ayala G., Marti-Bonmati L. “Evaluation of the registration of temporal series of contrast-enhanced perfusion magnetic resonance 3D images of the liver,” Computer Methods and Programs in Biomedicine, Vol. 108, Issue 3,pp.932-945, (2012).
[17]Cheung M.R. and Krishnan K., “Using manual prostate contours to enhance deformable registration of endorectal MRI,” Computer Methods and Programs in Biomedicine, Vol. 108, Issue 1,pp. 330-337, (2012).
[18] http://medical.nema.org/
[19]Loresen W. E., Cline H. E., “Marching Cubes: A high resolution 3D surface construction algorithm,” Computer Graphics, Vol.21, pp. 163-169, (1987).
[20]Watt A., 3D Computer Graphics, Addison Wesley, Chapter 9, (1994).
[21] http://www.drabba.com/avizo/
[22] http://bmeimage.be.cycu.edu.tw/lab/database/CT/CT.html
[23]李文胜，CT 機的物理原理，工科物理，1997年。
[24]Zollikofer, C. P. E., and Ponce de Leon, M. S., “Tools for Rapid Prototyping in the Biosciences,” IEEE, C.G.&A., pp.48-55, (1995).
[25]呂信億，應用X 射線電腦斷層掃描於BGA 檢測之最佳化，清華大學動力機械工程所，2006。
[26]莊克士，醫學影像物理學，合記圖書出版社，2003。
[27] http://web.it.nctu.edu.tw/~hcsci/hospital/ins/mri.htm
[28]Lorensen, W. E., and Cline, H. E., “Marching cubes: A high resolution 3D surface construction algorithm, Computer Graphics, Vol. 21, No. 4, pp. 163-169, (1987).
[29]Schroeder W., “A topology modifying progressive decimation algorithm,” Proceedings of Visualization, IEEE Computer Society Press, Los Alamitos, CA, (1997).
[30]Schroeder W., Zarge J. and Lorensen W. E., “Decimation of triangle meshes,” Computer Graphics, pp. 65-70, (1992).
[31]Van A., Kim K., “Direct volume rendering with shading via three-Dimensional Textures”, University of California, Santa Cruz, (1996).

[32]Cameron G. & Undrill E., “Rendering volumetric medical image data on a SIMD-architecture Computer”, (1992).
[33]Wittenbrink M., Kim K., “Data dependent optimizations for permutation volume rendering”, Computer Systems Laboratory, HPL-97-59 , (1988).
[34]Foley J. D., Van D., Feiner A., & Hughes F., “Computer graphics: principles and practice , ” New York, (1990).
[35]Brown L. G.,” A survey of image registration techniques” ACM Computing Surveys , Vol. 24, No. 4, (1992).
[36]Besl P. J., and Mckay N. D., “A method for registration of 3-D shapes,” IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol. 14,No. 2, pp. 239-256, (1992).
[37]Pluim J., Maintz, J., and Viergever M., “Mutual-information-based registration of medical images a survey”, IEEE Trans. Med. Image, Vol. 22, No.8, (2003).
[38] Shannon C. E., “A mathematical theory of communication,” Bell Syst. Tech. J., Vol.27, pp623–656, (1948).
[39] Maes F., Collignon A., Vandermeulen D., Marchal G., and Suetens P., “Multimodality image registration by maximization of mutual information,” IEEE Trans. Med. Image., Vol. 16, pp.187–198, (1997).