基于Zernike相衬法的相衬成像探索

２０ｏ６年　工程图学学报　ＪＯＵＲＮＡＬ　ＯＦ　ＥＮＧＩＮＥＥＲＩＮＧ　ＧＲＡＰＨＩＣＳ　２００６　Ｎｏ．４　第４期　基于Ｚｅｒｎｉｋｅ相衬法的相衬成像探索　刁麓弘　，２一，　李　华　２，　于　斌４，　彭　翔４，　牛憨笨４　（１．中国科学院计算技术研究所智能信息处理重点实验室，北京１０００８０；　２．中国科学院计算技术研究所智能中心，北京１０００８０；３．中国科学院研究生院，北京１０００８０；　４．深圳大学光电子学研究所，广东深圳５１８０６０）　摘　要：传统ＣＴ由于利用人体组织对ｘ光的吸收进行成像具有很大的局限性。利　用相位衬度对物体成像可以获得更高的对比度和分辨率，因此相衬成像成为近年来研究的热　点。但是目前相衬图像的获取十分困难，一般需要利用第三代同步光源。为了摆脱物理手段　的限制和探索图像相位信息的意义，提出了一种基于Ｚｅｍｉｋｅ相衬法处理医学图像的方法，　并对该方法进行了分析。实验表明，结果图像呈现出了相衬成像的特点。　关键词：计算机应甩；医学图像处理；相衬成像；Ｚｅｍｉｋｅ相衬法　文章编号：１００３．０１　５８（２００６）０４．０１００．０５　中图分类号：ＴＰ　３９１　文献标识码：Ａ　Ｅｘｐｌｏｒａｔｉｏｎ　ｏｆ　Ｐｈａｓｅ　Ｃｏｎｔｒａｓｔ　Ｉｍａｇｉｎｇ　Ｂａｓｅｄ　ｏｎ　Ｚｅｒｎｉｋｅ　Ｐｈａｓｅ　Ｃｏｎｔｒａｓｔ　Ｍｅｔｈｏｄ　ＤＩＡＯ　Ｌｕ．ｈｏｎｇ　，２　，ＬＩ　Ｈｕａ　一，ＹＵ　Ｂｉｎ　，ＰＥＮＧ　Ｘｉａｎｇ　，ＮＩＵ　Ｈａｎ．ｂｅｎ　（１．ＫｅｙＬａｂｏｒａｔｏｒｙｏｆＩｎｔｅｌｌｉｇｅｎｔＩｎｆｏｒｍａｔｉｏｎＰｒｏｃｅｓｓｉｎｇ，Ｉｎｓｔｉｔｕｔｅ　ｏｆＣｏｍｐｕｔｉｎｇ　Ｔｅｃｈｎｏｌ０ｇ　ＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅ，Ｂｅｉｊｉｎｇ１０００８０，Ｃｈｉｎａ；　２．ＮＣＩＣ，ＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｎｃｅ，Ｂｅｉｅｊｉｎｇ１０００８０，Ｃｈｉｎａ；３．ＧｒａｄｕａｔｅＵｎｉｖｅｒｓｉｔｙ，ＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｎｃｅ，Ｂｅｉｅｊｉｎｇ１０００８０，Ｃｈｉｎａ；　４．Ｉｎｓｉｔｕｔｅ　ｔｏｆＯｐｔｏｅｌｅｃｔｒｏｎｉｃ，Ｓｈｅｎｚｈｅｎ　Ｕｎｉｖｅｒｓｉｙ，Ｓｈｅｎｚｈｅｎ　ｔＧｕａｎｇｄｏｎｇ　５１８０６０，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ｔｒａｄｉｔｉｏｎａｌ　ＣＴ　ｈａｓ　ａ　ｇｒｅａｔ　ｌｉｍｉｍｔｉｏｎ　ｂｅｃａｕｓｅ　ｉｔｓ　ｉｍａｇｅ　ｉＳ　ｐｒｏｄｕｃｅｄ　ｂａｓｅｄ　ｏｎ　ｔｈｅ　Ｘ－ｒａｙ　ａｂｓｏｒｂｅｄ　ｂｙ　ｈｕｍａｎ　ｔｉｓｓｕｅ．Ｐｈａｓｅ　ｃｏｎｔｒａｓｔ　ｉｍａｇｉｎｇ　ｃａｎ　ｇｅｔ　ｈｉｇｈ－ｒｅｓｏｌｕｔｉｏｎ　ａｎｄ　ｈｉ曲－ｃｏｎｔｒａｓｔ　ｉｍａｇｅｓ　ｂｙ　ｕｓｉｎｇ　ｔｈｅ　ｃｈａｎｇｅ　ｏｆ　ｐｈａｓｅｓ　ｔｏ　ｄｉｓｐｌａｙ　ｔｈｅ　ｏｂｊｅｃｔ．Ｔｈｅｎ　ｐｈａｓｅ　ｃｏｎｔｒａｓｔ　ｉｍａｇｉｎｇ　ｈａｓ　ｂｅｅｎ　ｐａｉｄ　ｍｕｃｈ　ａｔｔｅｎｔｉｏｎ　ｔｏ　ｉｎ　ｒｅｃｅｎｔ　ｙｅａｒｓ．ＨｏｗｅｖｅＬ　ｐｈａｓｅ－ｃｏｎｔｒａｓｔ　ｉｍａｇｅｓ’ａｃｑｕｉｓｉｔｉｏｎ　ｉｓ　ｖｅｒｙ　ｄｉｉｃｕｌｆｔ　ｎｏｗ．Ｉｎ　ｇｅｎｅｒａｌ，ｉｔ　ｎｅｅｄｓ　ｔｏ　ｂｅ　ｒｅｓｏｒｔｅｄ　ｔｏ　ｔｈｉｒｄ－ｇｅｎｅｒａｔｉｏｎ　ｓｙｎｃｈｒｏｔｒｏｎ　ｌｉｇｈｔ　ｓｏｕｒｃｅｓ．Ｉｎ　ｏｒｄｅｒ　ｔｏ　ｂｒｅａｋ　ｔｈｅ　ｌｉｍｉｔ　ｏｆ　ｐｈｙｓｉｃｓ　ｉｍｐｌｅｍｅｎｔａｔｉｏｎ　ａｎｄ　ｄｉｓｃｕｓｓ　ｔｈｅ　ｍｅａｎｉｎｇ　ｏｆ　ｉｍａｇｅｓ　ｐｈａｓｅ　ｉｎｆｏｒｍａｔｉｏｎ，ａ　ｍｅｔｈｏｄ　ｔｏ　ｐｒｏｃｅｓｓ　ｍｅｄｉｃａｌ　ｉｍａｇｅｓ　ｂａｓｅｄ　ｏｎ　ｈｅ　Ｚｅｍｉｔｋｅ　ｐｈａｓｅ－ｃｏｎｔｒａｓｔ　ｍｅｔｈｏｄ　ｉｓ　ｐｒｅｓｅｎｔｅｄ．Ｅｘｐｅｒｉｍｅｎｔｓ　ｓｈｏｗ　ｔｈａｔ　ｔｈｅ　ｐｒｏｄｕｃｅｄ　ｉｍａｇｅｓ　ｈａｖｅ　ｔｈｅ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ　ｐｈａｓｅ　ｃｏｎｔｒａｓｔ　ｉｍａｇｉｎｇ．　Ｋｅｙ　ｗｏｒｄｓ：ｃｏｍｐｕｔｅｒ　ａｐｐｌｉｃａｔｉｏｎ；ｍｅｄｉｃａｌ　ｉｍａｇｅ　ｐｒｏｃｅｓｓｉｎｇ；ｐｈａｓｅ　ｃｏｎｔｒａｓｔ　ｉｍａｇｉｎｇ；　Ｚｅｍｉｋｅ　ｐｈａｓｅ－ｃｏｎｔｒａｓｔ　ｍｅｔｈｏｄ　收稿日期：２００６—０３—２９　基金项目：国家自然科学基金重点资助项目（６０５３２０９０）：中国科学院计算所青年基金资助项目（２００５６６００—１７）　作者简介：刁麓弘（１９７８一），男，山东济南人，博士研究生，主要研究方向为计算机图形学与图像处理。　维普资讯 http://www.cqvip.com

