基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法技术

一种基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法，通过对待测冷冻电镜单颗粒图片集进行功率谱高通滤波器参数估计，然后计算待配准图像的傅里叶功率谱，并提取特征点，根据特征点对待配准图像进行反旋转后对图像进行矫正；基于傅里叶功率谱的旋转周期性，再次根据特征点计算另一相对旋转角后进行相对反旋转和相对矫正，最后根据矫正以及相对矫正得到的矫正图像与目标图像进行相似性计算，取相似性最高的矫正图像作为配准参数结果。本发明专利技术使用傅里叶功率谱域的分析对噪声、颗粒以及背景进行了有效分离，同时基于快速傅里叶变换的有效实现大大提高了算法的计算性能，无需进行迭代修正可直接计算出配准参数，使得配准无论从计算速度还是精度与鲁棒性上都得到了提升。

【技术实现步骤摘要】
基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法
本专利技术涉及的是一种图像处理领域的技术，具体是一种基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法。
技术介绍
在冷冻电镜的检测过程中，为了保证生物大分子的活性尽量不受电子辐射影响，电子显微镜使用的电子剂量必须保证在一个很低的程度下，导致冷冻电镜单颗粒图像的噪声很大，信噪比极低(通常低至0.1以下)。冷冻电镜图像的分析技术包括二维分析技术与三维分析技术，二维分析技术主要对图像进行聚类来降低噪声获取较为清晰的电镜图像，而三维分析技术则通过冷冻电镜图像重建其三维结构。无论是二维分析技术还是三维分析技术都需要在实施过程中对冷冻电镜图像进行配准以便于更好计算图像距离、图像相似性以及相对方向关系，且最终结果严重依赖于配准算法的准确性。冷冻电镜图像配准的主要难点在于图像极低的信噪比以及颗粒部分只占据图像的一部分使得颗粒需要与背景进行分离。目前的冷冻电镜配准算法主要还是迭代算法，这样的方法无法保证收敛速度，使得算法非常耗时，且对噪声的抗性也不强。配准算法是整个冷冻电镜分析中的重要一环，也是重要耗时环节，因此目前需要一种全新的配准算法来对整个过程进行加速并且能保证算法的精度与鲁棒性。
技术实现思路
本专利技术针对现有技术存在的上述不足，提出一种基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法，使用傅里叶功率谱域的分析对噪声、颗粒以及背景进行了有效分离，同时基于快速傅里叶变换的有效实现大大提高了算法的计算性能，无需进行迭代修正可直接计算出配准参数，使得配准无论从计算速度还是精度与鲁棒性上都得到了提升。本专利技术是通过以下技术方案实现的：本专利技术涉及一种基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法，通过对待测冷冻电镜单颗粒图片集进行功率谱高通滤波器参数估计，然后计算待配准图像的傅里叶功率谱，并提取特征点，根据特征点对待配准图像进行反旋转后对图像进行矫正；基于傅里叶功率谱的旋转周期性，再次根据特征点计算另一相对旋转角后进行相对反旋转和相对矫正，最后根据矫正以及相对矫正得到的矫正图像与目标图像进行相似性计算，取相似性最高的矫正图像作为配准参数结果。所述的功率谱高通滤波器参数估计是指：利用配准算法在随机挑选的小样本数据集上进行配准精度评估，选取最佳的功率谱高通滤波器参数，即功率谱域滤波半径r作为后续完整配准过程的方法参数。所述的配准精度评估是指：将小样本数据集上的所有数据向指定图像进行配准，并计算所有配准过后的图像与指定图像的相似度的均值作为精度指标，该精度指标采取相关熵其中：P，Q分别为两幅冷冻电镜单颗粒图像，m，n分别为其高度与宽度，e为自然对数，Pij，Qij为两幅图像位于(i，j)处的像素值，σ为给定参数，一般设为1。所述的待配准图像的傅里叶功率谱P(u，v)＝|F(u，v)|2，其中：F(u，v)是图像的傅里叶变换u，v是频域坐标，|·|表示取幅值操作，P(i，v)是图像的功率谱，对于两张需要进行配准的冷冻电镜单颗粒图像I1、I2，分别计算其傅里叶功率谱并根据估计所得参数进行高通滤波得到对应功率谱P1、P2。所述的提取特征点是指：对两个傅里叶功率谱P1、P2，寻找满足条件的局部极大值，计算其邻域的重心作为特征点，其中：领域为位于局部极大值邻近区域且幅度不低于局部极大值的一定比例的傅里叶功率谱点集，邻域的重心其中：U(m)为局部最大值的邻域，pij为位于(i，j)处的像素值，gx、gy为重心的横纵坐标。所述的反旋转是指：对于两张需要进行配准的冷冻电镜单颗粒图像，根据特征点的相对方位，计算特征点的相对旋转角度α0，并据此对一张冷冻电镜单颗粒图像进行反旋转，得到两张仅存在平移变换的冷冻电镜单颗粒图像。所述的矫正是指：通过傅里叶空间的平移配准计算两张反旋转后的图片的平移参数Δx0，Δy0，并根据α0，Δx0，Δy0对图像进行矫正得矫正后图像I′1，其中：傅里叶空间的平移配准为F1(u，v)和F2(u，v)分别是两张图像的傅里叶变换，u，v是频域坐标，|·|表示取幅值操作，*表示复数共轭，表示傅里叶反变换，δ(x+Δx0，y+Δy0)是位于(Δx，Δy)处的冲激函数，通过求解冲击函数最大值位置即可找到两张冷冻电镜单颗粒图像间的平移量，通过平移实现矫正。所述的相对反旋转是指：对于两张需要进行配准的冷冻电镜单颗粒图像，选取相对旋转角为α1＝(α0+180°)％360°，并据此对一张冷冻电镜单颗粒图像进行反旋转，得到两张仅存在平移变换的冷冻电镜单颗粒图像，其中：％是指取余运算。所述的相对矫正是指：通过傅里叶空间的平移配准计算两张图片的平移参数Δx1，Δy1，并根据α1，Δx1，Δy1对图像进行矫正得矫正后图像I″1。所述的相似性计算是指：计算过两张矫正图像I′1、I″1与图像I2的相似性，采取相似性度量为相关熵。所述的配准参数结果，即选取相似性高的图像输出其配准参数α，Δx，Δy。本专利技术涉及一种实现上述方法的系统，包括：特征提取单元、结果计算单元以及结果决策单元，其中：特征提取单元指傅里叶功率谱特征提取与精修部分，对输入信息进行处理，与特征精修单元相连并传输特征信息；结果计算单元指根据特征信息计算多种可行配准结果，与决策单元相连并传输多种结果信息；结果决策单元对多种结果信息进行决策，最终输出最佳结果。技术效果本专利技术整体解决了现有技术的配准计算效率低下的问题，本专利技术通过傅里叶变换与功率谱分析算法对颗粒、噪声与背景进行分离，进而对颗粒进行方位配准，减小了噪声与背景对配准过程的影响，提升了配准算法的鲁棒性。采取快速傅里叶变换作为配准算法的基本步骤，同时无需迭代计算，可以显著提高算法的计算效率，减少配准时间。附图说明图1为本专利技术流程图；图2为实施例GroEL真实冷冻电镜单颗粒图像样本示意图；图3为实施例GroEL真实结构示意图；图4为实施例配准结果以及生成的参考图像；图5为实施例配准算法进行聚类的类平均结果示意图；图6为实施例配准算法进行重构的三维模型结果示意图。具体实施方式本实施例中傅里叶变换采取快速傅里叶变换，图像旋转平移变换中采取双线性插值，旋转中心均为图像中心。本实施例的配准算法对参考图像生成方法无特殊要求，此实例仅给出数据集配准中一种生成方法以便对整个过程能有更清晰的认识。如图2所示，为GroEL蛋白的真实冷冻电镜单颗粒图像。该数据集为4096张128×128的GroEL蛋白的冷冻电镜单颗粒图像。GroEL是一个具有D7对称性的结构，配准的目标是讲所有图像矫正到同一方位上。后续聚类的目标是将数据集中的单颗粒图像聚类到32个类别中，重构的目标是生成GroEL的三维结构。如图1所示，为本实施例涉及一种基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法，具体包括以下步骤：第一步，随机选取GroEL数据集中的100张图像作为子数据集进行高通滤波器参数估计，以子数据集第一张图像作为参考图像，应用一组高通滤波器参数本文档来自技高网...

