基于非均匀傅里叶变换的超声层析成像研究

【摘要】：超声波成像技术作为现代医学四大影像技术之一,有着不可替代的作用。X光成像可以准确获得人体组织的信息,是一种很好的检测手段,但其对软组织的成像不理想,且因其有电离辐射特性,对于人体的影响较大,所以它的应用有局限性。超声成像仪器因其价格便宜,无电离辐射,可以制造成便携式设备等特点,在医院的使用频率很高。常用的B超能应用于常规性的体检或胎儿的检查等情况。但是B超的成像较模糊,对人体情况的判断需要依靠医生的经验,个人的主观性会影响判断的结果,因此超声波用于人体的定量的成像成为一个重要的研究方向。层析成像(Computed Tomography)是通过从许多不同方向照射物体,根据透射或反射的数据来重建断面图像的成像技术。 本文讨论的是考虑散射的情况下采用超声波这种安全的能源对人体组织进行定量成像的超声衍射层析成像问题。所完成的主要工作和创新点如下: 基于傅里叶散射投影定理,频域插值(Gridding)算法是在频域对变换后的散射场数据进行插值,得到均匀网格上的值,从而进行二维傅里叶反变换得到介质成像平面的分布。 非均匀傅里叶变换(Non-uniform Fourier Transform)包括前向NUFT和后向NUFT。其中前向NUFT的快速算法可以采用对过采样的FFT变换的频域值进行插值得到。min-max方法对其插值参数进行优化,得到基于min-max的快速非均匀傅里叶变换的算法。本文讨论的是后向NUFT,即二维情况下从非均匀样点到均匀样点的前向NUFT的反变换问题。因为计算的规模很大,不能采用直接求解逆矩阵的方法实现,所以采用最小二乘法的思想,采用迭代算法(CG)逐步接近最优化解。采用插值方法的成像结果作为迭代的初始值,可以减少迭代次数,减小运算量。经试验分析,采用本文的方法,成像效果得到提高。 对仿真数据的获取方法进行了讨论分析,在弱散射情况下,基于傅里叶散射投影定理,可以对Shepp-Logan模型的频域任一点的值进行解析求解,基于此可以得到不引入任何误差的用于仿真的数据。对多频率投影的情况,分析其频域分布特点,在多频率入射的情况下,可以减少成像的时间,在医学图像的提取中是一个好的发展方向。在成像过程中会出现不完整投影的情况,分析其不完整程度对成像效果的影响,和数据冗余以减少投影或采样点,使信息的利用率变高。对透射和反射模式下散射场所提供的信息量进行分析,在不增大采样时间的情况下可同时获取反射场和透射场的数据,将其进行结合,可以提高成像质量。