一种基于压缩感知的光声图像重建方法
【专利摘要】本发明公开了一种基于压缩感知的光声图像重建方法，包括以下步骤：在目标组织周围放置M个传感器；用宽束短脉冲激光周期性照射组织W次，各传感器采集光声信号；对各传感器采集的W次光声信号求平均，得到M组平均光声信号；对每组平均光声信号做FFT，得到M组频谱数据；设定M个随机向量，根据各随机向量对各组频谱数据进行N点随机抽样；选取正交基矩阵，根据传感器脉冲响应及随机抽样得到的频率值计算测量矩阵；利用压缩感知求解目标生物组织的光吸收分布，得到光声重建图像。本发明结合压缩感知，利用少量传感器采集的光声信号完成目标生物组织内部的三维图像重建，具有传感器数量少、传感器摆放位置相对任意等优点。
【专利说明】一种基于压缩感知的光声图像重建方法
【技术领域】
[0001]本发明涉及图像处理领域，特别是一种基于压缩感知的光声图像重建方法。
【背景技术】
[0002]光声成像以超声波作为信息的载体，反映生物组织的空间光学吸收分布特性，它同时具备了纯光学成像和超声成像二者的优点，在保证一定成像深度的同时又具有高的对比度和空间分辨率，尤其适合于对生物软组织的成像。光声成像可根据待测目标的光学吸收特性选取合适的激励光波长，可定量的测量与光学吸收相关的生理参数，并可在完全无损伤的情况下对生物的高级功能活动进行观察并提供各方向高分辨率和高对比度的组织断层图像。尽管光声成像具备以上优势，但由于使用单传感器进行圆周扫描采集信号周期长或多传感器阵列受制于传感器的带宽、阵元数及间距而制造价格昂贵，该技术若要步入到临床实践阶段，仍需做出大量研究和改善。
[0003]压缩感知技术是一种新的采样理论，其认为只要信号是可压缩的或在某个变换域中是稀疏的，就可以利用随机采样直接得到压缩数据，然后通过求解最优化问题从采集到的少量压缩数据中重建出原始信号。利用压缩感知技术可以大大减少光声重建需要采集的数据，从而可以利用少量观测角度采集到的信息或较少数量传感器采集到的信息进行光声重建，大大节省传感器采样时间和制造成本。
【发明内容】
[0004]发明目的:本发明所要解决的技术问题是针对传统光声重建方法信号采样周期长或传感器阵列制造成本高的问题，提供一种基于压缩感知的光声图像重建方法。
[0005]为了解决上述技术问题，本发明公开了一种基于压缩感知技术的光声图像重建方法，包括以下步骤:
[0006]步骤一，在目标组织周围放置M个传感器；
[0007]步骤二，用宽束短脉冲激光周期性照射组织W次，各传感器采集光声信号；
[0008]步骤三，对各传感器采集的W次光声信号求平均，得到M组平均光声信号；
[0009]步骤四，对每组平均光声信号做FFT，得到M组频谱数据；
[0010]步骤五，设定M个随机向量，根据各随机向量对各组频谱数据进行N点随机抽样；
[0011]步骤六，选取正交基矩阵，根据传感器脉冲响应及随机抽样得到的频率值计算测量矩阵；
[0012]步骤六，利用压缩感知求解目标生物组织的光吸收分布，得到光声重建图像。
[0013]本发明中，优选地，所述步骤一中M个传感器可以放置在目标生物组织周围的任何位置，比如放置于组织一侧的矩形网格中的任意网格点上，或放置在包围组织的整个球体网格中的任意网格点上，M的取值与所要重建光声图像经过正交基矩阵变换后的图像的稀疏度有关，图像越稀疏，所需要的传感器个数就越少，通常可选择在30~150之间。
[0014]本发明中，优选地，所述步骤二中激光发射器周期性的发射宽束短脉冲激光，照射整个生物组织。每发射一次激光，所有传感器就采集组织发出的光声信号，重复发射W次，接收过程中不需改变传感器位置。
[0015]本发明中，步骤三对每个传感器采样得到的W组光声信号求平均，得到M组平均光声信号，以减少随机噪声干扰，其中W的值太小重建图像受噪声影响较大，W的值太大将耗费大量信号采样时间，通常可选择在I~128之间。
[0016]本发明中，优选地，所述步骤四中需要对每组平均光声信号进行快速傅里叶变换(FFT)，得到对应的频谱数据。
[0017]本发明中，优选地，所述步骤五中针对不同传感器计算得到的频谱数据产生不同的随机向量对该频谱数据进行N点随机抽样，将M组抽样得到的N点频谱数据合并为一个M行N列的频谱抽样矩阵Y，其中N即为重建图像宽度。 [0018]本发明中，优选地，所述步骤六中正交基矩阵Φ应根据需要重建图像在该正交基变换后的系数是稀疏的或可压缩的为前提进行选择，通常可以是傅里叶正交基、小波正交基、曲波正交基中的一种；测量矩阵K可根据传感器和声源的相对位置、随机抽样得到的频率值及传感器脉冲响应计算得到。
[0019]本发明中，优选地，所述步骤七中利用压缩感知技术，根据随机抽样得到的频谱抽样矩阵Y、用于稀疏信号表示的正交基矩阵Φ以及计算得到的测量矩阵K，通过求解最优化问题，求解得到目标生物组织的空间光吸收分布。
【专利附图】
【附图说明】
[0020]下面结合附图和【具体实施方式】对本发明做更进一步的具体说明，本发明的上述和/或其他方面的优点将会变得更加清楚。
[0021]图1是本发明方法的流程图。
[0022]图2是本发明方法的设备安放示意图。
【具体实施方式】:
[0023]本发明结合压缩感知技术，通过随机放置的少数传感器采集光声信号，并转换到频域进行随机抽样，利用随机抽样得到的频谱数据、计算得到的测量矩阵以及用于稀疏化重建图像的正交基矩阵，通过求解一个最优化过程重建出生物组织的空间光吸收分布。
[0024]如图1所示，本发明公开了一种基于压缩感知的光声图像重建方法，包括以下步骤:
[0025]步骤一，在目标组织周围放置M个传感器；
[0026]步骤二，用宽束短脉冲激光周期性照射组织W次，各传感器采集光声信号；
[0027]步骤三，对各传感器采集的W次光声信号求平均，得到M组平均光声信号；
[0028]步骤四，对每组平均光声信号做FFT，得到M组频谱数据；
[0029]步骤五，设定M个随机向量，根据各随机向量对各组频谱数据进行N点随机抽样；
[0030]步骤六，选取正交基矩阵，根据传感器脉冲响应及随机抽样得到的频率值计算测量矩阵；
[0031]步骤六，利用压缩感知求解目标生物组织的光吸收分布，得到光声重建图像。
[0032]本发明中，步骤一，将M(30~150)个传感器安放在目标生物组织周围的任何位置，比如放置于组织一侧的矩形网格中的任意网格点上，也可以放置在包围组织的整个球体网格中的任意网格点上。一种安放激光发生器、生物组织、球体网格传感器阵列等装置的示意图如图2所示，其中透镜用来将较窄的激光束进行扩束，以覆盖到整个生物组织。
[0033]本发明中，步骤二，激光发射器周期性的发射宽束短脉冲激光，照射整个生物组织。每发射一次激光，所有传感器就采集组织发出的光声信号，重复发射w(l~128)次，接收过程中不需改变传感器位置。
[0034]本发明中，步骤三对每个传感器采样得到的W组光声信号求平均，得到M组平均光
声信号Pm(rm，t）:
【权利要求】
1.一种基于压缩感知的光声图像重建方法，其特征在于，包括以下步骤: 步骤一，在目标组织周围放置M个传感器； 步骤二，用宽束短脉冲激光周期性照射组织W次，各传感器采集光声信号； 步骤三，对各传感器采集的W次光声信号求平均，得到M组平均光声信号； 步骤四，对每组平均光声信号做FFT，得到M组频谱数据； 步骤五，设定M个随机向量，根据各随机向量对各组频谱数据进行N点随机抽样； 步骤六，选取正交基矩阵，根据传感器脉冲响应及随机抽样得到的频率值计算测量矩阵； 步骤七，利用压缩感知求解目标生物组织的光吸收分布，得到光声重建图像。
2.根据权利要求1所述的一种基于压缩感知的光声图像重建方法，其特征在于，所述步骤一中放置传感器，可以任意放置在位于生物组织一侧的矩形网格中，或放置在包围生物组织的球体网格中。
3.根据权利要求1所述的一种基于压缩感知的光声图像重建方法，其特征在于，所述步骤二激光发生器每发射一次激光，M个传感器同时采集一段光声信号，并重复W次，测量中，不需改变传感器位置。
4.根据权利要求1所述的一种基于压缩感知的光声图像重建方法，其特征在于，所述步骤三分别对每个传感器采样得到的W组光声信号求平均，得到M组平均光声信号。
5.根据权利要求1所述的一种基于压缩感知的光声图像重建方法，其特征在于，所述步骤四需要分别对每组平均光声信号进行快速傅里叶变换(FFT)，得到对应的频谱数据。
6.根据权利要求1所述的一种基于压缩感知的光声图像重建方法，其特征在于，所述步骤五需要分别为各组频谱数据设定随机向量并对其进行N点随机抽样，将M组频谱抽样结果合并为一个MXN大小的频谱抽样矩阵Y。
7.根据权利要求1所述的一种基于压缩感知的光声图像重建方法，其特征在于，所述步骤六中选择的正交基矩阵Φ可以是傅里叶正交基、小波正交基、曲波正交基中的一种，也可以选取其它正交基矩阵；所述步骤六中测量矩阵K可根据传感器和声源的位置、随机抽样得到的频率值及传感器脉冲响应计算得到。
8.根据权利要求1所述的一种基于压缩感知的光声图像重建方法，其特征在于，所述步骤七中利用压缩感知理论，使用压缩感知信号重建时的最优化求解方法，结合选择的正交基Φ、计算得到的测量矩阵K以及随机抽样得到的频谱抽样矩阵Y，求解得到生物组织的光吸收分布，即为光声图像的重建结果。
【文档编号】A61B5/00GK103584835SQ201310436669
【公开日】2014年2月19日 申请日期:2013年9月24日 优先权日:2013年9月24日
【发明者】袁杰, 邵真天, 夏煦菁, 封婷, 朱毅, 李文超, 张星, 储哲琪, 温馨, 顾鹏, 沈庆宏, 都思丹 申请人:南京大学