第４期　刁麓弘等：基于Ｚｅｒｎｉｋｅ相衬法的相衬成像探索　ＣＴ作为可以探测人体内部信息的医疗设备　现在已经成为一种重要的辅助诊断手段。它根据　ｘ光穿过人体后的强度衰减，利用ｒａｄｏｎ逆变换　反求人体各组织器官的衰减系数，然后进行显　示。因此，对于衰减因子很小的器官和组织，传　统的ＣＴ图像很难清楚地对它们进行描述，这也　是ＣＴ成像的主要局限性之一。其他医学成像方　式，如ＭＲｊ、ＰＥＴ、超声等目前的分辨率和灵敏　度都相对比较低；而新发展的方法，如光学相干　层析成像，虽然精度较高，但穿透深度受到很大　限制。为了获得更高质量的图像，Ｄａｖｉｓ【１］和　Ｗｉｌｋｉｎｓ／２］分别于１９９５年和１９９６年发表了关于ｘ　射线相衬成像的奠基性工作。他们将相位信息引　入ｘ射线成像，为传统ｘ射线成像科学开辟了　一个崭新的研究领域Ｌ３Ｊ。光线在传播的过程中会　发生强度和相位的两种变化：强度变化是由传播　媒介的衰减因子造成的，而相位的改变则与媒介　的折射率有关。对于衰减因子很小的物质，由该　物质的折射率变化所引起的光线的相位变化幅　度仍然相对很大。一般认为，两者之间的差距会　达到３～５个数量级。利用相位衬度实现成像的最　大优点在于，即使物质对光线的吸收微弱到无法　探测，人们也有可能通过探测光线相位的变化来　实现对弱吸收物质的成像。这有望克服传统成像　方式中的图像对比度低，活体数据难以获得等困　难，因此也越来越成为当前研究的热点Ｌ４卜Ｌ７Ｊ。但　是目前ｘ射线相衬成像的理论与成像机制还不　成熟，对它的研究也大多停留于这两个阶　段嗍～【ｌ０Ｊ。因此，ｘ射线相衬图像的获取还十分困　难。图像处理中对于相衬成像的研究大都是针对　ＭＲ相衬图像的【¨Ｊ～【　Ｊ。　为了摆脱物理探测手段的限制，分析图像的　位相结构与相衬成像的关系，作者将傅立叶变换　与Ｚｅｍｉｋｅ相衬法相结合，提出了一种新的图像　处理方法并分析了处理后的结果。这种方法首先　利用傅立叶变换对医学图像的相位信息进行提　取，然后将获得的相位信息按频率进行滤波，最　后进行重建。实验表明：处理后的图像含有更加　丰富的信息并且出现了边缘增强的效果，这些均　符合相衬图像的特点，同时也说明合理改变图像　相位结构可以显示出更多的图像信息。　１　Ｚｅｒｎｉｋｅ相衬法　人眼和一般的成像系统都是根据光线的强度　成像的，光强的变化与振幅有关。因为位相物体　（只有空间的位相结构而透明度一样的物体，如　生物切片，油膜以及热塑等）只改变入射光的相　位而不影响其振幅。因此，为了显示位相物体，　Ｚｅｒｎｉｋｅ提出了一种将物的位相结构转换成像的　对比度的方法Ｌ１　。他使用一个使零频衰减同时相　移９０。的位相板进行空间滤波，使像平面上光强　分布近似正比于物体的位相分布。Ｚｅｒｎｉｋｅ相衬　法奠定了现代显微镜的基础，利用这种技术可以　直接观察生物切片中的活体组织，而无需对它们　染色。　位相物函数ｆ（ｘ）总是可以写成ｅ　的形　式（为简便本节所举的例子为一维函数，其中　是物体的坐标，　（　）是物体的位相分布）。当　（　）很小时，可以对函数做如下展开　ｆ（ｘ）＝１＋ｉｑ￣（ｘ）　若忽略相干光学的有限孔径的影响，则物的频谱　为　Ｆ（　）＝ｅ（１＋ｆ　））Ｐ　＝　（　）＋　（　）　其中　５（ｖ　）是狄拉克函数。Ｚｅｍｉｋｅ相板的光学　传递函数可写成　日　Ｖ　＝｛±　三　：　其中　，，是衰减因子，Ｖ０代表滤波的频谱平面上　的范围。由上可知，经滤波后物体的频谱变为　，’（１，　）＝Ｆ（ｖ　）　（１，）＝±ｆ　（１，　）＋ｉ￣ａ（ｖ　）　所以，物体的像函数可以写成　∞　…　、　ｆ’（　）＝Ｉ　Ｆ’（　）Ｐｕ　ｄ　＝±ｆ　＋ｆ　（　）　于是，像平面上可观察到的光强度为　ｌ（ｘ）＝ｆ’（　）　’木　）＝ｙ２±２　（　）＋一　（　）　厂’木　）是像函数厂’（　）的共轭。当　（】：）很小时，　维普资讯 http://www.cqvip.com

