1.一种荧光多分子定位方法,其特征在于,包括如下步骤:采集多个荧光分子通过光学系统的荧光图像;对所述荧光图像进行傅里叶变换,获得所述荧光图像的傅里叶频谱;对所述荧光图像的傅里叶频谱进行解卷积运算,并根据解卷积之后的图像构造相应的傅里叶字典;基于所述傅里叶字典和解卷积之后的图像,根据正交匹配追踪算法计算多个荧光分子的坐标位置;所述采集多个荧光分子通过光学系统的荧光图像的步骤中,所述荧光图像为:,其中Nq为单个荧光分子发出的光子数,σ为高斯函数的标准差,(x,y)为图像上的坐标,(xq,yq)为荧光分子的坐标位置,b为噪声的强度;所述对所述荧光图像进行傅里叶变换,获得所述荧光图像的傅里叶频谱的步骤中,所述荧光图像的傅里叶频谱为:其中,M,N为图像的像元数,F(k,l)为光学系统在像元(k,l)处的光学传递函数,B[k,l]为像元(k,l)处的噪声的傅里叶变换;所述对所述荧光图像的傅里叶频谱进行解卷积运算,并根据解卷积之后的图像构造相应的傅里叶字典的步骤包括:基于方程对所述荧光图像的傅里叶频谱进行解卷积运算,其中λ为根据图像的噪声大小预设的系数;将解卷积之后的图像矩阵转换为向量形式,并根据图像大小构造相应的傅里叶字典,得到g=Av,其中g是由解卷积之后的图像,A为傅里叶字典,v是由Nq组成的向量元素;基于所述傅里叶字典和解卷积之后的图像,根据正交匹配追踪算法计算多个荧光分子的坐标位置的步骤包括:对参数进行初始化,另残差r0=g,被选原子索引集 φ为空集,迭代次数i=0;基于方程 找出索引值ti,并添加索引值到被选原子索引集 其中d为矩阵A的列数;计算由被选原子 展成的正交投影Pi;计算出待求的近似解与残差:ai=Pig,ri=g-ai;判断i是否小于等于K,若是,则另i=i+1,并继续基于方程找 出索 引值 后 添加 索 引值 到被 选原 子 索引 集若否,计算在 处的 近似值, 并将当前的vi作为荧光分子的的坐标向量,其中K为预估分子数, 为第K次迭代得到的索引集,为定位后的像素值。2.一种超分辨成像方法,其特征在于,包括如下步骤:采集多幅由多个荧光分子通过光学系统的荧光图像;采用如权利要求1所述的荧光多分子定位方法对每幅荧光图像中的荧光分子进行定位,得到其坐标位置;对定位后的所有荧光图像进行累加,根据全部荧光分子的坐标位置描绘出超分辨荧光图像。3.一种荧光多分子定位装置,其特征在于,包括:采集模块,用于采集多个荧光分子通过光学系统的荧光图像;傅里叶变换模块,用于对所述荧光图像进行傅里叶变换,获得所述荧光图像的傅里叶频谱;解卷积模块,用于对所述荧光图像的傅里叶频谱进行解卷积运算,并根据解卷积之后的图像构造相应的傅里叶字典;定位模块,用于基于所述傅里叶字典和解卷积之后的图像,根据正交匹配追踪算法计算多个荧光分子的坐标位置;所述解卷积模块包括:解卷积单元,用于基于方程对所述荧光图像的傅里叶频谱进行解卷积运算,其中λ为根据图像的噪声大小预设的系数;转换单元,用于将解卷积之后的图像矩阵转换为向量形式,并根据图像大小构造相应的傅里叶字典,得到g=Av,其中g是由解卷积之后的图像,A为傅里叶字典,v是由Nq组成的向量元素;所述定位模块包括:初始化单元,用于对参数进行初始化,另残差r0=g,被选原子索引集 φ为空集,迭代次数i=0;搜寻单元,用于基于方程 找出索引值ti,并添加索引值到被选原子索引集 其中d为矩阵A的列数;第一计算单元,用于计算由被选原子 展成的正交投影Pi;第二计算单元,用于计算出待求的近似解与残差:ai=Pig,ri=g-ai;判断定位单元,用于判断i是否小于等于K,若是,则另i=i+1,并继续基于方程找出索引值后添加索引值到被选原子索引集 若否,计算在 处的 近似值,在其它地方,其值为0,并将当前的vi作为荧光分子的的坐标向量,其中K为预估分子数, 为第K次迭代得到的索引集,为定位后的像素值。4.一种超分辨成像系统,其特征在于,包括:图像采集器,用于采集多幅由多个荧光分子通过光学系统的荧光图像;如权利要求3所述的荧光多分子定位装置,用于对每幅荧光图像中的荧光分子进行定位,得到其坐标位置;重构模块,用于对定位后的所有荧光图像进行累加,根据全部荧光分子的坐标位置描绘出超分辨荧光图像。
|