【技术保护点】
1.一种基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法，其特征在于，通过对待测冷冻电镜单颗粒图片集进行功率谱高通滤波器参数估计，然后计算待配准图像的傅里叶功率谱，并提取特征点，根据特征点对待配准图像进行反旋转后对图像进行矫正；基于傅里叶功率谱的旋转周期性，再次根据特征点计算另一相对旋转角后进行相对反旋转和相对矫正，最后根据矫正以及相对矫正得到的矫正图像与目标图像进行相似性计算，取相似性最高的矫正图像作为配准参数结果。/n

【技术特征摘要】
1.一种基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法，其特征在于，通过对待测冷冻电镜单颗粒图片集进行功率谱高通滤波器参数估计，然后计算待配准图像的傅里叶功率谱，并提取特征点，根据特征点对待配准图像进行反旋转后对图像进行矫正；基于傅里叶功率谱的旋转周期性，再次根据特征点计算另一相对旋转角后进行相对反旋转和相对矫正，最后根据矫正以及相对矫正得到的矫正图像与目标图像进行相似性计算，取相似性最高的矫正图像作为配准参数结果。


2.根据权利要求1所述的基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法，其特征是，所述的功率谱高通滤波器参数估计是指：利用配准算法在随机挑选的小样本数据集上进行配准精度评估，选取最佳的功率谱高通滤波器参数，即功率谱域滤波半径r作为后续完整配准过程的方法参数。


3.根据权利要求1所述的基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法，其特征是，所述的配准精度评估是指：将小样本数据集上的所有数据向指定图像进行配准，并计算所有配准过后的图像与指定图像的相似度的均值作为精度指标，该精度指标采取相关熵其中：P，Q分别为两幅冷冻电镜单颗粒图像，m，n分别为其高度与宽度，e为自然对数，Pij，Qij为两幅图像位于(i，j)处的像素值，σ为给定参数。


4.根据权利要求1所述的基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法，其特征是，所述的待配准图像的傅里叶功率谱P(u，v)＝|F(u，v)|2，其中：F(u，v)是图像的傅里叶变换u，v是频域坐标，|·|表示取幅值操作，P(u，v)是图像的功率谱，对于两张需要进行配准的冷冻电镜单颗粒图像I1、I2，分别计算其傅里叶功率谱并根据估计所得参数进行高通滤波得到对应功率谱P1、P2。


5.根据权利要求1所述的基于傅里叶功率谱的冷冻电镜单颗粒图像配准方法，其特征是，所述的提取特征点是指：对两个傅里叶功率谱P1、P2，寻找满足条件的局部极大值，计算其邻域的重心作为特征点，其中：领域为位于局部极大值邻近区域且幅度不低于局部极大值的一定比例的傅里叶功率谱点集，邻域的重心其中：U(m)为局部最大值的邻域，pij为位于(i，j)处的像素值，gx、gy为重心的横纵坐标。


6.根据权利要求1所述的...

【专利技术属性】
技术研发人员：陈宇轩，沈红斌，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：上海;31