工程图学学报　２００６生　光强度可以近似看作与　（　）成ｉＥＬ￣，。　２该文对医学图像的处理　Ｚｅｍｉｋｅ相衬法的实现需要有相应的物理手　段的保证，这也给它的应用带来了限制。但是在　计算机处理中则没有这样的要求。本文设计了新　的传输函数，对原来的方法进行了改进，并充分　借鉴了其中的思想：通过数学变换突出了图像的　相位信息。　首先，对图像ｆ（ｘ，ｙ）做二维离散傅立叶变　换，提取图像的相位信息　Ｆ（　Ｖｙ）　赤　荟八　ｋ１　Ｎ　Ｎ　２）ＷＮｈⅣ＾和Ｎｗ是图像的高度和宽度。然后，设计一个　新的传输函数　（１，　，Ｖｙ），对图像进行处理　，　）＝月　，　）　，　）＝　ｌ砘，Ｖｙ）＋ｉＰｃ（ｂ－［ｖｘ，　，　Ｉ　盹，　）　＞Ｖｏ　传输函数Ｈ（ｖ　，Ｖｙ）的作用是将滤波范围１，０内频　谱函数的实部转化为虚部。最后将滤波后得到的　频谱函数做傅立叶逆变换重建新的图像　厂’（　），）　八　）：　ｌ　＂ｗ兰　，（　　３结果与分析　图ｌ至图４是利用第２节的方法对医学图像　处理的结果，左图均为原图，右图是处理过的图　像。滤波范围均选择为√１０。可以看出，处理后　的图像显示出更加丰富的信息，特别是包含衰减　因子较小的组织与器官的区域。另外，处理后的　图像也出现了边缘增强的效果。这些特点与相衬　成像的性质是相符的，这说明直接利用图像的相　位信息模拟相衬图像的生成是可行的。　经过分析可知：这种对图像处理的方法，其　本质是改变图像函数的对称部分和不对称部分　在低频分量的分布。图像函数ｆ（ｘ，Ｙ）可以写成　偶函数与奇函数的和　ｆ（ｘ，Ｙ）＝妄［厂（　，　）＋厂（一　，一　）】＋　去［厂（　，　）一厂（一　，一Ｙ］＝Ｐ（　，Ｙ）＋Ｄ（　，Ｙ）　偶函数ｅ（ｘ，Ｙ）和奇函数ｏ（ｘ，Ｙ）可以分别看作图　像的对称和不对称的部分。根据傅立叶变换的性　质可知，频谱函数Ｆ（ｖ　，Ｖｙ）的实部经傅立叶逆　变换后得到ｅ（ｘ，Ｙ），虚部的重构结果是ｏ（ｘ，Ｙ）。　设ｅ（ｘ，Ｙ）和ｏ（ｘ，Ｙ）在低频部分的分量分别为　ｅｚ（ｘ，Ｙ）和ｏｌ（　，Ｙ），对ｅｌ（　，Ｙ）和ｏｌ（　，Ｙ）作二　维离散傅立叶变换　巧（　，Ｖ　）＝　赢　（　（　，　）＋Ｄ』（　　），　　：　（　，Ｖｙ）＋ｉＯ￣（Ｖｘ，Ｖ　）　经传输函数Ｈ（ｖ　，Ｖｙ）滤波后，　’（Ｖｘ，１，　）＝　Ｆ（　，ｖｙ）Ｈ（ｖ￣，　）：ｉ（Ｅｚ（ｖ￣，　）＋０：ｒ（　，Ｖｙ））　对　’（１，　，Ｖｙ）进行逆傅立叶变换重构　’（　，Ｙ）＝　１＂∑Ｖｏ　窖　：　ｏｚ（ｘ，　）＋ｉｅ￣（ｘ，Ｊ，）　所以　＝　一　Ｉｆ，・Ｉ＝Ｄ』＋　一√　＝　＿ｅ　ｌ一　以上的实验结果与分析表明：改变图像的相位分　布，突出图像的相位信息，可以获得更多的图像　信息。　维普资讯 http://www.cqvip.com

第４期　刁麓弘等：基于ＺｅｒＭｋｅ相衬法的相衬成像探索　图１人体胸部的ＣＴ图像１　图２人体肺部的ＣＴ图像２　图３人体腹部ＣＴ图像１　图４人体腹部ＣＴ图像　维普资讯 http://www.cqvip.com

・１０４・　工程图学学报　２００６血　４未来的工作　相衬成像研究目前尚处在刚开始的阶段，作　者所作的工作主要是摆脱了物理探测手段的限　制，探索利用相位信息的生成具有相衬成像特点　图像的可行性。应该说，经过傅立叶变换提取的　相位信息与物体的相位分布的概念并不相同。如　何正确提取更加有用的相位信息，是下一步工作　的重点。另外，脱离了物理手段的限制，如何充　分利用数学变换的性质来构造较好的传输函数　和确定滤波范围，以取得更好的图像效果，也是　下一步的主要工作之一。　参考文献　Ｄａｖｉｓ　Ｔ　Ｊ，Ｇａｏ　Ｄ，Ｇｕｒｅｙｅｖ　Ｔ　Ｅ，ｅｔ　ａ１．Ｐｈａｓｅ－ｃｏｎｔｒａｓｔ　ｉｍａｇｉｎｇ　ｏｆ　ｗｅａｋｌｙ　ａｂｓｏｒｂｉｎｇ　ｍａｔｅｒｉａｌｓ　ｕｓｉｎｇ　ｈａｒｄ　Ｘ－ｒａｙｓ［Ｊ】．Ｎａｔｕｒｅ，１９９５，３７３：５９５—５９８．　Ｗｉｌｋｉｎｓ　Ｓ　Ｗ，Ｇｕｒｅｙｅｖ　Ｔ　Ｅ，Ｇａｏ　Ｄ，ｅｔ　ａ１．　Ｐｈａｓｅ—ｃｏｎｔｒａｓｔ　ｉｍａｇｉｎｇ　ｕｓｉｎｇ　ｐｏｌｙｃｈｒｏｍａｔｉｃ　ｈａｒｄ　Ｘ－ｒａｙｓ［Ｊ】．Ｎａｔｕｒｅ，１９９０，３８４：３３５—３３８．　Ｆｉｔｚｇｅｒａｌｄ　Ｒ．Ｐｈａｓｅ－ｓｅｎｓｉｔｉｖｅ　Ｘ－ｒａｙ　ｉｍａｇｉｎｇ［Ａ】．Ｉｎ：　Ｐｈｙｓｉｃｓ　Ｔｏｄａｙ［Ｃ】．２０００．２３—２６．　Ｎｕｇｅｎｔ　Ｋ　Ａ，Ｇｕｒｅｙｅｖ　Ｔ　Ｅ，Ｃｏｏｋｓｏｎ　Ｄ　Ｊ，ｅｔ　ａ１．　Ｑｕａｎｔｉｔａｔｉｖｅ　ｐｈａｓｅ　ｉｍａｇｉｎｇ　ｕｓｉｎｇ　ｈａｒｄ　Ｘ－ｒａｙｓ［Ｊ】．　Ｐｌａｙｓ．Ｒｅｖ．Ｌｅｔｔ．１９９６．７７：２９６１—２９６４．　Ｍｏｍｏｓｅ　Ａ，Ｆｕｋｕｄａ　Ｊ．Ｐｈａｓｅ　ｃｏｎｓｔｒａｓｔ　ｒａｄｉｏｇｒａｐｈｙ　ｏｆ　ｎｏｎｓｔａｉｎｅｄ　ｒａｔ　ｃｅｒｅｂｅｌｌｒａ　ｓｐｅｃｉｍｅｎ［Ｊ】．Ｍｅｄ．Ｐｈｙｓ．，　１　９９６，２２：３７５—３８０．　Ｃｌｏｅｔｅｎｓ　Ｐ，Ｌｕｄｗｉｇ　Ｗ，Ｂａｒｕｃｈｅｌ　Ｊ，ｅｔ　ａ１．　Ｈｏｌｏｔｏｍｏｇｒａｐｈｙ：Ｑｕａｎｔｉｔａｔｉｖｅ　ｐｈａｓｅ　ｔｏｍｏｇｒａｐｈｙ　ｗｉｔｈ　ｍｉｃｒｏｍｅｔｒｅ　ｒｅｓｏｌｕｔｉｏｎ　ｕｓｉｎｇ　ｈａｒｄ　ｓｙｎｃｈｒｏｔｒｏｎ　ｒａｄｉａｔｉｏｎ　Ｘ—ｒａｙｓ［Ｊ】．Ａｐｐ１．Ｐｈｙｓ．Ｌｅｔｔ．，１９９９，７５：　Ⅲ　２９１２－２９１４．　［７】　Ｃｈａｐｍａｎ　Ｄ，Ｔｈｏｍｌｉｎｓｏｎ　Ｗ，Ｊｏｈｎｓｔｏｎ　Ｒ　Ｅ，ｅｔ　ａ１．　Ｄｉｆｆｒａｃｔｉｏｎ　ｅｎｈａｎｃｅｄ　Ｘ－ｒａｙ　ｉｍａｇｉｎｇ［Ｊ］．Ｐｈｙｓ．Ｍｅｄ．　Ｂｉｏ１．，１９９７，４２：２０１５－２０２５．　［８】ＰａｇａｎｉｎＤ，Ｍａｙｏ　ＳＣ，ＧｕｒｅｙｅｖＴＥ，ｅｔａ／．Ｓｉｍｕｌｔａｎｅｏｕｓ　ｐｈａｓｅ　ｎａｄ　ａｍｐｌｉｕｔｄｅ　ｅｘｔｒａｃｔｉｏｎ　ｆｒｏｍ　ａ　ｓｉｎｇｌｅ　ｄｅｆｏｃｕｓｅｄ　ｉｍａｇｅ　ｏｆ　ａ　ｈｏｍｏｇｅｎｅｏｕｓ　ｏｂｊｅｃｔ［Ｊ】＿Ｊ．Ｍｉｃｒｏｓｃｏｐｙ，　２００２，２０６：３３—４０．　【９】Ｍｏｍｏｓｅ　Ａ．Ｎｕｃ１．Ｉｎｓｔｒｕｍ［Ｊ］．Ｍｅｔｈ．Ａ．，１９９５，３５２：　６２２．　［１０】Ｐａｓｃｈａｌ　ＣＢ，Ｃａｒｒｏｌｌ　ＦＥ，Ｗｏｒｒｅｌｌ　ＪＡ，ｅｔ　ａ１．Ｖｏｌｕｍｅｔｒｉｃ　ｍｏｎｏｃｈｒｏｍａｔｉｃ　Ｘ－ｒａｙ　ｔｏｍｏｇｒａｐｈｙ　ｏｆｔｈｅ　ｌｕｎｇｓ［Ａ］．Ｉｎ：　Ｃｏｍｐｕｔｅｒ　Ｔｏｍｏｇｒａｐｈｙ［Ｃ】．２００３．５—１０．　【ｌ１】　Ｗｉｎｋ　Ｏ，Ｎｉｅｓｓｅｎ　Ｗ　Ｊ，Ｖｉｅｒｇｅｒｖｅｒ　Ｍ　Ａ．Ｆａｓｔ　ｄｅｌｉｎｅａｔｉｏｎ　ａｎｄ　ｖｉｓｕａｌｉｚａｔｉｏｎ　ｏｆ　ｖｅｓｓｅｌｓ　ｉｎ　３－Ｄ　ｎａｇｉｏｇｒａｐｈｉｃ　ｉｍａｇｅｓ［Ｊ】．ＩＥＥＥ　Ｔｒａｎｓａｃｔｉｏｎｓ　ｏｎ　Ｍｅｄｉｃａｌ　Ｉｍａｇｉｎｇ，２０００，１　９（４）：３３７—３４６．　［１２】　Ｗｏｎｇ　Ａ　Ｌ　Ｎ，Ｌｉｕ　Ｈｕａｆｅｎｇ；Ｓｈｉ　Ｐｅｎｇｃｈｅｎｇ．　Ｓｅｇｍｅｎｔａｔｉｏｎ　ｏｆ　ｍｙｏｃａｒｄｉｕｍ　ｕｓｉｎｇ　ｖｅｌｏｃｉｔｙ　ｆｉｅｌｄ　ｃｏｎｓｔｒａｉｎｅｄ　ｆｒｏｎｔ　ｐｒｏｐａｇａｔｉｏｎ［Ａ】．Ｉｎ：Ｐｒｏｃｅｅｄｉｎｇｓ．　Ｓｉｘｔｈ　ＩＥＥＥ　Ｗｏｒｋｓｈｏｐ　ｏｎ　Ａｐｐｌｉｃａｔｉｏｎｓ　ｏｆ　Ｃｏｍｐｕｔｅｒ　Ｖｉｓｉｏｎ，２００２．（ＷＡＣＶ　２００２）［Ｃ】．２００２．８４—８９．　［１　３】Ｂｏｇｒｅｎ　Ｈ　Ｇ　Ｂｕｏｎｏｃｏｒｅ　Ｍ　Ｈ．Ｆｏｕｒ－Ｄ　ｂｌｏｏｄ　ｆｌｏｗ　ｐａｔｔｅｒｎｓ　ｉｎ　ｔｈｅ　ｔｈｏｒａｃｉｃ　ａｏｒｔａ　ｉｎ　ｐａｔｉｅｎｔｓ　ｗｉｈｔ　ＣＡＤ　ｃｏｍｐａｒｅｄ　ｔｏ　ａｇｅ－ｍａｔｃｈｅｄ　ｎｏｒｍａｌ　ｓｕｂｊｅｃｔｓ［Ａ】．Ｉｎ：　Ｐｒｏｃｅｅｄｉｎｇｓ．Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｗｏｒｋｓｈｏｐ　ｏｎ　Ｍｅｄｉｃａｌ　Ｉｍａｇｉｎｇ　ｎａｄＡｕｇｍｅｎｔｅｄ　Ｒｅａｌｉｔｙ［Ｃ】．２００１．７２—７６．　［１　４】Ｆａｔｏｕｒａｅｅ　Ｎ，Ａｍｉｎｉ　Ａ　Ａ．Ｒｅｇｕｌａｒｉｚａｔｉｏｎ　ｏｆ　ｆｌｏｗ　ｓｔｒｅａｍｌｉｎｅｓ　ｉｎ　ｍｕｌｔｉｓｌｉｃｅ　ｐｈａｓｅ—ｃｏｎｔｒａｓｔ　ＭＲ　ｉｍａｇｉｎｇ［Ｊ】．ＩＥＥＥ　Ｔｒｎａｓａｃｔｉｏｎｓ　ｏｎ　Ｍｅｄｉｃａｌ　Ｉｍａｇｉｎｇ，　２００３，２２（６）：６９９—７０９．　【１５】母国光，战元令．光学［Ｍ］．北京：人民教育出版社，　】９７８．３９５—３９７．